MUSCL 重构 | Codeboy's Blog

形式陈述

MUSCL 从单元平均 $U_i$ 构造单元内分片线性函数

p_i(x)=U_i+\sigma_i\frac{x-x_i}{\Delta x},

其左右界面迹为

U_{i-1/2}^{+}=U_i-\frac12\sigma_i, \qquad U_{i+1/2}^{-}=U_i+\frac12\sigma_i.

因线性项关于单元中心积分为零， $p_i$ 的单元平均仍是 $U_i$ ；这保证重构本身不改变守恒量。随后在界面调用数值通量

\widehat f_{i+1/2} =\widehat f\!\left(U_{i+1/2}^{-},U_{i+1/2}^{+}\right)

并作有限体积更新。MUSCL 重构本身是有限体积法调用的界面状态构造工具；只有再配上数值通量与时间推进，才形成完整的高分辨率有限体积格式，不能把 reconstruction 单独分类为整个有限体积 method 的特例。

最简单的 minmod 斜率为

\sigma_i=\operatorname{minmod}(\Delta_i^-,\Delta_i^+), \qquad \Delta_i^-=U_i-U_{i-1},\quad \Delta_i^+=U_{i+1}-U_i,

其中

\operatorname{minmod}(a,b)= \begin{cases} \operatorname{sgn}(a)\min(|a|,|b|),&ab>0,\\ 0,&ab\le0. \end{cases}

同号差分表示局部单调，取较小斜率防止重构越过邻值；异号表示极值或间断，斜率置零。更一般的 MC、van Leer、superbee 限制器在 Sweby 图中作不同压缩—耗散折中。

重构只给空间半离散通量。要得到完整TVD 格式，还需合适 CFL 和 TVD/SSP 时间推进；把二阶空间重构直接配任意高阶时间方法，不能仅凭限制器名称保证 TVD。

直觉

一阶有限体积把每个单元画成平台，界面跳跃包含了光滑斜坡本可提供的信息。MUSCL 在每个平台上恢复一条保持平均值的线段，让左右界面看到更接近真实解的迹。若三点趋势一致，线段可以倾斜；若趋势在中心反向，限制器宁可退回水平平台，也不让线段跨出邻域范围。

限制器像局部可信度判断，而不是统一滤波器。平滑单调区的左右差分接近，保留中心斜率给二阶精度；激波与极值处差分冲突，斜率归零并回到一阶鲁棒通量。精度与非振荡在不同位置分别取舍，正是非线性高分辨率格式绕开线性单调阶障碍的方式。

例子与边界

取周期单元平均

U=(0,1,3,2).

逐单元左右差分分别为 $(-2,1)$ 、 $(1,2)$ 、 $(2,-1)$ 、 $(-1,-2)$ ，故 minmod 斜率

\sigma=(0,1,0,-1).

对 $a=1$ 的右行平流，迎风通量只取每个单元的右迹。四个周期界面从“单元 4 到 1”开始，通量依次为

(\widehat f_{1/2},\widehat f_{3/2}, \widehat f_{5/2},\widehat f_{7/2}) =(1.5,0,1.5,3).

取 $\lambda=0.4$ ，更新得到

U^{n+1}=(0.6,0.4,2.4,2.6).

总和由 $6$ 保持为 $6$ 。周期总变差从

|1|+|2|+|-1|+|-2|=6

降为 $|{-0.2}|+|2|+|0.2|+|{-2}|=4.4$ 。该算例也显示局部极值处 $\sigma=0$ ，而单调区单元 2 保留斜率 $1$ 。

minmod 较耗散，superbee 较压缩；后者在某些问题中保持更锐利跳跃，也更容易暴露网格对齐与系统分量限制问题。对系统逐个保守分量应用 limiter 可能制造非物理压力或密度，常需转到特征变量并额外施加正性限制。

非均匀网格上， $\Delta_i^\pm$ 必须按中心距离缩放后比较，界面迹也需使用真实几何距离；直接复制均匀网格的半斜率公式会失去线性精确性。边界单元缺少一侧邻居时，应使用与流入数据一致的单侧重构或降低阶数，不能越界读取。

推论与应用

MUSCL 的验证可拆成三项独立检查：重构后每个单元平均不变，受限迹位于允许邻域或正性集合内，完整时间步在目标 CFL 下满足 TVD/熵诊断。上面的四单元算例覆盖周期斜率、界面左右顺序、守恒更新和总变差，适合作为最小回归检查。

在光滑解上，受限斜率若还原中心导数，空间误差可达二阶；在激波附近，局部退化到一阶使震荡受控。若应用需要更高阶，可采用 PPM 或 WENO，但仍必须保留共享通量、边界重构与稳定时间推进这些接口。仅增加更宽 stencil 而不处理非光滑区，不会自动得到可靠高阶守恒格式。

参考资料

Bram van Leer, “Towards the Ultimate Conservative Difference Scheme. V. A Second-Order Sequel to Godunov’s Method,” Journal of Computational Physics 32(1), 1979, pp. 101–136.
P. K. Sweby, “High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws,” SIAM Journal on Numerical Analysis 21(5), 1984, pp. 995–1011.
Randall J. LeVeque, Finite Volume Methods for Hyperbolic Problems, Cambridge University Press, 2002, Chs. 6 and 16.