一维黎曼问题数值解与计算程序.pdfVIP

一维 Riemann 问题数值解与计算程序 一维 Riemann问题，即激波管问题，是一个典型的一维可压缩无黏气体动力 学问题，并有 解析解。对它采用二阶精度 MacCormack两步差分格式进行数值求 解。同时，为了初学者入门和练习方便 , 这里给出了用 C 语言和 Fortran77编写的 计算一维 Riemann问题的计算程序，供大家学习参考。 MacCormack Riemann A-1 利用 两步差分格式求解一维 问题 1. 一维 Riemann问题 一维 Riemann问题实际上就是激波管问题。激波管是一根两端封闭、内部充 满气体的直管， 如图 A.1 所示。在直管中由一薄膜将激波管隔开， 在薄膜两侧充 有均匀理想气体（可以是同一种气体， 也可以是不同种气体） ，薄膜两侧气体 的压力、密度不同。当 t 0 时，气体 处于静止状态。当 t 0 时，薄膜瞬时 突然破裂，气体从高压端冲向低压端， 同时在管内形成激波、 稀疏波和接触间 断等复杂波系。 2. 基本方程组、初始条件和边界条件 图 A.1 激波管问题示意图 设气体是理想气体。一维 Riemann问题在数学上可以用一维可压缩无黏气体 Euler方程组来描述。在直角坐标系下量纲为一的一维 Euler方程组为： u f 0, 1 x 1 (A.1) t x u 其中 u u , f u 2 p (A.2) E ( E p)u 这里 、 u 、 p 、 E 分别是流体的密度、速度、压力和单位体积总能。理想气体 状态方程： 1 2 2 p 1 e 1 E u v （A.3 ） 2 初始条件： 1,u 0, p 1 ； 0.125, u 0, p 0.1 。 1 1 1 2 2 2 供参考 边界条件： x 1 和 x 1处为自由输出条件， u u , u u 。 0 1 N N 1 3. 二阶精度 MacCormack差分格式 MacCormack两步差分格式： 1 n u 2 un r f n f n j j