连续均匀热气流中液核/袋状破碎特性实验

为了获得热气流中单液滴液核/袋状破碎特性,采用高速摄像机对液滴变形、破碎过程进行了捕获.实验设计工况下所得结果表明:①液核/袋状破碎过程是液滴中液体流动与外界气动力共同作用下发生的,可以分为伞形破碎与勺形破碎;②液滴初始直径减小、气流温度增加,都可以降低液滴发生液核/袋状破碎所需最小气动力;③液滴破碎特征时间,随气动力增加呈负增长变化趋势,其值与变化梯度都随液滴初始直径增加而增大;④子液滴质量分数随气动力增加而增加,且随液滴初始直径增加而减小;⑤子液滴群质量分数的空间离散度都随气动力增加而增加.      

超声速气流中横向煤油射流的数值模拟

为了对超声速气流中液体横向射流进行数值模拟,采用Eulerian Lagrangian方法描述气液两相流动现象,用随机轨道模型追踪液滴的运动历程,并考虑气流可压缩性、流场不均匀性及液滴变形对液滴运动的影响,开发了相应的计算程序。计算结果发现,液体射流与气流之间存在着强烈的相互作用,液滴在进入气流中不久就破碎成很小的子液滴,受煤油液雾的影响气流速度、温度急剧下降,同时在喷孔上游出现一道弓形激波。与液雾结构的纹影图像和液雾穿透的实验结果对比显示,数值计算结果基本合理。     

超声速气流中液体喷雾流动数值模拟研究

超声速气流中液体喷雾流动特性对超声速燃烧基础研究及其工程应用具有重要意义。为了定量探索超声速气流中液体横向射流雾化特性，本文对超声速气流中液体喷雾流动进行了数值模拟研究。数值模拟方法基于Eulerian-Lagrangian两相流计算架构，考虑气液双向耦合，采用KH/RT液滴二次破碎模型计算液滴雾化过程，采用大涡模拟计算气相流动。结果表明，该数值模拟方法所获得的液雾场突起结构、穿透高度、液滴平均速度分布等液雾特性均能与试验结果较符合；初始液滴直径分布对破碎后液滴平均速度影响较小而对破碎后液滴平均直径及液相平均体积分数影响较大，初始液滴直径分布需在后续的建模与模拟中进行更多研究。     

超声速气流中液滴破碎的大涡模拟

为了研究超声速气流中的液滴破碎过程,采用基于CLSVOF界面追踪方法的大涡模拟方法。该方法中用可压流和不可压流求解器分别计算超声速气流流场和液体流场,并且在分别计算气相或液相时,另一相作为该相计算时的边界条件参与计算,从而实现可压流求解器和不可压求解的融合。通过去散度化液体速度外推解决界面处大液气密度比所带来的动量误差大的问题。利用该方法对来流Ma=1.358的超声速气流中的水滴破碎进行了数值模拟。模拟结果展示了超声速气流中液滴破碎的几何结构变化细节和基本破碎形态,且与实验吻合较好。数值模拟结果揭示了Rayleigh-Taylor不稳定和气动剪切力对液滴破碎的作用过程。破碎过程中,液滴在垂直来流方向上的最大无量纲宽度(Dmax/D0)约为4.5,无量纲破碎时间(tbU∞/D0)约为29.85。     

初始直径对单液滴破碎特性影响的试验

为了获得均匀横向气流中单液滴变形、破碎特性,采用高速摄像机对液滴的变形、破碎过程进行了测量.通过试验获得了液滴初始直径对液滴破碎模式、变形时间、变形与破碎总时间、无量纲索太尔平均直径的影响.在试验设计工况下所得结果表明:液滴破碎临界韦伯数随液滴初始直径的增加而增加;液滴初始直径增加使得液滴由单一的破碎模式转变成多种破碎模式共生的混合型破碎模式;相同韦伯数下,在初始直径尺寸较大的区域液滴变形时间、变形与破碎总时间都比尺寸较小的区域要大,变形时间随韦伯数增加呈先减小后不变的线性变化关系,而变形与破碎总时间随韦伯数增加呈先减小后增加再减小的线性变化关系;液滴初始直径对无量纲索太尔平均直径的影响能力要强于气液初始相对速度对其影响.      

连续均匀气流中单液滴破碎特性试验

为了获得均匀气流中单液滴变形、破碎特性,采用高速摄像机对液滴的变形、破碎过程进行了测量。试验设计工况下所得结果表明: (1)冷态气流中,液滴变形速率随[We]数的增加呈线性增加的变化趋势,[We]数在15～40液滴变形速率在0.2m/s～0.8m/s变化,[We]数在45～60液滴变形速率在0.6m/s～1.4m/s变化;(2)模糊数学模型计算结果表明,冷态气流中,液滴初始直径对液滴变形速率的影响能力弱于初始气/液相对速度,[We]数在15～40液滴直径、气/液相对速度的权重系数分别为0.29,0.71,[We]数在45～60液滴直径、气/液相对速度的权重系数分别为0.343,0.657;(3)液滴的初始直径增加使得液滴发生破碎的临界[We]数增大;(4)液滴的破碎时间分别随气流流量、气流温度的增加呈负指数的变化趋势,且变化趋势都随气流流量或气流温度的增加而减弱。     

超声速条件下亚毫米液滴的变形破碎模态

以超声速气流中液滴变形破碎行为为研究内容,对水平激波管中承受激波冲击的亚毫米水液滴(0.44~1.09 mm)变形破碎过程进行了观测,实验激波马赫数范围为1.07~2.11。利用纹影法,结合高分辨率高速相机对不同破碎模态下液滴的变形破碎特征进行了记录,得到了袋状、多模态、剪切和爆炸式等破碎模式下的液滴纹影图像,分析了液滴运动参数的时空关系。得出了液滴变形阶段,液滴无量纲横向变形宽度以及液滴无量纲迎风面位移随无量纲时间的变化发展规律,并且得出在液滴初始直径相同时,不同液滴破碎模式的无量纲最大横向变形宽度的变化,其中袋状、多模态、剪切破碎模式的无量纲横向最大变形宽度均在1.15~1.61范围内变化,爆炸式破碎模式的无量纲横向最大变形宽度均在0.21~0.68范围内变化。     

激波冲击下液滴变形破碎的黏性特征

为了获得不同黏度的液滴在高速气流中变形破碎的形态特征图像,定量分析黏性对变形破碎过程的影响,使用高速相机直接拍摄法在水平激波管中实验研究了液滴的变形破碎过程,测量了液滴迎风面位移、横向变形宽度、破碎时间等特征参数。所得结果表明:液滴的黏度较低时,液滴尾部形成的尾迹形状为细长的尖锥形,当液滴的黏度较高时,尾迹形状呈现波纹状,随着黏度的逐渐增大,液滴的尾迹更加复杂且紊乱;由于波后气流的作用,液滴迎风面发生变形失稳并演变出“尖钉”结构;黏性对液滴变形破碎过程起阻碍抑制作用;随着液滴黏度的增加,Rayleigh-Taylor(RT)不稳定性发展速率减缓,“尖钉”的数量先增加后减小。      

超声速冷态流场液体射流雾化实验研究

采用全息诊断和高速纹影,对不同来流总压、来流马赫数、喷孔直径和喷射压力等条件下超声速冷态流场液体射流雾化进行了研究。初步了解了超声速流场中液体射流的雾化过程和机理,得到了射流的Weber数和Oh数等雾化参数,比较了不同条件下射流穿透高度的差异,得到了液滴平均直径和数量密度的空间分布。研究表明：射流表面不稳定波的增长是超声速流场中射流破碎的主要原因;射流与气流的动量通量比和喷孔直径影响射流穿透高度,动量通量比和喷孔直径增加都会增加穿透高度;实验中液体射流的雾化过程非常迅速,在喷嘴下游20mm处,直径0.5mm的射流就破碎成平均直径10μm左右的液滴群,随着液滴向下游运动,平均直径逐渐减小,平均直径和数量密度分布逐渐均匀。     

剪切气流驱动液滴运动的受力分析

实验研究了剪切流驱动的液滴在固体表面上起始运动的受力机理。工作中使用一系列液体和固体表面来获得不同的液滴接触角,并在小型风洞中进行实验。实验中对液滴的启动气流速度进行了测量,并综合各种起始时刻的参数信息,建立了一个关于液滴接触线表面张力和剪切气流拖拽力平衡的数学模型,揭示了液滴脱落时刻的受力情况。所建立的模型更适合液滴1变形情况,但对于其它类似情况的剪切气流驱动液滴运动也能够进行合理的描述。     

不同液滴破碎模型的验证及柴油机在燃气混合过程中的应用

在射流不稳定性理论的基础上,将喷嘴内部的湍流流动以权重的形式加入初次破碎模型中;根据液滴的变形量采用TAB与KH-RT模型竞争的机制来计算液滴二次破碎;建立了Hybrid液滴破碎的数学模型。通过与实验数据、TAB和KH-RT破碎模型的计算结果对比,验证模型计算的准确性。Hybrid模型计算的液滴碰壁后飞溅半径和高度的最大相对误差分别为8%和6%,远小于其他两个模型的最大相对误差值。将验证过的三个破碎模型应用于柴油机中,比较了三个破碎模型对燃油分布的影响。研究发现,Hybrid模型的计算结果较比TAB和KH-RT模型准确度更高。      

过冷大水滴动力学特性对结冰影响数值研究

对比常规水滴,过冷大水滴结冰对飞机气动性能和飞行安全具有更大的危害性。在常规水滴结冰过程数值模拟基础上,针对过冷大水滴条件下各种动力学特性及其对成冰过程的影响进行了数值研究。采用结构化网格和中心有限体积法求解 N-S 方程获得空气流场,用拉格朗日法求解水滴流场,基于改进的 Messinger 热力学模型完成结冰过程模拟;对于过冷大水滴,采用泰勒类比理论和水滴碰撞模型,对比分析了水滴动力学特性及其对结冰过程的影响;通过数值模拟获得了水滴 Langmuir D 分布的撞击特性及结冰计算结果,研究了水滴多尺度分布对成冰过程的影响规律。通过一些典型结冰算例的数值模拟,并与参考文献结果及实验数据的对比分析,说明了本文计算方法及计算模型模拟过冷大水滴结冰是可行和正确的。     

变形破碎特性对SLD结冰过程影响

基于结冰过程数值模拟,针对过冷大水滴（SLD）条件下的变形破碎特性及其对成冰过程影响进行了研究。利用动态阻力计算模型,分析了水滴变形对水滴运动轨迹影响,采用基于泰勒类比理论的计算模型,研究了水滴破碎过程、子液滴粒径分布特性等。完成了NACA 0012翼型典型SLD结冰算例的数值分析,与参考文献计算结果和实验数据进行对比。结果表明:对于SLD结冰,变形和破碎主要改变了水滴运动轨迹和表面撞击水分布,使水滴撞击极限变小,在加入变形破碎特性的计算模型后,可以较好预测结冰极限位置和冰形轮廓,说明该方法对SLD变形破碎效应及其影响的模拟是可行和正确的。      

水滴动力学行为对其撞击特性影响

设计了一种压电式等径离散水滴发生器,利用高速摄影的手段对水滴变形现象进行了实验研究。建立了水滴变形后的非球形阻力模型,同时引入了水滴破碎和飞溅模型,采用欧拉法对不同条件下水滴撞击NACA23012翼型的过程进行了数值模拟,并与传统模型下水滴撞击特性计算进行了对比。研究表明:水滴变形会导致其阻力系数增大,水滴阻力系数曲线在雷诺数大于500时逐渐脱离球形阻力曲线,不断增大接近圆盘阻力系数;水滴阻力特性的改变会影响水滴运动轨迹,增强水滴的气流跟随性;水滴飞溅会造成水滴质量损失,减小翼型表面局部水收集系数,这与实验数据更加吻合。      

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

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

气液比对气液同轴离心式喷嘴缩进室内流动过程的影响

为了研究气液同轴离心式喷嘴缩进室内部非定常流动过程,采用Level Set和VOF相耦合的方法结合网格自适应技术对缩进长度为8 mm的液体中心式气液同轴离心式喷嘴流动过程进行了数值仿真研究,计算得到了较为精细的液膜一次破碎过程、流场结构和压力振荡特性。结果表明：液膜的破碎模式受气液比的影响较大,随气液比的增加,液膜破碎模式由曲张表面波主导的破碎变为穿孔破碎。此外,清晰获得了自激振荡过程,分析了缩进室内部压力场及速度场分布特征,研究发现随着气体压降增加,气体环缝出口会出现膨胀波和激波,形成一个“扇环形”的超声速流场区域；激波后气流分离,出现旋涡,形成局部高压区,旋涡中心随激波面周期性地上下移动,致使局部压力出现周期性振荡。     

超声悬浮乙醇液滴气动破碎特性研究

采用超声波悬浮液滴的方式,选取直径分别为1.0mm、1.5mm和2.0mm的乙醇液滴为研究对象,在韦伯数We=15~90的范围内,对其在脉冲气流中的破碎特性进行了实验研究。利用高速摄像机拍摄破碎现象,分析了液滴的破碎模式与运动特性。研究结果表明,随着韦伯数的增加,液滴先后产生袋状破碎、袋状/蕊心破碎、羽状/液膜稀释破碎。三种不同的破碎模式本质上都是袋状破碎——边缘袋状破碎与核心袋状破碎,低We数下主要为核心袋状破碎,高We数下主要为边缘袋状破碎。We数越大,破碎形成的袋直径越小,数量越多,破碎越剧烈。初始直径增大会使液滴转变破碎模式的临界We数增大。脉冲气流中乙醇液滴的迎风面运动遵循匀加速规律,其无量纲位移以We-1/2的比例缩小后是无量纲时间的简单二次函数。液滴在横向的扩散过程分为两个阶段,一个是惯性控制阶段,另一个是毛细效应阶段。     

过冷大水滴变形与破碎的影响因素

过冷大水滴（SLD）结冰过程中,影响水滴变形与破碎及其导致的效应的因素众多,多种影响因素耦合作用于水滴变形与破碎的规律还未得到清晰认识。采用数值计算的手段,模拟了SLD运动过程中的变形与破碎现象,分别研究了水滴直径、流场速度及绕流物体尺寸变化对水滴破碎与变形及其效应的影响。通过研究,提出了包含SLD的典型流场区域划分方法,得到了不同因素对水滴破碎区域面积、物面的水滴收集特性、破碎子水滴分布特性的影响规律,发现水滴直径、流场速度、翼型前缘半径与变形和破碎影响撞击特性的程度呈正相关关系。在此基础上,提出了一个综合影响参数X,线性回归分析的结果初步说明,X是表征水滴变形与破碎及其效应的合适参数。相关工作可以为SLD结冰的计算、试验以及防/除冰设计提供参考。