[论文解读] Hamiltonian Boundary Value Methods (Energy Conserving Discrete Line Integral Methods)

Luigi Brugnano, Felice Iavernaro+1 more|arXiv (Cornell University)|Oct 19, 2009

Numerical methods for differential equations参考文献 17被引用 171

一句话总结

本文提出哈密顿边界值方法（HBVMs）作为一类能量守恒的数值积分方法，用于哈密顿系统，通过离散线积分和静默阶段精确保持多项式哈密顿量。主要贡献在于推导出静默阶段趋于无穷时的极限，从而得到新型方法——无穷HBVMs，其达到最优阶数2s，并推广了能量守恒配置法。

ABSTRACT

Recently, a new family of integrators (Hamiltonian Boundary ValueMethods) has been introduced, which is able to precisely conserve the energy function of polynomial Hamiltonian systems and to provide a practical conservation of the energy in the non-polynomial case. We settle the definition and the theory of such methods in a more general framework. Our aim is on the one hand to give account of their good behavior when applied to general Hamiltonian systems and, on the other hand, to find out what are the optimal formulae, in relation to the choice of the polynomial basis and of the distribution of the nodes. Such analysis is based upon the notion of extended collocation conditions and the definition of discrete line integral, and is carried out by looking at the limit of such family of methods as the number of the so called silent stages tends to infinity.

研究动机与目标

• 建立基于离散线积分和扩展配点条件的哈密顿边界值方法（HBVMs）的统一框架。
• 研究当静默阶段数量趋于无穷时HBVMs的极限，阐明其与现有能量守恒方法的联系。
• 通过分析多项式基和节点分布对收敛性和阶数的影响，识别最优HBVM公式。
• 引入并表征一类基于移位勒让德多项式的极限方法——无穷HBVMs。
• 证明极限方法与所用节点分布无关，且收敛至相同的算子特征函数，确保各类公式间行为一致。

提出的方法

• 使用离散线积分形式化HBVMs，作为连续线积分的推广，确保对多项式哈密顿量实现精确能量守恒。
• 将静默阶段定义为基本阶段的线性组合，实现在不增加计算成本（超过s个阶段）的前提下保持能量守恒。
• 通过扩展配点条件推导HBVM族，其中静默阶段数r被选择为匹配哈密顿量多项式的次数ν。
• 取k = s + r → ∞的极限，推导极限方法，证明其收敛至唯一极限，且与节点分布无关。
• 建立使用拉格朗日基的HBVM极限与能量守恒配置法（EPCMs）之间的等价性，并基于移位勒让德基推导新的极限方法。
• 证明所得的无穷HBVMs（HBVM(∞,s)）达到阶数2s，与高斯-勒让德方法阶数一致，且与节点选择无关。

实验结果

研究问题

• RQ1 当静默阶段数量趋于无穷时，HBVMs的理论极限是什么？
• RQ2 多项式基的选择（如拉格朗日基与移位勒让德基）如何影响所得极限方法？
• RQ3 极限方法能否独立于配点过程中使用的节点分布进行表征？
• RQ4 极限方法是否保持与经典高斯-勒让德方法相同的阶数？
• RQ5 在实际应用中，极限方法与现有能量守恒积分器相比，在能量守恒性和稳定性方面表现如何？

主要发现

• 当静默阶段k → ∞时，HBVMs的极限为唯一极限方法，且与配点节点的分布无关。
• 使用拉格朗日基时，极限方法与能量守恒配置法（EPCMs）一致，证实了Hairer的猜想。
• 使用移位勒让德基可导出一类新的极限方法——无穷HBVMs（HBVM(∞,s)），其阶数为2s，与高斯-勒让德方法相同。
• 极限方法由某一算子的单位特征值对应的特征函数表征，该特征函数与节点位置无关。
• 数值实验表明，HBVM(10,2)在长时间积分中能量守恒精度可达10−12量级，优于标准方法如Itoh-Abe和四阶格式。
• 由于在机器精度下多项式对哈密顿量的高阶逼近，该方法即使对非多项式哈密顿量也表现出实际的能量守恒。