[论文解读] Enhancing stability of correction procedure via reconstruction using summation-by-parts operators I: Artificial dissipation

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

Medical Imaging Techniques and Applications参考文献 13被引用 6

一句话总结

该论文通过在求和分割（SBP）和同时近似项（SATs）框架内引入人工耗散，增强了双曲守恒律的校正过程重构（CPR）格式的稳定性。它使用显式欧拉时间离散化证明了全离散稳定性，并提出了一种自适应耗散策略，确保在足够小的时间步长下保持稳定性。

ABSTRACT

The correction procedure via reconstruction (CPR, also known as flux reconstruction) is a framework of high order semidiscretisations used for the numerical solution of hyperbolic conservation laws. Using a reformulation of these schemes relying on summation-by-parts (SBP) operators and simultaneous approximation terms (SATs), artificial dissipation / spectral viscosity operators are investigated in this first part of a series. Semidiscrete stability results for linear advection and Burgers' equation as model problems are extended to fully discrete stability by an explicit Euler method. As second part of this series, Glaubitz, Ranocha, Öffner, and Sonar (Enhancing stability of correction procedure via reconstruction using summation-by-parts operators II: Modal filtering, 2016) investigate connections to modal filters and their application instead of artificial dissipation.

研究动机与目标

• 解决高阶CPR格式在双曲守恒律中缺乏非线性稳定性的缺陷。
• 将半离散SBP-SAT稳定性结果扩展至使用显式时间积分的全离散格式。
• 设计一种自适应人工耗散策略，确保在显式时间推进下保持稳定性。
• 为CPR格式中使用SBP算子的人工耗散建立理论基础。
• 为系列论文的第二部分奠定基础，后者将探讨与模态滤波的联系。

提出的方法

• 通过求和分割（SBP）算子和同时近似项（SATs）重新表述CPR格式，以实现可证明的稳定性。
• 在连续设置中引入人工耗散算子，然后使用SBP算子进行离散化，以保持守恒性和稳定性。
• 在半离散设置中采用耗散算子的特殊表示形式，以维持稳定性并模拟分部积分。
• 应用显式欧拉时间离散化以获得全离散格式。
• 提出一种自适应算法，根据时间步长大小调整人工耗散的强度。
• 将耗散算子的离散特征值与连续勒让德多项式对应物进行比较，以验证稳定性特性。

实验结果

研究问题

• RQ1 能否在SBP-SAT框架中有效整合人工耗散至CPR格式，以确保全离散稳定性？
• RQ2 人工耗散强度的选择如何影响使用显式时间积分的全离散CPR格式的稳定性？
• RQ3 耗散算子的离散特征值与连续勒让德多项式谱之间存在何种关系？
• RQ4 能否设计一种自适应耗散策略，以在显式CPR格式中对足够小的时间步长保持稳定性？
• RQ5 基于SBP的人工耗散算子在如伯吉斯方程等非线性问题中，如何保持守恒性和稳定性？

主要发现

• 所提出的自适应耗散算子在SBP-SAT框架内保持了CPR格式的守恒性和稳定性。
• 线性对流和伯吉斯方程的半离散稳定性结果已通过显式欧拉方法扩展至全离散格式。
• 耗散算子的离散特征值与连续勒让德多项式谱高度一致，验证了算子设计的有效性。
• 自适应耗散策略确保了在时间步长足够小时的全离散稳定性。
• 该方法无需隐式时间积分即可实现稳定性，从而支持高效高阶模拟。
• 该理论框架为系列论文的第二部分奠定了基础，后者将研究与模态滤波的联系。