航空发动机燃烧室流动数值计算中湍流模型的比较

分别采用标准k-ε湍流模型、RNGk-ε湍流模型、Realizablek-ε湍流模型以及雷诺应力模型,对某型航空发动机燃烧室流动进行了数值计算.近壁处理采用标准壁面函数法,计算得到速度矢量分布以及质量流量、湍流粘度比和湍流强度等参数.四种湍流模型计算的总体流动差别较小,但流场的细节有较明显的不同.标准k-ε模型、Realizable k-ε模型和雷诺应力模型的湍流粘性计算结果较为接近,而RNGk-ε模型计算的湍流粘性较小.     

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

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

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

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

两种湍流模型在跨声速绕流计算的应用研究

湍流模型在对复杂流场的数值计算中起着非常重要的作用.本文采用SA一方程湍流模型和SST k-ω二方程湍流模型,通过解耦求解雷诺平均N-S方程和湍流模型方程实现对亚跨超声速湍流流场的数值模拟.对NACA0012翼型和ONERA-M6机翼跨声速绕流流场进行了计算,对压力分布和激波位置与实验结果进行了细致的比较,并分析了不同离散格式、不同网格疏密及壁面函数对计算结果的影响,在计算过程中这两个模型体现出了较好的简捷性和健壮性.     

后掠压缩斜板进气道进口流场数值研究

对一种后掠压缩斜板进气道(Caret进气道)的进口激波结构、进口流动对管道后部流场的影响进行了数值分析,对进口倒圆修形的影响、激波覆盖面也进行了讨论。数值模拟采用N-S方程和标准k-ε模型加壁面函数方法。计算结果表明:超音速情况下,由于进口存在较强的正激波,进口斜激波强度的不均匀性,对管道内部流动看不出有明显的不利影响;但两压缩面交角处的流动对管道后部流场有较大的影响。      

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

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

过膨胀状态下轴对称收-扩喷管内外流场计算及分析

轴对称收-扩喷管在过膨胀状态下工作时, 喷管内部将产生强烈的激波以及附面层分离, 流态非常复杂.计算实践表明, 应用广泛采用的重整化群(RNG)k-ε湍流模型, 并结合标准壁面函数对其进行数值模拟时, 误差较大.针对标准壁面函数的不足, 采用增强型壁面函数处理固体壁面附近的流动, 对轴对称收-扩喷管的内外流场进行了数值模拟, 并与试验数据进行了分析比较.结果表明:增强型壁面函数更加适合处理此类存在强逆压梯度导致附面层分离的问题.      

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

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

雷诺应力模型在三维湍流流场计算中的应用

从雷诺应力模型出发,通过求解雷诺平均N-S方程组获得三维湍流流场的数值解。计算中比较了多种湍流模式,并进行了相应的流场计算。本文完成了两个典型算例。从算例1的计算结果与实验值比较中发现:采用雷诺应力模型(RSM)计算的三维流场比采用k-ε模型更贴近实验值;算例2采用了RSM模型及三维非结构网格,对一典型内流问题进行了三维流场计算。算例的数值实践表明:采用雷诺应力模型可以有效的计算各向异性的湍流流场;另外,发展非结构网格有助于模拟壁面附近的流动,并节省计算机内存。      

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。     

偏心对双级旋流器出口流场影响数值模拟

采用Fluent软件,数值模拟了一级旋流器和二级旋流器偏心情况下,双级旋流器出口流场变化,湍流模型采用标准k-ε模型,计算了3种不同偏心量(0,0.5,1.7mm)下的双级旋流器出口流场,并将计算结果与试验结果进行了对比.研究表明:偏心对双级旋流器出口流场影响较大,0.5mm偏心下双级旋流器出口的回流区发生畸变,径向速度峰值区域位于X/D=0.5~1.0范围内;轴向速度和径向速度的计算结果与试验结果吻合良好.      

湍流模型对发动机排气系统流场和红外辐射特征计算的影响

计算分析了采用k-ω SST(shear-stress transport)模型、RNG(renormalization group)k-ε模型配合标准壁面函数(RNG k-ε+标准壁面函数)、RNG k-ε模型配合增强型壁面函数法(RNG k-ε+增强壁面函数)3种湍流模拟方法,对典型涡扇发动机轴对称收敛喷管排气系统流场、温度场、组分质量浓度场及红外辐射特征计算结果的影响.研究表明:对于这3种方法,排气系统腔体内壁面压力分布基本一致;但在温度场、组分质量浓度分布以及红外辐射特征方面的差别较大.RNG k-ε+增强壁面函数计算得到的内涵壁面与燃气流的温度、CO₂与H₂O的质量浓度和红外辐射强度值最大,RNG k-ε+标准壁面函数的计算结果其次,而k-ω SST 模型的计算结果最小.通过与试验结果对比,发现k-ω SST模型对排气系统红外辐射特征的模拟更加合理一些.      

三维亚声速冲击射流流场的数值模拟

利用有限体积法对三维可压缩的N－S方程进行离散，对喷嘴的亚声速垂直冲击射流和斜冲击射流进行了数值模拟。网络的划分采用非结构性网络，湍流模型为RNGκ-ε湍流模型，在近壁处采用壁面函数进行修正，流场计算结果与实验值吻合一致。计算可得到垂直冲击流场结构的三个区域：自由射流、滞止区和壁面射流区域，并显示出垂直冲击点附近区域的卷吸回流、速度分布。冲击射流近壁处存在着回流现象，由于边界层流动和近壁回流引起距离固壁2D内的速度梯度变化较大。斜冲击在近壁处平行壁面速度的最大值，以及壁面上的总压最大值都偏离冲击点，并且相对应。     

环形燃烧室冷态流场数值模拟中的数学方法

针对带双级径向旋流器环形燃烧室火焰筒,研究不同紊流模型(standardk-εmodel,RNG(renormalization group)k-εmodel和代数应力模型(algebraic stress model))及算法(SI MPLE(semi-i mplicitmethod for pressure-linked equations),PISO(pressure i mplicit split-operator))对复杂流场求解准确性及效率的影响.火焰筒三维贴体网格采用型线定点与偏微分方程相结合的方式生成,用壁面函数法处理近壁区域流动.计算结果表明,PISO算法提高求解效率;不同紊流模型预测结果均与试验数据符合较好,其中代数应力模型最为准确,更适用于复杂流场的数值预测.      

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-ω预测得更准确.      

结构参数对圆转矩形喷管换热的影响

为了解结构参数对圆转矩形内喷管再生冷却换热的影响，设计了多个圆转矩形喷管，考虑了三种结构参数：转方位置、出口高宽比和出口圆角大小的影响。采用有限体积法求解三维可压缩的N-S方程对其内部流动和换热进行了数值模拟。湍流模型采用标准的k-ε双方程模型，壁面附近的流动和传热采用壁面函数法处理，速度与压力的耦合采用SIMPLE算法求解。结果表明：在型面一阶导数连续的情况下，转方位置对圆转矩形内喷管的换热影响不大；出口高宽比对圆转矩形内喷管的换热影响较大，出口高宽比不能太小，否则影响内喷管流场和换热；出口圆角大小影响内喷管周向上的温度分布，圆角太小造成周向温度分布不均匀。     

收敛-扩张喷管中运用激波矢量控制法的计算研究

采用二阶TVD格式的有限体积法耦合RNG k-ε湍流模型和非平衡壁面函数求解二维守恒型雷诺平均N-S方程，数值模拟了运用激波矢量控制法的二维收敛一扩张喷管中的流动现象，很好的模拟出了由于次流喷射所产生的激波及上下壁面的压力分布，并将计算结果与国外的实验结果进行了对比。     

计算风工程中 k-ε模型的边界条件研究

从 k-ε湍流模型控制方程入手,对两类满足平衡大气边界层理论要求的入口边界条件进行了详细分析,比较了二者的理论差异。然后以我国公路桥梁抗风规范中建议的 A 类风场为来流条件,定义了相应的 k-ε模型常数与壁面条件,使用两类边界条件进行数值模拟。计算结果表明：两类边界条件结合相应的模型常数与壁面条件,均能在数值模拟中构建基本满足规范要求的平衡大气边界层,但各有其适用范围。在此基础上,对模型常数 C 1和湍动能耗散率与数值模拟结果之间的关系进行了研究,采用修改模型常数与边界条件的方法,改善了数值模拟结果。     

飞行器腹部空腔绕流的数值模拟

采用分块对接"O"网格技术,基于雷诺时均NS方程与两方程k-ε模型计算DLR-F6模型的三维粘性流场.模拟了三种几何构型在马赫数M=0.75、0.85、1.5、1.9下的流动情况.计算结果表明开空舱使升力减小,阻力增加;不同的空腔长深比气动特性变化不明显.     

粗糙壁面表面传热数值计算比较

为了计算冰形模拟中的表面传热系数,给出了三种粗糙壁面表面传热系数计算方法,即传统的边界层积分法、基于k-ε湍流模型及标准壁面函数的计算方法(k-ε法)和基于Spalart-Allmaras湍流模型在粗糙壁面扩展的计算法(S-A扩展法).计算了直径为15cm、表面等效沙粒粗糙高度为1.35mm的圆柱,在Re分别为8.8×...