[论文解读] Stability, convergence, and geometric properties of second-order-in-time space-time discretizations for linear and semilinear wave equations

Matteo Ferrari, Ilaria Perugia+1 more|arXiv (Cornell University)|Jan 6, 2026

Numerical methods for differential equations被引用 0

一句话总结

本论文建立了二阶时间的时空伽辽金离散化与一阶时间格式之间的等价性，适用于线性和半线性波方程，得到无条件稳定性、收敛性结果，并对能量守恒和辛结构性质提供见解。还通过高斯-勒让德和高斯-洛巴托时间积分引入了辛结构的变体。

ABSTRACT

We revisit second-order-in-time space-time discretizations of the linear and semilinear wave equations by establishing precise equivalences with first-order-in-time formulations. Focusing on schemes using continuous piecewise-polynomial trial functions in time, we analyze their stability, convergence, and geometric properties. We consider first a weak space-time formulation with test functions projected onto discontinuous polynomials of one degree lower in time, showing that it is equivalent to the scheme proposed in [French, Peterson 1996] in the linear case, and extended in [Karakashian, Makridakis 2005] to the semilinear case. In particular, this equivalence shows that this method conserves energy at mesh nodes but is not symplectic. We then introduce two symplectic variants, obtained through Gauss-Legendre and Gauss-Lobatto quadratures in time, and show that they correspond to specific Runge-Kutta time integrators. These connections clarify the geometric structure of the space-time methods considered.

研究动机与目标

• 为线性与半线性波方程的二阶时间-时空伽辽金方法提供动机与分析。
• 证明与一阶时间格式的精确等价性，以便转移已知性质。
• 研究稳定性、收敛性及几何性质，包括能量守恒与辛性。
• 提出实现技术，避免求解投影的辅助局部系统。

提出的方法

• 给出一个弱的时空变分形式，测试函数在时间上投影到一个较低次数的一致多项式。
• 证明对线性波方程，带稳定化的二阶时间方法等价于一阶时间 DG–CG 方法。
• 把该框架扩展到半线性波方程，在所有项中除了时间刚性项以外，进行一个次数较低的投影。
• 在时间上引入两种辛结构变体，利用高斯-勒让德和高斯-洛巴托积分，对应 Runge–Kutta 时间积分。
• 通过借助与一阶格式建立的等价性，证明无条件稳定性和收敛性。
• 在网格节点处 Demonstrate energy conservation（在非线性情况下近似于四舍五入误差，且讨论基于积分的辛结构保持）。

实验结果

研究问题

• RQ1 二阶时间的时空伽辽金离散是否可以严格地重写为一阶时间格式？
• RQ2 这些二阶方案是否继承了已知的一阶 DG–CG 方法的稳定性和收敛性？
• RQ3 降阶的时间离散空间投影对稳定性和能量守恒有何影响？
• RQ4 是否能通过高斯-勒让德和高斯-洛巴托时间积分实现辛性质，并且对稳定性有何影响？
• RQ5 在线性与半线性波方程中，在能量守恒与几何性质方面有哪些差异？

主要发现

• 二阶时间-时空离散化与一阶时间 DG–CG 公式等价，便于转移稳定性与收敛性结果。
• 带有时间方向投影到较低次数的一致多项式的稳定化变体，在张量积网格上具有无条件稳定性。
• 基于高斯-勒让德和高斯-洛巴托积分的辛变体对应于 Runge–Kutta 时间积分，并在四舍五入误差范围内保持能量等结构。
• 在 F=0 的线性情形下，能量在网格节点处被严格守恒，非线性情形下能量守恒延伸到四舍五入误差范围内。
• 对于半线性问题，该框架保持无条件稳定性与收敛性，同时将未知量减半，便于时间步进实现。
• 该方法通过时间方向的子积分实现，无需引入辅助变量即可实现。
• 数值结果验证线性情形下的能量保持，以及正弦-Gordon 型模型的收敛性速率。