非对称交叉激波和湍流边界层相互作用的数值研究

针对7°×11°双尖鳍外形的非对称交叉激波与湍流边界层相互作用,采用Navier-Stokes方程和5种湍流模型进行了计算.主要考察对壁面压强、热传导、绝热壁面温度和壁面摩擦力线分布的计算精度.计算结果表明:两道斜激波相交后的区域的壁面压强和热传导都比较高;计算的压强和壁面摩擦力线与试验吻合很好,绝热壁面温度次之,热传导最差,峰值高达试验的3倍左右.湍流模型对壁面压强和壁面摩擦力线影响很小,对绝热壁面温度和热传导影响很大.在5种湍流模型中,TNT(turbulent/non-turbulent)和SST(shear-stress transport)模型表现较好.      

应用GAO-YONG湍流模式数值模拟三维激波/湍流边界层干扰

应用GAO-YONG可压缩湍流模式数值模拟了三维激波/湍流边界层干扰算例之一——单鳍流动.攻角20°,来流马赫数2.93,雷诺数9.8×105.对流项和扩散项分别采用Roe格式和二阶中心差分格式计算.Runge-Kutta显示时间推进方法求解了半离散的控制方程.包括壁面压力分布,边界层内流动偏移角等在内的计算值与试验数据进行了比较.准确地预测出了三维激波/湍流边界层干扰流场的主要流动特性——λ波结构,主分离涡核,膨胀区,滑移线等.计算与Alvi等提出的单鳍流动的理论模型符合很好,得到了平板表面压力以及分离线、再附线等在单鳍流动中所独有的半圆锥特性.      

恒热流斜板自然对流边界层的稳定性研究

对恒热流条件下的斜板自然对流边界层的稳定性进行了分析及数值模拟.研究发现在一定的倾角和Prandtl数范围内,三维波可存在两种不稳定模态,其临界频率分别接近于重力内波和二维T-S(Tol-mien-Schlichting)波的频率.更为特殊的是,在适当条件下该三维波模态可比二维T-S波更不稳定.这一看似违反Squire定理的结果是由于控制方程中的浮力项造成的.数值模拟初步验证了三维波的临界特性.由于稳定性总是同层流至湍流的转捩过程密切相关,计算结果可为高效换热翅片的设计提供理论依据.      

直升机尾桨轴疲劳试验技术研究

尾桨轴为直升机传动系统中的关键部件,疲劳破坏为其主要的失效模式。本文结合某直升机尾桨轴的疲劳试验,介绍了疲劳试验设计、试验实施方案设计的工作内容和方法。考虑试验的考核要求,对锥齿面啮合处约束边界进行工程简化,设计了防扭结构。将试验载荷分解到作动器上进行协调加载,实现旋转弯矩和剪力绕尾桨轴轴线的旋转,模拟了真实载荷的传递和分布。根据飞行谱和使用载荷确定低周疲劳载荷谱,并给出了既满足加载要求,又易于实现的加载方法。     

开槽位置和槽道结构对叶栅性能的影响

基于利用小槽出口射流控制叶片吸力面尾缘分离气流的思想,为进一步研究开槽位置和槽道结构对叶栅性能的影响,设计了三种叶片开槽处理方案,用CFD方法对开槽前后的流场进行数值模拟分析.计算结果表明,小槽出口射流可以有效吹除叶片吸力面尾缘分离气流,控制和延缓附面层分离,从而增大了气流转折角,降低了总压损失,提高了叶栅流场稳定性,并且叶片开槽位置和槽道结构对叶栅性能也有较大的影响.     

局部抽吸对边界层流动稳定性的影响和控制研究

在边界层壁面上,设计局部抽吸结构,采用直接数值模拟的方法,获得稳定的三维基本流.在此基础上,研究稳定及最不稳定的二维扰动T-S波的时、空演化机制;进一步探讨了局部抽吸的形式、强度大小及分布结构对二维T-S波的非线性演化影响及其对增长率的贡献大小.结果表明,局部抽吸结构诱导产生的三维基本流是扰动波得以快速增长的一个关键性因素,这是由于平均流剖面的改变及展向速度的出现,增强了流体运动中的不稳定性、扩大了中性曲线的不稳定区域范围.在最不稳定的二维扰动T-S波的非线性演化过程中,由于非线性作用的不断增强,逐渐激发产生出三维扰动波及高次谐波,其三维扰动波的流向波数和频率与二维扰动波的流向波数和频率相同;同时展向速度的大小对二维扰动波的增长、流动的失稳、流向涡的形成等方面都起着激励的作用.随着时、空的不断发展和非线性作用的迅速加强,正、负相间的流向涡逐渐形成,强度逐渐增大,流向涡的影响区域也在不断扩大,涡的形状逐渐拉伸变长,并出现强的剪切层,流动开始失稳等其它机制;这些结论与文献[4、5]的结果相吻合.     

钝体头部边界层“逾越”型转捩机理研究

利用数值模拟的方法,研究了再入条件下,钝头体边界层中一种新的失稳机制.这种失稳机制无法用线性稳定性理论解释.最近的研究工作表明,流动中存在着一种与局部有限幅值扰动相关的失稳机制.它与常规的T-S波或二次失稳的过程不同.数值模拟了受空间模式的局部扰动影响的钝体头部超声速边界层的失稳过程.在模拟中,明确了不同扰动模型非线性相互作用的物理机制.      

带有压力梯度的平板边界层转捩数值模拟

对零压力梯度和带有压力梯度的平板边界层转捩实验T3A,T3C1和T3C2进行了数值模拟,计算结果与实验值吻合的较好,并得出以下的结论:压力梯度的存在对转捩起始的位置有较大的影响;M-L转捩模型能比较准确地预测转捩的发生和发展过程,但M-L一方程转捩模型因为要根据经验选取转捩起始动量厚度雷诺数,所以在通用性上相对M-L二方程转捩模型要弱;使用M-L转捩模型计算过程中,根据湍流的性质选取进口粘性比有助于正确地预测转捩的发展过程.      

平板附面层等离子体流动控制的数值模拟

通过求解电场中的拉普拉斯方程和赫姆霍兹方程的变形形式,成功地将等离子体激励对平板附面层流动的影响,以体积力向量的形式引入到NS方程之中.借助求解的电荷密度,通过将计算结果与实验数据的对比,可确定德拜长度、最大电荷密度、形状因子等可调参数的取值原则,建立起平板附面层等离子体流动控制的数值模拟方法,为将等离子体流动控制方法应用于外流及内流场中的强剪切流动控制,奠定关键性的技术支撑.      

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

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