[论文解读] Enhancing stability of correction procedure via reconstruction using summation-by-parts operators II: Modal filtering

Jan Glaubitz, Hendrik Ranocha+2 more|arXiv (Cornell University)|Jun 3, 2016

Medical Imaging Techniques and Applications参考文献 12被引用 5

一句话总结

该论文通过在求和分割（SBP）框架内引入模态滤波，增强了双曲守恒律的校正过程重构（CPR）格式的稳定性。它建立了模态滤波与谱黏性在首阶精度内的等价性，提出了自适应滤波强度策略，并证明了使用显式欧拉时间积分的全离散格式的可证明稳定性，克服了先前人工耗散方法中存在的时间步长限制。

ABSTRACT

A recently introduced framework of semidiscretisations for hyperbolic conservation laws known as correction procedure via reconstruction (CPR, also known as flux reconstruction) is considered in the extended setting of summation-by-parts (SBP) operators using simultaneous approximation terms (SATs). This reformulation can yield stable semidiscretisations for linear advection and Burgers' equation as model problems. In order to enhance these properties, modal filters are introduced to this framework. As a second part of a series, the results of Ranocha, Glaubitz, Öffner, and Sonar ("Enhancing stability of correction procedure via reconstruction using summation-by-parts operators I: Artificial dissipation", 2016) concerning artificial dissipation / spectral viscosity are extended, yielding fully discrete stable schemes. Additionally, a new adaptive strategy to compute the filter strength is introduced and different possible applications of modal filters are compared both theoretically and numerically.

研究动机与目标

• 通过在SBP-SAT框架内引入模态滤波，提升双曲守恒律中校正过程重构（CPR）格式的稳定性。
• 将第一部分中关于人工耗散的稳定性结果扩展至模态滤波，确保全离散格式的稳定性。
• 开发一种新颖的自适应策略，自动选择滤波强度，减少用户输入并提高鲁棒性。
• 比较并分析不同模态滤波应用策略：对解进行滤波、对时间导数进行滤波以及算子分裂方法。
• 在SBP-CPR格式背景下，建立模态滤波与谱黏性的理论与数值等价性。

提出的方法

• 通过使用带同时近似项（SATs）的求和分割（SBP）算子，重新表述CPR方法，以确保线性对流和Burgers方程的保守且稳定的半离散化。
• 通过指数滤波器对解或时间导数应用模态滤波，滤波矩阵记为 $\underline{\underline{F}}$。
• 采用算子分裂方法分析模态滤波与谱黏性的等价性，表明对时间导数进行滤波对应于在由 $\underline{\underline{M}}\underline{\underline{F}}^{-1}$ 诱导的范数下的稳定半离散化。
• 提出一种新的自适应滤波强度策略，以补偿 $(\Delta t)^2$ 误差项，从而实现使用显式欧拉时间积分的可证明稳定性。
• 时间离散化采用显式欧拉格式，滤波后的时间导数表示为 $\partial_t \underline{u} \mapsto \partial_t \underline{\underline{F}} \underline{u}$，从而得到 $\underline{u}_{+} = \underline{u} + \Delta t \, \partial_t \underline{\underline{F}} \underline{u}$。
• 数值实验比较了对解、时间导数以及算子分裂方法应用滤波的效果，采用能量范数评估稳定性。

实验结果

研究问题

• RQ1 模态滤波能否有效应用于SBP-CPR格式，以实现可证明稳定的全离散格式？
• RQ2 在稳定性与精度方面，模态滤波与人工耗散及谱黏性相比如何？
• RQ3 对解进行滤波与对时间导数进行滤波在数值稳定性与振荡控制方面有何影响？
• RQ4 能否设计一种自适应滤波强度策略，根据时间步长与解的行为自动调整滤波强度？
• RQ5 所提出的滤波策略是否能消除人工耗散方法中观察到的时间步长限制？

主要发现

• 模态滤波在时间首阶精度范围内与谱黏性理论等价，为人工耗散提供了一种稳定替代方案。
• 自适应滤波强度策略成功补偿了 $(\Delta t)^2$ 误差项，使得使用显式欧拉时间积分的全离散格式可实现可证明稳定性。
• 对时间导数进行滤波在由 $\underline{\underline{M}}\underline{\underline{F}}^{-1}$ 诱导的范数下可实现稳定半离散化，但完整时间步长中逆滤波的应用会破坏该稳定性，因此该方法劣于对解进行滤波。
• 数值实验表明，对解进行滤波可减少振荡并改善能量行为，而对时间导数进行滤波则导致能量范数振荡，与理论预期一致。
• 对于具有间断初值的线性对流问题，该格式保持稳定且准确，但在存在激波形成的非线性情况下，振荡仍持续存在。
• 所提出的框架通过自适应滤波机制，显著扩展了SBP-CPR格式在非线性问题中的适用性，提升了整体鲁棒性。