双曲型方程组省公共课一等奖全国赛课获奖课件.pptx

计算流体力学讲义计算流体力学讲义 第二讲第二讲 双曲型方程组及间断解双曲型方程组及间断解李新亮李新亮；力学所主楼；力学所主楼219；82543801 知识点：知识点：双曲型方程组边界条件提法双曲型方程组边界条件提法 双曲型方程特征方程双曲型方程特征方程 双曲型方程间断解及熵条件双曲型方程间断解及熵条件 Riemann间断解间断解1讲义、课件上传至讲义、课件上传至 (流体汉字网）流体汉字网）-“流体论坛流体论坛”-“CFD基础理论基础理论”第1页第第2 讲讲 双曲型方程组及其间断解双曲型方程组及其间断解2.4 双曲型方程及其数学性质双曲型方程及其数学性质 考虑方程组：考虑方程组：令：令：1.双曲方程边界条件提法双曲方程边界条件提法假如矩阵假如矩阵A 能经过相同变换对角化能经过相同变换对角化双曲型双曲型Copyright by Li Xinliang2第2页1)一阶常系数偏微方程组一阶常系数偏微方程组假如矩阵假如矩阵A 能够被对角化：能够被对角化：令：有即：m个方程完全个方程完全解耦，解耦，可独立求解可独立求解有有m 条特征线：条特征线：m个特征相容关系式：个特征相容关系式：假如矩阵假如矩阵A能够（相同变换）对角化，则原方程是能够（相同变换）对角化，则原方程是双曲型双曲型Copyright by Li Xinliang3第3页双曲方程边界条件提法双曲方程边界条件提法变换成为了变换成为了彼此独立彼此独立n个单波方程个单波方程方法：方法：独立独立给定给定j个方程边界条件个方程边界条件 假如假如 l lj0,则则在在左左端端给给定定vj边边界条件界条件 假如假如 l lj 难点:间断 20世纪 70年代(人工粘性，Steger-Warming,Roe)80 年代TVD,NND 90年代 ENO,WENO,AUSM,GVC 二十一世纪 WENO,AUSM+，DG b)精确解：利用空气动力学 （积分关系式）激波：激波：R-H 关系式关系式膨胀波：膨胀波：特征相容关系（特征相容关系（Riemann不变量）不变量）接触间断：接触间断：随流体传输，仅密度突变随流体传输，仅密度突变 间断条件：质量、动量、能量守恒初始值不满足间断关系，会分解成三个波独立传输初始值不满足间断关系，会分解成三个波独立传输Copyright by Li Xinliang18第18页膨胀波膨胀波 接触间断接触间断 激波激波xtSod 激波管起动后气流演化过程示意图激波管起动后气流演化过程示意图膨胀波膨胀波 接触间断接触间断 激波激波示意图普通情况：五种可能普通情况：五种可能xt激波激波 接触间断接触间断 激波激波膨胀波膨胀波 接触间断接触间断 激波激波激波激波 接触间断接触间断 膨胀波膨胀波膨胀波膨胀波 接触间断接触间断 膨胀波膨胀波 膨胀波膨胀波 膨胀波膨胀波（1）（2）（3）（4）（5）分析分析Copyright by Li Xinliang19动画演示：动画演示：密度演化密度演化第19页 2.求解方法求解方法 针对每种情况分别考虑；针对每种情况分别考虑；利用积分关系，将微分方程化成代数方程计算利用积分关系，将微分方程化成代数方程计算xt激波激波 接触间断接触间断 激波激波Zone:1 3 4 2积分关系式：积分关系式：1-3 两区两区2-4 两区3-4 两区8个方程，个方程，8个未知数。可解！个未知数。可解！其中：1)对于情况（对于情况（1）Copyright by Li Xinliang20Z1,Z2:以向右运动为正以向右运动为正第20页1-3 两区两区 积分关系式积分关系式一样，由一样，由2-4 两区关系式，可解出两区关系式，可解出其中：详细求解方法详细求解方法 Copyright by Li Xinliang21激波激波 接触间断接触间断 激波激波Zone:1 3 4 23个方程，个方程，4个未知数个未知数将压力设为已知，解出速度将压力设为已知，解出速度x（1）（2）（3）（4）利用（利用（3）(4)两式，有两式，有1个方程、1个未知数 ，可解（比如Netwon迭代法）思绪：消元法思绪：消元法 利用激波关系式，解出速利用激波关系式，解出速度对压力依赖关系度对压力依赖关系代入（代入（3）解出）解出再代入原式解出解出解出OK第21页膨胀波内部物理量计算膨胀波内部物理量计算波尾标准：标准：先计算（先计算（4）（5）区区,再计算稀疏波内部再计算稀疏波内部Copyright by Li Xinliang222)对于情况对于情况2 Sod 激波管问题属于该情况激波管问题属于该情况膨胀波区膨胀波区 接触间断接触间断 激波激波（1）（2）（3）（4）（5）未知数7个：方程：方程：2-5 两区关系式：两区关系式：4-5两区关系式两区关系式补充一个：补充一个：1-4两区关系（两区关系（Riemann不变量）不变量）：xt波头（1）（2）（3）(4)(5)x=07 个方程，个方程，7个未知数，可解个未知数，可解！详细解法与前面类似详细解法与前面类似再补充一个：再补充一个：膨胀波内等熵膨胀波内等熵依然是消元法，利用依然是消元法，利用1-4区关系区关系给出速度对压力依赖关系给出速度对压力依赖关系?为何未知数比情况1少1个？第22页Copyright by Li Xinliang23xt波头（1）（2）（3）(4)(5)x=0 x（1）（2）（4）（5）激波、膨胀波前后速度激波、膨胀波前后速度-压力依赖关系压力依赖关系可写成统一形式：可写成统一形式：左波左波（激波或膨胀波）：（激波或膨胀波）：右波（激波或膨胀波）右波（激波或膨胀波）（表示（表示（4）（）（5）区速度和压力）区速度和压力）其中：其中：激波稀疏波得到方程得到方程：（8）1 个方程，个方程，1个未知数，可解个未知数，可解求解求解(8)得到得到4,5两区压两区压力力然后，解出速度和密度然后，解出速度和密度第23页膨胀波内部物理量计算膨胀波内部物理量计算xtxt波头波尾处理方法：处理方法：1）计算膨胀波范围计算膨胀波范围 波头传输速度波头传输速度 波尾传输速度波尾传输速度（1）（2）（3）(4)(5)2）在膨胀波区内，利用特征相容关系计算）在膨胀波区内，利用特征相容关系计算 利用简单波特征，简化计算利用简单波特征，简化计算简单波简单波x=0特征线由特征线由x=0发出发出再利用另一条特征线信息：再利用另一条特征线信息：解出再利用等熵关系，计算Copyright by Li Xinliang24膨胀波膨胀波 接触间断接触间断 激波激波第24页Copyright by Li Xinliang25求解步骤求解步骤 针对情况针对情况1，2x（1）（2）（4）（5）xt（1）（2）（3）(4)(5)x=0step1.求解方程求解方程(8)，解出解出 4，5区压力区压力 单未知数代数方程；数值方法求解单未知数代数方程；数值方法求解其中：step 2.求出求出4，5区速度、密度、激波移动速度区速度、密度、激波移动速度 step 3.计算出稀疏波区量计算出稀疏波区量 针对情况1，求解完成；对于情况2 继续step 3其中各区范围以下（以情况2 讨论）：1 区：3 区：4 区：5 区：2 区：以上步骤完全适合用于以上步骤完全适合用于 情况情况3，4，5 （因为（因为8式同时适合用式同时适合用于激波和稀疏波）！于激波和稀疏波）！第25页 Riemann 问题五种可能情况问题五种可能情况xt激波激波 接触间断接触间断 激波激波膨胀波膨胀波 接触间断接触间断 激波激波激波激波 接触间断接触间断 膨胀波膨胀波膨胀波膨胀波 接触间断接触间断 膨胀波膨胀波 膨胀波膨胀波 膨胀波膨胀波（1）（2）（3）（4）（5）怎样区分这怎样区分这5种情况？种情况？Copyright by Li Xinliang26假设准则以下：情况1情况3情况4利用函数 （由 8式定义）函数性质很好函数性质很好 单调连续单调连续情况5 思索题?为何没有出现情况 2 情况5情况4情况3情况1Riemann 求解总步骤：1）依据上述判决区分情况 2）按照上一页步骤求解 真空区真空区第26页Copyright by Li Xinliang27思索题：思索题：上述求解方法要求间断两侧流场分布为常数，假如初上述求解方法要求间断两侧流场分布为常数，假如初始时刻流场分布是始时刻流场分布是x函数，怎样利用该理论解计算函数，怎样利用该理论解计算？提醒：提醒：把曲线离散化，看成折线把曲线离散化，看成折线第27页作业：作业：以下以下Sod 激波管问题激波管问题:求出理论解，求出理论解，并分别画出并分别画出t=0.14时刻时刻 分布曲线。分布曲线。Copyright by Li Xinliang28知识点，知识点，主要练习，熟练掌握！主要练习，熟练掌握！第28页