三维加力燃烧室湍流流场的数值模拟

采用SIMPLE方法,对涡扇航空发动机加力燃烧室内无化学反应的湍流流场进行了数值模拟。该加力燃烧室带有波瓣型混合器,并有多个环形稳定器及径向稳定器。湍流模型采用kε双方程模型。湍流流场数值计算结果合理。     

涡流冷壁旋流燃烧室的数值模拟

为了获得涡流冷壁旋流燃烧室的双旋涡结构对其燃烧流动性能的影响规律,改变切向入口的进口速度,对其内部流场进行了数值模拟。结果表明:在燃烧工况下,双旋涡结构较冷态时更难形成;随着切向入口的进口速度增加,涡流冷壁旋流燃烧室的燃烧效率较旋流燃烧室增大2%左右,壁面温度从1200K降低至600K,出口CO2的排放量降低;双旋涡结构会增加燃烧室的总压损失,但是增加的程度很小。     

渐扩通道内湍流流动与混合的数值模拟

应用标准形式的ｋ－ε湍流模型,对具有逆向压力梯度的渐扩通道内两股等温或不等温射流的湍流流动与混合进行了数值模拟。结果表明,时均速度场与温度场的计算结果与实验测量数据相符合。     

燃烧室中湍流涡旋流动的数值模拟

ＡＳＭ也进行了研究。研究的流场包括有约束的和自由的强涡旋流动。研究结果表明,二阶矩封闭模式在计算强涡旋湍流流动中很有效,代数模式则不能令人满意。计算结果可供燃烧室设计参考。     

旋流预混燃烧室流动混合的大涡模拟

采用大涡模拟的方法研究旋流预混燃烧室中的流动及燃料空气混合过程,选取亚格子湍动能模式为基准的涡黏模型,通过模拟涡核的运动过程定性的表征了拟序结构的发展;引入混合不均匀度的概率密度函数,研究入口边界条件尤其是入口速度对于涡团破碎及混合程度的影响规律;通过研究速度分量与燃料质量分数脉动的相关性,分析了不同方向的速度脉动对于混合过程的影响因子.这些均为进一步研究旋流预混燃烧室的燃烧不稳定性及污染物排放提供了理论上的指导.      

燃烧室三维湍流的数字模拟

文中以SIMPLE方法为基础,采用k-ε两方程湍流模型,发展了一种模拟燃烧室内三维湍流内流场的数值计算方法和计算程序。通过对四种冷流状态进行数字模拟,初步结果表明,理论计算与已有实验结果一致性良好。     

支板构型超燃冲压燃烧室流场数值模拟

数值模拟了带凹腔的支板构型超燃冲压燃烧室内部的流场结构。在计算过程中，控制方程采用了含组分守恒方程，并与k-ω双方程湍流模型紧耦合的质量平均Navier—Stokes方程，对方程中的对流项采用了空间为二阶精度的TVD格式，扩散项则采用了二阶中心差分离散。采用在流体力学时间步内的当地积分法克服了非平衡源项的刚性问题，保持了LU-SSOR隐式求解算法的高效率。通过数值模拟，对比研究了凹腔构型的位置、长深比(S／H)对燃料混合性能的影响。结果表明，长深比大的凹腔增加了燃料在腔内的停留时间，相应地改善了燃料的混合。位于燃料喷嘴之后的凹腔比位于下游的凹腔更有助于提高燃料的混合效率。     

高温升旋流燃烧室性能的数值分析

为了研究高温升燃烧室,采用数值研究的方法,对所设计的高温升多级旋流燃烧室和中心分级燃烧室与现有的单环腔燃烧室(SAC,原设计油气比为0.027)在设计油气比0.037条件下进行分析。研究结果表明:多旋流和中心分级的设计方案均可获得理想的燃烧性能参数,出口温度分布系数(OTDF)分别达到0.138和0.16,满足高温升燃烧室的设计指标,而SAC燃烧性能急剧恶化,不能满足设计指标。其中,中心分级燃烧室的技术优势十分明显,显示出作为高温升高热容燃烧室的发展前景;多旋流燃烧室则兼具高温升和低排放两方面优势。     

湍流燃烧模型在燃烧室数值模拟中的对比分析

以某单头部矩形燃烧室为研究对象,采用多种湍流模型和燃烧模型进行组合计算,模拟燃烧室内部的速度场和温度场,并对计算结果进行对比分析。结果表明:湍流模型主要影响火焰筒内部主燃孔横截面上游的速度分布,Standard k-ε和Realizable k-ε模型的速度场计算结果差异相对较小;湍流动能预测受湍流模型的影响较大,并具有一定规律性;不同组合模型对燃烧室内部和出口温度分布的局部细节模拟差异较大,燃烧模型影响最大;PDF模型计算的温度值较合理,另外三种燃烧模型在单步完全反应燃烧机理下的计算值偏高,计算获得的OTDF也相差较大。     

双旋流驻涡燃烧室湍流燃烧流动的数值分析

提出了带导流片的双旋流驻涡燃烧室以进一步发展低排放、高效燃烧技术。改变旋流器的结构参数以及旋向,对双旋流驻涡燃烧室内部燃烧湍流流场进行了数值模拟。结果表明:增大旋流器内外径比对燃烧效率有一定影响,NO_x排放降低了约8%。旋流器中心间距对燃烧性能的影响不如内外径比的大。不同旋向下燃烧室均可获得良好的双涡结构,但反旋时下游火焰筒的高温度梯度场会对筒壁造成强烈热冲击,同旋时NO_x排放降低了约22%,因此同向旋转结构总体上要优于反向旋转结构。     

基于CFD分析改进三旋流燃烧室头部设计

利用FLUENT,对一多级旋流燃烧室的冷态单相流场进行了数值模拟,基于该流场结果,对该燃烧室头部控制油雾混合的三级涡流器唇口处进行了改进设计。并对两种方案的燃烧室在相同试验条件下进行了贫油熄火油气比、大功率工作时的冒烟指数进行了实验研究。对方案1燃烧室单相冷态流场的结果表明,该燃烧室内部流场和预想的流动组织有较大差别,三级旋流器出口的气流,迅速向径向流去,主燃区的回流区结构不明显,且二级旋流出口和三级旋流出口之间有小回流区。该流场结构解释了方案1燃烧室有很好的贫油熄火性能,但是冒烟指数较高。方案2燃烧室的流场结构,改变了三级旋流出口气流方向,消除了与二级旋流器之间的小回流区。方案2燃烧室在与方案1燃烧室相同试验条件下的实验结果表明,贫油熄火油气比性能下降了7.3%,但发动机大功率时冒烟排放指数改善了29.4%。基于数值分析的燃烧室设计,使燃烧室扩大稳定工作范围的研究取得了进展。     

弯曲管道内湍流流动的数值模拟

为了对弯曲管道内的湍流流动特性进行有效的预测，用显示代数应力模式对－典型的弯曲管道内流动进行了数值模拟和分析，计算结果表明：（1）显式代数应力模式能对弯曲管道内流动的平均量（如速度、压力、摩阻等）进行较为准确的预测；（2）对模式中的ε方程做适当修正后，该模式能较好地模拟出流动曲率效尖对湍流特性的影响，是对雷诺应力模式的一个比较好的近似。此外，该模式在形式上与标准k－ε两方程模式类似，使用方便，具有较强的工程实用价值。     

双燃烧室冲压发动机超声速混合层混合增强数值研究

研究了具有双燃烧室冲压发动机超声速燃烧室工程背景的特定超声速混合层的混合增强措施——低速流入口加入沿流向振荡.研究结果表明:它对二维混合层的混合增强效果明显,加振荡后流场燃烧强度加大,但仍属于较弱燃烧强度;振荡频率决定着涡量聚集区域的位置和大小,振幅决定着涡量分配和聚集的均匀程度.      

燃气源复杂通道三维冷流数值模拟

采用两方程RNGκ-ε湍流模型，通过非结构网格上的SIMPLE算法，对某远程火箭弹控制系统燃气源复杂通道内流场进行了冷流数值模拟，获得了通道内各气流参量的三维分布，了解了燃气源内各处特别是控制气流流经的迷宫形多层导热套环结构内的气流流动特性，计算结果与实验结果基本吻合。     

涡流器燃烧室头部两相反应流数值模拟

在任意曲线坐标系中数值模拟带涡流器燃烧室头部的三维两相反应流场,其中,采用偏微分方程法生成块结构网格;气相在Euler坐标系下采用修正的κ-ε双方程模型处理紊流特性,用修正的EDM-Arrhenius紊流燃烧模型模拟化学反应速率,用六通量模型考虑热辐射的影响;液相处理采用在Lagrange坐标系下颗粒群轨道模型.在非交错网格系下采用SIMPLE算法及混分差分格式求解离散方程,采用单元内颗粒源法(PSIC)算法进行两相之间的耦合,计算所得的旋流数与试验数据相当接近.结果表明涡流器几何形状对燃烧室流场影响较大,涡流器的优化设计可以提高燃烧性能.     

离心式喷嘴内流场特性的数值模拟

基于单流体模型及气液两相混合分数的概念，建立了离心式喷嘴内流场的数学模型，并进行了数值计算，计算结果给出了离心式喷嘴内部液膜与气涡共存的流场结构及气液交界面的几何形状，计算得到的壁面压力分布、出口处液膜内轴向速度分布以及雾角大小都与实验结果吻合的很好，与常用的跟踪气液界面的ALE方法相比，本文建立的数值模拟方法具有简单，计算量小和易于实现的特点。     

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

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

Numerical simulation of flow in swirl injector for power-law fluid

In order to investigate the flow characteristics of swirl injectors for gelled propellants, which actually behaved as non-Newtonian power-law fluid, a swirl axisymmetric model was established to solve Navier-Stokes equations and VOF (volume of fluid) equation, and the power-law constitutive equation was used to describe the rheology characteristics of the gelled propellants. The film thickness and velocity distribution in the swirl injector under different flow conditions were studied numerically. The simulation results show that the increased geometry characteristic parameter of the swirl injector contributes to the decrease of liquid film thickness. The liquid film thickness is almost independent of the pressure drop. The rheologic parameters have great influences on the inner flow of swirl injector: by increasing the fluid consistency index K and power index n, both the axial and the swirl velocities decrease dramatically; higher fluid consistency index K and power index n make the liquid film thickness increase. When the viscosity is large enough, the air core in the injector would vanish.