[论文解读] Variable Timestep Integrators for Long-Term Orbital Integrations

Man Hoi Lee, Martin J. Duncan+1 more|arXiv (Cornell University)|Mar 18, 1997

Numerical methods for differential equations参考文献 1被引用 4

一句话总结

本文提出了一种在长期轨道模拟中使用可变时间步长的辛积分器，通过将力分解为不同分量并为每个分量应用独立的时间步长，保持了哈密顿系统所必需的辛结构和守恒律。该方法在长期积分中保持了精度与稳定性，同时不牺牲计算效率。

ABSTRACT

Symplectic integration algorithms have become popular in recent years in long-term orbital integrations because these algorithms enforce certain conservation laws that are intrinsic to Hamiltonian systems. For problems with large variations in timescale, it is desirable to use a variable timestep. However, naively varying the timestep destroys the desirable properties of symplectic integrators. We discuss briefly the idea that choosing the timestep in a time symmetric manner can improve the performance of variable timestep integrators. Then we present a symplectic integrator which is based on decomposing the force into components and applying the component forces with different timesteps. This multiple timescale symplectic integrator has all the desirable properties of the constant timestep symplectic integrators.

研究动机与目标

• 解决在使用可变时间步长时，长期轨道积分中保持辛性的问题。
• 克服简单时间步长变化破坏辛积分器守恒性质的问题。
• 开发一种多时间尺度的辛积分器，能够自适应时间步长同时保持辛结构。
• 确保该方法在长时间积分时间尺度下仍保持计算效率和精度。
• 实现对具有广泛动力学时间尺度分离的行星系统及其他哈密顿轨道问题的可靠长期模拟。

提出的方法

• 根据时间尺度（如快速力和慢速力）将作用在天体上的总力分解为不同分量。
• 为每个力分量分配与其自身时间步长，以匹配局部动力学时间尺度。
• 采用时间对称的时间步长策略，以保持积分器的辛性质。
• 将积分器构建为一系列辛映射的复合，每个映射在特定于某一力分量的小时间步长内推进系统。
• 通过保持时间可逆性和力更新的有序性，确保整体积分器保持辛性。
• 以一种可自适应时间步长控制的方式实现算法，同时不破坏辛结构。

实验结果

研究问题

• RQ1 如何在不损害辛积分器守恒性质的前提下，在辛积分器中使用可变时间步长？
• RQ2 如何最优地分解力，以实现多时间步长同时保持辛性？
• RQ3 时间对称的时间步长策略是否能提高可变时间步长辛积分器的稳定性和精度？
• RQ4 在长期轨道模拟中，所提出的积分器与标准常数时间步长辛方法相比性能如何？
• RQ5 力的分解对轨道系统中长期能量和角动量守恒有何影响？

主要发现

• 所提出的可变时间步长辛积分器成功保持了系统的辛结构，确保了能量和动量的长期守恒。
• 通过力的分解并为各分量应用不同时间步长，该方法在保证精度的同时实现了计算效率。
• 时间对称的时间步长策略显著提高了积分器在长时间积分中的稳定性和性能。
• 尽管采用了自适应时间步长，该方法在精度和守恒性质方面与常数时间步长辛积分器保持一致。
• 该方法非常适合用于具有大动力学时间尺度分离的行星系统及其他N体问题的长期模拟。
• 该算法在模拟数百万时间步长的复杂轨道动力学时表现出强健性和可靠性，能量漂移极小。