A Runge-kutta Type Boundary Value Ode Solver with Defect Control

A popular approach for the numerical solution of boundary value ODE problems involves the use of collocation methods. Such methods can be naturally implemented so as to provide a continuous approximation to the solution over the entire problem interval. On the other hand, several authors have suggested as an alternative, certain subclasses of the implicit Runge-Kutta formulas, known as mono-implicit Runge-Kutta (MIRK) formulas, which can be implemented substantially more eeciently than the collocation methods. These latter formulas do not have a natural implementation that provides a continuous approximation to the solution; rather only a discrete approximation at certain points within the problem interval is obtained. However recent work in the area of initial value problems has demonstrated the possibility of generating inexpensive interpolants for any explicit Runge-Kutta formula. These ideas have recently been extended to develop inter-polants for the MIRK formulas. In this paper, we describe our investigation of the use of continuous MIRK formulas in a code for the numerical solution of boundary value ODE problems. A primary thrust of this investigation is to consider defect control, based on these interpolants, as an alternative to the standard use of global error estimates, as the basis for termination and mesh redistribution criteria.

[1]  J. H. Verner,et al.  Differentiable interpolants for high-order Runge-Kutta methods , 1993 .

[2]  Uri M. Ascher,et al.  Collocation for Singular Perturbation Problems III: Nonlinear Problems without Turning Points , 1982 .

[3]  W. H. Enright A new error-control for initial value solvers , 1989 .

[4]  Marino Zennaro,et al.  Derivation of Efficient, Continuous, Explicit Runge-Kutta Methods , 1992, SIAM J. Sci. Comput..

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

[6]  J. Butcher The numerical analysis of ordinary differential equations: Runge-Kutta and general linear methods , 1987 .

[7]  Uri M. Ascher,et al.  Collocation for Singular Perturbation Problems I: First Order Systems with Constant Coefficients , 1981 .

[8]  W. H. Enright,et al.  Effective solution of discontinuous IVPs using a Runge-Kutta formula pair with interpolants , 1988 .

[9]  N. Nichols,et al.  Numerical methods for stiff two-point boundary value problems , 1986 .

[10]  Margaret H. Wright,et al.  A Deferred Correction Method for Nonlinear Two-Point Boundary Value Problems: Implementation and Numerical Evaluation , 1991, SIAM J. Sci. Comput..

[11]  Richard Weiss,et al.  The application of implicit Runge-Kutta and collection methods to boundary-value problems , 1974 .

[12]  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 .

[13]  Kevin Burrage,et al.  Order results for mono-implicit Runge-Kutta methods , 1994 .

[14]  Uri M. Ascher,et al.  Collocation for singular perturbation problems. II. Linear first order systems without turning points , 1984 .

[15]  J. Lorenz,et al.  A high-order method for stiff boubdary value problems with turning points , 1987 .

[16]  Suchitra Gupta,et al.  An Adaptive Boundary Value Runge–Kutta Solver for First Order Boundary Value Problems , 1985 .

[17]  Paul H. Muir,et al.  Order barriers and characterizations for continuous mono-implicit Runge-Kutta schemes , 1993 .

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

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

[20]  L. S. Baca,et al.  Practical aspects of interpolation in Runge-Kutta codes , 1987 .

[21]  Marino Zennaro,et al.  Order barriers for continuous explicit Runge-Kutta methods , 1991 .

[22]  Wayne H. Enright,et al.  Interpolants for Runge-Kutta formulas , 1986, TOMS.

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

[24]  W. H. Enright Analysis of error control strategies for continuous Runge-Kutta methods , 1989 .