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

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

锥柱形药柱环形槽内传热现象的数值研究

采用标准k-ε湍流模型和离散传递热辐射模型对不同来流雷诺数下锥柱形药柱前向环形槽内的流场与传热进行了数值研究。对槽内各表面与热气流之间的换热量进行了比较,并对努塞尔数沿环形槽各表面的变化进行了分析。结果表明来流雷诺数的大小影响环形槽内的旋涡数和旋涡区域的分布,进而影响槽内温度的分布,最终影响槽内装药表面与点火热气流之间的换热量。整个计算结果可为装药设计和分析推进剂的热解燃烧提供参考。     

多项式温度分布下的自由盘湍流换热数值研究

在航空燃气涡轮发动机的二次流空气系统中,涡轮盘腔的流动和换热问题具有复杂的热边界条件,其流场和换热特性对发动机的设计具有重要作用。基于单项式和多项式边界条件下的自由盘的湍流流动换热模型,应用FLUENT进行数值研究。通过与理论值的对比,验证了使用Relizible k-ε湍流模型用于研究自由盘面上湍流的适用性。研究发现,对于不可压的湍流,盘面的局部努塞尔数受温度分布的控制,盘面温度为任意n次多项式分布下的局部努塞尔数可由盘面温度按0-n次单项式分布时的结果导出。     

用一种非线性k-ε模型计算突扩台阶后的流动

对一种新的非线性k-ε模型的无量纲系数进行了重新拟合.用新的系数计算了槽道流动、方截面管流和等应变率,平均流场无旋流动的正应力分布.计算结果与实验资料符合很好.将该模型并入常规的k-ε模型的计算程序中,并对突扩台阶后的湍流流动进行了数值计算,同线性k-ε模型的计算结果和实验结果进行了比较.结果表明,这种非线性k-ε模型同线性模型相比有较大的改进,同已有的其它非线性模型相比,也有一定的优点.     

肋角度对流量系数影响的数值模拟研究

利用有限体积法对三维不可压缩的N-S方程进行离散,对同时带肋和气膜孔出流的内流冷却通道在45°,60°,90°,120°肋的结构下进行了数值模拟。网格划分采用非结构化网格,湍流模型为Realizablek-ε模型,近壁处湍流采用壁面函数,采用SIMPLE算法求解速度与压力的耦合。计算获得了同时带肋和气膜孔出流的内流冷却通道在不同肋角度时的三维流场分布和气膜孔的流量系数值。结果表明同时带肋和气膜孔出流的通道流场结构比较复杂,肋在通道中诱发的二次流改变了气体在进入气膜孔时偏转角度,从而影响了流量系数值。通道中气膜孔的流量系数值随肋角度的增大而呈增大的趋势。计算结果与实验结果进行了比较,两者规律符合良好。      

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

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

多尺度湍流模型在圆柱绕流中的应用

使用基于可变时间间隔平均方法的多尺度湍流模型,对高超临界区雷诺数下非定常圆柱绕流进行了数值模拟.结果表明:多尺度湍流模型在升阻力系数预测上与实验值一致,而标准k-ε模型和SST(shear stress transport)k-ω模型却只能兼顾其一;在旋涡脱落特性预测上与实验值相一致,其误差小于SST k-ω模型的5.8%和标准k-ε模型的37.6%;在表面平均压力系数分布预测上其误差仅为3.6%,明显优于SSTk-ω模型的13%和标准k-ε模型的53.7%;在表面摩擦因数分布及分离角度方面与实验结果相吻合,且分离角度误差仅为0.78%,结果优于标准k-ε模型的1.04%和SSTk-ω模型的1.83%,充分验证了该湍流模型应用于复杂湍流模拟预报的潜力.     

五种湍流涡粘模型在二维方柱绕流数值模拟中的对比研究

为研究不同雷诺时均Navier-Stokes(RANS)模型在求解钝体绕流场的差异,采用五种两方程湍流涡粘模型计算了二维方柱的绕流场(雷诺数Re=22 000),得到了不同湍流模型在不同网格离散方案下的方柱气动力与流场特性。结合以往文献数据,通过对比不同湍流模型计算结果的差异,揭示了不同湍流模型的特点。研究结果表明:网格离散方案对方柱绕流数值模拟结果有重要影响;RNGk-ε湍流模型的计算效率最高,SSTk-ω湍流模型的计算效率最低;Standard k-ε湍流模型的计算结果准确程度整体弱于其余湍流模型;RNGk-ε湍流模型、Realizable k-ε湍流模型与Standard k-ω湍流模型的计算结果大致相当,较接近大涡模拟结果;SSTk-ω湍流模型的模拟结果优于其余湍流模型,其尾流区速度场与试验结果吻合较好;两方程湍流模型计算的二维方柱绕流结果与试验结果以及大涡模拟结果相比,存在不可忽略的差异。     

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

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

两种k-ε模型用于扩压器内流数值模拟的比较

同时采用标准k-ε及非线性k-ε模型在附体坐标系下求解定常不可压N-S方程组计算扩压器内湍流流场, 得到了扩压器的压力恢复系数和速度分布, 与实验结果对比发现非线性k-ε模型预测的结果更吻合实验值, 且能测出临界扩压角下的局部分离流。文中着重研究了非对称扩压器, 并给出了用两种模型计算的湍流量。      

膨胀循环推力室再生冷却换热的数值模拟

为了解液体火箭发动机膨胀循环推力室再生冷却换热特性,采用数值模拟方法,研究了冷却剂流动方式、推力室圆柱段长度、圆柱段室壁加肋和气壁面粗糙度等因素对冷却通道压降、冷却剂温升、壁面热流密度和温度分布等换热特性的影响.计算过程中采用k-ε双方程湍流模型.计算结果表明:采取顺流冷却要比逆流冷却的冷却通道压降低,但同时冷却剂温升也低;对于室壁加肋结构,在肋个数相同而只改变肋高度的情况下,总换热量正比于总换热面积.      

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

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

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

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

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

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

涡轮叶片内部沿周边分布的竹节孔换热特性研究

根据发动机涡轮叶片内部沿周边分布的竹节孔冷却通道结构建立了简化的传热分析模型。首先通过无量纲分析得到了竹节孔的平均努塞尔数不仅与雷诺数和普朗特数有关,还与“竹节”形状、通道外壁形状以及流体和固体的导热系数比等参数有关。然后,对雷诺数和冷却通道几何参数对竹节孔平均努塞尔数的影响进行了数值模拟研究,计算并对比了上述参数对平均努塞尔数的影响规律,结果表明：（1）“竹节”结构使平均努塞尔数显著增大,而通道外壁形状主要影响了局部热流密度和局部努塞尔数的分布情况;（2）随着通道内的雷诺数和“竹节”的节高/孔径比的增加,平均努塞尔数单调增大。但以增大节高/孔径比的途径来提高平均努塞尔数时,阻力系数迅速增大;（3）可通过增大雷诺数和优化节距/节高比来提高竹节孔的平均努塞尔数,计算表明最佳节距/节高比约为10。     

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

直接供气预旋转静系流动和换热数值模拟

采用CFD商业软件FLUENT分别从维度模型和湍流模型两个角度对接近真实的高压涡轮旋转盘腔进行数值模拟,研究了一种直接供气预旋转静系旋转盘腔模型内的流动和换热特性.通过计算研究发现:低位腔内的流动仅由中心入流控制,二维计算不能反映高位腔内中心入流与预旋入流的掺混;二维和三维计算所得的腔内压力分布相似,实验值介于两者之间;主盘面的换热主要受中心入流控制,二维和三维计算差别不大,Realizable k-ε湍流模型在定性和定量上均能更好地反映换热实验结果.      

叶片内部蒸汽冷却的数值模拟

在某燃气轮机第二级导向叶片的基础上,对叶片内部进行了设计并采用内部蒸汽冷却方式.采用非结构化网格和Realizable k-ε紊流模型,对同时带15°,30°和45°肋的叶片内外流场以及叶片本身进行了三维热耦合的定常数值模拟.热耦合采用在流固耦合面的热对流和固体导热来实现.内部冷却介质选用比空气热容大,粘度低的水蒸汽.结果表明:带肋的通道流场非常复杂,扰流肋的存在使各壁面的换热都得到了增强,相对于光滑通道换热显著,有明显的冷却痕迹,冷却效果明显;采用蒸汽冷却比用空气冷却效果要好;压力面温度普遍低于吸力面温度.