第 3–4 章 · 弱解、激波与熵条件
激波出现了，"逐点满足方程"已经说不通。这一节建立间断解的数学框架：弱解、Rankine-Hugoniot 跳跃条件、熵条件。这是理解一切激波捕捉格式的理论底座。
弱解：积分意义的解
乘光滑试函数 $\phi$ 积分、分部积分把导数甩给 $\phi$：
满足此式的 $q$ 就是弱解——不再要求可微，激波合法化。代价：弱解不唯一。同一个黎曼问题，"物理的激波"和"非物理的膨胀激波"都是弱解。
Rankine-Hugoniot：激波的速度公式
跨过以速度 $s$ 移动的间断，守恒律要求：
激波速度 = 通量跳跃 / 状态跳跃。对伯格斯方程可得 $s = (q_L + q_R)/2$。这条公式是激波数值验证的裁判：格式算出的激波位置必须满足它。
熵条件：挑出唯一的物理解
物理激波两侧特征线进入激波（信息被吸收），膨胀激波特征线离开（凭空造信息）——后者违反热力学第二定律的数学影子。常用判据：Lax 熵条件 $f'(q_L) > s > f'(q_R)$；或熵不等式（如人工粘性消失极限）。数值方法的隐性义务：离散格式必须只收敛到熵解——一阶单调格式自动满足；高阶格式可能收敛到非物理弱解（如欠压缩问题），这是高级陷阱。
老师的提醒
- "格式对不对"的最低门槛：算出的激波速度与 Rankine-Hugoniot 公式一致、且收敛到熵解。
- 接触间断不会被任何熵机制"修陡"——它被抹平是正常现象，比较解时别把它当 bug。
自检
- 弱解为什么需要额外条件挑出物理解？
- 写出 Rankine-Hugoniot 条件并说出它的验证用途。
- 接触间断在数值上有什么特殊困难？
更深一层：弱解不唯一是"物理的失语"，熵条件是"物理的补票"
第 3 节的核心思想值得再往深处推一步：守恒律的积分形式本来就比微分形式"更基本"（激波处微分失效但积分守恒仍成立），所以弱解不是数学的妥协而是物理的还原。真正的问题是积分守恒律丢信息了：同一个守恒律可以对应多种微观机制（粘性不同、热传导不同），它们的极限解不同——熵条件就是把你想要的那个物理机制"补票"回来。数学等价物：给守恒律加一个熵对 $(\eta, F)$，要求 $\eta_t + F_x \le 0$——不同熵对挑出不同的物理极限。数值格式的"熵稳定性"证明 = 构造离散熵不等式——Tadmor 学派的熵稳定格式（熵守恒通量 + 耗散项）正是这条路线的工程结晶。
工程警示：多维系统的熵条件仍未完全解决（标量有完整理论，系统只有充分条件），跨声速、多相流的数值结果有时差异巨大——根源常在熵机制处理不同。看到文献里同类问题不同结论，先查各自的熵处理。
推导细读：黎曼问题是"最小完备实验室"
黎曼问题（左右常数初值）是双曲系统的氢原子：一维等熵气体、欧拉、浅水方程的黎曼解都有半解析形式，是格式开发与验证的最小完备考场。它的三大用途：验证（数值解与精确黎曼解叠加对比，位置/速度/波形三对）、构造（Godunov 型格式的通量就是界面黎曼问题的解——Roe/HLL/HLLC 都是对精确解的近似策略）、标定（限制器、熵修复的参数在黎曼问题上调定）。更进一步，任意初值可以局部化为黎曼问题的拼接——Godunov 方法的全局哲学：把解看成相邻状态间无数个微型黎曼问题的演化。理解了这一点，整个激波捕捉学派的世界观就通了。
在真实世界里：接触间断的"数值体质"问题
不同波族对数值耗散的敏感度天差地别：激波（非线性场）有熵机制帮忙陡化，接触间断（线性退化的场——密度界面、涡量面）没有任何物理机制阻止数值抹平，纯被格式耗散单方面侵蚀。这就是为什么 Sod 管里激波干净、接触间断糊——正常现象。工程对策：低耗散迎风格式（Roe 好于 HLL）、界面追踪（Level-Set/VOF 与守恒律耦合）、或多矩有限体积。多相流、可 mix 混合、天体界面不稳定性的模拟质量，很大程度上就是"对接触间断多温柔"的问题。
专题深入：熵条件的"三语言"与 Tadmor 学派的熵稳定格式
熵条件的三种表述（Lax 特征线、熵不等式、粘性消失极限）的数值翻译值得展开：单调格式自动满足（熵条件的离散版）；高阶格式需熵粘性或熵稳定通量（Tadmor 学派——构造离散熵不等式）；Roe 格式需熵修复（Harten 修复——把临界特征值的 $|\lambda|$ 换成光滑化版本）。熵稳定格式是现代的"收敛性保证"：对一般非线性守恒律系统（含多维），熵稳定性是目前最强的数学保证——它把"格式收敛到熵解"从"经验"升级为"可证"。
数学细节展开：接触间断的"数值体质"
接触间断（线性退化的场——密度界面、涡量面）的数值体质值得单列：没有任何物理机制阻止数值抹平——激波有熵机制帮忙陡化，接触间断纯被格式耗散单方面侵蚀。这就是为什么 Sod 管里激波干净、接触间断糊——正常现象。工程对策：低耗散迎风格式（Roe 好于 HLL）、界面追踪（Level-Set/VOF 与守恒律耦合）、多矩有限体积。多相流、可 mix 混合、天体界面不稳定性的模拟质量，很大程度上就是"对接触间断多温柔"的问题。
工程案例细节：Riemann 求解器的"工程体检"
黎曼求解器的工程体检（验证一个新求解器或新格式）：Sod 激波管（三波系——位置、速度、波形三对）、Lax 问题（更强的激波）、慢流（Mach 数趋于零时的精度——预条件需求）、强激波（Mach 数高时的鲁棒性）。"一套黎曼求解器的可靠性，就是它在四张考卷上的成绩单"——这四张考卷是双曲数值方法的工业验收标准。
例题
线性化输运问题：$h_t + (a h)_x = 0$（$a = 1\ \mathrm{m/s}$），上下游 $h_L = 4\ \mathrm{m}$、$h_R = 1\ \mathrm{m}$ 形成间断。用 Rankine-Hugoniot 条件求间断传播速度；再说明真实浅水方程（圣维南）里波速由什么决定。
解答
RH 条件：
该线性模型下间断以平流速度整体平移。
真实浅水方程：$h_t + (hu)_x = 0$、$(hu)_t + \big(hu^2 + \frac{g}{2}h^2\big)_x = 0$ 的特征值 $\lambda_\pm = u \pm \sqrt{gh}$——波速由重力波速 $\sqrt{gh}$ 决定：$h = 1\ \mathrm{m}$ 时 $\sqrt{9.8} \approx 3.1\ \mathrm{m/s}$，远快于平流速度；溃坝解的激波/稀疏波速度都由它派生。
- 要点：系统问题的波速 = 雅可比矩阵的特征值，不是对流速度；黎曼求解器必须逐特征场迎风。