CFL 条件 | Codeboy's Blog

形式陈述

Courant–Friedrichs–Lewy 条件的本质陈述是：若一族离散格式收敛到双曲问题的解，则连续问题在某点的依赖域必须包含在数值格式对该点的渐近依赖域内。它是收敛的必要几何条件，不是所有格式的统一充分稳定性定理。

对均匀网格上的常速平流

u_t+a u_x=0,

若显式一步 stencil 在 $a>0$ 时只连接 $U_i^n$ 与 $U_{i-1}^n$ ，一时间步的数值依赖域向上游延伸 $\Delta x$ ，连续特征脚却位于 $x_i-a\Delta t$ 。必要包含关系给

0\le \frac{a\Delta t}{\Delta x}\le1.

绝对值形式 $|a|\Delta t/\Delta x\le1$ 适用于左右对称陈述，但实际上游 stencil 还必须随 $a$ 的符号改变。

对非线性标量守恒律，局部波速为 $f'(u)$ ，常用充分操作限制为

\frac{\Delta t}{\Delta x} \max_{u\in\mathcal R^n}|f'(u)| \le C_{\mathrm{CFL}},

其中 $\mathcal R^n$ 覆盖本步可能状态，常数由通量、stencil 和时间积分决定。双曲系统以 Jacobian $DF(U)$ 的最大特征速度替换 $|f'(u)|$ 。非均匀或多维有限体积网格需按单元体积和各面谱半径形成局部限制，不能只用全局最大边长。

CFL 只约束网格族的传播尺度。中心空间配前向 Euler 的平流格式即使满足 $|a|\Delta t/\Delta x<1$ ，仍有放大因子模大于一；隐式格式可能无同样的稳定上限，却仍需时间分辨率保证精度。把“满足 CFL”直接写成“已证明收敛”会混淆必要条件、稳定估计和误差分析。

直觉

离散算法每一步只能从 stencil 访问有限邻居，许多步后这些访问形成一个数值光锥。真实双曲方程的信息沿特征光锥传播。若真实信息在一步内越过算法能看到的最远邻居，两个不同初值可能对数值输入完全相同，却对真实解产生不同影响；任何局部更新都不可能同时逼近二者。

Courant 数衡量一时间步内波移动了多少个网格间隔。它不是物理速度本身，也不是只由 PDE 决定的常数：改变 stencil 宽度、时间积分、维数或局部网格都会改变允许上界。显式算法通常用减小时间步换取依赖域覆盖，代价是更多步数。

例子与边界

取 $a=2$ 、 $\Delta x=0.1$ 的右行平流，并使用一步只看一个左邻的迎风格式。CFL 条件给

\Delta t\le\frac{0.1}{2}=0.05.

若 $\Delta t=0.04$ ，一时间步连续特征移动 $a\Delta t=0.08$ ，在数值可见的 $0.1$ 上游区间内；三步后真实脚位移 $0.24$ ，数值依赖域已扩展到 $0.3$ ，包含关系仍成立。若 $\Delta t=0.06$ ，真实一步移动 $0.12$ ，超过一个单元，数值值尚未读取到真正影响它的数据。

依赖域包含仍不保证稳定。前向时间、中心空间格式

U_i^{n+1}=U_i^n-\frac C2(U_{i+1}^n-U_{i-1}^n)

的 stencil 向两侧各延伸一格， $|C|\le1$ 时满足几何覆盖；但 Fourier 因子

G(\theta)=1-iC\sin\theta,\qquad |G(\theta)|^2=1+C^2\sin^2\theta

对大多数非零频率超过一，固定非零 $C$ 的细化族不稳定。这个反例把 CFL 必要性与 von Neumann 稳定性清楚分开。

扩散方程的显式限制常按 $\Delta t=O(\Delta x^2)$ ，来源是离散扩散算子的谱而非有限速双曲依赖域；习惯上也有人称其为“扩散 CFL”，但不能据此说热方程有速度 $\Delta x/\Delta t$ 的有限特征。局部加密网格则由最小单元或小 cut cell 控制全局显式步长，除非采用局部时间步或隐式处理。

推论与应用

CFL 数用于把 Riemann 求解器报告的最大波速转成时间步。稳健实现应在每步状态更新后重新估计最大谱半径，以

\Delta t=C_{\mathrm{target}} \min_i\frac{\Delta x_i}{a_{i,\max}}

或相应多维面公式选步，并为零波速、非有限速度和极小单元显式处理。把上一步的速度上界永久缓存，强非线性波出现时可能突然越界。

验证应同时画出理论 CFL 上界与实测放大，而不是只记录“程序没有崩”。满足条件后仍需证明数值通量一致、格式稳定并建立紧性或误差估计；违反条件时偶然一组平滑初值没有爆炸，也不能推翻依赖域必要性。CFL 是连接连续传播和离散 stencil 的接口，不是替代全部数值分析的认证章。

参考资料

R. Courant, K. Friedrichs, and H. Lewy, “Über die partiellen Differenzengleichungen der mathematischen Physik,” Mathematische Annalen 100, 1928, pp. 32–74.
Randall J. LeVeque, Finite Volume Methods for Hyperbolic Problems, Cambridge University Press, 2002, Ch. 4.
John C. Strikwerda, Finite Difference Schemes and Partial Differential Equations, 2nd ed., SIAM, 2004, Ch. 2.

形式陈述 ​

反向链接

形式陈述