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

基于可实现的 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的排放量.      

超声速横侧射流燃烧中可压缩性影响的数值研究

为了研究湍流的可压缩性效应对超声速燃烧的影响,对复杂的超声速横侧射流氢气/空气燃烧流场进行了数值模拟。采用修正的k-ε湍流模型(包括膨胀可压缩性和激波不稳定性两部分修正),考虑H2/Air详细化学反应机理(GRI-Mech2.11机理,10组分,28基元反应)。结果表明,可压缩性修正k-ε湍流模型能够预测复杂的超声速横侧射流流场结构。可压缩性使流体分离增强,湍流混合受到抑制,从而显著影响了超声速燃烧过程,数值模拟中需要考虑可压缩性效应的影响。     

应用修正的k-ε模型研究超声速H2/Air燃烧

为了研究超声速燃烧中流体可压缩性的影响,对标准k-ε湍流模型进行可压缩性修正(包括结构可压缩性修正和膨胀可压缩性修正两部分)。分别应用标准k-ε模型、修正的k-ε模型和雷诺应力模型(RSM),考虑氢气/空气详细化学反应机理(GR I-M ech 2.11机理,10组分,28基元反应),数值模拟有壁面限制的超声速混合层冷态及热态流场。结果表明:壁面和燃烧对湍流影响都很大;修正模型对冷态以及燃烧场的预测结果优于其它两个;修正模型预测的混合层厚度更薄,燃烧区域更窄,与实验结果吻合地更好。     

钝体稳定射流火焰的火焰面模型数值研究

基于稳态火焰面(SLFM)和交互式非稳态欧拉颗粒火焰面(EPFM)模型对Sydney大学CH4/H2钝体稳定湍流扩散火焰行了数值研究,采用修正的雷诺应力模型(RSM),同时对两种不同规模的甲烷详细化学反应动力学机理进行研究,比较了燃烧模型和反应机理对湍流火焰结构、活性自由基以及氮氧化物预测精度的影响,与实验数据对比结果表明:两种反应机理得到的温度场和主要组分分布基本相同;SLFM模型能对速度场和标量场的分布进行较为准确预测,采用EPFM模型修正后,部分区域OH预测结果更加靠近实验结果;采用EPFM模型对SLFM模型耦合GRI-Mech 211的计算结果修正后,NO量级降低近2倍,预测精度明显改善,与实验结果实现较好的符合,验证了化学反应动力学机理以及非稳态效应对氮氧化物预测精度的影响。     

自适应湍流耦合建表燃烧模型的振荡燃烧数值模拟

航空发动机及燃气轮机等动力装备燃烧室广泛采用的贫燃燃烧方式经常遇到破坏性的非定常热声耦合振荡燃烧问题。非定常振荡燃烧数值预测是一个长期的研究热点和难题。发展了针对振荡燃烧的耦合直接求解数值模拟方法，包括优化的动态模型参数的高精度自适应湍流模型（SATES），耦合可压缩的详细化学反应建表FGM燃烧模型。选取的3种湍流燃烧模型包括有限速率模型（W1）及火焰面密度封闭方法中的Zimont(W2)和Fureby(W3)2种褶皱因子模型。针对经典的LIMOUSINE燃烧室多个部分预混振荡燃烧工况开展了数值研究，发现自适应湍流模拟框架下的3种燃烧模型均准确预测到了振荡燃烧的振荡频率，与试验相比，误差<6.4%；对于振荡燃烧压力脉动振幅的预测结果，有限速率模型（W1）和Zimont(W2)模型结果显著大于试验值，误差>380%;Fureby(W3)模型结果与试验值吻合较好，误差<17.9%。表明振荡燃烧的数值预测对不同的湍流及湍流燃烧模型具有较高的敏感性。不同的工况结果表明，振荡燃烧存在完全振荡模态和过渡模态，完全振荡模态中数值预测的特征主频在燃烧室上下游多个位置趋于一致；过渡振荡...     

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

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

航空发动机燃烧室燃烧过程与排放物生成的反应动力学数值模拟

选用正癸烷作为航空煤油的替代燃料,建立了正癸烷的化学反应详细机理与简化机理(包括50种组分,118个基元反应).分别采用详细机理与简化机理对正癸烷在激波管中的着火延迟时间、在预混燃烧炉内的燃烧过程进行了数值计算,并与实验结果进行了对比分析.同时,耦合该简化机理与CFD计算软件Fluent,对某型航空发动机环管形燃烧室中单个火焰筒内流动特性与燃烧过程、排放物及活性中间组分生成的反应动力学特性进行了详细分析,并与采用C₁₂H₂₃为燃料的单步反应机理的计算结果进行了对比分析.结果表明:采用简化机理计算得到的着火延迟时间、反应物与各主要生成物摩尔分数的整体变化趋势与实验数据吻合较好;与采用C₁₂H₂₃为燃料的单步反应机理相比,采用正癸烷为替代燃料的简化反应机理计算得到的温度场分布更符合实际,其出口平均温度亦更为接近燃烧室出口设计温度;同时,能更为详细了解燃料低温裂解过程及裂解产物、中间产物及主要排放物的生成规律.      

输运概率密度函数中的小尺度标量混合建模

输运概率密度函数方法能精确地求解湍流燃烧中的有限化学反应速率,但小尺度上分子扩散引起的组分变化则需要通过标量小尺度混合模型来模拟。论文综述了当前湍流预混燃烧中混合模型研究的最新进展。首先介绍了现有标量混合模型在不同预混燃烧模式下的表现:对于火焰面燃烧模式,保持组分空间的临近性对于标量混合模型非常重要,标量混合频率受火焰结构影响显著;对于破碎反应区模式,组分空间保持临近性的重要性有所减弱,标量混合频率由湍流主导。接着介绍了最新基于直接数值模拟数据的标量小尺度混合机制和建模研究,包括湍流预混燃烧中反应标量小尺度耗散机制、反应标量混合时间尺度在不同燃烧模式下的演化规律、依据线性混合思路构建的标量混合时间尺度模型。针对不同湍流预混燃烧模式的模型验证表明:相比于已有湍流混合时间尺度模型,新模型显著提升了对湍流预混火焰中的标量耗散率和燃烧特性的预测精度。最后介绍了大涡模拟/输运概率密度函数方法在近极限湍流预混火焰,如甲烷预混值班火焰和湍流对冲预混火焰中的应用,验证了其对贫燃、高湍流度条件下的湍流-化学反应相互作用及局部熄火/再燃现象的预测能力。     

用准稳态方法建立碳氢燃料点火燃烧的简化化学反应动力学模型

基于“准稳态”方法建立了一套碳氢燃料点火燃烧的化学反应动力学模型简化方法和相应的软件SPARCK,并从甲烷点火燃烧的GRI2.11详细基元反应动力学模型出发简化得出了包含16个组分12步总包反应形式的简化化学反应动力学模型,从庚烷点火燃烧的详细基元反应动力学模型出发简化得出了包含25个组分21步总包反应形式的简化化学反应动力学模型。通过其计算结果与CARM软件导出简化模型的计算结果和典型激波管试验结果的对比可以看出,本文简化得到的简化反应动力学模型能较为有效地再现详细基元反应模型的反应机理,简化模型的计算精度与CARM软件导出简化模型的计算精度相当。与详细基元反应动力学模型相比,简化模型有效地减少了反应组分,在工程计算中有比较好的应用前景。     

Hyshot超燃冲压发动机的CFD数值模拟

以Hyshot超燃冲压发动机试验为研究对象,分别对冷流和燃烧工况进行了CFD数值模拟.冷流模拟计算得到的壁面压力与试验结果吻合良好,且数值计算压力分布对采用的湍流流动模型不敏感;但是采用不同的湍流流动模型计算的湍流参数(湍动能和耗散率)差别很大,会对氢气燃料与气流的掺混产生重要影响,进而影响燃烧模拟结果.采用SST(shear stress transport)k-ω湍流流动模型和EDM(eddy dissipation model)湍流燃烧模型得到的模拟结果与试验结果基本吻合,但是计算压力峰值略靠后.对8个不同氢气燃料喷射角度的工况进行数值模拟结果表明喷射角度为99°和114°时燃料和空气混合最好、燃烧效率可达到最高.     

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

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

大涡模拟模型环形燃烧室污染特性

在任意曲线坐标系下,采用大涡模拟模型环形燃烧室两相燃烧流场中污染物的生成.用亚网格涡破碎(EBU)燃烧模型估算化学反应速率;分别用NO亚网格动力学模型与CO亚网格模型来模拟NO和CO的生成;采用随机离散模型模拟气液两相湍流流动,两相之间耦合采用PSIC算法.计算结果与瞬态速度场、出口温度和污染物质量分数分布的实验数据比较吻合,表明大涡模拟方法和亚网格模型可以用来预估实际燃烧室两相喷雾燃烧流场和污染物生成过程.      

含硼贫氧燃气补燃室反应流场研究

以King模型模化硼粒子燃烧,以几率密度模型模化湍流燃烧,对含硼贫氧推进剂固体火箭冲压发动机补燃室流场进行数值研究,并与试验结果进行了对比。计算表明:环向燃料喷口布局可增加贫氧燃气气相组分与空气的混合、反应效率;突扩形的几何结构有利于燃气产生大涡运动。     

高雷诺数槽道湍流的壁面模化大涡模拟研究

选取基准壁湍流的槽道流动,研究了多种模型的壁面模化大涡模拟.模型包括经典的大涡模拟、Spalart-Allmaras、分离涡模拟和一种动态混合模型.基于摩擦速度的雷诺数范围为395~12000,采用3组粗糙网格,流向和展向维数分别同取37,49和65,法向维数保证y+(1)~1.主要研究平均速度、雷诺切应力分布、详细分析了各模型的特性差异并展示了相应的湍流结构.研究表明:在高雷诺数粗糙网格下,大涡模拟失去求解精度,分离涡模拟出现对数律不匹配,动态混合模型的计算接近直接数值模拟,其对数率区可解应力约占雷诺切应力的93%,边界层外层可解应力约占99%.这说明合适的混合模型可以在经济成本下保证计算精度,具有解决实际问题的潜力.      

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

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

一种改进的类DES湍流模拟方法

构造了一种基于一方程S-A(Spalart-Allmaras)模型和一方程Yoshizawa亚格子模型的混合RANS/LES(Reynolds-averaged Navier-Stokes/large eddy simulation)湍流模拟方法.在涡黏假设的基础上,将Yoshizawa亚格子湍动能方程转化为等效的亚格子湍流涡黏性输运方程,并采用混合函数将其与S-A模型方程进行混合,从而改进了DES(detached eddy simulation)模型的亚格子行为,同时克服了其依靠网格控制模型转换的缺点.模拟了超声速的带斜坡凹腔流动,并与相同网格下的LES及DES结果进行了比较,结果表明该混合RANS/LES方法在远离壁面的自由剪切流区域与LES行为一致,而在附着边界层区域表现优于LES和DES方法.      

湍流边界层各向异性张量模式的实验研究

使用高时间分辨率粒子图像测速技术在水槽内精细测量平板湍流边界层内不同流向、法向位置的瞬时速度矢量空间场的时间序列信号。利用实验数据对湍流各向异性的涡黏张量模式进行了实验研究,给出了湍流边界层不可压缩槽道湍流中各向异性的涡黏张量各分量的空间分布规律。各向异性的涡黏张量分量的空间分布揭示了流向平均速度的法向梯度对雷诺应力三个分量的不同贡献。涡黏张量三个分量空间上的不同规律说明对涡黏张量的各向异性的考虑是合理的,而各向同性的涡黏模型不适宜描述剪切湍流雷诺应力的物理本质。     

平板射流中大尺度涡的理论模型

根据流动稳定性理论,提出了一种解释平板射流中大尺度涡产生机理的二维不稳定波理论模型,用这种模型计算所得的流动结构和相干扰动量与实验结果吻合较好.说明不稳定波模型可以较好地描述湍流平板射流中的大尺度涡结构的主要特征.      

模型燃烧室紊流燃烧的大涡模拟

采用了两种不同的亚网格尺度燃烧模型对带 V型稳定器的模型燃烧室紊流化学反应流动进行了大涡模拟 ,用 k-ε方程亚网格尺度模型确定亚网格紊流粘性 ,为了考虑热辐射对燃烧室壁温和气流温度的影响 ,运用热流法辐射模型估算热辐射通量 ,用 SIMPLE算法和混合差分求解大涡模拟各守恒方程 ,通过对两种不同亚网格尺度燃烧模型数值模拟结果与实验值的比较表明 ,两种燃烧模型都与实验值较吻合 ,但 G方程小火焰模型要比亚网格 EBU燃烧模型符合得更好些      

湍流边界层相干结构空间拓扑形态的 层析TRPIV测量

用层析TRPIV(time-resolved paticle image velocimetry)技术精细测量了水洞中平板湍流边界层三分量速度的时空序列信号,提出了空间局部平均多尺度速度结构函数的新概念描述湍流多尺度涡结构的空间拉伸、压缩、剪切变形和旋转.用空间局部平均多尺度速度结构函数对湍流脉动速度进行了空间多尺度分解.用空间局部平均多尺度速度结构函数的新概念,根据湍流多尺度涡结构在空间流向的拉伸和压缩特征,提出了新的湍流相干结构条件采样方法,检测并提取了层析TRPIV数据中相干结构的“喷射”和“扫掠”事件中的速度、涡量等物理量的空间拓扑形态.发现在喷射和扫掠事件中均存在一对反向旋转的准流向“马蹄形”涡结构.