高波,周志威,倪丹,張寧,顧嘉嶸
(江蘇大學能源與動力工程學院, 江蘇 鎮(zhèn)江 212013)
離心泵是指靠葉輪旋轉時產(chǎn)生的離心力來輸送液體的泵,其廣泛應用于電力、冶金、煤炭、建材等輸送含有固體顆粒的漿體.離心泵內動靜干涉是產(chǎn)生壓力脈動的主要激勵源,而葉片尾跡是影響動靜干涉現(xiàn)象的重要因素,特別對葉頻處壓力脈動能量影響顯著[1-2].葉片尾跡與隔舌產(chǎn)生干涉作用后,葉片尾緣的脫落渦撞擊隔舌后產(chǎn)生的切割和變形可誘發(fā)強烈的流體激勵力,在該區(qū)域造成較大的壓力脈動幅值[3-5].通過改變葉片葉型,可有效降低葉片尾緣的脫落渦強度,改善流動分離現(xiàn)象,提高整泵水力效率[6-8],同時也可以顯著改善葉片出口區(qū)域的流動均勻性,進而降低壓力脈動強度[9].因此開展離心泵葉輪葉片改型設計,改善葉輪尾跡結構,進而優(yōu)化離心泵內部干涉流場,是降低動靜干涉作用影響的重要途徑,也可為低振動噪聲離心泵水力設計提供新方法[10-11].
翼型葉片對優(yōu)化葉片尾跡結構與降低葉片脫落渦強度效果顯著[12],為此提出了一種用于離心泵葉輪的翼型葉片優(yōu)化結構.另外,錯列葉柵可改變單翼型的尾跡結構[13],故文中在翼型葉片的基礎上增設偏置小翼,形成錯列葉柵,改善葉輪出口尾跡.文中以一臺低比轉數(shù)離心泵為研究對象,基于DDES方法,開展普通葉輪、翼型葉片葉輪和偏置小翼葉輪這3種方案的數(shù)值計算,對比分析3種方案下泵壓力脈動特性與尾跡結構,揭示偏置小翼對泵壓力脈動和尾跡結構的影響規(guī)律,以期為高效低噪聲離心泵設計提供參考.
模型泵設計流量Qd=55 m3/h,設計揚程Hd=20 m,轉速n=1 450 r/min,具體參數(shù)可見文獻[5].普通葉輪方案如圖1a所示.
圖1 葉輪設計方案
翼型葉片葉輪與偏置小翼葉輪在該普通葉輪基礎上修改,具體修改如下:將常規(guī)厚度變化的圓柱葉片改為NACA4418翼型葉片,如圖1b所示;為改善葉輪尾跡,優(yōu)化葉輪出口流動狀態(tài),在翼型葉片基礎上增設偏置小翼,翼型與長葉片相同,如圖1c所示.
圖2為離心泵計算域示意圖,其中葉輪為旋轉域,進口管、蝸殼、擴散段為靜止域.采用結構化網(wǎng)格對泵全流域進行劃分,為了保證結果的精確性,對邊界層進行加密處理,如圖3所示.經(jīng)過網(wǎng)格無關性驗證發(fā)現(xiàn),當網(wǎng)格數(shù)量大于1.2×107時,邊界層壁面y+值變化較小.為保證計算的準確性,最終選取總網(wǎng)格數(shù)為12 729 370.
圖2 離心泵計算域
圖3 離心泵網(wǎng)格
為探究葉片尾跡對壓力脈動的影響,在蝸殼流道內沿旋轉軸均勻布置20個監(jiān)測點,如圖4所示.采用DDES對泵設計工況下的內流場進行非定常計算[5],可以阻止模型應力損耗以及網(wǎng)格導致的分離,能夠精細地捕捉葉片尾跡,且對邊界層網(wǎng)格要求低.采用SIMPLEC算法以及二階迎風格式,葉輪為旋轉域,其他部件為固定域,采用Interface對部件結合面進行處理,其余壁面Wall采用標準壁面函數(shù),殘差精度設為1×10-6.進口采用速度邊界條件、出口采用壓力出口條件,壓力設置為101.325 kPa.將基于SSTk-ω湍流模型的定常計算結果作為非定常計算的初始條件,為保證數(shù)值模擬的計算精度,葉輪流道內庫朗數(shù)約為0.36,其他流域庫朗數(shù)為0.69,最終確定非定常模擬計算步長△t=0.000 114 9 s,即離心泵葉輪旋轉1°所需要的時間.
圖4 壓力脈動監(jiān)測點示意圖
首先對模型泵的流量和揚程進行了變換,用下式進行量綱一化
(1)
(2)
式中:u2為葉片出口的切向速度;R2為葉輪的外半徑;g為重力加速度.
研究發(fā)現(xiàn)普通葉輪方案的SSTk-ω,DDES數(shù)值計算結果和試驗結果的性能曲線基本吻合,在設計點揚程系數(shù)誤差小于0.025[5].圖5為不同計算方法下的速度云圖,可以看出,計算方法流動結構是相似的,但是DDES能夠更加精確地捕捉到內部的流動細節(jié),可認為該方法更準確可靠.
圖5 不同計算方法下的速度云圖
定義偏置角φ為長葉片尾緣至葉輪旋轉中心與短葉片尾緣至葉輪旋轉中心的夾角,L=R0/R2為偏置小翼葉片長度,其中R2為葉輪外半徑,R0為葉片截斷點到葉輪中心的距離,將偏置截斷后的圓柱葉片按翼型葉片厚度分布規(guī)律設計為翼型葉片.對不同偏置角度下的葉輪進行性能預測可知(如圖6所示):偏置小翼葉輪做功能力增強,揚程H均有所增加,當小翼置于長葉片工作面?zhèn)葧r(負角度),效率ηh顯著下降,而置于背面?zhèn)葧r(正角度),效率幾乎不變;對不同R0/R2偏置小翼進行性能預測發(fā)現(xiàn),當短葉片長度為0.85時揚程達到最高點.綜合以上因素,偏置小翼方案選取偏置角φ為10°,L=R0/R2=0.85.
圖6 設計工況不同偏置角及不同長度下小翼的離心泵外特性曲線
圖7為不同方案泵外特性對比,其中Φ為流量系數(shù),圖中所示的OR為普通葉輪方案,AF為翼型葉片方案,DAF為偏置小翼方案.
由圖7可知:與普通葉輪方案相比,翼型葉片方案的揚程在小流量工況較低,在設計工況以及大流量工況相差較小.偏置小翼方案揚程明顯高于其他2個方案,在設計工況下,揚程增加12%,在設計工況以及大流量工況效率提升,效率最高點較普通葉輪方案增加3.06%,設計工況及大流量工況效率明顯提高,高效工作區(qū)變寬.
圖7 不同方案離心泵外特性曲線
通過非定常計算提取監(jiān)測點的壓力脈動信號,選取葉片-隔舌動靜干涉作用較強的隔舌附近監(jiān)測點的壓力信號,經(jīng)過FFT變換得到設計工況下壓力脈動頻譜圖,如圖8所示.3種方案均捕捉到明顯的葉頻(145 Hz)及其倍頻.在普通葉輪方案下,可以捕捉到明顯的低頻雜亂信號;在隔舌附近的監(jiān)測點P2有大量低頻雜亂信號.改用翼型葉片方案后該現(xiàn)象有明顯改善,大量低頻雜亂信號得到了一定程度的抑制,且捕捉到了1.5倍軸頻(37 Hz),可能為周期性的脫落渦導致.偏置小翼方案低頻信號進一步被抑制,1.5倍軸頻消失,葉頻幅值明顯下降,在隔舌處葉頻幅值最高下降43%,由此可知,偏置小翼對降低離心泵壓力脈動、抑制低頻雜亂信號有明顯效果.
圖8 設計工況壓力脈動頻譜圖
圖9為壓力脈動葉頻幅值周向分布,在隔舌附近干涉作用最強,沿葉輪旋轉方向,動靜干涉作用逐漸減弱.在隔舌附近普通葉輪方案葉頻幅值最高,翼型葉片方案次之,偏置小翼方案最低.在隔舌下游處,翼型葉片方案與偏置小翼方案葉頻幅值相差較小.
圖9 壓力脈動葉頻幅值分布
圖10為蝸殼監(jiān)測點0~fBPF頻段壓力脈動分布,普通葉輪方案在低頻處存在大量低頻雜亂信號,隔舌處脈動能量最高,翼型葉片方案及偏置小翼方案壓力脈動能量明顯降低,翼型葉片方案存在能量較高的1.5倍軸頻,在其他低頻處脈動能量較低,因此總體能量值低于偏置小翼方案.
圖10 0~fBPF頻段壓力脈動幅值分布
圖11為尾跡渦量(Q準則)時空分布圖,圖中t0為初始時刻;Δt為時間步長.普通葉輪方案在葉片尾緣工作面存在大量高渦量區(qū),當葉輪掃掠隔舌時該高渦量區(qū)從葉片尾緣脫落,形成大面積高渦量尾跡區(qū),撞擊隔舌后呈現(xiàn)團狀、條帶狀混合分布,高渦量區(qū)集中且密集,并向下游移動,該尾跡周期性撞擊隔舌,這是引起離心泵壓力脈動的主因[1];翼型葉片方案葉片尾緣的脫落渦多呈現(xiàn)為條帶狀,葉片工作面尾緣的高渦量區(qū)消失,蝸殼流道中大面積的團狀高渦量區(qū)被抑制;偏置小翼方案中小翼吸力面存在明顯尾跡區(qū),該尾跡區(qū)與長葉片尾跡相互干涉,使尾跡能量耗散加快,高渦量尾跡區(qū)明顯減少,尾跡多為零星條帶狀高渦量區(qū),蝸殼流道中高渦量區(qū)明顯減少.由此可知,偏置小翼對改善葉片尾跡、優(yōu)化葉輪出口流動具有積極作用.
圖11 葉輪尾跡Q準則
三維渦結構可更加直觀地展現(xiàn)流道中渦的分布規(guī)律,如圖12所示,等值面采用Q準則,用壓力著色.由圖可知:普通葉輪方案蝸殼流道內存在大量高渦量區(qū),翼型葉片及偏置小翼方案相較模型泵方案團狀高渦量區(qū)明顯減少,多為條帶狀高渦量區(qū).通過高渦量區(qū)對比,說明壓力脈動低頻雜亂的信號多為葉輪尾跡中密集、碎裂的團狀高渦量區(qū)撞擊隔舌導致,葉頻能量值則與條帶狀高渦量區(qū)撞擊隔舌相關.翼型葉片和偏置小翼對離心泵葉片的尾跡渦具有明顯的抑制作用.
圖12 葉輪尾跡三維渦量分布
以低比轉數(shù)離心泵為研究對象,采用DDES方法,開展了泵內精細非定常數(shù)值計算,對比了普通葉輪、翼型葉片葉輪和偏置小翼葉輪3種方案下的泵壓力脈動頻譜特性及尾跡渦結構,獲得的結論如下:
1) 偏置小翼改善了離心泵能量性能,與普通葉輪相比,設計工況揚程提高12%,高效點效率提高3.06%,且高效工作區(qū)變寬.
2) 偏置小翼可明顯抑制低頻雜亂信號,有效降低壓力脈動葉頻幅值,在隔舌處葉頻幅值最高降低43%,改善了因壓力脈動造成的離心泵振動噪聲.
3) 翼型葉片及偏置小翼能改善葉輪尾跡,抑制大量團狀高渦量區(qū)的形成,這是改善低頻雜亂激勵信號的重要因素,尾跡中條帶狀及團狀渦量區(qū)明顯減少,撞擊隔舌的能量明顯下降,這是壓力脈動葉頻幅值降低的主因.