李 亮,张 斌,付 玲,肖庆麟,滕新科,张 劲(.中联重科股份有限公司,湖南长沙4003;2.建设机械关键技术国家重点实验室,湖南长沙4003)
扫路车用风机内部涡流特性与压力脉动分析
李 亮1,2,张 斌1,2,付 玲1,2,肖庆麟1,滕新科1,2,张 劲1,2
(1.中联重科股份有限公司,湖南长沙410013;2.建设机械关键技术国家重点实验室,湖南长沙410013)
以某型号扫路车用风机内部涡流和压力脉动特性为研究对象,建立适宜于扫路车用风机内部湍流流动的C F D数值模拟方法,仿真与试验结果的对比验证了所提方法的合理性和可靠性.分析表明:扫路车用风机靠近蜗舌区域的叶轮流道会产生强烈的漩涡流动乃至回流,该流动特征是风机气动噪声的主要来源;蜗舌靠近出口始终存在一个低速、涡量值较小的漩涡区,蜗舌结构形式及叶片与蜗舌间距是影响风机气动噪声特性的关键参数;叶片压力面比吸力面上的脉动周期性特征更强烈,正对叶轮流道的蜗舌表面压力脉动最强.
声学;扫路车风机;涡流;蜗舌;数值模拟
随着国家环境战略的强势推进,作为城镇环卫作业主力军的扫路车的高噪高能耗特性与日益提高的低噪节能环保要求之间的矛盾日益突出.扫路车用风机是扫路车气力系统的关键零部件、动力源和主要噪声源[1],其降噪需求也日益迫切.
目前,针对扫路车用风机气动噪声的发声机理、声源特点以及声传播特性等研究仍显不足.因此,开展扫路车用风机气动噪声分析是扫路车技术领域中亟待探索的关键研究方向.
本文以某型号扫路车用风机内部涡流和压力脉动特性为研究对象,对扫路车用风机在额定工况下的内部流动进行三维数值仿真.首先对风机内部涡流特性进行分析,然后结合时域分析方法对流场内部压力脉动进行分析,以揭示扫路车用风机气动噪声产生的机理、气动噪声源的位置及噪声类型,为扫路车用风机气动降噪提供有益的信息.
本文数值模拟的扫路车用风机主要由进口段、出口段、叶轮以及蜗壳构成,具体结构如图1所示,其主要参数详见表1.
图1 风机结构示意图Fig.1 Schematic diagramofcentrifugal fan
表1 风机主要结构参数Tab.1 Main structure parameters of centrifugal fan
在蜗壳圆周上每隔30°取一个观察点Pi,总共选取11个观察点;在蜗舌上选取3个测试点A (P1)、B(P2)、C(P3);在靠近蜗舌附近的3个叶片尾端上取6个观察点,压力面上的点为Di,吸力面上的点为di.所有点均处在离心风机叶轮中截面上,部分观察点具体位置如图2所示.
图2 部分测试点位置示意图Fig.2 Location drawing of some viewpoints
2.1控制方程及湍流模型
扫路车用风机的内部流动是三维黏性不可压缩流动,遵循物理守恒定律,控制方程包括连续性方程、动量守恒方程以及湍流输运方程.
连续性方程:
动量守恒方程:
上述式中:ρ为密度;p为压强;ν为运动黏度;fi为体积力;μe为湍流黏性系数;μ为分子黏性系数,μt为涡流黏性系数.
由于扫路车用风机叶轮高速旋转和强曲率效应的影响,其内部流场结构极其复杂,包含有叶片尾缘的分离尾迹流动、漩涡脱落以及泄漏流动等复杂流动现象,因此本文采用工程上应用广泛的R NGk-ε湍流模型.该模型通过对湍流黏性进行修正,考虑了旋转和曲率效应对流动的影响,与标准k-ε模型相比可以较好地处理高应变率及流线弯曲程度较大的流动[2].在R NGk-ε模型中,k方程和ε方程分别为:
2.2网格生成
采用三维建模软件进行实体建模,运用专业网格划分软件生成网格.为了使扫路车用风机内部流动充分发展,在风机进、出风口处各添加了一段延长管道.在划分网格时,考虑到扫路车用风机结构的复杂性以及网格的生成质量,采用分块划分网格,各个区域单独生成合适的网格,相邻的区域共用一个面.划分网格时首先进行了网格无关性计算,确保网格数量对计算精度不造成影响,确立了网格数约为1.14×106,2.4×106,6.33×106的3种网格.采用不同网格数进行计算所对应的性能曲线与试验性能曲线对比如图3所示,考虑到数值预估的精度,本文选取网格数量约为6.33×106.蜗壳和叶轮采用适应性较强的非结构化四面体网格,并对曲率变化较大的面进行了网格加密[3](如图4所示).模型各部分的网格数量分别为:进口段814339,叶轮区域2717480,蜗壳区域2067721,出口段729813,风机网格总数为6329353.图4叶片和蜗舌表面网格
图3 全压与流量关系曲线Fig.3 Relation of total pressure and flow
图4 叶片和蜗舌表面网格Fig.4 Grid distribution of blade and scroll
2.3近壁面处理方法
近壁面区域的流动由于其自身复杂特性及重要作用,一直是湍流模拟的热点和难点,近壁面的处理对于湍流计算的精度也具有重要影响[4].目前,应用较广的湍流数值计算方法是R AN S方法,它需要引入湍流模型才能使雷诺方程封闭可解.常用的R NGk-ε湍流模型是高雷诺数湍流模型,这些模型均是对于充分发展的湍流才有效.而近壁面附近区域的流动雷诺数较低,湍流的脉动影响不如分子黏性影响大,湍流发展并不充分,故高雷诺数湍流模型不能直接适用于近壁面附近区域的流场计算.为了使这些高雷诺数模型计算能够延伸到壁面,可使用两种方法为近壁区建模:标准壁面函数法和增强壁面处理法.
采用不同的壁面处理,对网格有不同的划分要求.考虑到扫路车用风机强旋转和强曲率效应,结合标准壁面函数法要求,本文y+的范围确定为11. 25≤y+≤300[5].
本文近壁面流动的处理方法采用标准壁面函数法,叶片和蜗舌表面的y+分布如图5所示,从图中可以看出本文数值计算的网格划分能满足上述要求.
图5 叶片和蜗舌表面y+分布Fig.5 y+distribution of blade and scroll
2.4方程离散与求解
计算时采用“多重坐标系”(MR F)模型耦合动静部分区域,即叶轮区域为旋转区域,采用旋转坐标系,给定相应的旋转速度;其余区域为静止区域,采用静止坐标系.离散方法采用有限体积法,压力速度耦合关系采用S I MP L E算法;湍动能、耗散率、动量方程的离散采用二阶迎风格式,使用速度入口及压力出口边界条件,其中进口速度为24.76m/s,出口压力为1个大气压.
叶片在叶轮内沿圆周均匀分布,且转速恒定,计算时间步长为6.41×10-5s,这个时间步长的选择与叶轮的转速有关,即每经过360个时间步长,叶轮旋转1周(时间为T).这个时间很短,能够获取叶轮所需要的旋转角度,每个时间步长迭代次数为20,叶轮转速为2600r/m i m.
对扫路车用风机内流场先作定常计算,取残差小于1×10-4时计算收敛.然后将定常计算的结果作为非定常计算的初始化数据,进行非定常计算.
3.1方法验证
本文对扫路车用风机在不同工况下的流动情况进行了预测,并与试验结果进行对比分析.按照G B/T 1236—1985《通风机空气动力性能试验方法》的要求对扫路车用风机进行试验,试验装置如图6所示.
图6 试验装置Fig.6 Testing device
本文采用无量纲流量系数和全压系数对结果进行处理[6],各工况下效率和全压系数的仿真与测试结果的对比如图7,8所示.从图中可以看出,数值模拟结果的与试验结果变化趋势一致,吻合良好,并且全压内效率曲线和全压系数曲线的最大误差在4%以内,准确预测了扫路车用风机的性能变化情况,为后续扫路车用风机内流场特性和压力脉动分析奠定了基础.
图7 试验和模拟的效率曲线Fig.7 Efficiency curve of experiment and simulation
图8 试验和模拟的全压系数曲线Fig.8 Total pressure coefficient curve of experimentand simulation
3.2叶轮内部涡流特征分析
由于流体的黏性作用以及旋转部件与静止部件的动静干涉,扫路车用风机内部流场呈现出全三维、非定常的流动特征.这种流动特征引起压力场的压力脉动,导致扫路车用风机气动噪声的产生.
图9为扫路车用离心叶轮1个叶片流道内同一截面、不同时刻的瞬态流线图.从图中可以看出,在长叶片进口吸力面处,由于进口冲角过大,使得长叶片吸力面上产生流动分离,整个流道流动情况恶化,该过程一直延续到叶轮出口,并在叶轮流道出口靠近蜗舌附近形成一个回流区,该回流区不断地与蜗舌表面相互作用,使得蜗舌表面的压力产生脉动,形成扫路车用风机气动噪声的一个来源.同时,在短叶片吸力面处,流道中一直存在一个大的漩涡区,几乎堵塞整个流道,该漩涡区域与回流区的流动相互干涉,进行能量交换,使得漩涡在1个叶轮旋转周期内,周期性地经历产生、发展、耗散过程,该过程加剧了扫路车用风机内部压力的脉动,成为扫路车用风机内部涡流噪声的主要来源[7].上述分析结论表明,在扫路车用风机内部,由于蜗壳结构的非对称性作用,靠近蜗舌区域的叶轮流道会产生强烈的漩涡流动乃至回流,该非定常流动增加了叶轮气流作用在蜗舌表面的脉动作用力,增大了扫路车用风机的气动噪声.
图9 叶轮与蜗舌相互作用瞬态流线图Fig.9 Transient streamlines of the impellerand volute tongue interaction
如图10所示,在蜗舌区域,始终存在一个低速漩涡区域,该漩涡周期性地经历产生、发展、壮大、脱落过程.漩涡的产生过程使得蜗舌表面受到非定常脉动力的作用,从而产生气动噪声.在叶片的尾部,由于漩涡的脱落,使得该处的速度场极其复杂,流体微团受到拉伸、挤压和扭曲变形作用[8-9],从而产生较大的涡量值.湍流涡声理论[10]认为,气动噪声与漩涡的拉压和变形密切相关,涡量是关联气动噪声源的关键湍流量[11].上述分析表明在扫路车用风机靠近出口的蜗舌区域,始终存在一个低速、涡量值较小的漩涡区,而叶片尾部较大涡量值的存在加剧了风机气动噪声的产生.因此,蜗舌结构是影响气动噪声特性的关键参数,在其设计中需要着重考虑才能使得其气动噪声性能达到最佳.
图10 蜗舌附近速度大小随时间变化Fig.10Changes of volocity near volutetongue with time
3.3计算结果的时域分析
由于实际采集得到的压力信号值较大,不利于分析观察,因此本文将采集到的360个压力脉动信号进行去平均处理,即用各时间点压力脉动值减去采集时间内压力平均值[12]:
式中:pi为原始采集压力信号每一时间点压力脉动值;Pi为去平均处理后每一时间点压力脉动值;Pj为原始采集压力信号每一时间点压力脉动值.
如图11所示,在蜗壳周向表面,蜗舌附近的压力脉动最强烈,且B点处的压力脉动最强,这是由于叶轮流道流出的流体垂直冲击到蜗舌表面,使得此处的压力脉动最强,因此是气动噪声的主要来源.图12中叶片与蜗舌表面噪声源分布(脉动压力大小分布)也说明了这点.从图12中可以看出,蜗舌表面和靠近蜗舌附近的叶轮流道噪声源分布值最大,其他位置处噪声源值要小得多.这主要是由于叶轮和蜗舌的动静干涉,使得该流道附近存在漩涡流动,流体涡团周期性地从叶片尾端脱落,因此导致此处压力脉动增强,从而产生较大的噪声源.上述分析说明在扫路车用风机内部,正是由于叶轮与蜗舌非定常作用产生的漩涡,导致该区域的噪声源值最大,成为风机气动噪声的主要来源.
图11 蜗壳周向压力脉动情况Fig.11 Pressure fluctuation of thecircumferential volute
图12 叶片与蜗舌表面噪声源分布Fig.12 Noise source distribution ofimpeller and volute tongue
从图13可以看出,在1个周期内,叶片压力面上的脉动峰值有30个,这与叶片数量(15片长叶片+15片短叶片)对蜗舌表面的非定常作用相吻合,呈明显的周期性.而吸力面上压力脉动的周期性则要微弱得多,这主要是由于吸力面上容易产生流动分离,形成漩涡从而产生涡流噪声.由于叶片与蜗舌的非定常作用存在周期性,说明可以通过增大叶片与蜗舌的间距来削弱它们之间的非定常作用,由此降低扫路车用风机的离散噪声.
从图14可以看出,与叶片压力面上的脉动值不同,蜗舌表面的压力脉动并没有产生周期性变化.蜗舌表面的3个观察点中,B点的压力脉动值最大,A、C两点的压力脉动相对较弱.这主要是由于B点正对着叶轮流道气流的打击作用,因此该点压力脉动最强[13].
图13 叶片表面压力脉动随时间变化图Fig.13 Pressure fluctuation of the blade with time
图14 蜗舌表面压力脉动随时间变化图Fig.14 Pressure fluctuation of the volute tongue with time
(1)扫路车用风机靠近蜗舌区域的叶轮流道会产生强烈的漩涡流动乃至回流,该非定常流动增加了叶轮气流作用在蜗舌表面的非定常脉动作用力,使得该区域的噪声源值最大,成为扫路车用风机气动噪声的主要来源.
(2)扫路车用风机蜗舌结构形式及叶片与蜗舌间距对其气动噪声的产生具有重要影响,在设计时需要着重考量才能使得风机气动噪声性能达到最佳.
(3)叶片压力面比吸力面上的脉动周期性特征更强烈;正对叶轮流道的蜗舌区域压力脉动最强.
[1] 张斌,李亮,付玲,等.清扫车专用风机内流仿真分析与优化设计[J].机械设计与制造,2014,9(9):246-248.ZHANG Bin,LI Liang,FU Ling,et al.Analysis of theoptimized design and numerical simulation on the centrifugalfan of sweepers[J].Machinery Design & Manufacture,2014,9(9):246-248.
[2] Ansys Inc.Ansys fluent theory guide[M].Canonsburg:Ansys Inc.,2011.
[3] BALLESTEROS-TAJADURA R, VELARDE-SUARREZ.Numerical calculation of pressure fluctuations in the volute of acentrifugal fan[J].Journal of Fluids Engineering,2006,128:359-369.
[4] 张建.采用修正来流条件和粗糙壁面处理方法的绕流问题研究[D].北京:北京交通大学,2011.
ZHANG Jian.Analysis of flow around bluff body based on modified inlet profiles and rough wall treatments[D].Beijing:Beijing Jiaotong University,2011.
[5] 张斌.大涡模拟滤波网格分析及网格自适应控制研究与应用[D].上海:上海交通大学,2011.
ZHANG Bin.Research and application of filtering grid scale and meshing adaptive-control strategy for large eddy simulation[D].Shanghai:Shanghai Jiaotong University,2011.
[6] 李亮,王瑞,金晗辉.基于非线性湍流模型的多翼离心风机内部流场三维数值模拟及性能预测[J].流体机械,2014,40(10):41-45.
LI Liang,WANG Rui,JIN Hanhui.Three dimensional numerical simulation of internal flow field in multi-blade centrifugal fan based on the non-linear turbulent model[J].Fluid Machinery,2014,40(10):41-45.
[7] 付双成,刘雪冬,邹鑫,等.倾斜蜗舌对多翼离心风机流场及噪声的影响[J].噪声与振动控制,2013,33(3):87-89.
FU Shuangcheng,LIU Xuedong,ZOU Xin,et al.Numerical simulation of the influence of inclined volute tongue on flow field and noise of multi-blade centrifugal fans[J].Noise and Vibration Control,2013,33(3):87-89.
[8] 陈庆光,房菲,张永超,等.对旋风机级间流场涡量和噪声分布规律的分析[J].噪声与振动控制,2012,2(1):89-93.
CHEN Qingguang,FANG Fei,ZHANG Yongchao,et al.Vortex and noise characteristics analysis of the flow field between the two impellers of a contra rotating axial fan[J]. Noise and Vibration Control,2012,2(1):89-93.
[9]JAATINEN A,GRONMAN A,TURUNEN-SAARESTI T.Effect of vaneless diffuser wid th on the overall performance of a centrifugal compressor[J].Journal of Power and Energy,2011,225(5):665-673.
[10] 张强.气动声学基础[M].北京:国防工业出版社,2012.
ZHANG Qiang.Aero-acoustic foundation [M ].Beijing:National Defense Industrial Press,2012.
[11]ZHANG Bin, WANG Tong, GU Chuanggang. Blade optimization design and performance investigations of anultralow specific speed centrifugal blowers[J].Science inChina Series E:Technological Sciences,2011,54(1):203-210.
[12] 卜锋斌,赵博,代学昌,等.船用离心风机壁面周向压力脉动测试及分析[J].船舶工程,2013,35(6):40-43.
BU Fengbin,ZHAO Bo,DAI Xuechang,et al.Dynamic test and analysis of circumferential distribution of static pressureon marine fan volute wall[J].Ship Engineering,2013,35(6):40-43.
[13] 王军,冉苗苗,姚姗姗.多翼离心风机的气动噪声预测[J].工程热物理学报,2014,29(7):1144-1146.
WANG Jun,RAN Miaomiao,YAO Shanshan.Aero-acoustics prediction of multi-blade centrifugal fan[J].Journal of Engineering Thermophysics,2014,29(7):1144-1146.
Analysis on internal vortex and pressure pulsationproperties of road sweeper fan
LI Liang1,2,ZHANG Bin1,2,FU Ling1,2,XIAO Qing-lin1,TENG Xin-ke1,ZHANG Jin1,2(1.Zoomlion Heavy Industry Science &Technology Co.,Ltd.,Changsha 410013,China;2.State Key Laboratory of Construction Machinery,Changsha 410013,China)
Pertaining to the internal vortex and pressure pulsation properties of a specific road sweeperfan,a CFD numerical simulation method is first established.Then,the rationality and reliability of theproposed method are verified in comparison with simulation and testing results.Finally,it is indicatedthat,due to the strong vortex flow even recirculation near the volute tongue region,the fluidity is themajor source of fan aerodynamic noise.Owing that a low-speed vortex zone with small vorticity valuesexists near volute tongue outlet,the structural form of volute tongue and distance between blade andvolute tongue significantly impact aerodynamic noise.In addition,the periodic property of blade pressuresurface is stronger than that of blade suction surface,whereas the pressure pulsation of volute tongueagainst impeller flow passage is strongest.
acoustics;road sweeper fan;vortex;volute tongue;numerical simulation
TH 442
A
1672-5581(2015)06-0486-07
湖南省自然科学基金资助项目(13J J 3130)
李 亮(1985-),男,工程师.E-m a i l:20924010@z j u.e d u.c n
张 斌(1983-),男,高级工程师.E-m a i l:s j t u 2009@163.c Om