高等气体动力学第4章－可压缩一维非定常流动（老版本-左特征矩阵符号未改）.doc

高等气体动力学 2008.2.19 第四章 可压缩一维非定常流动 §4－1 可压缩非定常流动的控制方程组－欧拉方程组 §4－2 一维非定常欧拉方程组的特征分析 §4－3 双曲型守恒律的弱解 §4－4 一维非定常流动的一般特征－微弱扰动波的传播 §4－5 一维非定常流动的一般特征－微弱扰动波的反射与相交 §4－6 一维非定常流动的一般特征－有限强度扰动波的传播 §4－7 特征线法的一般理论 §4－8 一维非定常均熵流动的特征线理论 §4－9 一维激波管（Riemann）问题 §4－1 可压缩非定常流动的控制方程组－欧拉方程组 欧拉方程组是描述无黏流的最高级或说是最完整的方程组，其解是N-S方程组当黏性系数趋于零时的极限； 欧拉方程组描述了N-S方程组的对流性质。在高雷诺数下，N-S方程组是由对流控制的，这也是大多数工程与自然界中的流动情形； 因此，研究Euler方程组就显得非常重要。从数值模拟角度而言，一般情况下能够适用于Euler方程组的数值方法对N-S方程组也同样有效，这时，只需要将黏性项按中心差分离散即可。只有对低雷诺数流动，由于它受黏性耗散控制，才需要使用特殊的方法； Euler方程组的物理意义在于它表述了质量、动量和能量的守恒定律，所以推导时的基本出发点是守恒形式的控制方程组； 在无黏流中可能存在各种间断，就数值模拟而言，为了在数值模拟中正确捕捉这些间断，必须使用守恒形式的控制方程组。只有当确信不会出现间断时，才可以使用非守恒方程组； Euler方程组是一阶非线性耦合方程组，可以写成各种不同的等价形式，这取决于流动变量的选取。在本课程中，定义了如下的因变量： ① 守恒变量——是服从守恒方程的变量，定义为 ② 原始变量——使用更为自然、方便，在实验中可以直接测量或控制的变量，定义为 ，，，， 在计算中，用这样的变量赋边界条件较直接和方便。 ③ 特征变量——由于Euler方程组在时间上是双曲型的，即存在特征线（或特征面），那些沿着特征线（特征面）传播的变量可定义为特征变量。 一、欧拉方程组的守恒形式 为方便起见，将欧拉方程组的各种形式再写一遍。 积分守恒形式 (4-1-1) (4-1-2) (4-1-3) 式中：－作用在控制体上的外力； －单位质量流体总能量； －单位质量流体总焓； 为将守恒方程写成紧凑形式，令 ，　， 则积分守恒方程组可写成 (4-1-4) 说明 列矢量包含的变量是守恒变量； 包含的是对流通量； 要求解上式，需补充描述流体热力学性质的方程，即状态方程； 一般地，可将内能定义 　　或　　 则状态方程的形式为 理想气体：内能；状态方程。 微分守恒形式 (4-1-5) (4-1-6) (4-1-7) 写成矩阵形式 (4-1-8) 或 (4-1-9) 式中 ，， (4-1-10) 二、欧拉方程组的拟线性形式 为了研究欧拉方程组的性质，必须将它写成拟线性形式； 欧拉方程只包含一阶导数。如果外力与流动变量无关，则是一阶偏微分方程组。 1. 守恒变量的Jacobian矩阵 定义 欧拉方程组(4-1-8)可以写成拟线性形式 (4-1-11) 矢通量对守恒变量的Jacobian矩阵定义为 (4-1-12) 将式(4-1-11)改写成 (4-1-13) 或显式地写出 (4-1-14) 各方向通量的Jacobian矩阵可写成 ，， (4-1-15) Jacobian矩阵的齐次性质 对理想气体，或更一般地，对具有如下性质 (4-1-16) 的流体，各方向通量是守恒变量的一阶齐次函数，即对任意的，有 (4-1-17) 对求偏导有 令有 (4-1-18) 或写成分量形式 ，　，　 (4-1-19) 将上式代入式(4-1-8)或(4-1-9)可得 (4-1-20) 将上式与式(4-1-14)相比，可见Jacobian矩阵既可写在微分外，也可写在微分内，它们在数学上是完全等价的。但数值计算时并不是一回事，它们将导致完全不同的离散形式，因为式(4-1-14)是非守恒形式，而式(4-1-20)是守恒形式。 Jacobian矩阵的具体形式 为求Jacobian矩阵，将和写成守恒变量的函数，令守恒变量为 ，　，　 其中 ，　，　，　 用上式可将和写成如下形式 ，，，(4-1-21) 在推导之前，还需将表示成的函数。对理想气体，有 即 (4-1-22) 用以上表达式可求出Jacobian矩阵，例如 (4-1-23) 式中每一项都代表一个有五个元素的列矢量，如 (4-1-24) 由式(4-1-22)得 (4-1-25) 所以有 (4-1-26) 将(4-1-23)中的每一项都求出，并还原，，和的符号后可得 (4-1-27) (4-1-28) (4-1-29) 可见，Jacobia