计算流体力学课件：双曲型方程组.pptVIP

双曲型方程组及间断解 知识点： 双曲型方程的特征方程 双曲型方程的弱解及熵条件 Riemann间断解—— 精确解、近似解初步 * * 知识回顾 1. 流体力学基本方程 概念： 连续介质假设； Euler描述/Lagrange描述 N-S方程—— 描述 质量、动量、能量守恒 的方程组 流通量： 单位时间内通过垂直于x/y/z 轴单位面积的 质量、动量、能量 无量纲量： 物理量与参考量（特征量）之比 2. 偏微分方程（组）及其类型 解耦成N个独立的方程 双曲型 有N个实特征根（含重根） N个独立特征向量 全部为复特征根 有1个N重特征根 独立特征变量数N 抛物型 椭圆型 特征线； 特征相容关系； 双曲方程边界条件提法 方法： 独立给定j个方程的边界条件 如果 lj0, 则在左端给定vj的边界条件 如果 lj0, 则在右端给定vj的边界条件 A B j=1 j=2 * 一维Euler方程 条件 描述 边界条件设定 超音速入口 给定3个边界条件 亚音速入口 给定2个边界条件 超音速出口 无需给定边界条件 亚音速出口 给定1个边界条件 * 变系数方程组的情况 令： 令 （行向量） 在x-t空间引入特征线： 1. 双曲型方程组的特征方程 * （变系数情况）虽然不能解耦，但能转换成常微方程组 2.1 双曲型方程组 x t 特征线 沿特征线 * 若不考虑粘性，流体微团运动过程中熵不变； 如果来流熵均匀分布，则全流场熵均匀分布 例：一维等（均）熵运动 预备知识： 完全气体中的热力学量 密度、压力、温度、熵、焓 内能、声速 只有两个独立变量 （完全气体）仅与温度有关 小常识： 等熵（绝热）关系 与等温相比，绝热气体更难压缩 等熵情况下，仅有一个独立的热力学变量； 给定任何一个都意味着给定全部热力学量； 矩阵A的特征值 若不考虑粘性，流体微团运动过程中熵不变； 如果来流熵均匀分布，则全流场熵均匀分布 均熵运动情况下，能量方程可用熵为常数替代 一维均熵流动控制方程（Euler方程简化版） 沿特征线1： 有： 沿特征线1： R不变 （1）转化为 x t 特征线 积分因子 定义： 注意：声速c 是温度的函数，可不是常数！ c2 ~ T ( c2 就是温度啊！） 绝热关系式 特征相容关系 * 知识点，牢记！ 一维均熵流动沿特征线Riemann不变量保持不变 x t 特征线1 特征线2 同理推导， 沿特征线2： 在（x,t）空间： 沿特征线1： 沿特征线2： A B C 扰动源 扰动向两侧传播 扰动波以当地声速向两侧传播 观测者 感受到两侧的扰动 例2.1： 有限区域内扰动波的传播问题 考虑一维无粘流动（Euler方程），初始时刻（t=0）流动状态如下： 试分析t=t0时刻的流动状态 （假设流场不出现间断） 不同时刻的速度分布（A=1） 不同时刻的速度分布（A=0.01） 思考题： 小扰动的传播情况？ 数值解 x t （1） （2） （3） （4） 利用特征线，分析不同区域的差异 等（均）熵情况下，同族特征线不会相交 * 目的： 学会如何运用Riemann不变量解题 * 一维扰动波的传播 （上： A=1; 下： A=0.01） 大扰动，非线性波 小扰动，线性波 基本解题思路： 利用特征关系 1 2 3 x t 解出 x1, x2 利用Riemann不变量得： 解出 x t （1） （2） （3） （4） 区域（2），（4） 未扰动 区域（1）内的流动使用基本方法计算 区域（3）内的计算可简化 A B D C E F G (3) 区内的波传播速度为常数，且在传播过程中物理量保持不变—— 简单波 ? 特征线为直线 因而方程是非线性的 给定x3,t3 利用 * （假设t3充分小） 解出t3时刻的流场，继续推进下个时刻 概念： 简单波 区域 （3） 内扰动波的传播特点 考虑 （3）区内的, 同属一条特征线M 上的任意两个点4 和5： 由于点1 和点3 均在未扰动区： x t （1） （2） （3） （4） 1 2 3 4 5 M 在（3）区内， 所有物理量（u,c）沿特征线M不变 特征保持直线，特征波传播速度不变 简单波 * * 各区物理含义 x t （1） （2） （3） （4） x 扰动区 t=0时刻 t=t1时刻 右行波 左行波 区域（1），感受到左、右波的影响 区域（3），仅感受到左行波的影响 简单波 区域（2），尚未感受到波 x t=t2时刻 区域（4），波已传播过去，恢复平静 波型、波速不变 3. 双曲型方程的间断解 （弱解） 双曲方