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

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

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

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

基于RANS/LES混合方法的箱梁气动导纳函数研究

为了研究RANS/LES(Reynolds-averaged navier-stokes/large eddy simulation)混合方法在桥梁非定常气动力的应用以及大积分尺度湍流中桥梁断面的气动导纳函数,首先选取3种典型的RANS/LES混合方法,即SBES(stress-blended eddy simulation)方法、SAS(scale-adaptive simulation)方法、WMLES(wall-modeled LES)方法,分别对3种方法的特点和适用性进行简要的概述,然后在均匀流中将其用于箱梁断面和分离式双箱梁断面气动力的计算,最后在提出的改进湍流合成方法基础上,实现了与实际桥位处积分尺度类似的大积分尺度湍流模拟,在此流场中识别了箱梁断面的气动导纳函数。数值计算结果表明,SBES方法适用于弱流动不稳定和强流动不稳定流场气动力的模拟,SAS方法仅适用于显著流动分离的断面,而WMLES方法对两种断面整体表现都较差。对于流线型箱梁断面,大积分尺度湍流场中识别的升力气动导纳在低频处大于Sears函数和小尺度湍流场中的结果,并显著大于节段模型风洞试验结果。在高频处大积分尺度湍流场的结果与小尺度湍流场以及风洞试验结果吻合较好。研究结果表明,将传统小尺度格栅湍流场风洞试验识别的气动导纳用于桥梁抖振分析时,会导致低频处的抖振力显著偏小。

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

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

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

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

亚临界区圆柱绕流相干结构壁面模化混合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方法研究

混合RANS/LES方法是近年来复杂湍流模拟的重要方法。简要回顾了混合RANS/LES模拟方法的发展历程,着重总结和分析了分离涡混合模型、类Menter SST混合模型以及湍流能量谱一致混合方法三类混合模型的构造方法和发展历程,对这三种方法的应用以及优缺点进行了简要评述,最后指出了混合RANS/LES方法发展应该兼顾的问题,为下一步的混合RANS/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方法的亚声速空腔流动主要影响因素的数值研究

随着计算流体力学(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混合方法的分离流动模拟

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

不同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混合模型更为准确和高效,消耗的计算资源更少。

联合RANS/LES方法数值模拟方柱绕流

在粗细两套网格上,采用联合RANS/LES数值模拟Re=2.2×104的非定常方柱绕流经典算例.以标准的k-ε模型为基础,采用两种亚格子湍流模化方法,分别为PRNS(partially resolved numerical simulation)和DES(detached eddy simulation).PRNS模型与计算网格无显式关联,而DES与网格有一定关联.计算结果表明,联合RANS/LES方法即使在较粗的网格上也能得到较好的结果,DES方法预测到了二维模拟中的"双周期"现象,整体表现优于PRNS方法;亚尺度湍流模化中与计算网格的关联程度对结果有较大影响.

二维激波/湍流边界层干扰的混合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模拟

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

基于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方法具有很好的优越性和工程实用性。

基于SDES方法的翼型动态失速数值模拟研究

动态失速是航空航天领域中一种普遍存在的流动现象,因气动特性很难预测,对飞行器的性能和飞行安全有着严重影响。为了对机翼动态失速过程中的复杂分离流动进行准确模拟,应用屏蔽分离涡模拟(shielded detached eddy simulation,SDES)方法对NACA 0012翼型的动态失速进行了数值模拟,并与延迟分离涡模拟(delayed detached eddy simulation,DDES)方法进行对比研究。首先验证基于静态翼型的CFD(computational fluid dynamics)模型准确性,随后对动态失速获得的升阻力系数和实验数据进行比较评估,同时分析和量化了动态失速的特征现象,并对动态失速的几个典型阶段进行研究。计算表明:对于静态翼型的分离流动,SDES、DDES方法计算结果相差不大;对于具有强烈流动非定常性的动态失速问题,SDES方法体现了较高的预测准确度,能精细的描述流场及漩涡结构并准确地捕捉到动态失速下复杂流动现象,表现出良好的对于复杂分离流动精确模拟的能力,这对于新型RANS/LES混合方法的应用以及研究翼型绕流机理及相关问题有重要意义。

基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG／FV 混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度 DG 方法计算量、存储量小等优点，自提出以来已成功应用于一维、二维标量方程和 Euler／N-S 方程的求解。综述了 DG／FV 混合方法的研究进展，重点介绍了 DG／FV 混合方法的空间重构算法、针对 RANS 方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析，并给出了系列粘性流算例的计算结果，包括用于验证混合方法数值精度的库埃特流，以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明 DG／FV 混合方法达到了设计的精度阶，且相比同阶 DG 方法计算量减少约40％，而隐式方法能大幅提高定常流的收敛历程，较显式 Runge-Kutta 的收敛速度提高1～2个量级。

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

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

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

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