气液针栓喷注单元雾化角模型

为了准确预测气液针栓喷注单元的雾化角,基于动量守恒建立了液束撞击气膜的雾化角理论模型,通过试验结果获得变形因子,并分析了结构参数和工况参数对雾化角的影响规律。结果表明:局部动量比对雾化角影响最大,其他结构参数和工况参数都是通过影响局部动量比而进一步决定雾化角;随着局部动量比增大,液束变形程度减小,液束变形导致有效撞击动量比小于几何动量比。根据试验结果分段给出了变形因子,当局部动量比范围为0～3时,变形因子推荐值为0.61,当局部动量比范围为3～4.5时,变形因子推荐值为0.70,当局部动量比范围为4.5～7.2时,变形因子推荐值为0.75,引入变形因子的理论预测值与试验结果吻合很好,该模型可为气液针栓喷注器的理论研究和工程设计提供参考。      

针栓式喷注器液膜下漏率预估模型

为了对针栓式喷注器液膜下漏率进行准确预测,基于针栓式喷注单元喷雾场结构分析,结合理论推导、数值仿真及试验研究3种方法建立了液膜液束各自变形的相对变形量模型;在考虑液膜液束变形的基础上,引入相互影响系数表征多喷注单元间相互影响,建立了实际阻塞率和实际下漏率模型。通过数值仿真及试验结果的多参数充分验证,结果表明:理论预估模型与数值仿真及试验结果一致性较好。液膜液束相互作用下,液膜绕液束流动和液束根部横截面前后缘位置移动不同步导致的展向变宽分别是液膜和液束发生变形的主因,且有效动量比越大,液膜相对变形越大,液束相对变形越小。对于一定阻塞率的几何结构,结果表明:下漏率随着有效动量比的增大而增大,增大趋势呈先快后缓,且实际下漏率均小于几何下漏率,这是由膜束变形导致的实际阻塞率比几何阻塞率更大造成的。另外,发现液膜下漏率仅与表征流场结构(有效动量比)及几何结构的无量纲参数(液膜厚度与液束直径之比和阻塞率)有关,与喷射速度的绝对值无关,并给出了模型中的常系数供工程设计预估使用,对从设计初期就考虑针栓头的热防护问题具有重要的指导意义。     

针栓式喷注单元膜束撞击雾化混合过程数值模拟

为了全面认识针栓式喷注器喷雾场结构，基于自适应网格加密技术和分三相计算的PLIC VOF（Piecewise Linear Interface Calculation Volume of Fluid）方法对针栓式喷注单元膜束撞击雾化混合过程进行了仿真分析，通过对两路推进剂分别进行界面追踪，获得了膜束撞击雾化混合过程的详细结构特征，与高速摄影试验结果定性定量对比均吻合较好，验证了数值方法的准确性。以此为基础对膜束撞击的喷雾场结构、撞击变形过程、流场涡结构、雾化破碎典型特征及破碎后的雾化混合分布特征进行了识别分析，结果表明：膜束撞击形成了液束未穿透液膜和液束穿透液膜2种不同的喷雾扇结构。膜束撞击形成的喷雾扇呈"Ω"形，膜束同时发生弯曲变形和横截面变形。另外，膜束撞击同时受到正压和剪切应力作用，导致了一系列复杂涡流现象，使得相互作用增强，雾化混合均增强，这也是膜束撞击喷注构型优于膜膜撞击的本质原因。最后，还发现膜束撞击喷雾场液滴分布呈现分区结构特征，分别是液束控制主导的上雾化区、液膜控制主导的下雾化区及夹在中间的混合区，实际中应兼顾雾化特性和混合特性，选取中等动量比膜束撞击，这可为针栓式喷注器的理论研究和工程设计提供重要参考。     

壁面边界对撞击合成动量角的影响研究

为了研究动量比对针栓式喷注器撞击雾化合成动量角的影响规律,基于三相的与水平集耦合流体容积法(CLSVOF)对贴壁液膜/自由液膜撞击的合成动量角进行了数值模拟,并与自由液膜/自由液膜撞击的合成动量角进行了对比,结合试验结果及理论预估结果详细考察了有无壁面条件下合成动量角与动量比的关系,深入分析了造成两者之间差异的根源,揭示了壁面边界对合成动量角产生影响的作用机制。结果表明:CLSVOF方法计算的合成动量角与试验结果一致,最大相对误差约为10%,大多数工况点的相对误差小于5%;有壁面边界的贴壁液膜/自由液膜撞击合成动量角显著大于无壁面边界的自由液膜/自由液膜撞击的;仅一路流体贴壁的撞击合成动量角与常用的入口动量比理论预测值最大相差20°以上,而两路均无贴壁或者均贴壁的撞击合成动量角与理论预测值吻合很好。这一显著差异的根源在于两者撞击形成的高压区分布显著不同。壁面的存在使得撞击点附近形成的高压区对两路流体作用不对称,有壁面边界时壁面承受高压迫使其对贴壁流体有强的作用力,导致垂直于壁面方向动量不守恒,根据入口动量比预测的理论不再适用。     

燃烧条件下凝胶自燃推进剂雾化特性试验研究

为了研究凝胶自燃推进剂的雾化特性及敏感因素,在单互击式喷嘴矩形燃烧室内进行了凝胶一甲基肼/四氧化二氮(MMH/NTO)喷雾燃烧过程可视化试验,采用光源后置消光法湮灭火焰自然辐射发光,采用彩色高速摄影获取了燃烧条件下的高质量雾场阴影图像,通过图像处理,有效提取了雾场的雾化锥角、破碎长度、液丝直径及液丝运动速度,分析了撞击角、射流速度和动量比的影响。结果表明,凝胶MMH/NTO稳态燃烧时可观察到液膜、贯穿视场的液丝和红棕色NO_2气体,推进剂混合燃烧不充分;撞击角从75°增大到105°,凝胶MMH/NTO撞击后的破碎长度、液丝直径减小,视场内可视红棕色NO_2气体变少,撞击角为105°时,推进剂会附着在喷注面上,从而影响液膜横向铺展,雾化锥角反而最小,建议撞击角选取90°。燃料射流速度从23m/s增大到45m/s,凝胶MMH/NTO撞击后的雾化锥角及液丝运动速度增大,破碎长度及液丝直径减小,雾化模式发生改变。动量比从1.04增大到1.52,凝胶MMH/NTO撞击后的雾化锥角及液丝运动速度增大,视场内红棕色NO_2气体变少。故一定量程内增加撞击角、射流速度、动量比有助于凝胶MMH/NTO推进剂混合燃烧。     

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

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

射流撞击雾化索太尔平均直径的计算

针对射流撞击雾化初始形成液膜的特征，推导出极坐标下自变量为极角的液滴尺寸分布函数的解析式，这种形式的分布函数与经典的Ｒ～Ｒ分布、上限对数正态分布等概率分布函数不同，可以计算任意撞击角下两股相同射流撞击或一股射流撞击壁面边缘在任意极角区域雾化所产生液滴的索太尔平均直径。计算结果表明，索太尔平均直径随撞击角增加而单调减小；射流直径、射流速度、射流温度变化时，索太尔平均直径也呈一定单调规律变化。本文还结合新一代液体火箭发动机层板式喷注器喷注单元雾化的特征进行了有关雾化机理的探讨。     

凝胶推进剂雾化特性研究进展

为了进一步认识凝胶推进剂雾化过程,促进凝胶发动机的设计和优化,综述了射流撞击式、离心式、气泡式以及燃烧条件下凝胶推进剂雾化特性的研究进展。综述结果表明,凝胶推进剂雾化性能明显差于牛顿流体推进剂,凝胶液膜尺寸显著增大,液丝较难破碎为小粒径液滴;射流撞击式喷嘴对凝胶推进剂的雾化效果优于离心式喷嘴;随着射流雷诺数和韦伯数增大,撞击式凝胶液膜的雾化模式依次为边缘闭合模式、边缘开放模式、无边缘射线模式、液丝分离模式和充分发展模式;采用锥形结构、减小喷嘴出口长径比、方形和椭圆形喷嘴出口皆有利于凝胶液膜破裂,且增大喷注压力、撞击角、温度、室压和减小撞击距离均能改善液膜雾化效果。燃烧条件下MMH/NTO凝胶推进剂撞击液膜破裂雾化机制在宏观上与冷模条件下凝胶推进剂模拟液撞击液膜较好地吻合。此外,对凝胶推进剂雾化特性的进一步研究工作提出了建议。     

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

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

互击式气-气喷注器混合场影响因素流动显示

虽然针对互击式气-液喷注器开展了大量研究,但对互击式气-气喷注器的了解甚少。为了深入理解互击式气-气喷注器气气混合过程的特点,利用粒子图像测速(PIV)系统对其开展了流动显示实验研究。实验结果表明:对于O-F-O互击式气-气喷注器,气流撞击角度和氧/氢动量比直接影响着气流撞击作用的大小,控制着气流混合过程;对于喷孔间距保持一定时,撞击角越大,气流间的相互作用越强,气流混和效果越好;氧/氢动量比越大,气流间撞击力越强,混合距离越短,对中心气流的挤压作用越厉害,其流场分布三维特征越明显;可以通过组合高的氧压降、低的氢压降和大撞击角来设计混合效果较好的互击式气-气喷注器。     

动量比对两股互击式喷嘴雾化特征的影响

基于CLSVOF（coupled level-set and volume of fluid）方法对两股互击式喷嘴雾化过程进行了数值模拟,仿真结果与实验结果进行对比,并详细考察了动量比对液膜特性的影响.结果表明:①CLSVOF方法能对射流撞击形成的液膜形态进行有效捕捉,计算获得的液膜形态与实验结果一致.②两股不同孔径的射流在撞击点撞击过程中,由于大孔径一侧的部分射流未参与撞击,而中心位置的液体在撞击之后速度迅速下降,周围的液体则保持撞击前的速度,使得射流撞击之后形成弧形状液膜.③在不同动量比条件下,射流撞击后轴线位置的合成角主要取决于参与撞击射流的动量比,数值计算结果与理论计算结果吻合较好,误差均小于10%.      

倾斜射流撞壁形成的液膜外形的理论建模

倾斜射流撞击壁面在燃烧室液膜冷却、溅板式喷注器雾化等领域均有广泛应用.为了研究倾斜射流撞壁形成的液膜的基本形态、液膜边界大小,开展理论建模研究.通过建立在液膜边界的守恒方程及液膜在壁面上的厚度及速度分布关联式,形成一套能够求解倾斜射流撞击壁面后液膜形态及边界的理论方法.模型计算结果与文献中的实验结果对比表明,建立的模型...     

幂律流体撞击射流破碎机理的实验

通过自行研制的撞击射流系统,采用高速摄像技术,对4种不同卡波姆质量分数剪切变稀的卡波姆凝胶非牛顿幂律流体在不同射流条件下撞击液膜的破碎模式、喷雾锥角、破碎长度进行实验与分析.结果表明:不同流变特性的幂律流体出现的破碎模式不同,黏度越高,破碎模式越复杂;4种卡波姆凝胶在韦伯数为8000左右时,达到完全发展模式,并在韦伯数为8000左右时表现出相似的破碎特性;随着韦伯数的增大,喷雾锥角逐渐增大并趋于稳定,稳定值小于撞击角度,当卡波姆质量分数增大时,最大喷雾锥角逐渐减小,最终收敛于60°;随着韦伯数的增大,4种卡波姆凝胶的破碎长度先小幅变大然后减小,理论结果与实验结果整体趋势相似.      

液气式针栓喷注单元喷雾场结构试验研究

径向孔型针栓喷注器相对于径向缝型针栓喷注器具有更复杂的喷雾场。为了研究径向孔型针栓喷注器的喷雾场结构，将径向孔型针栓喷注器简化为单个气体射流与液膜碰撞的针栓喷注单元，采用了背景光成像系统结合激光相位多普勒技术（PDA），以水和空气为模拟介质，对液气式针栓喷注单元的喷雾场进行了试验研究。试验结果表明，液气式针栓喷注单元喷雾的三维结构呈现“喇叭”状。根据喷雾的形成过程及液滴的分布，液气式针栓喷注单元喷雾可以划分为4个区域：碰撞区、液滴区、液雾区及液丝区。液气式针栓喷注单元喷雾的分布范围可由内边界角、外边界角、中线角及散布角表示，均随局部动量比的增大而增大。液滴区的粒子主要由碰撞过程产生，SMD较大；液雾区的粒子经碰撞过程产生后，在气动力作用下进一步雾化，SMD小。由于液雾区的速度和粒径同时受到气动力作用的影响，粒径分布与速度分布在空间上呈现负相关趋势。     

低Weber数射流撞击雾化的数学模型

通过分析低Weber数无粘射流撞击形成液膜雾化的特征, 由质量守恒和动量守恒, 并借鉴前人的理论和实验结果, 推导出Taylor心形波作用下计算液膜厚度和形状的近似解析式, 数值计算结果与文献所给实验数据基本吻合一致。结合矩形喷孔形成的扇形液膜的破碎理论, 进一步计算了不同撞击角下雾化形成液滴尺寸的大小及分布。