期刊
MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY
卷 491, 期 2, 页码 2885-2901出版社
OXFORD UNIV PRESS
DOI: 10.1093/mnras/stz2870
关键词
gravitation; methods: numerical; planets and satellites: dynamical evolution and stability
资金
- National Aeronautics and Space Administration (NASA) through the NASA Hubble Fellowship - Space Telescope Science Institute [HST-HF2-51423.001-A]
- Natural Sciences and Engineering Research Council (NSERC) [RGPIN-2014-04553]
- Centre for Planetary Sciences at the University of Toronto Scarborough
- NASA [NAS5-26555]
Symplectic methods, in particular the Wisdom-Holman map, have revolutionized our ability to model the long-term, conservative dynamics of planetary systems. However, many astrophysically important effects are dissipative. The consequences of incorporating such forces into otherwise symplectic schemes are not always clear. We show that moving to a general framework of non-commutative operators (dissipative or not) clarifies many of these questions, and that several important properties of symplectic schemes carry over to the general case. In particular, we show that explicit splitting schemes generically exploit symmetries in the applied external forces, which often strongly suppress integration errors. Furthermore, we demonstrate that so-called 'symplectic correctors' (which reduce energy errors by orders of magnitude at fixed computational cost) apply equally well to weakly dissipative systems and can thus be more generally thought of as 'weak splitting correctors'. Finally, we show that previously advocated approaches of incorporating additional forces into symplectic methods work well for dissipative forces, but give qualitatively wrong answers for conservative but velocity-dependent forces like post-Newtonian corrections. We release REBOUNDx, an open-source C library for incorporating additional effects into REBOUND N-body integrations, together with a convenient python wrapper. All effects are machine independent and we provide a binary format that interfaces with the SimulationArchive class in REBOUND to enable the sharing and reproducibility of results. Users can add effects from a list of pre-implemented astrophysical forces, or contribute new ones.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据