高超声速二维钝楔边界层参数分析

为了指导二元高超声速进气道型面设计,基于两类高超声速二维钝楔外形,采用高精度CFD方法计算了四种典型工况下二维钝楔边界层参数,分析了前缘半径、壁面温度、压力梯度对边界层参数的影响,得出了边界层位移厚度、边界层总压恢复及壁面压力梯度沿流向的变化规律。研究表明,前缘半径及壁面温度对边界层参数影响较大,在两者一定的情况下,逆压梯度使边界层位移厚度增速放缓,顺压梯度会使边界层位移厚度快速增长。     

一种基于湍流模式的转捩预测方法

通过对平板边界层进行计算分析,基于湍流模式来预测流动转捩是可行的,但预测出的转捩位置不准,比实际情况靠前.因此,本文从转捩流动的物理特征出发对湍流模式中的混合尺度作出了修正,提出了一种考虑间歇函数影响的k-ω SST两方程湍流模式.采用该模式来对平板湍流边界层和不同攻角下的NACA0012翼型绕流进行了转捩预测,计算结果表明,改进后的模式对转捩位置具有较好的预测能力.     

边界层稳定性的三维线性抛物化算法研究

通过把扰动波在流向上的变化分解成缓慢变化的形状函数部分和快速变化的波动部分,同时把扰动波在展向上的分布也分解成形状函数部分和波动部分,得到三维抛物化稳定性方程(3D PSE).然后在流向采用一阶向后差分,法向和展向采用4阶中心差分,发展出了迭代求解三维线性PSE的方法.通过求解二维扰动波和三维斜波在三维空间中的发展,发现二维波可以自动调制出展向扰动,且展向扰动的最大值点在边界层与外流的交接面附近.三维斜波的计算结果表明,网格密度应该与波矢方向上的波长相匹配,而不是与流向波长相匹配.     

近壁湍流边界层壁面局部扰动诱导形成相干结构的动力学模型

在壁面局部脉冲扰动作用下,采用直接数值模拟的方法,数值研究了湍流边界层近壁区两个相干结构形成的理论机制和非线性作用问题.然后,通过数值计算的手段,进一步探讨了近壁湍流边界层流中两个相干结构之间的非线性作用机理,分别确定它们在流向、展向方向上排列的不同结构类型分布的两个相干结构之间相互作用的效果有何区别以及它们之间间距的大小对相干结构的产生和发展有什么影响,从而可能合理地解释平面剪切湍流中低速条纹结构出现的原因.     

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

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

激波/湍流边界层干扰分离泡直接数值模拟

采用直接数值模拟(DNS)方法对来流马赫数为2.25、33.2°激波角的入射激波/平板湍流边界层干扰分离泡进行了数值研究。在验证了计算结果可靠性的基础上,通过分析比较3个不同展向站位分离泡的非定常运动特性、分离微团几何特征和相干结构等,定量考察了三维展向结构差异的影响规律。研究发现,分离泡存在复杂的三维结构,其流向长度明显大于法向高度和展向宽度,整体上沿展向呈现中间高两边低的扁平型单峰结构。分离泡面积脉动预乘功率谱结果表明,分离泡的非定常运动表征为大尺度低频膨胀/收缩过程,其展向三维结构对峰值频率的影响较小,且分离泡两侧略滞后于中间。采用经验模态分解(EMD)方法对分离泡低频膨胀/收缩过程进行了条件统计分析。统计结果表明,膨胀和收缩运动对分离微团几何特征没有实质影响,各展向站位分离微团高度/长度比值的概率峰值出现在0.1附近,同时分离微团面积和法向高度近似满足二次方分布。此外,流向速度脉动场的本征正交分解(POD)分析指出,分离泡的非定常运动与低阶模态密切相关,而高阶模态的贡献相对较小。采用前10个低阶模态可以准确重构出分离泡的低频膨胀/收缩过程。     

激波/边界层干扰及微型涡流发生器控制研究进展

激波/边界层干扰是一种发生在超声速/高超声速流动中的普遍现象。该现象将引起分离、流场结构振荡、局部高热通量和压力载荷。主要总结了近十年来激波/边界层干扰特性与微型涡流发生器及其组合体在流动控制中的最新进展。微型涡流发生器是目前研究最多、应用最广泛的控制方法,其流动机理和控制特性被大量挖掘。为了适应来流条件的变化、满足实际工况的需要,应开发定量评估和参数化设计方法。同时,应探索微型涡流发生器与其他控制方法的组合,实现更大程度、更广范围流场的控制。     

不同大气稳定度下台风\"莫拉克\"(2009)的近地风场特性

为了探究台风影响期间不同大气稳定度条件对近地层风场特性的影响,基于2009年台风\"莫拉克\"移动路径上的三座测风塔测得的高频三维风速时程,分析了台风影响期间中性、稳定、不稳定三种大气层结状况的分布频率,以及不同稳定度条件下近地层平均风场特性(平均风速和风向)和脉动风场特性(湍流强度、阵风因子、湍流积分尺度和脉动风速谱).研究结果表明,台风外围雨带区影响期间非中性层结发生概率约44％;在非中性层结影响下,平均风速出现频率最高的区间(10～12?m/s,14～16?m/s)不同于中性层结(12～14?m/s),且湍流强度最高频数分布较为集中(8％～10％);验证von-Karman谱描述不同大气稳定度下台风脉动风速谱的适用性,发现该经验谱对中性层结下的风速谱刻画较好,对非中性层结下的风速谱存在低频段高估、高频段低估谱能量的现象.进行工程结构抗风设计时,需注意非中性大气层结对风速、湍流强度、脉动风速谱等与结构风荷载确定密切相关的参数的影响.     

有/无尾迹作用下低压涡轮叶栅分离边界层转捩的大涡模拟

采用经过大量算例验证的可压缩大涡模拟求解器对雷诺数为60154、马赫数为0.402的低压涡轮叶栅T106D-EIZ进行了细致模拟,计算了定常来流和周期性尾迹来流两种工况.对计算结果的分析表明:定常来流工况下,叶片吸力面后部出现大尺寸的层流分离泡,分离剪切层的转捩过程受Kelvin-Helmholtz （K-H）不稳定性控制;尾迹来流工况下,由于来流尾迹的周期性扫掠,时均分离泡尺寸变小,叶栅总压损失降低.对相位平均和瞬态流场的分析表明,尾迹引起的逆射流使分离点后移,形成卷起涡结构,逆射流掠过卷起涡的过程中与其发生强烈的相互作用,产生大量气动损失,而后卷起涡破碎,流动转捩为湍流.      

三维坑湍流分离流动的粘流／无粘流干扰计算

本文对坑的三维分离流动做了低速粘流与无粘流的相互作用计算。对三维边界层反方法进行了分析和讨论。用数值试验的方法验证了在H和α作为已知量的情况下,三维边界层反方法的积分方程是双曲型的,并提出了一种近似数值特征线法进行求解。无粘流采用低速位流面元法。计算表明所用方法可计算出三维效应很强(即横向变化很大)的三维边界层分离流动。