大型泵站箱形進水流道內(nèi)部流動數(shù)值模擬

周濟人,金 燕,劉 超

(1.河海大學(xué)水利水電工程學(xué)院,江蘇 南京 210098;2.揚州大學(xué)水利科學(xué)與工程學(xué)院,江蘇 揚州 225009)

0 引 言

泵站進水流道是前池與水泵葉輪室之間過渡段,其作用是使水流由前池流向葉輪室的過程中能更好地轉(zhuǎn)向和收縮,使水泵葉輪進口具有較好的流速分布和壓力分布,進水流道內(nèi)流態(tài)的優(yōu)劣直接影響水泵的性能及安全運行[1]。某大型低揚程泵站為滿足自引、自排、抽引及抽排要求,泵站進出水流道設(shè)計為箱形雙向形式。目前對雙向流道泵裝置的研究多為模型試驗[2-6],對其內(nèi)部流場了解較少。本文采用CFD數(shù)值模擬方法[7-9]能夠獲取進水流道內(nèi)的流場,可以詳細(xì)了解流速分布情況,尤其是泵進口的流速分布,對改善進口水流條件,保證水泵性能的發(fā)揮和安全運行具有重要的意義。

1 計算模型

本文采用FLUENT軟件對大型低揚程泵站雙向箱形進水流道內(nèi)部三維紊流流動進行數(shù)值模擬。

1.1 計算區(qū)域

數(shù)值模擬以進水流道為計算對象,前池和進水流道外形輪廓為計算邊界,取整個進水流道包圍的水體為計算區(qū)域,如圖1所示。為了模擬水泵葉輪進口條件,進水流道出口設(shè)置了圓柱輪轂體。應(yīng)用有限體積法離散和SIMPLE算法計算求解速度場和壓力場。

1.2 控制方程

進水流道內(nèi)水流運動作為一種收縮型管道運動,只要設(shè)計得當(dāng)不易出現(xiàn)回流現(xiàn)象,因此采用標(biāo)準(zhǔn)k～ε紊流模型可滿足計算精度要求。

圖1 計算區(qū)域

標(biāo)準(zhǔn)k～ε紊流模型是半經(jīng)驗公式[10],主要是基于紊動能和紊動能耗散率。其中k方程是精確方程,ε方程是由經(jīng)驗公式推導(dǎo)出的方程。

紊動能k方程

紊動能耗散率ε方程

式中 :μt=為紊流粘性系數(shù);Cμ為常數(shù);Gk為由于速度梯度引起的應(yīng)力生成項:

標(biāo)準(zhǔn)k～ε紊流模型方程中的常數(shù)通過經(jīng)典實驗得到,各相關(guān)常數(shù)如下:

1.3 邊界條件

1.3.1 進口邊界

取在進水池上游較遠(yuǎn)處,認(rèn)為來流已是充分發(fā)展的紊流,比較均勻。進水池入口的流速場給定,在水深方向設(shè)為對數(shù)式分布。

1.3.2 出口邊界

取在泵吸水管出口處,認(rèn)為出口流動已成單向狀態(tài),計算區(qū)域的解不受出口下游影響。

1.3.3 固壁邊界

進水池中的固體壁面包括邊墻、吸水管及其圓柱輪轂體等,采用標(biāo)準(zhǔn)壁面函數(shù)法,并在固體邊壁處規(guī)定無滑移條件。

1.4 網(wǎng)格劃分

進水流道采用非結(jié)構(gòu)化四面體網(wǎng)格和六面體網(wǎng)格相結(jié)合的混合網(wǎng)格,以適應(yīng)隔墩、喇叭口等不規(guī)則部件。

2 流場計算結(jié)果及分析

2.1 流場計算工況點的選取

通過數(shù)值模擬,可以得到進水流道內(nèi)流場的速度、壓力、流線等信息。為了便于分析、說明和比較,選取三個典型的工況點進行分析,分別是最優(yōu)工況(Q=20 m3/s)、小流量工況(Q=15 m3/s)和大流量工況(Q=25 m3/s)。

以下各圖中壓強單位為Pa,流速單位為m/s。

2.2 粒子跡線圖

圖2為最優(yōu)工況下整個流道的粒子跡線圖,由圖可知,水流由進水池進入喇叭口過程中直線段的流線均較為平順,當(dāng)進入喇叭口下方時流線急劇彎曲,由于進水流道的后部空間較大,水流的自由度大,在來流的帶動下流道的后壁死水區(qū)流線紊亂,形成緩慢回流和漩渦,但該區(qū)流速較小。

圖2 進水流道內(nèi)的粒子跡線圖

2.3 壓力和流線分布

為了進一步研究進水流道內(nèi)的流動情況,在最優(yōu)工況下,截取了幾個具有代表性的斷面進行分析,分別是進水流道隔墩兩側(cè)的對稱面(Z=-0.5 m,Z=0.5 m),靠近進水流道面層的斷面(y=-0.5 m),靠近進水流道底層的斷面(y=-2.7 m),具體位置見圖3(a)和圖4(a)。圖3(b)、圖3(c)、圖4(b)、圖4(c)顯示的是不同斷面上的速度等值線和壓力云圖。由圖可知,從前池到進水流道,壓力分布沿著水流的方向逐漸降低,流速逐漸增大,流態(tài)平順。但在進水流道的后部是死水區(qū),流速很小,且有明顯的漩渦。從縱斷面來看,進水流道內(nèi)在隔墩的兩側(cè)斷面上的壓力分布和流速分布十分接近,說明水流基本上是對稱的,死水區(qū)的漩渦靠近后壁;從水平斷面來看,靠近流道底部的流速和壓力比靠近流道面層的流速和壓力分布更均勻些,靠近面層的斷面上在喇叭口后部出現(xiàn)了對稱的一對漩渦。

圖3 不同縱截面速度等值線圖和壓力云圖

圖4 不同水平截面速度等值線圖和壓力云圖

2.4 橫向流速

圖5是進水流道內(nèi)從進口到出口不同橫截面(從左到右依次編為 1-1、2-2、3-3、4-4、5-5)上的橫向流速分布云圖,由圖5可知,1-1至4-4斷面的橫向流速均較小,其中3-3斷面的橫向流速最小,斷面最大的橫向流速僅為0.08 m/s。由于靠近吸水喇叭口,5-5斷面表現(xiàn)為靠近隔墩橫向流速大,靠近流道底部和頂部橫向流速小的特征。

圖5 不同橫截面上的橫向流速云圖

2.5 水力損失計算

在計算中,可利用表面積分功能,得到流道進出口的總壓,流道水力損失 Δ h由下式算得:

式中:P0out為流道出口的總壓;P0in流道進口的總壓;PS為表面靜壓;u為流體速度。

各工況的數(shù)值計算水力損失值與流量關(guān)系曲線如圖6所示。由圖6可知,對不同工況的計算表明,進水流道內(nèi)的水力損失隨著流量的增大而逐漸增大,并與流量基本上成二次方關(guān)系,這與文獻[11]得到的結(jié)論相符。

圖6 進水流道內(nèi)水力損失

3 結(jié) 語

(1)基于CFD數(shù)值模擬,獲得了箱形進水流道內(nèi)部流動的規(guī)律,從前池到進水流道直段水流順直,流態(tài)平穩(wěn),當(dāng)進入喇叭口下方時流線急劇彎曲,在進水流道后部死水區(qū)有明顯的漩渦。

(2)分析了進水流道內(nèi)不同斷面的流速和壓力云圖,從前池到進水流道,壓力分布沿著水流的方向逐漸降低,流速逐漸增大,流態(tài)平順。但在進水流道的后部靠近后壁和喇叭口的位置出現(xiàn)漩渦。

(3)通過對不同工況內(nèi)水流流動的數(shù)值模擬,預(yù)測了不同工況下進水流道內(nèi)的水力損失,其規(guī)律是水力損失隨著流量的增大而逐漸增大,并基本呈現(xiàn)二次方的關(guān)系。

[1]劉超.水泵及水泵站[M].北京:中國水利水電出版社,2009.

[2]劉超,周濟人,湯方平,等.低揚程雙向流道泵裝置研究[J].農(nóng)業(yè)機械學(xué)報,2001,32(1):49-51.

[3]周濟人,劉超,袁家博,等.大型泵站箱涵式雙向進水流道試驗研究[J].揚州大學(xué)學(xué)報(自然科學(xué)版),1999,2(4):79-82.

[4]葛強,陳松山,王林鎖,等.涵箱型雙向流道在船塢泵站應(yīng)用研究[J].水泵技術(shù),2004,(4):40-43.

[5]葛強,陳松山,王林鎖,等.鐘形箱涵式進水流道泵裝置特性模型試驗研究[J].水力發(fā)電學(xué)報,2006,25(5):129-134.

[6]王林鎖,陸偉剛,陳松山,等.“工字型”雙向流道泵站裝置特性試驗研究[J].水力發(fā)電學(xué)報,2001,20(2):79-85.

[7]LI Yaojun,WANG Fujun.Numerical investigation of performance of an axial-flow pump with inducer[J].Journal of Hydrodynamics,Ser.B,2007,19(6):705-711.

[8]Constantinescu G S,PatelV C.A numericalmodel for simulation of pump-intake flow and vortices[J].Hydr Engrg,ASCE,1998,124(2):123-134.

[9]成 立,劉超,薛 堅,等.基于CFD流動分析的泵站肘形進水流道水力特性研究[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報,2008,16(6):891-899.

[10]Rodi W.Turbulence models and their application in hydraulics experimental and mathematical fluid dynamics[C]//Delft:IAHR Section on Fundamentals of DivisionⅡ,1980:44-46.

[11]湯方平,袁家博,周濟人.軸流泵站進出水流道水力損失的試驗研究[J].排灌機械,1995,(3):13-14.