回流燃烧室燃烧过程的三维数值模拟

在三维任意曲线坐标系下数值模拟回流燃烧室火焰筒内两相燃烧过程,采用RNG k-ε模型模拟紊流粘性,EBU-Arrhenius模型模拟燃烧反应速率、离散坐标模型以及六通量模型考虑辐射传热,液相采用颗粒轨道模型,气相采用SIMPLE算法求解,并用PSIC算法考虑气液两相之间的相互作用的影响,计算得到燃烧室内速度、温度等各气流参数分布.通过将计算与实验结果对比表明,计算方法可靠,离散坐标模型优于六通量辐射模型,更适用于模拟火焰筒内两相燃烧流场.      

RNG k-ε模型数值模拟油雾燃烧流场

采用RNG(R enorm alization G roup)k-ε紊流模型对某环形燃烧室火焰筒内三维两相燃烧流场进行数值模拟。运用偏微分方程法和区域法生成贴体网格,EBU-A rrhen ius紊流燃烧模型估算化学反应速率,六通量热辐射模型估算辐射通量。在非交错网格下求解气相采用S IM PLE算法,液相采用颗粒轨道模型和PS IC算法。近壁区处理分别采用两层和三层壁面函数。计算数据与试验值比较表明,采用三层壁面函数的RNG k-ε模型更适用于模拟三维两相燃烧流场。      

数值研究环形回流燃烧室紊流燃烧流场

在任意曲线坐标系下,采用四步反应机理数值研究回流环形燃烧室内燃烧过程中间及最终产物的生成.采用重整化群k-ε模型、涡团耗散模型及离散坐标法分别模拟紊流黏性、燃烧反应速率及辐射传热;采用随机颗粒轨道模型描述油珠温度、运动及其尺寸变化等过程;采用PISO(pressure implicits plit-oper-ator)算法对离散方程进行求解;在非交错网格体系下,预估该燃烧室三维两相燃烧流场以及燃烧产物分布.计算与实验结果较为相符,表明本文的数值方法及模型可靠,可为环形回流燃烧室的优化设计及污染物减排提供技术支持.      

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

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

燃烧室贫油熄火极限数值预测

在任意曲线坐标系下分别对双级轴向涡流器矩形燃烧室及涡流杯回流环形燃烧室三维两相反应流场进行数值模拟,采用RNG(renormalization group)k-ε,SOM-EBU(second order moment-eddy break up model)及六通量辐射等模型考虑紊流流动、化学反应及辐射传热对燃烧过程的影响;采用SI MPLE(semi-implicit method for pressure linked equation)算法求解流场,气液两相耦合采用颗粒源项法处理.采用提出的燃油稳态逐次逼近法对贫油熄火极限进行数值预测,所得计算值与实验数据较为相符,表明此数值预测方法可靠,可用于工程实践.      

数值分析二级涡流器环形燃烧室的燃烧性能

采用多维经验分析法在任意曲线坐标系下对包括二级突扩扩压器、双级轴向涡流器及火焰筒在内的单头部环形燃烧室的燃烧性能进行计算.并对RNGk-ε紊流模型加以改进, 使其更适用于数值模拟非交错贴体坐标系下的三维两相化学反应流.采用EBU-Arrhenius紊流燃烧模型、六通量热辐射模型以及颗粒群轨道模型对油雾燃烧进行模拟.同时, 数值分析了不同涡流器几何尺寸对燃烧室整体流场和燃烧性能的影响, 数值计算结果与实验数据相当一致, 表明计算方法合理, 计算程序可靠, 可用来估算环形燃烧室的燃烧性能.      

小发直流燃烧室扩压器和火焰筒匹配的数值模拟

采用数值模拟方法对某小发直流燃烧室扩压器和火焰筒匹配下的流场进行了计算和分析,用双方程k-ε模型描述紊流特性,用二步化学反应模型模拟化学反应,用空气流量分配作为火焰筒内外流场联合求解的桥梁,引入边界结点类型说明数组帮助复杂边界条件的处理。计算结果表明燃烧室环腔间形成两个强烈的回流区,内、外环腔流场极不对称,并且具有较强的三维性;前置扩压器倾斜角度是影响燃烧室流场和性能的最重要因素。      

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

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

R0110重型燃气轮机燃烧室三维数值模拟

采用SIMPLE算法,应用带有旋流修正的k-ε双方程湍流模型及有限速率／涡耗散化学反应模型,对R011O重型燃气轮机逆流环管型燃烧室单个火焰筒进行了三维数值模拟计算。将计算出的燃烧室出口温度场的分布、品质及火焰筒壁温与试验结果进行了对比分析。燃烧室进口流量、温度、压力等气动参数均与试验时保持一致,火焰筒各部分空气流量也均按火焰筒空气流量分配试验结果给定。计算和对比分析的结果表明,计算得到的燃烧室出口温度场的分布、品质及火焰筒壁温分布与试验结果比较接近。     

显式代数雷诺应力模型在超声速/高超声速流场计算中的应用

为了提高显式代数雷诺应力模型在高超声速流场中的计算精度,以修正后的CG k-ε模型作为湍流时间尺度决定方程,将Wallin和Johansson显式代数雷诺应力模型(WJ2000模型)应用于超声速/高超声速流场数值计算,并与线性湍流模型的计算结果进行了对比分析。结果表明WJ2000模型对高超声速压缩拐角附近壁面压力和热流分布的预测要优于线性湍流模型,同时网格依赖性明显低于后者;对CG k-ε湍流时间尺度决定方程所做修正不影响WJ2000模型对超声速流场的预测结果,但使其对高超声速流场的预测精度明显提高。     

高超声复杂流动中湍流模式应用的评估

本文选择几个高超声速基准流动：二维可压缩拐角、锥柱裙组合体绕流、斜激波与湍流边界层干扰，采用几个常见的湍流模式，BL模式，CHκ-ε模式，κ-ε模式，SST模式，CMOTTκ-ε模式，SHIHκ-ε模式，通过将数值计算结果和实验结果进行比较，对有关的湍流模式地 评估，得到一些有意义的结论。     

不可压缩湍流的多尺度模型

基于可变时间间隔平均方法,提出了一种不可压缩湍流多尺度模型及平均流动方程.与传统的单尺度湍流模型不同,该模型在建立的过程中,保留了湍流的多尺度特性,结合平均流动方程,可以更好地预测湍流流场特征.通过模拟平面后台阶流动和不对称平板扩压器流动,并将预测结果与标准k-ε模型的预测结果对比,初步验证了其可信性及优越性.结果表明:计算所得的平面后台阶流动的流向再附长度与台阶边压力系数比比标准k-ε模型提高精度约20%;平板扩压器流动的回流区位置误差约为7%、倾斜壁面摩擦因数误差约为5%,而标准k-ε模型未能预测出分离现象.可以看出该模型适用于典型的分离流动,在湍流流场的预测中表现优异,具有一定的工程应用价值      

Sajben扩压器复杂流动中湍流模型的性能评估

用8个常用的湍流模型对Sajben扩压器中跨声速流动进行了数值模拟,评估了Spalart-Allmaras, 标准k-ε, RNG (re-normalization group) k-ε,realizable k-ε,标准k-ω,SST(shear stress transport) k-ω,v2-f,Reynolds stress共8个湍流模型对激波/湍流边界层相互作用的模拟预测能力.通过与实验数据比较发现:SST k-ω模型和v2-f模型比其他模型模拟的更准确,其中SST k-ω模型比v2-f更能准确地预测壁面压力,然而对于分离点、再附点以及分离区长度v2-f比SST k-ω预测得更准确.      

高超声速飞行器-进气道一体化热流数值计算

采用CFD(计算流体动力学)技术, 开展了飞行器前体/发动机一体化气动热环境分析.对层流区、转捩区和湍流区分别采用计算模型, 在湍流区利用压缩性修正的SSGZ-Jk-ε湍流模型, 在转捩区引入代数型转捩因子模型描述边界层由层流逐渐过渡为完全湍流的流动过程.计算了前体和内通道的表面热流, 并与实验结果进行了对比.结果表明所采用的计算方法可以较好地预测前体及发动机内通道热流率, 流动状态、几何结构及激波入射对热流值影响较大.      

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

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

不可压缩湍流的色散模型

分析了广泛存在于湍流运动中的能量逆转现象,揭示了其产生的原因,提出了湍流的色散性质.在此基础上,修改了Boussinesq假设,建立了包含色散效应的新的雷诺应力封闭式和湍流色散系数,给出与不同模型相结合时,湍流色散系数所具有的不同形式,并阐述了湍流运动中能量的传递方向及条件.采用不可压缩平板边界层流动和平面后台阶流动验证了其可信性和优越性.平板摩擦阻力系数及边界层速度型与实验结果吻合良好,平面后台阶流动的流向再附长度、台阶边压力系数及湍流强度等参数均比标准k-ε模型更接近实验结果.结果表明:色散项的加入可以在不显著增加计算量的同时显著改善预测精度,模型具有一定的工程应用价值.      

单列叶栅紊流流场有限分析数值研究

本文将有限分析方法用于曲线座标系上紊流 N-S方程数值计算,研究了高雷诺数时叶栅粘性紊流流场。有限分析方法在网格单元上对非线性偏微分程进行线化处理,在解析边界条件下求出线化偏微分方程的解析解,以此解构造离散方程。有限分析方法能够根据对流的方向和大小自动改变格式系数,并具有数值扩散小、精度高和稳定性好等优点。本文以 k-ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。      

封闭转静盘腔中紊流的数值模拟

对一个旋转盘腔的流场进行了数值计算研究。该盘腔由一个旋转盘,一个静止盘及静止的外围屏组成。在盘腔的核心区采用标准k-ε模型,在近壁区分别采用了两种低雷诺数紊流模型,即(1)单方程模型和(2)Launder-Sharma低雷诺数k-ε模型,并且在第二个模型中考虑了固体旋转引起的耗散率降低的修正。