形式陈述
一维无源标量守恒律是在区域 ( x , t ) ∈ R × ( 0 , ∞ ) 上求
u t + ∂ x f ( u ) = 0 , u ( x , 0 ) = u 0 ( x ) , 其中 u 是一个守恒量,f : R → R 是通量。它是偏微分方程 公理库 偏微分方程 Partial differential equation · PDE 含未知多元函数及其偏导数的方程。 的一类,但“守恒”不只是把某个 PDE 改写成含导数的样子:对任意固定区间 [ a , b ] ,足够光滑的解都满足
d d t ∫ a b u ( x , t ) d x = f ( u ( a , t ) ) − f ( u ( b , t ) ) . 区间内总量的变化只由两端净流入决定。多维形式为 u t + ∇ ⋅ F ( u ) = 0 ;有源项时右端写成 s ( x , t , u ) ,此时体内生成与边界通量必须分别入账。
若 u 光滑,链式法则把方程写成
u t + f ′ ( u ) u x = 0. 沿满足 x ˙ = f ′ ( u ) 的特征,u 保持不变。非线性通量使速度依赖状态:后方特征若更快会相交,经典导数随之失效并形成间断;前方特征更快则会展开为连续扇形。因而即使 u 0 光滑,经典解通常也只能在有限时间内存在,之后必须转向弱解 公理库 守恒律的弱解 Weak solution of a conservation law · Distributional solution of conservation law · 守恒律分布解 通过测试函数积分恒等式定义允许跳跃的守恒律解,并保留初值、边界与跨间断通量平衡。 。
边界条件由特征方向决定。对区间上的平流型守恒律,只能在信息进入区域的边界给独立数据;流出边界的迹应由内部解决定。把左右两端都强行指定为 Dirichlet 数据,可能过度约束方程,也会破坏积分守恒。数值上则由有限体积法 公理库 有限体积法 Finite volume method · Finite-volume discretization · 控制体积法 对每个控制体积分守恒律,以共享界面通量更新单元平均,从离散结构上保证内部通量抵消。 直接维护每个控制体的通量账本。
直觉
可以把 u 看成沿管道分布的密度,把 f ( u ) 看成单位时间穿过截面的量。相邻小区间共享同一界面:一个区间的流出正是下一区间的流入,所以把局部平衡相加时内部通量成对抵消。这种抵消结构比点态可微性更基本;激波处虽然没有普通导数,跨越激波的总量仍能精确平衡。
特征线像携带状态标签的运输轨迹。线性通量 f ( u ) = a u 给所有标签相同速度 a ,波形只平移;凸非线性通量会让高状态与低状态速度不同。轨迹挤在一起时,守恒量不能凭空消失,只能集中成移动跳跃;轨迹拉开时,中间状态由一族不同速度连续填满。激波与稀疏波因此不是人为添加的现象,而是同一个局部守恒机制的两种几何结局。
例子与边界
取 Burgers 通量 f ( u ) = u 2 / 2 。Riemann 初值 u L = 2 , u R = 0 的左侧特征速度为 2 ,右侧为 0 ,特征相撞。守恒跳跃以速度
s = f ( 2 ) − f ( 0 ) 2 − 0 = 1 移动,所以 u ( x , t ) = 2 当 x < t ,u ( x , t ) = 0 当 x > t 。在 t = 0.6 时跳点位于 x = 0.6 ;这不是换一组数字的示意,而可由任意跨过该点的控制体直接核对净通量与总量变化。
若改为 u L = 0 , u R = 2 ,特征从速度 0 展开到 2 ,熵可接受解为
u ( x , t ) = { 0 , x / t ≤ 0 , x / t , 0 < x / t < 2 , 2 , x / t ≥ 2. 例如 ( x , t ) = ( 0.9 , 0.6 ) 给 x / t = 1.5 ,故 u = 1.5 。同一对端态也存在满足分布守恒的“膨胀激波”,说明方程和初值本身还不足以选择物理解;熵解 公理库 守恒律的熵解 Entropy solution of a conservation law · Kružkov entropy solution · 守恒律可容许解 在弱守恒之外加入单向熵耗散不等式,从多个分布解中选择满足唯一性与稳定性的可容许解。 补上方向性条件。
凸性不能被默认为永远成立。对非凸通量,u L < u R 不必对应单一稀疏扇,u L > u R 也不必对应单一激波;通量图的凸包或凹包可能产生激波与稀疏波的复合结构。多维中波面几何、边界迹和可能的奇异集中还会增加困难。本页只讨论标量守恒量,不能把结论直接搬到具有多个特征族的双曲方程组。
推论与应用
局部守恒给出两个后续接口。分析上,把方程乘测试函数并分部积分得到弱形式,移动间断则满足 Rankine–Hugoniot 跳跃关系;再加入熵不等式,才能在适当函数类中恢复唯一性与 L 1 收缩。数值上,离散更新若让公共界面只出现一个数值通量,所有内部贡献也会望远镜式抵消,从而在机器精度范围内守恒总量。
守恒形式与非守恒形式在光滑区可由链式法则互换,在间断处却可能给出不同传播速度。计算 u t + f ′ ( u ) u x = 0 时若直接离散乘积而不保持通量差,所得极限即使存在也未必满足正确跳跃平衡。模型验证至少应报告通量、初边值域、允许的解概念及边界流入规则;只给一条平滑数值曲线无法说明激波位置或总量是否正确。
参考资料
Lawrence C. Evans, Partial Differential Equations , 2nd ed., American Mathematical Society, 2010, §3.4.
Constantine M. Dafermos, Hyperbolic Conservation Laws in Continuum Physics , 4th ed., Springer, 2016, Chs. 1 and 6.
Randall J. LeVeque, Finite Volume Methods for Hyperbolic Problems , Cambridge University Press, 2002, Chs. 11–12.