EnckeHH: an integrator for gravitational dynamics with a dominant mass that achieves optimal error behaviour

Abstract

We present EnckeHH, a new, highly accurate code for orbital dynamics of perturbed Keplerian systems such as planetary systems or galactic centre systems. It solves Encke's equations of motion, which assume perturbed Keplerian orbits. By incorporating numerical techniques, we have made the code follow optimal roundoff error growth. In a $10^{12}$ day integration of the outer Solar System, EnckeHH was $3.5$ orders of magnitude more accurate than IAS15 in a fixed time step test. Adaptive steps are recommended for IAS15. Through study of efficiency plots, we show that EnckeHH reaches significantly higher accuracy than the Rebound integrators IAS15 and WHCKL for fixed step size.

arXiv comment8 pages, 5 figures, 1 table. Matches accepted MNRAS version. Instructions on code usage available at https://github.com/matthewholman/rebound/tree/master/examples/encke_hh

References in corpus (15)

• IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits
• WHFast: A fast and unbiased implementation of a symplectic Wisdom-Holman integrator for long term gravitational simulations
• REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
• Hybrid Symplectic Integrators for Planetary Dynamics
• High order symplectic integrators for planetary dynamics and their implementation in REBOUND
• Is the outer Solar System chaotic?
• A Fast and Accurate Universal Kepler Solver without Stumpff Series
• A study of symplectic integrators for planetary system problems: error analysis and comparisons
• Slimplectic Integrators: Variational Integrators for General Nonconservative Systems
• Time-symmetric integration in astrophysics
• Are long-term $N$-body simulations reliable?
• The use of Kepler solver in numerical integrations of quasi-Keplerian orbits
• Should $N$-body integrators be symplectic everywhere in phase space?
• On the accuracy of symplectic integrators for secularly evolving planetary systems
• Improving the accuracy of simulated chaotic $N$-body orbits using smoothness

Cited by in corpus (5)

• FROST: a momentum-conserving CUDA implementation of a hierarchical fourth-order forward symplectic integrator
• Revisiting high-order Taylor methods for astrodynamics and celestial mechanics
• Switching integrators reversibly in the astrophysical $N$-body problem
• Terrestrial Exoplanet Simulator (TES): an error optimal planetary systems integrator that permits close encounters
• Significant improvement in planetary system simulations from statistical averaging