“Godunov 格式在每个时间层把单元平均 $U i^n$ 解释为分片常函数。每个界面 $(U i^n,U {i+1}^n)$ 形成一个局部 Riemann 问题;设其熵解在 $\xi=0$…”
形式陈述 ​
守恒律的 Riemann 问题取分段常初值
方程与初值在缩放
因此解由常状态、满足
若
若
严格凸性保证
系统
直觉
Riemann 初值把复杂波形压缩成一个最小实验:界面左、右各是常状态,所有变化都从原点发出。由于没有固有长度尺度,任意时刻的图形只是在
这也是有限体积法采用 Riemann 问题的原因。一个时间步开始时,分片常单元在每个界面形成局部跳跃;只要相邻波在本步内不过度相互作用,界面处
例子与边界
对 Burgers 方程,取
故
改取
在
非凸通量是重要边界:端态次序不再单独决定波型,必须比较弦与切线,复合波中某一段激波速度还需与相邻稀疏扇边缘匹配。系统中的接触波允许等号特征速度,强激波可能离开局部小数据理论。真空、相变与共振还会让中间状态或严格双曲性失效。
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.