Revisiting high-order Taylor methods for astrodynamics and celestial mechanics

Abstract

We present heyoka, a new, modern and general-purpose implementation of Taylor's integration method for the numerical solution of ordinary differential equations. Detailed numerical tests focused on difficult high-precision gravitational problems in astrodynamics and celestial mechanics show how our general-purpose integrator is competitive with and often superior to state-of-the-art specialised symplectic and non-symplectic integrators in both speed and accuracy. In particular, we show how Taylor methods are capable of satisfying Brouwer's law for the conservation of energy in long-term integrations of planetary systems over billions of dynamical timescales. We also show how close encounters are modelled accurately during simulations of the formation of the Kirkwood gaps and of Apophis' 2029 close encounter with the Earth (where heyoka surpasses the speed and accuracy of domain-specific methods). heyoka can be used from both C++ and Python, and it is publicly available as an open-source project.

References in corpus (18)

• Seven temperate terrestrial planets around the nearby ultracool dwarf star TRAPPIST-1
• REBOUND: An open-source multi-purpose N-body code for collisional dynamics
• La2010: A new orbital solution for the long term motion of the Earth
• IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits
• Extrasolar planet population synthesis I: Method, formation tracks and mass-distance distribution
• WHFast: A fast and unbiased implementation of a symplectic Wisdom-Holman integrator for long term gravitational simulations
• The nature of the TRAPPIST-1 exoplanets
• High order symplectic integrators for perturbed Hamiltonian systems
• Odeint - Solving ordinary differential equations in C++
• Yarkovsky-driven impact risk analysis for asteroid (99942) Apophis
• Convergent Migration Renders TRAPPIST-1 Long-lived
• New families of symplectic splitting methods for numerical integration in dynamical astronomy
• Enhanced interplanetary panspermia in the TRAPPIST-1 system
• Planetesimal-driven planet migration in the presence of a gas disk
• On the minimal accuracy required for simulating self-gravitating systems by means of direct N-body methods
• High order symplectic integrators for planetary dynamics and their implementation in REBOUND
• SLEEF: A Portable Vectorized Library of C Standard Mathematical Functions
• EnckeHH: an integrator for gravitational dynamics with a dominant mass that achieves optimal error behaviour

Cited by in corpus (6)

• Neural representation of a time optimal, constant acceleration rendezvous
• Nonlinear Propagation of Non-Gaussian Uncertainties
• Reliable event detection for Taylor methods in astrodynamics
• Dark Matter reconstruction from stellar orbits in the Galactic Centre
• The Fellowship of the Dyson Ring: ACT&Friends' Results and Methods for GTOC 11
• An efficient approach for searching three-body periodic orbits passing through Eulerian configuration