[论文解读] Symplecticity-preserving continuous-stage Runge-Kutta-Nyström methods

Wensheng Tang, Jing-Jing Zhang|arXiv (Cornell University)|Feb 4, 2016

Numerical methods for differential equations参考文献 10被引用 5

一句话总结

本文提出了一种保持辛结构的连续阶Runge-Kutta-Nyström（csRKN）方法，用于求解由可分离哈密顿系统产生的二阶常微分方程。通过利用勒让德多项式展开并推导出充分的辛条件，作者构建了高阶、对角隐式及显式辛积分器，相较于经典方法（需求解复杂的非线性阶条件），能更高效地保持哈密顿动力学的几何结构。

ABSTRACT

We develop continuous-stage Runge-Kutta-Nyström (csRKN) methods for solving second order ordinary differential equations (ODEs) in this paper. The second order ODEs are commonly encountered in various fields and some of them can be reduced to the first order ODEs with the form of separable Hamiltonian systems. The symplecticity-preserving numerical algorithm is of interest for solving such special systems. We present a sufficient condition for a csRKN method to be symplecticity-preserving, and by using Legendre polynomial expansion we show a simple way to construct such symplectic RKN type method.

研究动机与目标

• 开发能够保持可分离哈密顿系统中二阶ODE辛结构的连续阶RKN方法。
• 建立csRKN方法中辛性的充分条件，以实现结构保持积分器的系统化构造。
• 通过避免求解经典方法中典型的复杂非线性阶条件，简化高阶辛积分器的构造过程。
• 证明csRKN方法可利用不同的积分规则，从同一连续阶公式中导出多个辛积分器。
• 提供具有自由参数的四阶显式及对角隐式辛积分器，包括已知的Störmer-Verlet方法等。

提出的方法

• 通过分析内部阶段的连续阶公式及其对Butcher表函数的依赖关系，推导出csRKN方法中辛性的充分条件。
• 采用勒让德多项式展开来参数化系数函数 $ A_{\tau,\sigma} $, $ \hat{A}_{\tau,\sigma} $, $ B_\tau $, 和 $ \hat{B}_\tau $，从而实现辛方法的解析构造。
• 对展开后的系数施加辛性与阶条件，导出可求解自由参数的代数方程组。
• 使用高斯积分规则将连续阶方法离散化为经典RKN格式，得到显式与对角隐式积分器。
• 通过在系数上施加特定代数约束，构造出两参数族的对角隐式辛积分器。
• 证明可通过在 $ \bar{A}_{ij} $ 矩阵上进一步施加上三角性（下三角为零）条件，从csRKN方法中导出显式辛积分器。

实验结果

研究问题

• RQ1 能否构造出保持可分离哈密顿系统辛结构的连续阶RKN方法？
• RQ2 何种充分条件可确保csRKN方法为辛方法？其在系数函数中的表达形式为何？
• RQ3 如何利用勒让德多项式展开系统化构造具有自由参数的辛csRKN方法？
• RQ4 能否通过使用不同的积分规则，从同一csRKN公式中导出不同的经典RKN格式？
• RQ5 积分规则的阶数与所得离散RKN方法的阶数之间存在何种关系？

主要发现

• 推导出csRKN方法中辛性的充分条件，使得无需求解完整的非线性阶条件系统，即可构造辛积分器。
• 构造出两参数族的四阶对角隐式辛积分器，其显式形式见表3及公式(5.6)–(5.8)。
• 通过在系数矩阵上施加额外约束，导出显式辛积分器，得到三个新显式格式，其中之一为著名的Störmer-Verlet方法。
• 当 $ \alpha = \frac{1}{8}, \beta = \frac{\sqrt{3}}{24}, \gamma = -\frac{1}{8} $ 时，Störmer-Verlet格式可作为特例被恢复，验证了与已有结果的一致性。
• 通过将方法扩展至更高阶积分规则并求解所得代数系统，可构造更高阶辛积分器（如三阶段、四阶）。
• 与经典方法相比，该方法显著简化了高阶辛RKN方法的构造过程，因其避免了求解来自阶条件的复杂非线性代数系统。