Filtered Implicit Second-Derivative Time-Stepping Methods for Stiff Initial Value Problems

Abstract

This work focuses on the construction and numerical analysis of novel time-filtered second-derivative methods for stiff equations. The procedure is based on the recent first-derivative time filters (DeCaria et al., 2022 [6]). Such methods are developed by incorporating inexpensive pre-filtering and post-filtering steps into existing second-derivative schemes. We show that applying these filtering steps to second-derivative multistep or multi-stage methods results in new methods that combine features of both multistep and multi-stage approaches, referred to as second-derivative general linear methods (SGLMs). The well-established properties of SGLMs are utilized to analyze the accuracy and stability of the filtered methods and to design optimal new filters for time-stepping schemes. Several new embedded families of high-accuracy methods with low cognitive complexity and excellent stability characteristics are introduced. Finally, numerical experiments validate the stability and efficiency of the proposed methods.

Keywords

Time filters / General linear methods / Second-derivative methods / Stability / Implicit methods / 65L05

References

[1] Abdi A, Hojjati G. An extension of general linear methods. Numer. Algorithms, 2011, 57(2): 149-167.

[2] Abdi A, Hojjati G. Maximal order for second derivative general linear methods with runge-kutta stability. Appl. Numer. Math., 2011, 61(10): 1046-1058.

[3] Butcher JC, Hojjati G. Second derivative methods with rk stability. Numer. Algorithms, 2005, 40(4): 415-429.

[4] Chen Y, Shen C. A jacobian-free newton-gmres(m) method with adaptive preconditioner and its application for power flow calculations. IEEE Trans. Power Syst., 2006, 21(3): 1096-1103.

[5] Chouchoulis J, Schütz J. Jacobian-free implicit mdrk methods for stiff systems of odes. Appl. Numer. Math., 2024, 196: 45-61.

[6] DeCaria V, Gottlieb S, Grant ZJ, Layton WJ. A general linear method approach to the design and optimization of efficient, accurate, and easily implemented time-stepping methods in cfd. J. Comput. Phys., 2022, 455110927

[7] DeCaria V, Guzel A, Layton W, Li Y. A variable stepsize, variable order family of low complexity. SIAM J. Sci. Comput., 2021, 43(3): A2130-A2160.

[8] DeCaria, V., Layton, W., Zhao, H.: A time-accurate, adaptive discretization for fluid flow problems. arXiv:1810.06705 (2018)

[9] DeCaria V, Schneier M. An embedded variable step imex scheme for the incompressible navier-stokes equations. Comput. Methods Appl. Mech. Eng., 2021, 376113661

[10] Gottlieb S, Grant ZJ, Hu J, Shu R. High order strong stability preserving multiderivative implicit and imex runge-kutta methods with asymptotic preserving properties. SIAM J. Numer. Anal., 2022, 60(1): 423-449.

[11] Hairer E, Wanner G. Multistep-multistage-multiderivative methods for ordinary differential equations. Computing, 1973, 11(3): 287-303.

[12] Hairer, E., Wanner, G.: Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems. Second edition. Springer-Verlag, Berlin (1996)

[13] Hojjati G, Ardabili MR, Hosseini S. A-ebdf: an adaptive method for numerical solution of stiff systems of odes. Math. Comput. Simul., 2004, 66(1): 33-41.

[14] Hojjati G, Ardabili MR, Hosseini S. New second derivative multistep methods for stiff systems. Appl. Math. Model., 2006, 30(5): 466-476.

[15] Knoll DA, Keyes DE. Jacobian-free newton-krylov methods: a survey of approaches and applications. J. Comput. Phys., 2004, 193(2): 357-397.

[16] Moradi A, Abdi A, Hojjati G. Strong stability preserving implicit and implicit-explicit second derivative general linear methods with rk stability. Comput. Appl. Math., 2022, 41(4): 135.

[17] Moradi A, Farzi J, Abdi A. Order conditions for second derivative general linear methods. J. Comput. Appl. Math., 2021, 387112488

[18] Niemeyer KE, Curtis NJ, Sung C-J. Pyjac: analytical jacobian generator for chemical kinetics. Comput. Phys. Commun., 2017, 215: 188-203.

[19] Perini F, Galligani E, Reitz RD. An analytical jacobian approach to sparse reaction kinetics for computationally efficient combustion modeling with large reaction mechanisms. Energy & Fuels, 2012, 26(8): 4804-4822.

[20] Schütz J, Seal DC, Zeifang J. Parallel-in-time high-order multiderivative imex solvers. J. Sci. Comput., 2022, 90(1): 54.

[21] Shampine LF, Reichelt MW. The matlab ode suite. SIAM J. Sci. Comput., 1997, 18(1):1-22.

Funding

Otto-von-Guericke-Universität Magdeburg (3121)

Rights & permissions

The Author(s)