冯·诺依曼稳定性分析
概述
冯·诺依曼稳定性分析（von Neumann stability analysis）由约翰·冯·诺依曼于 1940 年代在曼哈顿计划数值计算工作中发展。其核心思想是将有限差分格式的误差展开为傅里叶级数，分析每个傅里叶模式经一个时间步后的放大倍数（amplification factor）$g(\xi)$。稳定性等价于所有模式的放大因子满足 $|g(\xi)| \leq 1$。这一方法提供了比CFL条件更精确的稳定性判据，将定性的收敛条件转化为可精确计算的不等式。
关键内容
方法框架
考虑线性有限差分格式，将误差 $\epsilon_j^n = u_j^n - u(x_j, t_n)$ 展开为傅里叶模式：
$$\epsilon_j^n = \sum_\xi G^n(\xi) e^{i\xi j\Delta x}$$
每个模式 $e^{i\xi j\Delta x}$ 经一步时间推进后变为 $g(\xi) \cdot e^{i\xi j\Delta x}$，其中 $g(\xi)$ 是放大因子（amplification factor）。
稳定性判据：格式稳定 ⟺ 对所有波数 $\xi$，$|g(\xi)| \leq 1$（Fourier 模式不被放大）。
应用示例：一维波动方程显式格式
对格式 $u_j^{n+1} = u_j^n - \frac{c\Delta t}{2\Delta x}(u_{j+1}^n - u_{j-1}^n)$（一阶迎风格式），放大因子：
$$g(\xi) = 1 - \frac{c\Delta t}{\Delta x}(e^{i\xi\Delta x} - e^{-i\xi\Delta x})/2 = 1 - ic\nu\sin(\xi\Delta x)$$
$|g(\xi)|^2 = 1 + \nu^2\sin^2(\xi\Delta x)$，当 $\nu \neq 0$ 时恒 $> 1$——无条件不稳定！
对Lax-Friedrichs 格式（Lax 1954年修改为稳定格式），放大因子在 $\nu = c\Delta t/\Delta x \leq 1$ 时满足 $|g(\xi)| \leq 1$——这精确验证了CFL条件。
与 CFL 条件的关系
| CFL 条件（1928） | 冯·诺依曼分析（1940s） | |
|---|---|---|
| 来源 | 依赖域几何论证（物理直觉） | 傅里叶模式代数计算（定量工具） |
| 输出 | 收敛必要条件 | 精确稳定性充要条件（对线性格式） |
| 适用范围 | 显式格式，双曲方程 | 线性格式，任意类型方程 |
| 局限性 | 必要非充分 | 仅适用线性、常系数格式 |
对许多格式，两者条件完全一致；冯·诺依曼分析提供了 CFL 条件的数学"证明"。
Lax 等价定理（1956）
Peter Lax 和 R. D. Richtmyer 将冯·诺依曼分析的稳定性概念与收敛性正式联系起来：
Lax 等价定理：对于适定的线性初值问题的一致差分格式，稳定性 ⟺ 收敛性。
这将冯·诺依曼分析的稳定性判据升级为收敛性的充要工具。
扩散方程 FTCS 格式示例
一维扩散方程 $u_t = \nu u_{xx}$，前向时间中心空间（FTCS）格式：
$$g(\xi) = 1 - 4r\sin^2!\left(\frac{\xi\Delta x}{2}\right), \quad r = \frac{\nu\Delta t}{(\Delta x)^2}$$
稳定条件 $|g(\xi)| \leq 1$ 要求 $r \leq 1/2$，即：
$$\Delta t \leq \frac{(\Delta x)^2}{2\nu}$$
这是经典的抛物型方程时间步长限制：空间步长减半则时间步长须减至四分之一（比波动方程的 CFL 更严格）。
矩阵计算误差分析（1947 论文）
冯·诺依曼与戈尔茨坦（Herman Goldstine）的 1947 年矩阵求逆论文是另一重要贡献，奠定了数值线性代数误差分析基础：
- 矩阵条件数（condition number）概念首次系统提出：$\text{cond}(A) = |\lambda_{\max}|/|\lambda_{\min}|$（正规矩阵情形），衡量问题对输入扰动的敏感性
- 误差界定理：$t$ 位浮点精度下矩阵求逆的相对误差 $\sim \text{cond}(A) \cdot 2^{-t}$
- 统计误差模型：将舍入误差建模为随机变量，推导期望误差——此方法被后来的后向误差分析（Wilkinson, 1963）的确定性方法部分取代，但在现代概率数值分析中重新受到重视
局限性
- 仅适用线性格式：非线性方程（如 Navier-Stokes）需要线性化（"冻结系数"）或其他分析工具
- 常系数假设：变系数格式需要局部分析（冻结系数启发式方法）
- 不处理边界：傅里叶分析假设周期边界，实际边界条件的稳定性影响需要 GKS 理论（Gustafsson-Kreiss-Sundstrom）等更复杂工具
- 单步增长因子：不直接处理长时积分中相位误差的累积问题，需修正方程（modified equation）分析补充
- 矩阵误差：统计舍入误差假设比 Wilkinson 确定性后向误差分析更粗糙
历史背景
冯·诺依曼的稳定性分析主要形成于 1944–1947 年洛斯阿拉莫斯国家实验室内部报告，源于曼哈顿计划中流体力学与辐射传输方程的数值模拟需求——差分格式不稳定时误差指数增长，结果完全无意义。该方法通过口头交流迅速传播，与理查德·柯朗等人合作的冯·诺依曼1950年论文（Charney, Fjortoft, von Neumann）将其应用于首次成功的数值天气预报。
来源
- raw/books/数值分析/11_cfl_condition.md
- raw/books/数值分析/13_von_neumann_stability.md