期刊
MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY
卷 485, 期 4, 页码 5490-5497出版社
OXFORD UNIV PRESS
DOI: 10.1093/mnras/stz769
关键词
gravitation; methods: numerical; planets and satellites: dynamical evolution and stability
资金
- NSERC [RGPIN-2014-04553]
- Centre for Planetary Sciences at the University of Toronto Scarborough
- NASA through the NASA Hubble Fellowship - Space Telescope Science Institute [HST-HF2-51423.001-A]
- NASA [NAS5-26555]
- 2018 Fields Undergraduate Summer Research Program at the Fields Institute for Research in Mathematical Sciences in Toronto
Hybrid symplectic integrators such as MERCURY are widely used to simulate complex dynamical phenomena in planetary dynamics that could otherwise not be investigated. A hybrid integrator achieves high accuracy during close encounters by using a high-order integration scheme for the duration of the encounter while otherwise using a standard second-order Wisdom-Holman scheme, thereby optimizing both speed and accuracy. In this paper we reassess the criteria for choosing the switching function that determines which parts of the Hamiltonian are integrated with the high-order integrator. We show that the original motivation for choosing a polynomial switching function in MERCURY is not correct. We explain the nevertheless excellent performance of the MERCURY integrator and then explore a wide range of different switching functions including an infinitely differentiable function and a Heaviside function. We find that using a Heaviside function leads to a significantly simpler scheme compared to MERCURY, while maintaining the same accuracy in short-term simulations.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据