von Neumann 稳定性分析

形式陈述 ​

von Neumann 分析研究常系数、线性、平移不变的规则网格差分格式公理库偏微分方程有限差分法Finite-difference method for PDE · Finite-difference PDE discretization · Difference stencil在整张网格上用差分 stencil 替换偏微分算子，并连同边界闭合组装为稀疏离散系统。。先考虑一维无限网格或周期网格上的标量一步更新

$$U_j^{n+1}=\sum _{\ell =-p}^q{a}_{\ell }{U}_{j+\ell }^n,$$

其中系数可以依赖无量纲网格参数，但不依赖 $j$。线性格式的误差满足同一齐次更新，所以只需分析误差传播。

离散纯频率写成

$$e_j^n={\hat{e}}^n{e}^{ij\theta },\theta \in [-\pi ,\pi ].$$

$\theta$ 是无量纲网格频率；$\theta$ 与 $\theta +2\pi$ 在节点上不可区分。将模态代入格式得到

$${\hat{e}}^{n+1}=G(\theta ){\hat{e}}^n,G(\theta )=\sum _{\ell =-p}^q{a}_{\ell }{e}^{i\ell \theta }.$$

$G$ 称为放大因子或格式符号。平移不变使不同频率彼此不耦合，这正是Fourier 分解公理库Fourier 变换Fourier transform · 傅里叶变换把非周期函数分解为连续频率成分，并将卷积和平移不变算子转为频域乘法。在此处能够对角化更新的原因。

对无限网格上的序列，定义与 ${\ell }_h^2$ 范数匹配的离散 Fourier 变换

$${\hat{e}}^{,n}(\theta )=\sqrt{\frac{h}{2\pi }}\sum _{j\in \mathbb{Z}}{e}_j^n{e}^{-ij\theta }.$$

于是 Plancherel 等式给出

$$\Vert e^n{\Vert }_{{\ell }_h^2}^2=h\sum _j|e_j^n{|}^2={\int }_{-\pi }^{\pi }|G(\theta ){|}^{2n}|{\hat{e}}^0(\theta ){|}^{2}d\theta$$

（周期网格时积分替换为离散频率求和）。因此在这一特定框架中，

$$\underset{\theta \in [-\pi ,\pi ]}{sup}|G(\theta )|\le 1$$

足以保证 ${\ell }_h^2$ 范数不增长，并且也是相应收缩结论的必要条件。这里的“充要”依赖标量、常系数、平移不变、周期或无限网格以及离散 $L^2$ 范数，不能脱离这些条件移植。

有限时间稳定不一定要求每步严格收缩。若存在与网格无关的 $C$，使

$$\underset{\theta }{sup}|G(\theta )|\le 1+C\tau ,$$

则对 $n\tau \le T$ 有

$$\underset{\theta }{sup}|G(\theta ){|}^n\le (1+C\tau {)}^{T/\tau }\le e^{CT}.$$

反之，若某个可达频率每步放大至少 $1+c$，其中 $c>0$ 不随 $\tau$ 缩小，那么约 $T/\tau$ 步会产生随网格加密发散的增长。判断必须针对整个网格族，而不是只看一个固定离散矩阵。

向量值系统会得到矩阵符号

$${\hat{e}}^{n+1}(\theta )=\mathsf{G}(\theta ){\hat{e}}^n(\theta ).$$

正确的有限时间条件是

$$\underset{\theta }{sup}\underset{0\le n\tau \le T}{sup}\Vert \mathsf{G}(\theta {)}^n\Vert \le K_T,$$

并要求常数与网格无关。多步格式可以把历史分量组成 companion state，再检查相应放大矩阵的幂。即使验证了 $\rho (\mathsf{G}(\theta ))\le 1$，一般也不能推出这项幂界：非正规矩阵会在特征值不变时产生显著瞬态，单位圆上的非平凡 Jordan 块还会导致多项式增长。

von Neumann 分析是格式稳定性的一个工具，不是稳定性的另一个同义定义。它把一般算子幂界化成频率逐点计算；当 Fourier 模态不再对角化完整更新时，仍应回到离散能量、矩阵范数或边界算子估计。

直觉 ​

每个网格误差都可以分解成从缓慢起伏到逐点交替的一组离散波。常系数 stencil 对纯波只改变复振幅：$|G|$ 决定幅值增减，$\arg G$ 决定相位推进。逐个频率检查后再取最坏值，便能找出肉眼平滑初值中很小、却可能被格式迅速放大的隐藏网格振荡。

高频端 $\theta \approx \pi$ 尤其重要。它代表相邻节点近乎交替变号的模式，连续图像可能没有明显对应，却常最先触碰显式格式的稳定边界。只测试低频正弦或一条光滑解，容易漏掉真正控制 CFL 的模式。

例子与边界 ​

考虑热方程

$$u_t=u_{xx}$$

的前向时间、中心空间格式

$$U_j^{n+1}=U_j^n+r(U_{j+1}^n-2U_j^n+U_{j-1}^n),r=\frac{\tau }{h^2}.$$

对模态 $e^{ij\theta }$，二阶差分的符号为

$$e^{i\theta }-2+e^{-i\theta }=-4{\sin }^2\frac{\theta }{2},$$

所以

$$G(\theta )=1-4r{\sin }^2\frac{\theta }{2}.$$

当 $r\ge 0$ 时，$G$ 为实数且最大值为 $1$。最小值出现在 $\theta =\pi$，等于 $1-4r$；要求所有频率满足 $|G(\theta )|\le 1$，恰好得到

$$0\le r\le \frac{1}{2},\text{即}\tau \le \frac{h^2}{2}.$$

$r=1/2$ 时最高频模态的 $G=-1$，幅值不变但每步翻转符号；$0<r<1/2$ 才对除零频外的模式产生严格衰减。这个限制属于所选显式时间—空间组合，不是热方程连续解本身要求按 $h^2$ 采样。

再看 $a>0$ 的平流方程

$$u_t+a{u}_x=0.$$

迎风格式为

$$U_j^{n+1}=U_j^n-\nu (U_j^n-U_{j-1}^n),\nu =\frac{a\tau }{h}.$$

它的放大因子与模平方为

$$G(\theta )=1-\nu +\nu e^{-i\theta },$$$$|G(\theta ){|}^2=1-4\nu (1-\nu ){\sin }^2\frac{\theta }{2}.$$

因此全部频率不增长当且仅当

$$0\le \nu \le 1,\text{即}0\le a\tau \le h.$$

这条 CFL 条件有清楚的传播图像：一个时间步内连续特征移动的距离不能越过一个上风网格间隔。中心空间的前向时间格式虽局部一致，却在固定非零 $\nu$ 的网格族上放大某些频率；该反例由PDE 三分框架公理库PDE 离散的一致性、稳定性与收敛性Consistency stability and convergence for PDE · PDE scheme stability · PDE convergence framework用离散缺陷、扰动放大和重构误差区分 PDE 格式的一致性、稳定性与收敛性，并明确范数与网格族。承担，不在本页重复推导。

对空间半离散系统

$$\frac{dU}{dt}=L_{h}U,$$

若时间积分法的稳定函数公理库ODE 方法的绝对稳定域Absolute stability region · Region of absolute stability · A-stability以标量衰减测试方程的放大因子定义一步法稳定域，并区分 A-stability、L-stability、精度与其他稳定概念。为 $R(z)$，完整一步是 $U^{n+1}=R(\tau L_h)U^n$。当 $L_h$ 也是标量平移不变算子，设其 Fourier 符号为 ${\lambda }_h(\theta )$，就有

$$G(\theta )=R(\tau {\lambda }_h(\theta )).$$

于是检查所有 $\tau {\lambda }_h(\theta )$ 是否落入 ODE 绝对稳定域，可以重现这一标量正规情形的 von Neumann 条件。ODE 稳定域只描述时间方法的测试方程；空间符号、允许频率、边界和所选范数仍由 PDE 离散提供，两种分析不能合并成一个无上下文的“稳定”标签。

可变系数会在 Fourier 空间中卷积不同频率，物理边界会破坏循环平移不变性。冻结系数或只分析内部 stencil 往往给出必要诊断，却未必控制边界闭合引起的模态；此时应使用能量法、正常模分析或完整更新矩阵的统一幂界。

对非正规矩阵符号，所有特征值位于单位圆内仍可能有大瞬态。例如

$$\mathsf{G}=\left(\begin{array}{cc}1& 1\\ 0& 1\end{array}\right),{\mathsf{G}}^n=\left(\begin{array}{cc}1& n\\ 0& 1\end{array}\right).$$

两个特征值都等于 $1$，算子幂却线性增长。若矩阵随网格变化，特征向量条件数也可能随 $h$ 恶化；只画谱半径会漏掉这项不一致放大。

非线性格式没有固定放大因子。围绕常状态线性化或冻结系数，可以给出小扰动的必要信息，但有限振幅波、激波、限制器切换和熵条件不由线性符号完整决定。把线性化分析通过不代表原非线性更新全局稳定。

推论与应用 ​

von Neumann 计算应报告 PDE、完整内部 stencil、边界假设、无量纲参数、频率区间、所用范数以及最终取到最坏值的频率。只写一个 $G(\theta )$ 而不解释它来自哪个网格族，无法判断条件是必要、充分还是仅供诊断。

热方程的 $r\le 1/2$ 展示空间最高频率随 $h^{-2}$ 增长如何压缩显式时间步，平流迎风格式的 $\nu \le 1$ 则把离散传播速度与网格间隔连接起来。两种限制的尺度不同，来自各自空间算子符号；不能把所有 CFL 条件都概括成同一个固定常数。

参考资料

• John C. Strikwerda, Finite Difference Schemes and Partial Differential Equations, 2nd ed., SIAM, 2004, Chs. 2 and 4.
• Lloyd N. Trefethen, Finite Difference and Spectral Methods for Ordinary and Partial Differential Equations, unpublished text, 1996, Ch. 4.
• Randall J. LeVeque, Finite Difference Methods for Ordinary and Partial Differential Equations, SIAM, 2007, Chs. 9–10.
• MIT OpenCourseWare, 18.336 Numerical Methods for Partial Differential Equations, stability and Fourier analysis notes.