超声速气流中简化微小液滴横向喷射的汽化过程探索

为研究超声速气流中简化液滴的汽化过程问题,本文分析了两相流计算中已有的两相传热模型,并对简化液滴绕流开展数值计算.在来流Ma≤0.6的条件下,数值计算得到的简化液滴-气流之间的传热速率与已有模型得到的结果相一致,而在来流Ma≥0.9的条件下,数值计算得到的简化液滴-气流之间的传热速率与已有模型得到的结果存在很大偏差.由此建立了考虑简化液滴与气流相对超声速相互作用的两相传热模型.进一步,采用Charles B.Henderson阻力系数关系式与新建立的传热模型,对不同直径简化液滴的运动与汽化开展工程计算.在来流2.7Ma的二维平板超声速流场中选取一个截面,作为气相流场,结果显示,(1)简化液滴与主气流存在相对超声速作用.当简化液滴直径dk≤0.12mm时,作用区域约为0.1m～0.4m,当dk>0.12mm时,作用区域明显增大,(2)简化液滴的穿透尺度不超过0.011m/m(深度/长度),时间尺度约为0.28ms～3ms,(3)简化液滴完成汽化的空间尺度约为0.1m(dk>=0.03mm)、0.45m(dk>=0.05mm)与1.24m(dk>=0.075mm),而当dk>0.09mm时,简化液滴完成汽化的空间尺度则大于1.9m.使用考虑简化液滴与气流相对超声速相互作用的两相传热模型与使用传统的传热模型对简化液滴的运动轨迹没有影响,而对简化液滴的汽化过程有较大影响.     

基于双流体模型的超声速两相横向喷射流场研究

对两相超声速横向喷射试验进行了数值模拟.采用双流体两相流模型,在不同直径液滴和不同动压比条件下分别对两相超声速横向喷射流场进行了计算.结果表明:雾化液滴直径和喷射动压比对液体的横向穿透深度影响很大;雾化液滴直径的大小或者动压比的大小对喷口下游流场的影响区域是不同的.数值模拟结果与试验结果符合良好,为进一步的研究打下良好基础.      

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

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

有/无凹槽通道内煤油超燃雾化的测量研究

为研究超燃流场中雾化煤油液滴直径分布,针对有/无凹槽(L/D=3)直通道,对煤油在超声速气流中的雾化过程进行了测量。利用纹影法得到雾化流场波系结构,利用平面激光诱导荧光(PLIF)得到煤油扩散范围,基于激光前向散射原理的激光粒度仪,得到液滴直径分布。结果表明,当来流和喷射条件相同,平板通道内的煤油射流穿透深度略大于凹槽通道。液滴直径呈在截面中心小、外侧大趋势,沿外侧的液滴直径波动较中心区大。煤油距喷口30 mm处已完成雾化,雾化区的液滴Sauter直径为4~8μm,凹槽对煤油穿透深度和雾化尺寸无明显贡献。     

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

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

超声速气流中液体横向射流雾化破碎模型改进

为研究超声速气流中液体横向射流的雾化过程,根据可压缩条件下线性稳定性分析解推导出的不稳定增长率和波数间的色散关系,对KH/RT模型进行了改进,并调整了离散相液滴的追踪算法,利用改进的KH/RT模型对1.94Ma气流条件下水横向射流过程进行模拟。研究表明,数值模拟中液滴的实际追踪数对计算结果有较大影响,而改进的液滴追踪算法有效增加了液滴数,提高了计算准确性;改进的KH/RT模型扩大了Weber数的适用上限至4000,可用于超声速流场;利用改进的KH/RT破碎模型得到的液雾穿透深度与试验测量结果基本一致,而液雾展向和近壁区的分布与试验有一定差异。     

燃油直接喷射入横向气流中雾化机理数值研究

采用商业软件Fluent对燃油直接喷射入横向气流中的破碎雾化过程进行数值模拟,选用欧拉-欧拉方法中的VOF（volume of fluid）耦合level-set方法模拟气液两相流动.数值模拟了横向气流中燃油直接喷射射流破碎雾化过程及其相应流场信息,分析了流场结构对液柱破碎雾化过程的影响,获得了横向气流中射流破碎点位置以及射流穿透深度,并与试验数据进行了对比.数值计算结果表明:①采用VOF耦合level-set算法能够较好地模拟射流喷入横向气流中燃油破碎雾化过程;②由于反向旋转的涡对的存在,使得液膜和液滴沿着展向（y向）输运,对于液滴的二次雾化起到促进作用;③射流液柱总是在喷嘴下游约8倍直径处开始破碎,且破碎点位置取决于喷嘴几何特性.      

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

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

超声速气流中液体燃料雾化数值模拟

为了模拟超声速气流中液体射流的雾化,发展了一种液滴破碎混合模型。该模型将Kelvin-Helmhotz(K-H)模型和Rayleigh-Taylor(R-T)波动模型耦合在一起,首先计算R-T波增长模型,当R-T模型不能导致液滴破碎时,然后计算K-H波增长模型。利用该模型对来流Ma=1.94的超声速气流中的水射流的雾化进行了数值模拟。模拟结果再现了超声速气流中射流的雾化结构,计算得到的射流穿透深度、颗粒直径分布和速度分布与实验结果吻合较好。     

超声速气流中横向射流雾化实验和数值模拟

为了对超声速气流中的横向射流雾化过程进行研究,在直连式实验台上进行了超声速气流中的横向射流雾化实验,并对实验工况进行了数值模拟.超声速来流的马赫数为2,实验喷雾的动压比范围为1.～11.7,实验工质为水.采用纹影法对射流雾化过程进行了拍摄,得到了有雾化的超声速流场结构和穿透深度拟合公式.同时采用欧拉-拉格朗日两相流计算方法对实验工况进行了数值模拟,雾化模型采用了一种新型K-H(Kelvin-Helmholtz)和R-T(Rayleigh-Taylor)混合模型,数值模拟与实验结果符合较好.      

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

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

不同攻角条件下偏转头弹箭流场与气动特性

基于DES（detached eddy simulation,分离涡模拟）方法及AUSM（advection upstream splitting method,对流迎风分裂方法）格式,对不同偏角的偏转头弹箭在非零攻角下的超声速流场和气动特性进行了数值模拟.结果表明:在非零攻角情况下,与常规弹箭流场相比,偏转头弹箭头部迎风面激波更强;而当有效入射角相等时,即使弹箭的攻角不一样,其头部流场结构仍相同,且尾翼流场基本不受头部偏角影响.另外,由气动力系数变化曲线可知,偏转头弹箭具有良好的纵向静稳定性,且升力系数随攻角的增长率大于阻力系数,即升阻比大于常规弹箭,因此适当增大攻角可以提高偏转头弹箭的气动特性.      

自适应翼型的气动外形优化设计

 二维翼型自适应的研究是设计自适应机翼的基础。提出了在不同 Ma数、迎角下,用 Powell法优化二维翼型前、后缘襟翼的偏转角,以获得比常规翼型在亚声速时升阻比大而在超声速时阻力系数小的自适应翼型的研究方案。并与原始翼型以及气动双目标 (亚声速时,大升阻比;超声速时,小阻力系数 )的优化翼型进行了比较,证明了自适应翼型比气动双目标的优化翼型气动效率更高。初步探讨了二维翼型前、后缘襟翼的偏转位置对气动效率的影响。     

加热表面水珠运动特性研究

为研究加热表面的水珠运动特性,提出了加热表面水珠的几何参数、受力及其运动过程的计算方法。试验获得了水珠的表面阻滞力、黏性阻力和气动力计算关系式中的相关系数,给出了不同风速条件下水珠运动的临界直径,进行了加热表面水珠运动试验并对其过程进行了数值计算。试验结果表明：水珠的无量纲表面阻滞力保持恒定,运动时受到的黏性阻力与其运动速度和宽度有关,在外流场作用下所受的气动力可采用修正圆球阻力公式进行计算。将试验获得的相关系数加入水珠运动模型中,对加热表面水珠运动过程进行数值计算,计算结果与试验结果吻合,说明本文所述方法能够准确地模拟加热表面水珠大小和位置随时间的变化。     

曲面形栅格翼气动特性研究

与常规栅格翼布局相比，曲面形栅格翼布局减小了栅格翼翼元的等固壁通道的长度，从而有效降低了翼元内气流的壅塞，另外，这种栅格翼布局结构简单，易于折叠，减小了飞行器的轮廓尺寸，具有很好的工程应用前景。本文通过数值模拟方法，在亚、跨、超声速条件下，研究了曲面形栅格翼布局的减阻效果和减阻机理，对比了曲面形栅格翼和常规后掠形栅格翼布局的减阻效果，并对不同迎风方式的曲面形栅格翼气动特性进行了分析。     

超燃冲压发动机凹腔火焰稳定器阻力分析

为了研究凹腔火焰稳定器在超燃冲压发动机燃烧室内的阻力情况,采用了试验和数值仿真的方法,结合已有的研究结果、纹影图像和阻力测量等进行了分析。结果表明:在冷流中凹腔阻力受自由剪切层与后壁作用情况决定,来流总压高,作用面积大,凹腔阻力大;在反应流中,凹腔阻力小于冷流阻力,在高当量比条件下凹腔可能产生正推力;凹腔附近最强燃烧放热区对凹腔阻力影响较大,当最强放热区分布在凹腔前壁附近时凹腔可能产生推力;当最强燃烧放热区分布在凹腔后壁时,凹腔可能产生阻力。     

侧风环境下行驶的直背式轿车气动力计算

分别对直背式简化轿车模型在无侧风、稳态侧风、非稳态侧风三种条件下汽车周围流场进行了数值模拟,并对汽车受到的气动力进行了计算.结果表明,定常侧风会产生较大侧向力,阻力和升力也有一定增加;在非定常侧风作用下,气动力变化趋势与侧风速率变化趋势基本相同.阻力增大幅度较为平缓,侧向力增大较为显著,升力处于两者之间.当侧风增大到一定速率后,车外流场会发生较大变化,尾涡会转变为侧向涡.      

边界层抽吸效应对分布式推进系统气动性能影响数值研究

为探究边界层抽吸（BLI）效应对飞机的气动性能影响,采用基于沿流线体积力模型的飞机/发动机一体化数值模拟方法对某分布式推进系统进行了计算和分析。结果显示,BLI效应主要影响飞机中心体部分及发动机外整流罩的气体流动,对翼身融合体（BWB）的融合段及外翼段的升阻力影响较小。保持飞行条件和飞行攻角不变,飞机的升、阻力系数均随着无量纲化的发动机流量增加而增大,存在最佳无量纲化的发动机流量对应最大升阻比。发动机的安装位置直接影响机身表面的局部超声区和整流罩外的激波分布,布局在机身尾缘处会获得更好的升阻比,最大升阻比对应的无量纲化的发动机流量为0.65。      

低雷诺数下协同射流关键参数对翼型气动性能的影响

协同射流是一种高效的新型主动流动控制技术,至今缺乏关键参数对翼型气动性能影响规律的系统研究。通过引入螺旋桨激励盘模型,发展了一种新的协同射流翼型流动模拟方法,使得射流反作用力计算更符合实际。在低雷诺数条件下,以NACA6415为基准翼型开展了射流动量系数、开口尺寸和位置等关键参数对翼型气动性能的影响规律研究,并探讨了相应的物理机制。结果表明：大迎角分离流状态下,射流动量系数对翼型气动性能的影响规律比小迎角附着流状态更复杂;随着吹气口尺寸增加,气泵功率系数先减后增,有效升阻比先增后减;随着吸气口尺寸增加,气泵功率系数逐渐减小,有效升阻比先增加后趋于平稳;吹/吸气口位置对翼型气动性能和气泵功率系数的影响很小。