基于KhokhlovZabolotskayaKuznetsov方程的风扇激波噪声频域预测法.pdfVIP

中国工程热物理学会 热机气动热力学 学术会议论文 编号：092112 基于Khokhlov-Zabolotskaya-Kuznetsov 方程 的风扇激波噪声频域预测法 陈鹏，李晓东 （北京航空航天大学流体与声学实验室，北京 100191） （Tel：010E-mail:lixd@buaa.edu.cn） 摘 要：建立了一种快速预测风扇激波噪声传播的频域方法，基于传统Khokhlov-Zabolotskaya-Kuznetsov (KZK)方程，描述噪声在有均匀流动的圆环形管道内的传播。将该模型用于真实风扇噪声的预测，利用计 算流体力学(CFD)方法得到流场；将CFD 结果作傅立叶变换后代入模型方程，快速求解噪声传播特性。预 测结果与，研究表明：预测方法能够捕捉风扇激波噪声的非线性传播，声衍射项使得模型能够描述激波径 向分布的影响，传播过程中波形变化还与频散作用相关。 关键词：风扇激波噪声，非线性传播，Khohklov-Zabolotskaya-Kuznetsov （KZK ）方程，频域法 1 引言 由于高涵道比航空发动机的广泛使用，新的附加噪声也随之产生。当飞机起飞或爬升时， 由于风扇叶片的展弦比较大，高转速会导致叶尖附近径向区域之线速度与来流的合速度超过 音速，即会在叶片前缘形成脱体激波；由于激波与来流通过叶片吸力面产生的膨胀波相互作 用，并且每个风扇叶片的安装角误差引起相邻激波的强度存在微小诧异，这两个因素则会导 致与风扇转子叶片一起运动的激波系结构变得相当复杂。这一复杂的激波系随风扇叶片一起 运动便形成风扇激波噪声的声源[1] 。 从20 世纪70 年代初期开始，科研工作者期望建立理论模型以预测激波噪声的传播特性。 研究者们试图通过理论手段分析激波噪声的传播。根据 Whitham[2] 的弱激波理论， Hawkings[3]、Morfey 和Fisher[4] 以及Kurosaka[5]给出了预测跨音转子激波噪声的一维传播模 型，但无法考虑发动机进出口边界条件对激波噪声的影响，同时无法考虑叶片安装角误差带 来的影响。McAlpine 和Fisher[6]等人假设转子叶片仅在叶尖处超音，从二维激波结构出发， 利用噪声传播的螺旋坐标与轴线坐标的关系得到预测风扇激波噪声的一维传播模型；该模型 构造了非规则波形的激波噪声传播，得到不同频率的谐波的传播特性，能考虑敷设声衬的进 气道对噪声传播的影响；但由于模型方程建立在叶尖处，故无法描述激波分布对噪声传播的 影响。 在航空发动机的设计过程中，风扇设计是重要环节之一。在保证风扇工作效率及裕度的 同时，需考虑激波分布对叶片加功量的影响，现代风扇叶片除了增加弯／扭设计之外，还会 在叶尖附近采用特殊造型，如前掠、后掠等，这样构造出激波沿径向的分布；因此，在预测 噪声传播时，必须考虑激波强度在径向的变化情况。 本文基于Khokhlov-Zabo-lotskaya-Kuznetsov (KZK)方程[7] ，引入来流影响，建立了航空 发动机进气道内激波噪声传播的预测模型，其中包含了激波强度的径向分布。此外，对于真 实航空发动机风扇叶片，其安装角的误差会导致激波噪声呈现不规则锯齿波形，因此利用随 机数构造不规则的锯齿波形。在均匀圆环直管道内，本文计算了噪声的传播特性，分析了激 波径向分布对传播的影响。 基金项目：航空基础科学基金（2006ZB51023 ）和新世纪人才支持计划（NCET-07-0038 ）资助项目 2 预测模型 假设航空发动机进气道为半无限长圆柱形管道，不考虑向下游传播的声波反射作用的影 ρ 响，同时也忽略端口反射。若设管道内为无旋流场，其流速可定义为V (0,0, V ) ；在以