SAS与PANS模型在圆柱绕流中的应用比较

分别运用SST(Shear Stress Transport)、SST-SAS(Scale-Adaptive Simulation)、两种变fk(模化湍动能的比例)函数的SST-PANS(Partially Averaged Navier-Stokes)湍流模型对Re=3900的圆柱绕流进行了数值研究,重点从湍流结构捕捉、气动力计算、涡黏性控制等方面,比较了SAS与PANS两类RANS/LES混合模型的计算能力,并通过不同网格计算分析了模型的网格敏感性。数值结果表明:SAS及两种变fk方法的PANS模型均具有求解小尺度涡运动的能力,并能较好地反映出绕流尾迹的三维非定常特性,同时PANS模型能捕捉到更多的非定常结构;SAS模型中自适应尺度Lvk立足于当地流动,对网格依赖较小,计算的湍动粘度分布更合理,能够更好地计算剪切层及回流区;两种PANS模型网格独立性较差,出现了雷诺应力不足的现象;类DES可变fk函数构造相对简单,所得fk分布更准确,使用tanh函数计算的尾迹区fk值偏低,对流场调控能力稍差。     

基于RANS/LES方法的超声速底部流场数值模拟

分别采用基于两方程k-ω剪切应力输运(SST)湍流模型的延迟DES(DDES)、更改的DDES(MDDES)和改进的DDES(IDDES)方法,并引入可压缩修正,结合三阶MUSCL-Roe和五阶WENO-Roe两种空间离散格式,针对超声速底部的复杂流动现象,开展了数值模拟研究。计算结果表明本文方法能够捕捉到超声速底部流动中丰富的湍流结构,通过分析计算结果对超声速底部的流动机理有了进一步的认识,为下一步的超声速底部流动减阻改进和雷诺平均NavierStokes/大涡模拟(RANS/LES)方法在非定常高可压缩性流动中的应用提供了参考。通过对比分析不同空间离散格式的计算结果研究了数值耗散对计算的影响,五阶WENO-Roe格式的计算结果与实验结果吻合良好;对不同RANS/LES混合方法的计算结果进行了对比分析,结果表明IDDES方法在近壁区的表现优于DDES和MDDES方法。     

脱体涡模拟方法研究及其应用

基于现有的非定常雷诺平均N-S方程的求解方法,研究与发展了脱体涡模拟方法（DES）。在紊流附面层内用RANS方法,湍流模型采用S-A模型,在其他区域结合S-A模型运用Smagorinski的LES模型,并将S-A模型中的物面耗散项进行适当的改进,使得改进的模型既充当了RANS中的S-A紊流模型,又充当了Smagorinski大涡模拟模型。为验证所发展方法的有效性,将研究的DES方法应用到一些典型流场计算中,并与已有结果或实验结果进行比较。     

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

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

一种新型的动态RANS/LES混合方法

为解决RANS/LES混合方法过渡依赖网格疏密的问题,使用熵函数fs作为交界面函数构造一种新型的动态RANS/LES混合方法。同时采用脱体涡模拟方法 (DES97),对平板边界层流动、超声速凹腔-压缩拐角流动以及超声速底部流动进行数值模拟,并与实验数据进行对比。结果表明,无论对于附体流动区域还是分离区域,fs都能够动态地指示出边界层的范围。构造的动态RANS/LES混合方法能够显著缓解DES97由于网格因素导致的雷诺应力模化不足的现象,同时保证分离区域涡结构的求解精度。另外,熵函数fs是从能量耗散的角度进行构造,并不依赖于特定的湍流变量,因此更具有普适性。     

基于当地流动结构的RANS/LES混合模型

针对传统雷诺平均Navier-Stokes/大涡模拟(RANS/LES)混合模型固有的RANS和LES区域衔接难题,借鉴Vreman提出的亚格子(SGS)模型,提取其中描述流动状态的判断函数,利用不可压缩槽道湍流标定模型系数,构造了一个新的RANS/LES混合模型。评估算例包括不可压缩槽道湍流、稳态超声速平板边界层、NACA4412翼型绕流和亚临界雷诺数圆柱绕流。该模型作为壁面模化LES,解决了传统RANS/LES方法中速度型偏离对数律(LLM)的缺陷,且拓宽了原Vreman模型在粗糙网格上的适用性,具有较高精度。该模型同时解决了模化应力折损问题,且对大尺度分离流动的模拟精度相较传统的RANS/LES混合模型如脱体涡模拟(DES)有进一步提升。     

舰船空气流场数值模拟方法研究

分离涡数值模拟方法(DES)在物面附近求解雷诺平均Navier-Stokes方程,在其他区域采用大涡模拟方法,兼具前者计算量小的优点和后者能模拟大分离湍流流动、计算精度高的优势。论文分别采用RANS和DES方法对舰船空气流场进行数值模拟,通过所得流场结果与实测数据的对比分析,说明DES方法对大分离流动具有更强的模拟能力,更加适用于舰船空气流场预报。     

连续尺度变换方程在稳态及非稳态流动模拟中的应用

为了提高涡粘性假设的湍流模型对于非稳态流动的求解精度,同时兼顾其对于稳态流动的求解性能,将雷诺应力项与连续变换方程(CSSE)结合而形成新的应力项,使其根据流场尺度、网格尺度及Kolmogorov尺度来自动调节当地的应力雷诺应力模化水平,避免网格因素在流场模拟中产生不利影响,改正了混合RANS/LES方法的速度型偏离对数率问题;同时,该方法并未引入显式亚格子模型(SGS),因此回避了亚格子系数确定对于流场模拟精度产生的影响,改善了湍流模型对于流动不稳定性的辨识精度。在湍流平板算例中,CSSE方法计算的边界层速度型精度与雷诺平均方法(RANS)相当,而对于圆柱尾迹的模拟则证明了CSSE方法具有混合RANS/LES方法的优点,即能够准确模拟流动的不稳定性特征。     

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

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

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

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

黏性彻体力建模方法及其应用

研究了将Navier-Stokes方程中黏性项简化为黏性彻体力的建模方法,并将该模型加入已有的计算流体力学(CFD)程序中,实现了用三维欧拉方程加黏性彻体力模型快速模拟叶轮机三维非定常流动的方法.使用该方法对进气周向总压畸变条件下的某双级压气机非定常流场进行了模拟,并与雷诺平均的Navier-Stokes(RANS)程...     

超声速湍流边界层中横向声速喷流的混合LES/RANS模拟

通过对超声速湍流边界层中横向声速喷流的计算,对一种混合大涡/雷诺平均Navier-Stokes(LES/RANS)模拟方法进行了测试,该方法采用一个依赖于到壁面的距离及当地湍流参数的混合函数结合两方程k-ω SST(shear stress transport)湍流模型和混合尺度亚格子模型来封闭湍流项.计算结果表明:混合模拟方法能够捕捉到喷流/湍流边界层相互干扰的非定常大尺度结构,且对分离区长度、壁面静压峰值和膨胀区静压分布的计算精度要高于RANS(Reynolds-averaged Navier-Stokes)方法.      

不同淬熄结构下RQL燃烧室流场特性的数值研究

采用雷诺平均Navier-Stokes(RANS)方法对不同淬熄结构的富油/淬熄/贫油(RQL)燃烧室内定常掺混流动过程进行了数值模拟,分析了不同截面的流场结构,研究了动量通量比及淬熄孔排布方式的影响。结果表明:射流进入RQL燃烧室后发生偏转且偏转方向和旋流方向有关,回流区长度、穿透深度和高湍动能区随着动量通量比的增大而增加;当J≥80时,淬熄射流下游出现高湍动能区,逐渐与上游融为一体;当J≥120时,回流区的长度基本不再发生明显变化。此外,非中心对冲结构C3具有更大的高湍动能区,掺混更均匀。      

内管射流占优时轴对称共轴射流结构分析

以轴对称Reynolds-averaged Navier-Stokes(RANS)方程为控制方程模拟了内流占优条件下轴对称共轴喷嘴射流流场,并基于湍动能、速度和压强等流场参数的变化规律,将这类流场划分为交汇段、初始段和主流段等三个区段,给出三段流场的结构特征和流动机理,发现初始段内湍流动能大大低于单管射流,将导致共轴射流噪声低于单管射流;主流段内湍流动能则高于单管射流,将加强射流的掺混作用.将流量计算结果与试验结果进行了对比,证明使用的算法是有效的.      

Prediction of Boundary Layer Transition Based on Modeling of Laminar Fluctuations Using RANS Approach

This article presents a linear eddy-viscosity turbulence model for predicting bypass and natural transition in boundary layers by using Reynolds-averaged Navier-Stokes (RANS) equations. The model includes three transport equations, separately, to compute laminar kinetic energy, turbulent kinetic energy, and dissipation rate in a flow field. It needs neither correlations of intermittency factors nor knowledge of the transition onset. Two transition tests are carried out: flat plate boundary layer under zero ...     

基于CFD的带控制舵导弹的动导数计算

介绍了基于计算流体动力学(CFD)的阻尼导数预测方法以及其在复杂带控制舵导弹外形的俯仰/滚转阻尼导数预测中的应用.方法采用求解三维非定常ALE(arbitrary Lagrangian-Eulerian)形式RANS(Reynolds-averaged Navier-Stokes)方程加湍流模型的方法得到模型强迫运动的气动力时间历程,再通过绕模型非定常运动过程得到的气动力迟滞曲线计算得到俯仰及滚转阻尼导数.为了验证采用数值方法的正确性,针对基本带翼导弹外形,采用强迫模型按正弦形式振动,计算得到了超声速范围内马赫数从1.5到2.5的俯仰及滚转稳定性参数,通过与风洞试验及弹道数据结果的比较发现,两者吻合较好.