期刊
MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY
卷 465, 期 1, 页码 1201-1217出版社
OXFORD UNIV PRESS
DOI: 10.1093/mnras/stw2758
关键词
gravitation; methods: analytical; methods: numerical; celestial mechanics; planets and satellites: dynamical evolution and stability; globular clusters: general
资金
- BIS National E-Infrastructure capital grant [ST/K000373/1]
- STFC DiRAC Operations grant [ST/K0003259/1]
- STFC [ST/M006948/1, ST/N000757/1, ST/K000373/1] Funding Source: UKRI
We study analytically and experimentally certain symplectic and time-reversible N-body integrators which employ the Kepler solver for each pairwise interaction, including the method of Hernandez & Bertschinger. 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 generalize this map to an integrator which employs the 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 the Kepler solver contribute to the truncation error.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据