[论文解读] Numerical integration of dynamical systems with Lie series: Relativistic acceleration and non-gravitational forces

David Bancelin, Daniel Hestroffer+1 more|arXiv (Cornell University)|Aug 29, 2016

Astro and Planetary Science参考文献 17被引用 8

一句话总结

该论文将李级数数值积分方法扩展至包含相对论加速度和雅科夫斯基效应的动力系统，实现了包含非引力作用力的N体问题的高精度、辛积分。该方法在水星近日点进动计算中精度达到1.1×10⁻⁴弧秒/轨道，且在雅科夫斯基效应引起的半长轴漂移方面与Radau方法和解析结果具有良好的一致性（差异约5%）。

ABSTRACT

The integration of the equations of motion in gravitational dynamical systems -- either in our Solar System or for extra-solar planetary system -- being non integrable in the global case, is usually performed by means of numerical integration. Among the different numerical techniques available for solving ordinary differential equations, the numerical integration using Lie series has shown some advantages. In its original form (Hanslmeier 1984), it was limited to the N-body problem where only gravitational interactions are taken into account. We present in this paper a generalisation of the method by deriving an expression of the Lie-terms when other major forces are considered. As a matter of fact, previous studies had been made but only for objects moving under gravitational attraction. If other perturbations are added, the Lie integrator has to be re-built. In the present work we consider two cases involving position and position-velocity dependent perturbations: relativistic acceleration in the framework of General Relativity and a simplified force for the Yarkovsky effect. A general iteration procedure is applied to derive the Lie series to any order and precision. We then give an application to the integration of the equation of motions for typical Near-Earth objects and planet Mercury.

研究动机与目标

• 将李级数积分方法从纯引力N体问题扩展至包含非引力作用力（如相对论加速度和雅科夫斯基效应）的系统。
• 推导适用于位置与速度相关摄动的李项的一般递推公式。
• 通过与解析解及既有的数值方法（如Radau方法）对比，验证扩展后的积分器在关键动力系统中的性能。
• 实现对近地小行星（NEOs）和水星等行星的长期、高精度模拟，包含相对论与非引力摄动。
• 通过将积分器推广至任意位置与速度相关的作用力，为未来在撞击概率评估、稳定性研究及彗星动力学中的应用奠定基础。

提出的方法

• 推导出适用于任意位置与速度相关作用力的李级数项的广义递推公式，扩展了Hanslmeier与Dvorak（1984）原始框架。
• 应用李算子形式化方法，通过李级数展开对运动方程进行积分，其中解表示为李导数算子作用于初始条件的指数形式。
• 通过后牛顿近似引入相对论加速度，近日点进动率基于标准相对论公式推导。
• 将雅科夫斯基效应建模为与位置和速度相关的力，采用简化的横向分量模型，与轨道根数的高斯方程一致。
• 实现一种通用的迭代过程，以任意所需阶数和精度计算李级数项，确保长期积分的高精度。
• 通过将李积分器结果与解析解（如水星近日点进动）及Radau积分器在雅科夫斯基效应引起的半长轴漂移方面的结果进行对比，验证该方法。

实验结果

研究问题

• RQ1 李级数积分方法能否被推广以包含非引力作用力，如相对论加速度和雅科夫斯基效应？
• RQ2 扩展后的李积分器在重现已知的相对论进动率（如水星近日点进动）方面精度如何？
• RQ3 李积分器在重现雅科夫斯基效应引起的半长轴长期变化方面，与现有积分器和解析模型相比表现如何？
• RQ4 在模拟近地小行星等小天体的长期轨道演化时，扩展后的李积分器能否保持高精度与高效性？
• RQ5 将李积分器推广至处理天体力学中任意位置与速度相关作用力的可行性如何？

主要发现

• 扩展后的李积分器对水星相对论近日点进动的误差仅为1.1×10⁻⁴弧秒/轨道，与解析值0.103526弧秒/轨道高度一致。
• 在雅科夫斯基效应方面，李积分器计算的半长轴漂移与Radau积分器结果的差异在0.1–0.3公里/轨道之间（即约0.3–0.5%），与Chesley等（2008）的结果差异在~4–5%以内。
• 该方法可同时处理无质量天体与有质量天体，其能力与Radau积分器相当。
• 该算法可实现对近地小行星（包括受相对论与雅科夫斯基效应影响者）的精确长期积分，这对撞击概率与近距离飞掠研究至关重要。
• 李项的递推公式已推广至包含位置与速度相关力，为未来扩展至其他摄动类型奠定基础。
• 实现结果表明，即使在引入非引力作用力后，李级数方法仍保持高精度与辛性质，支持其在稳定性与动力学演化研究中的应用。