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

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

湍流大涡数值模拟进展

本文简要陈述湍流大涡数值模拟的原理、优点,着重讨论湍流大涡数值模拟方法的关键问题及其可能解决的途径,包括脉动的过滤、亚格子模型、近壁模型和标量湍流的大涡数值模拟中的特殊问题.文章强调大涡数值模拟中亚格子应力的本质是可解尺度湍流和不可解尺度湍流动量间的输运,并以作者最近提出的新型亚格子模型说明发展亚格子模型的正确途径.文章最后提出湍流大涡数值模拟近期需要迫切解决的问题和其他具有挑战性的方向.     

大涡模拟中亚格子模型的改进及其在槽道湍流中的应用

在用大涡模拟的方法计算具有强剪切的槽道湍流时, 常用的亚格子应力模型, 包括考虑壁面修正的模型, 都会出现平均流剖面偏高的现象.这表明在壁面附近亚格子应力模型还不能描述实际情况.针对这一问题, 修正了亚格子雷诺应力模型的壁面函数, 得到了较好的计算结果.用修正后的模型计算出的平均速度分布、均方根速度的分布以及雷诺应力的分布, 均与直接数值模拟(DNS)的结果吻合较好.      

旋转光滑直通道湍流流动二维热线实验

采用二维热线测量了旋转光滑直通道内不同流向位置上的平均速度和雷诺应力。实验结果表明:较高的当地旋转参数使得旋转对平均速度的影响区域扩大并且导致了前缘面附近湍流流动的再层流化。后缘面附近无量纲主流平均速度型严格按照旋转数顺序依次排列，且在对数律区呈现对数律分布。与此同时，所有无量纲雷诺应力分量在后缘面附近基本不受旋转影响。再层流化导致了前缘面附近无量纲主流平均速度型无法在对数律区维持对数律分布，且所有雷诺应力分量都随着旋转数和流向位置半径增大不断衰减；经过u-v象限分析，再层流化现象的直接原因被归于湍流脉动生成减少。      

桥梁断面静风荷载的格子Boltzmann方法数值计算

通过引入反映湍流涡粘性的湍流松弛,得到了模拟高雷诺数湍流的BGK方程.在速度相空间、物理空间和时间上对BGK方程进行离散得到了三维十九速离散速度模型;结合分区计算技术,设计了格子Boltzmann并行算法;根据亚格子Smagorinsky模型,提出了直接从粒子分布函数计算湍流松弛时间的方法.用开发的并行计算程序对分体双箱截面和闭口箱梁截面的静风荷载进行了数值识别,得到的静力三分力系数和流场压力分布与风洞试验结果及CFD宏观方法计算结果吻合,并从表面压力分布入手分析了两种桥梁截面的绕流特点.     

旋转状态气膜换热机理的数值研究

大涡模拟研究了旋转对单孔平板气膜冷却的影响,气膜出流的雷诺数为2600,吹风比为0.5,计算了静止和旋转数为0.02两种状态的流动和换热,从湍流结构演化的角度讨论了旋转状态气膜换热的机理.结果表明:①发卡涡是决定气膜换热的主要湍流结构,旋转状态发卡涡的形状和运动规律主要受哥氏力的影响;②气膜在哥氏力作用下向高半径方向偏...     

不同密度比下旋转方通道内流场和换热的数值模拟

在雷诺数为25000、旋转数为0.24、密度比为0.07~0.22的范围内,以数值计算的方法模拟了旋转方通道内三维流场及换热分布,与公开文献中的换热实验结果进行的对比表明,低雷诺数k-ω模型的计算结果与实验值吻合得相对较好;重点研究了旋转状态下冷气密度比对通道内流场和换热的影响,分析了哥氏力和浮升力在通道中的交互作用机理,结果表明,哥氏力引发的截面二次流是造成旋转通道前后缘换热差异的主要因素,浮升力加剧了通道内主流型的偏移,同时,其在前缘表面诱发的流体分离改变了局部湍流强度和换热分布。     

旋转U型通道流动与换热特性的数值研究

应用k-ω SST（shear stress transport）湍流模型,计算分析旋转U型通道在不同进口雷诺数（10000~60000）和高旋转数（0~2.013）范围内的流动与换热特性.结果表明:在静止和旋转状态下,进口雷诺数越大,努塞尔数越大.相比于同一工况下的静止状态,旋转显著增强了径向外流直通道的换热强度,径向内流直通道换热强度增大不明显.旋转数对U型通道换热的影响主要通过改变哥氏力和浮升力的大小.受哥氏力的影响,径向外流直通道后缘面换热增强,前缘面换热减弱.浮升力诱发了近壁面的流动分离,使得径向外流直通道前缘面不同位置处的换热强度随旋转数的增加而先减小后增大,计算得到的临界旋转数变化规律与实验测量结果保持一致,即无量纲距离参数与临界旋转数的乘积为定值.      

适用于旋转通道各向异性湍流模型

基于代数雷诺应力方程的简化模型,并结合充分发展旋转通道的直接数值模拟（DNS）湍流脉动数据，发展了一种适用于旋转通道的各向异性k ω模型。采用该模型对进口雷诺数为6000，旋转数为0~0.26的旋转直通道进行模拟，结果表明:将旋转修正系数乘以传统的湍流黏度发展的各向异性k ω模型，能够准确地描述旋转状态下前缘面和后缘面的换热情况，是一种有效的各向异性湍流模型修正方法；旋转修正系数是否合理的关键是对雷诺应力比值进行准确地描述；通道的换热效果与旋转数和流向沿程无量纲位置有关,前缘面的换热随旋转数和无量纲流向距离的增加而减小，而后缘面的换热随旋转数和无量纲流向距离的增加而增大。      

可压缩湍流边界层近壁区马蹄涡的演化

采用大涡模拟方法对可压缩槽道湍流进行了数值模拟.通过给定扰动波方法确定进口边界条件,很快得到了完全发展的湍流.利用动态亚格子模型、预处理技术以及高阶格式和方法,得到了精确的数值结果.计算得到的速度剖面和脉动分量与直接数值模拟的解进行了对比,验证了计算结果的可靠性.结果清楚地展示了湍流场中马蹄涡的生成和演化过程,包括对称涡腿、单侧涡腿的马蹄涡结构,以及流场中存在由亚谐波引起的拟序结构的交错现象.分析了湍流边界层近壁区可能出现的复杂的马蹄涡结构,及其涡结构形状也不严格对称的现象.     

旋转带肋U型通道内流动与换热的数值模拟

为了获得能够模拟高参数涡轮叶片内冷通道换热效果的模型，数值模拟了旋转状态下U型通道内的流场和温度场，比较了数值模拟与实验的结果。结果表明：所采用的计算程序和模型与实验结果吻合。旋转状态下，通道内各面换热的变化是和通道内流场的变化密切相关的；哥氏力在垂直于旋转半径截面上的不均匀分布引起流动较大变化，对通道内各面换热的影响比较大。哥氏力的作用较大幅度强化指向面换热，小幅强化两侧面换热，而弱化背向面换热。对于带肋通道，总体上阻力系数随着旋转速度的增加而升高。     

涡轮导向叶栅通道内紊流特性的实验研究

本文对一种高压涡轮导向叶栅中的紊流特性进行了实验研究。在雷诺数为100000的工况下, 采用二维热线对叶栅通道内垂直于叶高的平面内的平均速度, 脉动速度, 及雷诺切应力相关系数进行测量。并结合以前通过五孔针测量取得的流场进行了分析, 得出了该叶栅中紊流参数的变化规律。实验结果表明: 叶栅端壁区的紊流流动, 基本保持各向同性; 通道涡和马蹄涡内存在高紊流度和雷诺切应力的紊流核心。本实验结果为分析研究叶栅端壁换热特性奠定了基础。      

旋转带肋直通道湍流流动TR-PIV实验

采用高频粒子图像测速系统(TR-PIV)测量了旋转带肋通道内的主流平均速度、雷诺切应力、再附点等参数,并研究其沿程变化规律,通道高宽比为1,肋的阻塞比为01,雷诺数为10 000,旋转数从0变化到052,实验结果表明:静止时流动呈对称分布,但旋转后产生的哥氏力会极大地影响通道内的湍流流动,随着转速的增加,主流速度型偏向后缘面,前缘面涡系结构不断增大,再附点不断后移,而后缘面正好相反,并且这种趋势会沿程发展;前缘面附近的雷诺切应力变得越来越弱,而后缘面则越来越强,沿程雷诺切应力极值基本不变,但下游区域有所扩大。      

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

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

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

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

雷诺数对涡轮叶尖流场影响的数值研究

采用基于压力修正的三维计算流体力学程序，结合雷诺应力湍流模型和剪切应力传输湍流模型加壁面函数的方法，对某一轴流涡轮转子叶尖间隙流场进行了数值计算研究，详细计算了不同涡轮叶尖间隙高度和来流湍流度条件下雷诺数对涡轮转子间隙流场的影响，最后计算了转子效率。结果表明：当泄漏流流经叶尖时因为叶尖剪切力做功有块总压增大区；雷诺数带来的影响比湍流度和叶尖间隙高度带来的影响要大，湍流度的变化对流场影响不大；雷诺数对泄漏涡尺寸的影响不大，但低雷诺数会引起流动分离，带来损失，当雷诺数在文中的范围内变化时，效率下降近八个百分点。