[论文解读] Mimetic Explicit Time Discretiztions

Stanly Steinberg|arXiv (Cornell University)|May 27, 2016

Advanced Numerical Methods in Computational Mathematics参考文献 27被引用 3

一句话总结

本论文提出了一种拟谱有限差分方法，结合交错空间与显式时间离散化，以在时间步长约束下求解波动方程，实现二阶精度、能量守恒和无条件稳定性。该方法通过构建一个镜像连续系统结构的离散能量泛函，保持了关键连续体性质（如无散场和能量型不变量），在1D、2D和3D的标量方程与麦克斯韦方程中均得到验证，能量守恒精度达到机器精度。

ABSTRACT

This paper is part of a program to combine a staggered time and staggered spatial discretization of continuum mechanics problems so that any property of the continuum that is proved using vector calculus can be proven in an analogous way for the discretized system. We require that the discretizations be second order accurate and have a conserved quantity that approximates the energy for the system and guarantees stability for a reasonable constraint on the time step. We also require that the discretization is time explicit so as to avoid the solution of large system of possibly nonlinear algebraic equations. The well known Yee grid discretization of Maxwell's equations is the same as our discretization and is an early example of using a staggered space and time grid . To motivate our discussion we begin by studying the staggered time or leapfrog discretization of the harmonic oscillator and use this to introduce the modification of the energy that is conserved. Next we use systems of linear equations to motivate the definition of the modified energy for more complex systems of ordinary differential equations and then apply our ideas to the scalar wave equation in one spatial dimension. We finish by discretizing the three dimensional scalar wave and Maxwell's equations. Because the spatial discretization is mimetic, we obtain that the divergence of the electric and magnetic fields are constant when there are no sources. Using the mimetic properties the proof of this trivial and is essentially the same as in the continuum.

研究动机与目标

• 开发一种时间显式、二阶精度的有限差分格式，用于波动方程，以保持能量守恒和无散场等基本连续体性质。
• 构建一个离散能量泛函，其近似于连续能量并在时间演化下保持守恒，从而确保数值稳定性。
• 将先前用于椭圆与抛物问题的拟谱空间离散化方法，扩展至双曲波动系统，采用交错时空网格。
• 通过结合拟谱空间算子与跃迁时间推进法，实现非均匀与各向异性材料的稳定模拟。
• 通过图追逐法与离散星算子，证明离散系统继承了连续体的代数结构，包括精确序列性质与伴随关系。

提出的方法

• 对由二阶波动方程导出的一阶系统采用交错时间（跃迁）离散化，确保二阶精度与显式时间推进。
• 构造一个修正的离散能量泛函，其作为守恒量，类似于连续能量，并在足够小的时间步长下保持正定。
• 采用两个对偶网格（主网格与对偶网格），使一个网格的节点与另一网格的单元中心重合，从而实现拟谱空间算子（如梯度、旋度与散度）。
• 定义离散星算子（乘法算子），用于在主网格与对偶网格之间映射，保持对偶性并实现一致的内积。
• 应用离散精确序列，确保拟谱性质（如 div(curl) = 0 及算子兼容性）在离散层面得以保持。
• 通过混合（罗宾型）边界条件实现边界条件，且对角点值进行一致处理，避免编程不一致。

实验结果

研究问题

• RQ1 能否为波动方程构造一种时间显式、二阶精度的有限差分格式，使其守恒一个类似于连续能量的离散能量泛函？
• RQ2 如何将先前用于椭圆与抛物问题的拟谱空间离散化方法，扩展至双曲波动系统，同时保持关键连续体性质？
• RQ3 在交错网格上，确保波动方程时间二阶精度与稳定性的离散能量泛函具有何种结构？
• RQ4 在拟谱离散化下，离散微分算子（梯度、散度、旋度）及其伴随算子的行为如何？它们保留了哪些代数性质？
• RQ5 该方法在多大程度上可推广至非均匀与各向异性材料？这在离散算子与能量泛函中如何体现？

主要发现

• 在数值测试中，离散能量泛函的守恒精度达到至少 10^{-15}，证实了实际中的精确守恒。
• 该方法在二维标量波动方程中实现了时间与空间的二阶精度，部分情形在参数调优后甚至达到四阶精度或精确解。
• 拟谱空间离散化确保在无源条件下，电场与磁场的散度保持为零，与连续体一致。
• 离散系统继承了连续体的精确序列结构，离散算子因图追逐与星算子的正确定义，满足相同的代数恒等式（如 div(curl) = 0）。
• 通过使用对偶网格与离散星算子，实现了对各向异性与非均匀材料属性在波动方程中的一致处理。
• 边界条件通过混合（罗宾型）形式实现，角点值的处理不影响显式时间推进格式的解。