Stability analysis of Arbitrary-Lagrangian-Eulerian ADER-DG methods on classical and degenerate spacetime geometries

技术摘要：任意拉格朗日-欧拉（ALE）ADER-DG 方法的稳定性分析

问题陈述
任意拉格朗日-欧拉（ALE）框架广泛用于在移动网格上求解双曲型偏微分方程（PDE），特别是在发生拓扑结构变化的情况下。在此背景下，直接 ALE 方法通过时空控制体将不同时间步的网格连接起来。Gaburro 等人的近期研究引入了“退化”时空单元（狭缝单元/sliver elements），用于处理拓扑变化场景，即单元在时间步的开始或结束时空间尺寸可能为零，但具有非零的时空体积。虽然 ADER-DG（任意高阶导数间断伽辽金）方法在解决这些问题方面非常有效，但针对这些方案——特别是针对退化时空几何体——进行严谨的理论稳定性分析一直处于缺失状态。现有文献通常依赖于经验性的 CFL（Courant-Friederichs-Lewy）界限，或仅限于低多项式阶数的有限稳定性研究，这导致在理解高阶显式和隐式 ALE ADER-DG 方法在经典及退化几何体上的稳定性约束方面存在空白。

方法论
作者采用冯·诺依曼（von Neumann）稳定性分析来研究显式和隐式 ALE ADER-DG 方法。本研究以线性平流方程作为双曲系统的代理模型。

1. 公式化： 论文首先建立了显式和隐式 ALE ADER-DG 方案的数学框架。
   • 显式方案： 采用预测-校正（predictor-corrector）方法。预测步通过不动点迭代，在每个控制体（或狭缝单元）内构建局部时空多项式逼近。校正步利用移动基函数和数值通量（Rusanov 型）更新解。
   • 隐式方案： 通过 Newton-GMRES 迭代，在所有控制体内同时求解全局时空多项式，直接对 PDE 在时空控制体上进行积分。
2. 退化几何建模： 为了分析退化几何，作者引入了一个一维替代设置，其中标准控制体之间的界面被“狭缝单元”取代。这些单元模拟了用于 2D/3D 拓扑变化的类孔洞狭缝，其在 $t_n$ 和 $t_{n+1}$ 时空间宽度为零，但具有非零的时空体积。
3. 稳定性分析：
   • 经典几何： 作者推导了显式和隐式方案的放大矩阵。他们在一定范围内的 CFL 数和相位角上计算了矩阵的谱半径 $\rho$。
   • 退化几何： 分析被扩展到包含狭缝单元的退化几何。该域被建模为一个包含标准控制体和狭缝单元的周期性块。构造了一个依赖于 CFL 数、狭缝宽度参数 $\delta$ 和相位角的放大矩阵。
   • 数值验证： 通过在精细的 CFL 值和多项式阶数（$N=1$ 到 $9$）网格上计算离散放大因子，验证了理论稳定性界限。一致性阶数也通过数值进行了验证。

主要贡献与结果

• 完善显式稳定性界限： 对于经典几何，研究证实对于低多项式阶数（$N \le 3$），其稳定性限制与广泛使用的经验 CFL 界限一致。然而，对于高阶情况（$N \ge 4$），本文证明了常用的经验 CFL 值实际上违反了严格的冯·诺依曼稳定性条件（即 $\rho > 1$）。作者确定了高阶显式方案实现严格稳定性所需的显著更低的、严谨的 CFL 限制。
• 网格速度约束： 该分析为固定目标 CFL 下的可容许网格速度范围提供了理论特征，展示了网格运动如何影响稳定性约束。
• 隐式无条件稳定性： 冯·诺依曼分析和理论证明（定理 1）确认，对于线性平流方程，隐式 ALE ADER-DG 方法在经典几何上是无条件稳定的，无论时间步长大小如何。
• 退化几何上的稳定性：
  • 显式情况： 狭缝单元的引入并未降低稳定性。研究发现，退化设置下的离散放大因子小于或等于经典设置下的放大因子。因此，适用于经典几何的相同 CFL 界限同样适用于退化几何。作者指出，在显式全局方案中对狭缝预测步进行隐式处理可能会略微增加可接受的 CFL 上限。
  • 隐式情况： 只要狭缝宽度参数 $\delta$ 相对于时间步和网格大小满足特定的几何约束（定理 2），隐式方法即使在存在狭缝单元的情况下仍保持无条件稳定性。
• 一致性： 数值实验确认，两种方案在经典和退化几何上均保持了预期的一致性阶数（$N+1$）。

意义与主张
本文声称填补了关于退化时空几何上 ALE ADER-DG 方法的第一个严谨冯·诺依曼稳定性分析的理论空白。其主要意义在于验证了使用退化单元（狭缝）来连接发生拓扑变化的移动网格的可行性。结果表明，使用这些零尺寸的空间单元不会引入额外的稳定性限制，也不会降低与经典几何相比的可允许时间步长。这一发现为在涉及拓扑变化的复杂场景中实际应用直接 ALE 方法提供了理论基础，并为开发新的基于时空切割单元（spacetime cut cell）的方法铺平了道路。作者强调，虽然在实践中经常使用经验 CFL 值（通过粘性和限制器进行缓解），但所推导的严谨界限对于全面理解该方法的行为至关重要。