流体力学第八章气体的一元流动

  • 格式:docx
  • 大小:476.74 KB
  • 文档页数:17

第8章 气体的一元流动一、 学习的目的和任务1.掌握可压缩气体的伯努利方程 2.理解声速和马赫数这两个概念3.掌握一元气体的流动特性,能分析流速、流通面积、压强和马赫数等参数的相互关系 4.掌握气体在两种不同的热力管道(等温过程和绝热过程)的流动特性。

二、 重点、难点128.1(8.1-1)12,u u ——流动气体两点的平均流速在气体动力学中,常以g ρ乘以上式(8.1-1)后气体伯努利方程的各项表示称压强的形式,即2212112222u u p gz p gz ρρρρ++=++(8.1-2)由于气体的密度一般都很小,在大多数情况下1gz ρ和2gz ρ很相近,故上式(8.1-2)就可以表示为22121222u u p p ρρ+=+(8.1-3)前面已经提到,气体压缩性很大,在流动速度较快时,气体各点压强和密度都有很大的变化,式(8.1-3)就不能适用了。

必须综合考虑热力学等知识,重新导出可压缩流体的伯努利方程,推导如下。

如图8-1所示,设一维稳定流动的气体,在上面任取一段微小长度ds ,两边气流断面1、2的断面面积、流速、压强、密度和温度分别为A 、u 、p 、ρ、T ;A dA +、u du +、p dp +、d ρρ+、T dT +。

取流段1-2作为自由体,在时间dt 内,这段自由体所作的功为()()()W pAudt p dp A dA u du dt =-+++(8.1-4)根据恒流源的连续性方程式,有uA C ρ=(常数),所以上式(8.1-4)可写成 由于在微元内,可认为ρ和d ρρ+很相近,则上式可化简为()p p dpdpW Cdt Cdt ρρ--==-(8.1-5)又对1-2自由体进行动能分析,其动能变化量为222111()22E m u du m u ∆=+- (8.1-6)同样地根据恒流源的连续性方程式uA C ρ=(常数),故有12m m uA C ρ===上式就可以写成1(2)2E Cdt udu Cudtdu ∆==(8.1-7)根据功能原理有WE =∆,化简得0dpudu ρ+=(8.1-8)图8-1ds 微元流段该式就是一元气体恒定流的运动微分方程对上式(8.1-8)进行积分,就得一元气体恒定流的能量方程22dpu C ρ+=⎰(8.1-9)式中C 为常数。

上式表明了气体的密度不是常数,而是压强(和温度)的函数,气体流动密度的变化和热力学过程有关,对上式的研究取要用到热力学的知识。

下面简要介绍工程中常见的等温流动和绝热流动的方程。

(1) 等温过程等温过程是保持温度不变的热力学过程。

因pRT ρ=,其中T =定值,则有pC ρ=(常数),代入式(8.1-9)并积分,得2ln 2pu p C ρ+= (8.1-10)(2) 绝热过程绝热过程是指与外界没有热交换的热力学过程。

可逆、绝热过程称为等熵过程。

绝热过程方程pC γρ=(常数),代入式(8.1-9)并积分,得 212pu C γγρ+=-(8.1-11)式中γ为绝热指数。

8.2声速和马赫数8.2.1声速微小扰动波在介质中的传播速度称为声速。

如弹拨琴弦,使弦振动了空气,其压强和密度都发生了微弱的变化,并以波的形式在介质中传播。

由于人耳能接收到的振动频率有限,声速并不限于人耳能接收的声音传播速度。

凡在介质中的扰动传播速度都称为声速。

如图8-2所示,截面面积为A 的活塞在充满静止空气的等径长管内运动,0u =时(0t =),管内压强为p ,空气密度为ρ,温度为T ;若以微小速度du 向右推进时间dt ,压缩空气后,压强、密度和温度分别变成了p dp +,d ρρ+和T dT +。

活塞从右移动了dudt ,活塞微小扰动产生的声速传播了cdt ,c 就为声速。

图8-2 微小扰动波的传播取上面的控制体,列连续性方程得()()cdtA d c du dtA ρρρ=+-(8.2-1)化简并略去高阶无穷小项,得du cd ρρ=(8.2-2)又由动量定理,得()[()]pA p dp A cA c du c ρ-+=--(8.2-3)同样化简并略去高阶无穷小项,得(8.2-4)(8.2-5)c 越(8.2-6)式中p C γρ=(8.2-7)上式两边对ρ求导,得11dp p pC d γγγγργργρρρ--=== (8.2-8)又由理想气体状态方程g pR T ρ=和上式(8.2-8)、式(8.2-5)联立,得c==(8.2-9) 综合上述分析,有(1)由式(8.2-5)得,密度对压强的变化率ddpρ反映了流体的压缩性,ddpρ越大,则dpdρ越小,声速c也越小;反则声速c越大。

由此可知,声速c反映了流体的可压缩性,即声速c越小,流体越容易压缩;声速c越大,流体也越不易压缩。

(2(8.2-10) (3Ma表示(8.2-10)下面讨论微小扰动波的传播规律,可分为四种情况:(1)如图8-3()a所示,0u=,扰动源静止。

扰动波将以声速向四周对称传播,波面为一同心球面,不限时间,扰动波布满整个空间。

(2)如图8-3()b所示,u c<,扰动源以亚声速向右移动。

扰动波以声速向外传播,由于扰动源移动速度小于声速,只要时间足够,扰动波也能布满整个空间。

(3) 如图8-3()c 所示,u c =,扰动源以声速向右移动。

由于扰动源移动速度等于声速,所以扰动波只能传播到扰动源的下游半平面。

(4) 如图8-3()d 所示,u c >,扰动源以超声速向右移动。

由于扰动源移动速度大于声速,扰动波的球形波面被整个地带向扰动源的下游,所以扰动波只能传播到扰动源的下游区域,其区域为一个以扰动源为顶点的圆锥面内。

称该圆锥为马赫锥。

锥的半顶角θ称为马赫角,从图中可以看出1sin c u Maθ==(8.2-11)知,8.3性。

(8.3-1)(8.3-2)度增大,导致压强增大、气体压缩。

马赫数为两者相对变化量的系数。

因此,当时,即超声速流动,密度的相对变化量大于速度的相对变化量;当1Ma <时,即亚声速流动,密度的相对变化量小于速度的相对变化量。

以下再分析流速与断面积的关系8.3.2气体流速与流道断面积的关系对一元气流得连续性方程uA C ρ=(常数)两边取对数,得 对上式微分,得或d du dA u Aρρ=-- (8.3-3)将式(8.2-13)代入上式,得2(1)dA du Ma A u=- (8.3-4)从上式我们可以看到,1Ma =是一个临界点。

下面讨论其在亚声速和超声速流动下的情况。

(1) 亚声速流动时,即1Ma <。

面积相对变化量和速度相对变化量反向发展,说明了气体在亚声速加速流动时,过流断面逐渐收缩;减速流动时,过流断面积逐渐扩大。

(28.4(8.4-1)(8.4-2)式中根据连续性方程,有111222u A u A uA ρρρ==,对于等径管道因12A A A ==,得11u u ρρ=(8.4-3)又由热力学等温过程方程pC ρ= 即1C p ρ-=和111C p ρ-=,有或11p p ρρ=和11p uu p= (8.4-4)将式(8.4-4)代入式(8.4-2)并改写为211102pdp du dsp u u Dλρ++= (8.4-5)如图8-3所示,设在等温管道中,取一微小流段ds ,在1-2段对上式(8.4-5)进行定积分,得上式积分得222121111u l p p u D λ⎛⎫- ⎪⎝⎭(8.4-6)若管道较长,且气流速度变化不大,则可以认为212lnu l u Dλ<<,略去对数项,上式可写成 2221111l p p p u D λρ⎛⎫=- ⎪⎝⎭(8.4-7)2211211()D u p p p lρλ=-(8.4-8)质量流量公式为25222111121()416m D D Q u p p p lρππρλ==- (8.4-9)上面各项就是计算等温管道压强、流速和流量的计算公式。

8.4.2流动特征分析前面已经给出了气体连续性方程uA C ρ=,其中A 不变,则有u C ρ'=,对该式取对数并积分,得0d duuρρ+= (8.4-10)由热力学方程pRT C ρ==,积分得dp d p ρρ= (8.4-11)联立上面两式(8.4-10)和(8.4-11),以及声速公式pc γρ=uMa c=并整理。

得 22(1)2du Ma ds u Ma Dγλγ=-(8.4-11)图8-3 微元管流从上式我们可以看出,如果Ma >,,即210Ma γ-<,0ds >,则0du <;又对于大多数气体的指数常数1γ>,且实际工程等温管道中气流的速度不可能无限增大,21Ma γ-不可能等于或小于0,所以只有Ma <Ma >时,只能按Ma =(极限值)计算,该极限值计算的管长又称为最大管长,即实际管长超过最大管长时,进口断面的流速将受到阻滞,必须减小管长。

8.5 气体在绝热管道中的流动在实际的气体输送管道中,常常在管道外面包有良好的隔热材料,管内气流与外界不发生热交换,这样(8.5-1)(8.5-2)(8.5-3)(8.5-4)可得111221211111ln 2u l p p p u u D γγγγγγλργ++⎛⎫+-=+ ⎪⎝⎭(8.5-5)考虑到管道较长,流速变化也不大,21ln2u s u Dλ<<,略去对数项,可写成11122111112lp p p u Dγγγγγγλργ+++=-(8.5-6)1u =(8.5-7)质量流量为(8.5-8)(8.5-9)1<时,8.6也为零,驻点处的流动参数也可认为是滞止参数。

滞止参数常用下标“0”标识,如000,,p T ρ分别表示滞止压强、滞止密度、滞止温度。

由绝热过程方程式(8.1-11),按滞止参数的定义,可得滞止参数和某一断面的运动参数间的关系为20112p pu γγγργρ=+--(8.6-1)又由完全气体状态方程pRT ρ=得,上式可写为20112u RT RT γγγγ=+-- (8.6-2)即20112T u T RT γγ-=+ (8.6-3)又声速c =上式改写成马赫数的形式为(8.6-4)程pρ(8.6-5)(8.6-6)(8.6-2)得(8.6-7)式中例题:8.6.2临界状态参数气体从当地状态等熵地改变速度达到声速时(即1Ma =),所具有的状态称为与该当地状态对应的临界状态,相应的状态参数称为临界参数,与滞止状态一样,临界状态可以是流动中实际存在的,也可以是假想的状态。

临界状态参数常用下标“*”表示。

如*T 、p *分别称为临界温度、临界压强等。

在等熵流中所有的临界参数都是常数,因此可作为参考状态参数。