Variable stepsize variable formula methods based on predictor-corrector schemes

Abstract The use of fairly general predictor-corrector (PC) schemes of linear multistep (LM) formulae in the numerical solution of systems of ordinary differential equations (ODE's) is considered. It is assumed that both the stepsize and the PC scheme can be varied during the computational process. The numerical methods obtained under these two assumptions are called predictor-corrector linear multistep variable stepsize variable formula methods (PC LM VSVFM's). The consistency, zero-stability and convergence properties of the PC LM VSVFM's are studied. Several results concerning these fundamental properties of the numerical methods are established. It should be emphasized that all theorems are formulated and proved under very mild assumptions on the stepsize selection strategy. The extension of the results for the so-called one-leg methods is briefly discussed. The use of PC LM VSVFM's leads to a very efficient treatment of many mathematical models describing different phenomena in science and engineering. Such methods have successfully been used in the numerical solution of systems of ODE's arising after the space discretization of some air pollution models.

[1]  G. Dahlquist Convergence and stability in the numerical integration of ordinary differential equations , 1956 .

[2]  P. G. Thomsen,et al.  Two-parameter families of predictor-corrector methods for the solution of ordinary differential equations , 1979 .

[3]  P. G. Thomsen,et al.  Application of backward differentiation methods to the finite element solution of time‐dependent problems , 1979 .

[4]  Zahari Zlatev,et al.  Stability properties of variable stepsize variable formula methods , 1978 .

[5]  Charles William Gear,et al.  Numerical Solution of Ordinary Differential Equations: Is There Anything Left to Do? , 1981 .

[6]  O. Nevanlinna,et al.  Stability of explicit time discretizations for solving initial value problems , 1982 .

[7]  J. Lambert Computational Methods in Ordinary Differential Equations , 1973 .

[8]  Z. Zlatev Zero-stability properties of the three-ordinate variable stepsize variable formula methods , 1981 .

[9]  G. Dahlquist Error analysis for a class of methods for stiff non-linear initial value problems , 1976 .

[10]  Zahari Zlatev,et al.  Automatic Solution of Differential Equations Based on the User of Linear Multistep Methods , 1979, TOMS.

[11]  Zahari Zlatev,et al.  Implementation of a variable stepsize variable formula method in the time-integration part of a code for treatment of long-range transport of air pollutants , 1984 .

[12]  Fred T. Krogh Predictor-Corrector Methods of High Order With Improved Stability Characteristics , 1966, JACM.

[13]  Zahari Zlatev,et al.  Testing subroutines solving advection-diffusion equations in atmospheric environments , 1983 .

[14]  Zahari Zlatev,et al.  Stability Restrictions on Time-Stepsize for Numerical Integration of First-Order Partial Differential Equations* , 1983 .

[15]  L. Shampine,et al.  Computer solution of ordinary differential equations : the initial value problem , 1975 .

[16]  A. Sedgwick,et al.  An effective variable-order variable-step adams method. , 1973 .

[17]  H. Stetter Analysis of Discretization Methods for Ordinary Differential Equations , 1973 .

[18]  A. Hindmarsh LSODE and LSODI, two new initial value ordinary differential equation solvers , 1980, SGNM.

[19]  Fred T. Krogh,et al.  On Testing a Subroutine for the Numerical Integration of Ordinary Differential Equations , 1973, JACM.

[20]  C. W. Gear,et al.  Numerical initial value problem~ in ordinary differential eqttations , 1971 .

[21]  C. W. Gear,et al.  The automatic integration of ordinary differential equations , 1971, Commun. ACM.

[22]  C. William Gear Runge-Kutta Starters for Multistep Methods , 1980, TOMS.

[23]  T. E. Hull,et al.  Comparing Numerical Methods for Ordinary Differential Equations , 1972 .

[24]  P. Henrici Discrete Variable Methods in Ordinary Differential Equations , 1962 .

[25]  G. Dahlquist Stability and error bounds in the numerical integration of ordinary differential equations , 1961 .

[26]  Charles William Gear,et al.  Stability and convergence of variable order multistep methods , 1974 .

[27]  Peter Piotrowski Stability, consistency and convergence of variable K-step methods for numerical integration of large systems of ordinary differential equations , 1969 .

[28]  R. März Variable multistep methods , 1984 .