1、第第1010章章 一维气体动力学基础一维气体动力学基础 在前面的章节中，都将流体视为不可压缩流体，即流体的密度 =常量。在工程实际问题中，当气体的流速很高，压差很大时；温度变化很大且伴随热效应时，气体的密度会发生显著的变化，此时气体的运动规律和不可压缩流体大相径庭。研究此类问题时必须采用可压缩流体模型。本章主要讨论完全气体（在热力学中称理想气体）一维恒定流动。在学习本章的过程中不仅需要流体力学知识，还需要一定的热力学知识，在进行气体动力学计算时，压强只能用绝对压强，温度只能用开尔文温度。工程流体力学工程流体力学10.1 10.1 声速和马赫数声速和马赫数 10.1.1 10.1.1 声速声速
2、凡是微小扰动在流体介质中的传播速度都定义为声速，它是气体动力学的重要参数。dvdvv=0p(a)p+pd+dcc-vdp(b)p+pd+dcy图 0.1 1声速传播过程 对于小扰动波的传播过程，可通过下例说明。取等断面积为A，左端带活塞的直长管如图10.1（a）。工程流体力学工程流体力学 管中充满静止的可压缩气体，压强为p，密度为 。活塞在力的作用下，以微小速度dv向右移动，产生的一个微小扰动平面波不断地从左端波及到右端，波的传播速度即声速，以符号c表示。特别要注意的是声速c与气体受扰动后的速度dv是不同的。为了便于分析波阵面前后流体状态参数的变化关系，将坐标系固定在波阵面上图10.1（b），
3、这样，对位于该坐标系的观察者而言，流体的流动是恒定的 声速公式 ddpc拉普拉斯在1816年提出声音的传递是一个等熵过程。工程流体力学工程流体力学 在声波传递过程中，热力学参数的变化是无穷小量，忽略了黏性作用，因而整个过程可视为可逆的绝热过程（即等熵 s=常数），那末式（10.1）更精确的表达式为 ddspc 公式中下标s代表等熵过程。上式不仅适用于微小扰动平面波，也适用于球面波，对气体、液体均适用。对于完全气体等熵流体的状态参数方程式为 p常数工程流体力学工程流体力学式中 称为比热比（或称为绝热指数），对于空气 1.4于是可导出完全气体的理论声速公式 cRT式中 称为气体常数，空气的 。R2
4、87J(kg K)R 由以上声速公式可得出：1.是反映流体的压缩性，当越大，表示流体越易压缩，此时由式（10.1）越小；反之，当流体越不易压缩，则声速c越大，若流体为不可压缩流体，那么声速 。因而声速是反映流体压缩性大小的物理参数。ddpddpddpcc 工程流体力学工程流体力学2.由式（10.4）可得，不同的气体有不同的比热比，及不同的气体常数R，因而不同的气体声速是不同的。如在常压下，15空气中1.4 287(273 15)340m/scRT在相同的压强和温度下，氢气的声速为 。1 295m/sc 3.声速与气体热力学温度T有关，如在常压下空气中声速为20.1cT由于在气体动力学中，温度是
5、空间坐标的函数，所以声速也是空间坐标的函数，为此，常称为当地声速。工程流体力学工程流体力学4.对于液体，由式（1.10）液体的弹性模量E和压缩系数 k的关系为1ddpEk代入式（10.1）得声速公式的另一种形式Ec10.1.2 马赫数和马赫锥马赫数和马赫锥 1 1马赫数马赫数 马赫数是惯性力与由压缩引起的弹性力之比，它是气体动力学中最重要的相似准数，即定义马赫数 vMac式中 v当地气流速度；c当地声速。工程流体力学工程流体力学气体动力学中，依据马赫数对可压缩气流进行分类：，即 ，称为超声速流动；，即 ，称为声速流；，即 ，称为亚声速流动。1Mavc1Mavc1Mavc对于气体流动以 为界，对
6、于 ，为不可压缩流动，对于 ，为可压缩流动。0.3Ma 0.3Ma 0.3Ma【例10.1】用声纳探测仪，探测水下物体，已知水温20，水的弹性模量 ，密度为 ，今测得往返时间为6秒，求声源到该物体的距离。Pa1088.19E3kg/m2.998【解】由式（10.6）91.88 101 372.4m/s988.2Ec从声源到物体之间的距离为1 372.4 34 117mct 工程流体力学工程流体力学【例10.2】某飞机在海平面和11 000m高空均以速度 飞行，问这架飞机在这两个高度飞行时的马赫数相同吗？319.4m/s【解】由于海平面的声速为 。340m/sc 故海平面的飞行飞机，为亚声速飞行
7、。319.40.94340vMac 在11 000m高空飞行时，该处的温度为216.5K（见第二章），则由式（10.5）20.1295.8m/scT故该高度飞行的飞机319.41.08295.8vMac为超声速飞行。工程流体力学工程流体力学工程流体力学工程流体力学2.2.马赫锥马赫锥 图10.2是一小扰动波（例如点声源）在四种流动中的转播。c2cc3c4cO2cv=0c3c4cO2cv2v3v4vv cc3c4cO2cABv=c(a)(b)(c)3c4cOABv cv2v3v4v(d)图 0.2 1小扰动波在不同来速流场内的传播工程流体力学工程流体力学（1）当小扰动波在静止流场中传播（）图10
8、.2（a）。0,0vvMac（2）当小扰动波在亚声速流场中传播（）如图10.2（b）。0,01vcMa（3）当小扰动波在声速流场中传播（），此种情况同上面（2）相同。10.2图(c)AOB平面是所有小扰动波的包络面，称为马赫波。它是寂静区和扰动区的分界面。,1vc Ma（4）当小扰动波在超声速流场中传播（），此时的马赫波不再保持为平面，而是以固定点O为顶点向右扩张的旋转圆锥面，这个圆锥面称为马赫锥，圆锥顶角的一半称为马赫角 如10.2图（d）。,1vc Ma其中 1arcsinarcsincvMa对于完全气体的马赫数可表示为 vMaRT 由于温度是气体分子运动动能的度量，所以式（10.8）说明
9、马赫数是流体宏观运动动能和分子运动动能之比。【例10.3】飞机距地面的1 000m上空，飞过人所在位置600m时，才听到飞机的声音，当地气温为15，试求飞机的速度、马赫数及飞机的声音传到人耳所需的时间。v1000m600m图 0.3 1马赫锥工程流体力学工程流体力学解：当地声速为 1.4 287(273 15)340m scRT马赫角 为(如图10.3)o1 000arctan59600由式（10.7）o1arcsin59Ma故马赫数 1.167Ma 飞机速度 340 1.167397m svcMa所需要时间 6001.51s397t 工程流体力学工程流体力学10.2 10.2 气体一维恒定流
10、动的基本方程气体一维恒定流动的基本方程 基本方程主要是由连续方程、欧拉运动微分方程和能量方程等组成。10.2.1 连续性方程连续性方程 如图10.4为一维恒定气流。V1V2A1A2图 0.4 1一维气体流动或者对任一过流断面满足V AC式（10.9）即为一维恒定气流的连续性方程，它的微分形式为ddd0VAVA工程流体力学工程流体力学10.2.2 欧拉运动微分方程欧拉运动微分方程 在一维恒定气流中，取长度为ds微段，并沿轴线方向为s轴。vpp+pdssds图 0.5 1一维气流微段如图10.5，即可得到 2dd02pv上两式称为完全气体一元恒定流动的欧拉运动微分方程式。式中的密度 不再是常数。而
11、 三者之间的关系由微分方程式来确定。为求解此方程式除了要应用气流的连续性方程外，还必须补充气体状态方程，或者热力学过程方程。pv,工程流体力学工程流体力学 连续性方程和欧拉运动微分方程也可引用马赫数Ma来表示如下 连续性方程可表达为 2dd1AvMaAv欧拉运动微分方程为2ddvMav 10.2.3 不同形式的能量方程不同形式的能量方程（1）气体一维定容流动 Cgvp22工程流体力学工程流体力学 上式为不可压缩流体，不计质量力的能量方程。表示一维气流各断面上单位质量（或重量）具有的压能和动能之和守恒。（2）气体一维等温流动 在等温流动中，T=常数，则气体状态方程 CRTp2ln2pvp 常数或
12、者2ln2vRTp 常数（3）气体一维等熵流动 工程流体力学工程流体力学 在热力学中，无能量损失且与外界又无热量交换的情况下，为可逆的绝热过程，又称等熵过程。212pv常量212vRT常量2112ppv常量 ，在热力学中，这项正是在等熵过程中，单位质量气体所具有的内能e。表示为 11p工程流体力学工程流体力学22pveC 该式表明完全气体的等熵流中，沿流任意断面上，单位质量气体所具有的内能、压能和动能之和是不变的。【例10.4】用文丘里流量计来测量空气流量（图10.6），流量计进口直径 ，喉管直径 ，实测进口断面处压强 （相对压强），温度为20，喉管处压强 （相对压强），试求空气的质流量。（设
13、当地大气压 ）150mmd 220mmd 135kPap 215kPap a101.3kPap 工程流体力学工程流体力学d11122图 0.6 1文丘里流量计d2【解】气流通过文丘里流量计时，由于流速大，流程短，气流和壁面接触时间短，来不及进行热交换，且摩擦损失亦可不计，因此按一维恒定等熵流动来处理。先计算进口断面1处，喉管断面2处空气的密度由 式，进口断面空气的密度RTp33111(35 101.3)101.62kg/m287(27320)pRT由式10.3），即 1212pp工程流体力学工程流体力学111.432211101.3 151.621.446kg/m101.335pp由连续性方程
14、式（10.9），得 1 11222V AV A21121112221.620.05471.4460.024AVVVVA将以上量代入等熵能量方程式（10.16）221122121212pVpV工程流体力学工程流体力学22331171.4136.3 101.4116.3 101.4 11.6221.4 11.4462VV解得 123.25m/sV 故空气的质流量21 111.62 23.250.050.074kg/s4mQV A【例10.5】氦气（）作等熵流动，在管道截面1处参数为 ，测得截面2处的速度为 ，求该截面上的以 及 值。1.67,2077J kg KR161T 1,65m sV 218
15、0m sV 2T21pp解：由等熵流动能量方程式（10.16）工程流体力学工程流体力学212pv常数从截面122212121212VVRTRT 因此 2221.67651.671802 077273612 0771.67 121.67 12T解得 或 2331.28KT 58.28由等熵过程 1.6711.67 12211331.280.979 827361pTpT工程流体力学工程流体力学10.2.4 10.2.4 一维等熵流动气体动力学函数一维等熵流动气体动力学函数 1 1用滞止状态参数表示的气体动力学函数用滞止状态参数表示的气体动力学函数 当流体质点由某一个真实状态经等熵过程速度降为零。（
16、可以假想）这时流体质点的状态称为对应于真实状态的滞止状态。流体质点所具有的流体参数称为该真实状态的滞止参数，以下标“0”表示。例如以 分别表示滞止压强，滞止密度，滞止温度，以及滞止声速。在工程中，如气体从大体积的容器中流出（如煤气储气罐等），容器内气体的流速可视为零，那其它参数就是滞止参数；当气流绕过某物体，则驻点处气流的流动参数也是滞止参数。0,000,pT c工程流体力学工程流体力学 12012011201220122001110 19()21110 1921110 1921110 1921110 192TMaaTpMabpMaccMadcvv cMaMaecc c上述公式称为用滞止参数表
17、示的等熵流动气体动力学函数。工程流体力学工程流体力学图10.7表示气体动力学函数的曲线。为便于计算，气体动力学函数列表10.2如下。1.00.80.60.40.2001230.20.40.60.81.01.21.41.61.82.0Mvpp00TT0cc0vc0cc0TT00pp0图 0.7 1气体动力学函数曲线ac0,工程流体力学工程流体力学表表10.2 气体动力学函数数值关系（一维等熵关系）气体动力学函数数值关系（一维等熵关系）0.0 1.000 0 1.000 0 1.000 0 1.000 0 0.000 0 0.000 0 0.000 00.1 0.993 0 0.995 0 0.9
18、98 0 0.999 0 0.099 9 0.171 8 0.173 00.2 0.972 5 0.980 3 0.992 1 0.996 0 0.199 2 0.337 4 0.346 90.3 0.939 5 0.956 4 0.982 3 0.991 1 0.297 3 0.491 4 0.523 00.4 0.895 6 0.924 3 0.969 0 0.984 4 0.393 7 0.628 8 0.702 20.5 0.843 0 0.885 2 0.952 4 0.975 9 0.487 9 0.746 4 0.885 30.6 0.784 0 0.840 5 0.932 8
19、0.965 8 0.579 5 0.841 6 1.073 50.7 0.720 9 0.791 6 0.910 7 0.954 3 0.668 0 0.913 8 1.267 50.8 0.656 0 0.740 0 0.886 5 0.941 6 0.753 2 0.963 2 1.468 20.9 0.591 3 0.687 0 0.860 6 0.927 7 0.834 9 0.991 2 1.676 41.0 0.528 3 0.633 9 0.833 3 0.912 9 0.912 9 1.000 0 1.892 91.1 0.468 4 0.581 7 0.805 2 0.893
20、 1 0.987 0 0.992 1 2.118 41.5 0.272 4 0.395 0 0.689 7 0.830 5 1.245 7 0.850 2 3.121 22.0 0.127 8 0.230 0 0.555 6 0.745 4 1.490 7 0.592 0 4.636 72.5 0.058 5 0.131 7 0.444 4 0.666 7 1.666 7 0.379 3 6.480 03.0 0.027 2 0.076 2 0.357 1 0.597 6 1.792 8 0.236 2 8.674 55.0 0.001 9 0.011 3 0.166 7 0.408 2 2.
21、041 2 0.040 0 21.164 0Ma01p pfMa02fMa 03T TfMa04c cfMa0v cAAq1qyfMa工程流体力学工程流体力学2.临界状态和临界参数临界状态和临界参数 当气流速度 等熵地加速或减速到当地声速 的状态，也就是流动的马赫数等于1（可以假想）。此状态称为对应于真实状态的临界状态。临界状态的流动参量称为临界参数。以下标“”表示。例如记作 等。vc*,*,pT c12*21*121*122*211211211211M aTTM appM aM acc工程流体力学工程流体力学以 代入上式，可得到某真实状态所对应的滞止状态参数和临界状态参数之间的关系式:0Ma
22、*01*011*012*021212121TTppcc工程流体力学工程流体力学对空气（），具体数值为:1.4*0*0*0*00.8330.5280.6340.913TTppcc工程流体力学工程流体力学 在等熵的条件下，当温度降到绝对零度时，此时速度达到最大 的状态称为最大速度状态。由于真实温度不可能达到绝对零度，因此最大速度状态只具有理论意义，反映气流总能量的大小。maxv最大速度和滞止参数关系为:max002211vRTc最大速度与临界速度（声速）的关系max*1111vcv对于空气（），则 1.4max*2.45vv工程流体力学工程流体力学3.最大速度状态最大速度状态【例10.6】大容积压
23、缩空气罐中的压缩空气，经一收缩喷管向大气喷出，设喷嘴出口处的大气绝对压强为101.3kPa，温度为5，流速为 ，试求压缩空气罐中的压强和温度。234m/s【解】本流动可看作等熵流动 方法一：压缩空气罐中的空气速度可视为零，其流动参数为滞止参数。1.4 287(2735)334.2m/scRT2340.7334.2vMac由式10.19（b）喷口出口处声速马赫数120112ppMa1.41.4 121.4 1101.310.7140.5kPa2由式10.19（a）22011.4 1127810.7305.2K32.222TTMaC方法二：由式（10.16），压缩空气罐中的温度为0T2202343
24、.5 287(2735)122305.2K3.5 2871vRTTR由完全气体状态方程33111101.3 101.27kg/m287278pRT工程流体力学工程流体力学由等熵过方程式 ，pC111.4000111.27101.3ppp上式代入（10.16）式223427.1103.1015.33.10127.1105.3234.11030pp解得0140.5kPap 本题从以上两种解法可看出，应用滞止参数表示的等熵流动的气体动力学函数方法要简单得多。得：工程流体力学工程流体力学【例10.7】空气气流在收缩喷管进口截面上的参数为 ，在出口截面上 ，试求出口处的压强，温度和直径。5113 10
25、Pa,340KpT 1150m sV 146mmd 1Ma 解：本题也视为等熵流动 进口处 1111500.405 81.4 287 340VMaRT 现求滞止温度，由式 10.19（a）121101.4 112TMaT得 0351.2KT 出口处，,气流达到临界状态。1Ma 工程流体力学工程流体力学出口处的温度，由式10.23（a）*00.833292.5KTT由式10.19（a）和式10.19（b）3.50011pTpT503.360 4 10 Pap 得由式10.23（b）5*00.5281.774 10 Papp53*1.774 102.113kg m287292.5pRT*1.4 2
26、87 292.5342.82m sVcRT工程流体力学工程流体力学531113 103.074kg m287 340pRT由连续方程式10.9*1 11V AV A解得 1 1*1*3.074 1504636.7mm2.113 342.82VddV工程流体力学工程流体力学10.2.5 10.2.5 气流按不可压缩流体处理的限度气流按不可压缩流体处理的限度 对于低速气流，可忽略气体容易压缩的个性，而按照不可压缩流体处理。那么“低速”的限度就是下面要讨论的内容。完全气体一维流动，按不可压缩流体时，能量方程为:220vpp工程流体力学工程流体力学或者写成 1220vpp412220Mavpp 比较无
27、量纲式（a）和（c），气流按不可压缩流体处理时，能量方程式的计算相对误差为42Ma 在常温（15）下，空气的声速 ，倘若允许的相对误差为1，那么，相应的 ，此时相应的气流速度v为 340m/sc 2.0Ma工程流体力学工程流体力学340 0.268m/svcMa即，当气流速度小于 时，按不可压缩流体来处理时，其相对误差 。实质上气流按不可压缩流体来处理的限度是由计算要求的精度来决定。68m/s%1由式10.19（c）密度比式1120112Ma当 时，空气 ,代入上式,得2.0Ma1.402.12.0214.1114.1120密度的相对变化为%2102.10工程流体力学工程流体力学 在同样的气流
28、速度下，按不可压缩流体处理的话，其密度的相对变化较大。若要求气流密度的变化不超过1，相当的马赫数为 ，相应的气流速度v为 。141.0Ma48m/s【例10.8】某空气动力计算中，允许压强的相对误差 ，对于常温下的空气速度小于多少时可按不可压缩流体来处理；此时密度的相对变化为多少？%5.122vp【解】按式（10.27）压强的相对误差 42Ma工程流体力学工程流体力学据题意%5.142Ma故 245.0Ma常温下声速 340m/sc 气流速度 340 0.24583.3m/svcMa即，空气速度小于时 ，可按不可压缩流体来处理。83.3m/sv 030.1245.0214.1114.1120故
29、密度的相对变化为%3103.10由式10.19（c）工程流体力学工程流体力学10.3 10.3 喷管的等熵出流喷管的等熵出流 喷管是指在很短的流程内，通过改变断面的几何尺寸来控制气流速度的装置。由于高速气流在喷管内流动时来不及和外界进行热交换，同时摩擦阻力也可忽略不计，这样的流动过程可作为等熵流动。22dd1MaAMaA22dd1pMaApMaA22dd(1)1TMaATMaA工程流体力学工程流体力学 利用上述关系式，将断面面积A、气流速度v、压强p、密度 及单位面积的质流量 等与马赫数Ma之间的关系，能很清楚地表示如表10.3所列。v表表10-3 一维气体各流动参数随马赫数的变化关系一维气体
30、各流动参数随马赫数的变化关系1Ma 1Ma v,pT工程流体力学工程流体力学几点结论：1.亚声速气流（）在收缩管（）中，将加速（）和减压（）；在扩张管（）中，将减速（）和增压（）和不可压缩流体相似。2.超声速气流（）在收缩管（）中，将减速（）和增压（）；在扩张管（）中，将加速（）和减压（），与亚声速流恰好相反。1Mad0Ad0v d0p d0Ad0v d0p 1Mad0Ad0v d0p d0Ad0v d0p 图10.8为一收扩管，流体自左向右流动。若在该管道中达到声速，必定在最小截面处即喉部，称喉部的截面为临界截面记作。图 0.8 1收扩管 在扩张段流体被加速成超声速，并不断加速。这种流动成为
31、喷管流。瑞典工程师拉伐尔（Laval）将先收缩后扩大的喷管拉伐尔喷管（图10.9），用于蒸气涡轮机中。拉伐尔喷管在冲压式喷气发动机、超声速风洞等领域广泛地被应用。p0pb图 0.9 1拉伐尔喷管工程流体力学工程流体力学10.4 10.4 可压缩气体管道流动可压缩气体管道流动 对于可压缩气体的管道流动，有时要考虑摩擦阻力和热交换对压缩性的影响，需针对不同的热力过程进行分析计算。工程流体力学工程流体力学10.4.1 10.4.1 一维恒定等熵管流一维恒定等熵管流 下面主要对变截面管道内的流动进行分析。1 1管截面积和流动马赫数的关系管截面积和流动马赫数的关系122122121212121MaAMa
32、AMaMa 若已知 和 及某一处的马赫数，如 ，要求另一 ，却并不是很容易，为此假定一个参考截面 ，当流动至该截面马赫数 ，该截面流动参量就是临界参量，这个截面可实际存在于管流中，也可以是假想的。1A2A1Ma2MaA1Ma 工程流体力学工程流体力学121*2121AMaAMa 以上公式在一维管流的计算中被大量使用，公式算得的数据已列于表10.2之中。图10.10是根据上式绘制的 曲线。*AMaA1234123450AAM图 0.10 1 系曲线关AAM工程流体力学工程流体力学*e0.963 2AAAA2*0.963 2 0.0030.002 89mA 为假想的临界截面，即假想流体沿继续延伸的
33、喷管流动，在截面积 处达到声速，喷管其他截面上的参数与该假想临界截面上的参数关系符合等熵流气动函数关系。现*A*A*0.002 890.5780.005xAA由表10.2上，按 插入，查得 。*0.578AA 0.36Ma 工程流体力学工程流体力学【例10.10】设一喷管内为等熵流,出口截面积 ,出口 ，求喷管内截面积为 处的 。e0.8Ma 20.005mxA 2e0.003mA Ma解：由于 ，说明这是一个收缩喷管。由 ，查等熵流气动函数表可得 exAAe0.8Ma 2 2质流量的计算公式质流量的计算公式 质流量 121200112mQRT MaMaA121200112pMaMaART12
34、10,max*012mprQART 或者 1121,max*0021mQAp工程流体力学工程流体力学 当管道内存在临界截面 时，那末该处质流量达到最大值*1A Ma【例10.11】一个容积很大密闭容器中装满氮气，氮气的 ，容器中 ，氮气通过一收缩管向外流出，设出口处直径为 ，背压为 ，求流出氮气的质流量。解：工程上常称管外的环境压强为背压（或反压）用 表示，当封闭容器中压强 时管内无流动。当 时，在压差作用下产生流动。本题首先要判断在流动中管内是否会出现临界状态。bp0bpp0bpp工程流体力学工程流体力学1.4,297J/(kg K)R5004 10 Pa,298KpT5b10 Pap 50
35、mmd 由式10.22（a）*021TT故 *2298248.33K1.4 1T 由式10.22（b）55*00.5284 100.5282.112 10 Papp 由于 ，则说明在管道出口处前已出现临界状态，流量为最大，以后管内流动不再变化，通常称这种现象为壅塞现象。出口处的压强 ，不再等于 ，气流流出后经稀疏过程才降到 。b*pp*ppebpbp工程流体力学工程流体力学1121,max*0021mQAp其中 530004 104.52kg m297 298pRT故11.4 121.4 12520.051.4 4 1045241.4 1mQ 1.807kg s方法二:上面已分析由于 ，则收缩
36、管出口处b*pp*epp工程流体力学工程流体力学方法一:按式（10.37）质流量按 计算 *Qv A3*00.6340.634 4.522.866kg m*1.4 297 248.33321.33m svcRT2*2.866 321.330.051.807kg s4mQv A故由式10.22（c）10.4.2 绝热摩擦管流绝热摩擦管流 实际管流一般有两种，一种是在隔热的长管中流动，即具有摩擦但不考虑热交换的流动，如果这种流动在等截面管中流动被称为范诺（Fanno）流动。一种由于管道很长，气体与外界能够进行充分的热交换，使管道中气流与周围环境保持相同温度，这是等温管流。工程流体力学工程流体力学1
37、 1几个实用公式几个实用公式 范诺流的气体动力学函数 2*122*122*121211111121TTMaMaTMaTMappMaMa 若管流两截面压强分别是 和 ，那末可推导出管流的质流量 为1p2pmQ1225121 11111kg s81mppdQV AlRTp工程流体力学工程流体力学式中 为管道沿程摩阻因数，对于亚声速流，可直接利用穆迪图确定。对于超声速流（）通常取 。对于绝热摩擦管流也同样存在一个最大长度 ，若气流在该处已达到临界状态，当实际管长 时将发生壅塞现象。此时对亚声速流造成的压强扰动可向上游传播至入口，使入口处溢流而造成流量减小直至出口截面正好为临界截面。对超声速流，壅塞在
38、管中产生激波，从而使临界截面移至出口截面。13Ma0.002 0.003maxlmaxll最大管长为22max22111ln221MadMalMaMa工程流体力学工程流体力学管道进、出口马赫数的关系式为222122212211111112ln1212MaMaldMaMaMaMa2.2.计算实例计算实例 【例10.12】马赫数为 的空气超声速气流进入一个沿程损失因数 的绝热管道,管道的直径 ，若要求管道出口马赫数 ，试求管道长度l。13Ma 0.02200mmd 22Ma 解：由式（10.40），对于空气 ，按题意 1.4120.02,0.2m,3,2dMaMa则 1112.4910.2 4ln
39、1.4 942.8410.2 9ld得 2.171 6ml 工程流体力学工程流体力学解：（1）进口处马赫数 1112000.5551.4 287 323VMaRT现要求最大管长，即出口截面为临界状态 ，由式（10.40）21Ma 22max220.11 0.5551.4 1(1.4 1)0.555ln0.025 1.4 0.5552 1.42(1.4 1)0.555l2.80m工程流体力学工程流体力学【例10.13】空气流在等截面管道作绝热摩擦流动，进口处状态参数 ，。若管径 ，沿程摩阻因数 ，试求：（1）最大管长 及出口处压强和温度（2）若管长为 ,试计算进口的马赫数。512 10 Pap
40、1323KT 1200m sV 100mmd 0.025maxl3.5ml 出口处温度由式10.37（a）2*121TTMa2*2(1.4 1)0.555323285.75K1.4 1T 出口处压强由式10.38（c）125*21.4 12 100.5552(1.4 1)0.555p51.04 10 Pa（2）由于按照初始条件，最大管长为 ，现实际管长为 ，因此将发生壅塞现象。max2.8ml3.5ml 工程流体力学工程流体力学即出口处达到临界状态，进口处的 不再保持0.555，而将由式（10.41）确定，式中 代入该方程得 1Ma21,Ma 3.5m,d=0.1ml 0.025,1.4212
41、211112.40.8751ln1.21.42.81 0.2MaMaMa令 ，则上面代数方程经化简为 211xMa ln(0.2)1.671 801.2xfxx或者 1.2 ln 0.21.671 8xx工程流体力学工程流体力学Ma=0.555Ma2=1Ma1llm图 0.12 1气流进入绝热摩擦管设初始马赫数 ，即 ，应用简单迭代法解得 ,因此进口处马赫数 ,该流动的实际情况是当管长 时，进口处亚声速气流发生膨胀减速，马赫数由0.555减小到0.526后才进入绝热摩擦管 如图10.12。10.5Ma 04x 3.611 9x 10.526Ma 3.5m2.8mmll工程流体力学工程流体力学第第10章章 结结 束束
