Energy-preserving numerical schemes of high accuracy for one-dimensional Hamiltonian systems
arXiv:1009.2738 · doi:10.1088/1751-8113/44/15/155206
Abstract
We present a class of non-standard numerical schemes which are modifications of the discrete gradient method. They preserve the energy integral exactly (up to the round-off error). The considered class contains locally exact discrete gradient schemes and integrators of arbitrary high order. In numerical experiments we compare our integrators with some other numerical schemes, including the standard discrete gradient method, the leap-frog scheme and a symplectic scheme of 4th order. We study the error accumulation for very long time and the conservation of the energy integral.
15 pages, 6 figures. Presented at the conference "BIT 50 - Trends in Numerical Computing" (Lund, 17-20 June 2010)
References in corpus (5)
- A unified approach to Hamiltonian systems, Poisson systems, gradient systems, and systems with Lyapunov functions and/or first integrals
- How to improve the accuracy of the discrete gradient method in the one-dimensional case
- An orbit-preserving discretization of the classical Kepler problem
- Long-time behaviour of discretizations of the simple pendulum equation
- Comment on `conservative discretizations of the Kepler motion'