The jump filter in the discontinuous Galerkin method for hyperbolic conservation laws

技术摘要：用于双曲守恒律的间断伽辽金法中的跳跃滤波器 (Jump Filter)

问题陈述
高阶线性数值格式，特别是间断伽辽金（DG）方法，因其高精度和对非结构网格的灵活性，被广泛用于求解双曲守恒律。然而，在模拟包含间断（如冲击波）的解时，这些格式往往会产生不希望出现的伪振荡（吉布斯现象）。这些振荡可能导致求解器发散或产生人工结构。虽然存在各种捕捉冲击的技术——例如全变分有界（TVB）限制器、WENO 限制器和人工粘性方法——但许多方法都存在缺点，包括计算成本高、在非结构或曲面网格上实现困难、需要特征分解，或在光滑区域产生过度的数值耗散从而降低精度。

方法论
作者提出了一种新型的跳跃滤波器 (jump filter)，旨在 DG 框架内运行，以减轻这些振荡，同时保持该方法的高阶精度和紧凑性。

• 局部粘性公式化：该方法的核心是引入一种基于单元界面处跳跃信息的局部粘性项。与全局谱粘性不同，这种方法是局部的。粘性系数 $\sigma_j$ 是根据单元边界处解及其导数的跳跃幅度自适应确定的。
  • 在跳跃极小的光滑区域，增加的粘性可以忽略不计，从而保持高阶精度。
  • 在具有强间断的区域，粘性会增加，以抑制产生振荡的高阶多项式模态。
• 通过时间分裂进行谱滤波：为了避免求解扩散项隐式方程带来的刚性和计算开销（这会限制时间步长），作者采用了时间分裂策略 (time-splitting strategy)。在标准显式 Runge-Kutta DG 格式的每个阶段之后，将跳跃滤波器作为后处理步骤应用。这涉及求解关于多项式解的模态系数的简单常微分方程 (ODE)：
  $$\frac{d}{dt} u^l_j + \sigma^l_j u^l_j = 0$$
  解通过将系数乘以指数衰减因子 $e^{-\sigma^l_j \Delta t}$ 来更新。该操作是显式的、计算廉价的，且不需要特征分解。
• 混合限制器集成：跳跃滤波器还集成到一个混合限制器框架中。在这种方法中，使用跳跃指示器来识别“问题”单元。对于这些单元，使用跳跃滤波器来计算一个低阶多项式解，然后将其与高阶 DG 解进行混合。这旨在进一步减少与仅使用滤波器相比的数值耗散。
• 适用性：该方法适用于一维和二维的标量及系统守恒律（具体为 Euler 方程）。它直接在物理（守恒）变量空间中运行，避免了特征分解的复杂性。

主要贡献

1. 新型跳跃滤波器：引入了一种利用单元界面处跳跃信息的局部粘性滤波器，用以自适应地阻尼高阶模态。
2. 保持 DG 特性：理论分析确认，所提方案保持了守恒性、$L^2$ 稳定性以及最优收敛阶（高阶精度）对于光滑解。
3. 计算效率：该方法设计为低计算成本。它不需要特征分解，完全在物理空间中运行，并保持了 DG 方案的紧凑性（局部数据通信），使其非常适合并行计算和 GPU 环境。
4. 混合集成：展示了该滤波器如何无缝集成到混合限制器框架中，在有效抑制强冲击波中伪振荡的同时，最大限度地减少了光滑区域的数值耗散。

数值结果
论文通过在 1D 和 2D Euler 方程上的广泛数值实验验证了该方法，涵盖了：

• 光滑流场：精度测试（例如 Burgers 方程、等熵涡旋）确认，在细网格上，该滤波器不会降低标准 DG 方案的收敛阶。
• 强冲击波：涉及 Lax 问题、爆炸波和双重马赫反射的测试表明，该滤波器能有效控制伪振荡，并捕捉到尖锐的冲击波轮廓，而不会出现未滤波高阶方案中常见的过冲现象。
• 复杂相互作用：在涉及冲击波-涡旋相互作用、冲击波-气泡相互作用以及 Kelvin-Helmholtz 不稳定性等场景中，该方法表现出色，其中的高阶方案配合滤波器能够比低阶方案更好地解析复杂的流动结构（如滑移线和湍流混合层）。
• 鲁棒性：该方案在高温高马赫数天体物理喷流模拟中表现出鲁棒性，并且在结合保正性限制器时能保持正定性。

意义与主张
作者声称，跳跃滤波器为传统的限制器和人工粘性方法提供了一个鲁棒、精确且高效的替代方案。其主要意义在于其简单性和计算效率：

• 它避免了许多其他高阶限制器所需的昂贵的特征分解。
• 它保持了 DG 方法的紧凑性，有利于高效并行化。
• 它提供了一种在需要之处引入必要耗散的无缝方式，平衡了冲击波捕捉能力与光滑区域高阶精度的保持。

论文得出结论，该方法有效地解决了使用不连续解求解双曲守恒律的挑战，为高分辨率冲击波捕捉提供了一个实用的工具。未来的工作指出将探索其在非结构网格和稳态问题上的实现。