如何分析边界处 PDE 离散化的稳定性?

计算科学 pde 有限差分 稳定
2021-12-20 00:28:34

我有一个偏微分方程的数值离散化,该方程在边界点似乎不稳定或稳定,具体取决于我使用的有限差分方案。是否有调查边界稳定性的标准方法?

我只知道冯诺依曼稳定性分析,但这似乎考虑了整个区域,我对来自边界的不稳定性特别感兴趣。

1个回答

如果您的 PDE 和离散化是线性的,那么您的数值方法可以写为

Un+1=AUn

在哪里A是一个方阵。研究稳定性的一种方法是查看一些诱导矩阵范数A:如果大于1,离散化不稳定。这是完全通用的,可能是您正在寻找的。下面我解释一下冯诺依曼分析只是这种一般方法的一个特例。

如果A是一个正规矩阵(即,如果AAT=ATA),然后为了显示稳定性,检查的特征值就足够了A模数不大于 1。这是因为对于正态矩阵,诱导的 2 范数等于最大特征值的模。

如果问题具有周期性边界条件并使用标准有限差分离散化,则A将是一个循环矩阵。循环矩阵是正规的,循环矩阵的特征向量只是离散的傅里叶模式。这使得计算它们的特征值变得容易,而冯诺依曼分析只是一种确定这些特征值的方法。