扰动风作用下多段翼型流动的数值模拟及流动控制 Numerical simulation of gusts over multi-element airfoils and their flow control.pdfVIP

第31卷第6期 空气动力学学报 V01．31。No．6 2013年12月 ACTAAERoDYNAMICASINICA Dec．，2013 文章编号：0258—1825(2013)06—0710—08 扰动风作用下多段翼型流动的数值模拟及流动控制 李 亮，李孝伟 (上海大学上海市应用数学和力学研究所，上海200072) 摘要：运用有限体积法、双时间推进法、动态嵌套网格技术、网格速度法等数值技术实现了多段翼型在扰动风作 用下的非定常流场的数值计算，并研究了襟翼摆动对扰动风下多段翼型上的气动力的作用和控制规律。首先，运 用NACA0006翼型阵风响应的理论解对本文数值方法的可靠性进行了验证。然后，对带30％襟翼的Gw一2翼型 在具有垂直于来流方向的周期性扰动风环境下的气动力性能进行数值计算，结果表明：多段翼型的升力随时间亦 呈现周期性变化，其周期变化频率与扰动风频率相同，升力的波动幅值与扰动风幅值呈正比关系。进一步地，通过 改变襟翼的摆动幅值和摆动相位的办法，研究了襟翼运动对扰动风的作用和控制规律，结论是：针对一定频率的扰 动风，运用适当幅值和相位的襟翼摆动加以控制，能将有效地控制多段翼型上因扰动风引起的升力波动 关键词：多段翼型；扰动风；流动控制；双时间推进法；动态嵌套网格 中图分类号：V211．3 文献标识码：A 析，大大提高了阵风响应计算的效率。以上研究均表 0 引 言 明，在阵风和扰动环境下，飞行器的气动性能会发生 现代飞行器设计中必须考虑各种特殊的环境，扰 严重的波动，甚至引起飞行的不稳定性。 动风作为其中一种不可避免的因素，已经在当今飞行 从飞行的角度来说，人们期望这样的波动或不稳 器研究中引起了高度重视。 定性尽量少的发生，或者能采取一定的办法来抑制扰 从20世纪50年代到现在关于飞行器在扰动风动风的负面作用。从流体力学的角度来说，抑制扰动 和阵风作用的响应的研究已经很多[1。8]，主要集中在 风作用的最有效的办法是进行流动控制。 试验和CFD数值计算方法。在实验方面，对于模拟 对于目前常见流动控制的研究方法[1”12]：吹气和 阵风的形式是件很重要的事，目前一般是通过阵风发 吸气、微射流、零质量射流、等离子体、磁流体动力、涡 生器(或者加入一些能引起流场扰动作用的控制体) 流发生器、机翼上翼刀运动等；主要集中在对流场中 产生扰动风，由于受到阵风发生器参数与阵风形式对 出现流动分离进行控制及机翼增升减阻的研究；也有 应等问题，试验中采取的阵风形式的研究比较单一和 通过采取飞行器上机翼的运动对流场进行控制的研 笼统。在数值计算方面，文献[7—8]中发展了一种网究。从发表的文献来看，El前对飞行器在扰动风作用 格速度方法，该方法以网格速度代替扰动运动的方 下的流动的控制研究很少。 法，既避免了数值振荡导致的计算不稳定[9]，又得到 基于此，本文拟对多段翼型在扰动风作用下的非 了精确的数值结果。接着杨国伟等人发展了一种基 定常流动及其控制进行研究。由于机翼的机械运动 于CFD与ARMA降阶模型结合的阵风响应分析方相对于采取吹气和吸气、微射流、零质量射流等更具 法，并通过分析比较知道了正弦阵风响应的参数辨识 有改变流场区域的能力，因此本文采用襟翼摆动的方 模型与CFD模型直接计算结果拟合程度最好的结 法进行流动控制。在襟翼的摆动过程中，襟翼和主翼 论，建立的ARMA模型能用于任意形状阵风响应分之间存在相对运动，所以，在网格布局时采用了动态 *收稿日期：2012—07—06；修订日期：2012—12-21 基金项目：国家自然科学基金；上海市自然科学基金(13ZRl41570