Skip to content

双曲守恒律的 Riemann 问题

Riemann problem for hyperbolic conservation laws · Conservation-law Riemann problem · 守恒律黎曼问题

求解单个初始跳跃产生的自相似激波、稀疏波与接触波,并作为守恒格式的局部波传播模型。

条目类型
模型

形式陈述

守恒律的 Riemann 问题取分段常初值

u(x,0)={uL,x<0,uR,x>0.

方程与初值在缩放 (x,t)(αx,αt) 下不变,所以唯一熵解若存在,应为自相似形式 u(x,t)=U(ξ),其中 ξ=x/t。在光滑的非常值区代入方程得到

(f(U)ξ)U(ξ)=0.

因此解由常状态、满足 ξ=f(U) 的稀疏扇,以及沿固定 ξ=s 的跳跃拼成。跳跃速度由Rankine–Hugoniot 条件决定,多个弱拼法再由熵条件选择。

f 严格凸,uL>uR 给单一激波,

U(ξ)={uL,ξ<s,uR,ξ>s,s=f(uR)f(uL)uRuL.

uL<uR,给稀疏波

U(ξ)={uL,ξf(uL),(f)1(ξ),f(uL)<ξ<f(uR),uR,ξf(uR).

严格凸性保证 f 可逆。非凸通量须改用通量图在端态区间上的适当凸包或凹包,解可能是多段复合波。

系统 Ut+F(U)x=0 的局部 Riemann 解通常按特征值顺序包含 m 个波族,并由中间常状态连接。严格双曲性、真正非线性或线性退化等假设决定局部结构;“两边各给一个向量”并不保证任意大数据都有简单、唯一的全局波扇。

直觉

Riemann 初值把复杂波形压缩成一个最小实验:界面左、右各是常状态,所有变化都从原点发出。由于没有固有长度尺度,任意时刻的图形只是在 x 方向按 t 放大,同一条射线 x/t=ξ 上看到同一状态。求解问题因此转化为在“速度轴”上排列波。

这也是有限体积法采用 Riemann 问题的原因。一个时间步开始时,分片常单元在每个界面形成局部跳跃;只要相邻波在本步内不过度相互作用,界面处 x/t=0 的状态就给出物理通量。Godunov 方法不是猜一个平均状态,而是在每个面询问这个最小连续问题。

例子与边界

对 Burgers 方程,取 uL=2,uR=1。因为左态较大,解是激波,速度

s=f(1)f(2)12=1/223=12.

t=0.4 时跳点在 x=st=0.2x=0.1 处取值 2x=0.3 处取值 1。注意左右特征速度分别为 21,满足 1<1/2<2,所以这是压缩激波。

改取 uL=1,uR=2,稀疏解为

U(ξ)={1,ξ1,ξ,1<ξ<2,2,ξ2.

(x,t)=(0.2,0.4)ξ=0.5,故 u=0.5;在同一时刻扇区横跨 x=0.4x=0.8。这些数值由特征速度直接确定,不能把扇区边界任意画得更宽来“平滑”间断。

非凸通量是重要边界:端态次序不再单独决定波型,必须比较弦与切线,复合波中某一段激波速度还需与相邻稀疏扇边缘匹配。系统中的接触波允许等号特征速度,强激波可能离开局部小数据理论。真空、相变与共振还会让中间状态或严格双曲性失效。

Riemann 解也不等于边界条件。计算域端点的“边界 Riemann 问题”一侧是内部迹,另一侧是外部或边界状态;可进入的波数由特征方向限制。把所有外部状态分量强塞进域内会产生不相容反射。

推论与应用

精确 Riemann 解给出 Godunov 数值通量;近似求解器如 Roe、HLL 或 HLLC 则保留选定的波速与中间状态信息,以减少成本。无论精确或近似,接口必须报告左右状态、各波速度、界面状态和最大传播速度,后者进入 CFL 限制。只返回一个通量数字会隐藏失败原因。

Riemann 问题还充当解析基准:激波位置检验守恒速度,稀疏扇边缘检验特征速度,复合波检验非凸通量处理。网格解可能在有限分辨率下涂抹激波,但加密后应在适当范数中逼近同一熵解;出现持续过冲、错误波序或错误平台状态,分别指向耗散、熵修正或求解器中间态问题。

参考资料
  • Randall J. LeVeque, Finite Volume Methods for Hyperbolic Problems, Cambridge University Press, 2002, Chs. 13–15.
  • Eleuterio F. Toro, Riemann Solvers and Numerical Methods for Fluid Dynamics, 3rd ed., Springer, 2009, Chs. 4–6.
  • Constantine M. Dafermos, Hyperbolic Conservation Laws in Continuum Physics, 4th ed., Springer, 2016, Chs. 9–10.
关系图谱7 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系