应用修正的k-ε模型研究超声速H2/Air燃烧

为了研究超声速燃烧中流体可压缩性的影响,对标准k-ε湍流模型进行可压缩性修正(包括结构可压缩性修正和膨胀可压缩性修正两部分)。分别应用标准k-ε模型、修正的k-ε模型和雷诺应力模型(RSM),考虑氢气/空气详细化学反应机理(GR I-M ech 2.11机理,10组分,28基元反应),数值模拟有壁面限制的超声速混合层冷态及热态流场。结果表明:壁面和燃烧对湍流影响都很大;修正模型对冷态以及燃烧场的预测结果优于其它两个;修正模型预测的混合层厚度更薄,燃烧区域更窄,与实验结果吻合地更好。     

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

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

基于火焰面/进度变量模型的超声速燃烧IDDES模拟

为了清晰地描述超声速条件下的复杂流动与燃烧过程,建立了基于IDDES模型和火焰面/进度变量模型的湍流燃烧数值模拟方法,并应用于德国宇航中心(DLR)超燃冲压发动机的数值模拟。计算结果与实验数据符合良好,表明该方法可以较为准确地描述超声速湍流燃烧过程,也验证了经过可压缩性修正的火焰面/进度变量模型适用于超声速燃烧问题的求解。同时,利用k-ωSST模型对同一算例进行了计算,发现其对物理量分布的计算精度要低于IDDES模型,说明湍流输运过程的描述对火焰面/进度变量模型的计算精度有重要影响。     

一种超声速射流计算方法与实验验证

高温射流流场计算是尾焰辐射目标特性计算的前提,然而由于缺少可靠实验数据,针对湍流超声速射流的数值模拟多集中于低温射流,高温射流计算与实验的对比工作还很少见。利用k-ωSST双方程湍流模型,模拟了多个典型超声速射流实验的流场速度与温度分布,通过与实验结果进行对比,建立了一种超声速射流计算方法。首先,通过对比多个低温射流的实验与计算结果,探索了湍流模型中可压缩修正以及来流湍动粘性比对超声速射流计算结果的影响;进而,针对火箭发动机尾焰实验,计算尾焰流场与流场红外辐射,流场辐射计算结果与实验观测结果符合较好,进一步验证了计算方法。最终认为经过可压缩修正的k-ωSST双方程湍流模型结合湍动粘性比取值30可以作为超声射流计算中较为典型的湍流计算方法。     

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

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

液体煤油超声速燃烧数值分析的湍流模型

利用计算流体力学软件对煤油在所设计的双模态燃烧室内的喷雾燃烧进行了数值模拟。采用k-ε三种形式、k-ω两种形式和雷诺应力湍流模型,计算了在飞行马赫数为5,煤油与空气的当量比为0.551时的情况。通过将数值模拟结果与实验数据进行对比,证明这些湍流模型都可用于煤油超声速燃烧的数值计算。并且,k-εRNG模型和k-ωSST模型的计算结果与实验数据最为接近,可以认为这两个模型在模拟煤油超声速燃烧流场中是非常有效的;而雷诺应力模型效果较差。     

超声速燃烧的湍流流场数值模拟

选择了6个有代表性的湍流模型,对超声速燃烧湍流流场开展了数值模拟.采用的模型包括5个线性涡粘性模型和1个非线性模型,在非线性模型中同时引入了可压缩修正.数值结果表明,本文发展的数值方法可以成功地应用于求解可压缩、含化学反应的多组分N-S方程.同时引入可压缩性修正的二阶涡粘性模型给出的结果较好.     

用双流体模型模拟超声速气固两相流动

采用了双流体模型对JPL(Jet Propulsion Laboratory)喷管中气固两相流动以及超声速射流两相流动进行了数值模拟,并研究了可压缩两相流动中气相与颗粒的相互作用规律.自主开发的一般曲线坐标系下二维轴对称可压缩双流体程序Solve2D,对气相求解Navier-Stokes(N-S)方程组,采用k-ε湍流模型,颗粒相求解Euler方程组.对JPL喷管内的两相流场和湍流两相射流流场进行了数值模拟,研究了不同颗粒质量百分数以及不同颗粒直径时的气固两相流场的流动规律.      

等离子发生器燃烧流场数值模拟

采用漩涡破碎（EBU）燃烧模型、k-ε双方程湍流模型及SIMPLEC算法,对等离子发生器内部的燃烧流场进行了数值模拟,得到了关于等离子燃烧流场相关特性参数分布。该数值模拟对等离子发生器的研究具有指导作用。     

两种差分格式和两种湍流模型在轴对称冲击射流数值计算中的比较

用延迟修正QUICK和乘方两种差分格式，结合标准k-ε模型和一种基于重正化群（RNG）思想的k-ε模型，对半封闭轴对称涡流冲击射流场进行了数值模拟。通过与实验结果的比较，考察了两种差分格式对计算结果的影响，并讨论了两种湍流模型对冲击射流场的数值预测能力。研究表明：当都采用标准k-ε模型时，QUICK格式得到的结果比乘方格式更接近于实验值，尤其对湍流能分布；当都采用QUICK格式时，RNG k-ε模型的预测值比标准k-ε模型更准确。     

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

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

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

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

横向脉冲喷流超声速燃烧流场数值模拟与分析

为了研究横向脉冲喷流对燃烧室流场结构及燃烧特性的影响,对带后台阶的超燃冲压发动机燃烧室横向氢气喷流超声速燃烧流场进行了数值模拟。采用有限体积法求解多组元Navier-Stokes (NS)方程,分别对定常喷流与两种脉冲喷流超声速燃烧流场进行了数值计算,对比分析了三种氢气喷注方案下的流场结构、氢气掺混燃烧特性及燃烧室总压恢复特性。研究表明:脉冲干扰仅对燃烧室内局部流场产生周期性影响,且有利于氢气的横向扩散;采用与定常喷注喷流参数相同的脉冲喷注方案时,可在脉冲干扰区内保持氢气掺混量与定常喷注基本相同,氢气燃烧量与总压恢复系数整体上与定常喷注基本一致;采用特征时间内喷出氢气质量流量与定常喷注相同的脉冲喷注方案时,在燃烧室出口处氢气掺混量和燃烧量较定常喷注分别提高了21.94%和32.24%,总压恢复系数仅减小4.17%。采用脉冲喷注方案增加了燃料与空气接触面面积,对燃料掺混及燃烧起促进作用。     

隐式耦合求解燃烧室内喷雾两相流场

在任意非结构网格下,采用开源程序对某型涡喷发动机环形燃烧室内的喷雾燃烧过程进行数值模拟.首先选用k-ε双方程模型、EBU(eddy break up)-Arrhenius紊流燃烧模型、颗粒随机轨道模型等模拟流场中的物理过程.然后用隐式耦合算法求解气相流动、拉格朗日方法跟踪液滴相在流动区域的运动,并考虑气相与液滴相的完全双向耦合作用.仿真实验得到的 燃烧室内温度分布、马赫数分布、速度分布、压力分布均较合理,温度分布与文献中的基本一致,这表明隐式耦合方法能够较好地模拟燃烧室内流场,反映其内部流动特点.      

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

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

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

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

用Tam-Thies湍流模型研究带小突片的喷管流场特性

 首先,采用标准k-ε模型、可实现（realizable）k-ε模型、重正化群（RNG）k-ε模型和Tam-Thies模型4种湍流模型模拟了不同工况下的两种轴对称喷管和带有4片小突片的喷管流场,并将计算结果与实验值进行对比。结果发现：Tam-Thies模型所得出的模拟结果与实验值符合最好。然后,用Tam-Thies模型模拟了小突片后倾角不同的喷管的流场,和没有小突片的喷管的计算结果对比,发现：带有小突片的喷管的尾喷流的核心区长度变短,喷管出口下游气体的混合加剧;而当小突片后倾角增大时,核心区长度先减小后增大;在每个小突片下游产生一对方向相反、强度相同的流向涡,当小突片后倾角增大时,流向涡强度呈现出先增大后减小的趋势,而喷管的引射系数增益和推力损失系数都逐渐减小。     

平行平板内连续布置肋块流场与传热的数值模拟

为了研究平行平板内连续布置8个肋块的流场结构和传热现象, 利用计算流体力学软件FLU-ENT6.2内置的标准k-ε湍流模型、RNGk-ε湍流模型和Realizablek-ε湍流模型对此进行了数值模拟并对结果作了比较.流场结果显示在相同来流条件下, 由Realizablek-ε模型预测的肋块间逆时针方向回流区最大, 标准k-ε模型最小.数值结果表明从第2个肋块开始加热段表面压力系数Cp和努塞尔数Nu开始呈周期性变化.同时将第7个肋表面的努塞尔数数值解与实验值进行了比较, 结果为三种湍流模型都能准确模拟出努塞尔数的变化趋势.计算结果可为工程设计和使用提供参考.