高位进气、径向出流的旋转腔内流动与换热的数值研究

用标准k-ε模型模拟了航空发动机涡轮盘冷却系统的高位进气、径向出流盘腔模型的湍流流场和温度场,分析了其盘面平均Nu数随旋转雷诺数Reω和进气雷诺数Re的变化。     

旋转盘腔紊流流动的数值研究

为了比较近壁紊流模型对旋转盘腔内流动特征模拟精度的影响，采用在核心区选用标准k-ε模型、在近壁区选用单方程和Launder-Sharma两方程两种低雷诺数模型的区域紊流模化方法，对具有径向出流的转-静盘腔和具有径向出流的共旋盘腔结构的流动进行了数值计算研究，计算结果表明：在计算的两个算例中，前者采用近壁单方程模型要稍优于两方程模型，后者采用近壁两方程模型要明显优于单方程模型，因此近壁区紊流特性的处理对旋转盘腔内的流场模拟精度具有重要的影响。     

应用不同湍流模型计算叶栅流场的结果与比较

采用有限体积TVD-LW格式求解N-S方程,使用了B-L模型、高雷诺数k-ε模型和低雷诺数k-ε模型等多种湍流模型。并且用具体算例说明本文所用程序的可靠性,比较了3种的湍流模型的特点。B-L模型最大的优点是简单、使用方便。k-ε方程适应性强,分为低雷诺数k-ε方程和高雷诺数k-ε方程两种。高雷诺数k-ε方程比低雷诺数k-ε方程精度略低,但所需计算时间和存储量都比后者小,网格也可以粗些。k-ε模型对粘性层、分离区有比较高的分辨能力,高雷诺数k-ε模型比较适于一般数值计算。      

反旋进气盘腔内流动与换热的数值模拟

应用RNGk-ε湍流模型对围屏上带反旋喷嘴的径向内流旋转涡轮盘腔内的流动与换热进行了数值模拟,揭示了盘腔内的压力损失及冷气流量、旋转雷诺数、湍流参数等因素对盘腔内流动换热的影响.计算结果表明:盘腔内的流动结构主要由湍流参数决定;在某一旋转雷诺数下盘腔内压力损失随冷气流量的增大而呈现S型变化;反旋喷嘴的应用能有效地降低盘腔内的压力损失;转盘附近的努赛尔数随冷气流量及旋转雷诺数的增大而增大.      

高位垂直进气转静系旋转盘腔流场的实验

利用粒子图像测速技术(简称PIV)对高位垂直进气的转静系旋转盘腔流场进行了实验研究.该盘腔由一个旋转盘、一个静止盘及静止的外围盘罩组成.实验结果表明:PIV粒子成像技术可以应用于旋转盘腔结构的流场测试, 在本实验的工况范围内, 大部分区域的流动已为湍流, 两盘间的流动结构明显具有Batch-elor流型的特点, 转盘和静盘两个表面附近形成各自独立的边界层, 两边界层之间有一个旋转核心, 核心区内旋流系数β比大间隙封闭系统湍流流动时旋流系数的值0.43要小.流量系数和旋转雷诺数的变化对旋流系数β的影响较小, 但转静间隙的变化对β的影响较为显著.      

导流叶片宽度对预旋系统性能影响的数值研究法

为研究导流叶片结构尺寸对盘腔预旋性能的影响,采用RNG k-ε模型对导流叶片无量纲宽度为0、0.1、0.2、0.4、0.6、0.8和1.0的预旋系统在不同旋转雷诺数工况下进行了数值研究.结果表明:导流叶片能够提高喷嘴压比和温降系数,旋转雷诺数越大,提高的幅度越大;旋转雷诺数较大时,喷嘴压比和温降系数随导流叶片无量纲宽度增大而增大,无量纲宽度超过0.6后,趋于平稳.当导流叶片无量纲宽度在0~0.2范围内,总压损失系数随导流叶片无量纲宽度的增大而增大;当导流叶片无量纲宽度大于0.2时,总压损失系数随导流叶片无量纲宽度的增大基本不发生变化.     

导流叶片宽度对预旋系统性能影响的数值研究

为研究导流叶片结构尺寸对盘腔预旋性能的影响,采用RNG k-ε模型对导流叶片无量纲宽度为0、0.1、0.2、0.4、0.6、0.8和1.0的预旋系统在不同旋转雷诺数工况下进行了数值研究。结果表明:导流叶片能够提高喷嘴压比和温降系数,旋转雷诺数越大,提高的幅度越大;旋转雷诺数较大时,喷嘴压比和温降系数随导流叶片无量纲宽度增大而增大,无量纲宽度超过0.6后,趋于平稳。当导流叶片无量纲宽度在0~0.2范围内,总压损失系数随导流叶片无量纲宽度的增大而增大;当导流叶片无量纲宽度大于0.2时,总压损失系数随导流叶片无量纲宽度的增大基本不发生变化。     

反向旋转盘腔内部流动特性

对反向旋转盘腔系统内部流动特性进行了数值模拟,计算结果与已有文献中的实验结果对比表明,采用的数学模型和计算方法是准确的.数值模拟结果表明:反向旋转盘腔内流动结构较为复杂,转速比对流动结构有重要影响.根据流动结构特点研究了旋转雷诺数和转速比对盘腔内压力分布以及转盘壁面摩擦力矩的影响,旋转雷诺数越大,盘腔内的压力越低,转盘壁面的摩擦力矩越大.反向同速旋转有助于增大盘腔内的压力,但与其他转速比相比其摩擦力矩最大.      

径向进气轴向出流旋转盘腔总压损失特性研究

在径向进气轴向出流的旋转盘腔中,在哥氏力的影响下,流体速度的切向分量和流体的总压损失均增大,用数值计算方法研究其损失特性缺乏验证。采用试验方法研究了径向进气轴向出流旋转盘腔的压力损失特性,测试了不同工况下的旋转盘腔的进、出口总压,分析了流量系数和旋转雷诺数对径向进气旋转盘腔总压损失的影响规律。试验结果表明：旋转盘腔的总压损失随旋转雷诺数的增大而增大。随流量系数的变化规律较复杂,在较小旋转雷诺数下,总压损失随流量系数的增大而增大;在较大旋转雷诺数下,总压损失随流量系数的增大先减小后增大。     

轴向通流旋转盘腔流动换热的数值研究

为了探究旋转盘腔内的流动和换热规律,对轴向通流旋转盘腔进行了非稳态数值模拟,将计算结果与实验数据进行了对比,探究了流动不稳定性的发展过程,分析了盘腔内流动结构和盘面换热特性随旋转雷诺数的变化规律。结果表明:旋转引起的正旋涡从盘罩附近开始发展,随转速的增大而变大,挤压低半径区域的强迫对流区,最终扩展到整个盘腔,盘腔中轴面的涡对数与流动不稳定性的强度有关。上游盘和下游盘的高半径区域换热强度随转速的增大而增强,下游盘低半径区域的换热强度在低转速下由于冲击作用而较强,但该冲击作用随转速的增大而减弱,低半径区域的换热强度也就随之减弱。当旋转雷诺数增大到4.94×105时,下游盘低半径区域受到的冲击作用减小到可以忽略。     

一种近壁理论湍流模型初探

在传统的k-ε湍流模型基础上,在近壁区放弃求解ε方程而代之以根据共振三波理论模型在近壁区系统计算的ε值构成封闭,通过对NACA0012二维翼形流场的计算,初步验证并适当修正了新的理论模型,由于注意到了理论模型计算结果的无量纲性,从而改善了该流场的计算结果,得到了与实验更为相符的C_p分布。本文的分析还表明,逆压力梯度对ε的影响非常显著,验证了理论分析和计算得到的结论。      

不同去旋角度进气共转盘腔内流动与换热研究

以带有去旋喷嘴的共转盘腔为研究对象,运用RNG k-ε湍流模型,分别对3种喷嘴角度（90°、60°、30°）的3维模型进行了数值计算,揭示了去旋角度对盘腔内气体切向速度、盘腔内压降和转盘表面换热的影响规律。研究结果表明：随着去旋角度的减小,盘腔内气体的切向速度减小,盘腔内压降减小,转盘表面的局部努赛尔数Nur增大。     

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

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

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

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

为了研究平行平板内连续布置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-ε模型

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