ESERK5: A fifth-order extrapolated stabilized explicit Runge-Kutta method

Abstract A new algorithm is developed and analyzed for multi-dimensional non-linear parabolic partial differential equations (PDEs) which are semi-discretized in the spatial variables leading to a system of ordinary differential equations (ODEs). It is based on fifth-order extrapolated stabilized explicit Runge–Kutta schemes (ESERK). They are explicit methods, and therefore it is not necessary to employ complicated software for linear or non-linear system of equations. Additionally, they have extended stability regions along the negative real semi-axis, hence they can be considered to solve stiff problems coming from very common diffusion or reaction–diffusion problems. Previously, only lower-order codes (up to fourth-order) have been constructed and made available in the scientific literature. However, at the same time, higher-order codes were demonstrated to be very efficient to solve equations where it is necessary to have a high precision or they have transient zones that are very severe, and where functions change very fast. The new schemes allow changing the step length very easily and with a very small computational cost. Thus, a variable step length, with variable number of stages algorithm is constructed and compared with good numerical results in relation to other well-known ODE solvers.

[1]  L. Shampine,et al.  RKC: an explicit solver for parabolic PDEs , 1998 .

[2]  P. Gremaud,et al.  Super-time-stepping acceleration of explicit schemes for parabolic problems , 1996 .

[3]  L. M. Skvortsov.  Explicit stabilized Runge-Kutta methods , 2011 .

[4]  V. I. LEBEDEV.  A new method for determining the roots of polynomials of least deviation on a segment with weight and subject to additional conditions. Part I , 1993 .

[5]  Steven J. Ruuth,et al.  Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations , 1997 .

[6]  M. Yousuf,et al.  Smoothing schemes for reaction-diffusion systems with nonsmooth data , 2009 .

[7]  Marlis Hochbruck,et al.  Exponential Integrators for Large Systems of Differential Equations , 1998, SIAM J. Sci. Comput..

[8]  J. Martín-Vaquero,et al.  Extrapolated stabilized explicit Runge-Kutta methods , 2016, J. Comput. Phys..

[9]  Conall O'Sullivan,et al.  On the acceleration of explicit finite difference methods for option pricing , 2011 .

[10]  Marc Massot,et al.  A new numerical strategy with space-time adaptivity and error control for multi-scale streamer discharge simulations , 2011, J. Comput. Phys..

[11]  J. Verwer,et al.  Numerical solution of time-dependent advection-diffusion-reaction equations , 2003 .

[12]  J. Martín-Vaquero,et al.  SERK2v2: A new second‐order stabilized explicit Runge‐Kutta method for stiff problems , 2013 .

[13]  Laurent Stainier,et al.  Combined implicit/Explicit time-integration algorithms for the numerical simulation of sheet metal forming , 2004 .

[14]  Jeff Cash,et al.  The integration of stiff initial value problems in ODEs using modified extended backward differentiation formulae , 1983 .

[15]  Assyr Abdulle,et al.  S-ROCK: Chebyshev Methods for Stiff Stochastic Differential Equations , 2008, SIAM J. Sci. Comput..

[16]  T. Downes,et al.  A three-dimensional numerical method for modelling weakly ionized plasmas , 2006, astro-ph/0612580.

[17]  Richard C Aiken,et al.  Stiff computation , 1985 .

[18]  David I. Ketcheson,et al.  Optimal stability polynomials for numerical integration of initial value problems , 2012, 1201.3035.

[19]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[20]  V. I. Lebedev,et al.  Solution of the parameter ordering problem in Chebyshev iterative methods , 1974 .

[21]  A. Hindmarsh.  LSODE and LSODI, two new initial value ordinary differential equation solvers , 1980, SGNM.

[22]  J. Martín-Vaquero,et al.  Second-order stabilized explicit Runge-Kutta methods for stiff problems , 2009, Comput. Phys. Commun..

[23]  Assyr Abdulle,et al.  PIROCK: A swiss-knife partitioned implicit-explicit orthogonal Runge-Kutta Chebyshev integrator for stiff diffusion-advection-reaction problems with or without noise , 2013, J. Comput. Phys..

[24]  Abdul-Qayyum M. Khaliq,et al.  Stabilized explicit Runge-Kutta methods for multi-asset American options , 2014, Comput. Math. Appl..

[25]  J. Martín-Vaquero,et al.  SERK2v3: Solving mildly stiff nonlinear partial differential equations , 2016, J. Comput. Appl. Math..

[26]  Marc Massot,et al.  Time–space adaptive numerical methods for the simulation of combustion fronts , 2013 .

[27]  John C. Butcher,et al.  ARK methods for stiff problems , 2005 .

[28]  Abdul-Qayyum M. Khaliq,et al.  The locally extrapolated exponential time differencing LOD scheme for multidimensional reaction-diffusion systems , 2015, J. Comput. Appl. Math..

[29]  J. Verwer.  Explicit Runge-Kutta methods for parabolic partial differential equations , 1996 .

[30]  Stephen O'Sullivan,et al.  A class of high-order Runge-Kutta-Chebyshev stability polynomials , 2015, J. Comput. Phys..

[31]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[32]  Marc Massot,et al.  Simulation of human ischemic stroke in realistic 3D geometry , 2010, Commun. Nonlinear Sci. Numer. Simul..

[33]  Jesús Vigo-Aguiar,et al.  Exponential fitting BDF algorithms: explicit and implicit 0-stable methods , 2006 .

[34]  Manuel Torrilhon,et al.  Essentially optimal explicit Runge–Kutta methods with application to hyperbolic–parabolic equations , 2007, Numerische Mathematik.

[35]  A. Medovikov.  High order explicit methods for parabolic equations , 1998 .

[36]  H. De Meyer,et al.  Frequency evaluation in exponential fitting multistep algorithms for ODEs , 2002 .

[37]  Assyr Abdulle,et al.  Second order Chebyshev methods based on orthogonal polynomials , 2001, Numerische Mathematik.

[38]  Assyr Abdulle,et al.  Fourth Order Chebyshev Methods with Recurrence Relation , 2001, SIAM J. Sci. Comput..

[39]  David I. Ketcheson,et al.  A comparison of high-order explicit Runge–Kutta, extrapolation, and deferred correction methods in serial and parallel , 2013, 1305.6165.