阵风与翼型干涉的高阶数值模拟.pdfVIP

中国J：程热物理学会 热机气动热力学 学术会议论文 编号：062045 阵风与翼型干涉的高阶数值模拟宰 郝璇，李晓东 北京航空航天犬学能源与动力工程学院 北京航空航天大学407信箱100083 电话：010 E-mail：1ixd@buaa．edu．cn 摘要：通过求解二维线化欧拉方程对阵风与Joukowski翼型干涉发声问题进行数值模拟．分别采用频 散相关保持格式(DRP格式)和优化的4／6低频散低耗散龙格库塔法(LDDRK4／6)对方程进行空间离散 和时间推进：在进出口边界采用合适的无反射边界条件。将阵风与Joukowski翼型干涉的模拟结果与 Scott的精确解进行比较，符合较好。在此基础上研究了阵风与不同厚度Joukowski翼型干涉发声现象。 数值结果表明，在高频范围内，随着翼型厚度的增加，阵风／翼型干涉噪声的强度将减小． 关键词：计算气动声学，阵风／翼型干涉噪声，频散相关保持格式，龙格库塔法 1引言 一 在当今社会，噪声污染是危害社会的主要环境污染源之一。随着民用飞机日益普及， 航空噪声问题日益引起世界范围的关注。对于民用直升机来说，必须降低直升机近场与 远场噪声以满足乘客及周罔环境的要求。而对于军用直升机来说，则必须考虑降低直升 机的声学可探测性以满足战争需要。涡桨干涉(Blade．Vortex Interaction，BVI)噪声是一 种非常强的脉冲噪声，一旦出现就成为直升机的主要噪声。因此如何降低BVI噪声对直 升机设计十分重要。 BVI噪声是一种宽频噪声，．经过傅立叶分解后可得到连续频谱，对应的每一个单频 噪声分量都可以看作足由单频阵风与翼型干涉产生。因此在实际研究中，研究者多采用 阵风／翼型干涉模型作为BVI噪声的简化模型。除此之外，阵风／翼型1二涉可以作为航空 发动机内部来流扰动与进口导向器或风扇尾流与静子之间的干涉模型以及飞机飞行过 程中来流扰动与机翼之间的干涉发声问题的简化模型，具有广泛的应用背景。 CAA方法是一种高精度、低频散、低耗散格式，与传统CFD方法相比，可以快速准确 得获得计算结果，因此较为广泛的应用于求解声传播问题。1996年，Agarwal＆Huh 在涡／翼相互作用问题上对CAA进行了初步运用[1]。1997年，Lockard对二维零厚度平板 翼型与阵风干涉发声问题进行了数值模拟[2]。2003年，Golubev等采用空间一时问映射 分析的方法，对舴定常、非线性阵风／翼型干涉发声问题进行了数值研究[3]。2004年， Wang＆Himansu等通过空问一时问守恒元与求解元的方法对阵风／翼型干涉发声的标准 问题进行了数值模拟[4]。Hao＆Li采用CAA方法模拟了阵风与平板翼型干涉并研究了 翼型攻角对阵风／翼型干涉噪声的影响【5】。本文将在此基础上研究翼型厚度对阵风／翼型 干涉噪声的影响。 本文的内容主要仃以下两方面：首先使用阵风／翼型干涉标准篼例对数值模拟方法 进行校核，在此基础上研究翼型厚度以及阵风入射角度对阵风／翼型干涉噪声的影响。 。国家自然科学基金与航空基础科学基金(05C51030)资助项目 166 第二节介绍了所使用的数值方法；第一节给出阵风与不J一厚度Joukowski翼型十涉发声 的数值模拟结果及分析：第四节对本义的-T作进行了总结。 2数值方法 2．1主控方程 阵风／翼型干涉发声问题中所加入的阵风相对于平均流场为小量。采用小扰动假设， 叮得到二维线化欧拉方程。形式如下： 署+Ax罢+▲r万tgU+(Bx+Br)u=s+sr㈩ 其中：U= S= JD p ．p p p U ， p 0 一a万 ”瓦 ．．0F “一0x 0 AJ= 万 一P 0 0 o 万 。三 p O 0 万0 0万0历 ^=去 其中，歹，蟊玩矿分别为平均流的密度、 速度扰动。P、P分别密度扰动和压力扰动。 的声源滤波源项[6]。 小巨 0 万 矿 O O 矿 0 0佰 万0 0 万 O 0 矿 0 O 万0，y订 0 O 1 ●__- |ID 矿 ’By：昙 oy 矿0万0 0 O 订 0 0 0 矿 O 0 0万7矿 压力、x向的速度和y向的速度。U、1，为 方程右端Sf项是为了抑制不稳定波所加 2．2主控方程的离散 本文对主控方程空间导数项的离散采用频散相关保持(dispersion-relation·preserving scheme，DRP)格式E7]；时间导数离散采用四六步交替的优化的低频散低耗散的龙格库 塔法(optimized low-dissipation and dispersion Runge·Kutta scheme，LDDRK4／6)[8]。 2．3边界条件 本文所进行的数值模拟工作在计算域的出流边界选用了T柚和Webb[7]提出的出