posted on 2017-01-04, 15:32authored byWalter Dehnen, D. M. Hernandez
We study analytically and experimentally certain symplectic and time-reversible N-body integrators which employ a Kepler solver for each pair-wise interaction, including the method of Hernandez & Bertschinger (2015). Owing to the Kepler solver, these methods treat close two-body interactions correctly, while close three-body encounters contribute to the truncation error at second order and above. The second-order errors can be corrected to obtain a fourth-order scheme with little computational overhead. We generalise this map to an integrator which employs a Kepler solver only for selected interactions and yet retains fourth-order accuracy without backward steps. In this case, however, two-body encounters not treated via a Kepler solver contribute to the truncation error.
Funding
This work used the DiRAC Complexity system, operated by the
University of Leicester IT Services, which forms part of the STFC
DiRAC HPC Facility (www.dirac.ac.uk). This equipment is funded
by BIS National E-Infrastructure capital grant ST/K000373/1 and
STFC DiRAC Operations grant ST/K0003259/1. DiRAC is part of
the National E-Infrastructure.
History
Citation
Monthly Notices of the Royal Astronomical Society (February 11, 2017) 465 (1): 1201-1217.
Author affiliation
/Organisation/COLLEGE OF SCIENCE AND ENGINEERING/Department of Physics and Astronomy
Version
VoR (Version of Record)
Published in
Monthly Notices of the Royal Astronomical Society (February 11
Publisher
Oxford University Press (OUP), Royal Astronomical Society