湍流两相燃烧不稳定性分析

采用数值模拟的方法对液体火箭发动机燃烧室内湍流两相流动与燃烧过程进行了分析,同时研究了燃烧室内高频不稳定燃烧现象。气相控制方程用欧拉坐标系下的Navier-Stokes方程组描述,液相控制方程在Lagrangian坐标系下描述,湍流模型采用高雷诺数的k-ε双方程模型,并采用高压蒸发燃烧模型。通过数值模拟分析了当燃烧室处于稳态燃烧时从喷注面到喉部之间两相混合物的变化规律;当发生高频不稳定燃烧时所研究的控制体内混合气体和喷雾的量发生振荡,并且具有相同的频率和半个周期的相位差。     

用双流体—轨道模型模拟四角喷燃模型炉内三维湍流两相流动和煤粉燃烧

目前炉内两相流动和煤粉燃烧数值模拟中多半用颗粒随机轨道模型和单流体无滑移模型，这些模型都难以完整地给出三维空间内颗粒速度，浓度，湍流度分布的信息。主采用双流体－轨道模型（颗粒相连续介质－轨道模型）对一个四角喷燃模型炉内三维湍流两相流动及煤粉燃烧进行了模拟。此模型基于欧拉气相方程组、欧拉颗粒连续方程组和动量方程组以及拉氏颗粒能量方程和质量变化的方程，并使用k-ε-kp两相湍流模型，EBU－Arrhenius湍流燃烧模型，离散坐标辐射传热模型，煤粉颗粒的水分蒸发，热解挥发模型和焦炭燃烧的扩散－动力模型等。热态模拟中，为了减小为信散造成的影响，采用了扭转坐标法（将坐标扭转一定的角度使之与煤粉射流方程一致）。为了检验数值模拟，采用三维相位移普勒测速仪（PDPA）对于冷态模型炉内湍流两相流场进行了测量，得到了两相速度，湍流脉动及颗粒浓度的分布。分别对冷态模炉内两相流动和热态模型炉内三维两相流动和煤粉进行了模拟，冷态两相流动的计算与实验结果的对比表明预报的两相流场是合理的，热态模拟的结果给两相速度，气相温度、组分浓度及壁面热，显示出靠近出口处气相速度和温度分布不对称，造成一个局部高温区。     

折流燃烧室两相喷雾燃烧流场数值模拟

在任意曲线坐标系下对折流燃烧室两相喷雾燃烧流场进行数值模拟,采用偏微分方程法生成三维贴体网格,采用k-ε双方程模型模拟湍流黏性,分别采用EBU(eddy break-up)-Arrhenius燃烧模型、EDC(eddy díssipation concept)燃烧模型以及二阶矩(second-order-moment,...     

湍流燃烧模型在航空发动机喷雾燃烧中的应用

使用自主开发的计算软件,数值模拟航空发动机旋流杯环形燃烧室三维湍流煤油两相燃烧过程,采用多维经验分析法与多步反应火焰面模型预测旋流杯环形燃烧室污染物排放,对4种湍流燃烧模型预测喷雾燃烧流场与污染物排放能力进行评估.各模型所得的计算值与试验数据较为相符,表明这些模型都可用来预测旋流杯环形燃烧室两相燃烧流场和污染物排放能力,但是不同模型的计算精度不相同.其中多步反应火焰面模型所得的出口温度分布与污染物排放与试验结果最接近,其预测精度高于其他模型.因此合理选择湍流燃烧模型可以提高计算准确性和可靠性.      

大涡模拟模型燃烧室燃烧性能计算

对带双级扩压器的模型燃烧室气液两相瞬态喷雾燃烧过程,在三维贴体坐标系下采用欧拉-拉格朗日两相大涡模拟方法进行数值研究,同时采用多维经验分析法预估燃烧性能.采用 k 方程亚网格尺度模型模拟亚网格湍流黏性;亚网格EBU(eddy-break-up)燃烧模型预估化学反应速率;多维经验分析法计算燃烧性能;并在非交错网格体系下气相采用SIMPLE(semi-implicit method for pressure-linked equations)算法对控制方程进行求解,液相采用随机离散模型,两相之间的耦合采用PSIC(particle-source-in-cell)算法.通过大涡模拟瞬态及时均计算结果表明:与粒子图像测速仪(PIV)测量的瞬态速度场、出口温度分布试验数据吻合,表明在三维贴体坐标系下采用欧拉-拉格朗日两相大涡模拟方法,数值模拟模型燃烧室两相喷雾燃烧流场,所采用的亚网格模型可以用于燃烧室气液两相喷雾燃烧流场的大涡模拟;燃烧性能计算结果与试验测量结果基本一致,说明所采用多维经验分析法可以用来数值模拟航空发动机燃烧室燃烧性能的计算,特别是污染物的预估,为设计低污染高性能航空发动机燃烧室提供有用的设计依据.      

液氢/液氧火箭发动机喷雾燃烧过程三维数值模拟

用三维湍流Ｎ－Ｓ方程及拟流体模型描述液体火箭发动机内部喷雾两相燃烧流动过程,两相之间的质量交换采用蒸发模型计算,气相化学反应速率采用Ａｒｒｈｅｎｉｕｓ公式计算。采用ＳＩＭＰＬＥＲ与ＩＰＳＡ算法求解两相控制方程并对氢氧发动机燃烧流动过程进行了三维数值分析,得到了发动机推力室内燃气参数的详细分布情况。结果表明对氢氧发动机而言,蒸发过程非常快,推进剂混合过程是燃烧过程的决定因素。     

斜切径向旋流器环形燃烧室数值模拟

采用三维贴体坐标系统,对包括突扩扩压器、帽罩、旋流杯、火焰筒以及内、外环通道的环形燃烧室全流程2相燃烧流场进行数值研究。采用k-ε湍流模型、2阶矩-EBU湍流燃烧模型、6通量辐射模型和颗粒轨道模型等模拟2相湍流燃烧流动,分析了进口工况对全流程燃烧流场的影响。计算结果表明:随着进口工况的改变,燃烧室出口温度场也发生相应变化;计算结果与试验数据比较表明:所用的数学模型合理、计算方法可行,其结果可为某型燃烧室优化设计提供可靠依据。     

全状态两相湍流燃烧模型及其检验

在欧拉-拉氏体系中,在传统点源两相湍流燃烧模型(TM)中考虑液滴带火燃烧状态,提出全状态两相湍流燃烧模型(FSM),并用甲醇-空气两相湍流射流火焰实验数据对两种模型进行检验,对比结果表明:FSM模型能够给出平均温度的双峰分布,得到更合理的预报结果.      

环管燃烧室热态流场的数值计算及性能分析

借助流体计算软件对某环管燃烧室内的热态流场进行了数值分析。气流的湍流流动采用Realizab lek-ε湍流模型进行模拟,用拉格朗日法跟踪油滴的轨迹,非预混PDF燃烧模型模拟气液两相的燃烧过程,辐射传热通过离散坐标模型来考虑。分析了气流流动规律、燃油液滴的运动轨迹、温度场的分布,获得了总压恢复系数、燃烧效率、燃烧室出口的径向温度分布系数等性能参数。计算结果可为环管燃烧室的优化设计提供参考。     

基于曲线坐标系浸没边界方法的折流燃烧室模拟

航空发动机折流燃烧室几何结构复杂,其高保真数值模拟需要高效的网格与边界条件处理方法。采用曲线坐标系隐式浸没边界方法结合大涡模拟-概率密度函数输运方程湍流燃烧模型开发自研软件,并实现WP11中折流燃烧室的高保真模拟。流动模拟中准确解析了该燃烧室中的三股主要气流,且三股气流分别约占进口流量的75%,12.5%和12.5%。两相燃烧模拟中针对拉格朗日框架下的液滴运动和欧拉框架下的湍流燃烧采用不同网格标记,模拟得到的出口径向温度分布规律与实验一致,平均相对误差为17.95%,表明基于本方法开发的自研软件能准确模拟折流燃烧室中的两相湍流燃烧现象。     

加力燃烧室流场形态与振荡燃烧数值模拟

利用数值方法对发动机加力燃烧室内的流场形态和振荡燃烧进行模拟,获得流场平均参数及涡系结构,观察到燃烧室内火焰稳定器唇口附近剪切涡对选择模型的敏感性,发现火焰稳定器后的涡脱落和回流区的周期性变化是引起加力"蜂鸣"现象和振荡燃烧的重要因素,证明了Rayleigh准则在研究此类问题时的适用性.      

加力室液雾燃烧的模型与应用

提出了描写航空发动机加力燃烧室内气液两相湍流燃烧的控制方程,采用了双流体模型,高温辐射通量模型,双时间尺度湍流K-ε模型,应用SIMPLE方法解这组方程,对二维轴对称的涡喷加力燃烧室进行计算,证明应用本文的液雾燃烧模型可以得出合理的计算结果。      

模型加力室大涡模拟的并行计算

采用基于MPI(消息传递库)的并行算法,在贴体网格下对带V形槽稳定器模型加力燃烧室紊流化学反应流场进行数值模拟,湍流模型采用k方程亚网格尺度模型,燃烧模型采用亚网格EBU模型,采用热通量辐射模型估算辐射通量。在程序设计中,采用动态内存分配、分区算法和多点重合交错网格系统,并行计算的结果与单机运行结果的对比表明计算结果是正确的,可以明显的提高运算效率,是解决复杂燃烧流动大规模数值模拟的有效手段。      

接力喷嘴不同径向高度和方位角对加力燃烧室 热射流点火性能的影响

为了研究不同热射流点火状态下的燃烧性能,针对采用波瓣混合器的某型航空发动机加力燃烧室,基于N-S方程建立了3维数值计算模型,得到了接力喷嘴径向高度和方位角对加力燃烧室流场、燃烧特性和流阻特性的影响规律。结果表明:随着径向高度增加,热射流火焰传播距离逐渐减小,传播到稳定器下游区域从内涵逐渐向外涵移动,且稳定器壁面高温分布区域逐渐减小;随着方位角增大,热射流火焰径向穿透深度逐渐增大,且稳定器壁面高温区域逐渐减小,在方位角α=0°和α=5°时稳定器壁温最高,为1450 K左右;在加力燃烧室出口截面上,径向高度和方位角对无量纲总压影响不大,整体小于0.005。     

内外函混合加力燃烧室扩压器内液雾湍流两相流场和浓度场数值模拟

 本文用κ-ε湍流模型和液雾轨道模型,对二维内外函混合加力燃烧室扩压器内有蒸发的液雾两相流动和混合进行了数值预报。其中对气相速度场和温度场的计算结果与实验符合令人满意,液雾两相流的计算结果在趋势上也是合理的,表明本文的模拟方法可用于加力燃烧室性能估算及优化设计。     

加力燃烧室中缝式稳定器技术研究

提出了一种加力燃烧室新型中缝式稳定器.冷、热态试验研究表明中缝式稳定器具有结构简单、重量轻、低流阻、高燃烧效率、宽广的稳定工作范围等综合性能.利用PIV激光流场测试技术研究了中缝式稳定器尾迹流场, 揭示其良好的燃烧性能取决于回流区结构和特殊的旋涡脱落方式.      

数值模拟纵向隔热屏加力室热态流场

数值研究纵向隔热屏加力燃烧室的在三维贴体坐标系下燃烧整体流场,所用的数学模型有:k-ε模型、EBU-Arrehenius和二阶矩紊流燃烧模型及六通量模型等.数值分析两种不同几何形状、进口气流参数分布及紊流燃烧模型等对加力室内各气流参数、隔热屏和加力室筒体壁面温度分布的影响,计算结果与试验数据比较表明:不同几何形状对加力室紊流燃烧流场的影响要比进口气流参数分布大些,二阶矩模型更适用模拟紊流燃烧流动.      

一种一体化加力燃烧室的数值模拟

为适应新一代航空发动机高推质比的设计要求,设计出一种一体化加力燃烧室方案,利用截尾支板与带凹腔的分流环组合结构取代了传统火焰稳定器。对该方案进行了数值模拟研究及试验验证,结果表明:在该一体化加力燃烧室内涵中有3个低速回流区;截尾支板结构不仅起整流支板的作用,还能够起到火焰稳定器的作用;燃烧效率在90%～93%之间,流阻系数约为0.26;在所研究的工况下总压恢复系数均高于0.975,且主要的总压损失集中在截尾支板及分流环凹腔处。      

一体化加力燃烧室冷态流动特性数值研究

将混合扩压器与火焰稳定器融为一体设计,并引入外涵冷流对高温结构件进行强化冷却,采用数值仿真方法,对比分析入口涵道比（0.278~0.583）和总温比（0.418~0.464）对一体化加力燃烧室中流动特征及其冷却性能的影响。结果表明:在研究参数范围内,大部分外涵冷流冷却隔热屏或直接进入加力燃烧室参与混合,仅有少部分进入火焰稳定器和联焰器,对其下游回流区没有影响,便于组织燃烧与稳定并传播火焰;外涵气流的流量分配对涵道比和温比变化不敏感;涵道比增加,特征面的流量系数下降约6.2%,总压恢复系数下降约3.2%,但是沿程热混合效率增加约6.1%;火焰稳定器的冷却效率随涵道比增加而提升约37.2%,但是其表面峰值温度降低不明显;温比对于该结构的流动特征与冷却特性影响都很小。      

气体燃料气动火焰稳定器混合特性的初步研究

为了克服航空发动机加力燃烧室传统钝体火焰稳定器存在的流阻较大、燃烧效率较低和红外辐射偏大的缺点，设计了一种气体燃料气动火焰稳定器，通过喷射气体燃料射流形成气动屏障来产生回流区从而稳定火焰，并采用数值模拟计算和实验测试的方法研究了气体燃料火焰稳定器的混合特性。数值模拟计算表明气动火焰稳定器掺混速度快，可在回流区内形成余气系数比较均匀的混合物,且回流区内余气系数分布随来流和射流的参数变化基本保持恒定不变，实验结果证实了数值模拟的结果，并表明采用气动火焰稳定器的燃烧效率较高，部分工况可达98%以上，可为加力燃烧室火焰稳定器的研究和设计提供参考和依据。