强迫对流下硼颗粒燃烧特性影响因素研究

针对冲压发动机燃烧室内强迫对流下的硼颗粒燃烧特性展开了系统研究,考虑气相流动、扩散和表面单步有限化学反应动力作用,建立了强迫对流下硼颗粒燃烧过程的物理和数学模型;采用有限体积法求解含多组分反应流的二维轴对称Navier-Stokes方程,并验证了数值仿真方法的正确性。首先通过数值仿真研究了来流速度、颗粒半径、环境中氧气质量分数和环境压力等因素对硼颗粒燃烧特性的影响,并对其成因展开了详细分析。研究表明,在强迫对流作用下,硼颗粒总的燃烧质量流率和质量流率通量均随来流速度、颗粒半径、环境中氧气质量分数和环境压力的增加而增大。通过深入分析发现,强迫对流下硼颗粒的燃烧质量流率通量随着来流雷诺数的增加而增大。然后基于大量数值仿真结果,对相对静止气氛下的硼颗粒质量流率通量进行了修正,用于描述强迫对流下的硼颗粒燃烧特性。     

HTPB固体燃料冲压发动机流场仿真与燃速分析

基于守恒方程建立了固体燃料冲压发动机燃速仿真模型,采用二维轴对称模型和二方程化学反应模型开展了HTPB(端羟基聚丁二烯)固体燃料冲压发动机流场数值仿真,分析了不同空气来流条件对流场分布及燃速的影响.结果表明:火焰层在氧气和固体燃料壁面之间形成,随着来流空气流量和空气总温的增加,火焰层厚度变薄并向固体燃料壁面侧移动;随着发动机轴向位置的增加,燃速先迅速增加后缓慢增加,最后在补燃室附近快速减小,变化趋势与文献中试验结果吻合较好;固体燃料平均燃速随来流空气总温及发动内空气流率的增加而增大,并根据仿真结果拟合得到了燃速公式.      

激光辐射下AP/HTPB复合推进剂低压燃烧特性数值模拟

为了研究激光辐射对AP/HTPB复合推进剂低压燃烧的作用机理，建立二维三明治稳态燃烧模型进行数值仿真研究。针对激光辐射强度1.4W/mm²，压强30～90kPa的工况下推进剂的燃烧进行仿真计算，并与130～160kPa压强下自持燃烧的仿真结果进行差异分析。结果表明，激光助燃工况下，仿真计算所得燃速与实验基本一致，最大误差不超过4.2%。在激光助燃和自持燃烧工况下，组分扩散速度快，气相火焰均为预混结构。激光助燃工况下，激光辐射的能量使得近燃面处的反应比较充分，初扩散反应在气相反应中占主导地位，气相火焰高温区域离燃面近；自持燃烧工况下，AP预混反应和初扩散反应共同主导气相反应过程，气相火焰温度在垂直于燃面方向上增长缓慢且高温区域离燃面远。激光助燃比自持燃烧时的压强指数低，推进剂燃烧对压强的敏感性较弱。     

采用N_2O与含金属HTPB燃料固液火箭发动机燃速试验研究(英文)

为了研究含金属固体燃料固液火箭发动机的燃速及发动机的工作特征,采用N2O氧化剂与含Al,Mg和C的HTPB基固体燃料的固液发动机进行试验研究与分析,试验在3种尺寸药柱、多种氧化剂流量下进行,获得发动机的燃速、燃烧室压强及燃烧效率等参数。为分析辐射换热对燃速的影响,把燃速分成两部分并提出相应的分析式:一部分由对流换热控制的燃速,由通过药柱通道的总质量流率来衡量;另一部分由辐射换热控制,通过燃烧产物中凝相组分对固体燃料壁面的辐射换热量决定。对试验结果进行分析发现,除富氧的情况外,热辐射控制的燃速约占总燃速的30%~60%,辐射控制燃速与对流控制燃速的比例趋势与燃烧产物中凝相组分的质量分数随氧燃比的变化规律相似。同时,试验结果显示,试验燃速总体上随着固体药柱通道中氧化剂质量流率的增加而增大。发动机的燃烧效率绝大多数位于80%~97%的范围,在化学当量比附近和富氧范围内,随着燃烧温度的降低而降低,效率曲线与理论绝热燃烧温度值的平方相似。     

聚乙烯在固体燃料冲压发动机中燃烧特性的数值研究

为准确模拟聚乙烯在固体燃料冲压发动机中的燃烧状态,估算了乙烯-空气两步反应中,相关组分的粘性、导热系数随温度的变化,并进行多项式拟合;同时使用UDF进行质量、动量、能量添加.数值计算的结果表明:燃面退移速率的计算误差不超过11.3％,其变化规律与已有实验相同;补燃室温度的误差不大于5.2％;补燃室温度随来流空气质量流率的增大而减小;高来流空气总温或低来流空气质量流率的工况下,燃料充分燃烧所需的补燃室更长;几何相似,其它条件相同时,小尺寸发动机内流场传向燃料表面的热流密度更大,燃面退移速率更高.     

UDMH/NTO火箭发动机尾焰流场特性数值仿真

采用AUSM(advection upstream splitting method)空间离散格式,探究了常用的3种k-ε湍流模型和2阶、3阶迎风格式对激波捕捉效果的影响。采用标准k-ε双方程模型和2阶迎风格式对不同高度和不同来流马赫数下的偏二甲肼/四氧化二氮(UDMH/NTO)火箭发动机尾焰流场进行喷管-尾焰流场一体化仿真。复燃反应采用12组分18步化学反应模型,喷管入口参数由热力计算给出。结果表明:随着来流马赫数的增加,波节数逐渐降低;随着飞行高度的上升,尾焰影响区域逐渐扩大;复燃反应造成混合区中的O2、N2质量分数大幅下降,而混合区中O、OH、NO质量分数则有所上升。      

外侧面燃烧固体燃料冲压发动机燃烧室流场的数值研究

为研究外侧面燃烧固体燃料冲压发动机燃烧室流场及燃面退移速率的特点,在Fluent平台上完成了内孔、外侧面燃烧SFRJ的燃烧室内燃烧流场的数值计算。在所涉及的工况中,计算结果表明:外侧面燃烧的SFRJ中,再附着点之前,燃面退移速率较大,来流空气质量流率150g/s,总温600 K时,最大燃面退移速率比内孔燃烧增大43.5%;再附着点之后,燃面退移速率快速减小;随着来流空气总温的减小,固体燃料末端的燃面退移速率开始沿轴向增大;随着来流空气质量流率的增大,燃面退移速率开始增大的位置不断前移,而其增大速率不断减小;因大部分区域内燃面退移速率较小,导致其平均燃面退移速率比内孔燃烧减小21.9%至40.5%;外侧面燃烧的推力比内孔燃烧的小,但比冲相差不大;补燃室及喷管表面处流场温度比内孔燃烧低500~1000 K。     

固液火箭发动机工作过程三维数值仿真

根据固体燃料壁面与气相间的流固耦合得出了固体燃料燃速模型,对采用星形装药的H₂O₂/HTPB(hydroxyl-terminated polybutadiene)固液火箭发动机进行了燃烧流动三维数值仿真,得到了流场参数的分布及不同位置的固体燃料燃速,与二维轴对称仿真结果进行了对比.计算结果表明:装药截面的火焰层形状与装药星孔型面形状相似,但火焰层厚度与位置在星根与星角处存在差异;随着轴向位置的增加,氧化剂不断消耗,火焰层向通道中心移动;固体燃料燃速与氧化剂流率及不同装药位置有关,其大小随氧化剂流率的增加而增加,星根处燃速比星角处大;在相同氧化剂流率下,三维星形装药比二维轴对称装药的平均固体燃料燃速大.      

不同半径圆柱诱导CH4/空气预混燃烧的计算研究

本文研究不同半径圆柱诱导CH4/空气预混燃烧流场。采用保自由流5阶WENO格式求解贴体坐标变换后的多组分Euler方程,用基元反应模型描述CH4/空气燃烧。不同于标准WENO格式通量构造方法,该WENO格式数值通量由方程的解进行WENO插值得到,在曲线坐标系下具有保自由流性质。首先给出了保自由流WENO格式精度和保自由流的数值验证,然后计算圆柱诱导CH4/空气预混气燃烧流场,并考察不同半径圆柱的影响,给出燃烧流场压力与温度等值线和云图、压力和温度沿过驻点线分布。结果表明:在考核计算结果网格无关性基础上,该WENO格式可准确地捕捉激波和火焰阵面形状、激波和火焰面驻点距离,得到的计算结果和文献结果相符。增大圆柱半径,激波和火焰面被推向来流方向,激波和火焰面之间距离也减小。和TVD格式相比,5阶WENO格式采用四分之一的网格数可得到近似相同的计算结果。     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切。使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格。应用特征恢复技术恢复凹面特征并改善凸面处的网格质量。流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量。加入Spalart-Allmaras一方程湍流模型以获得湍流解。上下对称高斯 赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性。算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算。     

辅助动力装置环形回流燃烧室数值研究

利用Fluent软件计算辅助动力装置(auxiliary power unit,简称APU)环形回流燃烧室三维两相喷雾燃烧流场,研究不同进口温度和油气比对两相喷雾燃烧流场的影响,采用标准 k-ε 模型模拟湍流黏性,离散相模型(DPM)追踪油珠运动轨迹,燃烧模型采用非预混平衡化学反应模型.计算结果表明:随着进口温度和油气比的增加,燃烧室出口温度相应增加,但温度分布规律基本保持不变;计算结果与试验测量结果比较吻合,说明采用的数学模型和计算方法可用于预估实际APU回流燃烧室三维两相喷雾燃烧流场.      

氢气超音速燃烧非定常过程数值模拟

用有限差分法求解二维Ｎ－Ｓ方程和组分守恒方程，使用代数涡粘性湍流模型，两步化学反应模型和ＭａｃＣｏｒｍａｃｋ时间分裂法，模拟超音速气流中横喷氢气自动着火及火焰传播过程，清晰地描述了喷氢、火焰传播到主流、稳定燃烧过程中，温度、压力、速度和油气当量比的变化。计算结果表明，高温附面层和喷嘴前回流区能促使氢气自动着火，进口气流的静温是火焰向主流传播的重要影响因素。计算还证明了喷嘴后氢与空气两股流呈现平行流动，其间存在着强烈的湍流交换和化学反应边界层，超音速燃烧效率则主要取决于两股气流的混合速率。     

基于火焰面波动机理的热释放模型理论推导及其模态分析

针对一个中心体稳定的燃烧室通过G方程的方法推导建立了火焰面波动与声波速度波动之间的耦合关系,并进而得到了燃烧热释放模型的数学函数.推导过程中考虑了燃烧膨胀效应对火焰面运动的影响,弥补了前人研究中忽略了燃烧膨胀效应的不足.该函数表明在几何结构和燃料成分等参数保持不变的情况下,湍流燃烧速度是影响热声振荡的关键参数,能够同时强烈的影响燃烧热释放模型的振幅和相位,当量比和入口空气温度等参数是通过改变湍流燃烧速度最终影响了热声振荡现象;而只有当燃烧火焰传播速度相对较大时,密度膨胀比才会对热释放幅值产生影响,并且会对热释放的相位产生影响,从而指出了各参数影响热声振荡的具体途径和相对大小.通过将推导的燃烧热释放波动模型与模态分析方法相结合,系统连续地分析了压力振幅增长率、固有频率和声压模态等重要的热声振荡特性参数随着当量比和入口空气温度的变化关系,从而得到了热声振荡稳定区间和非稳定的区间的云图.对理解和控制热声振荡现象具有重要的实际指导意义.      

采用flamelet模型分析旋流杯燃烧室近贫油熄火点

针对进口温度388 K不变,不同进气速度下对旋流杯燃烧室常压模化实验得到的贫油熄火油气比数据,采用带有详细化学机理的非预混flamelet燃烧模型进行数值模拟,对近贫油熄火点的局部热态流场进行了分析.结果表明:近贫油熄火点时刻,随着油气比的下降,回流涡旋中心和高温度区分离,回流区内气流温度下降,着火点后移,点火距离增加;进气速度的增加和燃油流量的下降,使得回流区回流量增加,回流涡旋中心两侧温差加大,而着火点更靠后,增加了对这种变化的耐受性.      

固体燃料冲压发动机燃速理论分析与数值模拟

从理论和数值模拟两方面探讨了固体燃料冲压发动机 (SFRJ)燃速影响因素, 给出了平均燃速与入口空气流量、空气温度以及燃烧室压强之间的指数关系式, 同时考虑了燃速沿轴向的变化。与实验数据的比较表明, 本文提出的模拟方法是正确而且可靠的。所得结果对于固体燃料冲压发动机设计和实验均具有一定指导意义。      

支板火箭引射冲压发动机引射模态燃烧流动 (Ⅱ )二次燃烧及构型的影响

对比研究了两种构型的支板火箭引射冲压发动机引射模态的瞬时掺混燃烧(SMC)三维掺混和反应流场，详细分析了静态海平面零马赫数情况下燃烧及构型对引射流场的影响，发现几何构型和二次燃烧的综合影响决定了引射掺混后流体的速度、总温、总压及引射流量，从而也主要确定了发动机的性能，其中构型因素主要决定了掺混的质量，从而决定了低速模态的性能。     

LPP低污染燃烧室单头部燃烧性能试验

对贫油预混预蒸发(LPP)低污染燃烧室单头部三级旋流器进行燃烧性能试验,研究不同的油气比、进口空气流量和进口空气温度以及值班级喷嘴安装位置对燃烧室出口截面燃烧性能的影响,获得了燃烧室出口截面温度分布、燃烧效率以及污染物排放的规律.试验结果表明:①油气比增加,NOx排放相应增加;头部A燃烧性能稍优于头部B;②同一油气比下进口空气温度越高,其燃烧污染物排放越多;进口空气流量越大,污染物排放越少;③值班级喷嘴安装位置对LPP低污染燃烧室燃烧性能有一定影响.      

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

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

采用详细化学反应机理的火焰面模型模拟煤油两相燃烧流场

基于可实现的 k-ε湍流模型、颗粒随机轨道模型、火焰面模型和航空煤油详细化学反应机理对模型燃烧室内两相燃烧流场进行了数值模拟.其中详细反应机理由替代燃油(80%质量分数的正癸烷,20%质量分数的1,2,4-三甲基苯)的反应机理和NO_x的反应机理组合而成.通过与实验数据的比较,考察采用该详细化学反应机理的火焰面模型模拟RP-3航空煤油燃烧流场的准确性,特别是污染物排放计算的精度.结果表明:稳态火焰面模型模拟的温度场和CO₂生成量与实验数据吻合较好,而预测的NO排放量与实验值偏差较大;非稳态火焰面模型显著提高了NO的预测精度,在工况1(来流马赫数为0.16,进口温度为537K,油气比为0.0048,常压)条件下与实验数据吻合较好,但在工况2 (来流马赫数为0.155,进口温度为523K,油气比为0.010,常压)条件下仍过高估计了NO的排放量.