Energy-preserving numerical schemes of high accuracy for one-dimensional Hamiltonian systems
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 fourth order. We study the error accumulation for a very long time and the conservation of the energy integral.
- Publication:
-
Journal of Physics A Mathematical General
- Pub Date:
- April 2011
- DOI:
- 10.1088/1751-8113/44/15/155206
- arXiv:
- arXiv:1009.2738
- Bibcode:
- 2011JPhA...44o5206C
- Keywords:
-
- Computer Science - Numerical Analysis;
- Mathematics - Numerical Analysis;
- Physics - Computational Physics;
- 65P10;
- 65L12
- E-Print:
- 15 pages, 6 figures. Presented at the conference "BIT 50 - Trends in Numerical Computing" (Lund, 17-20 June 2010)