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

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

环形燃烧室两相燃烧数值研究

根据环形燃烧室设计需要,运用计算机模拟技术研究模型环形燃烧室三维有旋液雾燃烧流场,研究两相流动与混合,液雾蒸发对燃烧过程的影响。两相流动模型采用颗粒群轨道模型,气相紊流输运方程求解采用双方程k-ε模型和EBU-Arrhenius紊流燃烧模型,热辐射采用热通量法模型,本文提供计算方法可以用来研究不同性能参数对环形燃烧室两相燃烧特性影响。     

喷管面积调节精度对某型涡扇发动机加力性能影响的数值仿真

应用考虑混合室、加力燃烧室、主燃烧室和外涵道容积效应和变几何通道、主燃烧室供油量和加力燃烧室供油量等控制因素的发动机动态过程的仿真模型,对某型变几何涡扇发动机接通加力动态过程进行了仿真计算,并研究了喷管临界面积调节精度对发动机过渡过程的影响,得到了与实验数据符合较好的结果。仿真模型与方法可为加力系统调节提供理论基础。      

加力燃烧室三维两相化学反应流数值模拟

利用FLUENT软件对某型发动机加力燃烧室进行了三维两相化学反应流的数值计算;应用有限速率化学反应模型和离散相雾化模型描述了加力燃烧室中复杂的三维两相化学反应过程。计算中给出了冷热态的流场变量和化学组分分布,以及燃烧效率和总压恢复系数沿流程的变化关系,与试验结果符合较好。     

涡轮风扇发动机接通加力过程的数值模拟

发展了一种低涵道比混合排气加力涡轮风扇发动机接通加力时的数值模拟方法。该方法考虑混合室加力燃烧室、主燃烧室和外涵道的容积效应,风扇、压气机、高压和低压涡轮等部件的特性,加力燃烧室供油管的充油时间。以混合室加力燃烧室、主燃烧室和外涵道3个容积室的动态方程为基础,用很小的时间步长逐渐逼近,计算不需要迭代,速度快,精度满足工程要求。      

加力燃烧室用齿冠状混合器的数值研究与验证

基于CFD(computational fluid dynamics)技术,分别采用有限速率/涡耗散(finite-rate/eddy-dissipation)、涡耗散(eddy-dissipation)和涡耗散概念(EDC)燃烧模型,对一种曲面齿冠状混合器作用下的加力燃烧室进行数值计算.计算结果表明:曲面齿冠状混合器下游产生明显的流向涡对,流向涡的存在可以增强内外涵冷热气流掺混;与直面齿冠状混合器相比,曲面齿冠状混合器下游流向涡强度大于直面齿冠状混合器,其混合效率高于直面齿冠状混合器;加力燃烧室内氧气和二氧化碳质量分数分布结构表明了涡耗散模型的计算更为合理;计算结果可为加力燃烧室混合器选型和燃烧数值模拟提供依据.      

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

本文对涡扇发动机三维加力燃烧室内的两相湍流燃烧过程进行了数值模拟。两相流模型采用双流体模型,湍流模型采用标准 k-ε模型,湍流燃烧采用涡旋破碎 ( EBU)模型,数值方法采用 LEAGAP算法。计算结果定性合理。      

混合排气加力涡扇发动机过渡态的数值计算

考虑混合室加力燃烧室、主燃烧室和外涵道的容积效应，风扇、压气机、高压和低压涡轮等部件的特性，加力燃烧室供油管的充油时间等因素，提出了模拟低涵道比混合排气加力涡轮风扇发动机过渡过程的计算方法，以混合室加力燃烧室、主燃烧室和外涵道3个容积室的动态方程为基础，详细分析计算了苛发动机过渡态的性能，用很小的时间步长逐渐逼近，计算不需要迭代，速度快，精度满足工程要求。     

旋流杯套筒混合段长度对点火特性的影响

为研究旋流杯套筒过长的混合段长度对燃烧室点火特性的影响,设计了两种不同的套筒混合段长度试验件,在常温常压条件下,采用单头部模型燃烧室进行了起动点火性能试验,同时采用马尔文进行了燃油雾化颗粒度测量,并结合Fluent软件的数值模拟对结果进行分析.研究结果表明:套筒混合段长度对燃烧室头部流场和油雾分布有很大影响,过长的套筒混合段长度将造成液雾与套筒壁碰撞,限制旋转射流和液滴的径向运动,造成喷雾锥角减小了32%,贫油起动点火油气比增加了20%左右.      

加力涡扇发动机喘振与消喘模拟

基于气动耦合原理,发展了涡扇发动机风扇出口分流环计算模型,并进一步基于多级轴流压气机系统的逐级单元控制体模型,建立了混合加力涡扇发动机喘振及其喘振消除的数值模拟方法,并将其应用于某型全台加力涡扇发动机过失速特性模拟。在加力状态下采用收小尾喷管喉道面积诱导了风扇压气机过失速—喘振过程,并给出了预先对加力燃烧室实施切油消除风扇压气机喘振数值模拟结果。      

两相紊流燃烧过程的数值计算

 本文描写了作者发展的两相紊流燃烧理论及数值方法;简介了作者自编的计算机程序GEMCHIP;提供了将本方法应用于几种燃烧室流动的计算结果。     

一个新的可压缩性修正的k-ε模型

考虑结构可压缩性修正的影响,发展了一个同时考虑结构可压缩性修正和膨胀可压缩性修正的k-ε湍流模型,新模型包括Chang可实现性、Heinz湍流动能产生项以及Sarkar可压缩性三部分修正.新模型扩宽了以往发展的可压缩性修正模型的适用范围,适用于高超声速(M>5)复杂湍流流动中.通过对多个复杂超声速横侧射流工况的计算,验证了新模型的预测效果.与实验结果相比表明,几个工况下新模型的预测精度都显著高于标准k-ε模型.流体分离强度越大,新模型的修正效果越显著.与标准k-ε模型相比,新模型计算结果与实验更加接近.     

二元矢量喷管流场的N－S方程计算

对二元矢量喷管在不同工况下的流场进行了数值模拟。求解的控制方程为三维非定常Ｎ－Ｓ方程。空间方向用有限体积离散，时间方向用Ｒｕｎｇｅ－Ｋｕｔｔａ法推进。粘性项的计算上采用了一种有限体积框架下新的离散化方法。湍流模型采用了广泛使用的Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数模型。计算表明，本文方法能准确地计算二元矢量喷管、非矢量喷管在设计与大设计落压比下的流场。     

可压流双方程紊流模型及其应用

 本文提出了可压流紊流模型。所用的方法类似于K-KL双方程模型,但考虑了气流压缩性的影响。方程中的常数已由实验确定。用本模型计算了管内超音速混合流动,其理论计算曲线与实验数据相当符合,证明该模型较不可压紊流模型更为适用。     

突扩回流与大速差射流回流湍流气－固两相流动的数值模拟

本文建立了气－固两相流的多连续介质模型，并分别对突扩回流与大速差射流回流这两种与实际工程应用密切相关的复杂流动条件下的受限湍流气一固两相流动进行了数值模拟。单相湍流回流流场的模拟结果与实验符合良好。两相流动的模拟结果表明，与突扩回流相比，大速差射流回流更有利于实现较高的颗粒浓度与气相回流区及低速尾流区的匹配，增加两相滑移速度，延长颗粒停留时间及强化两相间的混合。     

Algebraic turbulent energy flux models for hypersonic shear flows

Algebraic truncation models were derived for the turbulent energy flux for high-speed shear flows with and without thermal and mechanical non-equilibrium effects. To facilitate closure, fluctuating dilatation moments were modeled via conservation of mass truncations and the inhomogeneous terms were neglected. The resulting energy flux model provided significant improvements, up to 20%, in the temperature predictions over the gradient diffusion model for equilibrium flows over a Mach number range of 0.02–11.8. Our experiments showed that vibrational relaxation significantly altered the basic decay rate of grid turbulence. Also, our models indicated that compressibility introduces an explicit pressure gradient dependency. Analyses of a governing parameter indicated that these terms are negligibly small for low speeds, but may be important, or even dominant, for high-speed flows.     

微分雷诺应力模型在激波分离流中的应用

由于Boussinesq涡粘性假设的湍流模型在一些流动中如激波分离流中预测失真,本文介绍一种基于完全雷诺应力方程的微分雷诺应力模型(DRSM),应用压缩性修正和改进耗散率方程构造出应用于可压流体的可实现DRSM,SSG压应力模型测试了两组参数。将其应用于ONERA M6机翼在跨声速大攻角的计算,结果表明DRSM能更好地模拟激波位置和流动分离。     

NUMERICAL STUDY OF TURBULENT SWIRLING FLOW IN COMBUSTOR

ＮＵＭＥＲＩＣＡＬＳＴＵＤＹＯＦＴＵＲＢＵＬＥＮＴＳＷＩＲＬＩＮＧＦＬＯＷＩＮＣＯＭＢＵＳＴＯＲＦｕＳｏｎｇ（ＤｅｐａｒｔｍｅｎｔｏｆＥｎｇｉｎｅｅｒｉｎｇＭｅｃｈａｎｉｃｓ，ＴｓｉｎｇｈｕａＵｎｉｖｅｒｓｉｔｙ，Ｂｅｉｊｉｎｇ，Ｃｈｉｎａ，１０００...     

High-order accurate,low numerical diffusion methods for aerodynamics

In recent years numerical methods have been widely used to effectively resolve complex flow features of aerodynamics flows with meshes that are reasonable for today's computers. High-order numerical methods were used mainly in direct numerical simulations and aeroacoustics. For many aeronautical applications, accurate computation of vortex-dominated flows is important because the vorticity in the flow field and the wake of swept wings at an incidence and rotor blades largely determines the distribution of loading. The main deficiency of widely available, second-order accurate methods for the accurate computation of these flows is the numerical diffusion of vorticity to unacceptable levels. Application of high-order accurate, low-diffusion numerical methods can significantly alleviate this deficiency of traditional second order methods. Furthermore, higher-order space discretizations have the potential to improve detached eddy simulation predictions of separated flows with significant unsteadiness. Recently developed high-order accurate finite-difference, finite-volume, and finite-element methods are reviewed. These methods can be used as an attractive alternative of traditional low-order central and upwind computational fluid dynamics methods for improved predictions of vortical and other complex, separated, unsteady flows. The main features of these methods are summarized, from a practical user's point of view, their applicability and relative strength is indicated, and examples from recent applications are presented to illustrate their performance on selected problems.     

Numerical simulation of airfoil flows with a turbulence model

Turbulent flows over AS240 and NACA4412 airfoil were simulated numerically using a two-equation turbulence model named k-ξ model. The predictions of velocity profiles and the pressure coefficient of airfoil AS240 at 8°/19° attack angle and NACA4412 at 13.87° attack angle were calculated. The results were compared with those using k-ε and k-ω models,as well as experimental data. It indicates that the new k-ξ model offers more realistic prediction than the other two models. The main finding shows that the new k-ξ model is good at predicting separated flows around airfoils, and it captures the flow feature of pressure-induced separation adequately. All calculations are implemented as per openFOAM 1.7.1(open source field operation and manipulation).