[论文解读] Locally exact modifications of numerical integrators

Jan L. Cieśliński|arXiv (Cornell University)|Jan 3, 2011

Numerical methods for differential equations被引用 2

一句话总结

本文提出了一类局部精确的数值积分格式改进方法，可在每一点精确保持常微分方程的线性化，从而确保更高的精度与 A-稳定性。通过改进现有格式（尤其是离散梯度法），作者为任意多维哈密顿系统在正则坐标下构造了能量守恒、A-稳定的积分器，实现了数个数量级的精度提升，尤其在一维情况下，精度提升最高可达 8 个数量级。

ABSTRACT

We present a new class of exponential integrators for ordinary differential equations. They are locally exact, i.e., they preserve the linearization of the original system at every point. Their construction consists in modifying existing numerical schemes in order to make them locally exact. The resulting schemes preserve all fixed points and are A-stable. The most promising results concern the discrete gradient method (modified implicit midpoint rule) where we succeeded to preserve essential geometric properties and the final results have a relatively simple form. In the case of one-dimensional Hamiltonian systems numerical experiments show that our modifications can increase the accuracy by several orders of magnitude. The main result of this paper is the construction of energy-preserving locally exact discrete gradient schemes for arbitrary multidimensional Hamiltonian systems in canonical coordinates.

研究动机与目标

• 开发一类新的数值积分器，使其在每一点均具有局部精确性，精确保持常微分方程的线性化。
• 以保持其几何性质（如精确能量守恒）的方式，对现有数值格式（尤其是离散梯度法）进行改进。
• 将局部精确方法推广至任意多维哈密顿系统在正则坐标下的情形。
• 在保持 A-稳定性和几何结构的同时，显著提升数值精度。
• 实现能量守恒积分器中的变时间步长，克服传统辛方法的局限性。

提出的方法

• 该方法通过在每一点对现有数值格式进行精确线性化系统的离散化，构造局部精确积分器。
• 对线性化常微分方程使用精确矩阵指数解，利用雅克比矩阵的矩阵指数定义修正时间步长：δₙ = 2Ω⁻¹ₙ tan(hₙΩₙ/2)。
• 该方法利用离散梯度法，将其重新表述为具有能量守恒性质的两步格式。
• 关键方程为 δₙ⁻¹(xₙ₊₁ − xₙ) = ∇ₛᵀ(pₙ, pₙ₊₁)，确保对任意 δₙ 均实现精确能量守恒。
• 通过将原系统的线性化动力学与线性化常微分方程的精确解匹配，利用雅克比矩阵的矩阵指数推导出该格式。
• 通过在平衡点处定义 Ωₙ² = TₚₚVₓₓ，将该方法推广至多维哈密顿系统。

实验结果

研究问题

• RQ1 能否对数值积分器进行改进，使其在保持能量守恒等几何性质的同时具备局部精确性？
• RQ2 如何系统性地构造适用于一般多维哈密顿系统的局部精确改进方法？
• RQ3 与标准离散梯度法相比，局部精确格式在精度上可提升多少？
• RQ4 能否在不破坏几何结构的前提下，将变时间步长引入能量守恒积分器？
• RQ5 是否能够构造出保持精确能量守恒的离散梯度格式的局部精确改进？

主要发现

• 所提出的离散梯度法的局部精确改进在任意多维哈密顿系统（正则坐标下）中，能精确保持能量积分。
• 在一维哈密顿系统中，该方法相比标准离散梯度格式，精度最高可提升 8 个数量级。
• 所得格式具有 A-稳定性，并精确保持原系统的所有不动点。
• 该方法支持变时间步长，可实现传统辛方法无法实现的自适应积分策略。
• 精确线性化系统的离散化作为基础，所得格式在线性情况下退化为已知的精确谐振子解。
• 该构造方法具有普适性，可通过在平衡点处计算雅克比矩阵并使用矩阵指数，应用于任意自治常微分方程系统。