离心风机叶片噪声数值研究I：声源分析.pdfVIP

中国工程热物理学会 流体机械 学术会议论文 编号：097002 离心风机叶片噪声数值研究 I:声源分析 毛义军 祁大同 (西安交通大学能源与动力工程学院，西安，710049) Tel：029Email: dtqi@mail.xjtu.edu.cn 摘要：采用混合计算气动声学方法分析了离心风机叶片辐射噪声。首先通过求解风机内部粘性非定常 流动获得叶片声源信息，然后采用频域形式的FW-H 方程对叶片辐射噪声进行了分析。结果表明：叶 片辐射噪声表现出典型的偶极子特性，其声源主要位于叶片尾缘附近区域；叶片声源强度主要受叶片- 蜗舌位势干涉影响，当叶片尾缘掠过蜗舌附近区域时，叶片辐射噪声最大。 关键词：离心风机，噪声，计算气动声学，声源 0 前言 随着国民经济和社会发展，工作和居住环境噪声越来越受到关注。2008 年 10 月 1 日国家颁布实施了三套新的强制性标准，对工业企业和居民生活环境噪声的排放要求做 了严格的规定[1-3]。离心风机被广泛应用在冶金、电力、汽车、空调、供暖等工业生产 和居民生活场合，是很多应用场合重要的噪声源。据作者利用百度新闻检索发现，仅在 2007 年10 月1 日至2008 年10 月1 日一年时间内国内新闻报道风机噪声扰民事件就达 到39 起。本文采用混合计算气动声学方法(HCAA)对某离心风机的叶片噪声产生机理和 声源位置进行了分析，本文工作为下文分析蜗壳散射对叶片噪声传播的影响提供准备。 1 非定常流动数值模拟 1.1 计算模型 以 T9-19No.6.3A 离心风机为例，数值预测该离心风机的叶片噪声。T9-19No.6.3A 离心风机结构如图1 所示，其主要结构参数如下：叶轮进口直径D1=246mm ，叶片出口 直径D =630mm ，无叶旋转扩压器出口直径为D =725mm ，叶片数为Z=12 ，叶片进口安 2 3 装角 β =38º ，出口安装角β =126º 。计算工况参数叶轮转速N ＝1450r/min，进口速度 1A 2A v=12.7m/s 。叶片通过频率f BPF ＝NZ/60=290Hz 。 图1 离心风机结构示意图 首先通过求解非定常RANS 方程，获得气动噪声源强度，然后通过FW-H 方程求解 叶片辐射的气动噪声。采用二维计算模型，求解非定常RANS 方程；湍流模型采用标准 k-ε模型，取标准壁面函数；由于风机内部流动马赫数很低，认为气体不可压缩且粘性系 数为常数；离散方程采用隐式分离方法求解，压力修正采用SIMPLE 算法；对流项采用 二阶迎风格式离散，扩散源项采用二阶中心格式离散，时间项采用二阶隐式格式离散。 非定常流动计算取时间步长 -4 ∆ t=1.1494 ×10 s，即叶轮每转 1°所需的时间，叶轮旋转 一周需要的时间步数T=360 。 本文采用二维计算模型，主要基于两点考虑：1) 简化计算模型，减少计算成本， 但能够模拟基本的流动和声传播现象；2) 离心风机的气动设计大多数都是采用一维设 计方法[4,5]，采用二维模型预测噪声能够比较方便的和目前的气动设计方法相衔接，以 便在设计阶段优化噪声控制方案。但是二维模型对某些物理信息的三维特征可能会遗 失，本文的方法也可以推广到三维计算模型。 1.2 叶片位置和气动力计算 为了节省计算成本，在风机内部非定常流场计算完毕后，提取叶片声源信息时作了 如下假设处理：1）叶轮叶片沿周向呈均匀分布；2 ）叶片在旋转一圈过程中，叶片上任 意位置处的静压呈周期性变化，且不同叶片在同一周向位置时的静压相等。 基于以上假设，则可以对声源信息做简化处理如下。程序中只存储了第1 个叶片上 Nnode 个网格节点在初始时刻的位置Y ( Nnode) 和单位法向矢量(UMV) N