Exponential Integrators

 

Exponential integrators constitute an efficient tool for the numerical solution of stiff and highly oscillatory problems. They rely on the variation-of-constants formula and require the evaluation of the action of certain matrix functions. Such computations can be carried out efficiently on massively parallel systems.


Introduction:


Many problems from science and engineering are modelled by partial differential equations. Their solutions describes the temporal evolution of the modelled processes. In most cases, however, the arising equations are too complex to be studied analytically. Consequently, their solutions have to be approximated by numerical methods. In this regard, topics like accuracy, numerical stability, efficiency and required computer memory are of paramount importance.

Exponential integrators constitute a very efficient class of numerical methods for the solution of high-dimensional stiff or highly oscillatory problems. Therefore, they are perfectly suited for carrying out complex simulations.


Exponential Integrators:


Consider the system of stiff or highly oscillatory differential equations u ' ( t ) = F ( u ( t ) ) , u ( t 0 ) = u 0 with a sufficiently smooth solution u . Such systems arise from spatial discretizations of parabolic or hyperbolic partial differential equations. For its numerical solution we consider discrete times t 1 , , t N and compute for each of these times t n an approximation u n to the exact solution u ( t n ) . Exponential integrators carry out these computations in an explicit way. A simple example is given by the exponential Rosenbrock-Euler method, which is a method of order 2: u n + 1 = u n + h n + ϕ 1 ( h n J n ) F ( u n ) with ϕ 1 = e z - 1 z . Here, J n denotes an appropriate approximation to the Jacobian matrix F ' ( u n ) . The implementation of this method requires the evaluation of the action of the matrix function ϕ 1 ( h n J n )  on a vector which we perform by interpolation based on Leja points. Note that exponential integrators of arbitrarily high order exist. Their properties are studied in this research project.

   
© Numerical Analysis - University of Innsbruck 2015