A Deferred Correction Method for Nonlinear Two-Point Boundary Value Problems: Implementation and Numerical Evaluation

A deferred correction method for the numerical solution of nonlinear two-point boundary value problems has been derived and analyzed in two recent papers by the first author. The method is based on mono-implicit Runge–Kutta formulas and is specially designed to deal efficiently with problems whose solutions contain nonsmooth parts—in particular, singular perturbation problems of boundary layer or turning point type. This paper briefly describes an implementation of the method and gives the results of extensive numerical testing on a set of nonlinear problems that includes both smooth and increasingly stiff (and difficult) problems. Results on the test set are also given using the available codes COLSYS and COLNEW. Although the intent is not to make a formal comparison, the code described appears to be competitive in speed and storage requirements on these problems.

[1]  R. O'Malley.  Topics in singular perturbations , 1968 .

[2]  Victor Pereyna.  Iterated deferred corrections for nonlinear boundary value problems , 1968 .

[3]  L. Segel,et al.  Introduction to Singular Perturbations. By R. E. O'MALLEY, JR. Academic Press, 1974. $ 16.50. , 1975, Journal of Fluid Mechanics.

[4]  V. Pereyra,et al.  A variable order finite difference method for nonlinear multipoint boundary value problems , 1974 .

[5]  H. J. Pesch,et al.  Comparing routines for the numerical solution of initial value problems of ordinary differential equations in multiple shooting , 1976 .

[6]  H. A. Watts,et al.  Computational Solution of Linear Two-Point Boundary Value Problems via Orthonormalization , 1977 .

[7]  Milan Kubicek,et al.  Test Examples for Comparison of Codes for Nonlinear Boundary - Value Problems in Ordinary Differential Equations , 1978, Codes for Boundary-Value Problems in Ordinary Differential Equations.

[8]  U. Ascher,et al.  A collocation solver for mixed order systems of boundary value problems , 1979 .

[9]  Robert D. Russell,et al.  Collocation Software for Boundary-Value ODEs , 1981, TOMS.

[10]  Philip E. Gill,et al.  Practical optimization , 1981 .

[11]  J. Cash,et al.  High order methods for the numerical solution of two-point boundary value problems , 1982 .

[12]  J. Cash,et al.  Mono-implicit Runge—Kutta Formulae for the Numerical Integration of Stiff Differential Systems , 1982 .

[13]  U. Ascher,et al.  On Spline Basis Selection for Solving Differential Equations , 1983 .

[14]  M. R. Maier,et al.  Numerical solution of singular perturbed boundary value problems using a collocation method with tension splines , 1985 .

[15]  S. Pruess.  Interpolation schemes for collocation solutions of two point boundary value problems , 1986 .

[16]  J. Cash.  Numerical integration of non-linear two-point boundary-value problems using iterated deferred corrections—I , 1986 .

[17]  U. Ascher,et al.  A New Basis Implementation for a Mixed Order Boundary Value ODE Solver , 1987 .

[18]  J. Cash.  On the Numerical Integration of Nonlinear Two-Point Boundary Value Problems Using Iterated Deferred Corrections. Part 2: The Development and Analysis of Highly Stable Deferred Correction Formulae , 1988 .

[19]  J. Cash.  A comparison of some global methods for solving two-point boundary value problems , 1989 .