两种亚网格燃烧模型的旋流扩散燃烧大涡模拟

用二阶矩亚网格（SOM—SGS）燃烧模型和文献中的涡旋破碎亚网格（EBU—SGS）燃烧模型，对甲烷一空气旋流扩散燃烧进行了大涡模拟，将二者得到的LES统计平均温度分布和实验结果以及用二阶矩燃烧模型的统观模拟（RANS—SOM）结果比较，表明LES—SOM和RANS—SOM的模拟结果都和实验符合较好，而LES—EBU的模拟结果和实验不符合，在不同区域内高估或者低估了燃烧温度。其原因是由于EBU模型不能有效地考虑有限反应动力学的作用。LES—SOM模拟的瞬态结果显示了旋流扩散火焰的湍流大涡结构不同于射流火焰的特点。     

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

对带双级扩压器的模型燃烧室气液两相瞬态喷雾燃烧过程,在三维贴体坐标系下采用欧拉-拉格朗日两相大涡模拟方法进行数值研究,同时采用多维经验分析法预估燃烧性能.采用 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算法.计算结果与瞬态速度场、出口温度和污染物质量分数分布的实验数据比较吻合,表明大涡模拟方法和亚网格模型可以用来预估实际燃烧室两相喷雾燃烧流场和污染物生成过程.      

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

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

五种湍流涡粘模型在二维方柱绕流数值模拟中的对比研究

为研究不同雷诺时均Navier-Stokes(RANS)模型在求解钝体绕流场的差异,采用五种两方程湍流涡粘模型计算了二维方柱的绕流场(雷诺数Re=22 000),得到了不同湍流模型在不同网格离散方案下的方柱气动力与流场特性。结合以往文献数据,通过对比不同湍流模型计算结果的差异,揭示了不同湍流模型的特点。研究结果表明:网格离散方案对方柱绕流数值模拟结果有重要影响;RNGk-ε湍流模型的计算效率最高,SSTk-ω湍流模型的计算效率最低;Standard k-ε湍流模型的计算结果准确程度整体弱于其余湍流模型;RNGk-ε湍流模型、Realizable k-ε湍流模型与Standard k-ω湍流模型的计算结果大致相当,较接近大涡模拟结果;SSTk-ω湍流模型的模拟结果优于其余湍流模型,其尾流区速度场与试验结果吻合较好;两方程湍流模型计算的二维方柱绕流结果与试验结果以及大涡模拟结果相比,存在不可忽略的差异。     

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

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

三维槽道湍流的大涡模拟

利用大涡模拟方法模拟了雷诺数为13000的充分发展槽道湍流流动。为了比较网格和亚格子模型对计算结果的影响,采用了三种不同的网格数目,以及Smagorinsky模型和动力模型等两种亚格子模型,并与直接数值模拟结果进行了比较。结果表明,大涡模拟可以比较准确地给出速度分布、脉动速度均方根分布等湍流统计参数。在此基础上,对近壁区小尺度湍流结构进行了研究。     

湍流模型对三维翼梢涡流场数值模拟的影响

 梢涡流场数值模拟对减少梢涡的不良影响意义深远。应用雷诺平均Navier-Stokes(RANS)方法建立了三维翼梢涡流场数值模型,选用的湍流模型包括两种经典两方程模型和一种代数雷诺应力模型。计算内容包括翼梢涡流场速度分布、压力分布及涡核位置,为了准确模拟梢涡涡核内系统旋转和流线曲率影响,将旋转和曲率修正方法应用于选用的湍流模型中。计算时对翼梢涡涡核处网格进行了加密处理。将计算结果与实验结果对比后发现,修正后的湍流模型较标准湍流模型明显提高了预报精度,与实验结果吻合较好,具有工程应用价值。     

环形燃烧室两相喷雾燃烧的大涡模拟

在三维任意曲线坐标系下采用欧拉-拉格朗日两相大涡模拟方法对环形燃烧室火焰筒气液两相紊流瞬态反应流进行数值模拟.采用椭圆偏微分方程生成三维贴体网格,计算中所采用的数学度模型有:k方程亚网格尺度模型估算亚网格紊流粘性;亚网格EBU燃烧模型估算化学反应速率;热通量辐射模型估算辐射换热.并在非交错网格体系下气相采用SIMPLE算法和混合差分格式求解,液相采用随机离散模型(Stachasttc Separated Flow,简称SSF),在拉格朗日坐标系下追踪各油珠群沿各自轨道运动、质量损失及能量变化.通过计算结果与实验数据相比较,表明在三维贴体坐标系下对燃烧室火焰筒两相紊流油雾燃烧流场进行大涡模拟,采用欧拉-拉格朗日两相大涡模拟方法能反映两相紊流化学反应流流动及实际燃烧过程.     

Large eddy simulation of hydrogen/air scramjet combustion using tabulated thermo-chemistry approach

Large eddy simulations(LES) have been performed to investigate the flow and combustion fields in the scramjet of the German Aerospace Center(DLR). Turbulent combustion is modeled by the tabulated thermo-chemistry approach in combination with the presumed probability density function(PDF). A b-function is used to model the distribution of the mixture fraction, while two different PDFs, d-function(Model I) and b-function(Model II), are applied to model the reaction progress. Temperature is obtained by solving filtered energy transport equation and the reaction rate of the progress variable is rescaled by pressure to consider the effects of compressibility. The adaptive mesh refinement(AMR) technique is used to properly capture shock waves,boundary layers, shear layers and flame structures. Statistical results of temperature and velocity predicted by Model II show better accuracy than that predicted by Model I. The results of scatter points and mixture fraction-conditional variables indicate the significant differences between Model I and Model II. It is concluded that second moment information in the presumed PDF of the reaction progress is very important in the simulation of supersonic combustion. It is also found that an unstable flame with extinction and ignition develops in the shear layers of bluff body and a fuelrich partially premixed flame stabilizes in the central recirculation bubble.     

三维钝体结构风压分布与绕流特性的大涡模拟研究

针对风场作用下三维钝体结构的绕流特性与空间风压分布特性,采用大涡模型对其进行了数值模拟研究,以平均风速和脉动风速作为入口边界条件,在稳态分析的基础上对9种不同来流方向的工况进行了瞬态动力分析,再现了钝体绕流的冲撞、分离、重附着等现象,揭示了钝体绕流主、次涡循环区的形成、发展机理,提出了空间风压场分布的一般规律。研究表明:钝体绕流对流场风速剖面影响的强弱程度依次为尾流区、屋顶区、来流区,且对中心区域流场风速剖面的影响要大于两侧流场,并由中心区域向两侧逐渐减小;来流区钝体绕流特性与空间风压分布特性较屋顶区和尾流区稳定,尾流区涡旋的非对称生成、脱落对钝体左右两侧流场的冲撞、分离、重附着等具有重要的影响;来流方向与钝体呈45°或135°夹角时,对尾流风场的影响范围最小,但随着来流方向逐渐垂直于钝体,尾流影响区域也逐渐增大。     

突扩燃烧室低频燃烧不稳定形成机理分析

突扩燃烧室在一定的工作条件下会出现燃烧不稳定现象。采用实验和数值模拟的方法对突扩燃烧室形成低频燃烧不稳定的机理进行了研究。通过实验研究发现突扩燃烧室压强振动过程中纵向振型占主导地位,但其振动频率并不与声学频率一致。建立了适合分析燃烧不稳定的多步化学反应动力学与大涡模拟耦合的数值分析方法,对实验发动机开展了非稳态数值模拟,获得了低频燃烧不稳定形成演化的详细过程和流场结构。实验和数值计算表明突扩截面形成的旋涡脱落,以及旋涡在燃烧室内的运动过程中引起燃烧面积、局部当量比和热释放率的脉动是激发低频压强振动的主要原因。压强振动引起上游速度脉动,进而形成旋涡脱落。大尺度旋涡在燃烧室内的运动又会引起热释放率的大幅度脉动,反过来又会促进压强振动。振动频率是由压强波和旋涡运动特征时间共同决定的。     

折叠V形火焰稳定器原型的冷态与稳燃特性研究

为了获得折叠V形钝体这一新型火焰稳定器原型的冷、热态特性,采用风洞实验与混合雷诺平均/大涡模拟结合的方法对10组不同折叠角的折叠V形钝体火焰稳定器模型进行了冷态实验与热态数值模拟研究。其中冷态实验风速10~50m/s,来流雷诺数10000～80000,热态数值模拟工况为进气温度700K、来流速度50~150m/s,来流雷诺数20000～60000,当量比0.2~1。通过实验测得了不同折叠角折叠V形钝体总压恢复系数,通过数值模拟进一步获得了阻力系数、吹熄性能与燃烧效率,并获得了折叠V形钝体下游火焰时均轮廓。研究表明,折叠V形钝体总压恢复系数与阻力系数优于标准钝体,总压恢复系数随钝体折叠角减小而增大,阻力系数随折叠角减小而减小;折叠V形钝体抗吹熄能力与燃烧效率均优于标准V形钝体;折叠V形钝体的下游火焰扩散特性与标准V形钝体有显著差异。上述特性揭示了折叠V形钝体具备改善加力燃烧室性能的较大潜力,在实际应用中,还需对稳定器的空间排布做进一步研究与优化。     

后视镜罩边缘结构对流场和气动噪声的影响

采用大涡模拟(LES)的计算方法,对车外后视镜不同边缘结构引起的外部流场和镜后车身表面监测点处的气动噪声进行了数值仿真.研究表明:不同的镜罩边缘结构在较大程度上影响了流经后视镜罩的气流速度和流线方向,对后视镜后部流场和监测点处声压级产生较大影响.相较于原模型的光滑边缘结构,模型1后部流场涡团更远离车身表面,有利于降低气动噪声,某监测点1/3倍频程中心频率处声压级最大降幅接近10dB;模型2后部流场涡团产生分离,并更靠近车身表面,反而使气动噪声增大.模型1的镜罩边缘结构改进方案对车外后视镜-A柱区域的流场和气动噪声都有较好的改善.     

模型燃烧室两相燃烧大涡模拟的并行计算

在三维任意曲线坐标系下采用欧拉-拉格朗日两相大涡模拟方法对模型燃烧室气液两相瞬态喷雾燃烧过程进行大涡模拟并行计算研究。通过大涡模拟计算结果与雷诺平均计算结果对比表明大涡模拟方法能更好的模拟流场细观结构。同时分别采用多个处理器和动态内存分配对大涡模拟并行计算程序进行数值模拟,并行计算结果表明:多个处理器之间并行计算结果与PIV测量的瞬态速度场以及出口温度分布实验数据相互吻合,表明采用并行计算和动态内存分配在保证计算结果正确的前提下,大幅度降低计算机内存和大涡模拟计算时间。     

超声速燃烧气动斜坡喷注器研究

为增强超声速燃烧过程的燃料混合,设计了通过喷嘴阵列形成流向涡强化混合的气动斜坡实验件.在直连式实验台上对气动斜坡实验件进行了纹影及油流谱等冷试实验研究.并采用大涡模拟(LES)方法数值模拟了超声速流场中气动斜坡的流动及燃烧特性.结果表明:气动斜坡喷注器的喷嘴阵列设计能有效加速喷流后流向涡的卷起,提高掺混特性;其燃烧过程也能使喷流形成的剪切层上抬,增进混合;由于喷流方向与主流方向夹角较小,总压损失也较物理斜坡、横向喷流等混合增强方式小.      

模型预混燃烧室燃烧不稳定性研究

针对钝体火焰稳定器结构的模型预混燃烧室的燃烧不稳定性现象进行了实验研究和数值模拟.实验中利用动态压力传感器测量了不同当量比时燃烧室内动态压力的频率和振幅.结果表明:随当量比增加,不稳定性频率从251Hz增加到258Hz.不稳定性振幅与大气压力比值从2.7%增加到6.1%.对燃烧室进行了定常和非定常数值模拟以及声学模态分析,得到了周期性旋涡脱落的频率为260Hz和燃烧室系统的前5阶本征声学模态频率.其中第3阶纵向声学模态频率与实验值基本一致,说明维持不稳定性的机理为钝体火焰稳定器后旋涡的周期性对称脱落和燃烧室系统的第3阶纵向声学模态形成了耦合.