杜媛英,高 強(qiáng),尚長(zhǎng)春, 郝昱宇, 趙天鵬, 王曉璐
(1.西安科技大學(xué)工程訓(xùn)練中心,陜西 西安 710054; 2.上海電力修造總廠有限公司,上海 201316)
?
·機(jī)電工程·
基于大渦模擬的導(dǎo)葉式混流泵葉輪內(nèi)壓力脈動(dòng)特性分析
杜媛英1,高強(qiáng)2,尚長(zhǎng)春1, 郝昱宇1, 趙天鵬1, 王曉璐1
(1.西安科技大學(xué)工程訓(xùn)練中心,陜西 西安 710054; 2.上海電力修造總廠有限公司,上海 201316)
為研究導(dǎo)葉式混流泵葉輪內(nèi)部非定常壓力脈動(dòng)特性,在其葉輪進(jìn)口截面及出口截面附近分別設(shè)置8個(gè)壓力脈動(dòng)監(jiān)測(cè)點(diǎn),采用大渦模擬方法(LES)對(duì)導(dǎo)葉式混流泵內(nèi)全三維流道(進(jìn)水管、葉輪、導(dǎo)葉及出水管)進(jìn)行模擬,并對(duì)8個(gè)監(jiān)測(cè)點(diǎn)進(jìn)行壓力脈動(dòng)時(shí)域圖和頻域圖的分析。結(jié)果表明:由于旋轉(zhuǎn)葉輪旋轉(zhuǎn)失速、脫流效應(yīng)及靜止導(dǎo)葉的干涉作用,葉輪出口處壓力脈動(dòng)系數(shù)幅值均大于進(jìn)口處脈動(dòng)系數(shù)幅值,且其最大壓力脈動(dòng)發(fā)生在葉輪出口處,脈動(dòng)波衰減較慢;葉輪進(jìn)、出口截面上監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)頻率以葉輪葉頻為主頻次,且壓力脈動(dòng)主要頻率為葉頻的倍數(shù)。
混流泵;葉輪;壓力脈動(dòng);大渦模擬;時(shí)域圖;頻域圖
混流泵是一種應(yīng)用非常廣泛的通用流體機(jī)械,常用于農(nóng)田灌溉、防澇排洪、水利工程、污水處理、電站(燃煤電站、核電站、蓄能電站、潮汐電站)冷卻系統(tǒng)等各個(gè)領(lǐng)域[1]。壓力脈動(dòng)是影響泵運(yùn)行穩(wěn)定性的關(guān)鍵因素,混流泵在運(yùn)行時(shí),葉輪高速旋轉(zhuǎn),葉輪和靜止導(dǎo)葉之間的動(dòng)靜干涉及流體的高黏性,使其內(nèi)部流動(dòng)呈現(xiàn)出復(fù)雜的非定常性,引起泵內(nèi)流場(chǎng)的壓力脈動(dòng),從而導(dǎo)致整個(gè)機(jī)組或設(shè)備的振動(dòng)及噪聲。目前張德勝等利用SIMPLEC算法研究了葉輪與導(dǎo)葉數(shù)對(duì)斜流泵壓力脈動(dòng)的影響[2];王福軍等利用大渦模擬方法對(duì)大型軸流泵內(nèi)部非定常流動(dòng)特性進(jìn)行了模擬和計(jì)算[3];黎義斌等對(duì)低比轉(zhuǎn)數(shù)混流泵壓力脈動(dòng)特性做了數(shù)值模擬[4];黎耀軍等采用LES方法和滑移網(wǎng)格技術(shù)對(duì)一低比轉(zhuǎn)數(shù)軸流泵不同工況下輪緣間隙泄漏渦形態(tài)及葉片近輪緣區(qū)壓力脈動(dòng)特性做了分析,揭示了軸流泵葉頂間隙流動(dòng)特征[5];靳栓寶等認(rèn)為混流泵內(nèi)最大壓力脈動(dòng)在葉輪進(jìn)口,從輪轂到輪緣脈動(dòng)幅值遞增,在葉輪進(jìn)口和葉輪出口壓力脈動(dòng)頻率主要為葉輪葉頻控制[6];姬晉廷等對(duì)考慮水壓力的混流式轉(zhuǎn)輪的振動(dòng)特性進(jìn)行了分析[7];袁建平等基于大渦模擬對(duì)離心泵蝸殼內(nèi)壓力脈動(dòng)特性進(jìn)行了分析[8]。盡管運(yùn)用CFD計(jì)算泵內(nèi)流場(chǎng)已經(jīng)比較廣泛,但運(yùn)用CFD預(yù)測(cè)混流泵性能目前還主要集中在泵的外特性預(yù)測(cè)方面,同時(shí)大渦模擬是目前CFD領(lǐng)域的熱點(diǎn)之一;因此,采用大渦模擬方法對(duì)混流泵內(nèi)非定常流動(dòng)進(jìn)行的模擬,可以為混流泵進(jìn)行合理的結(jié)構(gòu)優(yōu)化、減輕振動(dòng)與噪聲提供理論依據(jù)。
本文利用Fluent軟件,采用大渦模擬方法對(duì)導(dǎo)葉式混流泵吸入室、葉輪、導(dǎo)葉及出水管組成的全流道做了三維非定常數(shù)值模擬,對(duì)葉輪進(jìn)口和出口壓力脈動(dòng)進(jìn)行了比較深入的分析。
1.1模型泵主要參數(shù)
計(jì)算模型為一比轉(zhuǎn)速為503的導(dǎo)葉式混流泵,其吸入室進(jìn)口直徑為224 mm,葉輪出口寬度為74 mm,導(dǎo)葉長(zhǎng)寬比為0.3,葉輪數(shù)為4,導(dǎo)葉數(shù)為7,泵的設(shè)計(jì)轉(zhuǎn)速為2 900 r/min,設(shè)計(jì)流量為1 400m3/s,設(shè)計(jì)揚(yáng)程為31 m。
1.2大渦模擬的基本思想及運(yùn)動(dòng)方程
大渦模擬是介于直接數(shù)值模擬(DNS)與Reynolds平均法之間的一種湍流數(shù)值模擬方法,它將比網(wǎng)格尺度大的湍流運(yùn)動(dòng)通過瞬時(shí)Navier-Stokes方程計(jì)算出來,而小尺度的渦對(duì)大尺度渦運(yùn)動(dòng)的影響通過一定的模型在針對(duì)大尺度渦的瞬時(shí)Navier-Stokes方程中體現(xiàn)出來,從而形成LES。
其運(yùn)動(dòng)方程組[9]為:
(1)
(2)
1.3計(jì)算區(qū)域及網(wǎng)格劃分
在Pro/E中繪制混流泵全流道模型(如圖1所示),網(wǎng)格生成通過Gambit完成,整個(gè)計(jì)算域采用結(jié)構(gòu)化和非結(jié)構(gòu)化網(wǎng)格組成的混合網(wǎng)格,計(jì)算區(qū)域?yàn)槲胧?、葉輪、導(dǎo)葉及出水管組成的全三維流道,網(wǎng)格計(jì)算總數(shù)為117萬6 945。
一是全面深化水利改革。深化水利審批制度改革,理清各級(jí)各部門的權(quán)利和責(zé)任清單,切實(shí)提高審批效能。完善水利投融資機(jī)制,吸引更多的社會(huì)資金投入水利建設(shè)。積極探索公益性水利項(xiàng)目的市場(chǎng)引入機(jī)制,探索通過政府采購服務(wù)方式解決水利工程管養(yǎng)難題。健全民生水利工程建設(shè)督導(dǎo)機(jī)制,切實(shí)加快建設(shè)進(jìn)度。
圖1 三維計(jì)算模型圖
1.4三維非定常計(jì)算方法
在對(duì)混流泵進(jìn)行非定常湍流數(shù)值計(jì)算時(shí),計(jì)算域中靜止導(dǎo)葉與高速旋轉(zhuǎn)葉輪采用滑移網(wǎng)格技術(shù),湍流模型采用LES湍流模型,非定常計(jì)算采用SIMPLE算法,采用非定長(zhǎng)計(jì)算的時(shí)間步長(zhǎng)為Δt=0.223 ms,即時(shí)間步長(zhǎng)為葉輪旋轉(zhuǎn)周期的1/90。計(jì)算中先對(duì)混流泵進(jìn)行三維定常湍流計(jì)算,將定常計(jì)算收斂后的結(jié)果作為非定常計(jì)算的初始值。
為了檢測(cè)高速旋轉(zhuǎn)的葉輪處壓力脈動(dòng),在葉輪進(jìn)口截面處設(shè)置P11~P14共4個(gè)監(jiān)測(cè)點(diǎn),在出口截面處設(shè)置P21~P24共4個(gè)監(jiān)測(cè)點(diǎn),如圖2所示。
圖2 監(jiān)測(cè)點(diǎn)分布圖
2.1葉輪進(jìn)口附近及出口附近靜壓分布圖
為了研究混流泵內(nèi)葉輪處壓力脈動(dòng)情況,采用LES方法對(duì)混流泵全三維流道做了非定常模擬,取葉片的靜壓分布圖及速度分布圖,如圖3所示。由圖3(a)可以看出,葉片進(jìn)口工作面存在低壓區(qū),此位置易發(fā)生汽蝕,靜壓沿葉高方向逐漸增大,壓力分布比較均勻,沿圓周方向呈現(xiàn)周期性分布,這主要是由于旋轉(zhuǎn)葉片對(duì)流體做功的原因。由圖3(b)可以看出在葉片進(jìn)口處速度較大,沿葉高方向,速度呈流線分布,無明顯的二次流、脫流等現(xiàn)象。
(a)葉輪工作面靜壓圖
(b)葉場(chǎng)圖3 葉輪工作面靜壓及速度分布圖
為了形象、直觀地表征混流泵各監(jiān)測(cè)點(diǎn)壓力脈動(dòng)值,消除靜壓的影響,本文引用量綱一的壓力系數(shù)cp,其表達(dá)式[10]如下:
(3)
在設(shè)計(jì)工況下,采用大渦模擬方法對(duì)混流泵過流通道進(jìn)行全三維非定常模擬,對(duì)8個(gè)監(jiān)測(cè)點(diǎn)進(jìn)行壓力脈動(dòng)時(shí)域特性分析。
圖4為葉輪進(jìn)、出口處一個(gè)周期內(nèi)的壓力脈動(dòng)時(shí)域圖,圖4(a)為最優(yōu)工況下葉輪進(jìn)口處沿徑向不同檢測(cè)點(diǎn)P11、P12、P13、P14處非定常壓力脈動(dòng)時(shí)域波形圖,圖4(b)為最優(yōu)工況下葉輪出口沿徑向不同檢測(cè)點(diǎn)P21、P22、P23、P24處非定常壓力脈動(dòng)時(shí)域波形圖。由圖4(a)可以看出:葉輪進(jìn)口處沿徑向不同檢測(cè)點(diǎn)P11、P12、P13、P14在一個(gè)周期內(nèi)壓力脈動(dòng)基本一致,均呈現(xiàn)有規(guī)律的正弦曲線,且均有4個(gè)波峰和4個(gè)波谷,這主要是由于4枚葉片周期性旋轉(zhuǎn)掠過造成流體質(zhì)點(diǎn)波動(dòng)。P13、P14處壓力脈動(dòng)系數(shù)較P11、P12處大,這主要是由于P11、P12點(diǎn)處于層流和紊流區(qū)域,流體質(zhì)點(diǎn)運(yùn)動(dòng)較為規(guī)律,而P13、P14處于紊流和湍流區(qū)域,流體質(zhì)點(diǎn)流動(dòng)狀態(tài)較為復(fù)雜。由圖4(b)可以看出最優(yōu)工況下葉輪出口沿徑向不同檢測(cè)點(diǎn)P21、P22、P23、P24壓力脈動(dòng)基本一致,均有4個(gè)大波峰和4個(gè)波谷,相比較葉輪進(jìn)口處而言,葉輪出口處壓力脈動(dòng)較為復(fù)雜。這主要是由于旋轉(zhuǎn)葉輪旋轉(zhuǎn)失速、脫流效應(yīng)及靜止導(dǎo)葉的干涉作用。將圖4(a)、(b)比較可以發(fā)現(xiàn)出口處壓力脈動(dòng)系數(shù)幅值均大于進(jìn)口處脈動(dòng)系數(shù)的幅值,因此可以認(rèn)為葉輪出口附近區(qū)域內(nèi)的流動(dòng)更為復(fù)雜。
(b)葉輪出口時(shí)域圖圖4 葉輪進(jìn)、出口時(shí)域圖
2.3壓力脈動(dòng)頻域特性分析
根據(jù)各監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)系數(shù),采用快速傅里葉變換(FFT),對(duì)8個(gè)監(jiān)測(cè)點(diǎn)處進(jìn)行壓力脈動(dòng)頻譜分析,如圖5所示:圖5(a)為最優(yōu)工況下葉輪進(jìn)口處沿徑向不同檢測(cè)點(diǎn)P11、P12、P13、P14處非定常壓力脈動(dòng)頻域圖;圖5(b)為最優(yōu)工況下葉輪出口沿徑向不同檢測(cè)點(diǎn)P21、P22、P23、P24處非定常壓力脈動(dòng)頻域圖。由圖5(a)、(b)可以看出,壓力脈動(dòng)的最大幅值均位于4倍轉(zhuǎn)頻處,同時(shí)由于葉片數(shù)為4,可發(fā)現(xiàn)壓力脈動(dòng)主要頻率為葉頻的倍數(shù)。將圖5(a)、(b)比較可以發(fā)現(xiàn),壓力脈動(dòng)幅值最大點(diǎn)位于葉輪出口圖上,且出口處壓力脈動(dòng)比進(jìn)口處復(fù)雜,脈動(dòng)波衰減慢,這與高速旋轉(zhuǎn)的導(dǎo)葉與靜止葉輪之間的動(dòng)靜干涉造成過度流道內(nèi)流動(dòng)較葉輪入口處復(fù)雜的實(shí)際情況相符。
(a)葉輪進(jìn)口頻域
(b)葉輪出口頻域圖5 葉輪進(jìn)、出口頻域圖
2.4泵外特性試驗(yàn)
不同流量下混流泵外特性計(jì)算值與實(shí)驗(yàn)值的比對(duì)分析算如圖6所示??梢钥闯觯涸?.7Q~1.2Q之間,混流泵性能預(yù)測(cè)值與實(shí)驗(yàn)值趨勢(shì)基本一致,由此可以認(rèn)為,大渦模擬可以較準(zhǔn)確地預(yù)測(cè)最優(yōu)工況及其附近工況的壓力脈動(dòng)結(jié)果的正確性。
圖6 泵性能曲線計(jì)算值與實(shí)驗(yàn)值對(duì)比圖
本文以比轉(zhuǎn)速為503的導(dǎo)葉式混流泵為研究對(duì)象,基于大渦模擬方法(LES)對(duì)該混流泵整個(gè)過流通道進(jìn)行了非定常模擬計(jì)算,在葉輪進(jìn)口截面處及葉輪出口截面處設(shè)置了8個(gè)壓力脈動(dòng)檢測(cè)點(diǎn),得出以下混流泵內(nèi)葉輪進(jìn)、出口處壓力脈動(dòng)特性。
1)由于旋轉(zhuǎn)葉輪旋轉(zhuǎn)失速、脫流效應(yīng)及靜止導(dǎo)葉的干涉作用,葉輪出口處壓力脈動(dòng)系數(shù)幅值均大于進(jìn)口處脈動(dòng)系數(shù)的幅值,混流泵葉輪內(nèi)的最大壓力脈動(dòng)發(fā)生在葉輪出口處,且脈動(dòng)波衰減較慢。
2)葉輪進(jìn)、出口截面上監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)以葉輪葉頻為主頻次,且壓力脈動(dòng)主要頻率為葉頻的倍數(shù)。
[1]丁成偉.離心泵與軸流泵原理及水力設(shè)計(jì)[M].北京:機(jī)械工業(yè)出版社,1981:6.
[2]張德勝,施衛(wèi)東,王川,等.斜流泵葉輪和導(dǎo)葉葉片數(shù)對(duì)壓力脈動(dòng)的影響[J].排灌機(jī)械工程學(xué)報(bào),2012,30(2):167.
[3]王福軍,張玲,張志民. 軸流泵不穩(wěn)定流場(chǎng)的壓力脈動(dòng)特性研究[J].水利學(xué)報(bào),2007,38(8):1003.
[4]黎義斌,李仁年,王秀勇,等.低比轉(zhuǎn)數(shù)混流泵壓力脈動(dòng)特性的數(shù)值模擬[J].排灌機(jī)械工程學(xué)報(bào),2013,31(3):205.
[5]黎耀軍,沈金峰,劉竹青,等.軸流泵輪緣間隙非定常流動(dòng)的大渦模擬[J].排灌機(jī)械工程學(xué)報(bào),2013,10(1):113.
[6]靳栓寶,王永生,常書平,等.混流泵內(nèi)流場(chǎng)壓力脈動(dòng)特性研究 [J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2013,44(3):64.
[7]姬晉廷,羅興锜,鄭小波.考慮水壓力的混流式轉(zhuǎn)輪振動(dòng)特性分析[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2008,24(3):257.
[8]袁建平,付燕霞,劉陽,等.基于大渦模擬的離心泵蝸殼內(nèi)壓力脈動(dòng)特性分析[J].排灌機(jī)械工程學(xué)報(bào),2010,28(4):310.
[9]王福軍.計(jì)算流體動(dòng)力學(xué)分析:CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社,2004:139-140.
[10]黎義斌,李仁年,王秀勇,等.混流泵內(nèi)不穩(wěn)定特性的數(shù)值模擬[J].排灌機(jī)械工程學(xué)報(bào),2013,31(5):384.
[11]關(guān)醒凡.軸流泵和斜流泵[M].北京:中國宇航出版社,2008:64.
(編校:夏書林)
Analysis of Pressure Fluctuations in Impeller of Guide Vane Mixed-flow Pump Based on Large Eddy Simulation
DU Yuanying1, GAO Qiang2SHANG Changchun1, HAO Yuyu1, ZHAO Tianpeng1, WANG Xiaolu1
(1.EngineeringTrainingCenter,Xi’anUniversityofScienceandTechnology,Xi’an710054China;2.ShanghaiPowerEquipmentManufactureCO.,LTD,Shanghai201316China.)
In order to investigate the phenomenon of pressure fluctuations caused by unsteady flow in impeller of guide vane mixed-flow pump, 8 pressure fluctuation test points were set on impeller inlet-plane and impeller outlet-plane respectively. LES was adopted to calculate the whole flow field which includes inlet pipe, impeller, guide vane and outlet pipe, and analyze the pressure fluctuation time domain graph and spectrum graph of 8 pressure fluctuation test points. The results shows thatdue to the stall andseparate flow of rotating impeller and interaction with static guide vane, the pressure fluctuation coefficient amplitude at outlet of impeller are greater than that at inlet of impeller, and the maximum amplitude pressure fluctuation occurs at outlet of impeller, and the impulse wave decays slowly. The main frequency of pressure fluctuations at impeller inlet-plane and outlet-plane is blade frequency and blade frequency is a multiple of impeller frequency.
mixed-flow pump; impeller; pressure fluctuation; LES; time domain graph; spectrum graph
2015-07-16
西安科技大學(xué)培育基金(201632)
杜媛英(1988—),女,在讀博士,工程師,主要研究方向?yàn)榱黧w機(jī)械的設(shè)計(jì)與水潤(rùn)滑橡膠軸承潤(rùn)滑研究。
TH313
A
1673-159X(2016)04-0075-4
10.3969/j.issn.1673-159X.2016.04.016