[论文解读] Symplectic Mimetic Finite Difference methods for Hamiltonian wave equations in 2D

L. Beirão da Veiga, L. Lopez+1 more|arXiv (Cornell University)|May 5, 2015

Numerical methods for differential equations被引用 3

一句话总结

本文提出了一种用于求解二维空间中哈密顿波方程的辛仿射有限差分方法。通过结合仿射空间离散化（保持连续系统的关键结构特性）与辛时间积分，该方法在长时间模拟中保持能量和动量守恒，从而确保哈密顿系统中具有高精度和稳定性。

ABSTRACT

In this paper we consider the numerical solution of the Hamiltonian wave equation in two spatial dimension. We use the Mimetic Finite Difference (MFD) method to approximate the continuous problem combined with a symplectic integration in time to integrate the semi-discrete Hamiltonian system. The main characteristic of MFD methods, when applied to stationary problems, is to mimic important properties of the continuous system. This approach, associated with a symplectic method for the time integration yields a full numerical procedure suitable to integrate Hamiltonian problems. A complete theoretical analysis of the method and some numerical simulations are developed in the paper.

研究动机与目标

• 开发一种能够保持二维空间中哈密顿波方程几何结构的数值方法。
• 通过保持辛性与仿射性质，解决波方程在长时间模拟中的数值不稳定性问题。
• 将空间中的仿射有限差分与时间上的辛积分器相结合，以增强守恒性与精度。
• 提供所提方法稳定性与收敛性的理论分析及数值验证。

提出的方法

• 空间离散化采用仿射有限差分（MFD）方法，以保持连续哈密顿系统的根本特性，如守恒律与边界条件。
• 通过辛积分器推进由MFD空间近似得到的半离散系统，以保持系统的几何结构。
• 该方法确保离散系统继承关键特性，如反对称性以及与连续算子的一致性。
• 该方法的构型设计确保在长时间模拟中保持动量与能量守恒。
• 对完整格式进行了完整的理论分析，以验证其稳定性、收敛性及结构保持性。

实验结果

研究问题

• RQ1 将仿射有限差分与辛时间积分相结合，能否保持二维波方程的哈密顿结构？
• RQ2 该方法在长时间积分过程中对能量与动量的守恒性能如何？
• RQ3 与标准有限差分格式相比，该方法的收敛行为与稳定性如何？
• RQ4 该方法在具有复杂边界条件的问题中是否仍能保持精度？

主要发现

• 所提方法通过空间与时间的离散化，成功保持了二维波方程的哈密顿结构。
• 数值模拟表明，即使在长时间时间区间内，该方法仍表现出长期能量守恒，这是辛方法的典型特征。
• 仿射空间离散化确保了与连续系统守恒律及边界行为的一致性。
• 理论分析证实，在标准假设下，该方法具有稳定性与收敛性。