守恒律的 Lax–Wendroff 定理

形式陈述 ​

考虑守恒更新

$$U_i^{n+1}=U_i^n-\frac{\mathrm{\Delta }t}{\mathrm{\Delta }x}(\hat{f}(U_i^n,U_{i+1}^n)-\hat{f}(U_{i-1}^n,U_i^n)),$$

其中数值通量公理库数值通量Numerical flux · Conservative interface flux · 数值界面通量以左右离散状态给出共享界面通量，并用一致性、单调性与耗散量区分不同守恒格式。在有界集上连续并一致：

$$\hat{f}(u,u)=f(u).$$

令 $u_h$ 是由 $U_i^n$ 形成的时空分片常重构，网格尺度趋零，初值重构在 $L_{\mathrm{loc}}^1$ 中趋于 $u_0$。Lax–Wendroff 定理的典型标量版本断言：若 $u_h$ 一致有界，并且沿该网格族几乎处处收敛（或以足以通过非线性通量极限的方式在 $L_{\mathrm{loc}}^1$ 收敛）到 $u$，则 $u$ 是守恒律弱解公理库守恒律的弱解Weak solution of a conservation law · Distributional solution of conservation law · 守恒律分布解通过测试函数积分恒等式定义允许跳跃的守恒律解，并保留初值、边界与跨间断通量平衡。。

关键量词是“如果网格解已经收敛”。定理不从一致性推出收敛，不提供网格无关稳定界，也不判断极限是否为熵解。紧性可由单调性、BV/TVD、熵稳定或问题特定估计另行建立；这些属于一致性—稳定性—收敛性框架公理库PDE 离散的一致性、稳定性与收敛性Consistency stability and convergence for PDE · PDE scheme stability · PDE convergence framework用同一误差递推连接局部缺陷、扰动传播与全局误差，澄清归一化、网格范数和线性等价定理的条件。中的其他环节。

证明把离散方程乘采样测试函数 ${\phi }_i^n$，对 $i,n$ 求和并作离散分部求和。时间差转移到 ${\phi }_t$，相邻共享通量转移到 ${\phi }_x$；边界和终止项按测试函数支集处理。收敛与有界性使

$$\hat{f}(U_i^n,U_{i+1}^n)⟶f(u)$$

在相应积分意义中成立，一致性再把离散恒等式极限化为连续弱形式。

有限体积法公理库有限体积法Finite volume method · Finite-volume discretization · 控制体积法对每个控制体积分守恒律，以共享界面通量更新单元平均，从离散结构上保证内部通量抵消。提供这类共享面通量结构的标准实现，但定理本身针对一致守恒离散序列，不要求通量一定来自精确 Riemann 求解器。多维和非结构网格版本还需几何正则、面通量一致性与离散测试函数逼近条件。

直觉

连续弱形式的核心是分部积分，离散守恒式的核心是通量望远镜求和。Lax–Wendroff 定理说明这两种抵消在网格极限中相容：只要离散状态真的收敛、通量对相同左右状态回到物理通量，内部面就不会在极限里留下虚假源项。

它更像“极限身份证检查”，不是“把序列送到终点的发动机”。稳定性和紧性负责从无限多个网格解中找到收敛子列，定理负责核对该子列的极限满足哪条守恒方程，熵不等式再负责从多个弱解中选出唯一可容许者。把三步合成一句“守恒格式必收敛”会遗漏最困难的假设。

例子与边界

对周期线性平流 $u_t+(2u{)}_x=0$，取 $\mathrm{\Delta }x=1,\mathrm{\Delta }t=0.2$，迎风通量 $\hat{f}(u_L,u_R)=2u_L$，数据

$$(U_1^n,U_2^n,U_3^n)=(1,3,2).$$

周期界面通量从“单元 3 到 1”起为

$$({\hat{f}}_{1/2},{\hat{f}}_{3/2},{\hat{f}}_{5/2})=(4,2,6).$$

更新给

$$(U_1^{n+1},U_2^{n+1},U_3^{n+1})=(1.4,2.2,2.4),$$

总和仍为 $6$。对任意周期权重 ${\phi }_i$，通量项满足离散分部求和

$$\sum _i{\phi }_i({\hat{f}}_{i+1/2}-{\hat{f}}_{i-1/2})=-\sum _i{\hat{f}}_{i+1/2}({\phi }_{i+1}-{\phi }_i).$$

右侧差商在细化时逼近 ${\phi }_x$，这就是证明中可复算的代数核心；总量守恒只是取 ${\phi }_i\equiv 1$ 的特殊情形。

定理不会排除 Burgers 初值 $u_L=0,u_R=2$ 的膨胀激波，因为该跳跃确实是弱解。若某个一致守恒但缺少正确熵耗散的序列收敛到它，Lax–Wendroff 识别仍会通过。要得到稀疏扇，必须另外证明离散熵不等式或采用已知熵收敛格式。

仅有弱收敛也可能不足：非线性 $f(u_h)$ 不一定收敛到 $f(u)$，高频振荡可能产生 Young measure。边界问题则在离散分部求和后留下边界通量项，极限必须与流入边界熵条件匹配；把整线紧支集证明直接搬到开放边界会漏项。

推论与应用

一条完整收敛论证通常按顺序组织：证明 $L^{\mathrm{\infty }}$ 或不变量区域界，取得 BV/TVD 或补偿紧性，抽取强收敛子列，用本定理得到弱解，再用离散熵不等式识别唯一熵解。唯一性随后把“某个子列”升级为整个网格族收敛。每一步失败都对应不同诊断，不能由更细网格图像替代。

定理还解释了为何共享边界通量不可妥协。若相邻单元对同一面使用两个不同值，离散分部求和会留下内部源项，极限可能不再满足原守恒律。审计格式时应分别检查通量一致性、全局守恒、强收敛证据和熵残差；稳定、一致与收敛相互关联，却不是可交换的同义词。

参考资料

• Peter D. Lax and Burton Wendroff, “Systems of Conservation Laws,” Communications on Pure and Applied Mathematics 13(2), 1960, pp. 217–237.
• Randall J. LeVeque, Finite Volume Methods for Hyperbolic Problems, Cambridge University Press, 2002, §12.10.
• E. Godlewski and P.-A. Raviart, Numerical Approximation of Hyperbolic Systems of Conservation Laws, Springer, 1996, Ch. 3.