On overcoming Dahlquist’s second barrier for$A$-stable linear multistep methods

Document Type : Research Article

Authors

1 Faculty of Mathematics, Statistics and Computer Science, University of Tabriz, Tabriz, Iran.

2 Marand Technical College, University of Tabriz, Tabriz, Iran.

3 Institute of Analysis and Numerics, Otto von Guericke University Magdeburg, Universitätsplatz 2, 39106 Magdeburg, Germany.

Abstract

Dahlquist’s second barrier limits the order of $A$-stable linear multistep methods to at most two, posing significant challenges for achieving higher accuracy in the numerical solution of stiff ordinary differential equations. Leveraging various successful techniques, many efforts have been made to develop efficient methods that overcome this fundamental obstacle through different approaches. In this paper, we survey these techniques and analyze their impact on enhancing the stability and accuracy of the resulting methods. A comprehensive understanding of these advances can assist researchers in designing more effective algorithms for stiff problems.

Keywords

Main Subjects


[1] Abdi, A. and Hojjati, G. An extension of general linear methods, Numer. Algorithms, 57 (2011), 149–167.
[2] Butcher, J.C., and Hojjati, G. Second derivative methods with RK stability, Numer. Algorithms, 40 (2005), 415–429.
[3] Butcher, J.C. Numerical methods for ordinary differential equations, Wiley, Chichester, 2016.
[4] Cash, J.R. On the integration of stiff systems of ODEs using extended backward differentiation formulas, Numer. Math. 34(2) (1980), 235–246.
[5] Cash, J.R. Second derivative extended backward differentiation formulas for the numerical integration of stiff systems, SIAM J. Numer. Anal. 18(2) (1981), 21–36.
[6] Cash, J.R. The integration of stiff initial value problems in ODEs using modiffied extended backward differentiation formulae, Comput. Math. Appl. 9 (1983), 645–660.
[7] Cong, N.H. and Thuy, N.T. Stability of Two-Step-by-Two-Step IRK Methods Based on Gauss-Legendre Collocation Points and an Application, Vietnam J. Math. 40(1) (2012), 115–126.
[8] Curtiss, C.F. and Hirschfelder, J.O. Integration of stiff equations, Proceedings of the National Academy of Sciences, 38 (1952), 235–243.
[9] Dahlquist, G. A special stability problem for linear multistep methods, BIT Numer. Math. 3 (1963) 27–43.
[10] Enright, W.H. Second derivative multistep methods for stiff ordinary differential equations, SIAM J. Numer. Anal. 11 (1974), 321–331.
[11] Fazeli, S., Hojjati, G., and Shahmorad, S. Super implicit multistep collocation methods for nonlinear Volterra integral equations, Math. Comput. Model. 55 (2012) 590–607.
[12] Fredebeul, C. A-BDF: A generalization of the backward differentiation formulae, SIAM J. Numer. Anal. 35(5) (1998), 1917–1938.
[13] Gear, C.W. Numerical Initial Value Problems in Ordinary Differential Equations, Prentice-Hall, 1967.
[14] Hairer, E. and Wanner, G. Solving ordinary differential equations II: Stiff and differential–algebraic problems, Springer, Berlin, 2010.
[15] Hojjati, G. A class of parallel methods with superfuture points technique for the numerical solution of stiff systems, J. Mod. Meth. Numer. Math. 6 (2015), 57–63.
[16] Hojjati, G., Rahimi Ardabili, M.Y., and Hosseini, S.M. A-EBDF: an adaptive method for numerical solution of stiff systems of ODEs, Math. Comput. Simul. 66 (2004), 33–41.
[17] Hojjati, G., Rahimi Ardabili, M.Y. and Hosseini, S.M. New second derivative multistep methods for stiff systems, Appl. Math. Model. 30 (2006), 466–476.
[18] Hojjati, G. and Taheri Koltape, L. On the stability functions of second derivative implicit advanced-step point methods, J. Math. Model. 10 (2022), 203–212.
[19] Hindmarsh, A.C. ODEPACK, a systematized collection of ODE solvers, Scientific Computing, (1983), 55–64.
[20] Iserles, A. A first course in the numerical analysis of differential equations, Cambridge University Press, 1996.
[21] Jackiewicz, Z. General Linear Methods for Ordinary Differential Equations, Wiley, New Jersey, 2009.
[22] Psihoyios, G. Advanced step-point methods for the solution of initial value problems, Ph.D. Thesis, University of London, Imperial College, 1995.
[23] Psihoyios, G. A general formula for the stability functions of a group of implicit advanced step-point (IAS) methods, Math. Comput. Model. 46 (2007), 214–224.
[24] Shampine, L.F. and Reichelt, M.W. The MATLAB ODE suite, SIAM J. Sci. Comput. 18(1) (1997), 1–22.
[25] Skeel, R.D., Kong, A.K. Blended linear multistep methods, ACM TOMS 3 (1977), 326–343.
CAPTCHA Image