第 1 章 · 引言：守恒律为什么特殊

双曲守恒律 $\partial_t q + \partial_x f(q) = 0$ 看着简单，却是数值方法里最挑剔的一类问题。这一节先弄清"守恒"两个字的分量：它不是一种物理偏好，而是决定解存在性与唯一性的数学结构。

守恒律的积分形式才是本体

对区间 $[x_1, x_2]$ 积分：

文字版：区间内的总量变化 = 左端流入 - 右端流出。总量守恒是这套方程的身份证——质量、动量、能量守恒律都是它。数值方法若不能维持离散版的这条恒等式，激波位置就会错。

最小例子：对流方程与伯格斯方程

• 线性对流 $q_t + a q_x = 0$：形状以速度 $a$ 平移，特征线是直线。它是一切双曲问题的"氢原子"。
• 无粘伯格斯 $q_t + (q^2/2)_x = 0$：通量非线性后，波速 $f'(q) = q$ 依赖状态本身——高处跑得快，波形前缘越来越陡，有限时间形成激波。光滑初值自动演化出间断，这是与线性方程的本质区别。

特征线：双曲问题的"光线"

$qt + f(q)x = 0$ 沿曲线 $x'(t) = f'(q)$ 满足 $q$ 不变——特征线是信息传播的轨迹。特征线相交 = 同一点被要求有两个值 = 激波诞生；特征线张开 = 稀疏波。整个激波捕捉理论都建立在对特征线的追踪上。

老师的提醒

• 数值实验从对流方程开始：它有精确解（平移），任何格式的问题在对流方程上无所遁形。
• 记住因果链"非线性 → 陡化 → 间断"：它解释了为什么光滑初值也必须准备激波捕捉格式。

自检

• 守恒律的积分形式与数值方法的什么要求对应？
• 伯格斯方程的激波是怎么从光滑初值产生的？
• 特征线相交和张开分别对应什么波？

更深一层：双曲、抛物、椭圆——"时间之箭"的数学分野

守恒律属于双曲家族，这个分类决定了它数值处理的一切性格。三种类型的物理画像：双曲（波动方程）——信息以有限速度传播、有清晰"来路"（特征线）、解可以保有间断；抛物（热传导）——信息瞬时传遍全域（理论上）、解瞬间光滑化（$t > 0$ 后无限可微）；椭圆（拉普拉斯）——无时间之箭、全局耦合、任何局部扰动改写全场面貌。数值方法的性格必须匹配：双曲用显式推进 + 迎风（尊重有限传播速度）、抛物用大步长隐式（光滑化是朋友）、椭圆用全局迭代（多重网格是VIP待遇）。PDE 类型 = 数值方法的宪法——跨型抄方法（拿椭圆的中心差分去解双曲）必出乱子。

伯格斯方程的激波形成时间也有闭式答案值得记：$q_t + (q^2/2)_x = 0$，初值 $q_0(x)$，激波时间 $t_s = -1/\min(q_0'(x))$——初值导数最负处最先"撞车"。这给了一个免费的先验诊断：拿到初值先算 $\min q_0'$，就知道你的模拟什么时候需要激波捕捉能力上线。

推导细读：弱解的"测度论自由"与它的代价

弱解允许间断的代价是不唯一，这个不唯一性可以做一次直观实验：标量守恒律 $q_t + (q)_x = 0$（线性）只有唯一弱解（线性问题无病），但伯格斯方程的"左小右大"初值，膨胀激波（速度 $s = (q_L+q_R)/2$ 的间断）与扇形稀疏波都是合法弱解——积分恒等式分不出对错。熵条件充当裁判的物理理由：真实流体过激波熵增（第二定律），膨胀激波熵减被宇宙禁止。数值的转译：格式必须"知道"熵条件——单调格式自动满足；高阶格式靠熵粘性/熵稳定通量保证。第 5 节的 Roe 熵修复就是这个保证的工程补丁。

在真实世界里：双曲系统在哪里"营业"

双曲守恒律是几大工业的数学引擎：油气输运管道（一维欧拉/等温模型，水锤与管线断裂的模拟）、血管血流（一维 blood flow 模型，脉搏波传播 = 特征线）、交通流（LWR 模型，堵车向后传播的激波 = 交通激波）、浅水水利（溃坝、风暴潮）。共同点：都存在真实间断（水锤波前、交通堵车尾）且间断的位置与速度就是工程答案——这就是为什么这些行业对激波捕捉格式的研究投入如此认真。

专题深入：双曲守恒律的"交通流"原型

守恒律的最直观原型是交通流：密度 $\rho$（车流密度）守恒，通量 $f(\rho) = \rho v(\rho)$（流速随密度降）——堵车向后传播的"激波"是守恒律激波的日常化身。这个原型的教学价值：非线性通量 ⟹ 波速依赖状态 ⟹ 光滑初值自动演化出间断——刹车灯亮起、车流密度波向后传播，就是 Burgers 方程的现场直播。这个例子的工程延伸：交通流模型（LWR 模型）是智能交通的流量预测基础，激波捕捉格式在交通仿真里的应用是守恒律数值方法的直系出口。

数学细节展开：特征线的"信息图谱"

特征线是双曲问题的"光线"：$q_t + f'(q)q_x = 0$ 沿曲线 $x'(t) = f'(q)$ 满足 $q$ 不变——特征线是信息传播的轨迹。特征线相交 = 同一点被要求有两个值 = 激波诞生；特征线张开 = 稀疏波。整个激波捕捉理论都建立在对特征线的追踪上。这个"信息图谱"的思维方式（谁影响谁、沿哪条线、多快）是双曲问题的第一分析工具——边界条件的特征计数、CFL 条件的依赖域论证、Riemann 问题的波分解，全是它的应用。

工程案例细节：守恒律的"工业版图"

双曲守恒律的工业版图值得画出：油气输运管道（一维欧拉/等温模型，水锤与管线断裂的模拟）、血管血流（一维 blood flow 模型，脉搏波传播 = 特征线）、交通流（LWR 模型，堵车向后传播的激波）、浅水水利（溃坝、风暴潮）。共同点：都存在真实间断（水锤波前、交通堵车尾）且间断的位置与速度就是工程答案——这就是为什么这些行业对激波捕捉格式的研究投入如此认真。

常见误区与辨析（补充）

• "守恒律的间断是数值假象"：间断是物理的（激波、水锤、堵车），数值只是忠实再现它。
• "光滑初值不会有激波"：Burgers 方程的光滑初值有限时间形成激波——非线性通量的陡化是物理不是数值。