Lax-Wendroff格式在无粘Burgers方程间断初值下的二阶收敛性问询

咱们来逐个拆解你的问题：

问题1：Lax-Wendroff格式能否达到二阶收敛阶？ #

答案是不能。原因很直接：Lax-Wendroff是二阶线性守恒格式，在处理光滑解时确实能达到二阶精度，但无粘Burgers方程带间断初始条件时，解会演化出激波这类间断结构。在间断附近，Lax-Wendroff格式会产生明显的非物理振荡，这些振荡会导致间断区域的误差是一阶的，最终把整个计算的全局收敛阶拉低到一阶。简单说，线性的二阶格式没法适配拟线性方程的间断解特性，精度会被间断区域“拖后腿”。

这里是你提供的Lax-Wendroff格式Mathematica代码，我补全了主循环里的计算逻辑，方便直接运行：

(*初始数据*)
u0[x_] := 1 - UnitStep[x - 0.1];
Flux[u_] := 0.5 u^2;
u = u0[xTbl];
dt = Abs[\[Sigma]] / Abs[Max[u]] dx;
un = 0*u;
F = Flux[u];
(*Lax-Wendroff格式*)
t = 0;
While[t < tFin,
    (*主循环：Lax-Wendroff核心计算*)
    un[[2 ;; nx - 1]] = u[[2 ;; nx - 1]] - dt/(2 dx) (F[[3 ;; nx]] - F[[1 ;; nx - 2]]) + 
        (dt^2)/(2 dx^2) ( (u[[2 ;; nx - 1]] + u[[3 ;; nx]]) (F[[3 ;; nx]] - F[[2 ;; nx - 1]]) - 
        (u[[1 ;; nx - 2]] + u[[2 ;; nx - 1]]) (F[[2 ;; nx - 1]] - F[[1 ;; nx - 2]]) );
    (*边界条件：固定边界*)
    un[[1]] = u[[1]];
    un[[nx]] = u[[nx]];
    (*更新当前时刻的解和通量*)
    u = un;
    F = Flux[u];
    t += dt;
]

问题2：如何在拟线性问题中实现二阶收敛阶？ #

要在拟线性守恒律（比如带间断的无粘Burgers方程）中实现二阶收敛，核心是在保持二阶精度的同时，抑制间断附近的振荡，这就得用二阶TVD（总变差不增）格式，思路是“光滑区二阶，间断区单调”，具体做法包括：

• 二阶重构（Reconstruction）：对每个网格单元的左右界面进行二阶精度的状态重构，比如MUSCL重构。通过重构得到界面处的左右状态值，而非直接用网格平均值，这样能在光滑区域保留二阶精度。
• 搭配限制器（Limiter）：限制器是关键，它能在光滑区域允许二阶重构的斜率，在间断区域把斜率限制为一阶，避免振荡。常用的限制器有Minmod、Van Leer、Superbee等，不同限制器在精度和单调性上有不同权衡。
• 守恒型框架：必须采用守恒型格式，确保格式满足守恒律，这样才能正确捕捉激波等间断结构，这是拟线性守恒律求解的基础。

举个简单例子，基于MUSCL+Minmod限制器的二阶格式，会先对每个网格的左右界面进行重构，再用Riemann求解器计算界面通量，最后更新网格平均值，这样就能在光滑区保持二阶收敛，同时在间断区避免振荡，达到全局二阶精度。

内容的提问来源于stack exchange，提问作者Oleg Kravchenko