[论文解读] Analisys of Hamiltonian Boundary Value Methods (HBVMs) for the numerical solution of polynomial Hamiltonian dynamical systems

Luigi Brugnano, Felice Iavernaro+1 more|arXiv (Cornell University)|Sep 30, 2009

Numerical methods for differential equations被引用 18

一句话总结

本文提出了一类新型的对称、A-稳定数值积分器——哈密顿边界值方法（HBVMs），用于多项式哈密顿系统。HBVMs通过离散变分公式，精确保持任意次数的多项式哈密顿量，实现任意高阶精度，同时在长期模拟中保持超越经典辛方法的守恒性质。

ABSTRACT

One main issue, when numerically integrating autonomous systems, is the long-term conservation of some of its invariants, among which the function itself. For example, it is well known that classical symplectic methods can only exactly preserve, at most, quadratic Hamiltonians. In this paper, a new family of methods, called Hamiltonian Boundary Value Methods (HBVMs), is introduced and analyzed. HBVMs are able to exactly preserve, in the discrete solution, functions of polynomial type of arbitrarily high degree. These methods turn out to be symmetric, precisely A-stable, and can have arbitrarily high order. A few numerical tests confirm the theoretical results.

研究动机与目标

• 为解决经典辛方法仅能精确保持二次哈密顿量的局限性。
• 开发一种可精确保持任意高次多项式哈密顿量的数值积分方法。
• 确保所提出方法在长期模拟中保持理想的稳定性与对称性。
• 分析新方法的理论性质，包括阶数、A-稳定性与对称性。
• 通过基准动力系统的数值实验验证理论结果。

提出的方法

• HBVMs通过在特定高斯-勒让德点处进行配点的离散变分公式构建，以强制精确保持多项式哈密顿量。
• 该方法在时间分区上采用配点法，其中离散解满足修正的哈密顿边界条件。
• 该公式确保即使对于高次多项式，数值解仍位于哈密顿函数的精确水平集上。
• 该方法在设计上具有对称性，从而保证了有利的稳定性特性，包括A-稳定性。
• 通过增加配点数量并相应调整求积规则，可使精度阶数任意提高。
• 该方法通过将哈密顿结构直接嵌入离散变分原理，推广了经典配点方法。

实验结果

研究问题

• RQ1 能否设计一种数值积分器，以精确保持任意次数的多项式哈密顿量，突破经典辛方法仅限于二次哈密顿量的限制？
• RQ2 此类方法的稳定性和对称性特性如何？能否实现A-稳定性？
• RQ3 该方法能否在保持哈密顿结构的同时实现任意高阶精度？
• RQ4 离散变分公式如何确保数值解中哈密顿量的长期守恒？
• RQ5 HBVMs在实际模拟中的性能与标准辛积分器相比如何？

主要发现

• HBVMs在离散解中精确保持任意次数的多项式哈密顿量，克服了经典辛方法仅限于二次哈密顿量的局限性。
• 该方法具有对称性与A-稳定性，确保在长期模拟中具有鲁棒性与稳定性。
• 通过增加配点数量，HBVMs可实现任意高阶精度。
• 数值测试验证了理论预测，显示在长时间积分区间内哈密顿量的守恒性能极佳。
• 离散变分公式成功地将哈密顿结构嵌入数值格式，实现能量函数的精确保持。
• 该方法在保持高次多项式不变量方面优于经典辛积分器，后者通常无法精确保持此类不变量。