基于ANSYS-FLUENT模擬計算下泵站池流場特性分析研究

張世發(fā)

(南陵縣機電排灌管理站(養(yǎng)馬灘排澇站工程管理處)， 安徽 南陵 241300)

1 概述

由于地區(qū)水資源分布匱乏，需考慮從富水地區(qū)調(diào)度水資源[1-3]，因而需要修建泵站引水工程，服務(wù)于工農(nóng)業(yè)及生活需水。泵站作為重要水利工程，其長久穩(wěn)定運行能極大促進地區(qū)水資源供需，提高水資源開發(fā)度與高效利用，泵站穩(wěn)定運行與泵站各部分附屬結(jié)構(gòu)工程息息相關(guān)，其中泵站前池為聯(lián)通引水池與渠道的重要載體，前池內(nèi)流體運動穩(wěn)定關(guān)乎輸水安全，故而針對泵站前池流場分析具有重要意義[4-5]。已有基于水動力學理論研究了二維狀態(tài)下泵站前池運行狀態(tài)，準確預(yù)判泵站運行過程中前池流場特征[6-8]。另外還有學者采用水工模型試驗，如高傳昌等[9]研究了前池內(nèi)中心區(qū)域與末端壁面處流態(tài)特性。當然，亦可采用流場數(shù)值模擬手段[10-13]，研究不同狀態(tài)下的前池流場，為泵站安全運行提供重要參考，保障泵站引水工程安全。

2 流場分析理論

由于本文研究對象為泵站前池流場特性，而在前池流場中并不是僅存在流體介質(zhì)，另還存有懸浮泥沙等固體物質(zhì)，因而研究流場特性不應(yīng)忽視多場耦合狀態(tài)，故本文引入Mixture多相耦合模型方程。

Mixture模型連續(xù)性方程可表述為：

(1)

運動平均速度為：

(2)

式中αk為第k相體積分子。

故得到動量方程為：

(3)

流體質(zhì)點間相對速度可采用下式計算：

(4)

為保證Mixture模型在有限元插分計算時，不會出現(xiàn)多相場之間某一相場抵消與增強，因而將運動相對速度設(shè)定為代數(shù)關(guān)系，即為：

(5)

第二相場質(zhì)點運動時間τqp為：

(6)

式中dp為第二相場質(zhì)點尺寸。

擾力系數(shù)fdrag為：

(7)

即在不考慮相場間的運動干擾而致出現(xiàn)滑移速度，在流場外第二相場的體積微分方程為：

(8)

另外由于泵站前池特點，壁面邊界條件不同于流場與懸浮泥沙固體場，需引入近壁模型表述流場邊界面的壁面狀態(tài)，方便于求解邊界上的流場形態(tài)。近壁模型以修正湍流運動長度，如下式：

(9)

修正后湍流運動粘性系數(shù)：

(10)

式中fμ指雷諾修正系數(shù)。

流場邊界面上的運動阻尼函數(shù)為：

fμ=1-exp[-(a1R+a2R2+a3R3+a4R4+a5R5)]

(11)

壁面與多相場接觸系數(shù)：

(12)

通過上述分析，選擇合適的有限元單位體構(gòu)建模型，利用ANSYS-FLUENT有限元求解與流場后處理軟件，獲得泵站前池場中流態(tài)特征。

3 工程概況

西北某地區(qū)修建有一引水工程，承擔著區(qū)域蓄水調(diào)水水資源調(diào)度重要作用，該水工建筑包括有前池、出水池以及其他控制設(shè)備建筑。泵站設(shè)計提水流量為11 m3/s，每年提水流量指標為1.5×108m3，承擔農(nóng)田灌溉面積為2萬hm2，共建泵站10臺，裝機電量超過7×104kW，極大緩解了地區(qū)內(nèi)缺水，為當?shù)匦藿ňG色防護林提供重要水資源來源。設(shè)計輸水管直徑為1.4 m，泵站裝機高程達1 610 m，前池底部高程為1 613 m，池內(nèi)流量為15 m3/s，采用側(cè)向進水結(jié)構(gòu)形式，池內(nèi)設(shè)置有隔墩，另在中部區(qū)域設(shè)置有底坎，減弱水流對前池沖擊磨損。輸水渠道貫通至區(qū)域內(nèi)中度缺水狀態(tài)鄉(xiāng)鎮(zhèn)，緩和當?shù)剞r(nóng)業(yè)需水與生活需水缺口。由于泵站設(shè)計使用年限較久，不僅出水池有較多泥沙懸浮，在泵站前池中泥沙含量較高，淤泥沉積量已達池容積50%，長此以往對水道上的輸水管以及涵渠等水利設(shè)施均是較大的磨損，甚至會導(dǎo)致管路堵塞，影響泵站輸水效率。泵站前池在側(cè)向設(shè)計有進水口，水流方向平行于水泵，由于輸水管設(shè)計原因，泵站正向、側(cè)向視圖前池現(xiàn)狀如圖1所示。

圖1 泵站前池現(xiàn)狀示意

現(xiàn)場地質(zhì)踏勘及水質(zhì)檢測表明，池內(nèi)懸浮泥沙粒徑分布在0.02～1 mm，其中小于0.075 mm的粒徑占比達到90%以上，即懸浮泥沙基本上都屬細粒沙?，F(xiàn)場觀測泵站前池內(nèi)水流狀態(tài)，水流方向與輸水管產(chǎn)生交接，水流狀態(tài)較為錯亂，流場局部區(qū)域可見顯著回漩渦，加重水流的往返運動，造成輸水能耗多次消耗，增加輸水成本。另外，泵站場址內(nèi)以耕種栗鈣土和砂壤土為主，厚度最大約為7 m，分布普遍是較薄的砂壤土，由于水土流失嚴重，部分地表可見磨蝕的巖層，巖性主要為花崗巖，呈半風化特征。利用超聲波流速測試儀測試現(xiàn)場實地，水面層流速約為0.02～0.55 m/s，流速處于較動蕩分布狀態(tài)，池內(nèi)進水口左右兩側(cè)流速呈現(xiàn)一高一低，在中間地帶會有局部漩渦，考慮現(xiàn)場監(jiān)測儀器觀測僅僅停留在表面層，對池內(nèi)各深度流態(tài)特征觀測缺乏，因而本文將綜合工程資料，以數(shù)值手段分析泵站流態(tài)特征。

4 泵站前池流態(tài)特性

4.1 原設(shè)計前池與淤積后前池

前池結(jié)構(gòu)示意如圖2a所示，其基本形態(tài)及結(jié)構(gòu)尺寸均在結(jié)構(gòu)平面圖中標出，方便在ANSYS中建立數(shù)值模型(如圖2b所示)，由于前期實地觀測得知該泵站前池淤泥沉積已超過容積1/2，故建立淤積狀態(tài)下前池模型，并在ANSYS中劃分出網(wǎng)格單元(如圖3所示)。

圖3 淤積前池模型示意

a 前池剖面

原設(shè)計前池利用SOLID65為基本構(gòu)成微單元體，共劃分2 001 064個網(wǎng)格，節(jié)點數(shù)845 848個；淤積狀態(tài)下前池利用SOLID45與SOLID65模型構(gòu)成單元體，共劃分1 962 654個網(wǎng)格，節(jié)點數(shù)783 671個。模擬工況以1#、2#、3#、4#、5#、6#泵站開啟，7#、8#、9#、10#泵站關(guān)閉狀態(tài)，且含沙量為2.1 kg/m3、泵站流量為10.6 m3/s、攔污柵流速值為0.45 m/s時，開展模擬計算。在ANSYS-FLUENT求解中，以ANSYS計算兩相場基本特征參數(shù)，包括水力半徑、耗散能等，再利用FLUENT開展后求解分析，獲得側(cè)向前池內(nèi)水沙兩相場流場特征。

圖4為前池原設(shè)計模型不同深度處池內(nèi)流場分布示意，包括分布云圖與矢量圖。從圖4中可看出，原設(shè)計前池水流表面的流場內(nèi)存在回流現(xiàn)象，約占流場分布面積的50%，最大流速超過1.05 m/s，且表層流速分布范圍較廣，從0～1.05 m/s均有涉及，并未出現(xiàn)有主流速，即每個范圍流速均有一定流動性，故而流態(tài)分布較為紊亂；最大流速主要位于池壁邊緣處，漩渦回流主要位于池內(nèi)中間及部分與水流方向呈較大夾角之處，分析表明表面層由于水泵提水作用，攪亂了水流主向進入前池內(nèi)的平行方向，進而水流方向呈現(xiàn)一定傾斜，漩渦回流顯現(xiàn)，流場分布矢量圖中亦可看出。隨著計算水深增大至1 m、3 m，漩渦流速分布范圍有所降低，且均出現(xiàn)主流速，深度1 m處主流速范圍為0～0.22 m/s，大流速分布面積縮減，入口處流速相比水面表層降低，在各深度壁面處由于泵站運行影響，隨深度增大，流速值增大，且愈深出現(xiàn)漩渦回流范圍愈少，深度3 m處漩渦回流分布范圍最少，表明前池原設(shè)計考慮到了泵站運行與吸水管對池內(nèi)流速方向擾動的雙影響，通過控制吸水管流量與泵站運行頻率，減少了對池內(nèi)更深水流的影響，僅在水面表層出現(xiàn)回旋渦等現(xiàn)象。

(a)水流表面

圖5為淤積狀態(tài)下泵站前池流場分布云圖示意。從圖5中可看出，相比原設(shè)計模型水流表面層，淤積狀態(tài)下由于泥沙懸浮影響，在邊壁面處聚集較多漩渦回流，各個范圍流速階梯性顯著，形成過多小漩渦，且漩渦總分布面積超過池內(nèi)面積2/3，即淤積狀態(tài)下回流漩渦范圍更廣了。另外由于過多漩渦占據(jù)了水流表面層，導(dǎo)致主流流速只能傾向于左側(cè)邊壁流動，且愈向池內(nèi)中間區(qū)域，主流愈發(fā)無法占據(jù)控制流場的位置，極大影響了前池內(nèi)輸水效率。對比不同深度段前池內(nèi)流態(tài)特征可發(fā)現(xiàn)，隨深度增大，由于泥沙淤積在池底，在更深水流層中回流紊亂區(qū)域更多，且中部漩渦逐漸在深度3 m處時逐漸擴散至前池進水口處，造成進水口處流速愈大，但整體流速值相比淺水深處更小了，3 m處進水口最大流速為0.44 m/s，相比水流表面層降低了50%，較低水流在淤沙沖擊作用下會局部產(chǎn)生泥沙懸浮或沉積，進水口淤積較多泥沙即使如此，造成進水段面積降低，流速分布出現(xiàn)較為紊亂狀態(tài)，主流與回漩渦接觸面處流速分布較低，且進水口的紊亂水流逐步傾向于往水泵吸水管附近靠攏，并聚集著較大流速。對比原設(shè)計模型前池流態(tài)，由于淤積前池內(nèi)泥沙存在，進一步擠壓主向水流運動空間，造成主流速度降低，矢量方向上更靠近左側(cè)邊壁，回流漩渦矢量靠近水泵吸水管，分布范圍增多，經(jīng)計算得知吸水管附近水流方向與之平均夾角近90°，輸水效率嚴重降低，泵站安全運營受到極大挑戰(zhàn)。

(a)水流表層

4.2 改善后流態(tài)特征

前述分析已知淤積后前池內(nèi)流場處于較紊亂狀態(tài)，泵站運行極大受限，故而考慮對該前池采取一定改善措施，在前池內(nèi)布置壓水板與底坎相鉸接，夾角為45°，其中壓水板設(shè)置在距離前池入口23 m、20 m，深度取0.8 m、1.2 m、1.5 m，共6個對比方案，分別對該6種改善措施方案開展流態(tài)計算分析。限于篇幅，本文列出6種方案在水深2 m吸水管附近的流場分布特征(如圖6所示)。

從圖6可看出，架設(shè)45°壓水板后，不論壓水板距離入水口多遠，池內(nèi)吸水管附近回流渦旋數(shù)量顯著減少，流線平衡均勻，流速亦呈增大態(tài)勢，可帶動泥沙沖擊至出水口，提高了前池內(nèi)流體沖淤能力。對比同等水深不同壓水板與入水口距離流場特征可看出，當距離較近時，壓水板存在導(dǎo)致底坎前端流速值較大，導(dǎo)致主流過于控制流體運動方向，進而在底坎末端出現(xiàn)局部的小漩渦，流體翻滾轉(zhuǎn)動影響吸水管進水效率，故而考慮壓水板距離入水口23 m為最佳。對比不同深度布置壓水板可知，當壓水板距離23 m、深度為0.8 m時，未運行的吸水管水泵附近可見回流現(xiàn)象，泥沙懸浮顯著，易堵塞輸水管；當深度為1.5 m時，由于壓水板布置深度較深，入水口流速過大，直接擾亂了其他底坎前端小流速流體運動，形成一定數(shù)量的漩渦。綜上對比分析，表明設(shè)計深度為1.2 m時吸水管附近流場穩(wěn)定，既能裹挾泥沙，又不至于引起局部漩渦，故以壓水板距離入水口23 m、深度1.2 m為最優(yōu)改善方案。

(a)距離23 m、深度0.8 m (b)距離23 m、深度1.2 m

5 結(jié)語

1) 研究了原設(shè)計前池不同深度流場分布，水流表層分布約有1/2的回流漩渦，漩渦流速隨深度遞減，回流減弱，深度1 m處主流速范圍為0～0.22 m/s，池壁面處流速隨深度增大而增長。

2) 獲得淤積前池邊壁面存在較多漩渦回流，分布范圍超過池內(nèi)面積2/3，隨深度增大，回流紊亂區(qū)域更多，進水口流速隨深度增大逐漸降低，3 m處最大流速為0.44 m/s，相比水流表面層降低了50%，流速過小加劇了泥沙淤積；主流相比原設(shè)計前池更靠近左側(cè)邊壁，吸水管與水流方向夾角近90°，輸水效率較低。

3) 分析了壓水板與底坎6種組合的改善方案對流場分布影響，距離入水口過近，易導(dǎo)致底坎前端流速過大，在末端出現(xiàn)聚集漩渦；深度1.2 m時流場穩(wěn)定，流速適佳，既能裹挾泥沙，又不至于引起局部漩渦，故設(shè)置壓水板距離入水口23 m、深度1.2 m為改善方案。