湍流模型在单/双缝引射矢量喷管中的验证分析研究

利用数值模拟方法,分别求解以SA、SSTk-ω、EASMk-ω或k-ζ湍流模型封闭的RANS方程,针对单/双缝引射的激波控制射流推力矢量喷管展开研究,在多个主喷流压比NPR(4.6,7.0,8.78,10.0)和次主流压比SPR(0.7,1.0)下,系统考察了四种不同湍流模型对射流推力矢量喷管内外流场的流场结构,性能参数及主喷管内壁面压力分布的预测能力,探讨了由激波控制产生的推力矢量随不同参数的变化规律。数值计算结果表明,四种湍流模型都能比较准确地预测出射流推力矢量喷管的性能参数和喷管流场结构,双缝引射的推力矢量效率并不一定优于单缝引射。     

改进S-A湍流模型对横向射流的CFD模拟

通过在Spalart-Allmaras(S-A)湍流模型中加入应变率张量，来改善模型对横向射流(JICF)的计算精度。模拟了速度比分别为0.5和1.5工况下的横向射流。经实验值验证，结果表明:改进后的S-A模型预测精度有所提升；并对湍流模型改进部分对横向射流的作用机制做出了分析。将改进后S-A模型的计算结果与k-ε模型，雷诺应力模型(RSM)和大涡模拟(LES)的结果进行比较，发现S-A模型在速度比为0.5时精度最高，而在速度比为1.5工况下没有明显优势；并对这一现象做出了合理分析。      

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

影响冷气掺混数值模拟精度的若干问题分析

对冷气掺混数值模拟中的计算网格、湍流模型、射流边界条件等影响精度的若干问题进行分析。在吹风比0.5下,对不同网格分布、湍流模型、射流边界条件进行数值模拟,所得结果与Ajersch的实验结果进行对比。结果表明,出口下游的网格分布可适当稀疏以减少计算量;在所研究的几种湍流模型中,k-ε模型所得结果与实验结果吻合得最好;考虑射流通道内流动能提高精度,在射流出口给定流场分布也可保证计算结果的精度。     

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

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

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

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

湍流度和吹风比对叶片前缘双出口孔射流冷却效率影响

为了探讨湍流度对一种新型气膜孔射流气膜冷却影响,利用商业软件提供的有限体积法求解N-S方程,对湍流度分别为0.4%、10%和20%时的双出口孔射流冷却效率进行数值模拟。吹风比变化范围为0.5到2.0。首先将圆柱孔射流冷却效率计算结果和实验数据进行了对比,二者吻合较好。计算得到了双出口孔射流气膜冷却下的流场、径向平均冷却效率。结果表明,湍流度和吹风比对冷却效率都有较大影响。湍流度为0.4%和10%时,最高冷却效率在吹风比1.0时获得;湍流度30%时,最高冷却效率在吹风比2.0时获得。吹风比0.5和1.0时,冷却效率随着湍流度的增加而降低;吹风比2.0时,冷却效率随着湍流度的增加而增加。     

一种超声速射流计算方法与实验验证

高温射流流场计算是尾焰辐射目标特性计算的前提,然而由于缺少可靠实验数据,针对湍流超声速射流的数值模拟多集中于低温射流,高温射流计算与实验的对比工作还很少见。利用k-ωSST双方程湍流模型,模拟了多个典型超声速射流实验的流场速度与温度分布,通过与实验结果进行对比,建立了一种超声速射流计算方法。首先,通过对比多个低温射流的实验与计算结果,探索了湍流模型中可压缩修正以及来流湍动粘性比对超声速射流计算结果的影响;进而,针对火箭发动机尾焰实验,计算尾焰流场与流场红外辐射,流场辐射计算结果与实验观测结果符合较好,进一步验证了计算方法。最终认为经过可压缩修正的k-ωSST双方程湍流模型结合湍动粘性比取值30可以作为超声射流计算中较为典型的湍流计算方法。     

不同湍流模型对插板式进气畸变的数值模拟

为得出一种较为准确的湍流模型,应用三种湍流模型对插板式进气畸变进行数值模拟,通过与畸变实验实测数据的对比,研究不同湍流模型在插板式进气畸变计算中的预测能力。不同湍流模型计算出的出口流场畸变不同,重整化群k-ε模型得出高压区范围较大,S-A模型过渡区域范围大。结果表明:插板高度为畸变发生器内径的48%时,低压区范围角与实验值较为吻合,插板高度为畸变发生器内径的20%时,低压区范围角偏小。畸变指数的计算值以可实现化k-ε模型最为准确。     

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

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

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

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

两种k-ω型湍流模型在无网格方法中的应用研究

基于在高度各向异性点云结构中取得良好结果的点云重构技术,采用移动最小二乘无网格法解耦计算了k-ω型湍流模型方程和雷诺平均N-S方程(RANS),实现了对湍流流动的数值模拟。通过对NACA0012翼型、RAE2822翼型的绕流以及后向台阶大分离流动的数值模拟,比较分析了k-ωSST(Shear Stress Transport)和k-ωTNT(Turbulent/Non-Turbulent)两种模型的湍流预测能力。结果表明:对弱逆压梯度情况下的湍流,两者具有相当的预测能力;而k-ωSST较k-ωTNT具有更好的激波捕捉能力。对于大分离流动,两种模型较实验数据均存在一定差异,但总体趋势吻合良好。     

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

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

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

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

S弯扩压器中四种湍流模型的比较

采用Realizableκ-ε,S-A(Spalart-Allmaras),SST(Shear stress transport剪切应力输运),RSM(Reynolds stress equation model雷诺应力模型)四种湍流模型对S弯扩压器内的分离流动进行了模拟,分析了不同湍流模型对计算结果的影响.将不同轴向位置的表面压力系数、分离区的大小、马赫数等三个量的实验值与计算值进行了比较,结果表明:分离区湍流模型对表面压力分布影响较大,未分离区不同湍流模型模拟精度相当;S-A模型对分离区的模拟较保守,RSM模型对分离的模拟精度较高,SST与RSM模拟结果很接近;各湍流模型计算的高速区域均比实验值大.扩压器的平均总压恢复系数受湍流模型影响不大,和实验结果吻合较好.      

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

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

湍流模型源项的数值处理方法对计算结果的影响研究

涡粘性湍流模型广泛应用于各种工程湍流问题的计算。对于大多数湍流模型,在湍流控制方程的右端会出现源项,其刚性给数值计算带来很大的影响。从源项的物理意义出发分析了源项导致求解困难的原因,并在通用的求解湍流模型方程的点隐法基础上,以源项弱刚性的S-A湍流模型和源项强刚性的k-ωSST湍流模型为例,论述了生成项和耗散项之间的平衡关系对数值模拟的重要意义,给出了针对不同源项的具体处理方法。对RAE2882翼型跨声速流动算例的模拟结果表明,这些处理方法有效的提高了控制方程组的计算稳定性。     

涡轮叶片前缘气膜冷却数值模拟

对具有双排气膜孔的涡轮叶片前缘进行了气膜冷却数值模拟,在吹风比M为0.5,1.0的情况下得到了相邻两排孔的流场和气膜冷却效率的分布规律,并将计算结果与实验进行了对比分析.结果表明:①在不同的吹风比下,两排孔的冷气射流运动非常复杂,冷却效率沿流向的分布都是非线性的;在M为1.0时,25°的小角度喷射气膜冷却效果较好.②利用两层k-ε湍流模型比传统的k-ε湍流模型数值预测的精确度高,可较好地模拟叶片前缘的气膜冷却特性.      

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

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