[论文解读] Stability of explicit numerical schemes for convection-dominated problems

Erwan Deriaz|arXiv (Cornell University)|Jul 20, 2010

Advanced Numerical Methods in Computational Mathematics参考文献 12被引用 2

一句话总结

本文推导出对流主导问题中显式数值格式的新稳定性条件，表明时间步长必须满足 $\delta t \leq C\delta x^{2r/(2r-1)}$（$r$ 为整数），该条件强于经典的CFL条件。该条件通过数值验证，并在解光滑的假设下证明对非线性方程也成立，为高阶格式提供了更优的稳定性边界。

ABSTRACT

This paper presents original and close to optimal stability conditions linking the time step and the space step, stronger than the CFL criterion: $\delta t\leq C\delta x^\alpha$ with $\alpha=\frac{2r}{2r-1}$, $r$ an integer, for some numerical schemes we produce, when solving convection-dominated problems. We test this condition numerically and prove that it applies to nonlinear equations under smoothness assumptions.

研究动机与目标

• 建立比经典CFL条件更严格的显式时间积分在对流主导问题中的稳定性约束。
• 分析高阶数值格式中时间步长 $\delta t$ 与空间步长 $\delta x$ 之间的相互作用。
• 在解光滑的假设下，将稳定性结果从线性方程推广至非线性对流主导方程。
• 通过数值实验验证理论稳定性条件。

提出的方法

• 推导显式格式的稳定性条件 $\delta t \leq C\delta x^{2r/(2r-1)}$，其中 $r$ 为表示精度阶数的整数。
• 利用能量估计和傅里叶分析，分析格式中数值误差放大的因子。
• 构造特定的数值格式，以实现所推导的稳定性边界。
• 将稳定性条件应用于具有光滑初值的线性和非线性对流主导方程。
• 通过数值模拟验证，当 $\delta t$ 满足所提条件时，格式保持稳定。
• 在解满足适当光滑性假设的条件下，证明该稳定性条件对非线性方程也成立。

实验结果

研究问题

• RQ1 在对流主导问题中，显式格式的 $\delta t$ 与 $\delta x$ 之间最严格的可能稳定性条件是什么？
• RQ2 所提出的稳定性条件在严格性和适用性方面与经典CFL准则相比如何？
• RQ3 所推导的稳定性边界能否推广至非线性对流主导方程？
• RQ4 所提条件是否在实践中确保数值稳定性，如模拟结果所证实？
• RQ5 整数参数 $r$ 在决定稳定性条件中指数 $\alpha = 2r/(2r-1)$ 的作用是什么？

主要发现

• 本文建立的稳定性条件 $\delta t \leq C\delta x^{2r/(2r-1)}$ 对对流主导问题而言，严格强于经典CFL条件。
• 所推导的条件对所提出的数值格式而言接近最优，表明其在理论上具有紧密的稳定性边界。
• 数值实验表明，当时间步长遵循所提条件时，格式保持稳定。
• 在解光滑的假设下，该稳定性结果可推广至非线性对流主导方程。
• 指数 $\alpha = 2r/(2r-1)$ 随 $r$ 增大而增加，当 $r \to \infty$ 时趋近于1，表明高阶格式具有更好的稳定性。