第 4–6 章 · 有限体积法与 Godunov 方法

理论就位，开始造"正确"的方法。有限体积格式把守恒律的积分形式逐格离散：网格单元内平均值、单元边界上通量——激波捕捉的舞台正式搭好。

有限体积格式的骨架

把 $[x_{j-1/2}, x_{j+1/2}]$ 记为单元 $j$，其内平均 $Q_j = \frac{1}{\Delta x}\int q\,dx$。积分守恒律离散为：

$F$ 是数值通量。相邻单元共享同一 $F_{j+1/2}$，所以总量的离散守恒自动成立（内部通量两两相消）。设计格式 = 设计数值通量，全书的核心命题在这一行。

Godunov 方法：在界面解黎曼问题

最"物理"的通量：界面两侧状态 $Q_j, Q_{j+1}$ 构成局部黎曼问题，取其精确解在界面处的值作通量。优点：忠实于物理、天然熵条件、激波干净。缺点：精确黎曼解（欧拉方程）迭代求解昂贵。Godunov 定理再补一刀：一阶 Godunov 格式也逃不出"≤ 一阶精度"的单调性宿命——高阶化的路只能另辟。

高阶化第一招：分片线性重构

把单元内的常数分布升级为线性分布（斜率由相邻单元平均估计，限制器保证不产生新极值），界面两侧状态变成 $Q^- , Q^+$，再交给黎曼求解器——这就是 MUSCL 框架。时间上再配二阶（Hancock：半步预测），得到正式的二阶激波捕捉格式。斜率限制器（minmod、superbee、MC）的选择是精度与鲁棒性的旋钮。

老师的提醒

• 手写一个一维线性系统的 Godunov 格式（对流方程即可）是理解本章的最快路径，一百行以内。
• 限制器调参的体感：算例振荡 → 换 minmod 或降 CFL；激波过糊 → 换 MC/superbee。

自检

• 数值通量的设计如何决定格式的性质？
• Godunov 方法的优点与两条代价？
• MUSCL 的"重构-演化-投影"各对应什么操作？

更深一层：有限体积 = 弱形式的"守恒律特供版"

有限体积与有限元的关系值得说透：两者都是弱形式离散，差异在检验函数——有限元用光滑形函数（伽辽金），有限体积用"单元指示函数"（每格为 1 其余为 0）。这个选择让 FVM 天然输出单元平均的守恒律——与激波弱解的语言无缝对接（弱解本来就是积分守恒的意义下定义的）。这也解释了为什么守恒律领域 FVM 一统天下而 FEM 需要专门设计（ discontinuity-capturing / 弱伽辽金）：不是方法优劣，是语言匹配。同样，FVM 的边界条件必须走通量语言（边界通量 = 物理边界条件），直接"指值"会造成过约束与反射。

FVM 的几何自由度还是"任意多边形网格"——只要会算面积与形心、会算法向，任何网格拓扑都一视同仁。这是它在工业 CFD（复杂几何）压倒结构网格差分法的根本原因：几何灵活性与守恒性可以兼得。

推导细读：Godunov 方法的"以物理换精度"哲学

Godunov 格式的每一步是"把界面当作局部黎曼问题的初值，解出微小时间内的演化，取其界面通量"。它在一阶格式里精度垫底，却有两个高贵的品质：熵解自动满足（黎曼解是熵解）与单一间断精确捕捉（Roe 性质的前身）。Godunov 定理宣判它一阶封顶后，整条高分辨率路线（MUSCL 重构 + 精确/近似黎曼 + 限制器）本质上是"保持 Godunov 的物理内核、提高它的多项式表达力"。精确黎曼求解器（欧拉方程的非线性迭代）太贵 → 近似求解器（Roe 线性化、HLL 两波、HLLC 三波）接棒——"物理忠实度 vs 计算成本"的让步谈判构成了这一整章的算法谱系。

在真实世界里：从教科书格式到工业求解器的距离

手写 Godunov/MUSCL 求解器（几百行）与工业求解器（百万行）之间的距离不在核心格式而在配套工程：自适应网格细化（AMR，尖峰区加密、守恒的粗细界面通量修正）、并行分区（界面通量的跨进程一致性——守恒在并行下靠"界面通量单点计算"维持）、物性/真实气体效应、源项刚性处理（化学、重力——半隐式分步）、以及贯穿始终的验证套件。开源界的 OpenFOAM、AMReX、clawpack 把这些工程件全部开源化了——学完本书后最划算的下一步是读 AMReX 的双曲求解器文档与算例，看教科书算法如何被工程化到 $10^5$ 核心并行。

专题深入：从 Godunov 到 DG——"物理换精度"的谱系

从 Godunov 到 DG 的谱系值得画出：Godunov（界面解黎曼问题——一阶、物理忠实、熵保证）；MUSCL（重构 + 限制器——二阶、物理+精度）；WENO（多模板光滑度加权——高阶、放弃严格 TVD）；DG（单元内高阶多项式——任意阶、hp 自适应）。每一代都在"物理忠实度 vs 多项式表达力"的天平上加码精度——谱系的方向是从"低阶物理"到"高阶几何"。DG 的单元间断（界面通量用黎曼求解器）正是 Godunov 思想的高阶化——谱系的家谱第一页都写着 Godunov。

数学细节展开：边界条件的"ghost cell"工艺

双曲边界条件的 ghost cell 工艺值得展开：边界面上先解（或近似解）一个法向黎曼问题（内状态 vs 外 ghost 状态），通量取黎曼通量——ghost cell 的值不是"抄边界条件公式"而是"把特征相容关系解出来"。亚声速出流的完整实现：一条特征入流（指定压力）+ 两条出流（外推速度与密度）——特征变换后各方向的"指定/外推"分配严格按特征符号。检验边界的正确性有一招灵：圆形声波从域中心发出，看它穿过边界时是否产生反射波纹——有反射就是边界处理与特征结构失配。

工程案例细节：出师项目的"验收清单"

出师项目（二维浅水溃坝求解器：MUSCL-Hancock + HLLC + MC 限制器 + 干底处理）的验收清单值得做成可执行规格书：精确解对照（理想溃坝的相似解或高分辨率参考解，激波位置误差 < 半格）；静水还原（湖面初值算 1000 步，扰动 < 机器精度 × 安全系数——well-balanced 测试）；守恒审计（总量漂移 < 1e-12）；网格收敛（三套网格激波位置收敛）；干底鲁棒（下游干床溃坝不产生负水深）。这五个验收点每一个都对应本书的一章——把它们全跑通的那一刻，这本书就完成了从"读过"到"会用"的转化。