高等计算流体力学-05研讨.pptxVIP

第五讲迎风型有限体积格式的通量计算方法-Godunov格式 1 十六、一维Euler方程的Godunov格式及Riemann问题 2 一维问题的有限体积格式格式 考虑Euler方程 我们采用有限体积方法离散上述方程。在任意控制体 上积分上式,有: 3 其中 则 到目前为止整个过程是精确的。 采用零阶重构： 如何有从间断的初值出发计算通量? 4 Godunov的方案! 延拓! Riemann问题 初值 方程 计算数值通量须在每个时间步求解多个Riemann问题! 精确解 5 Riemann问题 所谓Riemann问题就是求解Euler方程 ? 在初值 下的解。 6 预备知识 一维一阶非线性守恒系统(守恒律)可写为下列一般形式 x∈R,t 0 守恒律可写为拟线性形式 7 流体力学无粘流动的Euler方程是典型的非线性守恒律，可以写为 其中： 这里ρ，u，p，E，H分别为密度、速度、压力、总能和总焓。 为内能， 为焓。γ为比热比，对于空气，γ=1.4。 对于完全气体， 写成拟线性形式，其Jacobi矩阵为 8 特征值与特征向量 左特征向量为行向量 右特征向量为列向量 具体到Euler方程，左特征向量和特征值为： 9 特征线与Riemann不变量 定义曲线 沿该曲线 特征线 特征相容关系 10 Euler方程特征相容关系为 其中 为熵 均熵流动 全场 Riemann不变量 一般情况 11 Euler方程解的结构 在均熵流动条件下，有 且全场 Euler方程的光滑解有如下几种可能 1）在求解域中，Riemann不变量 这是最一般的情况，Euler方程的解比较复杂，通常无解析解。 均不为常数。 2）均匀流：Riemann不变量 均为常数。此时，令 此时流动是均匀的。 12 3) 简单波：有一个Riemann不变量在某区域内为常数 考虑 的情况, ( 的情况类似 ) 全场 沿特征线 沿特征线 由均熵条件，密度 和压力 沿特征线 也为常数。 流线 从左侧穿过特征线 左简单波 向后简单波 右简单波 向前简单波 13 简单波可以分为压缩波和稀疏波（膨胀波）两类。 设流线与 交点处，流线的切线方向为 把 式沿 求方向导数，得： 当 时 压力密度沿流线减小，且特征线 发散 稀疏波 14 当 压力密度沿流线增加，且特征线 是收敛的的:压缩波 15 4）中心稀疏波 中心稀疏波是一种特殊的简单波。以向左中心稀疏波为例 对应的特征线 是通过某一点的中心直线族。 ,则特征线的方程为： 设这一点的坐标为 设左中心稀疏波的左边界（波头）的特征线为 右边界（波尾）特征线为 由简单波定义，还有 因此， 此外， 中心稀疏波的解: 相似解: 只与x/t有关 特征线是直线 压力密度可由等熵条件得到. 16 5）激波和接触间断 当流场中存在间断时，Rankine-Hogoniot关系为： 假定间断两侧速度是连续的，则必有 这种间断称为接触间断。由Rankine-Hogoniot关系易知，接触间断两侧压力也是连续的，发生间断的只有密度。 不是接触间断的间断称为激波。 对于激波而言，必有 17 结论： 与均匀流区相邻的区域，一定是简单波区或者间断线。 间断显然可以和均匀流区相邻。 如不存在间断，必有一族特征线从均匀流区进入这一区域。 这样，相应的Riemann不变量不仅沿特征线是常数， 而且在这个区域内保持常数，因此这一区域必然是简单波区。 18 Riemann问题 所谓Riemann问题就是求解Euler方程 ? 在初值 下的解。 19 Riemann问题的解法(参考Toro的书) Riemann问题解的结构 与均匀流区相邻的只能是简单波区或者间断。 如果是简单波区，其必为中心稀疏波。 中心的初始间断影响区和均匀流区通过左波（激波或中心稀疏波）和右波（激波或中心稀疏波）分开。 双波结构是超定的 Riemann问题的解是三波结构：左波和右波之间存在一个接触间断 4波或4波以上的解欠定，不可能是稳定的物理解 20 21 采用原始变量比较方便 注意 和 对应的压力和速度相同，密度不一定相同。即： 22 Riemann问题的求解过程 解法：牛顿迭代 23 ①当 时， 和 之间