[论文解读] High order symplectic integrators for perturbed Hamiltonian systems

J. Laskar, Philippe Robutel|arXiv (Cornell University)|May 4, 2000

Numerical methods for differential equations参考文献 15被引用 29

一句话总结

本文提出了一种构造性方法，用于为形如 $ H = A + \varepsilon B $ 的扰动哈密顿系统生成仅使用正时间步长的高阶辛积分器，其中 $ A $ 和 $ B $ 均为可积系统。该文证明了此类积分器的存在性，其余项为 $ O(\tau^p \varepsilon + \tau^2 \varepsilon^2) $，并表明通过引入校正步骤可将余项降低至 $ O(\tau^p \varepsilon + \tau^4 \varepsilon^2) $，显著提升了精度与稳定性，优于标准的跃迁法（leapfrog）方法。

ABSTRACT

We present a class of symplectic integrators adapted for the integration of perturbed Hamiltonian systems of the form $H=A+εB$. We give a constructive proof that for all integer $p$, there exists an integrator with positive steps with a remainder of order $O(τ^pε+τ^2ε^2)$, where $τ$ is the stepsize of the integrator. The analytical expressions of the leading terms of the remainders are given at all orders. In many cases, a corrector step can be performed such that the remainder becomes $O(τ^pε+τ^4ε^2)$. The performances of these integrators are compared for the simple pendulum and the planetary 3-Body problem of Sun-Jupiter-Saturn.

研究动机与目标

• 为扰动哈密顿系统 $ H = A + \varepsilon B $ 构建仅使用正时间步长的高阶辛积分器，克服现有高阶方法因使用负时间步长而引发的稳定性问题。
• 提供一种构造性证明，证明任意阶数 $ p $ 的辛积分器存在，且余项为 $ O(\tau^p \varepsilon + \tau^2 \varepsilon^2) $。
• 利用李代数形式与 Campbell-Baker-Hausdorff 定理，推导出所有阶数下主导余项的解析表达式。
• 提出一种校正步骤，将余项从 $ O(\tau^p \varepsilon + \tau^2 \varepsilon^2) $ 降低至 $ O(\tau^p \varepsilon + \tau^4 \varepsilon^2) $，在不增加步数的前提下提升精度。
• 在典型测试问题（如单摆和太阳-木星-土星三体系统）上评估性能，证明其在精度与稳定性方面显著优于标准跃迁法。

提出的方法

• 该方法采用李代数形式，将时间演化算符表示为 $ \exp(\tau L_H) $，其中 $ H = A + \varepsilon B $，并通过 $ \exp(\tau L_A) $ 和 $ \exp(\tau L_{\varepsilon B}) $ 的复合形式近似该算符。
• 辛积分器被构造为李导数 $ L_A $ 和 $ L_{\varepsilon B} $ 指数的复合，形成对称格式 $ \mathcal{S}\mathcal{A}\mathcal{B}\mathcal{A}_n $ 和 $ \mathcal{S}\mathcal{B}\mathcal{A}\mathcal{B}_n $。
• 通过求解由 Campbell-Baker-Hausdorff 展开导出的代数方程，确定复合中的系数，以确保形式哈密顿量 $ K $ 与 $ H $ 在 $ \tau^p $ 阶内一致。
• 使用自由李代数 $ \mathcal{L}(A,B) $ 的 Lyndon 基，对泊松括号进行无歧义分解，并定义 $ K $ 的形式级数。
• 在积分器输出后应用校正步骤，将余项从 $ O(\tau^2 \varepsilon^2) $ 降低至 $ O(\tau^4 \varepsilon^2) $，从而提升高阶格式的精度。
• 在单摆和三体行星系统上对方法进行数值验证，比较能量误差与计算成本。

实验结果

研究问题

• RQ1 能否为扰动哈密顿系统 $ H = A + \varepsilon B $ 构造出仅使用正时间步长的高阶辛积分器？
• RQ2 此类积分器在任意阶数 $ p $ 下的主导余项的解析形式为何？
• RQ3 能否设计一种校正步骤，将余项从 $ O(\tau^2 \varepsilon^2) $ 降低至 $ O(\tau^4 \varepsilon^2) $？
• RQ4 与常用于行星 N 体问题的跃迁法相比，这些积分器在精度与稳定性方面表现如何？
• RQ5 在何种条件下，低阶积分器的复合能提升性能？何时这种复合方式在成本上不具优势？

主要发现

• 本文构造了一类仅使用正时间步长的辛积分器 $ \mathcal{S}\mathcal{A}\mathcal{B}\mathcal{A}_n $ 和 $ \mathcal{S}\mathcal{B}\mathcal{A}\mathcal{B}_n $，其余项为 $ O(\tau^p \varepsilon + \tau^2 \varepsilon^2) $，适用于任意阶数 $ p $。
• 当 $ n=3 $ 或 $ n=4 $ 时，经校正的积分器 $ \mathcal{S}\mathcal{A}\mathcal{B}\mathcal{A}_{Cn} $ 和 $ \mathcal{S}\mathcal{B}\mathcal{A}\mathcal{B}_{Cn} $ 的余项达到 $ O(\tau^4 \varepsilon^2 + \tau^p \varepsilon) $，其中 $ p = n+2 $ 或 $ p = n+3 $。
• 校正步骤将余项从 $ O(\tau^2 \varepsilon^2) $ 降低至 $ O(\tau^4 \varepsilon^2) $，显著提升了大时间步长下的精度。
• 在单摆和太阳-木星-土星三体问题上的数值测试表明，经校正的积分器在精度-成本比方面优于标准跃迁法和未经校正的高阶方法。
• 积分器的复合通常不具成本效益，除非步长极小且需要极高精度。
• 由于角动量与开普勒部分 $ A $ 对易，其为辛积分器的精确积分，可作为有效的误差检验工具。