后向台阶层流分离剪切层特性研究

应用激光测速仪对生向台阶层流边界层分离产生的剪切层在再附前的发展进行了测量，得到了时均速度、湍流强度及剪切层厚度等的分布特性，并研究了来流湍流度对台阶后流场湍强度分布及剪切层发展的影响，实验中发现台阶后有一个流速降低的区域，对此从涡动力学的角度作了解释。     

来流湍流度对边界层发展的影响

应用ＬＤＶ（Ｌａｓｅｒ Ｄｏｐｐｌｅｒ Ｖｅｌｏｃｉｍｅｔｅｒ）测速技术及氢气泡流动显示技术对水流中平板边界层的发展进行了详细的观测，研究了来流湍流度对边界层转捩、湍流边界层特性等的影响，得到了层流边界层、转捩区、湍流边界层区的分区界限，在本实验条件下，来流湍流度对边界层转捩有非常明显的作用，但对湍流边界层的特性影响不大。     

基于Transition SST模型的高雷诺数圆柱绕流数值研究

为了研究高雷诺数下圆柱绕流边界层的转捩现象和圆柱尾迹近壁区的流动特征,首先通过在典型雷诺数下采用Transition SST四方程转捩模型模拟圆柱绕流得到的结果与实验结果及采用SST k-ω两方程湍流模型模拟结果的对比分析,验证了Transition SST模型在模拟高雷诺数下圆柱绕流的优越性,并较为准确地预测出了圆柱绕流边界层的转捩现象及尾迹近壁区的流动特征。然后分别对亚临界区、临界区、超临界区和过临界区的圆柱绕流问题进行了数值模拟,分析了不同雷诺数下圆柱绕流的流场结构及圆柱表面压力系数、摩擦力系数的变化规律,研究了圆柱绕流近壁区的流动特征、边界层转捩的流动机理、转捩位置及其随雷诺数的变化规律。结果表明,亚临界区二维圆柱绕流边界层发生层流分离,无分离泡和转捩现象;临界区和超临界区二维圆柱绕流边界层先产生了分离泡现象,之后流动发生了转捩并在转捩后发生湍流分离;过临界区二维圆柱绕流边界层流动在转捩之后发生湍流分离,无分离泡现象;在临界区、超临界区和过临界区,二维圆柱绕流边界层转捩位置随雷诺数增大向前驻点移动。     

用GAO-YONG湍流模型计算激波/边界层干扰

采用GAO-YONG可压缩湍流方程组,模拟了平板激波/湍流边界层干扰现象.运用SIMPLE算法求解方程组,并分别采用三阶精度的QUICK格式和中心格式离散对流项和扩散项.计算结果较好预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征: 分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波.对流场的时均参数与实验值进行了比较,计算得到的壁面压力分布、摩阻系数分布和速度型与实验值比较吻合很好.结果表明GAO-YONG可压缩湍流方程组能够高精度模拟平板激波/湍流边界层干扰流动.      

对流换热边界层内熵产率分布研究

本文对平板层流对流换热边界层内的熵产率分布做了一些理论分析。讨论了气流马赫数M_∞及无因次温差τ对熵产率分布的影响。在此基础上,对层流边界层内最大熵产率的位置进行了定性的分析。     

用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

轴对称钝体前缘分离剪切层的非线性行为

实验研究了锐前缘轴对称钝体前缘分离剪切层发展演化的动力学过程,以及分离剪切层与外加声激励之间的非线性作用.研究结果表明,自然状态下,存在着亚谐频现象,分离剪切层的演化发展是通过涡合并方式实现的.从分离到再附,层流分离剪切层一般经过2次涡合并过程,湍流分离剪切层一般经过5次涡合并过程.混沌动力学分析表明,层流分离剪切层是经过Hopf分岔和周期倍分岔途径进入混沌状态的.此外,讨论了分离剪切层对外加声激励的响应.     

G3风洞跨音速试验段流场M数和边界层特性测量

测定一个下吹、暂冲式三音速风洞的跨音速试验段流场M数分布和侧壁边界层特性。该试验段上下壁面开槽,截面尺寸为535×375mm。测定的试验M数范围为0.2～1.15。M=1.0时每米Re数为10×10~6。取国内规定的流场不均匀度ε_M作标准,所测各点未出现|ΔM_(max)|≥2ε_M情况。侧壁边界层的平均速度、湍流度和湍流应力分布的规律性良好,表明侧壁面可用作边界层实验。所测得的湍流度分布结果可作为今后制订跨音速风洞湍流度指标参考。     

粗糙度对涡轮叶片流动转捩及传热特性的影响

为研究表面粗糙度对涡轮叶片流动转捩以及传热特性的影响,在自行开发的CFD程序平台上提出了对γ-Re_θ转捩模型的粗糙度修正方法,并参考平板绕流和涡轮叶栅的实验数据对该方法进行验证。考虑粗糙度效应的影响,对Mark Ⅱ涡轮导叶5411工况进行数值模拟,得到如下结论:表面粗糙度对层流边界层换热系数影响不大,而对湍流边界层则有较大影响,进而显著改变壁面温度分布;与光滑壁面相比,5μm的等效沙粒粗糙度使吸力面湍流区域壁面温度升高约5.7K,100μm粗糙度使壁面温度升高28.4 K,增幅达5%左右;当壁面粗糙度较低时,激波干涉对吸力面边界层的转捩起主导作用,而当粗糙度大于某临界值时,其作用会使转捩位置突然变化,本算例中该临界值近似为150μm。      

超声速复杂流动中湍流模型的性能评估

对8°,16°,20°,24°超声速二维压缩拐角进行了数值模拟,系统评估了SA(Spalart-Allmaras), k-ω, SST(Shear Stress Transport)3种工程常用湍流模型对激波/边界层干扰复杂流动的模拟适用能力.对于"弱"干扰、小分离,工程常用湍流模型能够较准确地预测其压力、摩擦力、热流分布,而对于"强"干扰、大分离则会造成较大的偏差.同时,随着分离区的增大,各湍流模型的适用性能差别更加明显,数值方法也有一定的影响.      

用GAO-YONG湍流模式对分离流动的数值模拟

采用GAO-YONG湍流模式对二维管道内三角形钝体绕流问题以及槽道内后台阶分离流动进行了数值模拟.求解程序基于开源数值计算平台OpenFOAM (Open Field Operation and Manipulation),数值模拟结果很好地预测了钝体绕流问题中漩涡脱落的尾流的流动趋势以及后台阶流动中分离再附的流动结构,同时分析了速度分布以及摩擦系数等参数并与实验值进行了对比,结果符合很好.这表明GAO-YONG湍流模式对大分离流动有较好的预测能力,对工程实践具有指导性作用.     

三维斜掠后台阶流动数值模拟

介绍了开源计算流体力学软件OpenFOAM(Open Field Operation and Manip-ulation)及面向对象编程技术,并利用其对三维斜掠后台阶流动进行了数值模拟.通过对无斜掠角度情况下用OpenFOAM计算得到的流场速度型和展向涡强度场分布与实验数据的对比研究,证明其具有良好的计算精度和可靠性.通过对不同斜掠角度台阶下游流线与压力分布的详细研究发现随着斜掠角度从0°增加到60,°回流区涡团结构和再附距离均发生了显著变化.其总体趋势是再附距离在30°后会急剧缩小,旋涡强度也会减弱.     

沟槽面湍流减阻研究综述

对近20年来沟槽面湍流边界层特性、湍流拟序结构、湍流减阻及其机理的研究进展进行了综述.内容涉及沟槽面平板、旋成体、机翼等在压、跨、超音速流动情况下的实验研究成果;压力梯度、攻角、侧滑角等的影响;湍流猝发特性、紊动特性、近壁区带条结构的特征及减阻机理等方面的工作.为更有效地减少表面摩阻,必须深入开展对沟槽面湍流边界层特性、湍流拟序结构及湍流减阻机理等方面的研究.      

基于SA湍流模型的DES方法应用

使用基于SA(Spalart-Allmaras)湍流模型的分离涡模拟(DES,Detached-Eddy Simulation)方法,对圆柱绕流和超声速圆柱底部流动等高雷诺数大分离流动进行了数值模拟.无粘通量的计算采用二阶迎风格式,时间推进使用二阶隐式双时间步法.计算表明相对于雷诺平均N-S方程(RANS,Reynolds-Averaged Navier-Stokes)模拟,时空均为二阶精度的DES模拟能够分辨更为精细的分离区三维旋涡结构,对分离区主要湍流结构的捕捉和平均压力分布的预测取得了满意的结果.     

边界层转捩过程三维扰动的发展和演化

引入展向调制的二维人工扰动激发平板边界层转捩,利用水槽氢气泡流动显示技术观察到转捩边界层中典型三维扰动的产生和发展,包括Λ结构、发卡涡以及发卡涡演化生成环状涡的过程.实验结果显示出环状涡和尖峰结构的关系,直接证实了尖峰结构的产生机理.在流动紊乱化初始阶段发现常出现在湍流边界层中的暗斑结构,结果表明与暗斑相伴随的局部高剪切层不稳定可能是流动紊乱化的重要原因.     

沟槽面湍流边界层湍流度分布

应用激光测速(LDV)技术对沟槽面湍流边界层湍流度分布进行了精细的测量,实验结果表明,在湍流边界层区域沟槽面湍流度最大值比光滑面要小,但出现的位置比光滑面要远.此外,紧靠沟槽面,即y++<4内,沟槽面有较大的湍流度,说明沟槽面有使近壁区流动紊动增加的趋势.对于y/δ≥0.06的区域,湍流度可用u′\u=2.3 exp (-0.9y/δ)描述,而在y/δ<0.06内,则有u′/u =36.4y/δ.