亚临界区圆柱绕流相干结构壁面模化混合RANS/LES模型

本文发展了一种具有壁面模化大涡模拟能力的雷诺平均纳维-斯托克斯(RANS)和大涡模拟(LES)方法的混合模型(简称WM-HRL模型),致力于对亚临界区雷诺数钝体绕流相干结构这类复杂流动现象进行高置信度的CFD解析模拟研究.该方法通过一个仅与当地网格空间分布尺寸有关的湍动能解析度指标参数rk即可实现从RANS到LES的无缝快速转换,并且RANS/LES混合转换区的边界位置及其各个分区(包括RANS区、LES区及RANS/LES混合转换区)对湍动能的解析能力均可通过两个指标参数nrk1-Q和nrk2-Q准则进行预先设定.通过对雷诺数Re=3900下圆柱绕流场的系列数值模拟研究,获得了能够高置信度解析并捕捉其绕流场中三维时空瞬态发展相干结构特性的湍动能解析度指标参数nrk1-Q和nrk2-Q准则的组合条件.研究表明,该WM-HRL模型不仅能够准确获取圆柱绕流场中剪切层小尺度K-H不稳定性结构的精细谱结构,而且在同一套网格系统下通过变化湍动能解析度指标参数nrk2-Q和nrk1-Q准则的组合条件,还可以精细解析圆柱绕流场中两类不同回流区的长度结构特征,及其对应的圆柱尾部近壁面处V和U形两个平均流向速度剖面的分支结构特性.

采用RANS/LES混合方法研究分离流动

采用结合湍流模式和大涡模拟(LES)的混合方法数值模拟Aerospatial A-profile翼型以及细长椭球体大迎角分离流动。该类混合方法比较适于模拟分离流动:近壁区附体流动采用湍流模式,而远离物面分离区域内的大尺度运动则用LES方法。通过计算和实验结果的对比,混合方法切实可行并具有很大的发展潜力。

不同RANS/LES混合模型的汽车气动噪声分析

采用Fluent UDF构造Realizable k-ε/LES、SST k-ε/LES和Spalart-Allmaras/LES 3种RANS/LES混合模型,并对汽车的气动噪声进行数值模拟。在近壁面采用RANS湍流模型求解,而远离物面的分离区域内的大尺度运动则用LES求解,将3种混合模型的计算结果与全LES湍流模型和风洞试验结果进行了比较。结果表明,3种混合模型对气动噪声的计算精度高于LES湍流模型;Realizable k-ε/LES和SST k-ε/LES混合模型数值模拟结果与试验值的相对误差均小于5%,而其中Realizable k-ε/LES混合模型更为准确和高效,消耗的计算资源更少。

基于混合LES/RANS方法的引水隧洞施工通风数值模拟

为了获取准确的水电工程引水隧洞施工通风瞬态流场信息并且提高计算效率,基于Realizable k-ε两方程湍流模型建立引水隧洞施工通风三维非稳态混合LES/RANS模型,并结合某引水隧洞施工通风过程进行数值模拟,借助现场实测风速数据对数值模拟结果进行验证,同时将混合LES/RANS模型的模拟结果与RANS模型和LES模型的模拟结果进行对比。结果表明:混合LES/RANS数值模拟结果与现场实测数据较为吻合,其平均相对误差仅为4. 9%。混合LES/RANS方法仅需1 163 996个计算网格即可获得准确的瞬态风流场扩散规律,而LES方法所需网格数量近似为前者的2倍。混合LES/RANS方法可以为水电工程引水隧洞的施工通风数值模拟研究提供重要的理论与技术支撑。

复杂湍流流动的混合RANS/LES方法研究

混合RANS/LES方法是近年来复杂湍流模拟的重要方法。简要回顾了混合RANS/LES模拟方法的发展历程,着重总结和分析了分离涡混合模型、类Menter SST混合模型以及湍流能量谱一致混合方法三类混合模型的构造方法和发展历程,对这三种方法的应用以及优缺点进行了简要评述,最后指出了混合RANS/LES方法发展应该兼顾的问题,为下一步的混合RANS/LES模拟提供了参考。

用RANS/LES混合方法研究超声速底部流动

在剪切应力输运和弱非线性k-ω模式中引入压缩性修正,并在此基础上构造RANS/LES混合方法,对包含丰富流动结构和复杂流动机理的导弹类超声速底部流场进行数值分析,发现压缩性修正可更好地描述流场特征,适于具有较强压缩性效应流动的数值分析;混合方法捕捉到丰富的流动现象及非定常特性.

嵌入式RANS/LES混合方法在槽道流动数值模拟中的应用

与传统整体式RANS/LES混合方法相比,嵌入式RANS/LES混合方法不再依靠改进后的湍流模型本身或网格疏密自动划分流场中的RANS区域和LES区域,而是通过对流动现象的预估,人为将LES区域嵌入到RANS区域,从而使RANS区过渡到LES区的交界面尽可能与流向垂直,方便合理附加额外的湍流脉动信息。这样做不但确保了LES在某些复杂流动区域彻底支配地位,更重要的是能够避免RANS区对LES区解析湍流发展的污染和抑制。将槽道流场划分为上游RANS区和下游LES区,在经典的SST-DES模型基础上,采用合成涡方法生成交界面的湍流信息,取得了良好的模拟结果。

二维激波/湍流边界层干扰的混合LES/RANS模拟

为了降低高雷诺数条件下大涡模拟方法的计算量,将两方程k-ωSST湍流模型与Yoshizawa一方程亚格子模型通过一个衔接函数相结合,构造一种混合大涡/雷诺平均NS方程模拟方法(混合LES/RANS)。使用这种方法及AUSM+-up格式对20°压缩斜坡的马赫2.85流动进行模拟,并考察了固定入口和在入口添加白噪声两种湍流入口边界条件对于结果的影响。模拟结果再现了边界层的分离、再附以及分离激波等现象,计算得到的分离区要显著大于试验结果,对于这种混合模拟方法的缺点进行了分析,并提出了可能的改进方法。

多段翼型大迎角分离流动的DELAYED RANS/LES混合算法

多段翼型的大迎角绕流发生大范围附面层分离,具有明显的三维与非定常流动特性。RANS/LES混合算法继承了LES对流动分离区大尺度漩涡准确模拟的优点,避免了纯LES算法需求网格量巨大与亚格子模型壁面函数不成熟等问题,对分离流动的模拟效果优于RANS算法。以S-A湍流模型与Smagorinsky亚格子模型为基础,借鉴DDES的附面层延迟控制思想,构造了可用于对接网格、重叠网格的DELAYED RANS/LES混合算法。研究了GA(W)-1多段翼型的大迎角分离流动及其气动特性。

基于混合RANS/LES方法的亚声速空腔流动主要影响因素的数值研究

随着计算流体力学(Computational fluid dynamics,CFD)技术的发展,混合雷诺平均N⁃S方程(Reynolds averaged Navier⁃Stokes,RANS)/大涡模拟(Large eddy simulation,LES)类方法被越来越多地应用于航空领域的大分离流动及气动噪声的仿真分析中,并取得了很好的效果。本文选取M219空腔流动风洞试验模型,在Ma=0.85的情况下,进行了混合RANS/LES方法的验证,并进行了网格敏感性分析;然后针对影响空腔流动的几个关键因素进行了仿真研究,并给出了主要的研究结果;最后,针对L/D=7的大长深比空腔,研究了增加隔板形成前后两个L/D=3.5的串列空腔的流动特点。

后向台阶超声速流动的动力混合RANS/LES模拟

针对可压缩湍流流动,构建了1种动力混合RANS/LES模式。运用该模式模拟了超声速后向台阶流动,结果再现了二维后向台阶剪切层的拟序结构,计算得到了统计平均的压力分布,回流区长度与实验结果符合较好。与定系数的一般混合RNAS/LES模拟结果的对比表明,动力混合RANS/LES模式在对流场拟序结构的捕捉、流动时均压力分布方面表现较优。

基于RANS/LES混合方法的近岸水域波浪增阻影响的数值模拟

分析船舶在近岸水域航行时所受波浪增阻大小的影响因素,以寻求降低波浪增阻的方法。基于近岸水域中船舶阻力理论,采用雷诺平均(RANS)和大涡模拟(LES)相结合的方法进行数值求解,并运用两相流方法和造波消波技术建立数值水池模型。用Fluent求解器对船舶在不同工况下的波浪增阻进行数值模拟。试验结果表明:用该方法对近岸水域航行的船舶进行波浪增阻的数值模拟可有效提高计算机运算效率,在粗糙网格层面上得到的计算结果精度良好;船舶在近岸水域中航行时的波浪增阻随波高的增加而增加,并且与船体湿润面积呈正相关;船舶在近岸水域以中低速航行时所受波浪增阻随弗劳德数的增加而增加。该研究为近岸水域通航船舶减小航行阻力提供了可行的方法与思路。

超声速燃烧凹腔质量交换特性的混合RANS/LES模拟

对超声速冷流及燃烧流条件下超燃冲压发动机火焰稳定凹腔的质量交换特性进行了研究.采用混合RANS/LES(Reynolds-averaged Navier-Stokes/large-eddy simulation)方法对非定常流场进行数值模拟,系统研究了凹腔结构参数和横向喷流燃烧放热对凹腔质量交换特性的影响.首次采用"原子追踪法"对化学反应流条件下的凹腔质量交换过程进行了刻画.结果表明:在冷流条件和燃烧流条件下,凹腔驻留时间均随长深比的增大而增加;凹腔后缘倾角对驻留时间的影响在冷流条件下比较显著,而在燃烧流条件下明显减弱;相同结构的凹腔,在燃烧流条件下的驻留时间相对于冷流条件明显减小.

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

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

基于SA模型的改进延迟分离涡模拟方法研究

流动分离问题的精细化模拟需求日益增加,传统的URANS方法已经无法满足精细化模拟要求,以基于SA的RANS/LES混合方法为基础,对DES方法中存在的网格诱导分离(Grid Induced Separation,GIS)问题及DDES方法存在的对数层不匹配(Log Level Mismatch,LLM)问题进行了改进,进一步发展了IDDES方法。通过对单圆柱绕流、亚音速三角翼及CRM-HL高升力标模算例的仿真验证该方法的有效性。数值结果分析表明:与URANS方法相比,DES类混合方法模拟精度更高,能更好地反映出流场丰富的涡系结构;与DES和DDES方法相比,发展的IDDES方法由于克服了DES和DDES方法的缺陷,数值模拟结果更接近实验值,可见IDDES方法具有很好的优越性和工程实用性。

“回收/调节”方法在混合LES/RANS模拟方法中的应用

采用一种混合大涡模拟/雷诺平均Navier-Stokes(LES/RANS)模拟方法结合三阶加权基本无振荡(WENO)格式对马赫数为2.88的压缩斜坡流动进行了数值模拟,并采用"回收/调节"方法生成入口湍流脉动边界条件.当采用固定入口边界条件时,湍流边界层会缺乏合理的湍流能量,使之抵抗逆压梯度的能力减弱,会严重高估分离区的大小;而当采用"回收/调节"方法时,可以在10倍来流边界层厚度的距离内激励起湍流大尺度结构,计算得到的流场能够再现激波/湍流边界层干扰流场的非定常特性,时间及展向平均壁面压力和试验值吻合较好.

混合型三阶格式及关于两种DES算法的比较计算研究

为开展RANS/LES混合模拟,在传统三阶迎风偏置格式的基础上,提出一类混合型的三阶计算格式,格式在中高波数范围具有可调的耗散水平。在此基础上,文中通过计算比较研究两种RANS/LES混合算法:一种是基于Spalart-Allmaras一方程模型的DES模型(DES-SA),另一种是基于混合长度模型的DES模型(DES-ML),使用发展的三阶格式对圆柱绕流进行模拟,并将得到的结果与文献比较,对计算格式、DES-ML算法进行了初步探讨。

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

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

低速大雷诺数混合飞艇气动性能分析

针对混合飞艇体积巨大同时气动外形复杂使得现有条件的风洞试验很难精确测量其气动性能的问题,开展了适用于混合飞艇气动性能分析的计算流体力学(CFD)的数值分析方法研究。考虑混合飞艇低速大雷诺数的特点,将变分多尺度方法 (VMS)与动态Smagorinsky大涡模拟(LES)模型相结合,提出了组合的VMSLES湍流模型。将基于RANS方法和LES方法的其他三种湍流模型相对比,利用雷诺数相近、实验数据丰富的6:1长椭球飞艇对不同的湍流模型进行了对比验证。结果显示LES方法预测结果与实验结果吻合较好,优于RANS方法,并能显示更多流动细节,而组合的VMS-LES模型能够更精确地捕获实验研究中观察到的二次涡。利用组合的VMS-LES模型对有翼HAV与多囊瓣HAV进行了气动性能分析,并研究了不同部件对飞艇气动特性的影响。结果表明,由于尾翼表面产生的一次涡与二次涡相互作用,尾翼在增加气动升力的同时也增加了阻力。

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

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