[论文解读] On the stability of Forward in Time and Centred in Space (FTCS) scheme for scalar hyperbolic equation

Ritesh Kumar Dubey|arXiv (Cornell University)|Mar 6, 2015

Computational Fluid Dynamics and Aerodynamics参考文献 16被引用 4

一句话总结

本文挑战了传统观点，即对于标量双曲守恒律，时间向前空间中心（FTCS）格式在所有情况下均不稳定。通过引入基于连续空间梯度比（θ）的数据相关稳定性准则，作者证明：当初始数据分布位于由CFL数和波速定义的特定区域内时，FTCS是条件稳定且非振荡的。主要贡献在于提出了一套新的分析框架，用于识别依赖于数据的稳定性边界，并在线性和非线性情形下通过数值实验得到验证。

ABSTRACT

It is well known that Forward Time and Centred in Space (FTCS) scheme for scalar Hyperbolic Conservation Law (HCL) is unconditionally unstable. The main contribution of this work to show that FTCS is conditionally stable for HCL. A new approach is used to give bounds on the initial data profile by transforming FTCS into two point convex combination scheme. Numerical results are given in support of the claim.

研究动机与目标

• 挑战长期以来认为FTCS格式在标量双曲守恒律中无条件不稳定的信念。
• 识别尽管已知存在Von-Neumann不稳定性，FTCS格式仍可避免数值振荡的条件。
• 基于连续空间梯度比（θ）提出一种数据依赖的稳定性准则，以定义FTCS保持非振荡的区域。
• 通过线性和非线性双曲问题的数值实验验证理论边界。
• 通过引入初始数据的光滑性和梯度结构，将稳定性理解从传统的CFL准则扩展至更广泛范围。

提出的方法

• 利用梯度比θ将FTCS格式转化为两点凸组合形式，从而通过凸性约束实现稳定性分析。
• 将稳定性区域S_FTCS定义为θ空间中的一组区间，其取值取决于波速a的符号和CFL数C。
• 推导出保证凸组合系数非负且总和为一的θ解析边界，从而确保无虚假振荡。
• 采用混合数值策略：当θ位于理论稳定性区域S_FTCS内时切换至FTCS格式，否则使用一阶迎风格式。
• 通过正弦、间断和光滑初始条件的数值模拟，测试在不同CFL数下稳定性预测的有效性。
• 通过线性化通量，将分析扩展至非线性标量守恒律（如Burgers方程），并应用相同的θ准则。

实验结果

研究问题

• RQ1 尽管Von-Neumann分析已表明FTCS格式无条件不稳定，它是否可能在标量双曲方程中实现条件稳定？
• RQ2 初始数据分布——特别是连续空间梯度比——在决定FTCS中数值振荡发生时机方面起什么作用？
• RQ3 解空间中是否存在一个数据依赖的区域，使得即使满足CFL条件，FTCS仍保持非振荡？
• RQ4 不同类型的初始数据（光滑与间断）如何影响FTCS格式的稳定性与振荡行为？
• RQ5 θ的理论稳定性边界是否可在数值上验证，并推广至非线性双曲守恒律？

主要发现

• 当初始数据梯度比θ位于依赖于CFL数和波速符号的集合S_FTCS内时，FTCS格式呈现条件稳定且非振荡特性。
• 对于正弦等光滑初始数据，随着CFL数C趋近于零，数值振荡逐渐减弱，与理论预测一致：振荡区域随之缩小。
• 在间断初始数据下，FTCS立即引发强烈振荡，证实此类数据位于稳定区域S_FTCS之外。
• 采用混合FTCSUP格式（仅在S_FTCS内使用FTCS）的数值结果成功抑制了振荡，而FTUPCS格式（仅在S_FTCS外使用FTCS）则表现出虚假振荡，验证了理论边界的有效性。
• 该稳定性准则可推广至非线性标量守恒律（如Burgers方程），其中相同的θ区域可预测无振荡行为，前提是CFL数满足λ·max|g′(u)| ≤ 1。
• 当C → 0时，理论稳定性区域S_FTCS退化为(-1, 0)，与光滑数据下小时间步长可稳定格式的观察一致。