基于LES的非稳态火焰面/反应进度变量方法模拟部分预混抬举火焰

为了描述部分预混燃烧的火焰机制,捕捉燃烧现象中的非稳态效应,采用大涡模拟与非稳态火焰面/反应进度变量方法相结合来模拟湍流部分预混燃烧抬举火焰,并将计算结果与实验测量值进行了对比研究。结果表明基于大涡模拟的非稳态火焰面/反应进度变量方法捕捉到了湍流部分预混燃烧中的火焰抬举现象,计算得到的火焰抬举高度大致为35,和实验测量值完全一致。同时采用各热力学参数在混合分数空间的分布来研究部分预混抬燃烧中具有非稳态燃烧特性的局部熄火现象。各不同截面上的燃烧热力学参数与实验测量值吻合较好,说明采用该方法能比较真实地反映湍流部分预混燃烧的火焰机制。     

采用火焰面/反应进度变量方法模拟湍流燃烧

在自主开发的软件平台上添加了火焰面/反应进度变量方法,选取了两种反应进度变量分别对甲烷/空气的同轴射流燃烧室进行了数值模拟.将计算结果与稳态火焰面模型、文献参考值以及实验值进行了对比研究,结果表明:采用火焰面/反应进度变量方法能捕捉到稳态火焰面模型所不能预测到的火焰抬举等非稳态现象,得到的结果与文献参考值以及实验值出入很小,能准确反应真实的燃烧过程;另外不同反应进度变量的定义对计算结果有较大的影响,采用火焰面/反应进度变量方法对主要中间产物一氧化碳的预测比文献参考值更好,更接近于实验值.      

同轴射流燃烧室非预混湍流燃烧流场特性大涡模拟研究

为深入研究带突扩与回流的同轴剪切射流模型燃烧室流场结构与燃烧特性,采用单步快速化学反应、火焰面与反应进度变量三种燃烧模型结合动态亚网格模型在自有CFD平台上对甲烷-空气非预混燃烧流场进行了大涡模拟。计算检验了流场时空尺度及网格和时间步长设置,对流场0~2s内启动并发展进入统计定常状态的非定常流动过程进行了完整模拟。燃烧室内多处监测点湍动能谱统计结果表明,轴向速度在混合剪切层内的脉动存在从约1200Hz开始的多个特征频率序列,且随着测点位置向下游移动,其对应的特征峰从高频到低频逐渐衰减消失。流向多个截面上时间平均的混合分数及温度分布结果表明,反应进度变量模型能够模拟出该流场特有的抬举火焰,与实验数据吻合程度显著优于其余两种模型,即对于该类型强剪切非预混反应流,以燃烧产物质量分数替代采用梯度粘性输运模化的混合分数耗散率进行计算更为合适。     

超声速反应流火焰面/进度变量模型

基于Pierce低马赫数条件下的火焰面/进度变量模型,考虑可压缩性和激波的影响,发展了用于超声速湍流反应流的火焰面/进度变量模型.该模型通过求解混合分数和进度变量两个标量方程来描述化学反应,且采用火焰面方程解的完整的三个分支:稳态燃烧分支、具有局部熄火状态的非稳态分支和完全熄火分支.利用该模型对德国宇航中心(DLR)的超燃冲压发动机进行模拟,与稳态火焰面模型得到的附着火焰不同,计算结果捕捉到了熄火和火焰抬举现象;且数值阴影和平均温度分布与实验结果吻合较好.      

甲烷/空气射流抬举火焰的大涡模拟计算

为了研究中高温伴流中的自动着火特性,采用大涡模拟(LES)和均匀搅拌器(PSR)结合假定概率密度函数的建表湍流燃烧模型对甲烷/空气射流在高温伴流中的抬举火焰数值研究,并比较反应进度的不同概率密度分布对结果的影响。结果表明,计算结果与实验值符合得较好。PSR建表方法能够捕捉到高温伴流抬举火焰的点火过程,能够正确预测火焰的抬举高度,但在火焰稳定燃烧区域该建库燃烧模型还存在不足。同时,反应进度的亚格子脉动对准确模拟抬举高度十分重要。     

湍流扩散火焰局部熄火现象的大涡模拟研究

为了研究湍流扩散火焰的燃烧特性和局部熄火现象的机理,使用火焰面/反应进度模型对Sandia湍流值班射流火焰E进行大涡模拟。火焰面/反应进度模型适用于非预混燃烧,基于守恒标量和反应标量构建火焰面数据库,将详细的化学反应过程映射到这两个特征标量的平面上。数值模拟结果与实验数据符合得较好,火焰面/反应进度模型能较好地模拟湍流扩散火焰和预测火焰的不稳定性。研究发现,火焰在轴向距离与中心射流直径比值为7.5至15的截面处局部熄火概率较大,这与该位置标量耗散率偏大有关;火焰在轴向距离与中心射流直径比值为30以后的流场重新点燃,原因是混合物分数和当量比变小,火焰更多处于贫燃状态。     

超声速H_2/Air湍流扩散燃烧RANS数值模拟

为了研究修正的火焰面反应进度变量燃烧模型在超声速湍流扩散燃烧问题中的适用性,对德国宇航中心(DLR)超声速燃烧室开展RANS数值模拟。基于Open Foam软件平台中密度求解器分别对三维冷态场和燃烧场进行模拟分析。将网格自适应加密技术用于流场的计算;燃烧场计算中,通过分析不同压力下层流火焰面数据库,引入了反应进度变量源项的压力修正系数,压力修正系数α等于2.2。计算结果表明,冷态场中压力分布、波系分布、速度分布以及燃烧场中波系分布、速度分布、温度分布结果均与实验值符合较好。压力修正方法能够较好地解决超声速湍流扩散燃烧问题。湍流Schmidt数敏感性分析表明,湍流Schmidt数Sct对湍流火焰结构有较大影响,文中Sct等于0.7时能得到与实验值较为一致的分布。     

甲烷-液氧超临界压力非预混湍流燃烧的数值模拟

基于雷诺平均Navier-Stokes(RANS)方法和火焰面-反应进度变量湍流燃烧模型,在准确考虑热物性变化和详细化学反应机理的基础上,建立了适用于超临界压力非预混湍流燃烧的数值模型,开展了甲烷-液氧(LOx)超临界压力同轴喷射非预混湍流燃烧过程的数值模拟研究,着重探讨了推进剂混合比分别为1和3情况下超临界工作压力(6~15 MPa)对甲烷-液氧非预混湍流燃烧过程的影响。结果表明:不同混合比情况下压力对火焰温度和结构的影响会有显著的不同。在推进剂混合比为1时,随着压力的升高,火焰变得更长,且火焰温度升高;而在推进剂混合比为3时,随着压力的升高,火焰长度则会变短。在超临界压力下,火焰沿径向有突然的扩张现象(特别是在6 MPa压力下)。这主要与液氧物性在拟临界温度附近的突变所造成的拟沸腾现象有关,也会受到液氧喷射动量变化的影响。     

乙烯超声速燃烧部分预混火焰面模型数值研究

湍流涡结构的作用导致以扩散形式喷射燃料的超声速燃烧过程通常以部分预混的方式进行,为了使火焰面模型对这一燃烧过程描述的更准确,在模型中对扩散燃烧和预混燃烧状态均需要予以合理考虑。基于低Ma条件多区域火焰面模型(MRF)思想,发展得到了适用于超声速条件的部分预混火焰面模型。以碳氢燃料超声速双燃烧室结构作为验证算例,采用k-ωSST湍流模型、部分预混火焰面模型和乙烯(C2H4)的28组分72步化学反应机理对超声速湍流燃烧流场进行了数值研究。计算结果与实验数据吻合良好,验证了模型的准确性。加权系数计算结果表明,在中心剪切层和下壁面附近区域流场主要由预混燃烧控制,而其他大部分区域则由扩散燃烧控制。     

基于平衡化学反应和稳态火焰面模型的燃烧数值模拟

探讨了非预混燃烧模型在模拟甲醇燃烧室三维湍流两相燃烧数值模拟中的应用,简要介绍了非预混燃烧模型的优势.使用CFD软件对甲醇燃烧室的燃烧过程进行模拟,分别采用了平衡化学反应模型和稳态火焰面模型,模拟了甲醇燃烧室在试验参数状态下的温度场.通过对模拟结果比较后发现,两种模型在模型高温区域对温度的计算结果差别不大,但在低温区域差别较大.而平衡化学反应模型较稳态火焰模型达到平衡所需要的时间短.     

全流量补燃循环气气燃烧相似性缩尺试验研究

为指导全流量补燃循环发动机推力室全尺寸气气喷注器设计,采用气氢/气氧推进剂,在带可视化窗口的燃烧室中开展了气气燃烧流场相似性缩尺试验研究。采用高速摄影仪获得了不同流量工况下,同轴剪切喷嘴稳定燃烧和不稳定燃烧两种状态下近喷嘴区域的燃烧火焰结构,并分析了不稳定燃烧的频率特性。结果表明:在保持推进剂种类、推进剂混合比、推进剂温度、燃烧室及喷嘴结构尺寸不变的情况下,随着喷嘴流量的逐步增大,稳定燃烧和不稳定燃烧的喷嘴出口火焰结构均有一定的相似性,且不稳定燃烧的频率相同。     

等离子体助燃下超燃支板燃烧室燃烧特性数值研究

为了将支板喷注器与等离子体射流这两种促进超声速燃烧室燃烧的方式结合起来,设计了一种带有等离子体射流喷孔的支板燃烧室,并在超声速来流的条件下,针对燃料喷注总压、燃料喷注位置、等离子体射流介质、等离子体射流总压对燃烧室燃烧性能的影响进行了三维数值模拟。研究发现:增大燃料的喷注总压,燃烧室的燃烧范围明显增大,燃烧效率呈现出先增大后减小的趋势,在燃料喷注总压为2.0MPa时,燃烧效率达到最大值90.4%;不同的燃料喷注位置对燃烧室的燃烧范围影响较小;等离子体射流介质为O2时,燃烧效率最高,燃烧范围最广;提高等离子体射流的喷注总压,能够提升凹腔剪切层高度,有效促进燃烧,但同时也带来了更高的总压损失。     

某型航空发动机燃烧室煤油改柴油燃烧特性研究

燃料物性对燃烧室燃烧特性有着非常大的影响,是研究设计中首要考虑的参数。借助数值计算方法,采用Fluent稳态压力求解器、P1辐射模型和涡耗散模型(EDM)对某航空发动机燃烧室在巡航工况和最大工况下煤油与柴油两种燃料的燃烧特性进行计算及对比研究,得到该燃烧室使用航空煤油(RP-3)和0~#柴油的热态流场、空气流量分配、温度场、出口温度分布、污染物排放及头部燃油蒸发量。结果表明:在相同工况下,当该燃烧室的燃料由航空煤油改为0~#柴油后,其热态温度场分布基本一致,流量分配最大差异在0.45%之内,燃烧效率降低约4.3%,NO和碳烟粒子排放量相当,出口温度分布和总压损失差异分别在1.0%和4.1%之内。     

湍流燃烧模型对双旋流燃烧室喷雾燃烧的影响

采用数值模拟与试验测量相结合的方法,研究扩展旋涡破碎模型、扩展二阶矩模型和涡团耗散概念模型等三种湍流燃烧模型对双旋流湍流喷雾燃烧流场的影响.在任意曲线坐标系下数值研究双级轴向旋流器环形燃烧室全流程流场,采用粒子图像测速仪测量燃烧流场气流速度分布,热电偶测量燃烧室出口温度分布.计算结果与验证试验数据比较表明:不同湍流燃烧模型对双旋流湍流喷雾燃烧影响较大,所得的回流区形状、速度、温度场以及出口温度分布等都不太相同,其中扩展二阶矩模型所得的结果与试验值符合最好,更适用于模拟双旋流环形燃烧室湍流喷雾燃烧.     

基于时滞模型的纵向燃烧不稳定性分析

在集总燃烧时滞模型的基础上,研究了系统耗散对燃烧稳定性的影响。通过在流动控制方程中考虑剪应力引起的阻尼项,利用燃烧锋面处的匹配边界条件,推导了考虑耗散效应的集总燃烧时滞模型。其次,采用Tim Liuwen等的研究方法,通过对不稳定燃烧试车数据进行滤波操作和相关性计算,间接获得了各阶声学模态对应的耗散系数。将其作为阻尼项加入到新的模型中,重新计算了本次试车的稳定性曲线。结果显示,各阶耗散系数均小于0.1,且耗散系数会随着模态阶数的提高而降低。另一方面,考虑耗散后的系统增益较原始模型及现有模型降低了48%,说明系统耗散的加入影响了模型的预测精度。进一步研究证明,耗散的增加会降低系统增益,从而提高其稳定性。     

燃烧合成TiN陶瓷

采用燃烧合成工艺制备ＴｉＮ陶瓷件。压坯由７０μｍＴｉ加１５μｍＴｉＮ稀释剂经干压而成,孔隙率为４５％。压坯与８０ＭＰａＮ２反应合成ＴｉＮ陶瓷件。密度为理论密度的７５％。稀释剂为１５μｍ ＴｉＮ粉末,是由７０μｍＴｉ粉与８ＭＰａＮ２反应合成。对产物进行了ＸＲＤ相分析和ＳＥＭ微观组织分析。     

超燃冲压发动机燃烧室的燃烧特性

以一种低内阻光滑通道煤油超燃冲压发动机燃烧室为应用背景,采用有限差分法对燃烧室超声速流场进行了数值模拟.对流项采用3阶WENO（weighted essentially non-oscillatory）格式,湍流模型为SST（shear stress transport） k-ω模型,煤油（C₁₂H₂₃）/空气反应模型采用单步化学动力学模型.将燃烧室中沿侧壁的壁面静压的计算结果与实验结果进行了对比,结果符合良好,说明该算法适用于煤油超燃燃烧室计算. 研究了燃烧室来流静温、燃料/空气当量比和射流位置对煤油超声速流动与燃烧的影响.计算结果表明:燃烧集中在安装喷嘴一侧的壁面边界层附近,点火位置对当地静温非常敏感.随着来流静温降低、燃料/空气当量比减小和燃烧室扩张角增大,燃烧效率降低,燃烧性能下降,点火位置逐渐向燃烧室出口移动,燃烧放热形成的激波串结构消失.在燃烧室上、下壁面交错布置燃料喷嘴有利于提高燃烧效率.基于此,初步获得了光滑通道燃烧室内煤油点火燃烧的临界条件.      

爆震燃烧的特性研究

研究了爆震管中爆震燃烧的压力特性及爆燃到爆震转捩 (DDT)特性。在乙炔与氧气的预混气中 ,通过高频响压力传感器及自行研制的离子探针 ,测量了不同工况下爆震燃烧的压力与火焰传播速度的变化历程 ,获得了爆震波峰值压力、波后压力及 DDT距离随混气初始压力 (2×10 4 Pa～ 1× 10 5Pa)、混气当量比 (0 .3～ 1.0 )及混气浓度 (6 0 %～ 10 0 % )的变化规律。试验结果表明 :降低预混气的压力、混气当量比及浓度会使爆震波的峰值压力、波后压力不同程度的下降 ,DDT距离增大 ,其中 DDT距离对混气浓度最敏感     

含DNTF的改性双基推进剂燃烧性能

研究了3，4-二硝基呋咱基氧化呋咱（DNTF）的含量、不同催化剂体系、炭黑的含量、DNTF的粒度对DNTF改性双基推进剂燃速和压强指数的影响。实验表明：DNTF含量为30％时，以含能铅盐、芳香铜盐和炭黑（cB）所组成的复合催化剂可使推进剂压强指数降为0．37（8～14MPa）；而当DNTF的含量为50％时，所选用的三种催化剂体系（铅盐包括含能铅盐和芳香铅盐，铜盐为芳香铜盐）对推进剂的燃烧性能失去调节作用；炭黑含量的增加使推进剂的燃速增加，压强指数变大；DNTF的粒度对该推进剂的燃速影响明显，DNTF粒度较大时，燃速较高（53．19mm／s，16MPa），压强指数较大。