第 6–8 章 · MUSCL 与高分辨率格式

一阶格式稳定而糊，二阶格式锐而振。这一节讲"两全"的完整技术栈：高分辨率格式的构造原则——在光滑区高阶、在间断处自动降阶，判据就是限制器。

高分辨率格式的三项全能

一个合格的高分辨率格式： ① 单调保持（不产生新极值）；② 光滑区二阶以上；③ 间断处收敛到熵解（激波干净无振铃）。三项靠以下组件拼成：分片线性重构（MUSCL）+ 限制器 + 好的黎曼求解器 + 半步时间演化（Hancock 或 RK2）。

限制器的函数学

限制器 $\phi(r)$ 以相邻梯度比 $r = \theta_{j}/\theta_{j-1}$ 为输入，输出斜率缩放系数。全都在 TVD 区域内（Sweby 图）：

• minmod（$\phi=\max(0,\min(r,1))$）：最保守，最糊。
• MC limiter：精度与鲁棒的均衡，默认推荐。
• superbee：最锐，激波最陡但可能角化光滑波。
• TVD 的天花板：TVD 格式在极值点自动降为一阶（极值附近 $r\le 0$）——比 TVD 更弱的单调保持（如 WENO 的光滑指示器）能突破这个限制，在极值附近保高阶。

时间维度的精度

空间二阶 + 时间一阶 = 总体一阶。补齐时间：MUSCL-Hancock（预测-校正半步）或显式中点/RK2。CFL 数建议 0.4-0.6（二阶重构缩小了稳定域）。最终验收：Sod 问题 + 光滑对流问题（Isentropic vortex）双考卷——前者看激波锐度，后者看光滑区精度阶。

老师的提醒

• TVD 不是金标准而是底线：它保证不振荡，代价是极值点精度下降；WENO 是现代替代（极值附近五阶），但更贵更难调。
• 网格上看到"台阶状角化"的光滑波，是 superbee 类限制器过锐的签名——换 MC 先试。

自检

• 高分辨率格式的三项全能分别由哪个组件保证？
• 为什么 TVD 格式在光滑极值附近必降阶？
• 验收一个新格式的双考卷是什么？

更深一层：TVD 是"防线"而非"终点"

TVD（总变差不增）保证不振荡的证明干净利落，但它有一个鲜为人知的代价——精度在极值点系统性降为一阶（Sweby/Tadmor 的结果：二阶 TVD 格式在光滑极值处退化为 TVD 域边界上的低阶方案）。这不是实现瑕疵而是数学必然：极值点附近任何真正的二阶重构都会产生新极值。三条突破路线构成现代高阶方法的版图：WENO（用多个候选模板的光滑度加权，极值点附近自动切换到最光滑模板——保持高阶但放弃严格 TVD）；DG 有限元（单元内高阶多项式 + 限制器，几何与高阶兼得）；熵稳定格式（放弃 TVD 转而证明离散熵不等式——不振荡的另一种充分性）。TVD 是防线的底座，不是大厦的顶层——这个定位决定了你读到"XX 格式比 TVD 更优"时该有的正确期待。

限制器家族的性格对比也值得整理成卡片：minmod（最糊最稳，鲁棒性首选）、superbee（最锐最险，激波最陡但光滑波角化）、MC（均衡默认）、Van Leer/Van Albada（光滑性更好的中间派）。限制器选择 = 锐度与鲁棒性的连续旋钮，工业实践常提供 3-4 档让用户按算例脾气切换。

推导细读：WENO 的"加权民主制"

WENO（加权基本无振荡）的思路是"别押注单一模板，搞加权民主"：五阶 WENO 用三个三阶候选模板，各配一个光滑度指示器 $\beta_s$（模板上导数的离散积分），权重 $w_s \propto \frac{\alpha_s}{\sum\alpha}$、$\alpha_s = \frac{d_s}{(\beta_s + \epsilon)^2}$——模板越光滑权重越大，跨间断时含间断的模板权重自动趋零。光滑区权重回归理想值 ⟹ 五阶；间断区自动降级 ⟹ 无振荡。"用光滑度投票"是个可迁移的伟大思想：它以某种形式出现在自适应滤波、图像修复（全变分正则）、数据同化的混动方法里——哪里需要在"高阶"与"稳健"之间动态取舍，哪里就有它的变体。

在真实世界里：高分辨率格式的工业化调参

工业求解器里限制器/格式的调参有明确的决策路径：算例以激波为主（外流气动、爆炸）→ 锐度优先（superbee/Minmod-3 类）+ CFL 取 0.4；算例以光滑波与涡为主（声学、湍流 LES）→ 低耗散优先（WENO5/中心型 + 显式滤波）；不确定 → MC 限制器 + CFL 0.45 起步。调参的验证闭环永远不变：光滑波保真测试（相位与幅值误差）+ Sod 管锐度测试 + 长时间稳定性测试，三关全过才换下一个参数。格式的"脾气"是跑出来的，不是读出来的。

专题深入：TVD 区域的"Sweby 几何"与 WENO 的突破

Sweby 的 TVD 区域图（限制器 $\phi(r)$ 必须落在由 minmod 与 MC 围成的区域内）值得作为"限制器的几何宪法"记住：区域内的一切二阶 TVD 格式、区域外振荡。TVD 的天花板：极值点附近 $r \le 0$ 强制降阶——光滑极值处二阶 TVD 格式必降一阶（数学必然）。WENO 的突破路线：放弃严格 TVD，改用"光滑度加权"——三个候选模板按光滑度投票，间断处含间断的模板权重趋零。WENO5 在光滑区五阶、间断区无振荡——这是"用光滑度投票"对"TVD 区域"的超越。

限制器的工程选择口诀：算例以激波为主（外流气动、爆炸）→ 锐度优先（superbee/Minmod3 类）+ CFL 0.4；以光滑波与涡为主（声学、LES）→ 低耗散优先（WENO5/中心型 + 显式滤波）；不确定 → MC 限制器 + CFL 0.45 起步。

数学细节展开：WENO 的"光滑度指示器"设计

WENO 的光滑度指示器 $\beta_s$（模板上导数的离散积分）的设计值得一看：$\beta = \sum_{l=1}^{r}\int|\frac{d^lp}{dx^l}|^2 dx^{2l-2}$ 的离散版——模板越振荡 $\beta$ 越大。权重 $w_s \propto \frac{d_s}{(\beta_s + \epsilon)^2}$（$d_s$ 是理想权重、$\epsilon$ 防 除零）。"用光滑度投票"是一个可迁移的伟大思想：它以某种形式出现在自适应滤波、图像修复（全变分正则）、数据同化的混动方法里——哪里需要在"高阶"与"稳健"之间动态取舍，哪里就有它的变体。

工程案例细节：CFL 与限制器的"联动调参"

CFL 与限制器的联动调参值得展开：二阶重构缩小了稳定域（CFL 建议 0.4-0.6），限制器越锐稳定域越窄（superbee 比 minmod 的 CFL 上限低）。调参的验证闭环永远不变：光滑波保真测试（相位与幅值误差）+ Sod 管锐度测试 + 长时间稳定性测试，三关全过才换下一个参数。格式的"脾气"是跑出来的，不是读出来的——这句话对限制器选择尤其成立。