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

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

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

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

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

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

针栓式喷注单元雾化角模型分析

为了实现不同径向孔形的针栓式喷注器雾化角的准确预测，从动量守恒方程出发建立了液膜撞击液膜和液膜撞击液束的雾化角理论修正模型。对于液膜撞击液膜的喷注单元，模型中通过理论推导引入了2个变形因子，将撞击的几何变形效应与雾化角关联；对于液膜撞击液束，通过引入阻塞率定义有效撞击动量比，同时将液束入口孔形的影响隐含考虑在变形因子中，最后根据高速摄影试验结果和数值仿真结果获得了对应的变形因子组合系数，使得新的雾化角模型适应性更广、准确性更高。结果表明：引入变形因子和阻塞率的理论模型预测值与试验及数值仿真结果吻合很好；对于液膜撞击液膜，变形因子基本维持在0.9~1.1，根据试验结果及仿真结果，变形因子推荐值为C₁=0.99和C₂=1.06；对于液膜撞击液束，变形因子推荐值为C₁=0.75和C₂=1.25。该模型根据实际出口的轴向动量和合成总动量计算雾化角，隐含考虑了撞击作用造成的影响，较根据撞击前入口的轴向动量和合成总动量计算雾化角的常用模型预测值准确度显著提高，为针栓式喷注器的理论研究和工程设计提供了重要参考。     

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

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

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

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

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

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

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

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

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

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

倾斜射流撞壁实验研究及液膜几何参数建模

倾斜射流撞壁在液体火箭发动机液膜冷却、射流撞壁雾化等领域具有广泛的应用。为了研究倾斜射流撞击壁面后形成的液膜的关键特征，开展单股圆柱射流撞击壁面的实验研究。从实验中研究各射流参数对液膜外形的影响规律，继而开展理论建模，获取液膜外形的关键几何参数表达式。实验研究发现随着射流倾角增大，液膜长度减小而宽度增大，随着射流孔径和射流速度增大液膜长度和宽度均增大这一定性规律。理论分析得到了液膜最大宽度位置与液膜对称面的夹角近似等于射流倾角α，液膜的长宽比近似等于1+cot α。通过进一步分析得到液膜宽度、最大宽度距撞击点的距离、液膜长度这几个液膜关键参数的表达式。建立的液膜几何参数表达式预测结果与实验结果的误差均在20%以内。     

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

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

剪切稀化非牛顿射流撞击液膜破碎直接数值模拟

液态射流撞击是液体火箭推进系统中广泛采用的一种燃料雾化方法,其破碎特征直接影响燃料最终的掺混及燃烧效率。采用直接数值模拟（DNS）工具,研究了低雷诺数（Re_l=41）和中等韦伯数（We_l=163）条件下剪切稀化非牛顿射流撞击液膜破碎的问题,着重分析了对角液膜的三维结构、破碎特征和非牛顿特性等。研究结果表明：在所研究的射流参数下,该非牛顿撞击液膜破碎属于Open Rim类别,破碎过程具有三维特性并伴随液丝与边缘的融合、液丝向液滴的转变等时域流动特征。液体的总表面积随时间不断增长,但单位表面积随液膜破碎的发生而下降,液膜扩张半角随时间逐渐增加并趋于恒定值43°,而后部液膜的长度几乎不随时间发生变化。此外,撞击液膜表现出明显的剪切稀化特性,液体内部最低黏性系数仅为零剪切黏性系数的1/5。     

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

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

强迫扰动下的射流撞击雾化特性

为全面把握撞击式喷嘴的工作特性，进一步认识雾化在燃烧不稳定中所起的作用，采用试验结合数值模拟的方法开展了强迫扰动条件下撞击式喷嘴的非稳态雾化特性研究。试验方面，采用水力扰动装置产生喷前压力扰动，由脉动压力传感器记录喷前的脉动压力，由高速摄影对动态的喷雾场进行背光拍摄。数值模拟则是基于开源程序Gerris开展，通过给定周期性变化的速度入口来模拟前端压力扰动下的撞击雾化过程。首先验证了建立的数值模拟方案处理非稳态雾化的有效性，其次将自然雾化与强迫扰动雾化进行对比，分析了强迫扰动条件下的撞击雾化特性，最后研究了扰动频率与扰动幅值对于撞击雾化的影响。结果表明：强迫扰动下的射流撞击喷雾场出现了弓形液滴群局部聚集的现象，并且在时间上表现出周期性特征，雾化频率与强迫扰动的频率一致。在研究的频率范围（1 257~3 563 Hz）内，撞击式喷嘴的雾化对扰动都有响应。扰动频率主要影响相邻弓形液滴群之间的间距以及雾场与扰动压力之间的相位关系，扰动幅值则主要影响雾化Klystron效应的强度。随着扰动幅值的增大，液膜的破碎长度减小，撞击点下游的流量特性由线性向非线性转变，由正弦波形转变为陡峭前缘波形。     

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

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

平板式预膜喷嘴初次雾化特性试验

为深入了解平板式预膜喷嘴初次雾化特性,试验研究不同进口条件对平板式预膜喷嘴初次雾化特性的影响规律。试验采用了背光照明和片状激光照明相结合的高频拍摄手段分别获得液膜俯视和侧视破碎形态,同时运用本征正交分解(POD)法和液膜边缘定位等分析方法进行光学图像结果后处理,获得表征平板式预膜喷嘴初次雾化特性三个物理量:液膜波动频率、破碎距离和横向不稳定波长。试验结果表明:①通过分析侧视和俯视破碎形态,平板式预膜喷嘴液膜破碎形态可分为三类:末端破碎、波浪脱落和表面剥离,进口韦伯数对预膜喷嘴破碎形态的影响占主导地位;②把POD法和液膜边缘定位方法等相结合方法应用到高频非接触光学喷嘴雾化图像的后处理分析中,是一种非常有效的数据后处理方法;③液膜初次雾化特性主要受到进口韦伯数和气液动量比的影响,破碎距离和横向不稳定波长都随进口韦伯数的增加而降低,液膜表面不稳定波动频率随进口韦伯数的增加而增加,所获得的经验关系式与试验数据吻合较好。所获得的进口参数对液膜破碎形态和雾化特性的影响规律为喷嘴后续优化设计提供了依据。     

喷注角度对超声速气流中液体横向射流的影响

为探索超声速气流中变角度组合射流对喷注雾化过程的影响，采用高速摄影和基于KH-RT破碎模型的欧拉-拉格朗日数值仿真方法开展研究，仿真结果和实验吻合较好。研究表明，同一流量下逆流喷注，可以较为有效地增强混合效果，提升射流穿透深度。带扩张型面内液雾穿透深度呈现先增大后降低的趋势。等流量条件下顺流双孔喷注有效截面面积大体一致，略低于单孔直喷。等效流量下逆流喷注所得液雾截面面积收益最为明显。截面平均粒径方面，单孔喷注雾化粒径基本一致，顺流双孔喷雾粒径较小。等压降顺流对喷由于流量较大，雾化质量相对较差，平均粒径最大。     

低韦伯数非牛顿射流撞击破碎直接数值模拟

非牛顿射流的撞击破碎在液体火箭推进系统中被广泛用于燃料的喷注雾化,然而人们对其破碎物理机制却知之甚少。本文将采用基于液体体积法的直接数值模拟(DNS)工具,研究夹角为90°的2个等直径低韦伯数射流撞击现象,并分析二者形成的单一对角射流特征及其破碎机理。研究结果表明,撞击形成的单一对角射流直径较原射流直径大1.66倍,并在头部形成液滴诱导破碎的发生。除了头部破碎,在对角射流的发展过程中还观察到一类液柱破碎,表现为射流表面不稳定波不断发展形成新的弯曲波破碎,并产生卫星液滴及液滴的融合。伴随两股射流撞击的发生,气液两相交界面的面积也不断减小,同时,射流内部的黏性也不断变化,在本文的低雷诺数和低韦伯数条件下,流体内部黏性系数变化超过10%。     

非牛顿流体射流雾化特性研究进展

本文总结了有关非牛顿流体射流雾化特性的研究进展。首先,阐述了预测非牛顿液体射流初次雾化失稳特性的理论方法,介绍了有关非牛顿流体射流初次雾化的实验现象和特性参数。当射流初次雾化的过程结束后,破碎产生的液滴会在高速气流中发生二次雾化。随后,总结了国内外有关非牛顿流体液滴二次雾化实验研究的相关进展。分析了液滴二次雾化的实验现象,总结了不同种类液滴二次雾化过程中所研究特性参数,如破碎模态、临界韦伯数和初始变形时间等随来流气体参数之间的关系,并介绍了基于液滴二次雾化物理过程所建立的预测喷雾场液滴平均粒径的雾化模型。最后,基于目前的研究现状,给出了非牛顿液体射流初次雾化和二次雾化实验研究的后续重点研究方向及建议。     

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

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