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

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

超声速可压缩火焰面模型的研究进展

火焰面类模型最初是针对低速流动提出的湍流燃烧模型,本文针对其向超声速可压缩流动的推广研究进行了综述。对于火焰面模型,其数据库建立时环境压强为常数以及静焓与混合分数线性关系的假设在超声速流动中均不适用。数值实验发现:一方面,数据库中主要组分浓度对于环境压强的变化并不敏感,但中间产物浓度在不同压强下变化明显;另一方面,静焓与混合分数偏离线性关系这一因素所引起的火焰面模型获得的温度场差别很小。对于火焰面/进度变量模型,其数据库中进度变量的化学反应源项对于压强和静焓-混合分数关系的变化则很敏感,现有考虑这一影响的可压缩修正主要是通过针对利用未经修正的数据库PDF积分后的平均源项进行标度来实现。数值结果表明,可压缩标度方法可有效考虑高马赫数效应对进度变量源项的影响,从而改善火焰面/进度变量模型对超声速燃烧流动的模拟精度。     

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

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

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

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

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

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

基于火焰面模型的超声速湍流燃烧数值模拟研究

在非结构混合网格有限体积框架下,将RANS方程和层流火焰面模型相结合,开展超声速湍流燃烧流动数值模拟研究,时间和空间离散分别采用LU-SGS格式和HLLC格式。用SST k-ω模型模拟湍流,FlameMaster 3.9生成火焰面数据库,并采用β-PDF和δ-PDF分布来考察湍流和化学反应的相互作用。采用该数值方法,对DLR氢燃料超燃燃烧室和北航乙烯燃料气动斜坡超燃燃烧室开展了数值模拟研究,并将数值计算结果与实验数据进行对比,二者符合较好。研究结果表明:在主流中心支板顺流喷射超燃燃烧室内,燃烧主要发生在支板后的燃料/空气混合层内,并且火焰核心落在恰当混合物分数Zst附近区域。      

建库温度在超声速火焰面模型中的影响

在火焰面模型中,数据库计算结果的准确程度对之后湍流燃烧流场的预测精度有重要影响。为了研究建库温度在超声速火焰面模型中的影响作用,在建立火焰面数据库时增加了建库氧化剂温度这一维度,并提出了一种由流场各点平均温度确定其对应的建库氧化剂温度的查表方法。针对DLR结构的超声速燃烧室,以煤油作为燃料,在五个温度级别条件下建立火焰面数据库,考察了建库温度对一维火焰面数据库以及二维燃烧流场计算结果的影响。研究结果表明,建库温度对一维火焰面的温度和组分分布有较为明显的影响,增加建库温度这一查表维度有助于改善超声速湍流燃烧流场的预测精度。     

值班甲烷/空气射流火焰的小火焰模型模拟

将小火焰模型和假定的PDF方法相结合,模拟甲烷湍流射流扩散火焰,稳态火焰面结构由层流小火焰数据库得到,采用详细化学反应机理描述甲烷的氧化和NO以及碳烟的生成,数值模拟结果和实验数据吻合较好,说明小火焰模型可以较准确地计算出扩散湍流火焰中的燃烧过程。      

采用不同建表方法的火焰面模型在燃烧室中的应用研究

为加深对航空发动机燃烧室中湍流燃烧过程的理解,采用不同建表方法的火焰面模型对航空发动机模型燃烧室内的湍流燃烧过程进行数值模拟,包括层流火焰面数据库的构建和反应进度变量的PDF类型两个方面。其中,层流火焰面数据库的构造方法包括基于扩散火焰的FPV和基于预混火焰的FGM模型,反应进度变量的PDF类型包括δ和β分布。LISA和KHRT模型分别用于模拟液膜和液滴的破碎过程,非平衡Langmuir-Knudsen模型用于模拟液滴的蒸发过程。LISA模型得到的液膜破碎距离约为4.6mm,液滴直径在文氏管出口下游迅速减小到10μm左右,并在头部出口下游附近完全蒸发。通过与相干反斯托克斯喇曼散射(CARS)和可调谐二极管激光吸收光谱(TDLAS)测量温度的对比,验证了FPV和FGM模型的精度,并表明在流动变化较大的位置FPV模型具有更高的精度,而其他位置FGM模型具有更高的精度,采用β分布作为反应进度变量PDF的模型,可以有效提高温度的预测进度,而且主燃区内的误差基本都在5%以内。此外采用β分布作为反应进度变量PDF的FGM模型,可以更好地描述未燃混合物被回流燃气点火的过程,而且反应进度变量的PDF类型比层流火焰面数据库构建方法的影响更为显著。     

采用不同建表方法的火焰面模型在燃烧室中的应用研究（英文）

为加深对航空发动机燃烧室中湍流燃烧过程的理解,采用不同建表方法的火焰面模型对航空发动机模型燃烧室内的湍流燃烧过程进行数值模拟,包括层流火焰面数据库的构建和反应进度变量的PDF类型两个方面。其中,层流火焰面数据库的构造方法包括基于扩散火焰的FPV和基于预混火焰的FGM模型,反应进度变量的PDF类型包括δ和β分布。LISA和KHRT模型分别用于模拟液膜和液滴的破碎过程,非平衡Langmuir-Knudsen模型用于模拟液滴的蒸发过程。LISA模型得到的液膜破碎距离约为4.6mm,液滴直径在文氏管出口下游迅速减小到10μm左右,并在头部出口下游附近完全蒸发。通过与相干反斯托克斯喇曼散射(CARS)和可调谐二极管激光吸收光谱(TDLAS)测量温度的对比,验证了FPV和FGM模型的精度,并表明在流动变化较大的位置FPV模型具有更高的精度,而其他位置FGM模型具有更高的精度,采用β分布作为反应进度变量PDF的模型,可以有效提高温度的预测进度,而且主燃区内的误差基本都在5%以内。此外采用β分布作为反应进度变量PDF的FGM模型,可以更好地描述未燃混合物被回流燃气点火的过程,而且反应进度变量的PDF类型比层流火焰面数据库构建方法的影响更为显著。