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

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

三维贴体坐标系下燃烧室热态流场的大涡模拟

大涡模拟三维贴体坐标系下环形燃烧室火焰简热态紊流瞬态流场。利用椭圆方程方法生成三维贴体网格，计算中采用k方程亚网格尺度模型估算亚网格紊流粘性；亚网格EBU燃烧模型估算化学反应速率；热通量辐射模型估算辐射通量。并在非交错网格系下采用sIMPLE算法和混合差分格式求解离散方程，利用壁面函数处理固壁边界条件。计算结果与实验结果的比较表明，采用大涡模拟方法能更真实反映环形燃烧室火焰简内紊流化学反应流气流结构和燃烧过程。     

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

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

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

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

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

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

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

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

燃烧室污染物生成大涡模拟模型的研究

采用亚网格动力学模型和联合概率密度函数(PDF)亚网格模型对污染物NO的生成进行大涡模拟.计算中采用的数学模型有:k方程亚网格尺度模型来估算亚网格湍流粘性;亚网格EBU燃烧模型估算化学反应速率;在交错网格体系下, 采用SIMPLE算法和混合差分格式来求解差分方程.数值分析两种不同亚网格污染物生成模型和两种进口气流温度对NO排放的影响, 计算结果与实验数据比较表明:两种亚网格污染物生成模型都能较好的预估燃烧室污染物NO的生成, 但亚网格动力学模型的计算工作量要比联合概率密度函数亚网格模型少得多, 更适合工程应用.      

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

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

低热值气体燃料燃烧室数值模拟与试验研究

以某型航空发动机燃烧室为研究对象,对低热值燃料燃烧室进行了数值模拟和试验研究。采用RNGk-ε模型、小火焰紊流燃烧模型和P-1辐射模型,预估了紊流特性、化学反应速率和辐射通量;应用SIMPLE算法,对离散方程进行求解。计算结果与试验数据比较表明二者基本吻合,这说明计算方法合理,可用来估算低热值燃料燃烧室的燃烧性能。     

燃气轮机燃烧室中贴体网格的生成与应用

运用贴体坐标系统分别对燃气轮机燃烧室的冷态和热态流场进行了数值模拟。燃烧室的曲线壁面边界用Ｔｈｏｍｐｓｏｎ的非正交贴体坐标网格处理，并把ＳＩＭＰＬＥ算法应用到曲线坐标系下求解各守恒方程。采用双方程ｋ－ε模型来描述紊流特性；采用按简单的单步化学反应假设的Ａｒｒｈｅｎｉｕｓ－ＥＢＵ模型表述紊流燃烧；采用热通量法辐射模型来估算辐射通量和壁面温度。另外，还研究了冷态时不同旋流数和紊流强度对气流结构的影响。     

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

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

加力室液雾燃烧的模型与应用

提出了描写航空发动机加力燃烧室内气液两相湍流燃烧的控制方程,采用了双流体模型,高温辐射通量模型,双时间尺度湍流K-ε模型,应用SIMPLE方法解这组方程,对二维轴对称的涡喷加力燃烧室进行计算,证明应用本文的液雾燃烧模型可以得出合理的计算结果。      

内外函混合加力燃烧室扩压器内液雾湍流两相流场和浓度场数值模拟

 本文用κ-ε湍流模型和液雾轨道模型,对二维内外函混合加力燃烧室扩压器内有蒸发的液雾两相流动和混合进行了数值预报。其中对气相速度场和温度场的计算结果与实验符合令人满意,液雾两相流的计算结果在趋势上也是合理的,表明本文的模拟方法可用于加力燃烧室性能估算及优化设计。     

带隔热屏加力燃烧室热态流场计算

在任意曲线坐标系下对带有横向波纹隔热屏、外冷却通道和尾喷口的涡喷加力燃烧室的三维热态流场进行数值模拟.利用椭圆型微分方程和区域法生成三维贴体网格, k-ε模型预估紊流特性, EBU-Arre-henius紊流燃烧模型估算化学反应速率, 六通量模型预估辐射通量以考虑其对壁面和气流温度的影响.在非交错网格系统下应用混合差分格式离散控制方程, SIMPLE算法求解离散方程, 数值分析两种进口气流参数分布对加力室热态流场和壁温分布的影响.所得计算值与试验数据基本相符, 表明计算方法合理.      

用双流体—轨道模型模拟四角喷燃模型炉内三维湍流两相流动和煤粉燃烧

目前炉内两相流动和煤粉燃烧数值模拟中多半用颗粒随机轨道模型和单流体无滑移模型，这些模型都难以完整地给出三维空间内颗粒速度，浓度，湍流度分布的信息。主采用双流体－轨道模型（颗粒相连续介质－轨道模型）对一个四角喷燃模型炉内三维湍流两相流动及煤粉燃烧进行了模拟。此模型基于欧拉气相方程组、欧拉颗粒连续方程组和动量方程组以及拉氏颗粒能量方程和质量变化的方程，并使用k-ε-kp两相湍流模型，EBU－Arrhenius湍流燃烧模型，离散坐标辐射传热模型，煤粉颗粒的水分蒸发，热解挥发模型和焦炭燃烧的扩散－动力模型等。热态模拟中，为了减小为信散造成的影响，采用了扭转坐标法（将坐标扭转一定的角度使之与煤粉射流方程一致）。为了检验数值模拟，采用三维相位移普勒测速仪（PDPA）对于冷态模型炉内湍流两相流场进行了测量，得到了两相速度，湍流脉动及颗粒浓度的分布。分别对冷态模炉内两相流动和热态模型炉内三维两相流动和煤粉进行了模拟，冷态两相流动的计算与实验结果的对比表明预报的两相流场是合理的，热态模拟的结果给两相速度，气相温度、组分浓度及壁面热，显示出靠近出口处气相速度和温度分布不对称，造成一个局部高温区。     

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

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