High-Order Implicit Time-Marching Methods Based on Generalized Summation-By-Parts Operators

This article extends the theory of classical finite-difference summation-by-parts (FD-SBP) time-marching methods to the generalized summation-by-parts (GSBP) framework. Dual-consistent GSBP time-marching methods are shown to retain: A and L-stability, as well as superconvergence of integral functionals when integrated with the quadrature associated with the discretization. This also implies that the solution approximated at the end of each time step is superconvergent. In addition GSBP time-marching methods constructed with a diagonal norm are BN-stable. This article also formalizes the connection between FD-SBP/GSBP time-marching methods and implicit Runge-Kutta methods. Through this connection, the minimum accuracy of the solution approximated at the end of a time step is extended for nonlinear problems. It is also exploited to derive conditions under which nonlinearly stable GSBP time-marching methods can be constructed. The GSBP approach to time marching can simplify the construction of high-order fully-implicit Runge-Kutta methods with a particular set of properties favourable for stiff initial value problems, such as L-stability. It can facilitate the analysis of fully discrete approximations to PDEs and is amenable to to multi-dimensional spcae-time discretizations, in which case the explicit connection to Runge-Kutta methods is often lost. A few examples of known and novel Runge-Kutta methods associated with GSBP operators are presented. The novel methods, all of which are L-stable and BN-stable, include a four-stage seventh-order fully-implicit method, a three-stage third-order diagonally-implicit method, and a fourth-order four-stage diagonally-implicit method. The relative efficiency of the schemes is investigated and compared with a few popular non-GSBP Runge-Kutta methods.

[1]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[2]  David C. Del Rey Fernández,et al.  A generalized framework for nodal first derivative summation-by-parts operators , 2014, J. Comput. Phys..

[3]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[4]  D. Gottlieb,et al.  A Stable and Conservative Interface Treatment of Arbitrary Spatial Accuracy , 1999 .

[5]  John C. Butcher,et al.  A stability property of implicit Runge-Kutta methods , 1975 .

[6]  K. Burrage,et al.  Non-linear stability of a general class of differential equation methods , 1980 .

[7]  S. Abarbanel,et al.  Asymptotically Stable Fourth-Order Accurate Schemes for the Diffusion Equation on Complex Shapes , 1997 .

[8]  Jason E. Hicken,et al.  Summation-by-parts operators and high-order quadrature , 2011, J. Comput. Appl. Math..

[9]  David Gottlieb,et al.  Spectral Methods on Arbitrary Grids , 1995 .

[10]  John C. Butcher,et al.  General linear methods for ordinary differential equations , 2009, Math. Comput. Simul..

[11]  Ch. Hirsch,et al.  Fundamentals Of Computational Fluid Dynamics , 2016 .

[12]  P. Boom,et al.  High-order implicit time-marching methods for unsteady fluid flow simulation , 2015 .

[13]  David Gottlieb,et al.  CONVERGENCE RESULTS FOR PSEUDOSPECTRAL APPROXIMATIONS OF HYPERBOLIC SYSTEMS BY A PENALTY-TYPE BOUNDARY TREATMENT , 1991 .

[14]  E. Hairer,et al.  Solving Ordinary Differential Equations II , 2010 .

[15]  M. Carpenter,et al.  Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations , 2003 .

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

[17]  Daniele,et al.  CONVERGENCE RESULTS FOR PSEUDOSPECTRAL APPROXIMATIONS OF HYPERBOLIC SYSTEMS BY A PENALTY TYPE BOUNDARY TREATMENT , 1991 .

[18]  Ernst Hairer,et al.  Highest possible order of algebraically stable diagonally implicit runge-kutta methods , 1980 .

[19]  Jason E. Hicken,et al.  Superconvergent Functional Estimates from Summation-By-Parts Finite-Difference Discretizations , 2011, SIAM J. Sci. Comput..

[20]  H. Kreiss,et al.  Time-Dependent Problems and Difference Methods , 1996 .

[21]  Ken Mattsson,et al.  Optimal diagonal-norm SBP operators , 2014, J. Comput. Phys..

[22]  Gregor Gassner,et al.  A Skew-Symmetric Discontinuous Galerkin Spectral Element Discretization and Its Relation to SBP-SAT Finite Difference Methods , 2013, SIAM J. Sci. Comput..

[23]  G. Dahlquist.  A special stability problem for linear multistep methods , 1963 .

[24]  A. Prothero,et al.  On the stability and accuracy of one-step methods for solving stiff systems of ordinary differential equations , 1974 .

[25]  Jan Nordström,et al.  Summation-by-parts in time , 2013, J. Comput. Phys..

[26]  D. Gottlieb,et al.  A new method of imposing boundary conditions in pseudospectral approximations of hyperbolic equations , 1988 .

[27]  J. Butcher.  Implicit Runge-Kutta processes , 1964 .

[28]  D. Gottlieb,et al.  Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes , 1994 .

[29]  H. Kreiss,et al.  Finite Element and Finite Difference Methods for Hyperbolic Partial Differential Equations , 1974 .

[30]  David C. Del Rey Fernández,et al.  Multidimensional Summation-by-Parts Operators: General Theory and Application to Simplex Elements , 2015, SIAM J. Sci. Comput..

[31]  H. Kreiss,et al.  Initial-Boundary Value Problems and the Navier-Stokes Equations , 2004 .

[32]  P. Roache.  Fundamentals of computational fluid dynamics , 1998 .

[33]  J. Kuntzmann,et al.  Neuere Entwicklungen der Methode von Runge und Kutta , 1961 .

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

[35]  J. Dormand,et al.  High order embedded Runge-Kutta formulae , 1981 .

[36]  David W. Zingg,et al.  Runge-Kutta Characterization of the Generalized Summation-by-Parts Approach in Time , 2014, 1410.0202.

[37]  Jason E. Hicken,et al.  Dual consistency and functional accuracy: a finite-difference perspective , 2014, J. Comput. Phys..

[38]  F. Chipman.  A-stable Runge-Kutta processes , 1971 .

[39]  James Lu,et al.  An a posteriori Error Control Framework for Adaptive Precision Optimization using Discontinuous Galerkin Finite Element Method , 2005 .

[40]  Ernst Hairer,et al.  Error of Runge-Kutta methods for stiff problems studied via differential algebraic equations , 1988 .

[41]  K. Burrage,et al.  Stability Criteria for Implicit Runge–Kutta Methods , 1979 .

[42]  B. Gustafsson.  High Order Difference Methods for Time Dependent PDE , 2008 .

[43]  B. Strand.  Summation by parts for finite difference approximations for d/dx , 1994 .

[44]  M. Crouzeix.  Sur laB-stabilité des méthodes de Runge-Kutta , 1979 .

[45]  J. Verwer,et al.  Stability of Runge-Kutta Methods for Stiff Nonlinear Differential Equations , 1984 .

[46]  Jan Nordström,et al.  The SBP-SAT technique for initial value problems , 2014, J. Comput. Phys..

[47]  O. Axelsson.  A note on a class of stronglyA-stable methods , 1972 .