基于Gerris的离心式喷嘴锥形液膜破碎过程数值模拟

为了实现离心式喷嘴锥形液膜破碎过程的精细可视化和雾化特性的准确计算,基于自适应网格加密技术和VOF方法,建立了一种数值方法,可以捕捉到包含一次雾化和二次雾化的喷雾场结构特征,获得了全场液滴的空间分布和粒径分布,液滴捕捉效果逼真。针对给定结构尺寸的敞口型离心式喷嘴,计算了在不同流量和不同切向孔直径下的液膜锥角、液滴平均粒径SMD的变化,并与高速摄影拍摄的锥形液膜破碎过程图像和PDPA(Phase Doppler Particle Analyzer)测量的液滴SMD对比。结果表明,数值模拟与实验所获得的喷雾场吻合,获得的液膜锥角和液滴平均粒径最大相对误差分别为4.18%和11.82%,结果吻合较好,验证了数值方法的准确性,为离心式喷嘴的精细研究和设计应用提供一种新手段。研究表明,切向孔直径对液膜锥角和液滴平均粒径的影响较显著,在设计加工时,是一个比较重要的结构参数。     

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

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

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

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

内直外旋气液同轴式喷嘴流量及雾化特性

针对气体中心直流、液体涡流器离心式内缩进气液同轴喷嘴,采用高速摄影和PDPA测量了其流量和雾化特性。发现在常温常压下,喷嘴流量近似与喷注压降的0.477次方成正比。喷嘴喷雾场为空心锥形,雾化角约为110°,随液体压降增加,雾化锥角增大。液体压降不变,雾化锥角受中心气体卷吸作用随气体压降增加而减小。液膜破碎距离与雾化锥角有关,雾化锥角越大,液膜破碎距离越短,反之亦然。随径向距离增加,雾化液滴轴向速度、径向速度减小,SMD先增大然后略有减小,沿径向呈"马鞍"状双峰分布;随轴向距离增加,粒子轴向速度减小,SMD减小。     

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

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

喷嘴形式对幂律型非牛顿推进剂雾化特性的影响

为寻求幂律型非牛顿推进剂有效的雾化喷嘴形式,采用与该推进剂流变特性基本一致的水基模拟液对直流式喷嘴、撞击式喷嘴、剪切撞击式气动雾化喷嘴、气泡雾化喷嘴和离心式喷嘴进行了雾化特性实验,比较了喷嘴形式及其结构参数对雾化特性的影响.结果表明:撞击式喷嘴是一种有效的雾化形式,撞击角度和射流速度是主要影响参数;气泡雾化喷嘴是一种高效的雾化方法,在约为4%的气液比下使水基模拟液破碎并产生液滴.      

移动粒子半隐式方法GPU加速的双股射流撞击雾化模拟

为实现双股射流撞击雾化过程的高效数值求解并探究射流速度和撞击角度对雾化特性的影响规律，实现了移动粒子半隐式方法（MPS）GPU加速的双股射流撞击雾化模拟。GPU加速程序的最大加速比为16，取得了较好的加速效果。将GPU加速MPS方法应用于典型工况下的双股射流撞击雾化模拟，成功捕捉到了多尺度的液膜形成、液膜破碎成液丝继而破碎成液滴的瞬态过程，模拟得到的液膜破碎长度及雾化角度与试验较为吻合，误差分别为11.7%和0.5%，验证了GPU加速MPS方法在双股射流撞击雾化问题中处理能力。参数化分析了射流速度和撞击角度对液膜破碎长度、雾化角度及一次雾化液滴索尔直径的影响。结果表明撞击角度增加或者射流速度增加均会导致液膜破碎长度减小、雾化角度增加、一次雾化液滴索尔直径减小。     

射流自由长度对凝胶推进剂撞击雾化影响的实验

为研究射流自由长度对凝胶推进剂撞击雾化的影响,建立了撞击雾化实验台,制备了煤油凝胶和水基模拟液,测试了雾化装置的喷射系数及模拟液的黏性和稳定性。分析了3种射流速度不同射流自由长度条件下的凝胶撞击雾化特性,观测了射流和撞击喷雾图像。测量了液膜破碎长度、雾化液滴粒径分布和相应的SMD(索太尔平均直径)值。研究结果表明:低速时,随着射流自由长度增大,撞击液膜的喷雾形态会发生较大变化,而高速条件下,雾化形态则基本一致。3种射流速度下,破碎长度在45~9mm之间,并随射流自由长度逐渐减小。液滴分布服从Rosin Rammler规律,并具有较高的拟合精度。均匀度指数均在3~4之间,并随射流自由长度逐渐降低,粒径均匀度降低;较高射流速度下,SMD随自由长度逐渐增大。低射流速度时,SMD随射流长度先减小后增大,射流自由长度存在一个最优值,在设定研究条件下其值为25/3。因此,在设计撞击喷嘴时,根据射流速度选取适当的自由长度值可以获得更好的雾化效果。      

双层旋转锥形液膜一次破碎特性数值研究

为了全面加深对锥形液膜一次破碎机理的认识,对双层锥形液膜的雾化过程进行了数值模拟,重点研究了压降和同轴旋转空气对双层液膜宏观形态、液膜破碎模式、液膜破碎长度和喷雾锥角等液膜一次破碎特性的影响。数值计算的喷雾场宏观形态与试验结果接近,喷雾锥角和索特尔平均直径的计算最大误差分别为4.9%、7.4%。研究表明：同轴旋转空气参与雾化会改变喷雾场的整体形态;增大压降和空气速度会改变液膜的破碎模式和主导表面波模式;双层液膜的合并会在液膜表面产生剧烈的表面波动,同时会略微增大液膜的喷雾锥角;液膜的破碎长度会随着压降和同轴旋转空气轴向速度的增大而减小。该研究有助于进一步研究双层液膜一次破碎的机理,从而指导对双路离心式喷嘴的雾化认识。     

基于树形自适应网格的旋流液膜雾化过程仿真

为实现旋流液膜雾化过程的数值精确求解,基于Gerris采用的自适应网格技术和VOF方法,建立了一种模拟旋流液膜雾化过程的数值方法,分析了雾化破碎过程及三维雾场特征。研究结果表明:Gerris能够逼真地展示旋流液膜破碎成液丝、液丝进一步破碎成液滴全过程的细节特征,雾化破碎过程图像与实验拍摄的基本吻合;通过统计分析计算的带旋转速度的直射流雾化过程全场液滴粒径空间分布,与文献中实验测量值也吻合较好,分布曲线峰值对应的液滴直径的差值为1.8μm,相对误差为13.8%,表明建立的计算方法具有较高的准确性。另外,通过对旋流液膜破碎过程的精细仿真,对其有了更清楚的认识,液膜雾化过程中存在二次雾化现象,液丝在运动过程中受到气体力和表面张力的作用,开始断裂形成大液滴或液团,随着进一步运动收缩破碎成小液滴,液滴形状渐渐由不规则的柱形变成类球形。带旋转速度的直射流和空心旋流式锥形液膜的液滴空间分布存在不同,前者液滴在锥形区域内都有分布,而后者液滴只分布在锥形液膜两侧的环形区域。