第 10–15 章 · 推广与验证

一维标量问题练完，最后一段路走向真实：方程组、多维、边界、收敛验证。这一节是全书的"出师考核表"。

方程组与多维

• 方程组：特征分解代替标量斜率——每个特征场独立限制（波限制器）；欧拉方程的三波系各自处理。
• 多维：维分裂（算子分裂，简单但有分裂误差）或非分裂（角向通量 + 横向传播修正）。无结构网格上按面法向做一维黎曼问题，配合梯度重构（最小二乘/Green-Gauss）。
• 源项与平衡律：$q_t + f(q)_x = \psi(q)$——重力、化学反应源项会破坏"重构-演化"的平衡（静水湖问题：重力 + 压力梯度应精确抵消），well-balanced 格式专治此类。

边界条件的正确姿势

双曲问题的边界条件由特征方向决定：入流特征数 = 需指定的量数。激波捕捉格式的边界实现（ghost cell 填充）要与特征理论一致——照抄中心差分的边界写法会产生反射波污染全场。

收敛验证的正确流程

1. 精确解考卷：对流、Sod、Lax 问题——位置、速度、波形三对； 2. 网格收敛：加密检查光滑区精度阶是否达到理论值（二阶格式应观察到 $\Delta x^2$ 下降）； 3. 守恒检查：总量随时间的变化 = 边界通量积分（误差应机器精度级）； 4. 物理合理性：正密度正压力、熵不减少。

老师的提醒

• 出师项目建议：手写浅水方程溃坝求解器（MUSCL-Hancock + HLLC + MC 限制器），对照解析解与文献图形——做完这本的内化就完成了。
• 守恒检查脚本与求解器同时写：事后补的守恒检查从不被认真看。

自检

• well-balanced 格式解决什么"平衡被数值破坏"的问题？
• 双曲边界条件的特征数规则？
• 四步验证流程缺一不可的原因？

更深一层：从一维到多维的"推广税"

一维标量练熟的武器推广到方程组与多维时，每一站都要交"推广税"。方程组税：标量的斜率限制升级为逐特征场的波限制（特征分解 × 限制器）；黎曼求解器逐面解"法向一维问题"——但真实的多维流动沿切向也有结构，横向耗散不足是分维方法的通病（转角异常、Carbonou 问题）。多维税：维分裂简单但有分裂误差（交互激波算例会露馅）；非分裂格式（Corner transport + 横向传播）更准但模板更宽；无结构网格上还要补梯度重构（Green-Gauss/最小二乘）与限制器的方向化。源项税：重力、化学源项破坏离散平衡——静水大气、湖底压力梯度这类"应精确平衡"的解会被普通格式破坏，well-balanced 格式把平衡解硬编码进重构。每一条税则都对应一类真实模拟事故——推广不是抄写，是重新通过验证。

推导细读：边界条件的"特征语言"实现

双曲边界条件的正确实现在工程事故里排前三。特征语言的实现路径：边界面上先解（或近似解）一个法向黎曼问题（内状态 vs 外 ghost 状态），通量取黎曼通量——外状态的构造按特征计数：入流特征决定必须指定的物理量、出流特征外推。以亚声速出流为例：一条特征入流（$u - a < 0$ 方向看向域外的那条——取决于法向朝向）⟹ 指定压力、其余外推。 ghost cell 的值不是"抄边界条件公式"而是"把特征相容关系解出来"。检验边界的正确性有一招灵：圆形声波从域中心发出，看它穿过边界时是否产生反射波纹——有反射就是边界处理与特征结构失配。

在真实世界里：出师项目的完整规格书

第 7 节的"出师考核"值得给一份可执行的规格书。项目：二维浅水方程溃坝模拟（MUSCL-Hancock + HLLC + MC 限制器 + 干底处理）。验收清单：精确解对照（理想溃坝的相似解或高分辨率参考解，激波位置误差 < 半格）；静水还原（湖面初值算 1000 步，扰动 < 机器精度×安全系数——well-balanced 测试）；守恒审计（总量漂移 < 1e-12）；网格收敛（三套网格激波位置收敛）；干底鲁棒（下游干床溃坝不产生负水深）。这五个验收点每一个都对应本书的一章——把它们全跑通的那一刻，这本书就完成了从"读过"到"会用"的转化。后续进阶路线：Toro 的《Riemann Solvers and Numerical Methods for Fluid Dynamics》（黎曼求解器圣经）与 AMReX/clawpack 的开源代码精读。

专题深入：浅水方程的"well-balanced"革命

浅水方程的 well-balanced 格式值得单独展开：静水湖（$h$ 随地形变化、$u = 0$）是物理上最平凡的解——但普通格式把它算出虚假流动（源项 $-gh\partial b/\partial x$ 与压力梯度 $-g h\partial h/\partial x$ 的离散抵消不精确）。well-balanced 格式把"静水还原"硬编码进重构（对水面高度 $\eta = h + b$ 而非水深 $h$ 插值、源项用与通量同一的离散）——数值格式对"已知平衡态"的精确保持是 well-balanced 的定义。这个思想推广到一切"平衡律"系统：静水大气、化学平衡、 gravitational equilibrium——"平衡态被数值噪声破坏"是带源项守恒律的通病，well-balanced 是药方。

数学细节展开：自适应网格细化（AMR）的"守恒接口"

自适应网格细化（AMR）的工程核心是粗细界面的守恒通量修正：细网格块的边界通量必须与粗网格面的通量守恒（refluxing——通量差修正）——否则总量漂移。AMR 的三件套：聚类标记（误差大的细网格块聚类成新块）、细网格推进（时间子循环——粗一步 = 细多步）、界面同步（通量修正 + 边界插值）。AMReX/Berger-Colella 的 AMR 是双曲求解器的工业标准——它在保持守恒的同时把算力花在激波与界面。

工程案例细节：从教科书到 HPC 的"最后一公里"清单

从本书到 HPC 工业求解器的最后一公里清单：AMR（自适应网格细化——尖峰区加密、守恒的粗细界面通量修正）；并行分区（界面通量的跨进程一致性——守恒在并行下靠"界面通量单点计算"维持）；物性/真实气体效应（多组分输运、状态方程的扩展）；源项刚性处理（化学、重力——半隐式分步）；验证套件（本书的考卷 + 文献基准）。开源界的 AMReX、clawpack 把这些工程件全部开源化了——学完本书后最划算的下一步是读 AMReX 的双曲求解器文档与算例，看教科书算法如何被工程化到 $10^5$ 核心并行。