风扇进气畸变三维非定常数值模拟技术研究

采用三维欧拉方程加体积力模型的方法, 发展了一套三维非定常进气畸变数值模拟程序, 以实现风扇/压气机在进气畸变条件下的全环面三维非定常数值模拟.利用该程序模拟计算了某航空发动机风扇进口稳态总压畸变条件下特性和流场结构, 并与试验结果进行了对比, 结果表明该程序用较小的计算机资源即可以实现对风扇全环面三维非定场畸变流场的模拟, 能够以较高的精度预估气流畸变对压气机/风扇特性的影响.      

不同湍流模型对旋涡流动的数值模拟

针对某型旋流器,基于Navier-Stokes方程,湍流模型分别选择renormalization group(RNG)、κ-ε模型和Reynolds stress model(RSM)模型,采用计算流体动力学方法对旋流器内部旋流流场进行数值模拟,通过对计算结果的分析以及与试验结果的比较,进行了旋涡流动数值模拟中湍流模...     

一种强耦合Spalart-Allmaras湍流模型的RANS方程的高效数值计算方法

在工程实际中,一方程湍流模型或两方程湍流模型的求解通常和雷诺平均Navier-Stockes (RANS)方程的求解是解耦的,也称之为松耦合求解.在松耦合求解过程中,RANS方程和湍流模型方程通常采用不同的数值方法异步求解.这种求解方式很容易产生因两者计算精度不一致而引起的额外数值耗散.为了消除这种耗散,将RANS方程与Spalart-Allmaras模型方程耦合成一个系统方程——强耦合RANS方程,并发展了一种用于求解该系统方程的高效强耦合算法,其中对流项离散采用了Roe格式,时间项的离散采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式,为了提高计算效率,采用了三层V循环多重网格方法.通过翼型/机翼和振荡翼型/机翼等算例验证了本文发展的强耦合算法不仅具有较好的收敛性,而且计算精度明显优于松耦合算法,特别对于阻力的预测,强耦合算法更加准确.     

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

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

某涡扇发动机考虑级间引气的二维数值模拟

针对采用S₂流面流场模拟软件(AES-S2),对某涡扇发动机风扇、压气机(考虑级间引气)和整机进行地面0km及高空21km低雷诺数工况下转速特性的评估和流场分析.AES-S2数值仿真模型由带有黏性项的二维欧拉方程、燃烧模型和损失模型构成;采用任意曲线坐标系,可保证复杂几何边界的计算精度;采用隐式高阶精度Godunov求解格式,数值稳定性高,可自动捕捉激波.仿真结果表明:AES-S2计算稳定性好,收敛速度快;风扇和压气机的计算精度高于整机;整机二维计算满足工程精度需求,可为部件研制和整机部件匹配提供性能参数依据.      

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

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

跨音速低展弦比轴流风扇转子内部流场的数值研究

采用三维雷诺平均Ｎ－Ｓ方程程序求解跨音速压气机风扇转子内部流场。这个程序采用有限体积显式时间推进方法求解控制方程,方程自变量选取在控制体顶点,采用简单的Ｈ型网格离散方程,离散后的代数方程用两步Ｒｕｎｇｅ－Ｋｕｔａ法进行时间推进求解,使用多重网格和当地时间步长方法加速计算收敛。程序中使用了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ模型模拟湍流影响。用该方法计算了ＮＡＳＡ６７＃跨音速风扇转子的性能曲线,并重点计算最高效率点工况下内部流场。计算结果与实验结果的比较,证明该法可以用来模拟跨音速风扇转子内部流动。     

倾斜回流燃烧室燃烧性能的数值研究

采用Fluent软件对某倾斜回流燃烧室进行了三维数值模拟,研究了燃烧室内部流场结构、温度分布及压力损失等特性.数值计算采用PDF燃烧模型、Realizable k-ε湍流模型、SIMPLEC压力速度耦合算法以及二阶精度迎风差值格式.计算结果与试验结果吻合较好,对回流燃烧室的设计具有一定的指导意义.     

跨声速下迎风格式的比较研究

研究了不同迎风格式的耗散机理,通过求解N-S方程,比较了其耗散性能.采用SST,EASM k-ω湍流模型,分析了湍流数值模拟中不同格式的计算精度,重点比较了迎风格式抑制壁面振荡的能力.选取跨声速收敛-扩张喷管内流场算例,从不同角度考察了格式的计算性能,数值计算结果和试验结果吻合,获得了不同格式耗散性能的比较结果.     

基于混合RANS/LES方法与FW-H方程的气动声学计算研究

在气动噪声的数值计算过程中,非定常流动的求解精度对声学计算结果有着重要的影响.以机体噪声计算的标准算例双圆柱绕流气动噪声问题为研究对象,采用基于非线性k-ε湍流模型的限制数值尺度(LNS)方法对双圆柱绕流进行了数值模拟,将计算得到的气动特性和流动特征与相应的试验结果进行了对比分析.为了求解远场观测点处的气动噪声,在精确求解双圆柱绕流流动的基础上结合基于FW-H(Ffowcs Williams-Hawkings)方程的声类比方法进行数值计算,并通过圆柱体的展向相关性对计算结果进行了修正,将得到的最终结果与相应的声学试验结果进行了对比,两者吻合良好,表明该数值方法是准确、可靠的.     

应用各向异性湍流模型对气膜冷却的数值研究

应用引入各向异性系数的k-w湍流模型和标准k-ε模型,对平板气膜冷却经典实验进行了数值研究,分析了各向异性系数的引入对气膜孔下游流场与温度场的影响机理,探讨了该系数的适用性问题.研究表明:①展向各向异性系数的引入,增加了流体在展向的湍流脉动,强化了气膜冷气沿展向的扩散;②各向异性模型的计算结果在低吹风比(M=0.5)的工况下与实验值符合得较好;③各向异性系数在气膜脱离壁面与主流掺混的区域适用性较差.      

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

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

用Tam-Ganesan湍流模型研究排气系统红外辐射特性

通过三维流场、介质的质量浓度场与辐射场耦合求解方法,计算了1台有中心锥的轴对称喷管的空腔-喷流组合红外辐射特性.利用专门为自由剪切流动发展的Tam-Ganesan湍流模型,有效提高排气系统红外辐射模拟精度.采用有限体积法(FVM)结合窄带模型(NBM)计算了排气系统在3~5μm波段的辐射传输问题,并考虑大气衰减作用.对比了采用标准 k - ε 模型和Tam-Ganesan模型对排气系统核心区长度的影响,以及由此引起的排气系统红外辐射特性的变化.通过与试验数据对比证明了该结果的准确性,结果表明选择合适的湍流模型对计算排气系统红外辐射特性有重要的意义.      

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

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

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

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

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

用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个肋表面的努塞尔数数值解与实验值进行了比较, 结果为三种湍流模型都能准确模拟出努塞尔数的变化趋势.计算结果可为工程设计和使用提供参考.      

UDMH/NTO火箭发动机尾焰流场特性数值仿真

采用AUSM(advection upstream splitting method)空间离散格式,探究了常用的3种k-ε湍流模型和2阶、3阶迎风格式对激波捕捉效果的影响。采用标准k-ε双方程模型和2阶迎风格式对不同高度和不同来流马赫数下的偏二甲肼/四氧化二氮(UDMH/NTO)火箭发动机尾焰流场进行喷管-尾焰流场一体化仿真。复燃反应采用12组分18步化学反应模型,喷管入口参数由热力计算给出。结果表明:随着来流马赫数的增加,波节数逐渐降低;随着飞行高度的上升,尾焰影响区域逐渐扩大;复燃反应造成混合区中的O2、N2质量分数大幅下降,而混合区中O、OH、NO质量分数则有所上升。      

封闭转静盘腔中紊流的数值模拟

对一个旋转盘腔的流场进行了数值计算研究。该盘腔由一个旋转盘,一个静止盘及静止的外围屏组成。在盘腔的核心区采用标准k-ε模型,在近壁区分别采用了两种低雷诺数紊流模型,即(1)单方程模型和(2)Launder-Sharma低雷诺数k-ε模型,并且在第二个模型中考虑了固体旋转引起的耗散率降低的修正。