用DES数值模拟具有横向喷流的紊流流场

DES(DettachedEddySimulation)是把RANS(雷诺平均Navier Stokes方程)方法及LES(大涡模拟)方法结合起来的模拟有脱体涡的紊流流场的数值模拟方法。其主要思想是在物面附近解雷诺平均Navier Stokes方程、在其他区域采用Smagorinski大涡模拟方法。本文用DES及混合非结构网格数值模拟复杂紊流流场,算法采用Osher逆风格式,对超音速钝头体、三角翼、某型飞机及带喷流的导弹绕流用该程序进行了数值模拟。结果表明用DES模拟有分离的紊流流场是高效的。     

细长旋成体背风面非对称流场的DES模拟

利用脱体涡模拟DES方法对细长旋成体非对称绕流进行数值模拟,比较了基于S-A湍流模型的DES与RANS方法对背风面非对称分离涡的数值模拟能力,发现在大迎角非对称涡未破裂情况下,DES与RANS均能模拟出与实验相符的非对称绕流;当迎角增大,背风面非对称涡发生破裂时,RANS无法准确捕捉到背风面流场的非定常性,而DES能准确预测非对称涡的飘起与破裂,并与实验值接近。计算结果表明:与RANS方法相比,DES方法具有更好的模拟大攻角,大分离流动的能力,尤其在非对称流场的分离涡破裂模拟方面具有明显优势,能够更真实模拟出细长旋成体背风面分离涡破裂之后的非定常流动特征。     

一种改进的类DES湍流模拟方法

构造了一种基于一方程S-A(Spalart-Allmaras)模型和一方程Yoshizawa亚格子模型的混合RANS/LES(Reynolds-averaged Navier-Stokes/large eddy simulation)湍流模拟方法.在涡黏假设的基础上,将Yoshizawa亚格子湍动能方程转化为等效的亚格子湍流涡黏性输运方程,并采用混合函数将其与S-A模型方程进行混合,从而改进了DES(detached eddy simulation)模型的亚格子行为,同时克服了其依靠网格控制模型转换的缺点.模拟了超声速的带斜坡凹腔流动,并与相同网格下的LES及DES结果进行了比较,结果表明该混合RANS/LES方法在远离壁面的自由剪切流区域与LES行为一致,而在附着边界层区域表现优于LES和DES方法.      

基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

基于分离涡模拟的起落架气动噪声研究

用基于SA湍流模式的非定常RANS、分离涡模拟(DES)和延迟分离涡模拟(DDES)分别对四轮基本起落架模型进行了数值模拟,并根据所得的非定常流场计算了表面声压级分布和声压谱。三种方法所用的时间成本大致相同,URANS略低于其它两种,而DDES由于在附着流动区更好地保持RANS特性,故时间成本略低于DES。计算显示,非定常RANS在附着流动区能够得到合理的结果,但不能准确刻画分离流动的流动形态。DES和DDES都能较好地刻画起落架绕流的定常和非定常特性,DDES的结果略好于DES。因为起落架流动属于大分离流动,所以DES也能够得到相对正确的结果。研究结果验证了分离涡模拟在起落架大分离非定常流动预测与噪声预测中的可行性,对减小起落架噪声的方法研究具有一定的意义。     

部分平均Navier-Stokes方法模拟超声速斜面空腔流动

采用部分平均Navier-Stokes (PANS)方法对超声速斜面空腔流动进行了数值模拟研究,评估其对于超声速非定常湍流模拟的性能,并与雷诺平均N-S(RANS),脱体涡模拟(DES)的计算结果及实验数据进行对比.研究表明:模化湍动能比例全场可变的PANS方法预测的速度、壁面压力和摩擦力系数分布与DES的结果非常接近,都与实验值吻合得较好,优于RANS的计算结果;非定常雷诺平均N-S(URANS)计算的自由剪切层近乎为二维稳态的,而DES和PANS方法可以求解出更为丰富的流动结构;模化湍动能比例全场统一的PANS方法虽然可以求解出相比RANS更多尺度的流动结构,但在近壁区不能回归到RANS模型,预测的湍流边界层的速度型偏离对数律,后续的流动计算也偏离实验数据.      

改进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工况下没有明显优势；并对这一现象做出了合理分析。      

基于DES方法的轮胎部件外流场数值模拟研究

采用基于S-A湍流模型的脱体涡(DES)方法,数值模拟研究了三维轮胎部件外流场。通过非定常计算,将所得结果与相关文献对比分析。结果显示,该方法能够准确地模拟轮胎的非定常特性,后轮处能够很好地捕捉马蹄涡,其范围与文献接近,且起轮胎压力系数图能够和文献基本相似。     

分离涡模拟方法在内埋弹舱机弹分离模拟中的应用

为了研究内埋弹舱非定常流场对武器投放分离特性的影响,基于分离涡模拟（DES）方法和非结构重叠网格技术,建立了内埋武器机弹分离高精度数值模拟方法。在此基础上,首先通过典型空腔流动和机弹分离算例验证了数值方法的有效性,然后采用平均雷诺应力模拟（RANS）和DES两种方法探究了机弹分离过程中弹舱腔内涡结构的演化历程,分析了两种方法在机弹分离流场结构模拟中的差异,研究了内埋弹舱流场结构对导弹分离特性的影响。计算结果表明：机弹分离过程中,高强度涡结构逐渐向弹舱腔体后缘聚集,同时弹舱剪切层流动结构被导弹破坏,导弹出弹舱阶段俯仰方向所受气动力矩急剧变化。相比RANS方法,DES方法捕捉到的高强度小尺度涡结构对导弹压力分布影响较大,使得分离后期导弹俯仰角、俯仰角速度产生明显差异,DES方法更适合内埋武器舱分离特性精细化研究。该研究结果可为内埋武器机弹分离研究的计算方法选择和安全分离策略提供一定的参考依据。     

基于Ansys Fluent的近场翼尖涡数值模拟与分析

为了进一步研究飞机远场尾涡,提供网格分配及湍流模型的参考,并为整机模拟提供必要的参考依据,通过基于Ansys Fluent的数值模拟方法,研究了NACA0012机翼的近场翼尖涡流场,采用有限体积法求解不可压缩雷诺平均Navier-Stokes方程,其中雷诺应力项分别以S-A和Realizable k-ε模型封闭,模拟了近场翼尖涡卷起的过程,分析了机翼表面压力以及涡核参数,包括轴向涡量、涡核位置、涡核粘性等,并与风洞实验结果进行了对比。结果分析表明:基于局部O-网的六面体网格,RKE模型要优于S-A模型,与实验值更为吻合。     

高雷诺数槽道湍流的壁面模化大涡模拟研究

选取基准壁湍流的槽道流动,研究了多种模型的壁面模化大涡模拟.模型包括经典的大涡模拟、Spalart-Allmaras、分离涡模拟和一种动态混合模型.基于摩擦速度的雷诺数范围为395~12000,采用3组粗糙网格,流向和展向维数分别同取37,49和65,法向维数保证y+(1)~1.主要研究平均速度、雷诺切应力分布、详细分析了各模型的特性差异并展示了相应的湍流结构.研究表明:在高雷诺数粗糙网格下,大涡模拟失去求解精度,分离涡模拟出现对数律不匹配,动态混合模型的计算接近直接数值模拟,其对数率区可解应力约占雷诺切应力的93%,边界层外层可解应力约占99%.这说明合适的混合模型可以在经济成本下保证计算精度,具有解决实际问题的潜力.      

使用RANS/LES混合方法对钝前缘三角翼进行数值模拟(英文)

阐述了建立在一方程湍流模型基础上的分离涡模型(DES)和延迟脱体涡模型(DDES),它们很好地模拟了壁面附近的小尺度流动。该方法在大攻角分离流区域,使用了亚格子应力的Smagorinsky大涡数值模拟,模拟亚音速旋涡流动。使用大攻角下的钝前缘65°大后掠三角翼模型,研究了CFD对涡的发展、破裂以及复杂涡的发展和演变的模拟能力。研究结果说明了方法对模拟分离流动是可行的。     

湍流模型对喷流计算精度的影响研究(英文)

采用三维RANS代码MBNS3D对包含了内流、外流及后体干扰的三种简化尾喷流,进行了数值模拟研究。通过壁面压力分布及喷流截面速度型与实验数据的对比,研究了SA和SST两种湍流模型及其修正(SACC,SADES,SSTCC)对喷流计算精度的影响。本文计算表明,SSTCC模拟效果最为精确;利用SA模型进行脱体涡模拟,在同样的网格下优势不明显。湍流模型可压缩修正改善了模拟喷流的膨胀趋势,但对附面层分离位置的预测,在本文的工作中,还不如原始模型。两类模型相比,SST模型的结果好于SA模型。     

引入压力梯度的针对角区分离流动的Spalart-Allmaras模型的改进

针对Spalart-Allmaras（S-A）模型在角区分离计算中的问题,将无量纲的压力梯度引入其涡黏性输运方程的生成项,得到了改进的S-A模型.通过对两套含角区分离的低速压气机叶栅进行验证计算发现:与实验结果相比,原始S-A模型所得的分离区偏大,分离区内壁面压力偏低;而改进模型得到了与实验一致的分离区尺寸以及吸力面、压力面压力系数分布等结果.针对S-A模型涡黏性生成项和耗散项的分析表明:引入的无量纲压力梯度有效的识别了角区分离,在分离区内改变了涡黏性的生成、耗散关系,增大了涡黏性,从而缩小了计算所得分离区,同时在主流区保留了原始S-A模型的计算结果,进而带来了良好的改进效果.      

A study on turbulence transportation and modification of Spalart–Allmaras model for shock-wave/turbulent boundary layer interaction flow

It is of great significance to improve the accuracy of turbulence models in shock-wave/ boundary layer interaction flow. The relationship between the pressure gradient, as well as the shear layer, and the development of turbulent kinetic energy in impinging shock-wave/turbulent boundary layer interaction flow at Mach 2.25 is analyzed based on the data of direct numerical simulation(DNS). It is found that the turbulent kinetic energy is amplified by strong shear in the separation zone and the adverse pressure gradient near the separation point. The pressure gradient was non-dimensionalised with local density, velocity, and viscosity. Spalart–Allmaras(S–A) model is modified by introducing the non-dimensional pressure gradient into the production term of the eddy viscosity transportation equation. Simulation results show that the production and dissipation of eddy viscosity are strongly enhanced by the modification of S–A model. Compared with DNS and experimental data, the wall pressure and the wall skin friction coefficient as well as the velocity profile of the modified S–A model are obviously improved. Thus it can be concluded that the modification of S–A model with the pressure gradient can improve the predictive accuracy for simulating the shock-wave/turbulent boundary layer interaction.     

基于自适应混合网格的脱体涡模拟

基于混合网格和CGNS（CFD General Notation System）数据结构,建立了一种各向同性加密/稀疏的网格自适应方法。在悬空点的后处理中,让含有悬空点的单元转化为任意多面体,从而简化了自适应单元剖分模版,同时自适应网格单元之间可完全相容,自适应生成的网格能够直接用于可处理任意多面体的流场求解器。将该自适应方法与脱体涡模拟（DES）算法相结合,开展了65°后掠三角翼大迎角流动的数值模拟应用,并与初始网格的模拟结果进行了详细比较。对比表明：采用网格自适应方法适当增加局部网格量,能够以较小的成本迅速提高三角翼背风区的空间分辨率,增强数值模拟对小尺度涡系结构的解析能力,从而弥补了基于混合网格的脱体涡模拟中常用二阶格式计算的空间分辨率相对偏低、不利于湍流多尺度结构精细模拟的不足。     

基于DES模型的超燃气动斜坡/燃气发生器方案仿真研究

采用基于剪切应力输运(SST)k-ω两方程湍流模型的分离涡(DES)方法和雷诺时均Navier-Stokes(RANS)方法对超燃气动斜坡结合燃气发生器增进掺混方案进行数值仿真研究。通过对比DES仿真结果、RANS仿真结果与试验结果,发现DES仿真结果对流场内涡结构的捕捉和分辨能力强于RANS方法获得的结果。选取燃气发生器喷流形成的壁面涡流特征角为比较标准, DES仿真结果显示该角度为48°,RANS仿真结果为37°,而试验测得该角度为47°。可见DES仿真结果与试验结果更为接近,说明对于气动斜坡结合燃气发生器的超燃掺混方案,DES方法仿真结果是合理的。