基于谐振舵面的跨声速抖振抑制探究

跨声速抖振引起的非定常脉动载荷会造成飞行器结构疲劳甚至引发飞行事故,所以跨声速抖振的控制研究逐渐成为航空领域的热点。采用基于Spalart-Allmaras(S-A)湍流模型的非定常雷诺平均方程开展了基于谐振舵面的跨声速抖振抑制研究。首先验证静止NACA0012翼型的抖振边界和频率特性,然后分别从舵偏平衡位置、舵偏幅值、频率以及相角等角度研究了谐振舵面的控制效果。舵偏平衡位置等效于减小了翼型的有效迎角;幅值和频率对抖振抑制效果影响较大,当舵面振荡频率与抖振频率接近时发生共振现象;相角对控制效果有一定影响,在270°相角附近,升力系数幅值减小了60%。在合适的舵偏幅值、频率以及相角组合下,谐振舵面有可能成为跨声速抖振的有效开环控制策略。     

模型振动对风洞实验数据分散性的数值模拟研究

研究模型振动对风洞实验结果的影响对于修正风洞实验数据和加深空气动力学的理解具有重要的现实意义。基于κ－ω的SST两方程湍流模型,在时间域求解雷诺平均N—S方程,研究了模型振动在不同迎角下流场和气动力的变化规律。结果表明：模型振动在小迎角条件下对气动力测量结果几乎没有影响,但在大迎角存在一定分离的条件下,可以引起分离涡的非线性演化,导致气动力平均值的改变,从而引起风洞气动力测量值的分散性。结论从物理机理上解释了大迎角风洞实验的重复性和数据分散性。     

进口导流叶片对转子叶片颤振特性的影响

上游静子叶片的尾迹扰动是下游转子叶片发生强迫响应问题的主要原因,但通常认为不会对颤振特性产生影响。运用CFD技术求解非定常流场,用能量法对典型涡轮风扇NASA Rotor 67进行了气动弹性分析。结果表明:不加进口导流叶片时其叶片最有可能在第二阶模态、60°叶间振动相角的情况下发生颤振;加上IGV(进口导流叶片)后,IGV的尾迹会显著改变转子叶片的颤振特性,当转子叶片的最危险颤振模态频率接近IGV尾迹的扰动频率时,即使单独转子叶片十分稳定,但在尾迹的激励下颤振仍然可能发生;当该频率远离IGV尾迹扰动频率时,尾迹扰动同样会使转子叶片气动阻尼降低。该结果意味着尽管IGV的引入不会明显改变压气机性能,但有可能会对其颤振特性带来明显的负面影响,需要在涡轮机设计中予以考虑。     

悬停共轴双旋翼干扰流动数值模拟

基于非结构动态嵌套网格方法,通过求解三维非定常Euler方程对悬停共轴双旋翼直升机的复杂流场进行了数值模拟研究.双旋翼之间的干扰使得两者的拉力性能下降,上旋翼对下旋翼的干扰程度远大于下旋翼对上旋翼的影响程度.下旋翼处的下洗速度比单旋翼的要大得多,而上旋翼的下洗速度比单旋翼的略大一些,几乎相同.桨叶总星巨角增加,则下洗速...     

二维有限长度柔性壁面上T-S波演化的数值研究

受自然界鸟类羽毛的柔性特征启发,利用数值模拟的方式开展了柔性壁面对亚声速边界层中T-S波演化的影响研究。刚性壁面上的数值结果与线性理论吻合得很好,验证了数值方法的可靠性。在此基础上,将部分刚性壁面替换为柔性壁面,结果表明,柔性壁面可以抑制T-S波在空间上的增长,从而推迟边界层流动转捩。壁面的变形不只跟随T-S波的波形,还因为柔性段与刚性段相接的前缘和后缘引起与扰动源频率相同的更大尺度的壁面波动,壁面的实际变形由这几种波叠加而成。开展的参数研究结果表明,增大表面的质量密度对于柔性壁面衰减扰动的效果几乎没有影响;增大表面张力和增加底部支撑的弹性系数可以增加壁面的刚性,减小壁面变形的幅度;增加阻尼可以抑制柔性段前后缘产生的大尺度壁面波动的传播,而对跟随T-S波的变形影响不大。总体上,柔性壁面的变形程度越大,其扰动的抑制效果越强。      

弹性振动对翼型失速迎角附近流场的影响

通过求解雷诺平均Navier-Stokes方程得到气动力,结合翼型振动方程,计算了翼型不同迎角下的动态过程,分别通过层流和湍流情况的计算,重点研究了弹性振动对翼型失速迎角附近流场的影响,研究结果表明,在中低雷诺数、翼型具有弹性振动的情况下,翼型的失速迎角会比传统定常意义上的失速迎角提前出现,为长期以来数值计算得到的失速迎角与风洞实验、飞行试验结果的不同给出了一种物理解释。     

复杂外形下的无网格算法研究

研究了二维、三维复杂外形下的无网格算法。在一种布置点云方法的基础上,发展一种曲面拟合的重构方式构造流场物理量,并应用于AUSM+_up格式计算欧拉方程的数值通量;计算采用了隐式时间推进,并引入当地时间步长和残值光顺等加速收敛措施。通过对某三段翼型低速流动、M6机翼跨音速流动、某全机跨音速流动进行了数值模拟,表明本文发展的无网格算法能有效地模拟复杂外形的无粘绕流。     

基于嵌套网格的颤振激励系统数值计算研究

基于非结构动态嵌套网格方法,通过求解可压缩非定常Navier-Stokes方程,对一种新型颤振激励系统的工作特性进行了数值计算研究。旋转圆筒的开缝圆心角过大或过小、位置距离小翼后缘过远或过近都会使得激励效果下降,而适当增大旋转圆筒直径可以增加激励的效果,另外,开缝装置的截面形状也会对激振力大小产生一定的影响,在本文的计算状态下,采用月牙形截面时激励效果会较好一些。当开缝50%时,选择靠近翼根处的激励效果会较好一些。     

二维分离流中的NS方程数值计算方法

本文采用显式多步龙格－库塔方法求解了二维可压非常完全ＮＳ方程，为了使计算适用于一般的风格，采用有限体积法进行几何离散，研究定常问题计算中的一些加速收敛措施，如：当地时间步长、残值光顺技术等。讨论了模拟分离流问题中的技术问题，计算结果表明本文的计算结果与国外文献是吻合的。     

两种跨声速气动弹性问题分析研究

在非结构运动网格基础上,采用中心有限体积法进行空间离散和双时间方法进行时间推进求解非定常欧拉方程.通过与气动力方程的联立求解,在时域内用四步龙格-库塔方法求解结构运动方程.分析和研究了二维嗡鸣和三维机翼颤振这两种跨声速非线性气动弹性问题.二维嗡鸣问题的研究考虑了翼面-舵面系统的缝隙间网格运动、缝隙对嗡鸣的影响和扰流片对嗡鸣的抑制.耦合多自由度Lagrange结构运动方程数值模拟了三维机翼的颤振问题.通过跨声速标模算例AGARD445.6机翼的颤振计算,计算的颤振临界速度与实验值有5%左右的误差,验证本方法的正确性.由于本方法是在对外形具有良好普适性的非结构动网格基础上完成的,具有良好的工程实践价值.