第二讲双曲型方程组及间断解课件.ppt

Riemann 问题五种可能情况 x t 激波 接触间断 激波 膨胀波 接触间断 激波 激波 接触间断 膨胀波 膨胀波 接触间断 膨胀波 膨胀波 膨胀波 （1） （2） （3） （4） （5） 如何区分这5种情况？ Copyright by Li Xinliang * 假设 准则如下： 情况1 情况3 情况4 利用函数 （由 *式定义） 函数性质很好—— 单调连续 情况5 情况5 情况4 情况3 情况1 Riemann 求解总步骤： 1） 根据上述判决区分情况 2） 按照上一页的步骤求解 真空区 Copyright by Li Xinliang * Riemann问题的具体计算步骤 1. 判断可能会出现的情况（五种情形之一） a. 定义函数 b. 进行判断 情况5 情况4 情况3 情况1 情况5 情况4 情况2 情况1 单调增函数，性质很好 计算出 ， 根据 的大小进行判断，具体见下图： Copyright by Li Xinliang * 2. 求解中心区的压力和速度 单未知数的代数方程，迭代求解（例如Newton法，F(p)性质好，求解不困难） 3. 确定中心区接触间断两侧的密度 以及左、右波传播的速度 a. 左波为激波的情况（情况1,3） b. 左波为稀疏波的情况 （情况2，4，5） 中心区 接触间断左侧的物理量 膨胀波的波头及波尾速度 激波的传播速度 对于情况（5），波尾速度为： 中心区为真空，音速 无定义，改由该式计算 Copyright by Li Xinliang * c. 右波为激波的情况（情况1，2） 中心区 接触间断右侧的物理量 b. 右波为稀疏波的情况 （情况2，4，5） 4. 计算稀疏波区域的值（如果有稀疏波的话） a. 左稀疏波 b. 右稀疏波 情况2，4 情况5： Copyright by Li Xinliang * 思考题： 上述求解方法要求间断两侧流场分布为常数，如果初始时刻流场分布是x的函数，怎样利用该理论解计算 ？ 提示： 把曲线离散化，看成折线 Copyright by Li Xinliang * 4 近似Riemann解初步 精确Riemann解计算量较大； 近似Riemann解： 积分型 （HLL,HLLC）、微分型（Roe） 4.1 HLL 近似Riemann解 （Harten, Lax & van Leer） Ref.： E. F. Toro: Riemann Solvers and Numerical Methods for Fluid Dynamics, Springer, 2009 (Third Edition) 基本原理： 双激波近似 t=0 t=t0 激波1，速度Z1 激波2，速度Z2 假设间断面产生两道激波，速度分别为Z1,Z2 根据质量、动量、能量守恒，容易计算出图中控制体积内的总质量、总动量、总能量 t0 时刻激波才传到控制体边界，因此0 到t0时刻，控制体边界处物理量保持0时刻的值。 利用总量，求出图中控制体内的平均值，作为该区域物理量的近似值 Copyright by Li Xinliang * 4.2 HLLC 近似Riemann解 （Toro） 发展了HLL近似解，用三波模型来近似 （如图） 三波近似， 左、右波的速 T 时刻的流动状态 激波 激波 接触间断 模型： 左右两道激波，中间有接触间断 激波速度已知为: ZL, ZR 未知数（4个）： 方程（6个）：两道激波的RH关系式 方程多了两个？ （因为假设激波速度已知） 常用方法： 去掉两个方程 去掉两个能量方程， 4个未知数，4个方程，求解 求解过程简单，轻易可给出表达式 满足： 物理意义为平均增长率 Copyright by Li Xinliang * 4.3 Roe 近似Riemann解 u f(u) uL uR uRoe Roe平均 Riemann问题： 近似： 用平均增长率替代瞬时增长率 常系数线性方程组，求解简单（相似变换解耦求解） Copyright by Li Xinliang * 作业 2.1 公式推导 （1） 一激波从左向右传播。激波左侧物理量为 ; 激波右侧压力为 ， 试