喷管流动特性与管道截面变化规律的关系
小型拉瓦尔喷管的参数设计和模拟

㊀㊀辽宁大学学报㊀㊀㊀自然科学版第50卷㊀第2期㊀2023年JOURNALOFLIAONINGUNIVERSITYNaturalSciencesEditionVol.50㊀No.2㊀2023小型拉瓦尔喷管的参数设计和模拟范晓星ꎬ王㊀宇ꎬ王曦珠ꎬ张㊀硕(辽宁大学物理学院ꎬ辽宁沈阳110036)摘㊀要:冷喷涂相比于传统热喷涂技术有许多优势ꎬ能够制备传统喷涂技术难以制成的材料涂层ꎬ能更好地适应工作环境ꎬ使加工过程方便快捷ꎬ从而提高工作效率.本文以冷喷涂设备小型化为研究背景ꎬ首先在理论上对拉瓦尔喷管原理进行分析ꎬ建立喷管截面积变化与各流动性能参数间的关系ꎬ根据理论推导得出拉瓦尔喷管截面积与出口处流体速度等流动性能参数的关系ꎬ最后通过计算机模拟了各参数下拉瓦尔喷管的工作状态.关键词:拉瓦尔喷管ꎻFluent模拟ꎻ超音速流中图分类号:TH138.5㊀㊀㊀文献标志码:A㊀㊀㊀文章编号:1000-5846(2023)02-0146-08DesignofSmallSizeLavalNozzleandSimulationFANXiao ̄xingꎬWANGYuꎬWANGXi ̄zhuꎬZHANGShuo(SchoolofPhysicsꎬLiaoningUniversityꎬShenyang110036ꎬChina)Abstract:㊀Comparedwithtraditionalthermalsprayingtechnologyꎬcoldsprayinghasmanyadvantages.Itcanpreparematerialcoatingsthataredifficulttobemadebytraditionalsprayingtechnologyꎬbetteradapttotheworkingenvironmentꎬmaketheprocessingprocessconvenientandfastꎬandimprovetheworkefficiency.TakingtheminiaturizationofcoldsprayingequipmentastheresearchbackgroundꎬthispaperfirstanalyzestheprincipleofLavalnozzleintheoryꎬestablishestherelationshipbetweenthechangeofnozzlecross ̄sectionalareaandvariousflowperformanceparametersꎬdeducestherelationshipbetweenLavalnozzlecross ̄sectionalareaandflowperformanceparameterssuchasfluidvelocityattheoutletaccordingtothetheoryꎬandfinallysimulatestheworkingstateofLavalnozzleundervariousparametersbycomputer.Keywords:㊀LavalnozzleꎻFluentsimulationꎻsupersonicflow0㊀引言冷喷涂技术可实现涂料在超音速状态下与基底的牢固接触ꎬ实现对待喷涂表面的均匀顺滑的喷㊀收稿日期:2022-04-12作者简介:范晓星(1980-)ꎬ男ꎬ吉林长春人ꎬ博士ꎬ教授ꎬ研究方向:流体流场计算模拟ꎬ光催化气固反应器设计.㊀㊀涂覆盖[1].冷喷涂过程中温度较低ꎬ相对于传统的热喷涂技术ꎬ冷喷涂技术可以更好地实现对特殊材料的喷涂ꎬ有效降低喷涂材料的氧化ꎬ并且减少相变等情况的发生ꎬ具有诸多优势[2].冷喷涂技术的核心是采用拉瓦尔喷管对载送涂料的气体进行超音速加速ꎬ在喷管出口位置将颗粒物质加速到超音速[3].拉瓦尔喷管是利用出入口的压强差ꎬ通过截面积先减小后增大ꎬ来实现对气体的超音速加速ꎬ该过程中高压气体温度逐渐降低ꎬ压强逐渐减小ꎬ此过程中气体的内能转变为气体的定向动能ꎬ气体速度逐渐增大.冷喷涂设备在工业领域已经被广泛应用ꎬ但是其设备庞大ꎬ无法在实验室条件下应用.冷喷涂技术在实验室条件下可用于小型的实验研究和小尺寸的样品制备ꎬ从而降低实验的成本ꎬ因此冷喷涂设备的小型化具有一定的应用前景[4].本文针对冷喷涂设备小型化的需求ꎬ开展了拉瓦尔喷管的设计研究ꎬ介绍了拉瓦尔喷管的工作原理ꎬ对拉瓦尔喷管内流体运动过程进行简要分析ꎬ建立喷管截面积变化与各流动性能参数间的关系ꎬ计算得出拉瓦尔喷管的相关数据ꎬ并利用Fluent软件进行了计算机模拟.1㊀拉瓦尔喷管的工作原理分析为了设计针对气体的加速装置ꎬ本文首先研究一元气体的流动特性[5].流体在流管内的流动过程可以视为定常等熵流动[6]ꎬ气体状态满足如下的方程.连续性方程[7]:ρvA=C(常数)(1)能量方程:dH+dv22æèçöø÷=0(2)动量方程:dpρ+vdv=0(3)状态方程:pV=nRT(4)其中:H为气体的焓ꎻρ为气体密度ꎻp为压强ꎻv为气体速度ꎻV为气体体积ꎻR为气体常数ꎻT为温度.由式(1)~式(4)可以计算出流管内气体速度㊁压强㊁密度㊁温度随截面积的变化规律.dvv=1Ma2-1dAA(5)dpp=kMa21-Ma2dAA(6)dρρ=Ma21-Ma2dAA(7)dTT=(k-1)Ma21-Ma2dAA(8)其中:Ma为喉部位置的马赫数ꎻA为截面面积.式(5)~式(8)给出了流体管道中速度㊁压强㊁密度和温度的变化过程与截面面积变化的对应关系.通过其中速度与截面面积的变化关系ꎬ我们可以看出亚声速流体的速度随截面积的减小而增加ꎬ741㊀第2期㊀㊀㊀㊀㊀㊀范晓星ꎬ等:小型拉瓦尔喷管的参数设计和模拟㊀㊀而超声速流体的速度随截面积的增加而增加.在流体速度接近于声速时ꎬ应有dA=0ꎬ即截面积取最小值[8]ꎬ也就是说ꎬ在喉部应实现马赫数Ma=1.为了将气体从亚声速加速到超声速ꎬ人们通常采用拉瓦尔喷管实现超音速气体喷射.拉瓦尔喷管的形状如图1所示ꎬ形如细腰筒形鼓ꎬ两头大中间小ꎬ渐缩渐扩ꎬ以中间喉管为界限ꎬ分为渐缩段和渐扩段.由上述公式可知ꎬ流体通过渐缩段进入喉部的过程中ꎬ马赫数小于1㊁截面积减小㊁压强减小㊁密度减小㊁温度降低㊁气流速度增大ꎬ马赫数逐渐增大ꎬ到喉部时达到Ma=1ꎻ在进入渐扩段后ꎬ截面积增大㊁压强减小㊁密度减小㊁温度减小㊁气流速度进一步增大ꎬ预计马赫数在出口处将达到大于2的水平.气体流动参数与流管截面积变化的关系ꎬ管内气体压强㊁密度㊁温度㊁速度随截面积的变化情况如表1所示.图1㊀拉瓦尔喷管形状示意图表1㊀气体流动参数的变化参数渐缩段渐扩段压强渐减渐减密度渐减渐减温度渐减渐减速度渐增渐增2㊀拉瓦尔喷管参数设计2.1㊀计算过程拉瓦尔喷管的各处尺寸标记如图2所示.图2㊀拉瓦尔喷管尺寸示意图㊀㊀气体在拉瓦尔喷管内流动的过程中ꎬ具有很高的速度ꎬ可忽略与外界发生的热交换ꎬ依据绝热过程来处理分析.符号说明如下:假设气体为氮气ꎬ气体常数用R表示ꎬR=0.296kJ/(kg K)ꎻ绝热指数:k=1.41ꎻ滞止密度:ρ0ꎻ喉部密度:ρcrꎻ滞止压强:p0=p1T0T1æèçöø÷kk-1ꎻ出口压强:p2=101325Paꎻ滞止温度:T0=T1+v212cpꎬ其中ꎬcp为比定压热容.在入口处ꎬ气体的速度要远低于管内的速度ꎬ可近似为滞止状态[9].已知气体的能量方程为[10]841㊀㊀㊀辽宁大学学报㊀㊀自然科学版2023年㊀㊀㊀㊀H+v22=const(9)代入H=cpTꎬcp=kk-1Rꎬ式(9)可转化为kk-1RT+v22=const(10)又引入气体状态方程pρ=RTꎬ气体能量方程可转化为kk-1pρ+v22=const(11)气体在管内流动的能量损失忽略不计ꎬ可得出ꎬ入口处气体能量与流动过程中任意位置的能量守恒关系式ꎬ如下:kk-1p0ρ0=kk-1pρ+v22(12)由式(12)可得气体流动至各处的速度:v=㊀2kk-1p0ρ01-pp0ρ0ρæèçöø÷(13)考虑到绝热过程中ρ0ρcr=p0pcræèçöø÷1k(14)式(13)可转化为v=㊀2kk-1p0ρ01-pp0æèçöø÷k-1k[](15)根据式(15)ꎬ得到出口处速度:v2=㊀2kk-1p0ρ01-p2p0æèçöø÷k-1k[](16)再次代入气体状态方程pρ=RTꎬ有v2=㊀2kk-1RT01-p2p0æèçöø÷k-1k[](17)图3 出口处速度三维云图通过式(17)可以得出ꎬ拉瓦尔喷管出口处的速度ꎬ只和初始状态下的温度与压强相关.根据式(17)可绘制出初始状态温度在0~600ħꎬ压强在0.1~2.025MPa时ꎬ对应的出口处流速三维云图ꎬ如图3所示ꎬ可以看出高温高压的情况下气体速度更高ꎬ在低压段压强变化会导致气体速度增加得较快ꎬ高压段压强变化会导致气体速度增加得较慢.提高载气温度也可以提高喷口处气体的速度.考虑到实验室条件下拉瓦尔喷管小型化的需求ꎬ我们对拉瓦尔喷管的喉部尺寸进行了估算.实验室能够提供气体的流量一般在0.04~0.36m3/minꎬ根据941㊀第2期㊀㊀㊀㊀㊀㊀范晓星ꎬ等:小型拉瓦尔喷管的参数设计和模拟㊀㊀气源的气体流量数据ꎬ以及气体在拉瓦尔喷管喉部的密度数据可以算出喉部半径范围为0.80~2.38mm.当喉部直径为1.6mm时ꎬ我们对拉瓦尔喷管参数进行了估算.一般情况下ꎬ拉瓦尔喷管渐缩段的长度l1的取值为(3~5)dcrꎬ又由于转折角度过小ꎬ渐缩段的长度会过于长ꎬ影响实际应用ꎬ而转折角度过大ꎬ渐缩段的长度会过短ꎬ在转折处会使实际喉部截面变小ꎬ影响加速效果[11].综合考虑ꎬ可以取α=30ʎꎬ则渐缩管入口截面直径d1=2l1tan15ʎ+dcr.若喉部截面直径dcr=1.6mmꎬ可求得d1=2l1tan15ʎ+dcr=4.17~5.89mm.为了防止气体产生紊流ꎬ在喉部需建立过渡段[12]ꎬ该过渡段的计算公式为lcr=(0.5~1)ˑdcr.这里取lcr=1.6mm.对于渐扩段ꎬ当喷管渐扩角θ在8ʎ~12ʎ范围[13]内变化时对流场的影响不明显ꎬ可取θ=8ʎꎬ得到出口截面直径d2即可推出l2.对于d2的求解ꎬ可利用拉瓦尔喷管任意截面面积与喉管截面积之比与马赫数的关系求解:AAcr=1Ma1+k-12Ma2k+12éëêêêêùûúúúúk+12(k-1)(18)2.2㊀已知参数根据现有实验环境ꎬ本文的气体条件如下:入口压强:p1=1.5/0.8MPaꎻ入口温度:T1=873.15Kꎻ喷管背压:pb=101325Paꎻ出口压强:p2=pb=101325Pa.实验气体为氮气ꎬ则气体常数:R=0.296kJ/(kg K)ꎻ绝热指数:k=1.41(k=cpꎬmcvꎬmꎬ对于双原子分子气体ꎬk=1.41)ꎻ比定压热容:cp=1.0392kJ/(kg K).已知滞止压强:p0=p1T0T1æèçöø÷kk-1ꎻ滞止温度[14]:T0=T1+v212Cp.2.2.1㊀p1=0.8MPa时的数据计算当拉瓦尔喷管入口压强为0.8MPa时ꎬ根据式(17)可得到出口速度为v2=㊀2kk-1p0ρ01-p2p0æèçöø÷k-1k[]=896.04m/sꎬ代入Ma=v2a(a为当地声速)可得出口处马赫数Ma=2.64ꎬ进一步得AAcr=2.96.根据面积与直径关系d=㊀4Aπꎬ可得出口处截面直径d2=2.752mm.2.2.2㊀p1=1.5MPa时的数据计算当拉瓦尔喷管入口压强为1.5MPa时ꎬ根据式(17)可得到出口速度为051㊀㊀㊀辽宁大学学报㊀㊀自然科学版2023年㊀㊀㊀㊀v2=㊀2kk-1p0ρ01-p2p0æèçöø÷k-1k[]=982.72m/sꎬ代入Ma=v2a(a为当地声速)可得出口处马赫数Ma=2.89ꎬ进一步得AAcr=3.76.根据面积与直径关系d=㊀4Aπꎬ可得出口处截面直径d2=3.101mm.2.3㊀数据结果在已有实验条件基础上ꎬ可计算出拉瓦尔喷管的各项参数(表2).表2㊀拉瓦尔喷管几何结构参数入口压强/MPad1/mmd2/mmdcr/mmp1=0.85.032.751.60p1=1.55.033.101.60㊀㊀在实验室条件下ꎬ根据以上计算得出的数据ꎬ可制作出小型的拉瓦尔喷管以完成实验研究.图4㊀拉瓦尔喷管模型3㊀Fluent计算机模拟本文对已经获得的拉瓦尔喷管的参数使用流体仿真软件Fluent进行模拟仿真ꎬ如图4所示ꎬ首先建立拉瓦尔喷管模型ꎬ并将模型转移至前处理软件ICEMCFDꎬ仿真模拟中ꎬ流场计算模型的计算域应包括渐缩段㊁喉部㊁渐扩段[15].图5㊀拉瓦尔喷管网格划分模型3.1㊀喷管网格划分ICEMCFD软件是一种专业的前处理软件ꎬ能够为Fluent软件提供高效可靠的分析模型ꎬ利用其对喷管模型进行轴向方向上的 O 型网格划分ꎬ在对棱角处网格节点数进行调节等操作后ꎬ最终网格质量均在0.75以上[16]ꎬ如图5所示.3.2㊀Fluent仿真模拟Fluent是一款可以分析模拟压缩或不可压缩流体与离散相之间的耦合流动的大型流体仿真软件.本文对气体在喷管内的流动过程进行仿真模拟ꎬ采用线性压力-应力修正的雷诺应力模型(RSM)ꎬ辅以变尺度壁面函数处理近壁区[17]ꎬ据公式Re=ρvd/μ得到此次模拟适宜的雷诺数为5.本次模拟对象为可压缩且加速后达到超音速的气体ꎬ故选用密度场进行模拟ꎬ采用精密程度更高的二阶迎风方程ꎬ设置连续性残差小于10-6时认定收敛ꎬ停止计算[18].气体入口温度为873.15Kꎬ出口压强为101325Paꎬ选择压强边界条件ꎬ不计气体流动过程中与外界的热交换㊁摩擦ꎬ内壁选择无滑移㊁绝热边界[19]ꎬ分别对气体进口压强为0.8MPa㊁1.5MPa情况进行模拟[20].当入口压强为0.8MPaꎬ温度为600ħ时ꎬ气体速度㊁温度㊁压强㊁密度场模拟结果如图6所示ꎬ经过拉瓦尔喷管加速ꎬ气流到喉部达到当地音速ꎬ在喷嘴出口处达到953.04m/s.气体温度㊁压强㊁密度逐渐降低.当入口压强为1.5MPaꎬ温度为600ħ时ꎬ气体速度㊁温度㊁压强㊁密度场模拟结果如图7所示ꎬ经过拉瓦尔喷管加速ꎬ气流到喉部达到当地音速ꎬ在喷嘴出口处达到992.62m/s.气体温度㊁压151㊀第2期㊀㊀㊀㊀㊀㊀范晓星ꎬ等:小型拉瓦尔喷管的参数设计和模拟㊀㊀强㊁密度逐渐降低.图6㊀入口压强为0.8MPa时的模拟结果图7㊀入口压强为1.5MPa时的模拟结果4㊀结论本文对小型化拉瓦尔喷管的基本结构进行了设计ꎬ获得了拉瓦尔喷管的基本参数ꎬ可实现在家用气泵和低加热功率的条件下超声气体喷射.基于Fluent软件对气体在拉瓦尔喷管中的流动进行数251㊀㊀㊀辽宁大学学报㊀㊀自然科学版2023年㊀㊀㊀㊀值模拟仿真ꎬ结果显示ꎬ气体在喷嘴出口处速度能够达到超音速ꎬ且能达到2倍音速以上.喷枪设计具有可行性.数据表明:1)喷管模拟初始条件:入口温度为873.15Kꎬ出口压强为101325Paꎬ气体入口压强分别为0.8MPa㊁1.5MPa时ꎬ喷管出口处气体马赫数超过2ꎬ速度分别为953.04m/s㊁992.62m/s.2)通过数值模拟可以将气体在喷管内流动过程中各参数变化以图像形式直观展现ꎬ为喷枪的实物设计提供了参考数据.3)根据模拟结果可以发现ꎬ小型化冷喷涂喷枪对气体的加速效果与工业化冷喷涂技术的喷涂效果不相上下ꎬ可在一定的条件下使用ꎬ从而使冷喷涂过程变得更加便捷.参考文献:[1]㊀DeberneNꎬLeoneJFꎬDuqueAꎬetal.Amodelforcalculationofsteaminjectorperformance[J].InternationalJournalofMultiphaseFlowꎬ1999ꎬ25(5):841-855.[2]㊀赵国锋ꎬ王莹莹ꎬ张海龙ꎬ等.冷喷涂设备及冷喷涂技术应用研究进展[J].表面技术ꎬ2017ꎬ46(11):198-205.[3]㊀SherifSꎬLearWꎬSteadhamJꎬetal.Analysisandmodelingofatwo ̄phasejetpumpofathermalmanagementsystemforaerospaceapplications[C]//36thAIAAAerospaceSciencesMeetingandExhibit.Reno:AIAAꎬ1998:360.[4]㊀TakeuchiKꎬCunninghamJPꎬHochreiterLꎬetal.Steam ̄jetpenetrationintofreesurfaceofsubcooledwater[J].ThermalHydraulicsofAdvancedSteamGeneratorsandHeatExchangersꎬ1994ꎬ15:19-28.[5]㊀刘旭.L型双级谐振腔式低频超声雾化喷头的设计及试验[D].镇江:江苏大学ꎬ2018.[6]㊀焦峥辉.超音速喷管雾化器设计及理论研究[D].西安:西安石油大学ꎬ2017.[7]㊀钱翼稷.空气动力学[M].2版.北京:高等教育出版社ꎬ2008.[8]㊀文怀兴ꎬ刘永强ꎬ党新安.拉伐尔喷管曲线方程的建立及其数控加工[J].西北轻工业学报ꎬ1997ꎬ15(4):7-11.[9]㊀王冠群.拉瓦尔超音速雾化喷头的设计及试验[D].镇江:江苏大学ꎬ2020.[10]㊀任鹏宇.带拉瓦尔管的高压气体引射喷水泵特性实验研究[D].哈尔滨:哈尔滨工业大学ꎬ2008.[11]㊀CumoMꎬFarelloGEꎬFerrariG.Heattransferincondensingjetsofsteaminwater[C]//Proceedings6thInternationalHeatTransferConference(IHTC).Tornto:IHTCꎬ1978:101-106.[12]㊀MunsonBRꎬYoungDFꎬOkiishiTHꎬ等.工程流体力学[M].北京:电子工业出版社ꎬ2006.[13]㊀李文科.工程流体力学[M].合肥:中国科学技术大学出版社ꎬ2007.[14]㊀王新月ꎬ杨清真.热力学与气体动力学基础[M].西安:西北工业大学出版社ꎬ2004.[15]㊀ApteSVꎬGorokhovskiMꎬMoinP.LESofatomizingspraywithstochasticmodelingofsecondarybreakup[J].InternationalJournalofMultiphaseFlowꎬ2003ꎬ29(9):1503-1522.[16]㊀高全杰ꎬ汤红军ꎬ汪朝晖ꎬ等.基于Fluent的超音速喷嘴的数值模拟及结构优化[J].制造业自动化ꎬ2015ꎬ37(4):88-90ꎬ108.[17]㊀杨超ꎬ陈波ꎬ姜万录ꎬ等.基于拉瓦尔效应的超音速喷嘴雾化性能分析与试验[J].农业工程学报ꎬ2016ꎬ32(19):57-64.[18]㊀王冰川ꎬ张凯ꎬ张聃ꎬ等.基于Fluent的超音速气液混合喷嘴模拟仿真[J].清洗世界ꎬ2020ꎬ36(1):35-37.[19]㊀ArientiMꎬWangLꎬCornMꎬetal.Modelingwallfilmformationandbreakupusinganintegratedinterface ̄tracking/discrete ̄phaseapproach[C]//ASMETurboExpo2010:PowerforLandꎬSeaꎬandAir.Glasgow:ASMEꎬ2010:1059-1068.[20]㊀焦峥辉ꎬ王美妍ꎬ任雪娇ꎬ等.超音速雾化喷枪的设计与分析[J].现代制造技术与装备ꎬ2019(12):79-80ꎬ95.(责任编辑㊀郑绥乾)351㊀第2期㊀㊀㊀㊀㊀㊀范晓星ꎬ等:小型拉瓦尔喷管的参数设计和模拟。
液体火箭发动机喷管设计与流动特性分析

液体火箭发动机喷管设计与流动特性分析液体火箭发动机喷管是发动机中最关键的部件之一,它直接影响着发动机的性能和效果。
因此,对喷管的设计和流动特性分析十分重要。
本文将从设计原理、喷管构造、流动特性分析等方面进行分析,以期全面了解液体火箭发动机喷管的设计与流动特性。
液体火箭发动机喷管的设计原理是基于牛顿定律和热力学原理。
发动机喷管中的燃料和氧化剂通过压力调节阀进入喷管的燃烧室,然后在燃烧过程中产生大量的燃气。
这些燃气通过喷管的喷嘴部分喷出,产生反作用力推动火箭运行。
因此,喷管的设计需要考虑燃气喷出的速度、喷管形状以及材料的选用等因素。
喷管的构造是液体火箭发动机设计中的关键环节。
根据工作原理和流动特性的不同,液体火箭发动机喷管可分为喷嘴扩张段和喷管收缩段。
喷嘴扩张段是燃气通过喷管时的膨胀区域,它起到将高压燃气的内能转换成动能的作用。
喷嘴扩张段的形状通常为锥形或抛物线形,这样可以使燃气的速度逐渐增加,提高排放速度和推力。
另一方面,喷管收缩段是将喷嘴扩张段的高速燃气加速到超音速的区域。
喷管收缩段的形状通常为圆锥形或锥形,它能够有助于提高排气速度,并减少喷嘴周围的压力损失。
在流动特性分析方面,需要考虑的主要参数有流动速度、压力分布和喷管结构对流动的影响等。
首先,流动速度是液体火箭发动机设计中最重要的参数之一。
它直接影响着液体火箭的推力和效率。
因此,在喷管设计中需要确保喷气速度足够高,同时也不能超过材料的承受极限。
其次,压力分布是流动特性分析中需要关注的另一个重要因素。
喷气速度越高,压力降低越大,喷管收缩段的压力损失也越大。
因此,需要优化喷管的形状,合理分布燃气压力,以提高喷气速度并减小压力下降。
最后,喷管结构对流动的影响是设计中需要考虑的另一个重要因素。
喷管结构的形状、材料和尺寸等都会对流动产生影响,因此需要进行相关的模拟和计算,以确保喷管的工作状态达到最佳效果。
总的来说,液体火箭发动机喷管的设计与流动特性分析是液体火箭发动机设计中的重要任务。
拉伐尔喷管的设计

拉伐尔喷管的设计摘 要:本文针对拉伐尔喷管的几何条件和力学条件进行了推导。
建立了喷管截面积变化与流速、压强、密度、温度等流动性能参数间的关系,分析了喷管出口截面下游的外界反压对拉伐尔喷管工作过程的影响。
推导建立了拉伐尔喷管主要性能参数的计算方法。
针对实际流动损失的存在,为得到喷管的实际流动性能,对理论性能参数提出了修正方法。
本文研究内容为拉伐尔喷管的设计提供依据。
关键词:变截面;力学条件;性能参数;流动损失 1.引言拉伐尔喷管是火箭发动机和航空发动机最常用的构件,由两个锥形管构成,如图1所示,其中一个为收缩管,另一个为扩张管。
拉瓦尔喷管是推力室的重要组成部分。
喷管的前半部是由大变小向中间收缩至喷管喉部。
喉部之后又由小变大向外扩张。
燃烧室中的气体受高压流入喷嘴的前半部,穿过喉部后由后半部逸出。
这一架构可使气流的速度因喷截面积的变化而变化,使气流从亚音速到音速,直至加速至超音速。
所以,人们把这种喷管叫跨音速喷管。
瑞典工程师De Laval 在1883年首先将它用于高速汽轮机,现在这种喷管广泛应用于喷气发动机和火箭发动机。
图1 拉伐尔喷管结构图2.拉伐尔喷管的几何条件 2.1变截面一维定常等熵流动在变截面一维定常流动中只考虑截面积变化这一种驱动势,忽略摩擦、传热、重力等其他驱动势,因此流动是绝热无摩擦的,即等熵流动,变截面定常等熵流动模型如图2所示。
变截面一维定常等熵流动的控制方程组为:Const m VA ρ== (1) 0dp VdV ρ+= (2)2102d h V ⎛⎫+= ⎪⎝⎭ (3)2.2截面积变化对流动特性的影响管道的形状变化可以用截面积变化dA 来表示。
(a) 截面积变化对流速的影响图2 变截面一维定常等熵流动模型对连续方程(1)取对数微分,得0d dV dAV Aρρ++= (4) 将(2)两边同除以ρ,得20dV dp d V V d ρρρ+⋅= (5) 由声速公式及马赫数定义,得()21dV dAM V A-=(6) 这就是截面积变化与流速变化之间的关系。
喷管资料讲解

c2
2
k
k
1
R
gT
*
1
p2 p*
k
喷管的出口速度取决于滞止参数与滞止压力比,而当滞止参数一定 时,流体流速仅随出口压力与滞止压力比变化。特别地:当出口压 力与滞止压力的比值为1时,出口流速为零。
当出口压力与临界压力的比值减小时,流体流速增加,当这种比 值趋于零时,出口流速趋于一最大值:
cf2m a x2kk1p*v*2kk1Rg*T*
k
k
1
R
g
T
*
1
p cr p*
k
ca,cr
kRgTcr
k 1
kRgT *
p cr p*
k
k1
k1
2kk1RgT*1ppc*rk
k
gRT*ppc*rk
整理后有:
k
pcr p*
vc
r
2
k1
k1
临界压力比与绝热指数有关,当比热容为定值时:
对于单原子气体: 对于双原子气体:
对于三原子气体:
热效应时,
J
0
,温度升高;冷效应时, 0 J
,温度降低;
零效应时,
温度不变。
0
J
k1.6,7vc r0.1.67 k1.4,vc r0.528 k1.3,vc r0.546
(三)喷管的选型原则 由于气体在喷管中的流动工程中,参数变化为:压力降低,速度增 加,因此可以根据临界压力比确定喷管的形状。
Ma 1,Ma 1,Ma 1 ppcr,ppcr,ppcr
渐缩 临界 缩放
假定气体在喷管内实现完全膨胀,即喷管出口压力 p 2 等于背压 ,
一、一维稳定流动的基本方程 (一)稳定流动质量守恒方程(连续性方程)
喷管中气体流动特性分析_郑玉

第22卷 第7期2006年7月甘肃科技Gansu S cience and Technolo gyVol.22 N o.7J uly. 2006喷管中气体流动特性分析郑 玉1,张永恒2(1兰州交通大学数理与软件工程学院,甘肃兰州730030;2.兰州交通大学机电工程学院,甘肃兰州730070)摘 要:分析各种工况下气流流经渐缩喷管和缩放喷管各个截面压力和流量的变化规律,比较两种喷管的基本特性,确定工作条件对喷管中气体流动的影响。
并介绍了实际气流曲线偏离理想气流曲线的因素。
关键词:喷管;压力;流量中图分类号:TQ 0281 引言喷管是利用气体的压力降低使气流加速的特殊管道。
它广泛应用于汽轮机等动力设备中,也应用于通风、空调等热力设备中。
例如在燃气轮机和汽轮机里都需要装设喷管,利用气体流经喷管获得高速气流推动叶轮叶片而作功。
喷管是使流体增速的变截面流道,气流在管道中流动时的状态变化情况和管道截面积的变化情况有密切关系。
由于喷管的长度较短,流速较高,气流从进入喷管到流出喷管的时间很短,因此可以认为喷管中气体流动是可逆的绝热流动过程。
喷管一般分为渐缩喷管和缩放喷管。
渐缩喷管如图1所示,缩放喷管如图2所示。
渐缩喷管沿流动方向截面积是逐渐减小的,气流作亚音速流动。
缩放喷管中气流在渐缩部分作亚音速流动,在渐放部分作超音速流动,在最小截面,即喉部是亚音速流动向超音速流动的转折点,这时流速等于当地音速。
本文通过分析空气流经渐缩喷管和缩放喷管各个截面压力变化和流量变化的规律,气流在喷管内完全膨胀,膨胀过度和膨胀不足等现象,分析研究两种喷管的特性以及工作条件的改变对喷管中气体流动的影响。
2 实验方法实验中必须测量四个变量,测压孔在喷管内不同截面位置X ,气流在该截面上的压力P ,背压P b ,(喷管出口外介质的压力),流量m 。
这些量分别用位移指针位置,可移动真空表,背压真空表以及U 型管压差计来显示。
实验是在一喷管实验台上进行,采用真空泵为动力,大气为气源。
喷管特性实验指导书

喷管特性实验一、实验目的:1、验证和加深理解喷管中气体流动的基本理论。
2、观察气流在喷管中各截面的流速,流量,压力变化规律及掌握有关测试方法。
3、熟悉不同形式喷管的机理,加深对流动的临界状态基本概念的理解。
二、实验原理:1、喷管中气体流动的基本规律气体在喷管中作一元稳定等熵流动中,压力降低,流速增加。
气流速度C ,密度ρ及压力P 的变化与截面A 的变化及马赫数Ma (速度与音速之比)的大小有关。
它们的变化规律如下表: Ma 渐缩管 0<dx dA Ma渐扩管 0>dx dA dxdc dx d ρ dx dp dx dc dx d ρ dx dp <1>0 <0 <1 <0 >0 >1 <0 >0 >1 >0 <0(1)在亚音速(Ma<1)等熵流动中,气体在0<dx dA 的管道(渐缩管)里,速度C 增加,而密度ρ,压力P 降低,在0>dxdA 的管道(渐扩管)里,速度C 减小,而密度ρ,压力P 增大。
(2)在超音速(Ma>1)等熵流动中,气体在渐缩管中,速度C 减小,而压力P ,密度ρ增大,在渐扩管中,速度C 增加,压力P ,密度ρ降低。
(3)在Ma=1,即达到临界流动状态,此时,压力为临界压力,气流速度为音速。
2、喷管中流量的计算(1)理论流量根据气体一元稳定等熵流动中,任何截面上质量流量都相等,且不随时间变化。
流量大小由连续方程、动量方程、能量方程及绝热气体方程,等熵过程方程,得到气体在喷管中流量的计算式:⎥⎥⎦⎤⎢⎢⎣⎡-⋅-==+00011221211002222)()(12γγγγγP P P P V P A V C A q m (kg/s ) 式中:0γ—绝热指数 C 2—出口速度m/s A 2—出口截面积m 2V 2—出口比体积(m 3/kg ) P 2—出口压力(MPa )P 1—进口压力(MPa ) V 1—进口比体积(m 3/kg )若:P 1=P 2时 0=m q P 2=0时 0=m q ,即在0<P 2≤P c 渐缩喷管的出口压力P 2或缩放喷管的喉部压力P th 降至临界压力时,喷管中的流量达最大值,计算式如下:1112000minmax ,)12(12V P A q k m ⋅++=-γγγ 临界压力P c 为:11000)12(P P c ⋅+=+γγγ将0γ=1.4代入P c =0.528P 1 (2)、实测流量由于气流与管内壁间的摩擦产生的边界层,减少了流动截面,因为实际流量是小于理论流量,本实验台采用孔板流量计来测量喷管的流量。
拉伐尔喷管的设计

拉伐尔喷管的设计摘 要:本文针对拉伐尔喷管的几何条件和力学条件进行了推导。
建立了喷管截面积变化与流速、压强、密度、温度等流动性能参数间的关系,分析了喷管出口截面下游的外界反压对拉伐尔喷管工作过程的影响。
推导建立了拉伐尔喷管主要性能参数的计算方法。
针对实际流动损失的存在,为得到喷管的实际流动性能,对理论性能参数提出了修正方法。
本文研究内容为拉伐尔喷管的设计提供依据。
关键词:变截面;力学条件;性能参数;流动损失 1.引言拉伐尔喷管是火箭发动机和航空发动机最常用的构件,由两个锥形管构成,如图1所示,其中一个为收缩管,另一个为扩张管。
拉瓦尔喷管是推力室的重要组成部分。
喷管的前半部是由大变小向中间收缩至喷管喉部。
喉部之后又由小变大向外扩张。
燃烧室中的气体受高压流入喷嘴的前半部,穿过喉部后由后半部逸出。
这一架构可使气流的速度因喷截面积的变化而变化,使气流从亚音速到音速,直至加速至超音速。
所以,人们把这种喷管叫跨音速喷管。
瑞典工程师De Laval 在1883年首先将它用于高速汽轮机,现在这种喷管广泛应用于喷气发动机和火箭发动机。
图1 拉伐尔喷管结构图2.拉伐尔喷管的几何条件 2.1变截面一维定常等熵流动在变截面一维定常流动中只考虑截面积变化这一种驱动势,忽略摩擦、传热、重力等其他驱动势,因此流动是绝热无摩擦的,即等熵流动,变截面定常等熵流动模型如图2所示。
变截面一维定常等熵流动的控制方程组为:Const m VA ρ== (1) 0dp VdV ρ+= (2)2102d h V ⎛⎫+= ⎪⎝⎭ (3)2.2截面积变化对流动特性的影响管道的形状变化可以用截面积变化dA 来表示。
(a) 截面积变化对流速的影响对连续方程(1)取对数微分,得控制体p +dpdxρ+d ρV +dV T +dT A +dApT A图2 变截面一维定常等熵流动模型0d dV dAV Aρρ++= (4) 将(2)两边同除以ρ,得20dV dp d V V d ρρρ+⋅= (5) 由声速公式及马赫数定义,得()21dV dAM V A-=(6) 这就是截面积变化与流速变化之间的关系。
喷管和扩压管的流动特性及工程实例PPT资料优秀版

流动特性主要取决于流速是大于还是小于当地声速。
3、马赫数
• 定义:流体中任一确定点的速度与该介质
中的当地声速的比值
• 定义式:
Ma cf c
Ma<1 亚声速流 Ma=1 声速流 Ma>1 超声速流
管内流体流速与压力及比体积的关系
• 根据热力学第一定律能量方程,等熵过程
Ma=1 dA=0 临界截面
Ma>1 Ma<1 Ma>1 dA>0 dA<0 dA>0 渐扩 渐缩渐扩
注:扩压管dcf<0,故不同音速下的形状与喷管相反
喷管和扩压管流速变化与截面变化的关系
流动状态
管道种类
喷管 dcf>0 dp<0
扩压管dp>0 dcf<0
Ma<1
1
2
dA 0 A
p1 p2
1
2 dA 0
dh cf dcf 0
渐扩喷管: cf1 >c, cf2>c
通道截面渐扩 管内流体流速度与截面积的关系 节流后做功能力下降,即节流损失。
cf dcf vdp
Ma>1时,超声速流动,dA<0 ,
(1)加速气流dcf>0,dp<0; (2)减速气流dcf<0,dp>0 ;
管内流体流速度与截面积的关系
cf1 >c, cf2≥c
通道截面渐扩
节流后做功能力下降,即节流损失。 管内流体流速与压力及比体积的关系
0 dh vdp
渐扩喷管: cf1 >c, cf2>c
水中的c=1400m/s; 空气中的c=340m/s
超声速连续变为亚声速流动,通道截面先缩后扩,在喉部达到声速。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
喷管流动特性与管道截面变化规律的关系
摘要:针对管内流动规律的一般应用中存在的问题,着重讨论了喷管内工质流动特性与管道截面变化规律的关系,从而更准确更完整地反映了喷管内工质流动规律。
关键词:喷管;流动特性;变化规律
通常在研究喷管内工质流动特性时,只着重于对喷管外形的确定,所以总是以状态参数变化为前提,去探讨工质流动截面(即管道截面)的相应变化。
这时由可逆绝热流动的基本方程组,即连续性方程、能量方程和过程方程,整理出如下两个关系式:
很明显,式(1)、(2)反映了工质流速c、压力P、截面A之间的变化关系。
从数学角度而言,这几个量是可以互为变化前提的。
但对具体的管内流动来说,究竟谁是其中的决定性因素,从而控制着(导致)其它两个量的相应变化,这自然是一个非常重要的问题。
但这一问题在很多文献[1~3]中并无明确地阐述。
显然,要揭示清楚喷管内工质的流动规律,必须揭示清楚上式中各个量的决定与被决定关系,不然问题的实质就不会充分地显现出来,所得结论也是不完整的,也就无法满足实际应用的需要。
特别是个别文献还错误地强调了这种关系,从而让人产生各种疑惑甚至是误解。
这也是许多人在学习了喷管内流动特性之后,对一些管内流动现象还仍然解释不清,甚至出现概念上的错误的根本原因。
1对喷管内流动特性与管道截面变化规律关系的分析
任何一种流动都是在一定的外部条件作用下产生的。
随流动条件的不同,管内流动现象才是多种多样的。
就喷管流动而言,其流动条件应包括如下两个方面:(一)力学条件:即喷管前后的压差;(二)几何条件:即喷管长度L和喷管流动方向(设为x方向)的截面变化规律A=f(x)。
工质降压升速、升压减速等流动特性,即工质压力P、比容v、流速c包括流动截面A的相互变化关系,应属流体自身属性,这种属性不会自发地表现出来,它是从属于流动的外部条件而存在的。
这里的力学条件是工质流动和膨胀的动力,几何条件是工质连续降压增速的保证。
在流动产生前和流动过程中,其力学条件和几何条件都是客观的,两者共同确定了相应的流动特性,缺一不可。
比如,即使在力学条件完全具备的情况下,若没有几何条件的保证,流体降压升速等属性也不会自发地表现出来。
对此还可以用一个简单的例子来加以说明:设流动的
力学条件为初压P
1与背压P
b
,在流动产生之前,只有P
1
、P
b
是客观存在的,P
1
与P
b
之间的其它压力以及其它参数都不是客观的。
只有在流动产生之后才在各
截面处对应地产生相应的压力和其它参数。
另外,即使让流体在没有管道的情况下,由初压膨胀到背压环境,也将是典型的不可逆膨胀过程,肯定不会自发地出现任何种类的管内流动特性。
并且事实上,当力学条件确定后,喷管内流动特性(即P=f(x)和c=f(x)等)将完全决定于管道截面变化规律A=f(x)。
或者说,其只随管道截面变化而变化。
这些事实清楚地表明,任何稳定流动都是在流道中产生的。
“工质做稳定而连续的流动时,其截面随参数而变化”的陈述是错误的。
真实的流动特性恰恰与此相反,即工质截面变化在先,其它参数将随截面而发生变化。
这样,我们应该把工质流动特性与管道截面变化的关系做这样的理解:在一定的力学条件下,喷管内工质压力之所以能连续降低,形成相应的增速流动特性,是由于管道截面的连续变化形成的,其实质是在力学条件下几何条件发挥的控制作用,管道是流动产生的先决条件,并且决定着流动的特性。
有何种管道就产生何种流动。
也就是管道截面的变化导致了压力的相应变化,压力的变化又导致速度的改变。
2关于管道截面A、工质压力P及流速c之间关系式的建立
本文认为,要阐述清楚喷管内工质流动的客观规律,必须阐述清楚工质流动特性与管道截面变化规律的关系,即流体自身属性与流动的外部条件的关系。
也只有这样做,所得结论才是完整的。
这可以在文中(1)、(2)式不变的情况下,对式中各量的决定与被决定关系加以明确来实现,但这样做并不能改变(1)、(2)式没有直接反映出各量的决定与被决定关系的事实。
所以也可以完全按照管内流动的真实性,顺其自然,由具体的管道入手,仍然以可逆绝热为前提,借助可逆绝热流动的基本方程组,整理出如下两个关系式:
其结论自然是,在一定的力学条件(喷管前后压差)前提下,管道截面的变化将导致工质压力的变化,而压力的改变将导致速度的变化。
由此可对三种基本管形(直管、渐缩管、渐扩管)加以定性分析,从而得出喷管流动特性(分析过程暂略)。
3(1)、(2)式和(3)、(4)式以及相应分析过程的比较
比较式(1)、(2)和式(3)、(4)后看出,它们在结构上恰是前提和结论互为倒逆的。
从数学角度来说,两组关系式本身并不存在本质上的差别,因为它们都来自同一方程组。
但必须看到,数学上能够成立的关系不一定都是真的。
显然在具体流动中,式中各个量谁决定谁是不以人的意志为转移的。
从这个意义上讲,我们是必须要重视两组关系式在内涵上的重大差异的。
显然式(3)、(4)及相应的分析过程,是在明确了流体自身属性的同时,又更加准确地揭示了各个量之间的决定与被决定关系。
它不是将流动的条件抛开,而是当作导出结论的基础,其结论相对更完整、更准确,更切合管内流动的本质,应用起来也一定是可靠的、准确的。
完整的结论必有利于问题的正确求解和深入分析。
比如,对几个常让人疑惑和误解的流动现象可分析如下:
(1) 当设计工况下的拉伐尔喷管不断地提高和降低背压时(即改变力学条件),由于流动的几何条件未变,在喷管内部,将一定程度地维持原膨胀过程不变,实现一系列的“过度膨胀”和“膨胀不足”。
(2) 同理,将设计工况下的拉伐尔喷管在出口端截短或加长,在相同的力学条件下,随几何条件的变化,膨胀过程相应被截短或加长,产生相应的膨胀不足和过度膨胀,见图1。
(3) 引射器工作原理的热力学分析:
引射器属拉伐尔喷管和超音速扩压管的联合装置。
在工程上经常被采用。
例如电厂汽轮机的主要辅助设备之一抽气器就是如此。
其热力学原理可分析如下
(见图2):工作蒸汽压力为P
1,大气压力为P
b
,A
c
是最小截面,A
max
是最大截面,
A b 是L
1
段拉伐尔喷管在设计工况(相应于大气)下的出口截面。
其流动特性是在
P b <P
c
(临界压力)的力学条件下,L
1
段呈相应的拉伐尔喷管特性,产生超音速气
流;接着在不断扩展的几何条件作用下,气流在L
2段产生过度膨胀,在A
max
处产
生最低的压力P
min ,形成高度真空,L
3
段渐缩管中,超音速气流扩压后流入大气。
实际上,在管内流动特性的具体应用中,经常会遇到与上述各实例相类似的一些问题。
其主要特点是,从整体上已不属于管内基本流动现象。
这时就要求人们能运用所掌握的基本流动特性,对其加以正确的理论分析和求解才行。
不然,就无法满足实际问题的需要。
显然由(3)、(4)式提供的分析方法和过程能够很好
地保证这一点。
4热力学“可逆”条件下所得结论的局限性
对于喷管来说,在可逆前提下,其外形和出口(缩放喷管同时包括喉部)截面确定了唯一的可逆热力过程和工质流量,而喷管的长度以及处于进出口和喉部截面之间的所有其它截面大小,在一定外形前提下都是任意的,不论其大小如何,其所对应的理想热力过程及流量都是唯一的。
当然,对于实际流动来说,不同的喷管长度及不同的截面变化规律,其不可逆性是明显不等的。
但究竟哪一种情形属不可逆性相对较小,怎样才能设计出一只更加理想的喷管,这些实际结论的探讨已超出了热力学的研究范畴,这是需要由其它相关学科(如气体动力学)来回答的问题。
5结论
以上,本文着重分析和强调了喷管内流动中量与量的决定与被决定关系,从而提出了与以往不同的分析方法,重新建立了更切实际的流道截面、工质压力与流速之间的关系式。
界定了热力学课程和其它课程的分工与联系。
当然,喷管内实际流动过程是比较复杂的,许多实际结论的得出只靠热力学是远远不够的。
但不论如何,热力学可逆前提下的理想结论还是要由热力学本身来回答的,是其它学科无法替代的。
热力学应使本身所表达的内容更准确、更完整、更深入。
参考文献
[1]庞鹿鸣,汪孟乐,冯海仙.工程热力学(第二版)[M].高等教育出版社,1986.
[2]沈维道,郑佩芝,蒋淡安.工程热力学(第二版)[M].高等教育出版社,1983.[3]孔珑.工程流体力学(第二版)[M].中国电力出版社,1992.。