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 1012 d integration of the outer Solar system, ENCKEHH was 3.5 orders of magnitude more accurate than IAS15 in a fixed timestep 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.
- Publication:
-
Monthly Notices of the Royal Astronomical Society
- Pub Date:
- March 2021
- DOI:
- 10.1093/mnras/staa3945
- arXiv:
- arXiv:2010.13907
- Bibcode:
- 2021MNRAS.502..556H
- Keywords:
-
- methods: numerical;
- celestial mechanics;
- planets and satellites: dynamical evolution and stability;
- galaxies: evolution;
- Astrophysics - Earth and Planetary Astrophysics;
- Astrophysics - Astrophysics of Galaxies;
- Astrophysics - Instrumentation and Methods for Astrophysics
- E-Print:
- 8 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