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

基于SST k-ω的超声速气膜冷却湍流模型温度修正

针对现有湍流模型无法准确预测非等温超声速气膜冷却行为的问题,在现有SST k-ω可压缩修正模型基础上,以总温梯度为变量,完成了湍流模型的温度修正,并首先通过非等温可压缩自由剪切流动实验数据初步验证了其修正效果,在此基础上对温度修正模型预测超声速气膜冷却传热的准确性进行了验证。结果表明,基于剪切层总温变化的湍流模型修正效果显著,可准确预测大温度梯度下的自由剪切流动轴向速度分布;修正模型计算得到超声速气膜冷却壁面热流分布与对应的实验结果吻合;当用于剪切层温度大梯度变化的超声速气膜冷却数值模拟时,温度修正后的SST k-ω模型与可压缩修正的k-ω模型、SST k-ω模型相比,具有显著优越性。     

典型航空分离流动的雷诺应力模型数值模拟

先进飞行器设计对CFD方法的边界层分离模拟能力提出了更高要求,传统雷诺平均纳维-斯托克斯(RANS)涡黏性模型因建模和构造层面的理论缺陷导致其分离流动预测可信度较低。雷诺应力模型由于未对雷诺应力及其分量关系进行建模,对湍流非平衡、旋转以及雷诺应力各向异性等流动现象具有天然的理论优势。为验证与确认雷诺应力模型对典型航空分离流动的预测能力,基于SSG/LRR-g模型,以NACA4412翼型大攻角分离、M6机翼跨声速分离以及F6翼身结合区分离流动为例,探讨了雷诺应力模型对逆压梯度、激波诱导分离、二次流动分离等典型航空分离流动预测的适应性。通过与k-ω剪应力输运(SST)模型模拟结果对比发现,雷诺应力模型对分离泡大小、速度型分布、雷诺应力分布和激波位置等关键特征的模拟精度较涡黏性模型显著提升,基本验证了雷诺应力模型可在翼身接合区角区流动和三维强激波诱导分离等问题中得到正确的流动特征,而SST模型在此类流动中基本失效,显示了雷诺应力模型在典型航空分离流动中较涡黏模型的优势。同时,发现k-ω SST模型所包含的Bradshaw假设在三维激波诱导分离较强时严重影响了模型预测的准确性,是预测结果偏离试...     

Pω增强型k-ω湍流模型在三角翼旋涡流动的应用

在三角翼旋涡绕流数值模拟中,标准 Wilcox k-ω湍流模型生成项未考虑旋度的影响而导致预测的旋涡强度较弱。通过引入探测因子区分剪切层和涡核,在旋涡流动的高旋度区域增加ω方程生成项的方法,基于结构化网格上的 RANS 求解器,加入了 Pω增强型 k-ω湍流模型,对绕尖前缘三角翼亚声速和跨声速旋涡流场进行了数值模拟。计算结果与 NASA 的 NTF 风洞和 DLR 的 DNW-TWG 风洞试验数据进行了对比分析,结果表明：不论在亚声速还是跨声速自由来流条件下,Pω增强型 k-ω湍流模型计算的压力分布、涡破裂位置均与试验数据吻合良好,准确地预测出了三角翼上翼面的主涡、二次涡结构,特别是跨声速条件下激波干扰导致的涡破裂的临界迎角及涡破裂位置,表明 Pω增强型 k-ω湍流模型在绕三角翼旋涡流动数值模拟中具有良好的适用性。     

Sajben跨声速扩压器分离流动中湍流模式数值研究

为研究湍流模式对激波/湍流边界层干扰内流流动的影响,提高数值计算准确度,使用SA,SST k-ω,非线性EASM k-ω,Gao-Yong四个湍流模式对Sajben扩压器内激波/湍流边界层干扰流动进行了数值计算。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。文中展示了四个湍流模式计算得到的壁面压力、速度剖面、摩阻系数等分布。计算值与实验值符合很好,四个湍流模式总体上能够较好地模拟扩压器内激波/湍流边界层干扰复杂分离流动。Gao-Yong湍流模式对分离区内的压力、速度型的模拟更加准确,而非线性EASM k-ω模式对分离再附点位置计算最理想。     

翼伞平面形状对翼伞气动性能的影响

对带气室的展弦比为3的不同平面形状翼伞模型的流场进行了三维定常数值模拟,详细考察了平面形状对翼伞气动性能的影响.运用有限体积法对三维坐标系下不可压雷诺时均Navier-Stokes(RANS)方程进行了直接求解,采用剪切应力输运(SST)k-ω二方程湍流模型进行湍流模拟.数值模拟得出的原始翼伞的气动性能参数与试验数据在...     

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

使用基于可变时间间隔平均方法的多尺度湍流模型,对高超临界区雷诺数下非定常圆柱绕流进行了数值模拟.结果表明:多尺度湍流模型在升阻力系数预测上与实验值一致,而标准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%,充分验证了该湍流模型应用于复杂湍流模拟预报的潜力.     

无网格法耦合RNG k-ε湍流模型在亚、跨声速翼型黏性绕流中的数值模拟

基于移动最小二乘无网格方法,耦合RNG(Re-Normalisation Group)k-ε湍流模型求解雷诺平均Navier-Stokes方程。采用AUSM(Advection Upstream Splitting Method)+-up迎风格式求解数值通量,应用在高度各向异性点云结构中取得良好结果的点云重构技术结合移动最小二乘法拟合空间导数,并用三阶SSP(Strong Stability Preserving)型Runge-Kutta显式时间推进格式求解离散后的控制方程。在此基础之上,实现了对NACA0012、RAE2822翼型亚、跨声速黏性绕流的数值模拟,给出了翼型表面压力系数分布曲线、不同位置处的平均速度剖面、马赫数等值线等计算结果,并与实验值及相关文献数值模拟结果进行比较,结果吻合较好。表明所发展的结合点云重构技术的无网格方法耦合RNGk-ε湍流模型能够成功模拟翼型亚、跨声速黏性绕流,验证了所提算法的有效性,并拓展了无网格方法求解湍流流动的途径。     

湍流模型对湍流射流CFD模拟的影响

横流中的湍流射流涡旋结构复杂,湍流模型对其模拟结果影响很大。采用S-A,Realizable k-ε,k-ω和SST k-ω四种湍流模型对吹风比为0.5和1.5工况下的横流中射流进行了数值模拟研究。与实验值比较验证了数值模拟结果的可信性。分析了旋转张量和剪切张量对湍流生成的影响。发现四种湍流模型对下游速度分布的预测均较为准确,在小吹风比下S-A模型性能最好;在大吹风比下k-ω模型较为准确。     

基于湍流模型的转捩预测方法研究

阻力的精确计算,一直是低雷诺数飞行器设计的关键.而与附面层分离和转捩相关的摩擦阻力的精确计算则是此类飞行器阻力精确计算的关键.因此,开展了基于湍流模型的转捩预测方法的研究,应用几种湍流模型对平板边界层进行计算分析,研究各湍流模型预测流动转捩的能力.由于各模型对扰动的过于敏感,引起转捩提前发生,使预测与实验差距较大,由于k-ωSST两方程湍流模型预测能力优于其它模型,在此模型基础上引入了间隙函数进行粘性系数修正,对平板和NACA 0012翼型的转捩重新进行了预测,计算结果表明改进后的模型对转捩位置具有较好预测能力.     

燃烧室-涡轮耦合流动传热超大涡模拟研究

为了准确预测发动机燃烧室和涡轮耦合复杂气热环境下涡轮部件的流动和换热特性,应用基于BSL k-ω湍流模型的高精度超大涡模拟方法(VLES),以及SST k-ω雷诺平均湍流模型对雷诺数为Re=3.8×105的高压涡轮导叶在均匀进口条件以及燃烧室-涡轮耦合情况进口条件下的涡轮流动和传热进行了数值研究.耦合火焰面生成流型(F...     

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

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

超音速子母弹分离干扰流场数值模拟与试验研究

采用风洞试验和数值模拟相结合的方法对超音速子母弹分离干扰流场特性进行了研究。对比数值模拟与试验研究结果表明,结合k-ωSST两方程湍流模型求解雷诺平均N-S方程,并以AUSM+作为离散格式的数值方法,可较为细致地模拟子母弹分离干扰流场。通过对子母弹特定状态干扰流场结构的分析和对其流动机理的探讨,发现由于激波相互干扰使弹体表面压力分布多变、子弹攻角增大时纵向干扰量大影响子弹分离姿态。本文的工作为进一步研究子母弹分离干扰流场流动机理奠定了基础。     

二维地效翼波浪面失速特性研究

为了研究地效翼在波浪面下的失速特性,采用FLUENT商业软件求解非定常不可压N-S方程和k-ω/SST湍流模型,对二维NACA0012翼型在波浪地面下的流动进行了数值模拟分析。通过与平地面下地效翼失速特性进行对比分析,揭示了机翼在不同波浪地面下的失速特性。研究结果为地效飞行器的设计和优化提供了理论依据。     

进气道中激波-边界层干扰特性研究

针对高超进气道中激波/边界层干扰等现象,选取了24°压缩拐角及二维高超进气道DLR-GK01等算例,比较研究了混合网格并行数值模拟软件中SA、SST、TNT三种湍流模型对激波/边界层干扰等现象的模拟能力,主要考察了不同湍流模型对流场结构、壁面压强以及壁面摩阻的影响。结果表明:三种湍流模型在模拟大分离流场时存在不同的流动分离特性,TNT湍流模型较SST和SA湍流模型模拟的分离起始位置、分离区大小及边界层厚度与实验吻合度最好。     

前后翼型装药结构SRM气固两相流动特性

基于可压缩流动强守恒型Navier-Stokes方程,运用颗粒轨道模型(PTM)和shear stress transfer(SST)k-ω湍流模型,采用计算单元内颗粒源法(PSIC)进行气固两相耦合计算,建立了某型前后翼型装药结构固体火箭发动机(SRM)工作初期的三维两相内流场数值计算模型.对比分析了纯气相和气固两相...     

基于IDDES方法的翼型结冰失速分离流动数值模拟

应用基于k-ωSST湍流模型的IDDES(Improved Delayed Detached Eddy Simulation)方法,就失速点附近翼型前缘典型双角状积冰导致的复杂分离流动进行了数值模拟研究.通过与风洞试验结果进行对比,表明对于此类分离流动问题,IDDES方法能够在壁面附近取得良好的速度预测结果,有效解析分离区域内的中小尺度湍流结构,较为准确地描述大尺度时均分离泡的再附位置和形态特征,适用于翼型结冰后复杂流动的精细分析.同时计算结果显示当此带冰翼型位于失速点附近时,角状冰后方脱落剪切层内部的旋涡不稳定析出和输运过程促进了外部流动与回流区域流动间的掺混,将导致流动发生非定常再附现象.     

新型单方程湍流模型构造及其应用

为了提高雷诺应力本构关系式对于非平衡湍流的预测精度并且兼顾求解效率,发展了一种基于湍动能k的单方程(KDO)湍流模型。其主要思路为:采用平板直接数值模拟(DNS)数据对原始Bradshaw假设进行重新标定,使得当地湍动能和雷诺主应力之比能够根据当地流动条件进行自适应调节;同时,对标准k-ε模型中的湍流耗散率输运方程采用代数形式进行模化,进而形成一种一方程湍流模型。算例结果表明:KDO湍流模型对于对数率能够准确反馈,而在带有激波或部件干扰等流动现象的RAE-2822、ONERA-M6和DLR-F6算例中,KDO湍流模型能够准确控制湍动能的增长和衰减,相比于Spalart-Allmaras和Menter k-ω剪切应力输运(SST)模型,KDO湍流模型的计算结果有了较为明显的改善。     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

基于非结构/混合网格的脱体涡模拟算法

为了提高二阶精度有限体积算法的湍流数值模拟能力,在原始Roe格式基础上建立了与脱体涡模拟(DES)方法相匹配的二阶混合耗散自适应格式,根据流场信息自动调节格式耗散,并分别基于Spalart-Allmaras一方程湍流模型和k-ω剪切应力输运(SST)两方程湍流模型,发展了基于非结构/混合网格的DES方法。采用该方法计算了雷诺数为3 900的圆柱绕流和NACA 0021翼型60°大迎角分离流两个典型算例,通过与试验数据以及其他数值结果的对比验证了该方法的可行性。同时开展了不同数值格式、湍流模型的对比分析,研究结果表明:采用混合格式的DES算法能够解析更小尺度的湍流涡结构、计算数据更接近试验值;本文的DES类算法受其基准湍流模型影响较小。