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] AbdiA, HojjatiG. An extension of general linear methods. Numer. Algorithms, 2011, 57(2): 149-167.

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

[3] ButcherJC, HojjatiG. Second derivative methods with rk stability. Numer. Algorithms, 2005, 40(4): 415-429.

[4] ChenY, ShenC. 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] ChouchoulisJ, SchützJ. Jacobian-free implicit mdrk methods for stiff systems of odes. Appl. Numer. Math., 2024, 196: 45-61.

[6] DeCariaV, GottliebS, GrantZJ, LaytonWJ. 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, 455. 110927

[7] DeCariaV, GuzelA, LaytonW, LiY. 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] DeCariaV, SchneierM. An embedded variable step imex scheme for the incompressible navier-stokes equations. Comput. Methods Appl. Mech. Eng., 2021, 376. 113661

[10] GottliebS, GrantZJ, HuJ, ShuR. 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] HairerE, WannerG. 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] HojjatiG, ArdabiliMR, HosseiniS. A-ebdf: an adaptive method for numerical solution of stiff systems of odes. Math. Comput. Simul., 2004, 66(1): 33-41.

[14] HojjatiG, ArdabiliMR, HosseiniS. New second derivative multistep methods for stiff systems. Appl. Math. Model., 2006, 30(5): 466-476.

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

[16] MoradiA, AbdiA, HojjatiG. Strong stability preserving implicit and implicit-explicit second derivative general linear methods with rk stability. Comput. Appl. Math., 2022, 414135.

[17] MoradiA, FarziJ, AbdiA. Order conditions for second derivative general linear methods. J. Comput. Appl. Math., 2021, 387. 112488

[18] NiemeyerKE, CurtisNJ, SungC-J. Pyjac: analytical jacobian generator for chemical kinetics. Comput. Phys. Commun., 2017, 215: 188-203.

[19] PeriniF, GalliganiE, ReitzRD. 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ützJ, SealDC, ZeifangJ. Parallel-in-time high-order multiderivative imex solvers. J. Sci. Comput., 2022, 90154.

[21] ShampineLF, ReicheltMW. 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)

Just Accepted

This article has successfully passed peer review and final editorial review, and will soon enter typesetting, proofreading and other publishing processes. The currently displayed version is the accepted final manuscript. The officially published version will be updated with format, DOI and citation information upon launch. We recommend that you pay attention to subsequent journal notifications and preferentially cite the officially published version. Thank you for your support and cooperation.