旋转槽道湍流的格子Boltzmann方法模拟

基于多松弛格子Boltzmann方法的大涡模拟对雷诺数为194,旋转数从0~5.0的旋转槽道湍流进行数值模拟.采用动态亚格子应力模型模化滤波后的不封闭项,修正二阶矩作用力模型计算压力梯度、哥氏力.对平均速度、均方根脉动速度、雷诺应力以及湍流结构进行了计算.结果显示哥氏力使流场平均速度呈现不对称性:在压力面,随着旋转数的增加,湍流度增强;而在吸力面湍流脉动减弱,具有层流化的趋势.将格子Boltzmann模型与直接数值模拟求解进行对比,结果验证了格子Boltzmann方法在旋转湍流模拟中的可行性.      

一种剪切修正的涡黏性亚格子模型及其在槽道湍流中的应用

提出了一种利用壁面附近强剪切对涡黏系数进行修正的亚格子模型,称为剪切修正的涡黏模型.该模型能使涡黏系数在壁面附近自动满足壁面修正,在脉动为零的地方,涡黏系数自动为零,保证了壁面附近涡黏系数趋于零的特性.用该模型对槽道湍流壁面无吹吸、有吹吸,较低和较高雷诺数的情况进行了数值模拟,与直接数值模拟的结果进行了对比,结果吻合较好.      

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

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

湍流模型对三维翼梢涡流场数值模拟的影响

 梢涡流场数值模拟对减少梢涡的不良影响意义深远。应用雷诺平均Navier-Stokes(RANS)方法建立了三维翼梢涡流场数值模型,选用的湍流模型包括两种经典两方程模型和一种代数雷诺应力模型。计算内容包括翼梢涡流场速度分布、压力分布及涡核位置,为了准确模拟梢涡涡核内系统旋转和流线曲率影响,将旋转和曲率修正方法应用于选用的湍流模型中。计算时对翼梢涡涡核处网格进行了加密处理。将计算结果与实验结果对比后发现,修正后的湍流模型较标准湍流模型明显提高了预报精度,与实验结果吻合较好,具有工程应用价值。     

低雷诺数翼型层流分离现象大涡模拟方法

为准确预测翼型低雷诺数条件下出现的层流分离流动现象,发展了基于结构化拼接网格的大涡模拟方法。控制方程为Favre滤波Naver-Stokes方程,并选用多种亚格子模型。空间离散采用AUSM格式以及高阶WENO格式,时间推进采用显式方法和隐式方法。以SD7003翼型在雷诺数60 000及4°迎角下的层流分离流动为研究对象,对比分析了数值格式、亚格子模型、网格尺度对流场预测结果的影响。数值结果表明:划分的计算网格能够有效解析小尺度流动结构,基于隐式亚格子模型、采用AUSM格式和双时间步推进的大涡模拟方法能够准确预测流动分离、转捩、再附等复杂流动现象,计算得到的平均压强系数与雷诺应力分布与文献数据吻合较好。与转捩模式计算结果对比进一步表明:发展的大涡模拟方法能够准确预测翼型低雷诺数层流分离流动非定常演化过程,为下一步的研究工作建立了有效的数值模拟手段。     

三维槽道湍流的大涡模拟

利用大涡模拟方法模拟了雷诺数为13000的充分发展槽道湍流流动。为了比较网格和亚格子模型对计算结果的影响,采用了三种不同的网格数目,以及Smagorinsky模型和动力模型等两种亚格子模型,并与直接数值模拟结果进行了比较。结果表明,大涡模拟可以比较准确地给出速度分布、脉动速度均方根分布等湍流统计参数。在此基础上,对近壁区小尺度湍流结构进行了研究。     

基于格子Boltzmann方法的方通道湍流的大涡模拟

基于格子Boltzmann方程的大涡模拟方法,对以摩擦速度、方通道水利直径为特征尺度,雷诺数为300的直方通道内湍流流动进行数值计算.利用多松弛时间格子Boltzmann方法来模拟流场的流动,切应力改善亚格子应力模型来模化滤波后的非封闭项.将模化后的亚格子应力与格子Boltzmann方法中的松弛时间相关联,使得松弛时间当地化,从而能够准确地模拟湍流.对湍流的平均流向速度、平均二次流速度以表征湍流强度的均方根速度以及不同截面流向瞬时涡做了计算和评估.计算结果与直接数值模拟、实验数据相吻合,证明了格子Boltzmann方法在计算通道湍流中的精度.      

室火轰燃的大涡模拟

采用一套修正的适用于低马赫数浮力流的N-S方程组来描述火灾场,运用Smagorinsky亚格子模型的大涡模拟方法对室内池火的轰燃发展过程进行了数值模拟,把模拟与实验结果进行了对比和分析。结果表明,该方法能较好地预测轰燃。同时模拟结果显示火焰区向房间后部偏转稍大,近壁处温度预测的误差大一些,随着亚格子模型的不断完善和计算能力的提高,用该技术来预测轰燃现象会更加准确。     

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

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

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

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

高超声速钝楔边界层转捩大涡模拟

 以五阶迎风和八阶对称格式混合差分格式求解三维可压缩滤波Navier-Stokes方程,对来流Mach数为6.0、半锥角5°的高超声速空间发展钝楔边界层转捩至完全湍流进行了大涡模拟。时间推进采用紧致存储三阶Runge-Kutta方法,亚格子尺度模型为Driest因子修正的Smagorinsky涡黏性模型。通过定常流场入口边界附近吹／吸引入不稳定扰动斜波的方法数值模拟得到了层流失稳转捩直至完全湍流的空间发展全过程。对扰动的线性、非线性增长以及湍流斑的形成和发展进行了分析,给出了转捩及完全湍流下的速度相关量统计并与实验、DNS结果进行了对比分析,计算结果与理论及实验吻合。     

超声速小攻角下的圆锥转捩攻角效应

采用高阶差分格式求解三维可压缩滤波Navier-Stokes方程,对超声速零攻角和小攻角圆锥边界层转捩进行了空间大涡模拟研究.在三维Tollmien-Schlichting(T-S)波扰动作用下,计算得到了流场非定常转捩过程.计算结果清晰地展现了非定常大尺度结构演化特征和攻角效应引起的非对称转捩现象,流场发卡涡、横流涡等转捩主控大尺度结构与实验一致.     

二维沟槽粗糙槽道内的湍流特性研究

采用激光多普勒测速技术对光滑和粗糙槽道湍流特性进行了实验研究。粗糙元为二维横向V型沟槽,沟槽深度为0.8mm,沟槽间距为6.4mm,对应的槽道半高度与沟槽深度比为12.5。基于中线时均速度和槽道半高度的流动雷诺数范围为2740～17400。实验测量了包括时均速度、湍流强度、雷诺切应力和速度脉动偏斜因子和平坦因子在内的湍流统计量,结果表明沟槽型粗糙度对湍流的影响不仅局限于边界层内区,而是延伸到整个边界层范围。粗糙壁面上的粗糙度函数随雷诺数的增大而增大,时均亏损速度也较光滑壁面高。沟槽抑制了内区的流向湍流强度,同时增大了外区的湍流强度。粗糙壁面上的雷诺切应力高于光滑壁面,且与湍流强度一样表现出对雷诺数的依赖性。尽管沟槽型粗糙度对流向平坦因子影响不大,但对流向偏斜因子有显著影响。     

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

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

NACA4412翼型低速绕流数值计算中湍流模型对比

使用Spalart-Allmaras(S-A)、SSTk-ω、Gao-Yong 3种湍流模型对NACA4412翼型低速绕流进行了数值计算,研究了尾迹流动松弛效应。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。计算中对翼型尾缘流松弛效应进行了分析,比较了翼型表面压力系数、速度剖面、雷诺应力等的分布,3种湍流模型总体上能够较好地模拟NACA4412翼型低速绕流。SSTk-ω模型对流动细节及升力系数计算最好,Gao-Yong模型对翼型平均速度剖面及雷诺剪切应力分布计算最准确。     

平板边界层壁面局部脉冲激发大涡结构起因的数值模拟

构建了边界层壁面局部脉冲扰动诱导大涡结构的初始理论模型,采用直接数值模拟的方法,研究了边界层壁面局部脉冲激发大涡结构起因的问题.计算结果表明:局部脉冲激发大涡结构成因的许多特点,如幅值演化、雷诺切应力功率、高低速条纹结构、流向涡的产生、上喷和下扫过程、平均速度剖面存在的拐点和严重扭曲以及类似于发卡涡的特征等现象与实验结果和某些数值结果非常相似;另外,壁面局部脉冲的强度、影响区域、加载时间、分布方式对边界层流中大涡结构起因的特性起着决定性的作用,壁面脉冲也是边界层流中激发大涡结构的重要途径之一.