SERK2v2: A new second‐order stabilized explicit Runge‐Kutta method for stiff problems

Traditionally, explicit numerical algorithms have not been used with stiff ordinary differential equations (ODEs) due to their stability. Implicit schemes are usually very expensive when used to solve systems of ODEs with very large dimension. Stabilized Runge-Kutta methods (also called Runge–Kutta–Chebyshev methods) were proposed to try to avoid these difficulties. The Runge–Kutta methods are explicit methods with extended stability domains, usually along the negative real axis. They can easily be applied to large problem classes with low memory demand, they do not require algebra routines or the solution of large and complicated systems of nonlinear equations, and they are especially suited for discretizations using the method of lines of two and three dimensional parabolic partial differential equations. In Martin-Vaquero and Janssen [Comput Phys Commun 180 (2009), 1802–1810], we showed that previous codes based on stabilized Runge–Kutta algorithms have some difficulties in solving problems with very large eigenvalues and we derived a new code, SERK2, based on sixth-order polynomials. Here, we develop a new method based on second-order polynomials with up to 250 stages and good stability properties. These methods are efficient numerical integrators of very stiff ODEs. Numerical experiments with both smooth and nonsmooth data support the efficiency and accuracy of the new algorithms when compared to other well-known second-order methods such as RKC and ROCK2. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013

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

[2]  K. Burrage,et al.  Adaptive stepsize based on control theory for stochastic differential equations , 2004 .

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

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

[5]  Linda R. Petzold,et al.  Runge-Kutta-Chebyshev projection method , 2006, J. Comput. Phys..

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

[7]  Willem Hundsdorfer,et al.  RKC time-stepping for advection-diffusion-reaction problems , 2004 .

[8]  Kjell Gustafsson,et al.  Control-theoretic techniques for stepsize selection in implicit Runge-Kutta methods , 1991, TOMS.

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

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

[11]  K. Gustafsson,et al.  API stepsize control for the numerical solution of ordinary differential equations , 1988 .

[12]  Kjell Gustafsson,et al.  Control Strategies for the Iterative Solution of Nonlinear Equations in ODE Solvers , 1997, SIAM J. Sci. Comput..

[13]  Desmond J. Higham,et al.  Analysis of stepsize selection schemes for Runge-Kutta codes , 1988 .

[14]  A. Abdulle,et al.  S-ROCK methods for stiff Ito SDEs , 2008 .

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

[16]  C F Curtiss,et al.  Integration of Stiff Equations. , 1952, Proceedings of the National Academy of Sciences of the United States of America.

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

[18]  Jan G. Verwer,et al.  An Implicit-Explicit Runge-Kutta-Chebyshev Scheme for Diffusion-Reaction Equations , 2004, SIAM J. Sci. Comput..

[19]  P. Houwen,et al.  On the Internal Stability of Explicit, m‐Stage Runge‐Kutta Methods for Large m‐Values , 1979 .

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

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