冻结系数法及其与冯诺依曼稳定性分析的关系

计算科学 pde 有限差分 稳定
2021-12-04 21:43:04

我正在考虑两个方程

ut=a(x)uxx
vt=b(x)vx
作为抛物线和双曲线方程组的经典代表。如果a(x)=ab(x)=b是常数,为了显示我用于这些方程的任何有限差分方案的稳定性,我可以使用傅里叶或所谓的冯诺依曼稳定性分析,即必须计算放大因子的大小。

但是,在可变系数的情况下,其中a=a(x)b=b(x),我有些困惑。从 Strikwerda 的书(有限差分方案和偏微分方程)第 59 页,它说可以使用冻结系数的方法,即冻结每个网格点的值,使用傅里叶分析并推断变量系数的稳定性。另一方面,在 Ascher (进化微分方程的数值方法)一书中,第 153 页它说即使该方法在常数系数下是稳定的,那么它并不一定意味着变量系数的稳定性。

因此,我看到了两本书之间的矛盾。我的问题是我看到了冯诺依曼分析是如何工作的,我想知道我是否也可以将它与可变系数 PDE 一起使用。如果可以的话,我正在查看的两个方程是否相同,或者两者在冻结系数的方法方面存在一些差异?基本上,冯诺依曼稳定性分析对抛物线和双曲线方程的限制在哪里?

1个回答

一位作者务实,而另一位作者严谨。简短的回答是,如果系数和解足够平滑,则冻结系数问题的稳定性被证明可以保证可变系数问题的稳定性。

如果您想知道需要多少光滑度,请准备好弄脏您的手。抛物线侧的基本参考是

弗里茨,约翰。 用差分法积分抛物方程: I. 无限区间的线性和拟线性方程。纯粹和应用数学通讯,5(2):155--211 (1952)。

对于双曲线问题,请参见

斯特朗,吉尔伯特。 准确的偏差法Ⅱ数学数学 6, 37-46 (1964)。

附加说明:Ascher 声称他的示例 5.1 和 5.6 展示了在光滑情况下冯诺依曼稳定的方案如何在非光滑情况下不稳定。但这两个例子都具有误导性。例 5.1 涉及 KdV 方程,实际上一直有一个平滑解。同时,例 5.6 仍然是一个线性常数系数问题;只有解是不连续的。那里出现的振荡是由于数值分散而不是不稳定性。所以这两个例子并没有说明可能出现的病态。不幸的是,我不知道有一个这样做。