基于大涡模拟的离心式喷嘴雾化过程研究

为实现离心式喷嘴雾化过程的精确数值仿真,探究喷嘴内部流动特性与外部液膜破碎形式,采用基于大涡模拟的仿真方法,对一种典型的四进口离心式喷嘴进行研究,仿真结果揭示了喷嘴内部相界面的振荡现象与外部液膜的破碎细节,并通过耦合流体体积法（VOF）与离散相模型（DPM）,获得液滴粒径的空间分布特征。研究结果表明：在液体填充过程中,喷嘴内的气液相界面存在波动与褶皱,形状并不稳定,内部的空气芯直径呈现正弦模式的振荡变化,喷嘴出口液膜厚度沿周向分布不均,这些因素导致出口附近的液膜表面出现扰动。在不同的进口条件下,不稳定性导致液膜表面上的扰动波形式不同。进口压力为0.3MPa时,液膜破碎由开尔文-亥姆霍兹（K-H）不稳定性产生的轴向正弦波所导致,产生沿周向分布的环形液带;在0.7MPa下,液膜表面开始出现由瑞利-泰勒（R-T）不稳定性引发的周向扰动波;随着压力增加至1.1MPa,液膜的破碎则由R-T不稳定性主导,产生沿轴向分布的液带结构,随后在气动力与表面张力的作用下破碎成液滴。二次雾化破碎后,喷嘴外部截面内的粒径呈“单谷”分布,液滴平均粒径计算结果与实验的最大相对误差为5.1%,与实验数据吻合度较高。     

低温环境下RP-3航空煤油雾化特性试验研究

为了获取低温对RP-3航空煤油雾化特性的影响,在综合试验台上完成了供油压力为300～700 kPa、燃油温度为25、-15、 -35 ℃的某喷嘴雾化特性试验。使用高速摄像机拍摄雾化图像,利用Matlab/GUI图像处理功能编制的程序提取雾化锥角边缘,并 获得雾化锥角,使用相位多普勒粒子分析仪测得距喷口7 mm和17 mm处截面的雾化粒度和油滴平均速度。结果表明：在压力一 定、温度降低时,雾化锥角变小,油滴雾化粒度增大,油滴平均速度升高,在同测量截面上的油滴平均速度变化量减小;随着压力升 高,-35 ℃和25 ℃的雾化锥角相差不超过1°,-35 ℃的油滴雾化粒度增幅减小到8.36%。     

基于压电纤维复合材料的离心叶轮叶片高阶模态测量

为准确获取离心叶轮叶片的高阶振动特性及其模态,利用压电纤维复合材料(Marco fiber composite,MFC)对离心叶轮叶片进行模态测试。设计了MFC激励系统,提出了MFC选型、激励和测点位置的确定方法,研究了不同信号的高阶模态激励效果,并分析了MFC附加质量和刚度对叶片固有频率和模态的影响;结合有限元仿真结果,对比分析了前20阶振动模态的异同。试验结果表明:MFC具有操作简单、重复性好、信噪比高、能够稳定激励出高阶模态等优势,计算模态与测量模态吻合较好,在第6～8阶与第12～15阶之间存在模态密集区域;MFC附加质量和刚度对实验固有频率的影响低于4%。提出的MFC激励方法可用于具有较高固有频率的航空发动机部件的模态测量。     

文氏管出口张角对旋流杯性能的影响研究

为了研究文氏管出口张角对旋流杯综合性能的影响,对不同文氏管出口张角的双级旋流杯开展了流量特性、下游流场和雾化性能试验,并借助仿真对试验结果进行了分析。结果表明:文氏管出口张角对旋流杯流量以及流量系数无显著影响。回流区外侧扩张锥面轴向速度随出口张角的增大先增后减。出口张角的存在可增大文氏管出口湍动能,强化两级旋流气体之间相互剪切作用进而改善雾化性能,存在一个最佳的角度(本文研究参数范围内,该值为56°)使得液雾平均粒径最小且液滴尺寸分布最均匀。拟合了关于文氏管结构的可用于预估双级旋流杯SMD值的模型公式,预估值与试验值吻合较好,相对误差小于20%。     

三旋流燃烧室喷嘴雾化性能试验研究

为了获得三旋流燃烧室喷嘴雾化性能与喷嘴几何尺寸的相互关系,针对供油压力以及旋流槽长宽比和旋流槽角度等关 键结构参数对燃油流量、喷雾锥角和雾化性能的影响进行了试验研究。采用3 维相位多普勒粒子分析仪测量了某一直线上各点的 索太尔平均直径和数密度分布,以及Rosin-Rammler 分布的特征直径和均匀度指数。结果表明：当供油压力提高时,燃油流量和喷 雾锥角增大;旋流槽几何尺寸的变化对燃油流量和喷雾锥角有不同的影响,当旋流槽长宽比和旋流槽角度增大时,燃油流量减小, 喷雾锥角增大。研究所获得的规律为三旋流高温升燃烧室的喷嘴优化设计提供了重要的理论依据。     

煤油两相连续旋转爆震燃烧室工作特性试验研究

为了深入分析煤油燃料两相连续旋转爆震燃烧室的工作特性,采用富氧空气或氧气为氧化剂,通过试验得到爆震波的时域、频域特征,对两相连续旋转爆震燃烧室中爆震波的起爆过程和稳定后的传播过程进行研究。利用基于激光散射相位多普勒分析(PDA)技术对雾化流场进行了测量,得到喷注器出口不同平面处煤油液滴速度与直径的统计分布。试验结果表明当煤油流量为78g/s,氧气流量为224.0g/s,空气流量为72.5g/s,当量比为1.083时,燃烧室在单波模态下工作,爆震波传播频率为0.904k Hz,平均转速为649m/s。使用氧气作为氧化剂,当煤油流量为81.8g/s,氧气流量为231.8g/s,当量比为1.222时,燃烧室在双波模态下工作,爆震波传播频率为5.882k Hz,平均转速为1848m/s,传播过程中表现出很强的非定常性。在当量比为0.805~0.908的富氧工况下,随着氧化剂中含氧量的增加,爆震波的速度逐步增大,最终达到2440m/s;在当量比为1.057~1.220的富燃工况下,随含氧量的增加爆震波速度呈现线性增长的特征。     

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

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

离心压缩机密封及空腔流动的一维/三维耦合建模研究

为了在保证数值预测精度的同时降低离心压缩机的建模复杂度,以一维密封泄漏计算替代完整的密封及空腔建模,开展了离心压缩机盖盘密封及空腔流动的源项建模方法研究。结果表明:基于源项方法的密封及空腔简化建模仅需在主流道建模基础上施加泄漏质量流量、总温和摩擦损失源项边界条件,即可较高精度地还原完整的盖盘密封及空腔建模的性能和关键流动细节。与完整建模相比,在全工况下源项建模预测的等熵效率的绝对偏差小于0.65%,总压比和等熵效率的相对性能偏差小于1%。该方法可在满足下游部件匹配设计精度需求的同时,减少50%以上的网格数量,并显著提高计算的经济性和收敛性。泄漏流量的计算精确度是影响源项建模法性能预测准确性的最关键因素,耦合建模可以采用一维计算中推荐的密封泄漏系数值。      

偏心撞击对撞击式喷嘴雾化特性的影响

为研究偏心撞击对撞击式喷嘴雾化特性的影响,建立了求解自燃推进剂冷态射流撞击雾化过程的数值模拟方案,计算了不同偏心度条件下的射流撞击雾化过程。采用树形自适应加密算法直接求解不可压Navier-Stokes方程组,由分段线性的流体体积(VOF)方法对流体界面进行捕捉。结果表明偏心撞击会导致雾场发生偏转,当无量纲偏心度E为1/8时,雾场偏转角度约为9.2°,应控制加工偏差小于该值。随着偏心度的增大,液膜的偏转角度增大,理论推导得到的液膜偏转角度要小于数值计算得到的液膜偏转角度。正心撞击时燃料与氧化剂流强峰值接近,雾场的流强分布呈单峰分布。当发生偏心撞击时,由于燃料与氧化剂部分射流未参与撞击导致流强峰值出现交错,雾场的流强分布呈双峰分布,混合比的空间分布发生较大改变。正心撞击时撞击点下游液滴的速度分布近似呈轴对称分布,而偏心撞击之后的速度分布则呈中心对称分布。偏心撞击导致的射流动量损失使得雾化性能变差,当无量纲偏心度E为1/8时,一甲基肼(MMH)的Sauter平均直径增大约4.8%,四氧化二氮(NTO)的Sauter平均直径增大约5.8%。      

液体火箭发动机推进剂泵诱导轮与离心轮的匹配

为获得诱导轮离心轮周向匹配的时序效应对离心泵外特性以及压力脉动的影响规律,阐释相关作用机制,采用基于分离涡仿真（DES）的离心泵三维全流道数值仿真方法,引入熵产理论以及压力脉动强度系数等先进分析方法对不同匹配角度下离心泵内能量损失机制及压力脉动特性进行了研究。结果表明：离心轮诱导轮的时序效应对泵外特性有一定的影响,随着匹配角度的增加,扬程和效率均呈现先减小后缓慢增大的趋势,扬程变化为0.8%,效率变化为1.2%,其影响机制由不同匹配角度下叶轮通道分离涡、叶轮叶片尾迹以及靠近隔舌处扩压器通道回流涡变化决定;时序效应对离心轮扩压器动静干涉效应影响显著,当诱导轮叶片尾缘位于离心轮相邻主叶片中间位置时,能够有效消除3倍频成分,显著降低泵内压力脉动水平,其中动静干涉区域以及隔舌处扩压器叶片表面压力脉动平均降幅分别达到14.5%和16.7%。