直流互击式喷注单元雾化特性准直接数值模拟

为实现直流互击式喷注单元雾化过程的数值求解并探究结构及工作参数对雾化特性的影响规律,基于一种树形自适应加密算法与流体容积(Volume-of-Fluid,VOF)方法实现了雾化过程的准直接数值模拟。计算得到了两股射流由喷射到撞击形成液膜,液膜进一步破碎形成液丝、液滴的全过程,获得了液膜的破碎长度、液滴的Sauter平均直径、雾化频率等雾化特性参数。通过将典型算例的计算结果与试验数据进行对比验证了计算的有效性,给出了数值求解精度。探讨了撞击波的形成机理,分析了雾场液滴的尺寸分布规律,撞击夹角、孔径比、射流速度、动量比对雾化特性的影响规律。结果表明:所采用的算法可以实现多相、多尺度雾化过程的数值求解;撞击波的形成是由于两股射流撞击时惯性力的不完全对称导致的;雾场液滴的尺寸分布近似服从Rosin-Rammler分布;撞击夹角的增大与射流速度的提高导致液膜的破碎长度减小,液滴的平均粒径减小,撞击夹角增大雾化频率呈减小的趋势,射流速度提高雾化频率呈增大的趋势;动量比为1而孔径比不为1时会形成凹形液膜,雾场存在一定程度的偏斜;动量比主要影响雾场的偏斜程度。     

背压对撞击式喷嘴雾化特性影响研究

为研究背压对撞击式喷嘴雾化特性的影响,将压力的变化等效为气体密度的变化,基于一种树形自适应加密算法,通过直接数值求解不可压Navier-Stokes方程组实现了不同背压条件下射流撞击雾化的数值模拟。首先将数值模拟结果与试验数据进行对比,验证数值模拟的有效性,在此基础上开展了高背压条件下雾化过程的数值模拟。结果表明,随着背压的提高,气动力相应增强,液膜的破碎更加剧烈,一次雾化区域的液滴数密度增大,雾场由稀疏向稠密发展;液膜在向下游运动过程中波动速度的幅值逐渐增大,并且随着背压的提高,液膜的波动由线性向非线性转变;背压增大导致破碎长度减小,经过参数修正得到了液膜破碎长度的经验公式,并与试验数据进行了对比;背压对液滴尺寸分布规律没有显著影响,但随着背压的提高,同一时刻大液滴所占的比例提高,整个雾场的Sauter平均直径有增加的趋势,当背压从0.1MPa增大到1MPa时,雾场的Sauter平均直径由155.5μm增大到166.9μm;背压增大,液滴粒径分布的均匀度指数减小,液滴尺寸分布更加不均匀。     

液体射流撞击液膜振荡行为的实验研究

液体射流撞击过程是液体火箭发动机常见的高效雾化方式。为了进一步探究液体射流撞击雾化过程,采用高速摄像仪和图像处理技术对两股液体射流撞击液膜的振荡行为进行了实验研究,考察了射流韦伯数(21≤We≤1356)和撞击角(2θ=60°,90°,120°)对液膜振荡的影响。结果表明,随着射流韦伯数增大(We250),撞击液膜会从竖直稳定模式转变为振荡模式,促进液膜破裂和雾化。增大撞击角能加剧液膜振荡和增加雾化角。2θ=60°和90°的液膜振荡区间为250We400,2θ=120°的液膜振荡区间为We250,无量纲振幅的时均极大值约为5。结果揭示了随着液体射流速度的增大,射流不稳定性的发展会引起撞击点的动量不平衡,进而形成液膜振荡。此外,液膜向下传播时与空气界面间的Kelvin-Helmholtz不稳定性也促进液膜振荡。     

两股湍流射流撞击雾化过程的数值研究

采用新提出的雾化过程仿真算法,对两股湍流射流撞击雾化过程进行了数值研究,重点关注了撞击角和动量比的影响。为精细地模拟湍流对雾化的影响,采用了大涡模拟。对雾场瞬态和统计结果的分析表明,喷雾角随着两股射流撞击角的减小而减小,但液膜的破碎长度反而增加;动量比主要影响雾场的偏斜程度,动量比由1.0增大到2.25时,喷雾偏斜角由0°增大到约23°。仿真结果与试验结果吻合较好,表明新算法能够较好地反映两股湍流射流的雾化过程。     

撞击角对撞击式喷嘴雾化特性影响研究

为了研究液体火箭发动机中撞击式喷嘴的雾化问题,采用CLSVOF(Coupled Level-Set and Volume-of-Fluid Method)方法对撞击式喷嘴的雾化过程进行了数值仿真模拟,重点考察了撞击角对撞击式喷嘴雾化特性的影响。结果表明,通过CLSVOF方法能对射流撞击雾化形态进行较好的捕捉;分析整个雾化区域的液相分布,将射流撞击雾化过程中的液相分布大致分为三个区域:撞击前两股独立射流区域、撞击后形成的液膜区域、液膜破碎后液丝和液滴的生成区域;射流撞击雾化除了气液界面上速度差引起的不稳定之外,还存在由射流湍流或撞击波引起的其它不稳定因素;液膜破碎长度随撞击角的增大而减小,液膜表面上的表面波幅值及液膜破碎程度则随撞击角的增大而增大;射流撞击在撞击点位置处所形成的速度差对液膜的破碎和液膜上表面波幅值的大小起到了关键作用。     

凝胶推进剂模拟液撞击雾化特性研究

为了寻求凝胶推进剂雾化机理、提高其雾化效果,开展幂律流体双股圆柱射流撞击的理论和实验研究。首先推导了幂律型流体射流撞击形成液膜的理论模型,在其中引入粘性力和能量损失。进而通过设计凝胶模拟液射流撞击实验,使用高速摄像机拍摄图像的方法进行了相关的实验研究,以验证理论的正确性。在此基础上,对撞击角度、射流速度分布及液体部分物性参数对液膜形状、厚度及速度分布的影响进行了分析,得到了相关参数的影响规律。研究结果表明,射流撞击角、速度型、液体表面张力系数、流变特性等均对液膜特性有明显影响,且流体本身的物理性质对撞击形成的液膜的影响更甚。     

基于SIFT算法的含碳颗粒凝胶推进剂雾化场速度分析

雾化是凝胶推进技术的关键问题之一,雾化过程中的液膜、液丝及液滴等的速度分布、液膜厚度等参数对于判断雾化效果、揭示雾化机理具有重要作用。为了对含碳颗粒凝胶推进剂雾化场的速度进行定量分析,提出了一种基于SIFT关键点匹配的雾化场速度计算方法,并以雾化场速度分析为基础,提出了一种新的液膜厚度估算方法。结果表明:雾化场速度随射流速度的增大而增大,随撞击角度的增大而减小;雾化场平均速度与射流速度的比值vato/vjet在0.6~0.9,可用于表征雾化效果,vato/vjet越小,雾化效果越好;射流撞击形成液膜的厚度在0.04~0.13mm,液膜厚度随着射流速度及撞击角度的增大而减小;凝胶推进剂的雾化效果随着碳颗粒浓度的增大而降低,随着碳颗粒粒径的增大而略有改善。     

相邻离心式喷嘴液膜撞击雾化过程仿真

基于Gerris软件建立的锥形液膜雾化破碎过程数值仿真方法,对相邻多个离心式喷嘴液膜撞击雾化过程进行了数值仿真,可视化展示了喷雾场三维形态和结构特征,并做了流场分析,讨论了喷嘴之间的相互干扰作用,获得了液滴空间分布,对相邻不同数目的喷嘴雾化效果进行了比较。结果表明:双喷嘴液膜撞击会形成一个类似互击式射流撞击形成的扇形区域,但又与其不同;锥形液膜的撞击效果主要依赖径向速度而不是切向速度,液滴平均粒径与撞击液膜间的旋向几乎无关。多喷嘴液膜撞击会使得液滴空间分布发生大的变化,液膜破碎长度会缩短。多喷嘴液膜撞击后的雾化特性是否改善与单喷嘴原先的雾化特性以及喷嘴间的距离密切相关。对于文中特定的喷嘴结构、喷嘴间距及排列方式,双喷嘴的液滴SMD均比单喷嘴的增大约4.8%~6.1%;四喷嘴的液滴SMD比单喷嘴的增大约7.3%,比双喷嘴的增大约2.4%。     

高压环境下剪切稀化非牛顿撞击射流直接数值模拟

高压环境下的撞击射流是液体火箭推进系统中广泛采用的一种燃料雾化方法，其雾化效果将直接决定最终燃烧效率。采用直接数值模拟（DNS）工具，对10 MPa高压环境下剪切稀化非牛顿直角撞击射流的三维非定常雾化特征、机理及非牛顿特性进行了研究。结果表明：高压环境下该撞击射流的雾化流场呈圆形辐射状分布并形成Mushroom头部和Ω状的局部凸起，气体中的涡量分布表现为有序贴附区和无序爆炸区两类，液膜向液丝的破碎主要受平均气体力和平均黏性力作用影响，而液丝向液滴的破碎则主要受局部流场参数影响，撞击雾化过程中液体无量纲表面积不断增长并可分为5个阶段。此外，撞击射流头部在局部强剪切力作用下其无量纲黏性系数最低降至仅0.7。     

幂律型流体雾化SPH方法数值分析

幂律型流体的雾化过程存在复杂的界面运动,用传统网格法很难精确追踪运动界面.为研究幂律型流体的雾化特性,运用SPH方法对典型的双股幂律型流体撞击雾化问题进行三维数值分析.根据文献实验,采用SPH方法模拟获得与实验条件相应的数值结果,对比后表明,二者雾化角相当吻合,数值结果还成功捕捉到液膜向液丝的破碎过程及网状的液丝分布状态,验证了方法的有效性.分析了射流速度和撞击角度对雾化角的影响,得到雾化角随着射流速度和撞击角度的增加而增大;对雾化后的速度场进行数值分析后表明,撞击点附近,惯性力的作用使速度场变化剧烈;在流体远离撞击点的过程中,粘性耗散作用使速度场趋于稳定,但速度大小小于初始撞击速度.     

喷油杆和凹腔支板稳定器近距匹配的液雾分布可视化

液雾分布与稳定器的燃油布置方式、油气分配、火焰稳定及火焰传播密切相关。以RP-3为雾化介质,在来流马赫数为0.2及来流温度为10~400℃的条件下,采用高速摄影法和激光片光/照相法,对喷油杆与凹腔支板稳定器间隔31.5 mm且顺喷时的液雾分布特点进行了可视化研究,并探讨了来流温度及油气动量比对其液雾分布轨迹的影响。结果表明:顺喷喷油杆与凹腔支板稳定器近距匹配方式有利于燃油在支板前缘形成挡溅雾化;挡溅雾化后的一部分燃油在支板表面形成油膜,并在凹腔前缘与尾缘进行二次雾化,来流温度较高时高温支板表面也有利于燃油的蒸发雾化;另一部分燃油则以类似横向射流的形式进行雾化。当来流温度一定时,油气动量比增大,液雾轴向分布距离和横向穿透深度均增大;来流温度升高,液雾穿透深度增加,油气动量比对液雾分布的影响更明显。      

通道高度对受限长通道流动特性的影响

为更好地掌握涡轮叶片错排射流冲击受限长通道内的流动结构,了解横流对射流作用的影响,获取沿流动方向通道内的流动特性,对包含12个射流孔和1个出流孔的受限通道流场进行了详细的测量,着重研究受限通道和出流孔内的流动特性,同时获得沿通道密流比、侧壁静压以及孔流量系数的变化规律。实验结果表明:通道高度主导着通道和出流孔内的流动特性,但其影响随着通道高度的进一步增加而减弱;通道高度的增加削弱了横流强度,但也增加了射流流程,使其在下游位置无法对靶面形成有效冲击;随通道高度的增加,沿通道流量分配趋于均匀,密流比呈线性分布,沿程压降趋缓,出流孔流量系数提高。     

幂律型流体撞击射流破碎特性直接数值模拟

凝胶推进剂的双股射流撞击雾化广泛应用于液体火箭发动机的燃烧室中,其破碎特征及雾化效果直接影响燃烧效率。为探究雾化特性的发展规律,采用直接数值模拟DNS方法,对射流速度为100m/s的剪切稀化非牛顿液体正交撞击产生的雾化特征、液体表面积、表面波、涡特性以及非牛顿特性开展研究。结果表明,射流下形成的雾化流场迅速扩张形成液膜,液膜两侧边缘破碎成大量的液丝与液滴,核心部分产生撞击波后在气体力的作用下逐步发展为带有凸起和褶皱的不稳定表面波,其撞击波波长最大可达2.46倍射流直径。液体表面积不断增长,但无量纲表面积总体呈现先下降再上升的趋势。气体中的涡量分布则分为有序附着区和无序爆炸区两类,并且涡量主要集中分布于气相区域。此外,射流撞击时产生强剪切使该液体内部的粘性系数下降,最低仅为初始粘性系数的0.3倍。     

Modelling of the atomization of a plain liquid fuel jet in crossflow at gas turbine conditions

The penetration and atomization of a plain jet of kerosene fuel in air crossflow were modelled and compared with the experiment at test conditions relevant to lean premixed prevaporized combustion in gas turbines.Tests were conducted at an air velocity of 100 and 75 m/s and an air pressure of 6 and 9 bar at ambient temperature. The fuel nozzle diameter was 0.45 mm and the liquid-to-air momentum flux ratio was between 2 and 18. Measurement techniques employed include time-resolved shadowgraphs, Mie-scattering laser light sheets, and Phase-Doppler-Anemometry. Hence qualitative information on jet disintegration as well as quantitative data regarding jet penetration, fuel placement and dropsize distribution were obtained.The numerical model predicts the jet penetration by treating the liquid fuel jet conceptually as a cylinder deflected by the airflow. Fracture of the liquid column is modelled by a characteristic time criterion. Shear breakup of the liquid column and its fragments is modelled as liquid boundary layer stripping. The engineering model has been implemented as an extension into a 3d CFD-code with Lagrangian particle tracking. Within its application limits, the model fills the gap of lacking droplet initial conditions in numerical spray simulation, while at the same time requiring low computational effort.     

Experimental investigation and correlation development of jet impingement heat transfer with two rows of aligned jet holes on an internal surface of a wing leading edge

Extensive experimental studies on the heat transfer characteristics of two rows of aligned jet holes impinging on a concave surface in a wing leading edge were conducted, where50000 Rej 90000, 1.74 H/d 27.5, 66° a 90°, and 13.2 r/d 42.03. The finding was that the heat transfer performance at the jet-impingement stagnation point with two rows of aligned jet holes was the same as that with a single row of jet holes or the middle row of three-row configurations when the circumferential angle of the two jet holes was larger than 30°. The attenuation coefficient distribution of the jet impingement heat transfer in the chordwise direction was so complicated that two zones were divided for a better analysis. It indicated that: the attenuation coefficient curve in the jet impingement zone exhibited an approximate upside-down bell shape with double peaks and a single valley; the attenuation coefficient curve in the non-jet impingement zone was like a half-bell shape, which was similar to that with three rows of aligned jet holes; the factors,including Rej, H/d and r/d, affected the attenuation coefficient value at the valley significantly.When r/d was increased from 30.75 to 42.03, the attenuation rates of attenuation coefficient increased only by 1.8%. Consequently, experimental data-based correlation equations of the Nusselt number for the heat transfer at the jet-impingement stagnation point and the distributionof the attenuation coefficient in the chordwise direction were acquired, which play an important role in designing the wing leading edge anti-icing system with two rows of aligned jet holes.     

不同淬熄结构下RQL燃烧室流场特性的数值研究

采用雷诺平均Navier-Stokes(RANS)方法对不同淬熄结构的富油/淬熄/贫油(RQL)燃烧室内定常掺混流动过程进行了数值模拟,分析了不同截面的流场结构,研究了动量通量比及淬熄孔排布方式的影响。结果表明:射流进入RQL燃烧室后发生偏转且偏转方向和旋流方向有关,回流区长度、穿透深度和高湍动能区随着动量通量比的增大而增加;当J≥80时,淬熄射流下游出现高湍动能区,逐渐与上游融为一体;当J≥120时,回流区的长度基本不再发生明显变化。此外,非中心对冲结构C3具有更大的高湍动能区,掺混更均匀。      

煤油在超声速气流中非定常横向喷射的实验观察

本文目的是研究煤油在超声速气流中非定常横向喷射问题,并与静止和亚声速气流中的结果进行对比。改造普通油嘴并采用单次高压燃油喷射系统,得到不同延时的流场阴影照片。结果表明：在超声速气流中,射流柱破碎是由迎风面的表面波引起的,破碎点位于表面波的波谷。雾化过程过分四个区域。射流激波由近射流柱的迎风面,并在人壁反射,非定常喷射的波系演化复杂。煤油在亚声速气流中喷射,不存在射流激波,破碎点沿射流柱上移,穿透深     

有/无凹槽通道内煤油超燃雾化的测量研究

为研究超燃流场中雾化煤油液滴直径分布,针对有/无凹槽(L/D=3)直通道,对煤油在超声速气流中的雾化过程进行了测量。利用纹影法得到雾化流场波系结构,利用平面激光诱导荧光(PLIF)得到煤油扩散范围,基于激光前向散射原理的激光粒度仪,得到液滴直径分布。结果表明,当来流和喷射条件相同,平板通道内的煤油射流穿透深度略大于凹槽通道。液滴直径呈在截面中心小、外侧大趋势,沿外侧的液滴直径波动较中心区大。煤油距喷口30 mm处已完成雾化,雾化区的液滴Sauter直径为4~8μm,凹槽对煤油穿透深度和雾化尺寸无明显贡献。     

冠齿喷嘴射流冲击平直靶面对流换热实验

利用红外热像仪测试了冠齿喷嘴射流冲击平直靶面的对流换热特性,在射流雷诺数为5 000~20 000、冲击间距比为1~8范围内,与普通圆管射流进行了对比,并对冠齿数和冠齿长度的影响进行了初步分析。研究结果表明,冠齿射流冲击对流换热显著高于圆形射流,在小冲击间距下,冠齿射流冲击的局部对流换热系数分布在冲击驻点附近呈现明显的梅花瓣状特征,当射流冲击间距比达到4以后,冠齿射流的局部对流换热系数分布则呈现出常规圆形射流冲击的特征;以2倍或4倍射流直径作为区域平均范围,冠齿射流的区域平均努塞尔数相对圆形射流的增加幅度在15%~30%之间,相对增加幅度与射流雷诺数和射流冲击间距比相关;在本文的冠齿结构参数范围内,冠齿伸出长径比为0.6的6-冠齿结构取得的射流冲击强化传热效果较优。     

某航空齿轮润滑喷嘴结构对喷嘴射流特性的影响研究

为了研究喷嘴结构对喷嘴射流扩散角和射流流量的影响,采用基于FLUENT的大涡数值模拟方法对喷嘴三维射流流场进行了二相流仿真和单相流仿真,分别通过高速可视化和喷射滑油称重的方法测量了喷嘴射流扩散角和射流流量,并用测量结果对二相流仿真和单相流仿真进行了合理性验证.当喷嘴长径比大于2且喷孔轴线与喷孔上游管道轴线的夹角在0°~90°之间时沿喷孔轴线距喷孔出口100mm范围内射流扩散角近似为零,且增加喷孔直径、喷孔轴线与喷孔上游管道轴线的夹角或者减小喷孔轴向长度均会引起射流流量的增加.结果表明:当喷嘴长径比大于2且喷孔轴线与喷孔上游管道轴线的夹角在0°~90°的范围内变化时对射流扩散角的影响很小;喷孔直径、喷孔轴向长度和喷孔轴线与上游管道轴线的夹角通过影响局部损失和沿程损失,从而影响喷嘴射流流量.