葉片出口角影響離心泵噪聲輻射數(shù)值研究

丁 劍， 劉厚林， 王 勇， 談明高， 崔建保

(江蘇大學(xué) 流體機械工程技術(shù)研究中心，江蘇 鎮(zhèn)江 212013)

離心泵噪聲包括機械振動噪聲與流體誘導(dǎo)噪聲[1]。機械振動噪聲主要源于機械設(shè)計、制造，此已通過主動控制技術(shù)解決；而離心泵內(nèi)部流體誘導(dǎo)噪聲機理更復(fù)雜，涉及頻率范圍更廣泛。水泵誘導(dǎo)噪聲已成研究熱點。

目前國內(nèi)外對流體誘導(dǎo)噪聲研究主要集中于試驗及數(shù)值計算。Choi等[2]采用熱線風(fēng)速儀及壓力傳感器測量離心泵葉輪速度、壓力脈動分量，發(fā)現(xiàn)葉輪出口存在射流/尾跡現(xiàn)象，并分析誘導(dǎo)噪聲機理；Srivastav等[3]通過測定振動速度、噪聲聲壓研究離心泵在不同工況下葉輪與隔舌間隙對噪聲影響，結(jié)果表明噪聲隨該間隙的增大而降低；袁壽其等[4-5]通過實驗研究離心泵內(nèi)壓力脈動與流動噪聲關(guān)系；劉厚林[6-7]基于實驗研究不同葉片數(shù)、(不同葉輪出口寬度對離心泵流體誘導(dǎo)噪聲影響。伴隨CFD技術(shù)的快速發(fā)展，CAA(計算氣動聲學(xué))技術(shù)不斷完善。工程計算聲學(xué)使用較廣的混合方法為：CFD+Kirchhoff方法、CFD+Lighthill聲比擬方法[8]。文獻[9-14]基于CFD結(jié)合聲類比方法研究泵內(nèi)水動力噪聲；Jiang等[15-17]通過LES結(jié)合聲振耦合求解方法對離心泵流動誘導(dǎo)噪聲進行數(shù)值模擬，結(jié)果表明泵殼結(jié)構(gòu)振動模態(tài)對外聲場計算結(jié)果影響顯著。

上述研究僅在數(shù)值計算方面對泵噪聲進行預(yù)測，缺少實驗驗證。離心泵內(nèi)非定常流動激勵力為泵流體誘導(dǎo)噪聲的主要聲源[18]，可用水聽器測量聲波在泵內(nèi)流體介質(zhì)中傳播形成的內(nèi)場噪聲；聲波向外界傳播，引起結(jié)構(gòu)振動，聲場-結(jié)構(gòu)振動協(xié)同作用引起空氣輻射噪聲，即外場噪聲。

本文借助文獻[15,19]，用大渦模擬方法計算離心泵內(nèi)部瞬態(tài)流場獲得蝸殼壁面偶極子聲源。在對泵殼體結(jié)構(gòu)模態(tài)分析基礎(chǔ)上，用LMS Virtual Lab間接邊界元IBEM聲振耦合模塊計算非定常流動引起的離心泵內(nèi)場噪聲，進行試驗驗證，并計算離心泵外場噪聲及聲輻射，研究葉片出口角對離心泵外場噪聲輻射影響。

1 泵結(jié)構(gòu)及實驗方案

選轉(zhuǎn)速ns=96.9的單級離心泵作為研究對象，見圖1(a)。葉輪為封閉式，采用三種不同葉片出口角度布置的葉輪，分別為葉輪29、33及37。泵設(shè)計轉(zhuǎn)速n=2 900 r/min，設(shè)計流量Q=50 m3/h，設(shè)計揚程H=30.18 m。離心泵結(jié)構(gòu)參數(shù)及測試工況見表1。

圖1 泵結(jié)構(gòu)及水體計算域示意圖

表1 離心泵主要結(jié)構(gòu)參數(shù)及測試工況

2 流場計算方法

2.1 控制方程

用CFX-12求解離心泵內(nèi)部三維流場，基于LES方法模擬湍流作用。用濾波函數(shù)對N-S方程進行濾波獲得LES控制方程。濾波后不可壓縮連續(xù)性方程及N-S方程為：

(1)

式中：Bi為哥氏力作用產(chǎn)生的源項；τij為亞格子尺度應(yīng)力，考慮小尺度渦作用，定義為：

(3)

通過渦粘度νsgs建立亞格子應(yīng)力τij與大尺度應(yīng)變率張量關(guān)系為：

(4)

用WALE(Wall-Adapting Local Eddy-viscosity)模型[20]求解渦粘度。該模型中含無量綱系數(shù)Cw，模型定義為：

(5)

(6)

式中：Δ為濾波長度，取決于網(wǎng)格體積(Δ=(ΔxΔyΔz)1/3)；Cw=0.1為模型系數(shù)。

2.2 計算網(wǎng)格及邊界條件

計算域分三部分，見圖1(b)。用多重坐標系，葉輪流場在旋轉(zhuǎn)坐標系中計算，蝸殼與泄漏流道在靜止坐標系中計算。在動靜部件間用交界面進行數(shù)據(jù)交換。穩(wěn)態(tài)計算用凍結(jié)轉(zhuǎn)子交界面(Fronzen rotor interface)；瞬態(tài)計算用瞬態(tài)動靜交界面(Transient rotor/stator interface)。計算域進口設(shè)為質(zhì)量流量邊界，出口表面靜壓設(shè)為0。計算域中所有表面均采用無滑移壁面條件。用GridPro 5.1生成高質(zhì)量結(jié)構(gòu)網(wǎng)格，見圖2。葉片上第一層節(jié)點距壁面距離定為1 μm，其它表面上該距離定為5 μm。計算獲得葉片表面y+值小于1，滿足LES對y+的要求。用三套不同密度網(wǎng)格研究大渦模擬網(wǎng)格無關(guān)性，三套網(wǎng)格定義為CM，MM，F(xiàn)M，見表2。

圖2 計算域網(wǎng)格

對三套網(wǎng)格進行定常RANS計算，獲得大尺度(L=k3/2/ε)、泰勒尺度(λ=(10νk/ε)1/2)及Kolmogorov尺度(η=(ν3/ε)1/4)在葉輪子午面(0≤mi/m≤1)及蝸殼圓周環(huán)面(0°≤α≤350°)的面積平均值，并與當(dāng)?shù)貫V波長度(Δ=V1/3)比較，見圖3、圖4。由二圖看出，MM，F(xiàn)M網(wǎng)格在大尺度與泰勒尺度之間，與泰勒尺度接近且位于同數(shù)量級；若繼續(xù)增加網(wǎng)格數(shù)量，即可使濾波長度小于泰勒尺度，但對計算機性能要求更高。考慮計算機能力，計算中采用FM網(wǎng)格，計算平臺為曙光TC3600刀片服務(wù)器，處理器為Intel Xeon CPU 5600。時間步長設(shè)為1.15×10-5s以保證平均庫朗特數(shù)小于5，該時間步長與葉輪旋轉(zhuǎn)一周的1 800步相當(dāng)。

表2 三套網(wǎng)格細節(jié)

圖3 葉輪網(wǎng)格空間分辨率

圖4 蝸殼網(wǎng)格空間分辨率

3 聲振耦合計算方法

邊界元BEM方法為基于波動方程的數(shù)值模擬方法，多用于計算結(jié)構(gòu)低模態(tài)密度頻域內(nèi)聲學(xué)問題。與有限元方法、差分法相比具有輸入數(shù)據(jù)少、計算時間短等優(yōu)點，適用于無限域及有限域。FEM-BEM耦合計算方法結(jié)合FEM法與BEM法各自優(yōu)勢，能計算聲學(xué)流體域與薄壁結(jié)構(gòu)間聲固耦合效應(yīng)。本文采用程序LMS Virtual.lab中間接邊界元IBEM聲振耦合模塊計算離心泵噪聲?；贚ES非定常計算輸出兩完整旋轉(zhuǎn)周期的蝸殼內(nèi)表面偶極子聲源后，對泵殼體結(jié)構(gòu)進行模態(tài)分析獲得結(jié)構(gòu)模態(tài)響應(yīng)，在LMS Virtual.Lab中進行聲振耦合求解，獲得聲壓及指向性信息。計算流程見圖5。

圖5 聲振耦合計算流程

采用有限元方法分析泵殼體模態(tài)響應(yīng)。泵殼體材料各向同性，彈性模量E=200 GPa，密度ρ=7 800 kg/m3，泊松比μ=0.3，阻尼系數(shù)比設(shè)置為0.01。泵體通過地腳螺栓與地基連接，模態(tài)分析邊界條件為與地基接觸底面施加3個位移約束。計算獲得前7階模態(tài)固有頻率見表3。

表3 泵殼體前7階固有頻率

4 離心泵內(nèi)聲場計算及實驗驗證

離心泵外場噪聲受電機噪聲、管道噪聲及背景噪聲干擾，精確測量有一定難度。為驗證數(shù)值計算方法的可信度，通過離心泵內(nèi)部噪聲信號進行驗證。離心泵內(nèi)部噪聲信號用水聽器測量，水聽器采用齊平式安裝方式，直接安裝在管壁上，使傳感器探頭與測壓點周圍壁面齊平，直接測量管內(nèi)流體噪聲，水聽器測點布置在模型泵出口4倍管徑處。實驗裝置見圖6。

由于蝸殼表面偶極子為泵內(nèi)主要噪聲源，葉輪轉(zhuǎn)動偶極子影響可忽略[14]，故本文采用大渦模擬僅提取蝸殼壁面偶極子聲源，基于FEM/BEM方法計算泵內(nèi)場噪聲，介質(zhì)為水，密度1 000 kg/m3，聲波傳播速度1 500 m/s，參考聲壓1×10-6Pa。

1.真空泵 2.汽蝕筒 3.4.8.10.碟閥 5.渦輪流量計 6.壓力變送器 7.模型泵 9.電機 11.穩(wěn)壓罐 12.13.球閥

圖7 各工況點計算與實驗結(jié)果對比

圖7為葉輪37模型泵不同流量的噪聲實驗與計算結(jié)果對比。由圖7可見，耦合計算結(jié)果與實驗值趨勢基本一致，耦合計算結(jié)果較非耦合情況更合理；實驗值與耦合計算值在葉片通過頻率BPF(241 Hz)處吻合良好。對定性、定量而言，基于FEMBEM聲振耦合方法計算BPF噪聲可行。

5 葉片出口角對外聲場影響

在實驗驗證對比基礎(chǔ)上，用FEMBEM方法計算離心泵外場噪聲輻射，介質(zhì)為空氣，密度1.225 kg/m3，聲波傳播速度340 m/s，在泵殼體外距幾何中心1 m處創(chuàng)建一球面場點網(wǎng)格，以觀察聲壓分布及計算透過的聲功率。表4為不同葉片出口角、不同工況在葉片通過頻率BPF處輻射聲功率。由表4看出，0.8Qopt及1.2Qopt工況下，隨葉片出口角β2的增加，離心泵在葉片通過頻率BPF處輻射聲功率增加；Qopt工況下，輻射聲功率隨出口角β2的增加先增大后減小。觀察整體趨勢，離心泵輻射聲功率隨葉片出口角β2的增大而增大。葉輪29、葉輪33模型泵輻射聲功率在Qopt處均遠大于0.8Qopt及1.2Qopt；葉輪37模型泵輻射聲功率在三工況下相差不大。各工況輻射聲功率數(shù)值比較顯示，葉輪29模型泵較優(yōu)。

表4 不同葉輪出口角離心泵BPF處輻射聲功率

圖8 外部監(jiān)測點布置

聲音傳播有明顯方向性，不同空間點相對聲源位置，方向不同噪聲頻譜存在差異。為獲得離心泵聲壓級周向分布，在葉輪中心回轉(zhuǎn)面的泵殼外部一圈距葉輪軸1 000 mm處設(shè)置36個監(jiān)測點，每個監(jiān)測點之間夾角為10°，見圖8。圖9為不同工況、不同葉片出口角在葉片通過頻率BPF處的外場輻射噪聲指向性分布。

由圖9(a)看出，0.8Qopt工況下葉輪29及葉輪33的噪聲指向分布曲線幾乎重合，而葉輪37的噪聲值明顯較大，該工況下聲壓級指向性分布極值位于0°、150°；由圖9(b)看出，Qopt工況下聲壓級隨葉片出口角增大而遞增，遞增幅度較小，聲壓級極值位于30°、210°；由圖9(c)看出，1.2Qopt工況下葉輪33、葉輪37指向分布曲線幾乎重合，葉輪29聲壓級最小，聲壓級極值位于30°、210°。對比圖9看出，相同流量下，不同葉片出口角度噪聲指向分布在各方向趨勢相同，極值位置一致；隨流量增加，其極大值方向在30°以內(nèi)變化，因此出口角度對聲壓級指向變化影響不大。而幅值，小流量工況下噪聲聲壓級受葉片出口角影響較明顯，葉片出口角由29°增至33°，噪聲聲壓級無明顯變化，由33°增至37°時，噪聲聲壓級增幅較大，較33°增大約28%；Qopt工況下各向聲壓級隨角度的增加均勻增大，遞增幅值約6%；大流量工況下，葉片出口角由33°增至37°，噪聲聲壓級增幅不明顯，由29°增至33°時，聲壓級增幅明顯，較29°增大約11% 。因此，葉片出口角增加超過一定值時，小流量工況下聲壓級增大顯著；而大流量下其聲壓級增加不明顯。表明葉片出口角存在合適范圍，保證泵性能在滿足使用要求的同時，使外場噪聲聲壓級在小流量工況下較小。

圖9 BPF處輻射噪聲指向性分布

6 結(jié) 論

基于LES結(jié)合FEMBEM聲振耦合計算方法分析葉片出口角對離心泵在水動力激勵下泵殼振動輻射噪聲影響，結(jié)論如下：

(1) 通過離心泵內(nèi)場噪聲實驗驗證對比，基于FEMBEM的聲振耦合計算方法可行，耦合計算較非耦合計算更合理。

(2) 離心泵葉片通過頻率BPF處的輻射聲功率隨葉片出口角β2的增大而增大。

(3) 外場噪聲聲壓級指向性分布顯示，葉片出口角存在合適范圍，得泵運行在小流量工況下時噪聲較小。該分析結(jié)果有利于拓寬離心泵設(shè)計思路，為新泵設(shè)計提供參考。

參 考 文 獻

[1]Langthjem M A，Olhoff N.A numerical study of flow-induced noise in a two-dimensional centrifugal pump[J].Part I: Hydrodynamics Journal of Fluids and Structures,2004, 19(3):349-368.

[2]Choi J S，McLaughlin D K，Thompson D E．Experiments on the unsteady flow field and noise generation in a centrifugal pump impeller[J]．Journal of Sound and Vibration,2003,263(3): 493-514.

[3]Srivastav O P，Pandu K R，Gupta K．Effect of radial gap between impeller and diffuser on vibration and noise in a centrifugal pump[C]//Journal of the Institution of Engineers (India)：Mechanical Engineering Division, 2003:36-39.

[4]Yuan S Q,Yang J,Yuan J P,et al.Experimental investigation on the flow-induced noise under variable conditions for centrifugal pumps[J]. Chinese Journal of Mechanical Engineering,2012,25(3):456-462.

[5]袁壽其,薛 菲,袁建平,等. 離心泵壓力脈動對流動噪聲影響的試驗研究[J]. 排灌機械,2009, 27(5): 287-290.

YUAN Shou-qi,XUE Fei,YUAN Jiang-ping,et al. Experimental study on impact of pressure fluctuation on flow-noise in centrifugal pump[J]. Drainage and Irrigation Machinery, 2009，27(5):287-290.

[6]劉厚林,王 勇,袁壽其,等.葉片出口寬度對離心泵流動誘導(dǎo)振動噪聲的影響[J].華中科技大學(xué)學(xué)報(自然科學(xué)版)，2012，40(1):123-127.

LIU Hou-lin，WANG Yong,YUAN Shou-qi,et al.Effects of impeller outlet width on the vibration and noise from centrifugal pumps induced by flow[J].J. Huazhong Univ. of Sci. Tech.(Natural Science Edition)，2012,40(1):123-127.

[7]談明高,王 勇,劉厚林,等.葉片數(shù)對離心泵內(nèi)流誘導(dǎo)振動噪聲的影響[J].排灌機械工程學(xué)報,2012,30(2):131-135.

TAN Ming-gao,WANG Yong,LIU Hou-lin,et al. Effects of number of blades on flow induced vibration and noise of centrifugal pumps[J].Journal of Drainage and Irrigation Machinery Engineering,2012,30(2):131-135.

[8]Wang M，F(xiàn)reund J B，Lele S K.Computational prediction of flow-generated sound[J]. Annual Review of Fluid Mechanics，2005,38(1): 483-512.

[9]Kato C,Yamade Y,Wang H,et al.Numerical prediction of sound generated from flows with a low Mach number[J].Computers & Fluids,2007,36(1):53-68.

[10]黃俊雄,耿少娟,吳 瑞,等.不同葉輪形式下離心泵噪聲特性對比研究[J].聲學(xué)學(xué)報,2010,35(2):113-118.

HUANG Jun-xiong,GENG Shao-juan,WU Rui,et al. Comparion of noise characteristics in centrifugal pumps with different types of impellers[J].Acta Acustica,2010, 35(2):113-118.

[11]袁壽其,司喬瑞,薛 菲,等. 離心泵蝸殼內(nèi)部流動誘導(dǎo)噪聲的數(shù)值計算[J]. 排灌機械工程學(xué)報,2011, 29(2): 93-98.

YUAN Shou-qi,SI Qiao-rui,Xue Fei,et al. Numerical calculation of internal flow-induced noise in centrifugal pump volute[J]. Journal of Drainage and Irrigation Machinery Engineering, 2011,29(2):93-98.

[12]王宏光,徐小龍,楊愛玲,等. 軸流泵流動噪聲數(shù)值模擬[J]. 排灌機械工程學(xué)報,2011, 29(3):199-203.

WANG Hong-guang,XU Xiao-long,YANG Ai-ling,et al. Numerical simulation of flow noise in axial-flow pump[J]. Journal of Drainage and Irrigation Machinery Engineering, 2011,29(3):199-203.

[13]袁建平,金 榮,陳紅亮,等. 離心泵用赫姆霍茲水消聲器聲學(xué)特性數(shù)值模擬[J]. 排灌機械工程學(xué)報,2012,30(2):141-146.

YUAN Jian-ping,JIN Rong,CHEN Hong-liang,et al. Numerical simulation of acoustic performance of Helmholtz water-muffler for centrifugal pumps[J]. Journal of Drainage and Irrigation Machinery Engineering, 2012,30(2):141-146.

[14]何 濤,尹志勇,孫玉東.離心泵流動誘發(fā)振動特性數(shù)值計算分析[J].振動與沖擊,2012,31(12):96-102.

HE Tao,YIN Zhi-yong,SUN Yu-dong.Numerical analysis for flow induced vibration of a centrifugal pump[J].Journal of Vibration and Shock,2012,31(12):96-102.

[15]Jiang Y Y,Yoshimura S,Imai R,et al．Quantitative evaluation of flow-induced structural vibration and noise in turbomachinery by full-scale weakly coupled simulation[J]．Journal of Fluids and Structures,2007,23(3):531-544.

[16]譚永學(xué),王宏光,楊愛玲,等. 離心泵水動力噪聲預(yù)測[J].上海理工大學(xué)學(xué)報,2011,33(1): 89-94.

TAN Yong-xue,WANG Hong-guang,YANG Ai-ling,et al. Numerical prediction of hydrodynamic noise for a centrifugal pump[J]. Journal University of Shanghai for Science and Technology,2011,33(1):89-94.

[17]王鎮(zhèn)宇,楊愛玲,戴 韌. 離心泵流動誘導(dǎo)噪聲的數(shù)值預(yù)測[J].機械工程學(xué)報,2012,48(6):162-167.

WANG Zhen-yu,YANG Ai-ling,DAI Ren.Numerical prediction of the flow-induced noise of centrifugal pump[J]. Journal of Mechanical Engineering,2012,48(6):162-167.

[18]蔣愛華,章 藝,靳思宇,等.離心泵葉輪流體激勵力的研究[J].振動與沖擊,2012,31(22):123-127.

JIANG Ai-hua,ZHANG Yi,JIN Si-yu,et al.Fluid exciting forces on centrifugual pump part Ⅰ: force on volute[J].Journal of Vibration and Shock,2012,31(22):123-127.

[19]Wagner C,Hüttl T,Sagaut P.Large-eddy simulation for acoustics[M].NewYork：Cambridge University Press,2007.

[20]Nicoud F,Ducros F.Subgrid-scale stress modeling based on the square of velocity gradient tensor[J].Flow Turbulence Combustion,1999,62(3):183-200.