[论文解读] Continuous-stage Runge-Kutta-NystrÖm methods

Wensheng Tang|arXiv (Cornell University)|Jul 9, 2018

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

一句话总结

该论文通过引入权函数并将积分区间推广至任意有限或无限区间，提出了一种连续阶段Runge-Kutta-Nyström（csRKN）方法的新框架。通过利用正交多项式展开并施加简化假设，该方法实现了对二阶微分方程的高阶辛和对称积分器的系统构造，数值结果证实了在哈密顿系统中能量和不变量的近似守恒。

ABSTRACT

We develop continuous-stage Runge-Kutta-NystrÖm (csRKN) methods in this paper. By leading weight function into the formalism of csRKN methods and modifying the original pattern of continuous-stage methods, we establish a new and larger framework for csRKN methods and it enables us to derive more effective RKN-type methods. Particularly, a variety of classical weighted orthogonal polynomials can be used in the construction of RKN-type methods. As an important application, new families of symmetric and symplectic integrators can be easily acquired in such framework. Numerical experiments have verified the effectiveness of the new integrators presented in this paper.

研究动机与目标

• 通过引入权函数并将积分区间推广至任意有限或无限区间，扩展经典csRKN框架。
• 使经典加权正交多项式可用于构造RKN型方法，以提升精度与稳定性。
• 通过简化假设与正交展开，建立系统化推导高阶辛和对称积分器的方法。
• 通过数值实验展示新积分器在哈密顿系统中保持几何不变量（如能量、角动量、Runge-Lenz-Pauli向量）的有效性。

提出的方法

• 将权函数引入csRKN形式化体系，将离散阶段集转化为由权函数定义测度的连续区间。
• 使用一般区间 $ I $ 定义新的csRKN方法，允许使用经典加权正交多项式（如Hermite、Legendre和Jacobi多项式）。
• 对阶条件施加简化假设，实现无需求解复杂非线性代数方程的高阶方法系统推导。
• 利用正交多项式展开将Butcher系数表示为光滑函数，便于解析处理并整合几何性质。
• 通过确保系数函数满足由简化假设导出的特定代数条件，构造辛和对称积分器。
• 采用配点法与Galerkin型方法进行数值实现，误差分析基于最大范数与欧几里得范数。

实验结果

研究问题

• RQ1 经典csRKN框架如何推广以允许更广泛的正交多项式基？
• RQ2 在csRKN公式中引入权函数是否能带来更高效、更灵活的RKN型积分器？
• RQ3 在何种程度上可借助正交多项式展开与简化假设系统化推导辛和对称积分器？
• RQ4 新积分器在哈密顿系统中对能量、角动量及Runge-Lenz-Pauli向量等几何不变量的保持效果如何？
• RQ5 新积分器在混沌与非混沌哈密顿问题的长时间模拟中表现出怎样的数值性能？

主要发现

• 新csRKN框架成功将经典加权正交多项式（如Hermite、Legendre和Jacobi）融入RKN型方法的构造中。
• Hermite-4方法在测试积分器中表现出最优精度，而Hermite-3方法因阶数较低导致性能较差。
• 所有辛积分器均表现出哈密顿量与Runge-Lenz-Pauli不变量的近似守恒，二次角动量则保持至机器精度。
• 数值轨道在Hénon-Heiles问题中保持在预期区域内，确认了正确的定性行为，未发生逃逸至等边三角形区域之外。
• 位置与动量变量的解误差随时间线性增长，符合辛积分器的预期行为。
• 借助简化假设与正交展开，该框架可在无需求解复杂非线性阶条件的情况下，系统推导出高阶辛与对称积分器。