三角翼受迫滚转运动气动特性的数值模拟研究.pdfVIP

第26卷第 4期 气 动 研 究 与 实 验 Vo1．26．No．4 2008年 l2月 AERODYNAMICRESEARCH EXPERIMENT Dec．．2008 三角翼受迫滚转运动气动特性的数值模拟研究 高 娜 陈 宝 赵庆贺 (哈尔滨 150001) 摘 要：基于动网格的技术，利用有限体积法对N—s方程进行数值离散，选用湍流 模型。对 65。大后掠三角翼的受迫滚转运动进行了数值模拟，在此基础上，对滚转运动的 气动力特性和流场结构进行了分析。计算结果表明：滚转运动时．三角翼上表面涡分布的 非对称性产生横侧方向的偏航力矩和滚转力矩 ，偏航力矩、滚转力矩系数滞回曲线呈现 典型的 “单环”形态，俯仰力矩、法向力滞回曲线呈现典型的 “单 8”形态。 关键词 ：三角翼 气动特性 受迫滚转运动 涡结构 1．引言 在三角翼绕流这个典型流场环境中，包含了丰富的空气动力学现象。这些现象的研 究对空气动力学理论的发展和完善具有重大意义，对提高现代飞行器性能具有重要价 值。因此，国内外大量学者对其进行了研究。Abdelhamid研究了缩减频率对俯仰运动的影 响，发现缩减频率增加，升力和阻力的峰值增大。Menzies发现：65。三角翼在固定攻角下做 有限振幅的滚转运动时升力系数有较大提高。Kandil在模拟65。三角翼受迫俯仰滚转耦 合运动时发现：当滚转运动的缩减频率提高到2盯时，涡的破裂将消除，从而使升力有所 提高。国内，阎超等人对三角翼的数值模拟进行了综述，杨立芝和高正红研究了76。后掠 三角翼做俯仰振荡时缩减频率对气动力的影响，发现随着机翼俯仰加快，机翼大迎角抵 抗漩涡非对称破裂的能力增强。但对于大后掠三角翼，随着迎角的增大，三角翼将会出现 摇“滚”这一现象。本文针对65。大后掠三角翼进行了数值模拟，研究受迫滚转运动时的气 动特性、涡结构、表面压力分布及表面流态。 2．数值方法、边界条件和计算网格 物体做非定常运动，且网格刚性地固连于物体，此部分随物体一起运动，另一部分不 动，两部分通过 interface连接。在惯性cartesian坐标下，偏微分形式的非定常N—S方程 组可写作： 2 三角翼受迫滚转运动气动特性的数值模拟研究 第4期 +v( )：v( )+ f1、 、 其中：西为通用变量；r为广义扩散系数；S为广义源项，其详细表达式参考文献 [2】。式 l中各项依次为瞬态项、对流项、扩散项和源项。 对流项采用改进离散方程截差的QuICK格式。扩散项采用中心差分格式离散，湍流 模型采用 一o【5s模型。 三角翼表面边界条件取无滑移的绝热固壁，远场运用压力远场条件，即给出来流的 静压、静温与马赫数，值分别为92579．165Pa、298．05K、0．29。 本文所用的外形为65。后掠三角翼，根弦长为622mm，展长为580mm，厚度为9．5mm， 其详细构型见文献 5【】。文献 3【、4】的研究发现：计算网格的拓扑结构对计算结果有明显 的影响。0型网格由于引进了奇性轴，可能使得解的质量下降，故本文采用c—H型网格， 同时在翼面和涡核附近适当地进行了加密，以便能更好地模拟流场。第一层网格距物面 的距离为 l0—5倍根弦长。沿弦向、周向及法向的网格点分布为 88×99×80。 O2 妄0 、- ， )--