Kepler’s first law states that the orbit of a point mass with negative energy in a classical gravitational potential is an ellipse with one of its foci at the gravitational center. In numerical simulations of this system one often observes a slight precession of the ellipse around the gravitational center. Using the Lagrangian structure of modified equations and a perturbative version of Noether’s theorem, we provide leading order estimates of this precession for the implicit MidPoint rule (MP) and the Störmer-Verlet method (SV). Based on those estimates we construct some new numerical integrators that perform significantly better than MP and SV on the Kepler problem.
CITATION STYLE
Vermeeren, M. (2018). Numerical precession in variational discretizations of the Kepler problem. In Springer Proceedings in Mathematics and Statistics (Vol. 267, pp. 333–348). Springer New York LLC. https://doi.org/10.1007/978-3-030-01397-4_10
Mendeley helps you to discover research relevant for your work.