[论文解读] Energy-preserving continuous-stage partitioned Runge-Kutta methods

Wensheng Tang|arXiv (Cornell University)|Aug 7, 2018

Numerical methods for differential equations参考文献 42被引用 3

一句话总结

本文通过推导能量保持的充分条件，提出了一类用于哈密顿系统的能量保持连续阶分区 Runge-Kutta（csPRK）方法。该方法在使用归一化移位勒让德多项式及简化假设时，通过要求 Butcher 权系数 $B_\tau$ 和 $\widehat{B}_\tau$ 均等于 1，实现精确的能量守恒，从而构造出高阶积分器，并通过非线性系统上的数值实验验证了其能量守恒精度接近机器精度。

ABSTRACT

In this paper, we present continuous-stage partitioned Runge-Kutta (csPRK) methods for energy-preserving integration of Hamiltonian systems. A sufficient condition for the energy preservation of the csPRK methods is derived. It is shown that the presented condition contains the existing condition for energy-preserving continuous-stage Runge-Kutta methods as a special case. A noticeable and interesting result is that when we use the simplifying assumptions of order conditions and the normalized shifted Legendre polynomials for constructing high-order energy-preserving csPRK methods, both the Butcher "weight" coefficients $B_τ$ and $\widehat{B}_τ$ must be equal to $1$. As illustrative examples, new energy-preserving integrators are acquired by virtue of the presented condition, and for the sake of verifying our theoretical results, some numerical experiments are reported.

研究动机与目标

• 开发用于哈密顿系统的能量保持数值积分器，使其在长时间模拟中精确保持总能量。
• 将连续阶 Runge-Kutta 方法的理论扩展至分区形式，以更好地保持几何结构。
• 为 csPRK 方法建立能量保持的充分条件，推广现有连续阶 Runge-Kutta 方法的结果。
• 通过简化假设与归一化移位勒让德多项式构造高阶能量保持积分器，确保 $B_\tau = 1$ 与 $\widehat{B}_\tau = 1$ 同时成立。
• 通过非线性哈密顿系统上的数值实验验证理论结果，比较能量与不变量守恒性能与辛方法的表现。

提出的方法

• 利用级数展开技术推导连续阶分区 Runge-Kutta（csPRK）方法中能量保持的充分条件。
• 应用简化假设以满足阶条件，并利用归一化移位勒让德多项式构造高阶积分器。
• 施加约束条件，要求 Both Butcher 权系数 $B_\tau$ 与 $\widehat{B}_\tau$ 必须等于 1，以确保能量保持。
• 将方法重新表述为非线性哈密顿系统，以支持数值测试与验证。
• 采用三点高斯积分进行数值积分，并与二阶四阶高斯-勒让德 RK 方法（GLRK-4）进行比较。
• 通过 $H$、$I$ 与 $L$ 的全局解误差和不变量误差的最大范数评估性能。

实验结果

研究问题

• RQ1 能否为哈密顿系统中的连续阶分区 Runge-Kutta 方法推导出能量保持的充分条件？
• RQ2 在简化假设与基于勒让德多项式的构造下，Butcher 权系数 $B_\tau$ 与 $\widehat{B}_\tau$ 在高阶能量保持方法中的行为如何？
• RQ3 与辛方法（如 GLRK-4）相比，所提出的 csPRK 方法在能量与其他不变量守恒方面表现如何？
• RQ4 所推导的能量保持 csPRK 方法是否能在混沌或刚性系统中实现线性误差增长并保持准确的长期动力学？
• RQ5 所提出的方法是否可超越标准 Runge-Kutta 框架，特别是在 P 系列积分器的背景下适用？

主要发现

• 推导出的 csPRK 方法能量保持的充分条件，推广了现有连续阶 Runge-Kutta 方法的条件。
• 当使用简化假设与归一化移位勒让德多项式时，为实现能量保持，$B_\tau$ 与 $\widehat{B}_\tau$ 必须精确等于 1。
• 数值实验表明，所有提出的能量保持方法将哈密顿量 $H$ 守恒至机器精度，而辛方法 GLRK-4 精确守恒二次角动量 $I$。
• 能量保持方法在全局解误差与不变量误差中表现出线性误差增长，与辛方法相当。
• 所有能量保持方法的数值轨道均紧密逼近精确的椭圆轨迹，表明其具有高精度与稳定性。
• 尽管非辛，所提出的 csPRK 方法在保持几何与动力学特性方面与辛方法表现相当。