[论文解读] Numerical comparisons between Gauss-Legendre methods and Hamiltonian BVMs defined over Gauss points

Luigi Brugnano, Felice Iavernaro+1 more|arXiv (Cornell University)|Feb 14, 2010

Numerical methods for differential equations被引用 6

一句话总结

该论文提出了一类能量保持型数值积分方法——哈密顿边界值方法（HBVMs），适用于多项式哈密顿系统。HBVMs通过施加扩展的配点条件，将高斯-勒让德配点法推广至任意次数的多项式哈密顿量，从而实现哈密顿量的精确守恒。与仅能保持二次哈密顿量的常规高斯方法不同，HBVMs在任意多项式次数下均保持能量守恒，数值实验表明，即使采用较大的步长，其哈密顿量守恒精度也达到机器精度，而高斯方法则累积能量误差，需减小步长或提高阶数才能降低误差。

ABSTRACT

Hamiltonian Boundary Value Methods are a new class of energy preserving one step methods for the solution of polynomial Hamiltonian dynamical systems. They can be thought of as a generalization of collocation methods in that they may be defined by imposing a suitable set of extended collocation conditions. In particular, in the way they are described in this note, they are related to Gauss collocation methods with the difference that they are able to precisely conserve the Hamiltonian function in the case where this is a polynomial of any high degree in the momenta and in the generalized coordinates. A description of these new formulas is followed by a few test problems showing how, in many relevant situations, the precise conservation of the Hamiltonian is crucial to simulate on a computer the correct behavior of the theoretical solutions.

研究动机与目标

• 开发能够精确保持多项式哈密顿系统哈密顿函数的数值方法，尤其适用于哈密顿量为高次多项式的情形。
• 解决经典高斯-勒让德配点法的局限性，后者仅能保持二次哈密顿量，对更高次多项式则失效。
• 证明能量守恒对于正确模拟长期动力学至关重要，特别是在周期环带边界附近。
• 在长积分区间内，比较HBVMs与高斯-勒让德方法在能量守恒精度和轨道保真度方面的表现。
• 验证HBVMs可实现机器精度的周期环带边界定位，而高斯方法则引入显著的能量漂移。

提出的方法

• HBVMs采用离散线积分公式构建，该公式推广了保守向量场的连续线积分。
• 方法在相空间的多项式路径上施加扩展的配点条件，确保离散能量积分精确为零。
• 基于移位勒让德多项式构造正交多项式，用以定义离散解路径的基函数。
• 该方法被实现为具有 $k$ 个隐式阶段和 $s$ 个阶段的BVM，其中 $k \geq 2s-1$ 可确保对次数不超过 $2k-1$ 的多项式实现哈密顿量的精确守恒。
• 通过求解由扩展配点条件导出的非线性系统推进解，从而精确保持哈密顿量。
• 通过在相同步长和容差下积分同一哈密顿系统，将HBVMs与高斯-勒让德方法进行比较，测量能量误差和轨道精度。

实验结果

研究问题

• RQ1 能否设计一种数值方法，使其对任意次数的多项式哈密顿量（而不仅限于二次）实现精确的能量守恒？
• RQ2 在积分高次多项式哈密顿量系统时，HBVMs的能量守恒性能与高斯-勒让德方法相比如何？
• RQ3 能量漂移对数值轨道长期行为有何影响，特别是在周期环带边界附近？
• RQ4 HBVMs是否能以机器精度定位周期环带边界，而与步长或阶数无关？
• RQ5 HBVMs对能量的精确守恒是否导致其轨道轨迹比高斯方法更精确、更稳定，特别是在鞍点区域？

主要发现

• HBVMs（阶数4, 6, 8, 10，对应 $k=5,8,10,13$，$s=2,3,4,5$）计算周期环带边界点时，相对能量误差 $\leq 4.6 \times 10^{-16}$，达到机器精度。
• 相比之下，高斯方法（HBVM($s$,$s$)）在 $h=1$ 时相对能量误差高达 $2.15 \times 10^{-2}$，仅通过减小步长或提高阶数才能降低。
• HBVM(5,2)计算的数值轨道精确通过鞍点 $P_1$，而同阶高斯方法因能量漂移导致轨道明显不规则。
• 即使高斯方法中哈密顿量误差极小，也导致周期环带边界附近出现明显的轨道不规则性，如图11所示。
• HBVMs计算的点 $y_0^{(k,s)}$ 在所有仿真中保持一致，表明边界检测具有高度一致性和准确性。
• HBVMs在任意步长（$h=1$ 或 $h=0.5$）下均保持机器精度的能量守恒，而高斯方法需减小步长才能提高精度。