1、拉伐尔喷管的设计摘 要：本文针对拉伐尔喷管的几何条件和力学条件进行了推导。建立了喷管截面积变化与流速、压强、密度、温度等流动性能参数间的关系，分析了喷管出口截面下游的外界反压对拉伐尔喷管工作过程的影响。推导建立了拉伐尔喷管主要性能参数的计算方法。针对实际流动损失的存在，为得到喷管的实际流动性能，对理论性能参数提出了修正方法。提出了拉伐尔喷管的设计方法。本文研究内容为拉伐尔喷管的设计提供依据。关键词：变截面；力学条件；性能参数；流动损失1引言拉伐尔喷管是火箭发动机和航空发动机最常用的构件，由两个锥形管构成，如图 1 所示，其中一个为收缩管，另一个为扩张管。拉瓦尔喷管是推力室的重要组成部分。喷管的
2、前半部是由大变小向中间收缩至喷管喉部。喉部之后又由小变大向外扩张。燃烧室中的气体受高压流入喷嘴的前半部，穿过喉部后由后半部逸出。这一架构可使气流的速度因喷截面积的变化而变化，使气流从亚音速到音速，直至加速至超音速。所以，人们把这种喷管叫跨音速喷管。瑞典工程师 De Laval 在 1883 年首先将它用于高速汽轮机，现在这种喷管广泛应用于喷气发动机和火箭发动机。图 1.1 拉伐尔喷管结构图2拉伐尔喷管的几何条件21 变截面一维定常等熵流动在变截面一维定常流动中只考虑截面积变化这一种驱动势，忽略摩擦、传热、重力等其他驱动势，因此流动是绝热无摩擦的，即等熵流动，变截面定常等熵流动模型如图 2 所示
3、。控制体p+dpdx +dV+dVT+dTA+dApVTA图2.1 变截面一维定常等熵流动模型变截面一维定常等熵流动的控制方程组为：(2.1.1)ConstmVA(2.1.2)0dp(2.1.3)210dhV22 截面积变化对流动特性的影响管道的形状变化可以用截面积变化 dA 来表示。(a) 截面积变化对流速的影响对连续方程(1)取对数微分，得(2.2.1)0dVA将(2.1.2)两边同除以 ，得(2.2.2)2pd由声速公式及马赫数定义，得(2.2.3)21VAM这就是截面积变化与流速变化之间的关系。(b) 截面积变化对压强的影响将(2.1.2)代入(2.2.3)，由理想声速公式得到(2.2
4、.4)21dpdA(c) 截面积变化对密度、温度、声速、马赫数的影响联立(2.2.1)式与(2.2.3)式，消去速度项，得(2.2.5)21dMdA联立(2.1.2)式与(2.1.3)式，并将(2.2.4)式代入，得(2.2.6)2T将理想气体声速公式求对数微分，并将(2.2.6)式代入，得到(2.2.7)21dadAM对马赫数定义取对数微分，并将(2.2.3)式和(2.2.7)式代入，得(2.2.8)21ddA通过分析所得结果，截面积变化对各流动特性的影响可概括为：一维定常等熵流动具有膨胀加速或压缩减速额流动特性。收敛管道中的亚声速流和扩张管道中的超声速流是膨胀加速的，沿管道流速不断增加，而
5、压强、密度和温度不断减小；扩张管道中的亚声速流和收敛管道中的超声速流是压缩减速的，沿流道流速不断降低，而压强、密度和温度却不断增加。23 流动极限状态壅塞状态收敛管道中的一维定常等熵流动流速只能连续变化到 M=1，即达到临界状态，这是它的极限。在此之后，流速既不可能增大，也不可能减小，收敛管道中的这种现象称为流动壅塞。同样，超声速流也不可能通过收敛管道连续减速到亚声速流。如果在临界截面之后使管道扩张，则当管道出口截面处的下游物理边界条件满足一定要求时，流动能够从声速流变为超声速流。这种先收敛后扩张的管道即为拉伐尔喷管。这种先收敛后扩张的管道形状是从初始亚声速流获得超声速流的必要条件，称为拉伐尔
6、喷管的几何条件。3拉伐尔喷管的力学条件拉伐尔喷管为实现亚声速流向超声速流的连续变化，除几何条件外，必须对喷管出口截面下游的环境压强（外界反压）做出限制，即拉伐尔喷管的力学条件。为了分析外界反压对拉伐尔喷管流动的影响，假设出口截面外的环境压强保持不变，而喷管进口截面的滞止压强 可变。当总压 变化时，喷管出ap 0p0p口截面上的气体压强 随之变化。根据 和 的相对大小，气体在喷管中的流epae动状态分为以下三种情况。(1) 最佳膨胀状态 ea气体在喷管中得到了完全膨胀，这就是喷管的最佳膨胀状态，又称为设计状态，如图 3.1 所示。这种流动的主要特点是：喷管喉部达到了临界状态，出口流动为超声速，即
7、 Me1；流体流出喷管后，既不膨胀，也不压缩，而是一平行射流；由于管内流动为超声速，当外界环境发生微小扰动时，扰动的传播速度（即声速）小于流动速度，扰动不能传进喷管内部，即喷管中的流动觉察不到外界反压的变化。出口截面进口截面pt Pe=papepepepaMae1ptpa进口截面 出口截面图 3.1 喷管最佳膨胀时的流动 图 3.2 欠膨胀状态时的喷管流动(2) 欠膨胀状态 eap如果在最佳膨胀状态下提高喷管进口总压 ，则出口 同时增大，有0pe。气体没有得到完全膨胀，其能量未充分发挥，即气体热能没有最大限eap度地转变成定向流动动能。这种流动称为欠膨胀状态或膨胀不足状态，如图 3.2所示。欠
8、膨胀状态流动主要特点是：喷管喉部达到了临界状态，出口仍为超声速 M1；气体在喷管外继续膨胀，直到压强等于 时为止，因此喷管出口处有一ap系列膨胀波；喷管外的压强扰动也不能逆向传入喷管。(3) 过膨胀状态 eap如果在最佳膨胀状态下减小喷管进口总压 ，则喷管出口的气体压强也将0p减小，即 。气体在喷管中作了过分的膨胀。这种流动称过膨胀状态。根ea据 小于 的程度大小，气体在喷管中的流动状态又可分为下述四种情况。epa 稍小于 ap喷管出口的气体流动为超声速。在喷管外气体由于受到反压的突然压缩而产生不连续的压强增加，形成激波。因为 稍小于 ，激波是附着在扩张段出epa口截面上的激波，如图 3.3
9、所示。气体经过斜激波后，压强升高到 。ap 比 小于一定值epa随着压强差 的增大，喷管外的斜激波逐渐向喷管口收拢，并最终在ep小于 一定值时演变成覆盖在喷管出口截面上的正激波，如图 3.4 所示。气ea体压强 经过正激波压缩后升高到 ，这时的外界反压 称为第二临界反压。e apap 进一步减小ep当 比 小很多时，正激波从喷管出口截面向喷管内部移动，喷管扩张段a内的流动以正激波为分界线。激波后的流动就是扩张管道中的亚声速流动，流动的马赫数将逐渐减小，压强逐渐升高，并在喷管出口截面升高到 。ap eap=如果 ，则正激波最终移动到喉部。此时正激波消失，流动不再壅塞，全部喷管内的流动均为亚声速流
10、，气体的压强、流速和质量流率都为外界反压所控制。这种流动状态称为亚临界流动状态，喷管喉部达不到临界状态。pepaM1斜激波pt进口截面 出口截面 进口截面 出口截面正激波pe1ptpe2=paMe21pe2=pa图 3.3 过膨胀状态的喷管流动 图 3.4 正激波位于喷管出口截面时的流动综上所述，若要在拉伐尔喷管出口截面获得超声速气流，喷管出口截面的气体压强必须达到或超过反压值，这一条件称为力学条件。由此可知，拉伐尔喷管中的流动受几何条件和力学条件两方面的影响，在拉伐尔喷管的设计过程中必须同时考虑。4拉伐尔喷管的性能参数计算拉伐尔喷管的性能参数主要包括喷管出口速度(排气速度)、质量流率、推力等
11、。(1) 流速 和排气速度Ve根据式(2.1.3) 可求出喷管内任一截面流速，即(4.1)0022()(1)p TVcTR将气体动力学函数代入(4.1)式，得到(4.2)102VRT令流速函数为(4.3)1VF得到流速：(4.4)0vRT排气速度：(4.5),etVeF排气速度衡量火箭发动机性能高低的一个重要参数。火箭是利用燃气的高速喷出获得推动力的，所以排气速度越高，获得的推力就越大。(2) 质量流率当喷管的喉部截面达到临界状态时，临界截面积就是喉部截面积 At。质量流率为：(4.6)tpAmRT(3) 推力 F对于给定的固体火箭推进剂， 和 RT 均为常数，当喷管进口总压 和喷 0p管喉部
12、面积一定时，质量流率是确定的，火箭发动机的推力仅是扩张比的函数。对于一定值的环境压强，由于喷管处于最佳膨胀状态所具有的最大推力为最佳推力，用 表示，即opt(4.7)toptetVtpAFmVRTFpA (4) 推力系数推力系数 定义为FC(4.8)0FtCpA5拉伐尔喷管中的流动损失上述对拉伐尔喷管流动的讨论基于一维定常等熵流动假设，实际流动过程与这种理想情况存在一定差别，如喷管热损失、摩擦损失、非理想气体效应等。在拉伐尔喷管设计中，为得到喷管的实际性能，必须对理论性能参数进行修正。通常需要考虑的流动损失主要包括：两相流损失、流量损失、边界层损失和喷管扩张损失等，这些流动损失一般用修正系数来
13、表示，定义为 实 际 性 能 参 数理 论 性 能 参 数(1) 两相流损失在理想性能参数计算中，假设燃气为理想气体，实际上，应当考虑凝聚相影响。计算两相流损失的修正系数用如下经验公式(5.1)32410.5.8cstptetnd式中 凝聚相微粒浓度；sn凝聚相微粒直径；d燃烧室压强；tp喷喉直径；面积扩张比。2e(2) 流量损失燃气从喷管收敛段流到喷管喉部时，由于气流的惯性作用，流线不能完全适应流道截面变化，使实际喷喉直径减小，使流量下降。流量损失系数为(5.2)210.2.3tmgA 式中 喷喉面积；tA挡药板通气面积，无挡药板时可使用燃烧室末端面积；g喷管的收敛半角。(3) 边界层损失边
14、界层损失是指喷管壁面摩擦及散热所造成的损失，修正系数为(5.3)0.820.8121.69tcpdtbl ece式中 t发动机工作时间；c1、c2与喷管有关的常数，一般取 c1=0.2357，c2=0.0605。(4) 喷管扩张损失在理想喷管流动中，一维流动的方向平行于轴线，流动参数在垂直于轴线截面上均匀一致。实际上，燃气沿锥形向外扩张流动，这种流动更接近源流，即所有流线从源点出发向外扩张流动。火箭喷管流动计算主要关心出口截面即排气面上的流动参数。采用源流假设时，排气参数只有在球面 上才是均匀一致的，分别为压强 、密度 、速sAsps度 等，而喷管出口截面 上的参数仍用 、 、 等表示。sve
15、epevrAe AsdFO cb图 5.1 喷管源流示意图考虑喷管扩张损失，推力公式可改写成(5.4)1cos2eeaeeaFmVApmVAp 修正系数(5.5)1cos26. 拉伐尔喷管设计拉伐尔喷管设计包括亚声速收缩段、喉部和超声速扩张段 3 个部分的型线设计。(1)亚声速收缩段 亚声速收缩段的作用是使气流加速，同时要保证收缩段的出口气流均匀、平直而且稳定。收缩段的性能取决于收缩段进口面积和出口面积的比值及收缩段曲线形状。将收缩段设计成维托辛思基曲线 5，收敛段上任意截面 半径为：,22223111()(/)/()RRxLxL式中：R1 、R2、R 分别为收缩段进口、出口及任意 x 处 的
16、截面半径，L1为收缩段长度（如图 6.1 所示）图 6.1 收缩段型线图（2）喉部 喉部是气流从亚声速转变为超声速的过渡段，喉部直径的选取受到气流流量的限制。喉部曲线变化不能太快，这里选用了一段圆弧作为过渡曲线。（3）超声速扩张段 扩张段曲线采用基于特征线法的富尔士法进行设计 【5】 。超声速扩张段曲线包括 3 段曲线（如图 6.2） 。喉部过渡段、直线段和消波段，其中喉部过渡段和直线段是使气流加速的，消波段是设法将膨胀波在壁面的反射消灭，以保证实现出口气流均匀。6.2 扩张段曲线形状图结束语本文从原理上对拉伐尔喷管设计中的几何条件以及力学条件分别进行了推导与分析。在火箭发动机拉伐尔喷管的设计
17、过程中，必须同时满足先收敛后扩张的几何条件以及喷管出口截面的气体压强必须达到或超过反压值的力学条件。列举了拉伐尔喷管各性能参数的计算方法。由于实际运用中拉伐尔喷管存在流动损失，本文最后引入针对各流动损失的修正方法。参考文献：1 吴望一. 流体力学M. 北京：北京大学出版社， 20102 武晓松，陈军，王栋. 固体火箭发动机气体动力学M. 北京：国防工业出版社，20053 董师颜等. 固体火箭发动机原理M. 北京: 国防工业出版社,19834 刘志友等. 环境压力大偏差条件下拉瓦尔喷管发动机高空推力的确定 J. 航空动力学报，2006，21（3）：610-6145 刘建华，陈一坚，段军，等激光切割超音速喷嘴设计J. 激光技术，2000,24（1）：46-50
