高等气体动力学第7章－超声速流动.docVIP

- 1 - - PAGE 9 - 高等气体动力学 2008.3.04 第七章 超声速流动 7 超声速流动（6学时） 7.1 二维定常无黏、无旋、等熵超声速流动的特征线法 7.2 二维定常无黏、有旋、等熵超声速流动的特征线法 7.3 三维定常无黏、等熵超声速流动的特征线法 7.4 二维非定常无黏、等熵超声速流动的特征线法 7.5 特征线数值计算的预估－校正算法 本章以三维定常无黏、等熵超声速流动的特征线法为例，讨论超声速流动的特征线理论、求解方法以及在微推偏喷管中的应用。 有关特征线的一般理论已在第4章引入，不再重复。 对其它流动情况，如二维定常无黏无旋等熵超声速流动、二维定常无黏有旋等熵超声速流动以及二维非定常无黏等熵超声速流动的特征线法，请参阅霍夫曼气体动力学。 §7－3三维定常无黏等熵超声速流动的特征线法 一、原始控制方程组 在笛卡尔坐标系中，无黏可压流的三维定常等熵流动的控制方程组由动量方程（3个）、连续方程和声速方程（代替能量方程）所组成，即 (7-3-1 式中 (7-3-2 为流体中的声速，?、p、u、v和w分别是流体的密度、压强和三个速度分量。 方程组(7-3-1 二、特征方程和相容性方程 因为五个未知量?、p、u、v和w是耦合在一起的，所以必须联立求解。 仿照含两个方程的方程组情形，作一个线性组合的微分算子，即 (7-3-3 式中?1，?2、…，?5均为待定的任意函数。 将(7-3-1 (7-3-4 以上式中偏导数的系数作分量，构造如下的矢量 (7-3-5 于是，式(7-3-4)中的各项可以看成?、p、u、v和w的梯度在相应矢量上的投影，或称它的方向导数。即(7-3-4) (7-3-6 这里表示u在矢量方向上的方向导数，等等。上式就是三维定常等熵流动的相容性方程。 现在提出这样的问题，即三维空间中是否存在这样的面，使所有的(j=1, 2, …, 5)都落在这个面上？如果这样的面存在，则称它为特征面，其法线称为特征法线，而方程(7-3-6) 特征法线是一个三维矢量，定义为 (7-3-7 因为所有的矢量，都在特征面上，故有 (7-3-8 将5个矢量(7-3-5)式代入上式，展开后 (7-3-9 式中， (7-3-10 是流体的速度矢量。 为了使方程组(7-3-9 (7-3-11 这就是确定特征面的方程。 三、特征面 1. 流面 令方程(7-3-11 (三重根) (7-3-12 因为特征法线n的实际大小并不重要，所以可将其当成单位矢量，即 (7-3-13 以上两个方程，式(7-3-12)和(7-3-13 方程(7-3-13)是一个球面方程，其球心在流线上。方程(7-3-12)则是一个过球心的空间平面方程，平面的取向取决于三个速度分量u、v和 平面方程(7-3-12 图7-3-1 三维定常流中流面的法线平面 图7-3-2 流面的包络是流线，若设t表示流体微团沿流线运动的时间，则流线的方程可写成 ，， (7-3-14) 在dt时间内，流体微团沿流线运动的微元距离矢量为 (7-3-15 2. 波面 令方程(7-3-11 (7-3-16 式中的正负号是任选的。 由此可见，确定n的条件也只有两个，即方程(7-3-16)和(7-3-13 和流面一样，方程(7-3-13)仍是一个球面方程，方程(7-3-16)仍是一个空间平面方程，但该平面不通过球心，而且平面的取向不仅取决于三个速度分量u、v和w的大小，还取决于声速a 从球心发出、并指向此圆的每一条半径代表一个特征法线，所以全部法线组成了以球心为顶点的锥，即法线锥。这说明，法线锥的每一条母线都是特征法线，而每一条特征法线又对应一个特征面，这个特征面称为波面。 显然，无数条特征法线就对应于无数个特征面，即无数个波面。这些波面的包络称为马赫劈锥，它的顶点也在流线上。与马赫劈锥在锥顶点相切的锥是马赫锥。马赫锥代表均匀流场中的马赫劈锥，所以它是一个直立圆锥，并且是法线锥的倒锥。 方程(7-3-16)表明，垂直于特征面的流动速度的大小等于声速 图7-3-3 波面的法线锥 图7-3-4 以微分形式表示的马赫劈锥方程可写成 (7-3-17 式中，。 不难验证，只有在超声速流动时，即V2-a2>0时，由上式给出的马赫劈锥才是实的。 马赫劈锥与波面相切的曲线称为双特征线。图7-3-5表示在类空间曲线上任一点P(x,y,z)上的两族波面、流面、马赫劈锥和双特征线之间的关系。类空间曲线是指马赫劈锥以外的任意曲线。如果此类空间曲线旋转360° 马赫锥与波面相切的线称为双特征射线，它是直线。显然，双特征射线是双特征线的切线。 如图7-3-6所示。若以表示流体微团沿流线在dt时间内运动的距