支板喷射超声速湍流燃烧的大涡模拟

为了发展可行的超声速湍流燃烧大涡模拟方法,将设定型PDF(Probability Density Function)模型与LES(Large Eddy Simulation)相结合以封闭亚格子湍流-燃烧相互作用,并将模型用于支板喷射超声速湍流燃烧流场的数值模拟。分别对冷流及燃烧流场进行了模拟,计算结果与实验测量符合较好,表明了所采用方法及模型的可行性。冷流条件下,大尺度湍流涡通过卷吸、拉伸运动主导支板尾迹区的近场混合,并通过破碎过程影响远场混合。燃烧条件下,回流区尺度扩大,剪切层中形成的高温燃烧产物通过大涡卷吸以及回流区末端对流作用进入回流区并与其中的燃料喷流相互作用,使部分燃料预热升温并进入回流区两侧剪切层与主流空气混合、燃烧,从而实现火焰稳定。在薄反应剪切层及大尺度反应涡的边界区域,LES网格不足以直接求解湍流与燃烧的相互作用,PDF模型给出了较强的亚格子脉动。     

超燃冲压发动机支板喷射燃料的燃烧过程试验

利用高速摄影和高速纹影对超燃冲压发动机支板喷注燃料的燃烧流动过程进行了试验研究,结果表明:燃料喷注器与火焰稳定器之间的相对位置,对流场的火焰结构和发动机性能有较大的影响,支板侧喷较尾喷方案具有更好的燃烧效果.燃烧与流动强烈耦合,燃烧放热使得燃烧室压力升高,改变了上游来流流场结构,流场结构的改变同时又影响到了燃料的喷注、扩散和混合过程.      

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

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

超声速和高超声速横向喷流数值模拟

为了解横向喷流的干扰影响和喷口附近的流场结构,采用分块对接网格和“0”型网格技术,数值求解N—S方程来模拟超声速和高超声速流场中横向喷流的干扰流场。对两种尖拱弹身外形的超声速和高超声速喷流干扰流场进行了数值计算,计算结果与风洞实验数据吻合一致。在此基础上,开展了某型导弹多喷构型干扰流场的数值模拟,得出一些横向喷流数值模拟的结论。     

气液喷流与超声速凹腔流场的相互作用

为了解凹腔火焰稳定器的工作过程,通过试验和数值仿真手段,对超声速条件下凹腔的流动特性进行了深入细致的研究,探讨了气/液喷流与超声速凹腔流动的相互作用机理。研究结果表明,无喷流时超声速凹腔流场具有五个典型的特征;引入喷流会引入新的流场特征,同时流场结构会发生机理性的变化;气/液喷流对凹腔整体流动特征的影响是一致的;有气体喷流时,不同压降下凹腔流场结构是类似的;而对液体喷流而言,提高喷注压降会增加雾化距离、射流穿透度、喷流厚度,提高来流速度则会使其减小。     

燃料引射方式对超燃冲压燃烧室混合及燃烧效率的影响

应用含组分守恒方程的质量平均Navier Stokes方程和B L代数湍流模型,数值模拟了后台阶构型燃烧室在采用台阶上游支板引射和壁面垂直引射燃料时的内部流场。在计算过程中,对方程的对流项采用空间为二阶精度的TVD格式,扩散项则采用二阶中心差分离散。通过流场计算,对比研究了引射方式对燃料混合性能的影响。结果表明,台阶上游的支板在燃烧室的流场中产生了一对相对稳定的大尺度轴向旋涡,该旋涡不利于燃料的混合。采用壁面垂直引射时,在喷嘴下游的燃料流场中产生了小尺度轴向旋涡,该旋涡是提高燃料混合及燃烧效率的关键。     

支板喷射超声速燃烧火焰结构实验

为明晰支板喷射超声速燃烧火焰精细结构,并探索其影响因素,利用自发辐射成像技术和平面激光诱导荧光技术(PLIF)对支板喷射超声速燃烧流场进行研究。结果表明,燃烧室内形成了典型的扩散火焰,且火焰稳定在支板尾部的回流区,宽度与支板厚度处于同等量级;支板厚度和燃料喷注当量比对火焰结构有显著影响,增加支板厚度会使火焰基底变宽,火焰长度变短,有利于燃烧室内部的快速放热和缩短燃烧室长度,增加燃料喷注当量比会使火焰长度增加,燃烧强度增大。     

横向脉冲喷流超声速燃烧流场数值模拟与分析

为了研究横向脉冲喷流对燃烧室流场结构及燃烧特性的影响,对带后台阶的超燃冲压发动机燃烧室横向氢气喷流超声速燃烧流场进行了数值模拟。采用有限体积法求解多组元Navier-Stokes (NS)方程,分别对定常喷流与两种脉冲喷流超声速燃烧流场进行了数值计算,对比分析了三种氢气喷注方案下的流场结构、氢气掺混燃烧特性及燃烧室总压恢复特性。研究表明:脉冲干扰仅对燃烧室内局部流场产生周期性影响,且有利于氢气的横向扩散;采用与定常喷注喷流参数相同的脉冲喷注方案时,可在脉冲干扰区内保持氢气掺混量与定常喷注基本相同,氢气燃烧量与总压恢复系数整体上与定常喷注基本一致;采用特征时间内喷出氢气质量流量与定常喷注相同的脉冲喷注方案时,在燃烧室出口处氢气掺混量和燃烧量较定常喷注分别提高了21.94%和32.24%,总压恢复系数仅减小4.17%。采用脉冲喷注方案增加了燃料与空气接触面面积,对燃料掺混及燃烧起促进作用。     

冲压燃烧室内的燃料扩散性能研究

应用含组分守恒方程的质量平均Navier-Stokes方程和B-L代数湍流模型,本文数值模拟了带支板引射器的模型燃烧室内部的流场结构,对比研究了支板引射器数量以及燃料注入角度对燃料混合性能的影响.在计算过程中,对方程中的对流项采用了空间为二阶精度的TVD和AUSMPW格式,扩散项则采用了二阶中心差分离散.采用流体力学时间步内的当地积分法克服了非平衡源项的刚性问题,保持了LU-SSOR隐式算法的高效率.研究结果表明,燃料有角度注入及采用多支板构型可提高燃料的混合性能.     

钝头双锥喷流致冷流场结构及密度脉动特性

采用大涡模拟方法对钝头双锥喷流致冷流场开展了数值模拟,研究了超声速喷流混合流场结构特征及密度脉动特性。大涡模拟方法基于隐式亚格子模型,空间离散采用高精度通量限制型紧致格式,时间推进采用显式Runger-Kutta方法。数值模拟清晰地捕捉到了流场波系结构,精细地预测了流动发生失稳、转捩以及发展为充分发展湍流的物理过程,直接获得了流场密度脉动特性。通过有、无喷流状态对称面流场的对比,发现超声速喷流能够有效冷却光学窗口;喷流与主流形成的混合层不稳定,很快发生失稳和转捩,形成大尺度湍流结构,进而引起强烈的密度脉动。此外,获得了钝头双锥整体模型喷流致冷流场的空间发展形态特征。      

基于火焰面模型的超声速燃烧混合LES/RANS模拟

为了明晰超燃冲压发动机燃烧室内部燃烧过程的细节,建立了超声速湍流燃烧稳态火焰面亚格子模型,并采用混合LES/RANS方法对氢燃料超燃冲压发动机进行算例验证.控制方程对流项用五阶精度WENO格式离散,时间方向采用二阶Runge-Kutta方法.研究表明:(1)冷流流场中燃料分布与大尺度结构分布相似,说明混合过程受大涡控制;(2)燃烧流场中涡的尺寸明显变大,且仅存在于火焰面上,另外温度分布和主要生成物分布与涡量云图基本相同,说明燃烧过程也由大涡控制;(3)时均计算结果与实验阴影基本符合,速度剖面和温度剖面与实验测量值定性一致,说明本文的数值模拟方法和燃烧模型可以较好地描述和预测超声速流动燃烧过程.     

带支板的冲压燃烧室的燃烧性能研究

应用含组分守恒方程的质量平均Navier-Stokes方程和B-L代数湍流模型,数值模拟了带支板引射器的模型燃烧室内部的流场结构,对比研究了燃料注入角度对燃烧及发动机性能的影响.在计算过程中,对方程中的对流项采用了空间为二阶精度的AUSMPW格式,对扩散项采用了二阶中心差分离散.研究结果表明,燃料有角度注入改善了燃料的混合性能,强化了燃料的燃烧效率,提高了发动机燃烧室的动量增量.此外,燃烧室的性能研究需考虑燃烧效率、总压损失及动量增量等因素.     

基于acetone-PLIF技术的超燃冲压发动机燃烧室燃料/空气掺混特性

针对超燃冲压发动机燃烧室内燃料/空气掺混特性,设计了一套以acetone-PLIF（planar laser-in-duced fluorescence）为技术手段的流动显示测量系统,得到了4种不同结构交错楔形支板在不同特征截面处的丙酮荧光图像,研究了喷嘴位置和尾缘角这两个变量对燃料/空气掺混效果的影响。结果表明:“基准”支板NL1的流场结构主体上表现为横向“S”形曲线,燃料/空气的掺混效果相对普通; NL2支板虽然增大了一些丙酮的扩散范围,但有荧光信号的掺混区域面积却反而有所下降;而引入扩张型尾缘角变量的EL1支板则是大幅改善了燃料/空气的掺混效果,54 mm截面的丙酮分子充满了较大片区域,掺混区域超过截面的1/2。      

甲烷/空气超声速燃烧流动数值模拟

用流动项与化学反应生成源项解耦处理的化学非平衡流动计算方法，从薄层近似三维N-S方程出发，采用ENO差分格式数值模拟了超声速冲压加速器简化模型中，高速甲烷气流从后体多个喷口射入高超声速空气流形成流场，研究了高温异质气体效应和19组分65反应模型的非平衡效应对冲压加速器表面压力分布的影响。计算表明新型解耦方法适合反应机理复杂的碳氢燃料超燃冲压发动机内部流动模拟，为开发应用软件系统打下基础。     

支板阻力特性实验

实验在超燃冲压发动机直连式试验台进行,燃烧式加热器出口气流马赫数为2.7,总温为1 337 K.利用推力测量系统对不同结构尺寸支板的冷流阻力和热试阻力进行了研究.比较了高度分别为37.5 mm和71 mm,后掠角分别为0°,30°和45°的支板冷流阻力,结果表明支板的冷流阻力随支板高度的增大而增大,随后掠角的增大而减小;实验还在燃烧模态下对不同当量比时支板的阻力进行了对比,结果表明支板的热试阻力比冷流阻力小,且支板的热试阻力随当量比的增加而减小.      

超声速燃烧凹槽火焰稳定的研究动态

高超声速气流在燃烧室的停留时间非常短，使得超燃冲压发动机燃烧室内燃料与气流的混合及其燃烧变得非常困难。目前集燃料喷射、混合及火焰稳定为一体的凹槽有望改进这一情况而受到了普遍关注。介绍了超声速气流流过凹槽的自激振荡及其控制、停留时间和阻力等特性，总结、分析了近几年凹槽在超声速气流中增强混合及燃烧的实验与数值研究，指出了在把凹槽用于提高超燃冲压发动机性能方面存在的问题和进一步研究的方向。     

基于气动斜坡的超燃冲压发动机双燃烧室方案研究

为提高超燃冲压发动机工作稳定性,提出了基于气动斜坡的超声速燃烧冲压发动机双燃烧室方案,该方案属于高超声速飞行器动力装置新方案。超燃主燃烧室采用基于气动斜坡的燃料喷注方式,并以小型燃气发生器作为亚燃燃烧室布置于气动斜坡喷嘴下游。超声速来流空气经进气道分流,96%左右进入超燃主燃烧室,4%左右经燃料电池驱动的离心式压气机增压后进入亚燃燃烧室。亚燃燃烧室在富油工况下工作,其出口布置在超燃主燃烧室气动斜坡喷注模块的下游(距气动斜坡第1排喷孔10倍喷孔直径处),此模块在主燃烧室中高效、低损失地形成流向涡。亚燃燃烧室喷流位于流向涡之后,起到点火、增强掺混和稳定火焰的作用。在直连式试验台上进行了该方案燃烧室部分的燃烧试验,结果表明:该方案成功实现了碳氢燃料大当量比范围内的稳定燃烧,以燃料比冲为评判标准,初步证明了该方案的可行性。     

Numerical simulations of turbulent flows in aeroramp injector/gas-pilot flame scramjet

To uncover the internal flow characteristics in an ethylene-fueled aeroramp injector/gaspilot(ARI/G-P)flame scramjet,a Reynolds-averaged Navier-Stokes(RANS)solver is constructed under a hybrid polyhedral cell finite volume frame.The shear stress transport(SST)k-x model is used to predict the turbulence,while the Overmann’s compressibility corrected laminar flamelet model is adopted to simulate the turbulent combustion.Nonreactive computations for Case 1(G-P jet on),Case 2(ARI jets on),and Case 3(both ARI and G-P jets on)were conducted to analyze the mixing mechanism,while reactive Cases 4–7 at equivalent ratios of 0.380,0.278,0.199 and0.167 respectively were calculated to investigate the flame structure and combustion modes.The numerical results are compared well to those of the experiments.It is shown that the G-P jet plays significant role in both the fuel/air mixing and flame holding processes;the combustion for the four reactive cases takes place intensively in the regions downstream of the ARI/G-P unit;Cases 4 and 5are under subsonic combustion mode,whereas Cases 6 and 7 are mode transition critical and supersonic combustion cases,respectively;the mode transition equivalent ratio is approximately 0.20.