第 6 章 · 近似 Riemann 求解器

精确黎曼求解器太贵（欧拉方程要迭代解非线性方程组）。这一节的思想：通量不必精确，只需满足正确结构——近似黎曼求解器用线性化问题替代原问题，成本大降而质量几无损。

Roe 求解器：线性化的王者

在界面状态 $Q_L, Q_R$ 附近找一个"平均矩阵" $\hat A$，使 $\hat A(Q_R - Q_L) = f(Q_R) - f(Q_L)$（Roe 性质：任何跳跃都被线性问题精确复现）。数值通量就是线性黎曼解的通量：

三段结构：中心部分（平均通量）+ 按特征值绝对值的迎风耗散 + 特征分解系数。每一个近似的激波捕捉格式都是这个公式的变体——差别只在 $\hat\lambda$ 与特征分解怎么算。

波分解的统一视角

黎曼问题的解是一组波（特征场）的叠加：激波/稀疏波 + 接触间断。近似求解器就是把每个波简化为一个跳跃。由此立即得到两个著名变体：

• HLL：只保留最外两条波速，中间用一个平均状态——极稳健、极耗散（接触间断全糊）。
• HLLC：在 HLL 中间补回接触间断——稳健与保真兼顾，工程常用。

熵修复

Roe 格式有一处理论裂缝：跨声速稀疏波内它可能收敛到非物理膨胀激波（违反熵条件）。修补是熵修复：把临界特征值的 $|\lambda|$ 换成光滑化的 $|\lambda|_{smooth}$（Harten 修复），加一点点耗散保住熵解。这是"格式正确性"的最后一道保险。

老师的提醒

• 选求解器的体感顺序：Roe+熵修复（基准）→ HLLC（更稳）→ HLL（极端算例兜底）；HLLE 的过度耗散常被误诊为网格问题。
• 检查膨胀激波 bug 的标准算例：跨声速稀疏波（Sod 变体 Lax 问题的某些参数）——没有熵修复的 Roe 会在中间算出一个"假激波"。

自检

• Roe 性质是什么？它带来什么精度红利？
• HLL 为什么糊？HLLC 补了什么？
• 熵修复加在哪里、防什么病？

更深一层：Roe 平均的"魔法"与它的适用边界

Roe 求解器的核心是一个看似天真的要求：找一个常数矩阵 $\hat A(Q_L, Q_R)$，使得 $\hat A(Q_R - Q_L) = f(Q_R) - f(Q_L)$——任何跳跃都被线性问题精确复现。满足 Roe 性质的"平均状态"对欧拉方程有闭式解（Roe 平均：含 $\sqrt{\rho}$ 加权的速度、焓平均）。魔法有三重：单一激波被精确捕捉（位置、速度零误差）、常数状态零耗散、通量差分裂的完整特征结构。适用边界同样要背熟：跨声速稀疏的膨胀激波病（熵修复治）、壁面附近的强非物理状态（需裁剪）、以及多维"转角污染"（一维黎曼逐面求解的固有限制——横向耗散不足，跨分量扩散需要辅助）。

从更广的视角看，Roe/HLL/HLLC/AUSM 这些家族的差异可以总结成一句话：对黎曼解"简化到什么程度"的不同报价——Roe 全波精确但每面要做特征分解，HLL 只保留最外波速便宜而糊，HLLC 补回接触间断是性价比甜点。选型的决策输入只有两个：你的问题里接触间断/物质界面重要吗？你的算力预算允许每面多少次三角函数与开方？

推导细读：波分解的"守恒 + 迎风"合体

把 Roe 通量公式逐项读出味道：

第一项是中心通量（平均，无耗散）；第二项按特征场逐个加迎风耗散（特征值绝对值 × 跳跃幅度 × 特征向量）。这个结构是现代激波捕捉格式的"标准骨架"——所有格式的差异都在" $\hat\lambda$ 与 $\hat\alpha$ 怎么算"：Roe 用 Roe 平均、HLL 用极端波速夹逼、AUSM 用通量分裂的代数变体。识别出这个骨架，阅读任何新格式的论文都变成"它对哪一项做了什么改动"的填空题——文献阅读效率的质变。

在真实世界里：格式性能的"一张图"验证法

工程上快速鉴定黎曼求解器，有一张经典验证图：Sod 激波管的密度剖面叠加比较（Roe vs HLLC vs HLL vs 精确解）。看三个位置：激波（都该干净——它们都对激波友好）、接触间断（Roe 锐、HLLC 可、HLL 糊到十格宽）、稀疏波内（无熵修复的 Roe 可能现"膨胀激波"平台——假解的签名）。这张图 + 跨声速稀疏管一张图，十五分钟内可完成求解器的核心鉴定。开源验证参考：一维欧拉精确解的程序（Toro 书附带的 EXACT 算法）值得自己实现一份，它是所有双曲求解器的终身裁判。

专题深入：Roe 平均的"数学配方"与它的代数性质

Roe 平均的配方值得单独解剖：对欧拉方程，Roe 平均状态由 $\sqrt{\rho}$ 加权的速度、焓平均给出——$u_{roe} = \frac{\sqrt{\rho_L}u_L + \sqrt{\rho_R}u_R}{\sqrt{\rho_L} + \sqrt{\rho_R}}$ 型。这个"sqrt 加权"不是任意的：它恰好让 Roe 性质（$\hat A(Q_R - Q_L) = f(Q_R) - f(Q_L)$）成立。Roe 平均是"为满足一个代数条件而反推的平均"——它的每一项都是 Roe 性质的解。这个思路的推广：任何守恒系统都可以问"什么平均满足 Roe 性质"——有的系统有闭式（欧拉、浅水）、有的没有（磁流体 MHD 的 Roe 平均带约束问题）。

Roe 格式的代数性质还有一条值得记：常数状态零耗散（$Q_L = Q_R$ 时通量恰为物理通量——均匀流不被扰动）。这是"格式不产生虚假物理"的最低要求，但一些高阶格式在非均匀网格上会违反它（自由流保持性 freestream preservation——坐标变换度量项与通量散度的离散不一致）。

数学细节展开：波分解的"骨架"与格式谱系

Roe 通量公式的骨架 $F = \frac12(f_L + f_R) - \frac12\sum_p|\hat\lambda_p|\hat\alpha_p\hat r_p$ 值得作为"所有激波捕捉格式的标准骨架"记住：中心部分（平均通量）+ 按特征场逐个加的迎风耗散（特征值绝对值 × 跳跃幅度 × 特征向量）。每一个近似的激波捕捉格式都是这个公式的变体——差别只在 $\hat\lambda$ 与特征分解怎么算：Roe 用 Roe 平均、HLL 用极端波速夹逼、HLLC 补回接触间断、AUSM 用通量分裂的代数变体。识别出这个骨架，阅读任何新格式的论文都变成"它对哪一项做了什么改动"的填空题。

工程案例细节：MHD 与多物理场的"黎曼求解器扩展"

黎曼求解器向多物理场的扩展清单：磁流体力学（MHD）——八个波（快/慢磁声波、Alfvén 波、熵波），Roe 平均存在但需要 $\nabla\cdot\mathbf{B} = 0$ 的保持技术（约束输运 CT/求解投影）；两相流——界面处的相变与物性跳跃，五方程/七方程模型；化学反应流——多组分输运 + 刚性源项的分步处理。每一次物理扩展都要求黎曼求解器"长出"新的波族——MHD 的 Alfvén 波与化学的组分波都是双曲谱系的新成员。