Revisiting high-order Taylor methods for astrodynamics and celestial mechanics

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...

Full description

Saved in:
Bibliographic Details
Main Authors: Biscani, Francesco (Author) , Izzo, Dario (Author)
Format: Article (Journal)
Language:English
Published: 15 April 2021
In: Monthly notices of the Royal Astronomical Society
Year: 2021, Volume: 504, Issue: 2, Pages: 2614-2628
ISSN:1365-2966
DOI:10.1093/mnras/stab1032
Online Access:Verlag, lizenzpflichtig, Volltext: https://doi.org/10.1093/mnras/stab1032
Verlag, lizenzpflichtig, Volltext: https://gateway.webofknowledge.com/gateway/Gateway.cgi?GWVersion=2&SrcAuth=DynamicDOIArticle&SrcApp=WOS&KeyAID=10.1093%2Fmnras%2Fstab1032&DestApp=DOI&SrcAppSID=E1Ds9d5qwGMRH1DGxTT&SrcJTitle=MONTHLY+NOTICES+OF+THE+ROYAL+ASTRONOMICAL+SOCIETY&DestDOIRegistrantName=Oxford+University+Press
Get full text
Author Notes:Francesco Biscani and Dario Izzo
Description
Summary: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 specialized 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 time-scales. 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.
Item Description:Gesehen am 10.08.2021
Physical Description:Online Resource
ISSN:1365-2966
DOI:10.1093/mnras/stab1032