清华工程流体力学课件第七章气体一维高速流动

1、2021-6-13清华工程流体力学课件第七章气体一维高速流动 第七章第七章 气体一维高速流动气体一维高速流动 u 第一节 微弱扰动波的传播 u 第二节 气体一维定常等熵流动 u 第三节 气体一维定常等熵变截面管流 u 第四节 正激波 2021-6-13清华工程流体力学课件第七章气体一维高速流动 前几章讨论的是不可压缩流体的流动，例如对于液体， 即使在较高的压强下密度的变化也很微小，所以在一般情况 下，可以把液体看成是不可压缩流体。对于气体来说，可压 缩的程度比液体要大得多。但是当气体流动的速度远小于在 该气体中声音传播的速度（即声速）时，密度的变化也很小。 例如空气的速度等于50m/s，这数值

2、比常温20下空气中的 声速343m/s要小得多，这时空气密度的相对变化仅百分之一。 所以为简化问题起见，通常也可忽略密度的变化，将密度近 似地看作是常数，即在理论上把气体按不可压缩流体处理。 当气体流动的速度或物体在气体中运动的速度接近甚至超过 声速时，如果气体受到扰动，必然会引起很大的压强变化， 以致密度和温度也会发生显著的变化，气体的流动状态和流 动图形都会有根本性的变化，这时就必须考虑压缩性的影响。 气体动力学就是研究可压缩流体运动规律以及在工程实际中 应用的一门科学。本章中仅主要讨论气体动力学中一些最基 本的知识。 2021-6-13清华工程流体力学课件第七章气体一维高速流动 第一节

3、微弱扰动波的传播 一一. 微弱扰动波的一维传播微弱扰动波的一维传播 如图7-1所示，在一个截面积为A、足够长的直圆管中充满 了静止的气体，将圆管左端的活塞以微小速度 向右轻微地 推动一下，使活塞右侧的气体压强升高一个微小增量 ， 所产生的微弱压强扰动向右传播。活塞将首先压缩紧贴活塞的 那一层气体，这层气体受压后，又传及下一层气体，这样依次 一层一层地传下去，就在圆管中形成一个不连续的微弱的压强 突跃，就是压缩波mn，它以速度 向右推进。压缩波面mn是 受活塞微小推移的影响而被扰动过的气体与未被扰动过的静止 气体的分界面。设在压缩波前未被扰动过的静止气体的压强 为 、密度为 、温度为 ，波后已被

4、扰动过的气体以与 活塞的微小运动同样的微小速度 向右运动，其压强增高 到 ，密度和温度也相应增加到 和 。 Vd pdpd p T Vd ppdd TTd 2021-6-13清华工程流体力学课件第七章气体一维高速流动 图7-1 微弱扰动波的一维传播 2021-6-13清华工程流体力学课件第七章气体一维高速流动 显然，这是不定常流动。为了得到定常流动，可以 设想观察者随波面mn一起以速度c向右运动。气体相对 于观察者定常地从右向左流动，经过波面速度由c降为 c-dv，而压强由p升高到p+dp，密度和温度由 、 增 加到 、 。如图7-1(b)所示，取包围压缩波的 控制面，根据连续性条件，在 时间

5、内流入和流出该 控制面的气体质量应该相等，即 化简后，得 （7-1） 由于压缩波很薄，作用在该波上的摩擦力可以忽略不计。 于是对于控制面，根据动量定理，沿气体流动的方向，质 量为 的气体的动量变化率等于作用在该气体上的压力 之和，即 或 （7-2） T dTT d t d tAVctAcd)d)(d(d d d d c V Ac Appp t cVc tAc)d( d )()d( d p c Vd 1 d 2021-6-13清华工程流体力学课件第七章气体一维高速流动 由式（7-1）和式（7-2）得 由于是微弱扰动， 远小于 ，即 ，所以 （7-3） 式（7-3）与物理学中计算声音在弹性介质中传

6、播速度 （即声速）的拉普拉斯公式完全相同。可见气体中微弱扰 动波的传播速度就是声速。 在式（7-3）的推导过程中，并未对介质提出特殊要求，故 该式既适用于气体，也适用于液体，乃至适用于一切弹性 连续介质。不同介质的压缩性不同，压缩性小的扰动波传 播速度高，压缩性大的扰动波传播速度低，因此声速值反 映了流体可压缩性的大小。 式(7-3)是声速的通用表达式，要计算某种流体中具有的声 速值，尚需确定 和 的关系，以求出 的值。 d dd 1 2 p c d 1 d d d p c pdd d dp 2021-6-13清华工程流体力学课件第七章气体一维高速流动 由于微弱扰动波的传播过程进行得很迅速，与

7、外界来 不及进行热交换，而且其中的压强、密度和温度变化极为 微小，所以这个传播过程可以近似地认为是一个可逆的绝 热过程，即等熵过程。假定气体是热力学中的完全气体， 则根据等熵过程关系式 =常数和完全气体状态方 程 ，可得 代入式（7-3），得 （7-4） p RTp RT pp d d RT p c 为绝热指数 为气体常数，J/(kgK) 为热力学绝对温度，K 对于空气， , R= 287 J/(kgK)。4.1 2021-6-13清华工程流体力学课件第七章气体一维高速流动 由式（7-4）可知，气体中的声速随气体的状态参数 的变化而变化。于是在同一流场中，各点的状态参数若 不同，则各点的声速也

8、不同。所以声速指的是流场中某 一点在某一瞬时的声速，称为当地声速。 在实际计算中，通常用气体速度 与当地声速 的比值 来作为判断气体压缩性对流动影响的一个标准，即 （7-5） 称为马赫数，是一个无量纲数，也是气体动力学中一 个重要参数。 我们常根据马赫数的大小,把气流分为亚声速流 1， 跨声速流 1，超声速流1 3等 几类。亚声速流动和超声速流动有许多显著的差别，我 们将在以后各节中逐一介绍。 VcMa Ma Ma Ma Ma Ma c V Ma 2021-6-13清华工程流体力学课件第七章气体一维高速流动 二二 微弱扰动波的空间传播微弱扰动波的空间传播 前面讨论了微弱扰动波的一维传播，下面进

9、一步讨论 微弱扰动波在空间流场中的传播。 为了便于分析问题，假设流场中某点有一固定的扰动源， 每隔1s发生一次微弱扰动，现在分析前3s产生的微弱扰动 波在空间的传播情况。由于不论流场是静止的还是运动的， 是亚声速的还是超声速的，都将对微弱扰动波在空间的传 播情况产生影响，所以下面分四种情况来讨论。 1静止流场(V=0) 在静止流场中，扰动源产生的微弱扰动波以声速c向四周 传播，形成以扰动源所在位置为中心的同心球面波，微弱 扰动波在3s末的传播情况如图7-2(a)所示。如果不考虑微弱 扰动波在传播过程中的损失，随着时间的延续，扰动必将 传遍整个流场。也就是说，微弱扰动波在静止气体中的传 播是无界

10、的。 2021-6-13清华工程流体力学课件第七章气体一维高速流动 2亚声速流场(Vc) 在亚声速流场中，扰动源产生的微弱扰动波在3s 末的传播情况如图7-2(b)所示。由于扰动源本身以 速度运动，故微弱扰动波在各个方向上传播的绝对 速度不再是当地声速c，而是这两个速度的矢量 和。这样，球面扰动波在顺流和逆流方向上的传播 就不对称了。但是由于Vc） 在超声速流场中，扰动源产生的微弱扰动波在3s末的 传播情况如图7-2(d)所示。由图可见，由于vc，所以相 对气流传播的扰动波不仅不能向上游传播，反而被气流带 向扰动源的下游，所有扰动波面是自扰动源点出发的圆锥 面的一系列内切球面，这个圆锥面就是马

11、赫锥。随着时间 的延续，球面扰动波不断向外扩大，但也只能在马赫锥内 传播，永远不会传播到马赫锥以外的空间。也就是说，微 弱扰动波在超声速气流中的传播也是有界的，界限就是马 赫锥。 2021-6-13清华工程流体力学课件第七章气体一维高速流动 马赫锥的半顶角，即圆锥的母线与气流速度方向之间 的夹角，称为马赫角，用 表示。由图7-2(d)可以容易地 看出，马赫角 与马赫数 之间的关系为 （7-6） 马赫角从90这时相当于扰动源以声速V=c流动的情况， 如图7-2(c)所示 开始，随着马赫数的增大而逐渐减小。由 于圆锥顶就是扰动源，所以当物体以超声速运动时，它所 引起的扰动不能传到物体的前面。马赫锥

12、外面的气体不受 扰动的影响，微弱扰动波的影响仅在马赫锥内部，即微弱 扰动波不能向马赫锥外传播。这就说明了，为什么以超声 速飞行的弹丸在附着于它头部的波未到达观察者的耳朵以 前听不到声音的原故。 Ma MaV c1 sin 2021-6-13清华工程流体力学课件第七章气体一维高速流动 上述关系也适用于气流流过一静止微小障碍 物时的情况。假如气体以与上述扰动源的运动速 度数值相等而方向相反的速度作等速直线运动， 则扰动源就成为静止微小障碍物，即图7-2中的3 点就是静止扰动源，而扰动源所发出的扰动波 （图中的各圆）不断地被气流以速度-V带走。很 明显，在 （即 ）的亚声速流动时，带 走的各扰动波在

13、一定时间后可达到空间中的任何 一点。也就是说，扰动波不仅能顺流传播，而且 也能逆流传播。但在 （即 ）的超声速 流动时，带走的各扰动波只能在马赫锥内顺流传 播，不能逆流传播，也就是说在超声速流动中的 微弱扰动不能传播到整个空间。这就是超声速流 动和亚声速流动的一个重要差别，从而使这两种 流动的图形有着根本的不同。 cV 1Ma cV 1Ma 2021-6-13清华工程流体力学课件第七章气体一维高速流动 第二节 气体一维定常等熵流动 在讨论不可压缩流体流动时，应用连续性方程和伯努 利方程就可以对许多问题求解。但是对于可压缩流体 气体流动仅仅应用上面两个基本方程还不足以求解，因为 由于气体密度的变

14、化必然会引起热力学状态发生相应的变 化。就是说在气流流动中，不仅它的力学状态在发生变化， 而且热力学状态也在随着改变。因此必须把热力学中的状 态方程和过程方程一并考虑，才能解决气体流动问题。 本节将只讨论气体的一维定常等熵流动，即假定气体 是完全气体，在流动过程中与外界无热交换，摩擦影响很 小可以忽略不计。在一般情况下还认为各参数仅在一个方 向上有显著的变化，而且变化是连续的、不随时间而变化， 这就是一维定常等熵流动。在许多实际流动问题中，例如 气体在喷管、扩压管和短叶栅中的流动都可以近似地认为 是一维定常等熵流动。 2021-6-13清华工程流体力学课件第七章气体一维高速流动 一、气体一维定

15、常流动的基本方程一、气体一维定常流动的基本方程 1连续性方程 由于气体的密度在流动中是发生变化的，所以它的连续 性方程不能像不可压缩流体那样按体积流量来计算，而需 要用质量流量来计算，即气体在流管中流动时，每单位时 间内流过流管中任意两个有效截面的质量流量必定相等， 即 (7-7) 也可以把连续性方程写成微分形式，即对式(7-7)取对数后 微分，得 (7-8) 常数VA AVAV 222111 0 ddd A A V V 2021-6-13清华工程流体力学课件第七章气体一维高速流动 2能量方程 由于气体的密度很小，所以质量力可以忽略不计。气 体是一维定常流动，并令 ， ，则欧拉运动微 分方程可

16、写成 或 （7-9） 将式（7-9）沿流管（或流线）进行积分，得 对于等熵流动，将等熵过程关系式 常数，代入上式， 得完全气体一维定常等熵流动的能量方程为 （7-10） 显然，这个方程只能用于可逆的绝热流动。 Vu 0wv x p x V V d d1 d d 0d 1 dpVV 常数 2 d 2 Vp p 常数 21 2 Vp 2021-6-13清华工程流体力学课件第七章气体一维高速流动 热力学第一定律用于流体流动的能量关系式为 在绝热流动的条件下， ，上式可写成 ，积分 可得能量方程的另一表达式 （7-11） 这个方程可用于可逆的绝热流动，也可用于不可逆的绝热 流动，即式（7-11）在熵有

17、增加（有摩擦或其他不可逆因 素）的绝热流动中也是正确的。因为在与外界无热交换的 绝热过程中，消耗于抵抗摩擦所作的功完全转换为热能， 该热能重又加入气流中，使气流中的熵增加。所以在绝热 流动中总能量不变，摩擦损失的存在只会使气流中不同形 式的能量重新分配，即一部分机械能不可逆地转化为热能， 因而能量方程（7-11）的形式不变。 VVhqddd 0d q0ddVVh 常数 2 2 V h 2021-6-13清华工程流体力学课件第七章气体一维高速流动 对于完全气体，存在下列关系 代入式（7-11），也可得到与式（7-10）同一形式的完全 气体一维定常等熵流动的能量方程。现在来分析一下这个 方程中各项

18、的物理意义，可将式（7-10）改写成 （7-12） 根据热力学可知，对于完全气体上式第一项是单位质量气 体所具有的内能u，即 pp cc c p R c Tch Vp pp p 1 uTc p R cp cc cp V V Vp V 1 常数 21 2 Vpp 2021-6-13清华工程流体力学课件第七章气体一维高速流动 而式（7-12）的后两项是单位质量气体的压强势能和动能。 所以完全气体一维定常等熵流动的能量方程的物理意义是： 在完全气体一维定常等熵流动中，气流流管任一有效截面 （或流线的任一点）上单位质量气体的压强势能、动能和 内能之和保持不变。 由于 ，代入式（7-10）得到完全气体能

19、量方程 的又一个表达式 （7-13） 2 c p 常数 21 22 Vc 2021-6-13清华工程流体力学课件第七章气体一维高速流动 二、滞止参数二、滞止参数 在实际工程上，为了分析和计算流动问题方便起见， 常使用滞止参数这个概念，而且由于它比较容易测量，所以 滞止参数得到广泛的应用。设想气体流过流管的两个有效截 面时，在一个截面上完全滞止下来，也就是说，在这个截面 上的气流速度等于零。则这个截面上的气流状态称为滞止状 态，滞止状态下各相应参数称为滞止参数，分别 以 、 、 、 等表示之。气体绕过一个物体时，在驻点处气 流受到阻滞，速度等于零，这一点的气流状态也是滞止状 态。在滞止状态下式（

20、7-10）、式（7-11）和式（7-13）可 写成 （7-14） （7-15） （7-16） 0 p 0 0 T 0 c 常数 00 0 0 2 1121 TcRT pVp p 常数 0 2 2 h V h 常数 121 2 0 22 cVc 2021-6-13清华工程流体力学课件第七章气体一维高速流动 由式（7-14）和式（7-15）可知，在滞止状态下气流 的动能全部转变为热能，可以用滞止焓 表示之， 它表示单位质量的气流所具有的总能量，称为总焓。 式（7-15）又可改写成 （7-17） 上式表明，滞止温度要比气流的温度T高出 ，对于 J/（kgK）的空气，则高出 0 2 2 T c V T

21、 p 00 Tch p p c V 2 2 1005 p c 2010 2 0 V TTT 2021-6-13清华工程流体力学课件第七章气体一维高速流动 例如速度为100m/s的空气流，滞止温度超过气流的温度约 5K，也即约5。可见，将一个带小玻璃球的普通水银温 度计或热电偶温度计放在气流中来测量气流的温度，读出 的温度比气流的温度T要高。但小玻璃球上驻点处的温度虽 达到滞止温度，但其上的其他各点的温度升高要小一些， 所以普通水银温度计上读出的平均温度比滞止温度稍低一 些。因此用任何静止温度计都不能直接测得气流的真实温 度了，只有用与气流同样速度运动的温度计才能直接测得 利用关系式 和 可将式

22、（7-17） 改写为 或 （7-18） 对于等熵气流 和 Rc p 1 RT V c V Ma 2 2 2 2 0 2 2 1 1TMaT 20 2 1 1Ma T T 1 00 T T p p 1 1 00 T T 2021-6-13清华工程流体力学课件第七章气体一维高速流动 将式（7-18）代入上两式，得 （7-19） （7-20） 这样，只要知道气流的滞止参数和 值，就可由式（7- 18）、式（7-19）和式（7-20）以及 ，求得流管 内气流在某指定截面上的温度 、压强 、密度 和速 度 。反之，若已知 、 、 和 也可求得 、滞止参 数 、 和 。所以这三个公式是计算气体一维定常等熵

23、流 动问题的基本公式。 1 20 2 1 1 Ma p p 1 1 20 2 1 1 Ma Ma RT V Ma 2 2 T p V Tp VMa 0 T 0 p 0 2021-6-13清华工程流体力学课件第七章气体一维高速流动 第三节 气体一维定常等熵变截面管流 一、气流速度与密度的关系一、气流速度与密度的关系 由一维流动的运动微分方程式 （7-9） 得 变形得 （7-21） 由式（7-21）和能量方程式（7-9）可看出： （1）不管Ma1，或 Ma 0，则 0， 0； 反之 0， 0。 这说明加速气流（ 0），必引起压强降低（ 0）和气 体膨胀（ 0）；而减速气流（ 0）和气体压缩（ 0）

24、，即气体流动伴随着密度的 变化。亚声速气流和超声速气流都具有上述性质，但当 不 同时， 与 的变化值不同。 0d 1 dpVV dd d d d 2 c p VV V V Ma V V c V c VVdddd 2 2 2 2 Vdpdd Vdpdd Vd Vd pd pd d d V Vd d Ma 2021-6-13清华工程流体力学课件第七章气体一维高速流动 （2）Mal时密度相对变化量是小于速度的相对变 化量，即 1时，密度的相对变化量大于速度的相对 变化量，即 。 这种亚声速和超声速在变化数量上的差别，导致 了亚声速和超声速在速度与通道截面形状关系上 本质的差别。 d V Vd d V

25、 Vd 2021-6-13清华工程流体力学课件第七章气体一维高速流动 二、气流速度与通道截面的关系二、气流速度与通道截面的关系 由一维流动的运动微分方程（7-9），并利用 和 得 以上两式与连续性方程的微分形式（7-8）各联立一次， 消去 和 ，得到气体的压强变化率和速度变化率与通 道截面变化率的两个关系式，即 （7-22） （7-23） p MaMacV d d 2 p c 0 d1ddd 22 p p MaV V V p V V 0 d1dd d d1d 22 MaV Vp VV V V Vd d A A Ma Ma p pd 1 d 2 2 A A MaV Vd 1 1d 2 2021-

26、6-13清华工程流体力学课件第七章气体一维高速流动 （1）Ma1，超声速流动。 与 异号，而 与 同 号。当压强降低时，通道截面积随着气流速度的增加而扩 大，这就是超声速喷管。这是由于超声速气体在压强下降 时，密度剧烈减小、体积迅速增大，这时通道截面积必须 扩大，才能使剧烈膨胀的加速气流通过。反之，当压强升 高时，通道截面积随着气流速度的减小而缩小，这就是超 声速扩压管。 由式（由式（7-22）和式（）和式（7-23）可以得到三个重要结论：）可以得到三个重要结论： p pd p pd A Ad A Ad A Ad A Ad V Vd V Vd 2021-6-13清华工程流体力学课件第七章气体一

27、维高速流动 （3） =1，这时 。从以上两种情况知道，当降压加 速的气流由亚声速连续变为超声速时，通道截面先收缩后扩 大，在最小截面（ ）处速度达到声速（ ），该最 小截面称为临界截面，也称为喉部截面，简称喉部。当升压 减速的气流由超声速连续地变为亚声速时，通道截面也是先 收缩后扩大，在最小截面处速度达到声速。 在临界截面上的相应参数称为临界参数，分别以 、 和 等表示之。可将 代入式（7-18）、式（7-19）和 式（7-20），得到临界截面上气流的临界温度、临界压强和 临界密度各与滞止温度、滞止压强和滞止密度之间的关系式 （7-24） （7-25） （7-26） 常用气体的物理性质见表7-

28、2。 Ma 0dA cV 0dA T p 1Ma 0 1 2 TT 0 1 1 2 pp 0 1 1 2 2021-6-13清华工程流体力学课件第七章气体一维高速流动 表表7-1 气流参数变化与通道截面变化之间关系气流参数变化与通道截面变化之间关系 2021-6-13清华工程流体力学课件第七章气体一维高速流动 表表7-2 常用气体的物理性质（标准大气压强、常用气体的物理性质（标准大气压强、20） 2021-6-13清华工程流体力学课件第七章气体一维高速流动 三、气体经渐缩喷管和缩放喷管的流动三、气体经渐缩喷管和缩放喷管的流动 由上面可知，要使气流加速，当流速尚未达到当地声 速时，喷管截面应逐渐

29、收缩，直至流速达到当地声速时， 截面收缩到最小值，这种喷管称为渐缩喷管。渐缩喷管 出口处的流速最大只能达到当地声速。要使气流从亚声 速加速到超声速，必须将喷管做成先逐渐收缩而后逐渐 扩大形（在最小截面处流速达到当地声速），这种喷管 称为缩放喷管。缩放喷管是瑞典工程师拉伐尔（de Laval）在研制汽轮机时发明的，所以又称为拉伐尔喷管。 这种利用管道截面的变化来加速气流的几何喷管，在汽 轮机、燃气轮机、喷气发动机和流量测量中被广泛地应 用，本节以完全气体为对象，来讨论渐缩喷管和缩放喷 管基本设计关系式。 2021-6-13清华工程流体力学课件第七章气体一维高速流动 1、渐缩喷管 假定气体在等熵条

30、件下从大容器中经渐缩喷管流出，如 图7-3所示。由于容器的容量很大，可近似地把容器中的 气体速度看作是零（ ），即容器中的气体处于滞止 状态（ 、 、 ），而喷管出口截面上的气流参 数为 、 和 。对0-0，2-2截面列一维定常等熵流动的 能量方程（7-10），得 或 0 0 V 0 p 0 0 T 2 p 2 2 T 211 2 2 2 2 0 0 Vpp 2 0 0 2 0 0 2 1 1 2 p pp V 2021-6-13清华工程流体力学课件第七章气体一维高速流动 图7-3 渐缩喷管 2021-6-13清华工程流体力学课件第七章气体一维高速流动 将等熵过程关系式 代入上式，得出口截面处

31、的流速为 或 又 ，则出口截面上的马赫数为 （7-29） 2 2 0 0 pp 或 1 2 0 2 0 p p 1 0 2 0 0 2 1 1 2 p pp V 1 1 2 1 2 0 2 2 2 p pp V 2 2 2 p c 1 1 2 1 2 0 2 2 2 p pp Ma 2021-6-13清华工程流体力学课件第七章气体一维高速流动 222 AVq m 2 A 通过喷管的质量流量 式中 喷管出口截面积。 将 和式（7-27）代入上式，得 (7-30） 1 0 2 02 p p 1 0 2 0 0 2 1 0 2 0 1 1 2 p pp A p p qm 1 0 2 2 0 2 0

32、0 20 1 2 p p p pp A 2021-6-13清华工程流体力学课件第七章气体一维高速流动 由上式可知，当气体的滞止参数和喷管的出口截面积保持不 变时.质量流量 仅随压强比 而变化，由式（7-30） 描绘出的 与 的关系曲线如图7-4（a）所示。 m q 0 2 p p m q 0 2 p p 图7-4 气体流过渐缩喷管时流量与设计出口压强和环境压强的变化 2021-6-13清华工程流体力学课件第七章气体一维高速流动 当气体经过设计成的渐缩喷管时，实际上质量流量 随着 而变化, 为喷管出口截面外的气流压强， 称为环境压强。 与 的关系曲线如图7-4（b） 所示， 与图7-4（a）中的

33、 - 曲线相比，两者有 明显的差异。从中我们得到如下结论： m q 0 amb p p amb p m q 0 amb p p m q （1） = =1到 最大值时对应的压强与 相比，两曲线ab完全吻合。 0 2 p p 0 amb p p m q （2） =1时， =0，即当喷管的进、出口压强相等 时，气体不流动， 出口马赫数 。 0 amb p p m q 0 2 Ma 0 2 p p m q 2021-6-13清华工程流体力学课件第七章气体一维高速流动 （3） l时， 逐渐降低，出口马赫数 逐渐增加， 沿曲 线ab逐渐增加，当出口截面上的流速增加到声速时，即 时，流量达到最大值 ，此时的

34、压强比 称为临界压强比 。 可由下法求得，即将 代入式（7-29），得临界压强比 （7-31） 也可通过直接对式（7-30）求导，并令 的方法求得。 再将式（7-31）代入式（7-30）和式（7-27）中，即分别 得到临界流量 ，也就是最大流量 和临界速度 （也称 为临界声速 ）. （7-32） （7-33） 0 amb p p 0 amb p p 2 Ma 1 2 Ma m q m q 0 amb p p 1 2 Mamax q 1 00 2 1 2 p p p p p0 d d 2 p qm m q V c maxm q )1(2 1 002max 1 2 pAqq mm 0 0 1 2

35、p cV 2021-6-13清华工程流体力学课件第七章气体一维高速流动 （4） 从 再继续降低，即 ，则在整个喷管内部是亚声速气 流，如图7-5中ODE曲线所示，这时的缩放喷管相当于文 丘里管。 amb p 2 p 0 p amb p 0 p amb p 0 p amb p 0 p 1 p p1 p p 2021-6-13清华工程流体力学课件第七章气体一维高速流动 图7-5 缩放喷管内的压强和流量变化 返回(2) 返回(4) 2021-6-13清华工程流体力学课件第七章气体一维高速流动 （3）如果环境压强 继续下降到使最小截面上的压强达 到临界压强 ，则流量达到如式（7-34）所示的最大值 这

36、时在喷管扩张部分可能有两种流动状况：当 = （ 为喷管中气流只在喉部达到声速其余全为亚声速时出口 截面的压强）时，在整个喷管扩张部分中仍然都是亚声速气 流，如图7-5中OCF曲线所示；而当环境压强等于喷管出口 截面上的设计压强时，即 = ，在整个喷管的扩张部分 中都是超声速气流，如图705中OCJ曲线所示，即气流在缩 放喷管内压强从 下降到 （即亚声速连续变到超声速） 的连续变化曲线。 amb p p maxm q amb p m p 2 m p 2 amb p 2 p 2 p 0 p 2021-6-13清华工程流体力学课件第七章气体一维高速流动 （4）当环境压强在 和 之间，即 ，气流在扩

37、张部分会出现压强的不连续变化，也就是形成一个所谓正 激波，正激波的位置随着 的下降，从最小截面处移到 喷管出口处， 就是正激波移到喷管最小截面时的出口 压强，气流通过正激波从超声速变成亚声速，一直到出口 截面处，如图7-5中OCS1S2H线所示。 显然，对于缩放喷管，只要 ，不论环境压强 怎 样变化，气流通过缩放喷管的流量将始终保持为最大流量， 这是由于喷管最小截面处的临界参数没有变化；当 时，流量将减小；当 = 时，流量等于零. 如图7-5的右图所示。 m p2 m p 2 m p2 m p2 m p2 2 p 2 p amb p amb p amb p amb p amb p amb p

38、0 p 2021-6-13清华工程流体力学课件第七章气体一维高速流动 第四节 正 激 波 一、正激波形成一、正激波形成 本节以气体中的微弱扰动波在直圆管中传播的情况为例来 说明正激波形成的物理过程。 如图7-6所示，在一个充满静止气体的直圆管中，活塞向 右作加速运动，活塞右侧的静止气体受压后被扰动形成一 个压缩波向右移动，已被扰动的气体的压强从 升高 到 ，设 - 是一个有限的压强量。为了分析方便起见， 假定把这个有限的压强增量看作是无数个无限小压强增量 dp的总和。于是，可认为在活塞右侧形成的压缩波是一系 列微弱扰动波连接而成的。每一个微弱扰动波压强增加dp。 当活塞开始运动时，第一个微弱扰

39、动波以声速 传到未被 扰动的静止气体中去，紧跟着第二个微弱扰动波以声速 传到已被第一个微弱扰动波扰动过的气体中去。 1 p 2 p2 p 1 p 1 c 2 c 2021-6-13清华工程流体力学课件第七章气体一维高速流动 图7-6 在圆管中正激波的形成过程 2021-6-13清华工程流体力学课件第七章气体一维高速流动 显然，被第二个微弱扰动波扰动过的气体中的压强、密度 和温度都比被第一个微弱扰动波扰动过的气体中的相应参 数略大一些，根据 ，因此 ，也就是说第二个微 弱扰动波的声速比第一个微弱扰动波的声速略快一些。与 此相类似，第三个微弱扰动波又以比第二个略快一些的声 速（ ）向右传播，。如果

40、在某一时刻波形如图7-6 （a）所示，经过一段时间后，后面的微弱扰动波一个一 个追赶上前面的波，波形变得愈来愈陡，最后叠加成一个 垂直于流动方向的具有压强不连续面的压缩波，这就是正 激波，如图7-6（c）所示。气流通过激波除压强突跃地升 高外，密度和温度也同样突跃地增加，而速度则下降。发 生这种突跃地不连续变化是在与气体分子平均自由行程同 一数量级（在空气中约310-4mm左右）内完成的，也可 以说，各气流参数是在一个极小的激波厚度内连续地进行 变化的。当然也可以认为，是在一个几何面上突然变化的； 这就是说，可以把激波看作是一个不连续的间断面，气流 参数通过激波的变化是突跃的，不连续的。 RT

41、c 2 c 1 c 3 c 2 c 2021-6-13清华工程流体力学课件第七章气体一维高速流动 二、正激波前后气流参数二、正激波前后气流参数 如图7-7所示，正激波前和正激波后各气流参数的下标分 别为1和2。由于圆管的截面积不变，所以连续性方程可写成 （a） 若忽略摩擦的影响，则动量方程可写成 或 （b） 气流通过激波时受到急剧地压缩，由于其时间极短，所产 生的热量来不及外传，故使气流的熵增加。所以气流通过激 波时的突跃压缩过程是一个不可逆的绝热过程。于是，气流 在激波前后的总能量相等，并保持不变，对于完全气体能量 方程可写成 （c） 或 （d） 式中临界声速 也保持不变。 2211 VV

42、)( 121121 VVVpp 2 222 2 111 VpVp 常数 02 2 2 1 2 1 22 TcTc V Tc V ppp 常数 21 1 11212 2 0 0 2 2 2 2 1 1 2 1 cppVpV c 2021-6-13清华工程流体力学课件第七章气体一维高速流动 图7-7 正激波 2021-6-13清华工程流体力学课件第七章气体一维高速流动 将气体状态方程应用与正激波前、后的状态，得 （e） 将式（b）的两边各除以式（a）的两边，得 （f） 由能量方程（d）可得 （g） （h） 将式（g）和（h）代入式（f），简化后得 由于 ，所以 （7-35） 这就是著名的普朗特公式

43、，再由动量方程（b）和连续性方 程（a）可知 )( 112212 TTRpp 2 12 2 2 22 1 1 1 V V V p V p 2 1 2 1 1 1 1 2 1 Vc p 2 2 2 2 2 1 1 2 1 Vc p 2 122112 cVVVVVV 12 VV 2 21 cVV 1 22 11 2 22 2 1112 1 V V VVVpp 2021-6-13清华工程流体力学课件第七章气体一维高速流动 由于激波是压缩波，即 ，因此 。所以由式（7- 35）可得重要结论：若正激波前是超音速流，则在正激波后 必定是亚音速气流。 由于 和 ，则式（b）可改写成 所以 （7-36） 又由

44、 代入式（c）得 （i） 再由于 ，所以式（i）可写成 （7-37） 2 p 1 p 2 V 1 V p c 2 c V Ma 2 11 2 22 2 11 2 2221 MapMapVVpp 2 2 2 1 1 2 1 1 Ma Ma p p RTMacMaV 22 2 211 2 1 2 1 2 1 TcRTMaTcRTMa pp Rcp 1 2 2 2 1 2 2 2 1 1 2 1 2 1 1 1 2 1 1 2 1 2 1 Ma Ma Mac Mac T T p p 2021-6-13清华工程流体力学课件第七章气体一维高速流动 由状态方程和式（a） 所以 （7-38） 现将式（7-3

45、6）和式（7-37）代入式（7-38）得 或 1 2 11 22 11 22 21 12 1 2 RT RT Map Map Vp Vp p p T T 2 11 22 1 2 Map Map T T 2 1 2 2 2 2 2 1 2 2 2 1 1 1 1 2 1 1 1 2 1 1 Ma Ma Ma Ma Ma Ma 4 2 2 1 2 1 22 1 2 1 11 2 1 MaMaMaMa 2 1 2 2 2 1 2 1 2 2 1 2 1 121MaMaMaMaMa 0 2 1 1 2 1 Ma 2021-6-13清华工程流体力学课件第七章气体一维高速流动 简化成 或 （7-39） 式

46、（7-39）最简单但无意义的解是 ，即上、下游的 马赫数相等，无正激波存在的情况。式（7-39）的另一个解 就是正激波前、后马赫数的关系式 （7-40） 0 2 1 11 2 1 2 1 2 1 2 2 4 1 4 2 2 1 MaMaMaMaMaMa 0 2 1 1 2 1 2 1 2 2 2 1 2 1 2 2 MaMaMaMaMa 21 MaMa 2 1 2 1 1 2 1 2 1 2 2 Ma Ma Ma 2021-6-13清华工程流体力学课件第七章气体一维高速流动 将式（7-40）代入式（7-36）和式（7-37），得 （7-41） （7-42） 再将式（7-41）和式（7-42）代

47、入式（7-38），得 （7-43） 式（7-40）至（7-43）表示正激波前、后各气流参数之比都 是波前马赫数的函数。所以，当波前各气流参数已知时，就 可以从这些公式求得波后各气流参数之值。 1 12 2 1 1 2 Ma p p 1 1 1 12 1 2 1 2 1 2 1 2 1 2 Ma Ma Ma T T 2 1 2 1 2 1 1 2 1 12 Ma Ma V V 2021-6-13清华工程流体力学课件第七章气体一维高速流动 第五节 应 用 举 例 【例例7-1】 空气从大容器经喉部直径为25mm的缩放喷管流 向大气，设大容器中的计示压强和温度各为690kPa和 40，大气压强 10

48、1.3kPa（abs.）。若不计空气流 过喷管时的损失，试求临界流速、出口流速、出口马赫数 和出口截面的直径以及 。空气 J/（kgK）， 1.4。 解析 【例例7-2】 在亚声速飞行的飞机上，装有皮托管和静压管， 用来测量飞机的飞行速度。今测得大气的静压 0.75105Pa，气温 -10，动压 0.25105Pa。 试计算飞机的飞行速度。 解析 pt 0 p m p2 287R pp 0 2021-6-13清华工程流体力学课件第七章气体一维高速流动 【例例7-3】 核爆炸产生的爆震波（激波）以16000m/s的速 度在静止的大气中传播，原来空气的压强为 1.0133105Pa，温度为300K

49、，试计算： (1)激波相对静止空气的马赫数 ； (2)激波后相对静止观察者的压强和温度的滞止值。 解析 2021-6-13清华工程流体力学课件第七章气体一维高速流动 【解解】 空气的临界压强由式（7-25） 0.528（690+101.3）=417.7（kPa） 而 101.3kPa，于是 ，故采用缩放喷管。 空气的气体常数 J/（kgK），则 （kg/m3） 临界流速由式（7-33） （m/s） 出口速度为 返回【例例7-1】 * p 2 p * 2 pp 287R 81. 8 )40273(287 10) 3 .101690( 3 0 0 0 RT p 324 81. 8 10) 3 .1

50、01690( 08. 108. 1 3 0 0 * p cV 1 0 2 0 0 2 1 1 2 p pp V 4 . 1 14 . 1 3 3 .101690 3 .101 1 81. 8 10) 3 .101690( 14 . 1 4 . 1 2 528（m/s） 2021-6-13清华工程流体力学课件第七章气体一维高速流动 空气的定压比热 J/（kgK），则 ，由式（7-11） 所以 K=-99 出口马赫数为 通过喷管的流量为 1005 p c Tch p 2 00 2 2 2 2 1 2 VTc V Tc pp )(2 20 2 0 2 2 TTcVV p 0 0 V)40273(10052528 2 2 T 174 2 T 2 17