Runge–Kutta Defect Control Using Hermite–Birkhoff Interpolation
algorithmdefect controlHermite-Birkhoff interpolationnumerical testsRunge-Kutta methodstepsize selection
Nonlinear ordinary differential equations and systems (34A34) Numerical methods for initial value problems involving ordinary differential equations (65L05) Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations (65L06) Mesh generation, refinement, and adaptive methods for ordinary differential equations (65L50)
In a single step of a Runge-Kutta method applied to an initial value problem \(y'(x)=f(x,y(x))\), which takes the solution approximation from \(x_ n\) to \(x_{n+1}=x_ n+h\), the defect is defined as the value of \(P'(x)-f(x,P(x))\), where \(P\) is a suitable interpolant to the numerical solution computed within the step. Controlling the defect by controlling the value of \(h\) is a possible means of stepsize selection. What is required is an appropriate choice of \(P\). Amongst the requirements laid down for \(P\) are its global continuity (continuity at least \(C^ 1\) is suggested) as well as accuracy and efficiency. The current paper proposes the use of the Hermite-Birkhoff interpolation scheme as a means of achieving these goals. For an \(s\)-stage method of order \(p\), in addition to the usual stage values, the solution is approximated at values \(\gamma_ 3,\gamma_ 4,\dots,\gamma_ r\) of the way through the step. From the derivatives computed at \(\gamma_ 1=0\) (the initial point of the step) and \(\gamma_ 2=1\) (the end of the step), and at the \(r-2\) additional values, an interpolation scheme of the form \(H(x_ n+\tau h)=d_ 1(\tau)y(x_ n)+d_ 2(\tau)y(x_{n+1})+h\sum^ r_{i=1} e_ i(\tau)y'(x_ n+h\gamma_ i)\) is proposed. The paper discusses a number of practical details concerning the selection of \(\gamma_ 3,\gamma_ 4,\dots,\gamma_ r\). Even for such a low order as \(p=5\), where \(r=5\) is the appropriate choice, the relative merits of various options turn out to be exceedingly complicated. Using values of \(\gamma_ 3,\gamma_ 4,\dots,\gamma_ r\) which are in some sense optimal, and a standard embedded Runge-Kutta method [cf. \textit{J. R. Dormand} and \textit{P. J. Prince}, J. Comput. Appl. Math. 6, 19-26 (1980; Zbl 0448.65045)], a specific algorithm is developed. This compares favourably on numerical tests with a known method based on Hermite interpolation at the points \(x_ n\), \(x_{n+1/2}\) and \(x_{n+1}\).
- Parallel implementation of interpolants for Runge-Kutta pairs
- Solving ODEs and DDEs with residual control
- Continuous numerical methods for ODEs with defect control
- An error term and uniqueness for Hermite-Birkhoff interpolation involving only function values and/or first derivative values
- On the minimization of the defect in continuous Runge-Kutta methods
- Birkhoff Polynomial Basis
- Robust Defect Control with Runge–Kutta Schemes
- A Posteriori Analysis of Fully Discrete Method of Lines Discontinuous Galerkin Schemes for Systems of Conservation Laws
- Runge-Kutta research at Toronto
- Error estimation and control for ODEs
- Polynomial algebra for Birkhoff interpolants
- New interpolants for asymptotically correct defect control of BVODEs
- The reliability/cost trade-off for a class of ODE solvers
This page was built for publication: Runge–Kutta Defect Control Using Hermite–Birkhoff Interpolation
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q3977641)
