超声速燃烧气动斜坡喷注器研究

为增强超声速燃烧过程的燃料混合,设计了通过喷嘴阵列形成流向涡强化混合的气动斜坡实验件.在直连式实验台上对气动斜坡实验件进行了纹影及油流谱等冷试实验研究.并采用大涡模拟(LES)方法数值模拟了超声速流场中气动斜坡的流动及燃烧特性.结果表明:气动斜坡喷注器的喷嘴阵列设计能有效加速喷流后流向涡的卷起,提高掺混特性;其燃烧过程也能使喷流形成的剪切层上抬,增进混合;由于喷流方向与主流方向夹角较小,总压损失也较物理斜坡、横向喷流等混合增强方式小.      

子母弹不同舱段分离流场特性及运动特性研究

为研究子弹药在不同装配舱段下抛撒分离过程的干扰气动特性,以2层弹舱轴向排布子母弹系统为模型,基于有限体积法,结合非结构动网格技术,通过耦合求解任意拉格朗日欧拉(ALE)描述下的三维非定常可压缩N-S流体控制方程及六自由度刚体运动方程,对时序抛撒方式下子弹在不同舱段分离的三维非定常流场进行了数值模拟。得到了不同舱段下子母弹分离流场的干扰特性及子弹药气动参数变化曲线,分析了子弹分离过程的运动特性,揭示了子弹与母弹激波在不同分离阶段的相互作用过程。为进一步研究子母弹分离干扰流场机理提供了依据。     

火箭级间分离喷流干扰数值模拟与风洞试验研究

针对典型风洞试验状态,采用有限体积法求解三维雷诺平均Navier-Stokes(RANS)方程,数值模拟了给定级间距下火箭级间分离喷流干扰流动情况,获得了清晰的流场结构和有/无喷流干扰时两级气动特性变化规律,并对仅主机喷流和主机与四游机同时喷流的情况进行了计算,数值模拟研究结果为风洞试验方案的设计提供了参考依据。在典型状态数值模拟研究的基础上采用了冷喷流模拟技术对多个状态下火箭级间分离喷流干扰特性进行了风洞试验研究,获得了两级在有/无喷流、同轴变攻角情况下的气动力系数,研究结果表明,级间喷流干扰会对两级气动特性带来较大影响。相同状态下计算得到的两级气动特性变化与试验符合较好,表明了针对风洞试验状态的预先数值模拟对试验方案的设计验证可提供有益的参考。     

微型飞行器低雷诺数矩形扑翼非定常气动特性的数值模拟

采用双时间步方法求解三维可压缩非定常N-S方程,数值模拟了微型飞行器低雷诺数矩形扑翼的非定常绕流,首先将得到的结果与文献进行了对比,数据间具有较好的一致性.然后针对不同的展弦比、减缩频率及初始攻角,计算了矩形扑翼的非定常气动特性及表面流态和动态压力分布,并分析了翼尖涡对扑翼非定常气动特性的影响.     

求解可压缩N－S方程的无粘－粘性分数步法

用按物理过程作时间分裂的无粘－粘性分数步法求解可压缩Ｎ－Ｓ方程，无粘步用高效欧拉算法解欧拉方程，粘性步用全隐格式解抛物方程；并采取了有效措施消除由分步法产生的时间和空间不相容性，文中给出了轴对称喷管内外流场的数值模拟结果。     

机翼静气动弹性数值模拟研究

通过求解基于三维非结构网格的欧拉方程计算机翼气动特性,采用结构柔度影响系数方法进行机翼静弹性变形的计算,利用三维非结构弹性网格技术进行几何外形变形下的网格修正,将以上三部分进行耦合迭代.实现了机翼的静弹性特性的数值模拟.采用开发的数值分析程序,进行了某支线飞机机翼静弹性特性的计算分析,计算结果合理,表明所发展的数值模拟方法可用于机翼静弹性气动特性的分析研究.     

基于涡限制法的旋翼气动噪声计算

在直升机旋翼流动的数值模拟中,流动求解精度对旋翼气动噪声的预测有着重要影响.采用中心格式有限体积法求解可压缩Euler方程数值模拟旋翼绕流,并采用基于Ffowcs Williams- Hawkings( FW-H)方程的声类比方法计算旋翼气动噪声.为了抑制中心格式的数值耗散,保持机翼翼尖以及旋翼桨尖的尾涡结构,引入了涡...     

求解可压缩N—S方程的无粘—粘性分数步法

用按物理过程作时间分裂的无粘－粘性分散步法求解可压缩Ｎ－Ｓ方程，无粘步用高效欧拉算法解欧拉方程，粘性步用全隐格式解抛物方程；并采取了有效措施消除由分步法产生的时间和空间不相容性，文中给出了轴对称喷管内外流场的数值模拟结果。     

用DES数值模拟具有横向喷流的紊流流场

DES(DettachedEddySimulation)是把RANS(雷诺平均Navier Stokes方程)方法及LES(大涡模拟)方法结合起来的模拟有脱体涡的紊流流场的数值模拟方法。其主要思想是在物面附近解雷诺平均Navier Stokes方程、在其他区域采用Smagorinski大涡模拟方法。本文用DES及混合非结构网格数值模拟复杂紊流流场,算法采用Osher逆风格式,对超音速钝头体、三角翼、某型飞机及带喷流的导弹绕流用该程序进行了数值模拟。结果表明用DES模拟有分离的紊流流场是高效的。     

加油机喷流对受油机的气动干扰机理

针对软式平台(HDB)、软式吊舱(HDP)和硬式平台(FB)3种典型加油方案研究了喷流对加油流场的影响。采用结构化多块网格和雷诺平均Navier-Stokes(RANS)方程,通过带涡轮动力模拟器(TPS)和翼身组合体构型DLR-F6对数值方法进行了验证。通过对3种典型的空中加油方案进行数值模拟,并和不考虑喷流效应的空中加油流场进行对比,详细研究了不同加油方案中加油机喷流对受油机气动特性的干扰作用。结果表明:受油机升力系数、阻力系数以及低头力矩与无喷流相比均有增大。受油机在硬式平台加油方案中受喷流影响最小,但其横航向气动特性在软式吊舱加油方案中受加油机喷流影响较大。喷流对受油机的气动干扰主要表现在提高受油机来流动压、改变局部迎角以及对周边气流的引射作用,这些因素会导致受油机表面压力与无喷流相比发生较大变化。其中,局部迎角的改变是由喷流的加速效应及其对尾涡的耗散作用导致的。     

二维喷管的初始流动

基于可压缩Navier-Stokes方程,采用大涡模拟方法与高精度混合WENO/TCD格式,对Ma=1.4的超声速平面射流初始流场进行了数值研究.数值结果清晰地描述了超声速平面射流初始流场的结构特征,包括主涡环与激波结构以及它们演变过程.因主涡环内存在涡导激波对,激波与涡相互作用加速射流剪切层失稳,使剪切层首次卷起形成小涡的位置出现在涡导激波对后,此与亚声速射流情况不同.小涡串卷起成后,相继与涡导激波对相互作用,使激波出现明显的变形并加速主涡环失稳.      

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

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

超声速弱欠膨胀冲击射流流场结构

为了对超声速弱欠膨胀冲击射流的流场结构细节进行研究,使用大涡模拟方法对其进行了数值模拟。利用三阶迎风和四阶对称紧致格式对无量纲化轴对称可压缩滤波N-S方程进行空间离散,时间上推进采用的是三阶精度的TVD型Rugge-kutta法。亚格子尺度模型采用的是修正Sm agorinsky涡粘性模型。通过与经典的冲击射流实验比较,证明了程序的可靠性。数值模拟得到了剪切层以及壁面射流中的涡结构和主射流中的激波结构,并且在此基础上对涡合并和板前激波和涡干扰现象进行了深入研究。发现涡合并现象主要出现在流场的上游,越往下游出现的几率越小;涡和板前激波的相互作用会引起激波位置和强度以及冲击平板上冲击区的压强的显著变化,同时也会导致涡的变形。     

新型斜波瓣超声速混合器设计研究及其数值计算

根据超声速两股流掺混的特点, 在波瓣喷管的基础上提出了一种新的斜波瓣型超声速混合器设计方法, 设计的原则是在较小流动损失的前提下尽可能增加混合程度.利用有限体积法对三维可压粘流N-S方程进行离散, 紊流模型为k-ε模型, 分别对欠膨胀、接近完全膨胀和过膨胀3种工况进行了计算研究.计算结果表明, 当内通道瓣腔顶端与外通道相交处成一定夹角, 射流产生鞍形激波和旋涡可以增加掺混效率, 混合器内两股超声速流的流态基本符合设计要求.在欠膨胀、接近完全膨胀的工况, 效果是令人满意的, 但对严重过膨胀的工况, 到出口截面的混合效果不理想.      

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

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

基于自适应网格的超声速液体射流破碎过程研究

为研究超声速气流中液体横向射流的破碎过程,采用脉冲背景光方法和VOF方法开展了实验和数值研究。为提高液体横向射流中气液界面和气流场特征捕捉的精确性,采用自适应网格技术对于气液界面、激波出现位置进行网格细化,计算得到了较为精细的气液界面、激波特征及涡系结构。研究结果表明:在低成本仿真模拟条件下,利用自适应网格计算得到的射流轨迹和轮廓与实验吻合较好,射流轨迹的最大误差为10%;射流初始段在超声速气流条件下,仍然存在一段高度约为1.9倍喷孔直径且圆柱形态保持较好的连续光滑液柱。随着喷注压降的升高,液柱的长度逐渐增大;主流气体流经液柱发生三维绕流,在射流附近和近壁面区域形成不断演化的反转涡对,反转涡对的形成加速了液体射流一次破碎过程。     

超声速X形鸭翼-弹身组合体涡迹发展

应用蒸汽屏方法研究超声速X形鸭翼-弹身组合体涡迹发展。观察了起源于鸭翼后缘的四个翼涡在横截面上形成的“蛙跃”和上反角二翼涡与弹身一对对称脱体涡形成的“混合式蛙跃”现象。在临近蛙跃距离时,有不稳定特性发生。文中还给出了细长体理论计算的涡迹路径跟实验数据比较,结果表明:如果各个旋涡的初始位置和相对强度适合,这种数学模型可计算导弹上的各个旋涡路径,二者存在的差别,可能是由于计算未能模拟涡面和涡量耗散的缘故。为了有助于理解导弹的气动特性,用少量的离散涡计算涡迹路径,作为工程估算是适宜的。     

Review of atomization mechanism and spray characteristics of a liquid jet in supersonic crossflow

The injection and atomization process of a liquid fuel jet is critical for an ignition start of a scramjet engine. Airwall-mounted crossflow injection strategy is widely used in scramjet combustors, avoiding high total pressure loss and allowing the liquid fuel to rapidly undergo atomization, mixing, and evaporation. In this review, research progress on a liquid jet in supersonic crossflow was evaluated from aspects of atomization mechanism and spray characteristics. When a liquid jet is injected into a supersonic crossflow, primary and secondary breakups occur successively. The surface instability of liquid can significantly affect the breakup process. This review discusses the current understanding of the breakup process and spray characteristics of a liquid jet in supersonic crossflow including the mechanism of atomization and the characteristics of distribution and atomization. The development of windward Rayleigh-Taylor (R-T) unstable waves is the main factor in column breakup. The development of Kelvin-Helmholtz (K-H) unstable waves along the circumferential direction of the jet or droplets is the main factor of surface and droplet breakups. The liquid–gas momentum ratio is the most important factor affecting the penetration depth. The span width of the liquid jet is affected by the windward area. Breakup and coalescence lead to a transformation of the size distribution of droplets from S- or C-shaped to I-shaped, and the velocity distribution of the droplets on the central symmetry plane has a mirrored S-shape. The droplet distribution on the spanwise cross-section retains a structure similar to an “Ω” shape. At last, some promising recommendations have been proposed, namely a theoretical predictive model which can describe the breakup mechanism of a liquid jet, the distribution characteristics and droplets size distribution of a liquid jet under a cavity combustion chamber, especially for enthalpy flows with complex wave structures.     

超声速来流稳焰凹腔上游气体燃料横向喷注的流动混合特征

对Ma=1.7的空气超声速来流稳焰凹腔上游气体燃料横向喷注的流动混合特征进行了实验和数值模拟研究。在喷注的模拟燃料氦气中添加丙酮,采用平面激光诱导荧光技术(PLIF)观测了设置开式凹腔时燃料在展向和流向的空间分布。针对下游布置不同凹腔时燃料分布相似的实验结果,选取L/D=7的凹腔对其燃料混合流场进行了大涡模拟。实验和仿真结果表明:实验来流条件下喷注燃料大部分直接由喷流引起的上升反转旋涡对带走,只有小部分燃料经由喷流与凹腔流动涡结构的相互作用进入凹腔剪切层,并随剪切层运动对流进入凹腔内部。     

超声速高度欠膨胀冲击射流的大涡模拟

冲击射流广泛应用于短距起飞垂直降落飞行器(SVTOL)等航空航天领域.本文采用大涡模拟方法对高度欠膨胀的超声速冲击射流的流场进行了数值模拟.本文数值模拟得到了高度欠膨胀冲击射流流场中的激波结构和内外剪切层中不同尺度的涡结构.数值结果观察到了马赫盘的振荡,以及在斜激波、马赫盘及大尺度涡结构的共同作用下,射流内外剪切层之间的环形激波的生成与消失的周期过程.并对流场内剪切层的涡结构的演化进行了研究,数值结果显示内剪切层的大尺度涡结构的形成与马赫盘的振荡相关,在内外剪切层的作用下形成了壁射流区内外交错的涡结构.