形式陈述
考虑守恒更新
U i n + 1 = U i n − Δ t Δ x ( f ^ ( U i n , U i + 1 n ) − f ^ ( U i − 1 n , U i n ) ) , 其中数值通量 公理库 数值通量 Numerical flux · Conservative interface flux · 数值界面通量 以左右离散状态给出共享界面通量,并用一致性、单调性与耗散量区分不同守恒格式。 在有界集上连续并一致:
f ^ ( u , u ) = f ( u ) . 令 u h 是由 U i n 形成的时空分片常重构,网格尺度趋零,初值重构在 L loc 1 中趋于 u 0 。Lax–Wendroff 定理的典型标量版本断言:若 u h 一致有界,并且沿该网格族几乎处处收敛(或以足以通过非线性通量极限的方式在 L 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 用离散缺陷、扰动放大和重构误差区分 PDE 格式的一致性、稳定性与收敛性,并明确范数与网格族。 中的其他环节。
证明把离散方程乘采样测试函数 φ i n ,对 i , n 求和并作离散分部求和。时间差转移到 φ t ,相邻共享通量转移到 φ x ;边界和终止项按测试函数支集处理。收敛与有界性使
f ^ ( U i n , U i + 1 n ) ⟶ f ( u ) 在相应积分意义中成立,一致性再把离散恒等式极限化为连续弱形式。
有限体积法 公理库 有限体积法 Finite volume method · Finite-volume discretization · 控制体积法 对每个控制体积分守恒律,以共享界面通量更新单元平均,从离散结构上保证内部通量抵消。 提供这类共享面通量结构的标准实现,但定理本身针对一致守恒离散序列,不要求通量一定来自精确 Riemann 求解器。多维和非结构网格版本还需几何正则、面通量一致性与离散测试函数逼近条件。
直觉
连续弱形式的核心是分部积分,离散守恒式的核心是通量望远镜求和。Lax–Wendroff 定理说明这两种抵消在网格极限中相容:只要离散状态真的收敛、通量对相同左右状态回到物理通量,内部面就不会在极限里留下虚假源项。
它更像“极限身份证检查”,不是“把序列送到终点的发动机”。稳定性和紧性负责从无限多个网格解中找到收敛子列,定理负责核对该子列的极限满足哪条守恒方程,熵不等式再负责从多个弱解中选出唯一可容许者。把三步合成一句“守恒格式必收敛”会遗漏最困难的假设。
例子与边界
对周期线性平流 u t + ( 2 u ) x = 0 ,取 Δ x = 1 , Δ t = 0.2 ,迎风通量 f ^ ( u L , u R ) = 2 u L ,数据
( U 1 n , U 2 n , U 3 n ) = ( 1 , 3 , 2 ) . 周期界面通量从“单元 3 到 1”起为
( f ^ 1 / 2 , f ^ 3 / 2 , 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 。对任意周期权重 φ i ,通量项满足离散分部求和
∑ i φ i ( f ^ i + 1 / 2 − f ^ i − 1 / 2 ) = − ∑ i f ^ i + 1 / 2 ( φ i + 1 − φ i ) . 右侧差商在细化时逼近 φ x ,这就是证明中可复算的代数核心;总量守恒只是取 φ i ≡ 1 的特殊情形。
定理不会排除 Burgers 初值 u L = 0 , u R = 2 的膨胀激波,因为该跳跃确实是弱解。若某个一致守恒但缺少正确熵耗散的序列收敛到它,Lax–Wendroff 识别仍会通过。要得到稀疏扇,必须另外证明离散熵不等式或采用已知熵收敛格式。
仅有弱收敛也可能不足:非线性 f ( u h ) 不一定收敛到 f ( u ) ,高频振荡可能产生 Young measure。边界问题则在离散分部求和后留下边界通量项,极限必须与流入边界熵条件匹配;把整线紧支集证明直接搬到开放边界会漏项。
推论与应用
一条完整收敛论证通常按顺序组织:证明 L ∞ 或不变量区域界,取得 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.