×

High-order implicit time-marching methods based on generalized summation-by-parts operators. (English) Zbl 1359.65127

Summary: 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 favorable for stiff initial value problems, such as L-stability. It can facilitate the analysis of fully discrete approximations to PDEs and is amenable to multidimensional space-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.

MSC:

65L12 Finite difference and finite volume methods for ordinary differential equations
65L04 Numerical methods for stiff equations
65L05 Numerical methods for initial value problems involving ordinary differential equations
65L06 Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations
65L20 Stability and convergence of numerical methods for ordinary differential equations
65L60 Finite element, Rayleigh-Ritz, Galerkin and collocation methods for ordinary differential equations
PDFBibTeX XMLCite
Full Text: DOI arXiv

References:

[1] O. Axelsson, {\it A note on a class of strongly A-stable methods}, BIT, 12 (1972), pp. 1-4. · Zbl 0236.65047
[2] P. D. Boom, {\it High-order Implicit Numerical Methods for Unsteady Fluid Simulation}, Ph.D. thesis, University of Toronto Institute for Aerospace Studies, Toronto, 2015.
[3] K. Burrage and J. C. Butcher, {\it Non-linear stability of a general class of differential equation methods}, BIT, 20 (1980), pp. 185-203. · Zbl 0431.65051
[4] J. C. Butcher, {\it Implicit Runge-Kutta processes}, Math. Comp., 18 (1964), pp. 50-64. · Zbl 0123.11701
[5] J. C. Butcher, {\it A stability property of implicit Runge-Kutta methods}, BIT, 15 (1975), pp. 358-361. · Zbl 0333.65031
[6] J. C. Butcher, {\it Numerical Methods for Ordinary Differential Equations}, Wiley, Chichester, 2003. · Zbl 1040.65057
[7] M. H. Carpenter and D. Gottlieb, {\it Spectral methods on arbitrary grids}, J. Comput. Phys., 129 (1996), pp. 74-86. · Zbl 0862.65054
[8] M. H. Carpenter, D. Gottlieb, and S. S. Abarbanel, {\it Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: Methodology and application to high-order compact schemes}, J. Comput. Phys., 111 (1994), pp. 220-236. · Zbl 0832.65098
[9] M. H. Carpenter, J. Nordström, and D. Gottlieb, {\it A stable and conservative interface treatment of arbitrary spatial accuracy}, J. Comput. Phys., 148 (1999), pp. 341-365. · Zbl 0921.65059
[10] F. Chipman, {\it A-stable Runge-Kutta processes}, BIT, 11 (1971), pp. 384-388. · Zbl 0265.65035
[11] M. Crouzeix, {\it Sur la B-stabilité des méthodes de Runge-Kutta}, Numer. Math., 32 (1979), pp. 75-82. · Zbl 0431.65052
[12] G. G. Dahlquist, {\it A special stability problem for linear multistep methods}, BIT, 3 (1963), pp. 27-43. · Zbl 0123.11703
[13] K. Dekker and J. G. Verwer, {\it Stability of Runge-Kutta Methods for Stiff Nonlinear Differential Equations}, CWI Mongr., 1984. · Zbl 0571.65057
[14] D. C. Del Rey Fernández, P. D. Boom, and D. W. Zingg, {\it A generalized framework for nodal first derivative summation-by-parts operators}, J. Comput. Phys., 266 (2014), pp. 214-239. · Zbl 1311.65002
[15] B. L. Ehle, {\it On Padé Approximations to the Exponential Function and A-stable Methods for the Numerical Solution of Initial Value Problems}, Ph.D. thesis, University of Waterloo, Waterloo, Ontario, 1969.
[16] D. Funaro and D. Gottlieb, {\it A new method of imposing boundary conditions in pseudospectral approximations of hyperbolic equations}, Math. Comp., 51 (1988), pp. 599-613. · Zbl 0699.65079
[17] D. Funaro and D. Gottlieb, {\it Convergence results for pseudospectral approximations of hyperbolic systems by a penalty-type boundary treatment}, Math. Comp., 57 (1991), pp. 585-596. · Zbl 0736.65074
[18] G. J. Gassner, {\it A Skew-symmetric discontinuous Galerkin spectral element discretization and its relation to SBP-SAT finite difference methods}, SIAM J. Sci. Comput., 35 (2013). · Zbl 1275.65065
[19] B. Gustafsson, {\it High Order Difference Methods for Time Dependent PDE}, Springer, Berlin, 2008. · Zbl 1146.65064
[20] B. Gustafsson, H.-O. Kreiss, and J. Oliger, {\it Time Dependent Problems and Difference Methods}, Wiley, Chichester, 1996. · Zbl 0843.65061
[21] E. Hairer, {\it Highest possible order of algebraically stable diagonally implicit Runge-Kutta methods}, BIT, 20 (1980), pp. 254-256. · Zbl 0425.65037
[22] E. Hairer, C. Lubich, and M. Roche, {\it Error of Runge-Kutta methods for stiff problems studied via differential algebraic equations}, BIT, 28 (1988), pp. 678-700. · Zbl 0657.65093
[23] E. Hairer, S. Nørsett, and G. Wanner, {\it Solving Ordinary Differential Equations I: Nonstiff Problems}, 2nd ed., Springer, Berlin, 2000. · Zbl 1185.65115
[24] E. Hairer and G. Wanner, {\it Solving Ordinary Differential Equations II}, Springer, Berlin, 1991. · Zbl 0729.65051
[25] E. Hairer, G. Wanner, and C. Lubich, {\it Geometric Numerical Integration–Structure-Preserving Algorithms for Ordinary Differential Equations}, Springer Ser. Comput. Math. 31, Springer, Berlin, 2006. · Zbl 1094.65125
[26] J. E. Hicken, D. C. Del Rey Fernándex, and D. W. Zingg, {\it Multidimensional Summation-By-Parts Operators: General Theory and Application to Simplex Elements}, J. Comput. Phys., submitted. · Zbl 1382.65355
[27] J. E. Hicken and D. W. Zingg, {\it Superconvergent functional estimates from summation-by-parts finite-difference discretizations}, SIAM J. Sci. Comput., 33 (2011), pp. 893-922. · Zbl 1227.65102
[28] J. E. Hicken and D. W. Zingg, {\it Summation-by-parts operators and high-order quadrature}, J. Comput. Appl. Math., 237 (2013), pp. 111-125. · Zbl 1263.65025
[29] J. E. Hicken and D. W. Zingg, {\it Dual consistency and functional accuracy: A finite-difference perspective}, J. Comput. Phys., 256 (2014), pp. 161-182. · Zbl 1349.65559
[30] Z. Jackiewicz, {\it General Linear Methods for Ordinary Differential Equations}, Wiley, Chichester, 2009. · Zbl 1211.65095
[31] H.-O. Kreiss and J. Lorenz, {\it Initial-Boundary Value Problems and the Navier-Stokes Equations}, Classics Appl. Math. 47, SIAM, Philadelphia 2004. · Zbl 1097.35113
[32] H.-O. Kreiss and G. Scherer, {\it Finite element and finite difference methods for hyperbolic partial differential equations}, in Mathematical Aspects of Finite Elements in Partial Differential Equations, Academic Press, New York, 1974, pp. 195-212. · Zbl 0355.65085
[33] J. Kuntzmann, {\it Neuere entwicklungen der methode von Runge und Kutta}, ZAMM Z. Angew. Math. Mech., 41 (1961), pp. T28-T31. · Zbl 0106.10403
[34] H. Lomax, T. H. Pulliam, and D. W. Zingg, {\it Fundamentals of Computational Fluid Dynamics}, Sci. Comput., Springer, Berlin, 2001. · Zbl 0970.76002
[35] J. Lu, {\it An {\it A Posteriori} Error Control Framework for Adaptive Precision Optimization Using Discontinuous Galerkin Finite Element Method}, Ph.D. thesis, Massachusetts Institute of Technology, Cambridge, MA, 2005.
[36] T. Lundquist and J. Nordström, {\it The SBP-SAT technique for initial value problems}, J. Comput. Phys., 270 (2014), pp. 86-104. · Zbl 1349.65550
[37] J. Nordström and T. Lundquist, {\it Summation-by-parts in time}, J. Comput. Phys., 251 (2013), pp. 487-499. · Zbl 1349.65399
[38] P. Prince and J. Dormand, {\it High order embedded Runge-Kutta formulae}, J. Comput. Appl. Math., 7 (1981), pp. 67-75. · Zbl 0449.65048
[39] A. Prothero and A. Robinson, {\it On the stability and accuracy of one-step methods for solving stiff systems of ordinary differential equations}, Math. Comp., 28 (1974), pp. 145-162. · Zbl 0309.65034
[40] B. Strand, {\it Summation by parts for finite difference approximations for d/dx}, J. Comput. Phys., 110 (1994), pp. 47-67. · Zbl 0792.65011
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.