Stability with large step sizes for multistep discretizations of stiff ordinary differential equations

Author:
George Majda

Journal:
Math. Comp. **47** (1986), 473-502, S41

MSC:
Primary 65L20

DOI:
https://doi.org/10.1090/S0025-5718-1986-0856698-8

MathSciNet review:
856698

Full-text PDF

Abstract | References | Similar Articles | Additional Information

Abstract: In this paper we consider a large set of variable coefficient linear systems of ordinary differential equations which possess two different time scales, a slow one and a fast one. A small parameter characterizes the stiffness of these systems. We approximate a system of ODE's in this set by a general class of multistep discretizations which includes both one-leg and linear multistep methods. We determine sufficient conditions under which each solution of a multistep method is uniformly bounded, with a bound which is independent of the stiffness of the system of ODE's, when the step size resolves the slow time scale but not the fast one. We call this property stability with large step sizes.

The theory presented in this paper lets us compare properties of one-leg methods and linear multistep methods when they approximate variable coefficient systems of stiff ODE's. In particular, we show that one-leg methods have better stability properties with large step sizes than their linear multistep counterparts. This observation is consistent with results obtained by Dahlquist and Lindberg [11], Nevanlinna and Liniger [33] and van Veldhuizen [42]. Our theory also allows us to relate the concept of *D*-stability (van Veldhuizen [42]) to the usual notions of stability and stability domains and to the propagation of errors for multistep methods which use large step sizes.

**[1]**L. R. Abrahamsson, H. B. Keller & H.-O. Kreiss, "Difference approximations for singular perturbations of systems of ordinary differential equations,"*Numer Math.*, v. 22, 1974, pp. 367-391. MR**0388784 (52:9618)****[2]**G. Bjurel, G. Dahlquist, B. Lindberg, S. Linde & L. Oden,*Survey of Stiff Ordinary Differential Equations*, Report NA 70.11, Dept. of Information Processing, Royal Institute of Technology, S-100 44 Stockholm 70, Sweden, 1970.**[3]**K. Burrage, "High order algebraically stable Runge-Kutta methods,"*BIT*, v. 18, 1978, pp. 373-383. MR**520749 (80e:65073)****[4]**K. Burrage & J. C. Butcher, "Stability criteria for implicit Runge-Kutta methods,"*SIAM J. Numer. Anal.*, v. 16, 1979, pp. 46-57. MR**518683 (80b:65096)****[5]**J. C. Butcher, "A stability property of implicit Runge-Kutta methods,"*BIT*, v. 15, 1975, pp. 358-361.**[6]**M. Crouzeix, "Sur la*B*-stabilité des méthodes de Runge-Kutta,"*Numer. Math.*, v. 32, 1979, pp. 75-82. MR**525638 (80f:65081)****[7]**G. G. Dahlquist, "A special stability problem for linear multistep methods,"*BIT*, v. 3, 1963, pp. 27-43. MR**0170477 (30:715)****[8]**G. G. Dahlquist,*On Stability and Error Analysis for Stiff Non-linear Problems*, Part 1, Report TRITA-NA-7508, Dept. of Information Processing, Royal Institute of Technology, Stockholm, Sweden, 1975.**[9]**G. G. Dahlquist, "The sets of smooth solutions of differential and difference equations," in*Stiff Differential Systems*(R. A. Willoughby, ed.), Plenum Press, New York, 1974, pp. 67-80. MR**0405860 (53:9652)****[10]**G. G. Dahlquist,*Some Properties of Linear Multistep and One-leg Methods for Ordinary Differential Equations*, Report TRITA-NA-7904, Dept. of Numerical Analysis and Computing Science, Royal Institute of Technology, Stockholm, Sweden, 1979.**[11]**G. G. Dahlquist & B. Lindberg,*On Some Implicit One-step Methods for Stiff Differential Equations*, Report TRITA-NA-7302, Dept. of Information Processing, Royal Institute of Technology, Stockholm, Sweden, 1973.**[12]**G. G. Dahlquist & G. Söderland,*Error Propagation in Stiff Differential Systems of Singular Perturbation Type*, Report TRITA-NA-8108, Dept. of Information Processing, Royal Institute of Technology, Stockholm, Sweden, 1981.**[13]**B. L. Ehle,*On Padé Approximations to the Exponential Function and A-stable Methods for the Numerical Solution of Initial Value Problems*, Dept. of Applied Analysis and Computer Science, University of Waterloo, Research Report No. CSRR 2010, 1969.**[14]**W. H. Enright & T. E. Hull, "Comparing numerical methods for the solution of stiff systems of ODE's arising in chemistry," in*Numerical Methods for Differential Systems*(L. Lapidus and W. E. Schiesser, eds.), Academic Press, New York, 1976, pp. 45-66.**[15]**W. H. Enright, T. E. Hull & B. Lindberg, "Comparing numerical methods for stiff systems of ODE's,"*BIT*, v. 15, 1975, pp. 10-48.**[16]**S. F. Feschchenko, N. I. Shkil' & L. D. Nikolenko,*Asymptotic Methods in the Theory of Linear Differential Equations*, American Elsevier, New York, 1967. MR**0221029 (36:4081)****[17]**R. Frank, J. Schneid & C. W. Ueberhuber, "The concept of*B*-convergence,"*SIAM J. Numer. Anal.*, v. 18, 1981, pp. 753-780. MR**629662 (82h:65054)****[18]**C. W. Gear,*Numerical Initial Value Problems in Ordinary Differential Equations*, Prentice-Hall, Englewood Cliffs, N. J., 1971. MR**0315898 (47:4447)****[19]**A. R. Gourlay, "A note on trapezoidal methods for the solution of initial value problems,"*Math. Comp.*, v. 24, 1970, pp. 629-633. MR**0275680 (43:1433)****[20]**F. Hoppensteadt, "Properties of solutions of ordinary differential equations with small parameters,"*Comm. Pure Appl. Math.*, v. 24, 1971, pp. 807-840. MR**0288378 (44:5576)****[21]**F. Hoppensteadt & W. L. Miranker,*Numerical Solution of Differential Equations with Rapidly Changing Solutions*, IBM Report RC4792, IBM Research Center, Yorktown Heights, New York, 1974.**[22]**H.-O. Kreiss, "Difference methods for stiff ordinary differential equations,"*SIAM J. Numer. Anal.*, v. 15, 1978, pp. 21-58. MR**486570 (80a:65149)****[23]**H.-O. Kreiss, "Problems with different time scales for ordinary differential equations,"*SIAM J. Numer. Anal.*, v. 16, 1979, pp. 980-998. MR**551320 (81a:65087)****[24]**J. D. Lambert,*Computational Methods in Ordinary Differential Equations*, Wiley, London, 1973. MR**0423815 (54:11789)****[25]**L. Lapidus & W. E. Schiesser,*Numerical Methods for Differential Systems*, Academic Press, New York, 1976. MR**0438721 (55:11629)****[26]**B. Lindberg,*IMPEX--A Program Package for the Solution of Systems of Stiff Differential Equations*, Report TRITA-NA-72.50, Royal Institute of Technology, Stockholm, Sweden, 1972.**[27]**B. Lindberg,*IMPEX 2--A Procedure for Solution of Systems of Stiff Differential Equations*, Report TRITA-NA-7303, Royal Institute of Technology, Stockholm, Sweden, 1973.**[28]**G. Majda, "Filtering techniques for systems of stiff ordinary differential equations I,"*SIAM J. Numer. Anal.*, v. 21, 1984, pp. 535-566. MR**744172 (85i:65089)****[29]**G. Majda, "Filtering techniques for systems of stiff ordinary differential equations II. Error estimates,"*SIAM J. Numer. Anal.*, v. 22, 1985, pp. 1116-1134. MR**811187 (87d:65074)****[30]**G. Majda,*A New Theory for Multistep Discretizations of Stiff Ordinary Differential Equations*I:*Stability with Large Step Sizes*, Report NA-3, Division of Applied Mathematics, Brown University, Providence, R.I., 1983.**[31]**W. L. Miranker, "Numerical methods of boundary layer type for stiff systems of differential equations,"*Computing*, v. 11, 1973, pp. 221-234. MR**0386276 (52:7134)****[32]**W. L. Miranker,*Numerical Methods for Stiff Equations and Singular Perturbation Problems*, Reidel, Boston, 1981. MR**603627 (82h:65058)****[33]**O. Nevanlinna & W. Liniger, "Contractive methods for stiff differential equations, part 1,"*BIT*, v. 18, 1978, pp. 457-474. MR**520755 (80h:65053a)****[34]**F. Odeh & W. Liniger, "A note on unconditional fixed-*h*stability of linear multistep formulae,"*Computing*, v. 7, 1971, pp. 240-253. MR**0298957 (45:8006)****[35]**F. Odeh & W. Liniger, "Nonlinear fixed*h*-stability of linear multistep formulas,"*J. Math. Anal. Appl.*, v. 61, 1977, pp. 691-712. MR**0468197 (57:8035)****[36]**A. Prothero & A. Robinson, "On the stability and accuracy of one-step methods for solving stiff systems of ordinary differential equations,"*Math. Comp.*, v. 28, 1974, pp. 145-162. MR**0331793 (48:10125)****[37]**H. J. Stetter,*Towards a Theory for Discretizations of Stiff Differential Systems*, in Lecture Notes in Math., vol. 506, Springer-Verlag, New York, 1975, pp. 190-201. MR**0455414 (56:13652)****[38]**H. J. Stetter,*Analysis of Discretization Methods for Ordinary Differential Equations*, Springer-Verlag, Berlin, 1973. MR**0426438 (54:14381)****[39]**G. W. Strang, "Accurate partial difference methods II. Non-linear problems,"*Numer. Math.*, v. 6, 1964, pp. 37-46. MR**0166942 (29:4215)****[40]**M. van Veldhuizen,*Convergence of One-step Discretizations for Stiff Differential Equations*, Ph.D. Thesis, Mathematical Institute, University of Utrecht, Netherlands, 1973.**[41]**M. van Veldhuizen, "Consistency and stability for one-step discretizations of stiff differential equations," in*Stiff Differential Systems*(R. Willoughby, ed.), Plenum Press, New York, 1974, pp. 259-270. MR**0431694 (55:4689)****[42]**M. van Veldhuizen, "*D*-stability,"*SIAM J. Numer. Anal.*, v. 18, 1981, pp. 45-64. MR**603430 (82e:65089)****[43]**W. Wasow,*Asymptotic Expansions for Ordinary Differential Equations*, Interscience, Wiley, New York, 1965. MR**0203188 (34:3041)****[44]**O. Widlund, "A note on unconditionally stable linear multistep methods,"*BIT*, v. 7, 1967, pp. 65-70. MR**0215533 (35:6373)****[45]**R. A. Willoughby,*Stiff Differential Systems*, Plenum Press, New York, 1974. MR**0343619 (49:8359)**

Retrieve articles in *Mathematics of Computation*
with MSC:
65L20

Retrieve articles in all journals with MSC: 65L20

Additional Information

DOI:
https://doi.org/10.1090/S0025-5718-1986-0856698-8

Article copyright:
© Copyright 1986
American Mathematical Society