泵站側(cè)向進水前池整流措施數(shù)值模擬

張亞莉，宋世露，陳 勇，劉 霞，符向前

(1.甘肅省水利廳水利管理局，蘭州 730000；2.武漢大學動力與機械學院，武漢 430000；3.肇慶市景豐聯(lián)圍管理局，廣東 肇慶 526000；4.湖北省水利水電科學研究院，武漢 430000)

在泵站工程的進水建筑物中，前池、進水池是十分重要的組成部分。前池設置于引渠和進水池之間，它的作用是將引渠和進水池合理地銜接起來，使水流平穩(wěn)且均勻地流入進水池，為進水池提供良好的流態(tài)[1]。目前，國內(nèi)外學者對泵站進水結(jié)構(gòu)做了大量的研究，但對前池的研究則較為少見，尤其是對側(cè)向進水前池的研究甚少。前池內(nèi)的不良流態(tài)將引起池內(nèi)泥沙淤積，引發(fā)水泵機組的振動、氣蝕并降低水泵的效率[2]。

文中針對泵站側(cè)向進水前池的進水流態(tài)，采用數(shù)值模擬方法對前池尺寸、隔墩[3]對前池流態(tài)的影響進行計算，并將計算結(jié)果與試驗結(jié)果進行對比，同時比較不同隔墩形態(tài)對前池整流效果的影響，為改善側(cè)向擴散前池流態(tài)提供一定參考。

1 數(shù)學模型及計算方法

1.1 數(shù)學模型

本文中的泵站前池及進水池直接與大氣相通，可視為不可壓黏性流體的湍流流動進行三維模擬。選擇Realizablek-ε湍流模型模擬該流動。前池內(nèi)流體的流動狀態(tài)寫成張量形式，用連續(xù)性方程和動量方程來描述。

連續(xù)性方程和動量方程：

(2)

式中：ui,uj(i=1,2,3;j=1,2,3)為流體的速度分量；ρ為流體的密度；P為壓強；μ為流體的動力黏性系數(shù)；fi單位質(zhì)量的摩擦阻力。

1.2 模型布置

某泵站屬于側(cè)向進水泵站。泵站總裝機臺數(shù)3臺，3臺水泵同時運行。泵站總設計流量6.6 m3/s，單泵設計流量為2.2 m3/s。泵站進水池設計水位為6.5 m，最高水位7.0 m，最低水位6.0 m，停機水位5.8 m。

從取水管道到前池、進水池部分按照比例全部模擬，見圖1。每臺泵后接一管道，每條管道各安裝一個電磁流量計和流量調(diào)節(jié)閥，用以調(diào)節(jié)和監(jiān)測流量。管道連接到貯水池中，形成一個水循環(huán)系統(tǒng)。水泵為臥式離心泵，安裝在地面上，試驗時先用一小泵將吸水管抽真空，待吸水管內(nèi)充滿水即可開機進行試驗。

1-貯水池；2-引水管道；3-進水前池；4-進水池 ；5-循環(huán)管道；6-電磁流量計；7-流量調(diào)節(jié)閥圖1 模型布置圖

1.3 計算區(qū)域與網(wǎng)格

本文的主要研究對象為泵站前池與進水池，因為模型結(jié)構(gòu)比較復雜，故采用非結(jié)構(gòu)化網(wǎng)格對模型進行劃分。原方案網(wǎng)格劃分結(jié)果如圖2所示。

圖2 計算區(qū)域的網(wǎng)格劃分

通過網(wǎng)格無關(guān)性驗證，將網(wǎng)格劃分從90萬到150萬進行試算，之間結(jié)果相差不超過1%。最終網(wǎng)格總數(shù)在取100萬左右，隨著前池寬度和長度的變化，渦流防止裝置的布置，以及水位的變化而略有不同。

1.4 邊界條件

進口位置取在進水前池上游的半徑為300 mm的取水管道入口位置，采用速度進口條件。出口位置取在泵吸水管的入口斷面處。本文中，因未知出口斷面的速度和壓力，并假設出口處流動已完全發(fā)展，故在出口位置延長直徑100 mm，長1 m的管道。此時，可認定出口邊界對計算區(qū)域的影響較小，速度梯度為0。出口邊界條件設置為自由出流。固壁表面符合無滑移壁面條件，設為無滑移邊界條件(x，y，z方向上的速度分量均為0)。假設前池及進水池流動達到穩(wěn)定，即水體自由表面保持不變，自由表面上邊界條件指定為對稱面邊界條件，符合靜水壓力假定，沒有變量穿越該表面。

2 方案設計及結(jié)果

泵站進水前池未經(jīng)處理，開啟1～3號水泵，水位為最高水位7.0 m時，0.6H截面處的流速分布圖，見圖3。

圖3 水位7.0 m，開啟1～3號水泵，0.6H截面流速分布圖

2.1 方案設計

改善前池流態(tài)的方法[4]可以分為2大類。一類是改進前池的尺寸，一類是在前池尺寸確定的情況下增加改善流態(tài)的措施[5]。為了全面的研究這兩類情況對前池流態(tài)的影響，共設置6個方案，見表1。方案1～4主要改變前池的尺寸，方案5設置隔墩，方案6調(diào)整隔墩為翼型，觀察以下情況下該泵站的流速分布圖。

表1 各方案對前池流態(tài)的影響

2.2 方案效果分析

比較圖4(a)～4(d)，可得出在加寬、減寬、加長、減短前池的情況下，從截面流速分布圖4中可以看出，前池內(nèi)依然出現(xiàn)明顯的大回流漩渦。1號機組、2號機組進水池進口的水流出現(xiàn)明顯的偏流，只有3號機組水流是平行、均勻地進入前池的。

比較圖4(e)、4(f)，可以看出，增設隔墩之后1號、2號機組進水池進口的水流發(fā)生了變化，方案5中在與引水管道出口相連的泵站前池設置兩個隔墩。隔墩寬0.5 m，兩隔墩與原隔墩成60°角，以期望達到分流及導流的作用。隔墩端部與原隔墩對齊，布置圖如圖5。但是1號機組、2號機組進水池進口的水流出現(xiàn)明顯的偏流，只有3號機組水流是平行、均勻、無渦地進入進水池的。方案5比原方案在前池流態(tài)上改善不大。圖4(f)中，1號、2號機組進水口水流流態(tài)發(fā)生了明顯變化，流態(tài)基本平穩(wěn)無渦，起到了明顯的調(diào)節(jié)作用，方案6方法有效。

圖4 各方案在水位7.0 m，開啟1～3號水泵，0.6H截面流速分布圖

圖5 方案5隔墩布置圖

繞孤立翼流動是改善流場結(jié)構(gòu)的有效措施之一[6,7]。根據(jù)方案5的試驗結(jié)果分析，以及翼型原理，設計了方案6。設計兩個翼型，具體形狀及布置見圖6。兩翼型隔墩將來流均分為3等份。

圖6 方案6翼型隔墩布置圖

方案6實現(xiàn)了改善側(cè)向進水前池水流流態(tài)，因此在方案6的基礎(chǔ)上，調(diào)整隔墩尺寸，以期達到更好的調(diào)節(jié)作用。調(diào)整后速度分布圖如圖7所示。

圖7 調(diào)整翼型隔墩尺寸后速度分布圖

3 結(jié) 語

基于CFD技術(shù)，分析了湍流模型對泵站的前池內(nèi)流體流態(tài)的適用性，并繪制相關(guān)模型，分析了不同方案下對前池流態(tài)的影響。得出了以下結(jié)論。

(1)相對于標準k-ε模型、RNGk-ε模型，Realizablek-ε模型能夠更加精確的對某泵站前池流態(tài)模擬。

(2)引水管道、進水前池與進水池三者的流態(tài)息息相關(guān)，不良的引水管道形式會對前池流態(tài)產(chǎn)生不良影響，進而影響泵站水力穩(wěn)定性。因此，要改善進水流道的流態(tài)，就必須針對引水管道出口水流分布的特點，設置一定的導流整流措施改善前池流態(tài)，使進水前池及進水流道的水流均勻、順直、平穩(wěn)、無渦。

(3)對于某泵站來說，前池寬度對流態(tài)影響不大；在地形地勢及經(jīng)濟條件允許的情況下，加長前池的確能夠獲得更平穩(wěn)的流態(tài)。

(4)針對某泵站引水管道不良形態(tài)來流，采用翼型隔墩的整流的方式是可行的。翼型隔墩的長短、形狀、位置應根據(jù)來流形態(tài)，漩渦位置等綜合考慮，并不斷嘗試以最終獲得良好的整流效果。

□

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

[2] 張賢明，吉慶豐．泵站前池流態(tài)的數(shù)值模擬[J]．灌溉排水，2001,20(1):35-38.

[3] 周龍才．泵站前池隔墩整流的數(shù)值分析[J]．長江科學院院報，2010,27(2):31-33.

[4] 高傳昌，劉新陽，石禮文，等.泵站前池與進水池整流方案數(shù)值模擬[J].水力發(fā)電學報，2011,30(2):54-59.

[5] 成 立，劉 超，周濟人，等.泵站側(cè)向進水流態(tài)及其改善措施機理探討[J].排灌機械，2001,19(1):31-34.

[6] 谷傳綱，陳 西，苗永淼.不可壓縮實際黏性流體繞流孤立翼型的數(shù)值研究[J].空氣動力學報，1988,6(2):190-196.

[7] 關(guān) 鍵，郭 正.繞翼型低雷諾數(shù)流動的數(shù)值仿真[J].科學技術(shù)與工程，1988，13(24):7 275-7 281.