第十章 膨胀波和激波
折角与波后气流的马赫数的关系。
称为Prandtl-Mayer函数,或Prandtl-Mayer角。它表示对于来流马赫数为1时,经过膨
胀后气流速度达到M时,所能偏转的角度。这个式子已经制成表格,在知道M数或 后,从这个表格中可以查到另一个值。
(三)对于来流马赫数为M1而最后速度为M2,气流总偏转角为:
㈣如果壁面转折是朝上的,膨胀波将沿逆时针方向,此时普郎特-梅耶角取正数。否则为负数。(也可以理解为取绝对值。)
㈤当最终马赫数为∞,达到普郎特-梅耶角的最大可能值
但这只是一个理论值,因为早在达到这个速度前,气流就会冷凝了。
[例1]马赫数1.4的空气,绕一外钝角偏转了20o。已知来流的初始静压和静温分别是p1=101325N/m2,T1=288K,求 膨胀后气流的马赫数、静压和静温。
[解]由来流马赫数M1,可以查表或者根据Prandtl-Mayer角求得 v1=8.987 ,这个角度表示音速的来流经过膨胀后气流马赫数为M1后的偏转角。
这样,从音速的来流膨胀到M2的总偏转角为v2=8.987+20=28.987 再查表或计算得到M2=2.096。
因为气流经过膨胀波是绝能等熵的,所以总温总压不变,借此可以计算出波后静压和静温。
㈥在连续转折或凸曲面处的膨胀波。 不论多道转折,还是曲面转折,在已知来流马赫数后,只要知道气流膨胀之后的马赫数、或者总的折转角,便可求得另一个。
[例2]拉瓦尔喷管的出口处,Me=1.2 ,气流出口处的总压为3个大气压,问: ①出口处气流是膨胀还是压缩?
②膨胀或者压缩的气流偏转角多大?绝热指数取1.4。
[解] ①根据拉瓦尔喷管的出口总压和马赫数,可求得出口气流的静压为: p=1.237
因为这个压力比环境压力高,所以气流必须继续膨胀减压
②为了求得气流膨胀后的转角,必须求得气流膨胀后的马赫数M2。根据膨胀波是绝能等熵的过程,而膨胀后气流压力必须达到环境压力,因此可以求出气流膨胀后的马赫数为:
[说明]拉瓦尔喷管的作用是依靠气流加速产生推力的,当出口处产生膨胀时,表明气流在喷管内膨胀不足,因此这个工作状态下,喷管是损失掉能量了。常说这是喷管被截短了。
㈦气流绕凹曲面等熵压缩是膨胀波的逆过程,也可以用普郎特-梅耶膨胀波理论解决。
参见课本例题2。
需要特别指出,等熵压缩也称为微弱压缩,而且等熵压缩波最后会汇总在一起,形成强烈的压缩,那里 压缩不再是等熵的。
2、激波(shock wave)理论 1.激波的特征
当超音速气流流过大的障碍物时,气流在障碍物前受到急剧的压缩,压力和密度突然显著增加。所产生的压力扰动波以比音速大得多
的速度传播,波面所至之处气流参数发生突然的变化,这种强压力扰动波称为激波或冲波。
气流通过激波时,速度下降,而压力、密度和温度增加。总温不变但总压降低。因此这是一个熵增过程,换言之是一个不可逆过程。
激波是由许多微弱扰动波迭加而成的,有一定强度的压缩波。
2、激波的分类 (1)正激波
激波面与气流来流方向垂直,气流经正激波后不改变来流方向。 (2)斜激波
激波面与气流来流方向不垂直,气流经斜激波后要改变流动方向。 (3)曲激波(或称弓形激波、脱体激波)
由正激波(在中间部分)和斜激波系组成。
3.激波的厚度
在无粘性不导热的理想气体中,激波成为无厚度的数学上的间断面。
在实际气体中,由于粘性和热传导,激波是一个极薄的过渡区,厚度与气体分子的平均自由行程(≈10-5mm)同一数量级。各气体参数是在这个极小厚度内连续变化,所以也可以把激波看作是一个不连续的间断面。 4.波阻
超音速气流经过激波后,部分动能不可逆地转变为热能而损失掉,因而产生一种超音速气流所特有的阻力损失,称为波阻。气流通过正激波时波阻最大。
5、正激波参数计算 ㈠基本方程
取相对坐标系 ,把直管内的激波当作静止的,这样是超音速来流穿过激波。波前参数是p1、T1、ρ1、V1;波后参数为p2、T2、ρ2、V2。 取包围激波的一个
控制体,可以使用积分形 式方程建立波前波后气流 的基本方程。
思考:为什么不能用 微分方程?
连续方程: 动量定理:
能量方程:
状态方程:
借助于上面这个一维流的气体动力学方程组,可以得到关于激波参数计算的四组公式。 它们分别是:Prandtl激波关系式;波前波后气流参数关系;Ranking-Hugioniut关系;总压损失关系式(或称熵增关系式)。
㈡正激波基本关系式之一——普郎特激波公式
这个关系表明,正激波后气流总是亚音速的。
当 M1->00 ,那么正激波后最小的马赫数M2=0.378 。(γ=1.4) 【证明】用动量方程除以连续方程,得到
由于能量方程可写成
㈢正激波的激波关系之二——波前波后气流参数比与波前马赫数的关系 可见,波前波后参数之比只取决于波前M数
从上面的关系式可以分析出:
①如果波前马赫数无限大,则压强比也是无限大。这意味着,来流越强,激波也越强。可以通过激波实现任意大压强的压缩。
②因为波前马赫数无限大时,密度之比有一个极限: (r+1)/(r-1) 。所以通过激波不可能实现无限高密度的压缩。
③
还可以得到总压比的关系式(请自己思考)。
㈣激波关系之三——Rankine-Hugoniut关系
这组关系描述了波前波后气流参数比之间的相互联系。
Rankine-Hugoniut关系式的重要性在于,它对于任何激波都是正确的。
㈤激波关系之四——激波的熵增关系
因为经过激波后,气流总压变化,而总温不变。所以
由于波后,气流总压下降,所以气流的熵必然增加。
6、斜激波
㈠斜激波与正激波的关系
超音速气流经过一个楔角为 S 的楔形时,气流被迫转向了 S ,如果超音速气流足够强,则在楔形的顶部产生一道斜激波,把斜激波与来流方向的夹角 B 称为 激波角,而 B 称为气流 偏转角(气流角)。
按动量方程,可以证明气流沿切向的分速在波前波后是不变的(自己分析)。
而法向速度穿过激波时是变化的。从气体动力学基本的积分方程可知,气流参数的变化可以如下考虑:即在正激波的关系式中,如果遇到速度(包括马赫数或速度系数),则用斜激波的法向分速代替之,就得到穿过斜激波后气流参数的变化。 因为波前波后气流的法向速度分别为:
㈡斜激波的激波关系式: 1、波前波后气流速度
这表明,斜激波后气流的法向速度必定是亚音速的,因为波前气流法向速度必定是超音速的。但是波后气流的总速度可以是亚音速或者可以是超音速的,因为还取决于切向分速。
至于为什么斜激波前气流法向速度必定是超音速的,这可以从下面的压强比关系可以看出
2、斜激波前后气流参数比与来流马赫数的关系
即 V1n>C1 ,因此斜激波前法向速度必定是超音速的。
㈢气流偏转角与激波角的关系
对于斜激波的计算,最重要的就是确定这两个角度。通常知道一个而需要求另一个。 因为通过激波面的流量与沿波面的分速Vτ无关,故连续性方程可写为:ρ1V1n=ρ2V2n
这就是斜激波的激波角 与气流偏转角 、及波前马赫数M1之间的关系。
【说明】
1、这个公式计算起来是困难的,但是可以通过查斜激波表或者斜激波图得到。
2、最重要的在于,这个 表中一个M1和一个 对应了两个 。这两 个解分别称为强解和弱 解,分别对应了大的和 小的 。
到底取哪个解,必
须根据具体问题而定。
实验指出,这两道激波中较弱的一道最常发生,特别是尖头体前面的附体激波。但是如果下游压力极高,就有可能产生强激波解。
如果不能肯定到底是哪道激波,那么可以先按弱激波计算,再看看这个计算的结果跟下游的参数是否符合。
【例2】空气流过一个无摩擦的
表面并偏转 S=15 ,已知M1=3, p1=101325N/m2,T1=300K。假定产生一道弱激波,试求:波后气流的马赫数、压强、密度和速度。 【解】根据来流M数及气流偏转角,可以查得激波角为:B=32.2404 利用激波角和来流马赫数,可以得到波前波后气流参数比,比如
㈣脱体激波
对于一个给定的来流马赫数M1,气流的偏转角δ有一个极限值δmax。超过了这个极限值后,激波不再附体而离开物体形成脱体激波。因为气流角对应了激波 角,所以对应有一个
极限的激波角βmax
对给定的M1,按此式可求出最大的激波角,再按激波角与气流角关系求出最大的气流偏转角。
【说明】
1、一旦激波脱体,那么便不能再使用斜激波关系。
2、激波脱体后,物体正前方一段的激波可以当作正激波,因此,这个区域后方的气流是亚音速。而激波的外伸段是斜激波,因此,其后面的气流可能是超音速也可能是亚音速。 3、即使来流是理想无旋的,但是脱体激波后面的气流也是有旋的,这是由著名的Crocoo定理表明的。