Hostname: page-component-586b7cd67f-dlnhk Total loading time: 0 Render date: 2024-11-23T14:53:10.636Z Has data issue: false hasContentIssue false
  • English
  • Français

A Priori Estimates for the Global Error Committed by Runge-Kutta Methods for a Nonlinear Oscillator

Published online by Cambridge University Press:  01 February 2010

Jitse Niesen
Jitse Niesen
Affiliation:
Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road, Cambridge CB3 OWA J.Niesen@damtp.cam.ac.uk, http://www.damtp.cam.ac.uk/user/na/people/Jitse/index.html

Abstract

Core share and HTML view are not available for this content. However, as you have access to this content, a full PDF is available via the ‘Save PDF’ action button.

The Alekseev–Gröbner lemma is combined with the theory of modified equations to obtain an a priori estimate for the global error of numerical integrators. This estimate is correct up to a remainder term of order h2p, where h denotes the step size and p the order of the method. It is applied to nonlinear oscillators whose behaviour is described by the Emden–Fowler equation y″+tνyn=0. The result shows explicitly that later terms sometimes blow up faster than the leading term of order hp, necessitating the whole computation. This is supported by numerical experiments.

Type
Research Article
Information
LMS Journal of Computation and Mathematics , Volume 6 , 2003 , pp. 18 - 28
Copyright
Copyright © London Mathematical Society 2003

References

1Benettin, G. and Giorgilli, A., ‘On the Hamiltonian interpolation of near-to-the identity symplectic mappings with application to symplectic integration algorithms’, J. Statist. Phys. 74 (1994) 1117–1143.CrossRefGoogle Scholar
2Calvo, M.P. and Hairer, E., ‘Accurate long-term integration of dynamical systems’, Appl. Numer. Math. 18 (1995) 95–105.CrossRefGoogle Scholar
3Cano, B. and Sanz-Serna, J.M., ‘Error growth in the numerical integration of periodic orbits by multistep methods, with application to reversible systems’, IMA J. Numer. Anal. 18(1998) 57–75.CrossRefGoogle Scholar
4Chandrasekhar, S., An introduction to the study of stellar structure(University of Chicago Press, 1939).Google Scholar
5Gragg, W., ‘Repeated extrapolation to the limit in the numerical solution of ordinary differential equations’, Ph.D. thesis, University of California, Los Angeles, 1964.Google Scholar
6Hairer, E. and Lubich, Ch., ‘Asymptotic expansions of the global error of fixed-stepsize methods’, Numer. Math. 45 (1984) 345–360.CrossRefGoogle Scholar
7Hairer, E. and Lubich, Ch., ‘Asymptotic expansions and backward analysis for numerical integrators’, Dynamics of algorithms {Minneapolis, MN, 1997), IMA Vol. Math. Appl. 118 (ed.de la Llave, R., Petzold, L.R. and Lorenz, J., Springer, New York, 2000) 91–106.CrossRefGoogle Scholar
8Hairer, E., Lubich, Ch. and Wanner, G., Geometric numerical integration. Structure-preserving algorithms for ordinary differential equations, Springer Ser. Comput.Math. 31 (Springer, Berlin, 2002).Google Scholar
9Hairer, E., Nørsett, S. and Wanner, G., Solving ordinary differential equations I. Nonstiff problems, 2nd edn (Springer, Berlin, 1993).Google Scholar
10Iserles, A., ‘On the global error of discretization methods for highly-oscillatory ordinary differential equations’, BIT 42(2002) 561–599.CrossRefGoogle Scholar
11Lambert, J., Numerical methods for ordinary differential equations: the initial value problem (John Wiley & Sons, Chichester, 1991).Google Scholar
12Neville, E., Jacobian elliptic functions (Clarendon Press, Oxford, 1944).Google Scholar
13Wong, J., ‘On the generalized Emden-Fowler equation’, SIAM Rev. 17(1975) 339–360.CrossRefGoogle Scholar