前向离心风机叶轮通道中的流动扩压度较大,容易产生气流分离, 特别是叶轮后部的吸力面附近分离现象严重, 导致叶轮出口周向气流的分布很不均匀。由于叶轮转速较高, 不均匀的叶轮出口气流会与周围非轴对称分布的静止部件产生强烈的非定常干涉, 这是产生离心风机气动噪声, 特别是离散噪声的主要原因之一, 并对风机的气动性能带来不利影响。国内外对离心风机气动性能和气动噪声的研究大都着眼于改变蜗舌的倾角、蜗壳形状、叶片形状及数量等方面 。近年来国内外的一些研究表明, 合理的布置串列叶片能使离心叶轮出口流动参数沿周向分布比采用单列叶片更加均匀, 对降低叶轮与静止部件间的 非定常冲击损 失及气动 噪声有利 。到目前为止, 串列叶片对离心风机整机的气动性能和气动噪声二者影响的研究还比较匮乏。另外, 气动性能和气动噪声还常常存在着相互制约的关系,一方的提升有时会带来另一方的下降 , 能够同时改善离心风机气动性能和气动噪声的方法是当前研究所追求的主要目标。因此, 本文在综合考虑气动性能和气动噪声的前提下, 采用串列叶片结构对某单列叶片式前向离心风机进行了改进研究。本文采用数值方法对串列叶片的相对长度和相对周向位置两个参数对该前向离心风机气动性能和气动噪声的影响同时进行了研究。使用 F L UENT软件计算风机内部的非定常流动来获取流场和声源信息, 分别利用效率公式和 FW H 方程求得风机的-效率和 离 散 噪声, 并 采 用响 应 面 方 法 ( ResponseSurface Met hodo logy, RSM ) 拟合得到两个参数与风机气动性能和气动噪声的函数关系式, 通过优化分析, 寻找两个参数的*佳组合点, 在保持风机原有良好气动性能的情况下, 使气动噪声降低, 为串列叶片式离心风机的优化设计提供参考。

前向离心风机的几何参数

被改进的原单列叶片式前向离心风机( 后文中简称原风机) 叶轮进口直径 D1 = 156m m, 叶片出口直径 D2 = 400mm, 无叶旋转扩压器出口直径 D3 =460m m, 叶片数 Z = 12, 叶片进口安装角 B = 38b,1A出口安装角 B = 126 叶轮进口宽度 b1 = 70m m,2Ab,叶轮出口宽度 b2 = 36mm , 蜗壳宽度 B = 64mm, 设计转速 n = 2900rpm。图 1 给出改进后的串列叶片式离心风机结构示意图。串列叶片结 构及相关参数如图 2 所示。 2图中, 若把后排短叶片去掉, 将前排长叶片加上图中虚线部分, 就得到原风机的叶型。后排短叶片的叶型与改进前叶片后部的叶型完全相同。符号 a 表示后排叶片的弦长, b 表示改进前原单列叶片的弦长, ls 表示前后排叶片 重叠形成的缝道长度, H为后 排叶片相对前排叶片沿叶轮旋转方向转动的角度。定义叶片相对长度因子€ a/ b; 叶片相对周向位置因子 €l=H= H $, $/ H H= 30 为叶轮单个叶道所占的周向角度。b文献[ 3 4] 指出短叶片在靠近长叶片压力面一侧,-对改善离心风机气动性能和气动噪声更有利。本文在保持 ls / b = 1/ 20 不变, 分别选取相对长 度因子l为€ 2/ 20、 20、 20、 20, 相对周向位置因子 €3/4/5/H为1/ 15、 15、 15、 15, 并根据样本点的设计方案,2/3/4/计算在不同组合下的风机气动性能和气动噪声。

数值计算方法

三维非定常流场的数值计算方法

采用 F L UENT 软件计算风机内部的非定常流动来获取流场和声源信息。风机计算区域主要分为三个部分: 进气段、叶轮( 包括旋转扩压器) 和蜗壳。为了提高收敛速度, 忽略了进风口和旋转叶轮套接部分的内泄漏流动。为适应风机内部结构的复杂性,采用多块网格生成方法生成高质量网格。考虑到不同流动区域的不同重要性, 对叶轮内部、蜗舌壁面附近的网格节点进行了加密控制和非等距处理。为方便, 取叶轮旋转轴中心线为 Z 坐标轴, 进口气流沿 Z轴正向流向 Z 轴负向, 轮盘内壁面为坐标 Z = 0 平面。各部位的具体网格节点数分别为: 进气部分约 6万个, 叶轮部分约 61 万个, 蜗壳部分约 23 万个。整机网格节点总数约 90 万个。风机的计算网格模型如图 3、 4 所示。图

- S 方程, 湍流模型采用 RealizableJ E模型, 取标准-壁面函数。由于风机内部流动马赫数很低, 认为气体不可压缩且粘性系数为常数。离散方程采用隐式分离方法求解, 压力修正采用 SIMP LE 算法。对流项采用二阶迎风格式离散, 扩散项采用二阶中心格式离散, 时间项离散采用二阶隐式格式。由于叶片在叶轮 内沿圆周均布且认为转速恒定, 所以叶轮流道和蜗壳之间的耦合呈周期性关系,在计算中时间步长确定为$t =60/nK Z( 1)式中: K 为一个非定常计算周 期的时间步数, K =30; n 为叶轮转速, n = 2900rpm; Z 为叶片排数, Z =12。通过上式计算得到时间步长为$t = 5 747 @ 10 s1计算给定进口速度边界条件, 依据流量计算管-5道进口速度, 且假设进口速度均匀分布; 给定蜗壳出口压力边界为标准大气压; 壁面采用无滑移条件; 计算中取进口体积流量 qv = 27 7m 3 # min- 1 。1321气动性能和气动噪声的计算方法效率是衡量风机气动性能的主要参数之一, 风机效率的计算可确定为

             ( P2 - P1) qv( 2)G=MX其中: P1 和 P 2 为风机进口和出口的总压; M 表示叶轮施加给流体的作用力对 Z 转轴的力矩; X 表示叶轮旋转的角速度。计算中上述各参数均采用国际标准单位。由于在不同时刻( 叶片和蜗舌处在不同相对位置) , 计算所得的效率不同, 为了便于比较, 下文中所给出的效率为 360 个非定常时间步( 叶轮旋转一周) 的时均值。声压级是衡量风机噪声大小的常用物 理量之一, 本文运用三维非定常流场数值计算的结果, 采用FW H 方程对风机的离散噪声进行计算 [ 5] , 其表达-式为1 52255( 2 2 - ý )p =[ QD f ) ] -([ F iD f ) ] +(c0 5t5t5x i5 2 [ T ij H ( f ) ]5x i 5x j利用自由空间格林函数求解 FW H 方程, 其完-备解是由面积分和体积分组成的。面积分部分表示的是单极子和偶极子声源, 以及部分四极子声源的作用, 而体积分部分表示的是在声源面以外的四极子声源。对于所研究的前向离心风机, 其内部流动属于低亚音速不可压缩流动, 当声源区包含在声源面内时, 单极子和四极子声源对声场的贡献可以忽略。

常州市（无锡）文顺风机有限公司产品齐全，主要产品有：高、中、低离心风机、引风机；各种轴流风机；水泥立窑专用高压离心风机；化铁炉高压离心风机；化铁企业替代进口的不锈钢风机；空调专用风机；各种除尘、耐磨、防腐、非标等特种风机近三十个系列，三百多种规格的产品。

公司各类通用、专业设备齐全，工艺先进，技术力量雄厚，性能测试运用微机电脑处理。

常州市文顺风机有限公司

Changzhou Wenshun Wind Machine Co., Ltd.

电话：0519-88726638

传真：0519-88726638

手机：13601517069

网址：www.wsfengji.com

邮箱：ws@wsfengji.com

地址：江苏省常州市横林镇新东方村工业园（殷坂）