[论文解读] How to regularize a symplectic-energy-momentum integrator

Yosi Shibberu|ArXiv.org|Jul 22, 2005

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

一句话总结

本论文提出了一种针对辛-能量-动量（SEM）积分器的正则化方法，通过解决离散哈密顿表述中奇点引起的非物理解析轨迹（鬼魅轨迹），恢复了时间可逆性。通过引入对偶变分原理，并利用带括号的根求解方法求解能量守恒的中点，该方法在保持辛结构、能量和动量守恒的同时，实现了即使在相空间临界点（如单摆转折点）附近也稳定的、时间可逆的积分。

ABSTRACT

We identify ghost trajectories of symplectic-energy-momentum (SEM) integration and show that the ghost trajectories are not time reversible. We explain how SEM integration can be regularized, in a SEM preserving manner, so that it is time reversible. We describe an algorithm for implementing the regularized SEM integrator. Simulation results for the pendulum are given. Coordinate invariance of the regularized SEM integrator is briefly discussed.

研究动机与目标

• 解决在奇点（如单摆转折点）附近，辛-能量-动量（SEM）积分器出现的时间不可逆性与数值不稳定性问题。
• 开发一种正则化技术，使SEM积分器在保持辛结构、能量和动量守恒的同时，恢复时间可逆性。
• 提供一种计算可行的算法，利用改进的中点格式与自适应时间步长，计算正则化轨迹。
• 通过单摆和开普勒单体问题的数值模拟，验证该方法的有效性。

提出的方法

• 该方法引入扩展相空间表述，将时间视为具有共轭动量 ℘ 的动力学变量，从而基于最小作用量原理建立离散变分原理。
• 鬼魅轨迹被识别为在中点满足哈密顿约束 H(z) = 0，但因时间演化非单调而违反时间可逆性的解。
• 采用对偶变分原理对轨迹进行正则化：一个用于鬼魅路径（ψ ≠ 0），一个用于常规路径（ψ = 0），其中 ψ 为与时间步长符号相关的约束函数。
• 算法使用牛顿法求解 H(𝑧̄(λ, zₖ)) = 0 以确定中点 𝑧̄ₖ，并采用带括号的根求解策略以处理在 ψ = 0 附近的病态问题。
• 当穿过 ψ = 0 时，算法切换至使用 μₖ ≠ 0 的正则化路径，确保时间可逆性并避免鬼魅解。
• 应用重连程序以防止 ψ = 0 的立即再次穿越，从而保持轨迹连续性与数值稳定性。

实验结果

研究问题

• RQ1 辛-能量-动量积分器能否通过正则化方法在存在奇点（如单摆转折点）的情况下恢复时间可逆性？
• RQ2 SEM积分器中非时间可逆鬼魅轨迹的产生机制是什么？
• RQ3 如何修改离散变分原理，以在确保时间可逆性的同时，保持辛结构、能量和动量守恒？
• RQ4 何种数值算法能够实现正则化SEM轨迹在标准积分器失效区域的鲁棒计算？
• RQ5 负时间步长与多值时间依赖关系如何影响正则化SEM积分器的长期行为？

主要发现

• 正则化过程成功通过消除违反时间可逆性的鬼魅轨迹，使SEM积分器恢复了时间可逆性。
• 在每一步的中点，能量守恒至舍入误差水平；在顶点处动量守恒，证实了SEM性质的精确保持。
• 通过验证映射 zₙ(z₀) 的雅可比矩阵满足 (dzₙ/dz₀)ᵀJ(dzₙ/dz₀) = J 至舍入精度，数值验证了辛结构的保持。
• 该算法可处理负时间步长，并通过重连机制确保最终输出的时间演化保持单调性。
• 在单摆分离子附近，当 ψₖ = 0 时观察到类似混沌的行为，提示可通过引入小的 μₖ ≠ 0 进行进一步正则化。
• 单摆与开普勒单体问题的数值结果表明，正则化积分器保持了所有SEM性质，并在相空间临界区域表现出鲁棒性。