□耿鑫琪

（武漢中糧食品科技有限公司 湖北 武漢 430000）

離心泵小流量工況下的內(nèi)部流動特性

□耿鑫琪

（武漢中糧食品科技有限公司 湖北 武漢 430000）

隨著科學技術(shù)的發(fā)展和離心泵使用量的不斷增加，人們越來越重視離心泵運行的穩(wěn)定性和高效性，離心泵內(nèi)部流動特性決定著其性能和運行狀態(tài)，一般說來，離心泵內(nèi)部流動性越好，速度和壓力的分布更具規(guī)律性，則運行狀態(tài)更穩(wěn)定。不同的離心泵系統(tǒng)設(shè)置不同，且離心泵運行狀態(tài)也不同，多是在小流量工況下運行，一般情況下，離心泵在小流量工況下運行會出現(xiàn)不穩(wěn)定現(xiàn)象，比如會出現(xiàn)回流、失速、流動分離等情況的發(fā)生，這時會嚴重影響離心泵運行的穩(wěn)定性和高效性，因此，為更好地研究離心泵小流量工況下的內(nèi)部流動特性，文中針對IS160-50-65型離心泵在小流量工況下進行數(shù)值模擬分析，并為離心泵小流量工況下的內(nèi)部流動特性提供一定的數(shù)值參考。

離心泵；小流量工況；內(nèi)部流動特性

泵是現(xiàn)代工業(yè)中能量轉(zhuǎn)換和液體輸送的重要動力裝置，而離心泵是應(yīng)用最為廣泛的裝置，被廣泛的應(yīng)用在能源行業(yè)、農(nóng)業(yè)灌溉以及交通運輸行業(yè)等多個領(lǐng)域，通過對離心泵的研究可以發(fā)現(xiàn)，離心泵的運行狀態(tài)取決于其內(nèi)部流動特性。離心泵大部分情況下是在小流量工況下運行，但是，其在小流量工況下運行容易產(chǎn)生回流現(xiàn)象，且葉輪流道內(nèi)可能會發(fā)生分離漩渦，從而造成離心泵系統(tǒng)運行的不穩(wěn)定性，因此，本文結(jié)合數(shù)值模擬方法對小流量工況下的內(nèi)部流動性進行分析，得出流動規(guī)律，從而為離心泵運行的穩(wěn)定性提供支撐。

1 計算模型與計算方法

1.1 計算模型

本文的研究對象是低速率離心泵，其設(shè)計參數(shù)見表1，采用Pro/E軟件進行三維建模，且為了最大程度的排除干擾，在葉輪進口增加了長度為五倍的進口直徑，在出口增加了長度為五倍的出口延伸段。

表1 離心泵主要設(shè)計參數(shù)

1.2 網(wǎng)格劃分

模型建立完成后要進行數(shù)值計算，而網(wǎng)格劃分是數(shù)值計算的前提，網(wǎng)格劃分質(zhì)量影響著數(shù)值計算的準確性，本文進行離心泵內(nèi)部流動特性的研究采用ICEM對網(wǎng)格進行劃分，共劃分為四個部分，進口管道、葉輪流道、蝸殼流道和擴壓室，考慮到離心泵設(shè)計復雜且工況不定，進行建模時則采用性能比較好的四面體非結(jié)構(gòu)網(wǎng)絡(luò)。

1.3 網(wǎng)格無關(guān)性驗證

網(wǎng)格的劃分對模型數(shù)值的準確性影響很大，因此，對離心泵研究建模之前需要對網(wǎng)格做無關(guān)性驗證，便于選取合理的模型結(jié)果。網(wǎng)格無關(guān)性驗證指的是當網(wǎng)格數(shù)量有一個恒定值時，數(shù)值計算結(jié)果不再隨數(shù)量的增加而改變數(shù)值，計算結(jié)果也會得出一個恒定值。網(wǎng)格無關(guān)性驗證是數(shù)值計算基礎(chǔ)，也是保證模型穩(wěn)定性的前提，其中，揚程可當做參考標準，即網(wǎng)格數(shù)量達到恒定值，揚程變化比較小的情況下，可以驗證網(wǎng)格數(shù)量對數(shù)值驗證結(jié)果是可信的，且數(shù)值的準確度與網(wǎng)格數(shù)量無關(guān)?；诖耍疚倪x取一定范圍的網(wǎng)格數(shù)，分成八組網(wǎng)格，進行數(shù)值分析，見表2。

通過表2可以得出，網(wǎng)格數(shù)達到108萬時，揚程最高，而當網(wǎng)格數(shù)繼續(xù)增大達到200萬時，揚程則降低到20.657，網(wǎng)格數(shù)繼續(xù)增大達到220萬時，揚程數(shù)為20.652，變化不大，而且隨著網(wǎng)格數(shù)的增大，揚程數(shù)基本沒有變化，這說明，網(wǎng)格無關(guān)性驗證指的是當網(wǎng)格數(shù)量有一個恒定值時，數(shù)值計算結(jié)果不再隨數(shù)量的增加而改變數(shù)值，計算結(jié)果也會得出一個恒定值。

2 離心泵內(nèi)部流場定常數(shù)值研究

離心泵內(nèi)部流動特性會隨著流量的變化不斷地發(fā)生改變，本文主要針對離心泵小流量工況下的內(nèi)部流動流動特性進行定常數(shù)值模型研究。

2.1 離心泵外特性預(yù)測分析

離心泵的揚程和速率是內(nèi)部流動特性的表現(xiàn)形式，通過對1.25m3/h到16.25m3/h中多個工況進行數(shù)值計算，得出離心泵內(nèi)部流動特性曲線。

揚程計算公式：H=P2-P1/Pg，其中，P2為離心泵出口總壓，P1為進口總壓，單位為Pa；離心泵速率計算公式：ηh=pgQH/Mω,其中，M為葉輪扭矩，單位為N·s,ω為葉輪角速度，單位為rda/s。根據(jù)揚程、速率的計算公式，可以得出小流量工況下的各種數(shù)據(jù)，進而得出離心泵小流量工況下的內(nèi)部流動特性。根據(jù)不同流量下離心泵的數(shù)值結(jié)果，可以得出小流量工況下離心泵的揚程和速率的計算數(shù)值，進而通過模擬計算，得出不同工況下的流動特性，但是，考慮到設(shè)備條件的不同，對流動性能的模擬可能會存在誤差，將非全流場計算得出的揚程值的5%作為損失補充部分。

表2 離心泵不同網(wǎng)格數(shù)的預(yù)測揚程

2.2 離心泵內(nèi)部流場計算分析

為更好地測算離心泵小流量工況下的內(nèi)部流動特性，特分析流量工況Q/Qd=1.0和小流量工況Q/Qd=0.4，Q/Qd=0.2時內(nèi)部流動情況的數(shù)值分析，通過不同流量下的方案設(shè)計，通過計算可以得出，葉輪軸面分布均勻的情況下，離心泵的內(nèi)部流量穩(wěn)定；小流量工況下，離心泵出現(xiàn)回流情況，而且隨著流量越來越小，回流現(xiàn)象會變得越來越嚴重，而且根據(jù)分析可以的得出，小流量Q/Qd=0.2的工況情況下，離心泵流道中的漩渦是最多的，也就是說回流最為明顯，因此，可以得出結(jié)論，在液體回流的情況下，流體逆主流方向回流至進口管的離心泵進口管道中，甚至可能會回流到吸水管中。

3 離心泵內(nèi)部流場非定常數(shù)值研究

上一部分通過定常數(shù)值分析了離心泵不同流量工況下的數(shù)值分析，并獲得了內(nèi)部流動特性，但是離心泵的內(nèi)部流動特性是一個復雜的三維模型，僅憑定常數(shù)值分析無法真實反映離心泵小流量工況下的內(nèi)部流動特性，因此，還需要對離心泵內(nèi)部流場非定常數(shù)值進行分析。

為了對蝸殼中的速度進行分析，選取低比轉(zhuǎn)速泵，通過徑向尺寸的對比，軸向速度基本無變化，通過分析可以得出結(jié)論，隨著離心泵流量的降低，周向速度和徑向速度也是不斷地降低的，使得進入蝸殼流道中的液體不斷地減少，進而使得蝸殼內(nèi)的速度較低。本部分采用非定常數(shù)值分析法對離心泵小流量工況下的內(nèi)部流動性進行分析，通過能夠監(jiān)測點的設(shè)置分析了速度和壓力，并通過快速傅里葉變換獲得壓力脈動的頻率特性。通過模型的建立可以得出，小流量工況下，液體出現(xiàn)分離和回流情況，且高能流體會對進口管道內(nèi)的主流產(chǎn)生干擾，使得進口管道的流體變得沒有規(guī)律性，甚至可能會發(fā)生漩渦現(xiàn)象，而且隨著流量的不斷減小，離心泵內(nèi)部的流線也是越來越?jīng)]有規(guī)律，進而占據(jù)了整個空間，極易引起內(nèi)部流量的不穩(wěn)定性，情況嚴重時可能會發(fā)生內(nèi)部斷流現(xiàn)象。

結(jié)束語

通過分析可以得出結(jié)論：不同流量的設(shè)計工況下，離心泵揚程預(yù)測偏差為1.47%時，則偏差為3.61%；在小流量工況下，流量越小，離心泵的預(yù)測值與實際值偏差趨勢越小，而效率則不斷地增大；離心泵在小流量工況下，進口管道會出現(xiàn)回流現(xiàn)象，且隨著流量的減小，回流現(xiàn)象越明顯，甚至可能會出現(xiàn)堵塞葉輪流道的情況。

[1]黃建德.開式和閉式離心泵進口回流的研究[J].工程熱物理學報,1997年01期.

[2]祝磊,袁壽其,袁建平,裴吉.階梯隔舌對離心泵壓力脈動和徑向力影響的數(shù)值模擬[J].農(nóng)業(yè)機械學報,2010年S1期.

[3]張金鳳,梁赟,袁建平,袁壽其.離心泵進口回流流場及其控制方法的數(shù)值模擬[J].江蘇大學學報(自然科學版),2012年04期.

[4]李曉俊,袁壽其,潘中永,李亞林,陽君.離心泵邊界層網(wǎng)格的實現(xiàn)及應(yīng)用評價[J].農(nóng)業(yè)工程學報,2012年20期.

1004-7026（2016）18-0089-02

S277.9;TH311

A

DOI：10.16675/j.cnki.cn14-1065/f.2016.18.069