On Richardson extrapolation for low-dissipation low-dispersion diagonally implicit Runge-Kutta schemes

Abstract In the modeling of wave propagation phenomena it is necessary to use time integration methods which are not only sufficiently accurate, but also properly describe the amplitude and phase of the propagating waves. It is not clear if amending the developed schemes by extrapolation methods to obtain a high order of accuracy preserves the qualitative properties of these schemes in the perspective of dissipation, dispersion and stability analysis. It is illustrated that the combination of various optimized schemes with Richardson extrapolation is not optimal for minimal dissipation and dispersion errors. Optimized third-order and fourth-order methods are obtained, and it is shown that the proposed methods combined with Richardson extrapolation result in fourth and fifth orders of accuracy correspondingly, while preserving optimality and stability. The numerical applications include the linear wave equation, a stiff system of reaction–diffusion equations and the nonlinear Euler equations with oscillatory initial conditions. It is demonstrated that the extrapolated third-order scheme outperforms the recently developed fourth-order diagonally implicit Runge–Kutta scheme in terms of accuracy and stability.

[1]  Ben P. Sommeijer,et al.  Diagonally implicit Runge-Kutta-Nystrm methods for oscillatory problems , 1989 .

[2]  Manuel Calvo,et al.  Short note: a new minimum storage Runge-Kutta scheme for computational acoustics , 2004 .

[3]  C. Bogey,et al.  Low-dissipation and low-dispersion fourth-order Runge–Kutta algorithm , 2006 .

[4]  W. Hundsdorfer,et al.  A note on iterated splitting schemes , 2007 .

[5]  Zahari Zlatev,et al.  The convergence of diagonally implicit Runge-Kutta methods combined with Richardson extrapolation , 2013, Comput. Math. Appl..

[6]  M.Y. Hussaini,et al.  Low-Dissipation and Low-Dispersion Runge-Kutta Schemes for Computational Acoustics , 1994 .

[7]  J. Butcher Numerical methods for ordinary differential equations , 2003 .

[8]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[9]  Miguel R. Visbal,et al.  An implicit discontinuous Galerkin method for the unsteady compressible Navier–Stokes equations , 2009 .

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

[11]  W. Habashi,et al.  2N-Storage Low Dissipation and Dispersion Runge-Kutta Schemes for Computational Acoustics , 1998 .

[12]  Christophe Bailly,et al.  Optimized explicit schemes: matching and boundary schemes and 4th-order Runge-Kutta algorithm. , 2004 .

[13]  Xiaolong Gou,et al.  A dynamic multi-timescale method for combustion modeling with detailed and reduced chemical kinetic mechanisms , 2010 .

[14]  Z. Zlatev,et al.  Richardson Extrapolation combined with the sequential splitting procedure and the θ-method , 2012 .

[15]  Cornelis Vuik,et al.  Numerical Methods for Reacting Gas Flow Simulations , 2007 .

[16]  J. Stoer,et al.  Numerical treatment of ordinary differential equations by extrapolation methods , 1966 .

[17]  Zhen-Xing Yao,et al.  Optimized explicit finite-difference schemes for spatial derivatives using maximum norm , 2013, J. Comput. Phys..

[18]  S. Lele Compact finite difference schemes with spectral-like resolution , 1992 .

[19]  Farshid Nazari,et al.  High-order low-dissipation low-dispersion diagonally implicit Runge-Kutta schemes , 2015, J. Comput. Phys..

[20]  C. Bogey,et al.  A family of low dispersive and low dissipative explicit schemes for flow and noise computations , 2004 .

[21]  P. Deuflhard,et al.  A semi-implicit mid-point rule for stiff systems of ordinary differential equations , 1983 .

[22]  M H Chawla,et al.  A Noumerov-type method with minimal phase-lag for the integration of second order periodic initial-value , 1986 .

[23]  Christophe Bailly,et al.  Semi-Implicit Runge-Kutta Schemes: Development and Application to Compressible Channel Flow , 2014 .

[24]  Luc Mongeau,et al.  A low-dispersion and low-dissipation implicit Runge-Kutta scheme , 2013, J. Comput. Phys..

[25]  Li Wang,et al.  Implicit solution of the unsteady Euler equations for high-order accurate discontinuous Galerkin discretizations , 2006, J. Comput. Phys..

[26]  L. Richardson The Approximate Arithmetical Solution by Finite Differences of Physical Problems Involving Differential Equations, with an Application to the Stresses in a Masonry Dam , 1911 .

[27]  Tom E. Simos,et al.  Runge-Kutta methods with minimal dispersion and dissipation for problems arising from computational acoustics , 2005 .

[28]  Chi-Wang Shu,et al.  High Order Weighted Essentially Nonoscillatory Schemes for Convection Dominated Problems , 2009, SIAM Rev..

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

[30]  C. Tam,et al.  Dispersion-relation-preserving finite difference schemes for computational acoustics , 1993 .

[31]  G. Marchuk,et al.  Difference Methods and Their Extrapolations , 1983 .