基于Realizablek-ε湍流模型的泵站進水池水流流態研究

陶東 周冬蒙 劉子金

摘 要：泵站進水池進水喇叭口附近因抽吸力的作用而產生大量旋渦，進而嚴重影響水泵吸水效率。為探究進水池內部水流流態特性，基于Realizable k-ε湍流模型對泵站進水池水流流態進行數值模擬，并與HXH03-1S型超聲波多普勒流速流向儀在泵站進水池的原型監測結果進行對比，驗證數值模擬結果的可靠性。在此基礎上，對泵站進水池水流流動特性進行分析。結果表明：Realizable k-ε湍流模型對泵站進水池水流流態的模擬與HXH03-1S型超聲波多普勒流速流向儀在泵站進水池的原型監測結果吻合度較高，表明該數學模型模擬泵站進水池水流流態具有較高可靠性;喇叭口兩側出現相反方向的渦流，與吸入喇叭口內的水流形成新的渦帶，致使管道內的流態受到了前池流態的影響。

關鍵詞：泵站進水池;Realizable k-ε湍流模型;旋渦;數值模擬

中圖分類號：TV675;TV698.1+6 ? 文獻標志碼：A

doi：10.3969/j.issn.1000-1379.2020.11.006

Abstract：In the vicinity of the inlet bell of the pumping station， a large number of vortices are generated due to the suction force， which seriously affects the water absorption efficiency of the pump. In order to explore the flow characteristics of the water in the inlet pool， the flow regime of the pumping station was numerically simulated based on the Realizable k-ε turbulence model and compared with the prototype monitoring results of the HXH03-1S ultrasonic Doppler flow rate meter in the pumping station inlet tank to verify the reliability of the numerical simulation results. On this basis， the flow characteristics of the pumping station inlet water were analyzed. The results show that the coincidence degree of the flow pattern of the pumping station inlet and the prototype monitoring results of the HXH03-1S ultrasonic Doppler flow rate meter used in the pumping station is high， showing that has higher reliability. The vortex in the opposite direction on both sides of the bell mouth and the water flow in the suction bell mouth form a new vortex belt， so that the flow regime in the pipeline is affected by the flow state of the front pool. The research results have important guiding significance for guiding the optimization design of the front pool of large pumping stations and the flow analysis during the structural transformation process.

Key words： pump sumps; Realizable k-ε Turbulence Model; vortex; numerical simulation

泵站是灌區最主要的建筑物，而進水前池是泵站的重要組成部分，其主要作用是引導水流平穩地進入泵站進水池，保證水泵良好的進水條件。泵站運行中，前池流態是否平穩直接影響水泵吸水管的進水流態，進水池中的水流流態對泵站提水效率及水泵的氣蝕性能具有重要影響[1]。目前該領域學者針對泵站進水池內部渦流特性開展了大量研究，取得了豐富的研究成果。LI等[2]采用二維PIV系統，從PIV數據中得到流速、流線、湍流強度和渦度，通過對這些參數的分析，找出了具有特殊結構的集水池的流動特性，并在實驗的基礎上確定了最佳的流動條件;SUERICH-GULICK等[3]通過對進近流和渦旋內部流的速度測量，提出了一種半經驗模型，粗略地預測了進近氣流和進氣口幾何形狀如何決定核心渦特性;TEMMERMAN等[4]利用大渦模擬湍流，考慮流體黏度的影響，討論了速度分布的數值計算結果，并與ANSAR和NAKATO的實驗數據進行了比較，取得了較為理想的模擬結果;JAKIRLI等[5]采用大渦模擬技術對某小型水電設施進水口三維非定常流場進行了數值模擬，基于針對自由表面流動改進的高度簡化的MAC迭代法，采用近壁面阻尼修正的標準Smagorinsky模型對亞網格尺度應力進行了計算，利用數值網格對實驗室尺度雷諾數的近壁面流場進行了數值求解，較好地再現了渦旋流場的復雜性和非定常性;YU等[6]首次通過對黏性子層的充分求解，對SST k-ω湍流模型進行改進，并與標準的SST k-ω湍流模型進行比較，改進后的模型對氣體污染物在孤立建筑物周圍的擴散模擬有很好的適用性;為保障泵站穩定、高效運行，學者們運用多種湍流模型對前池的水流流動特性進行仿真模擬，并提出了多種整流措施[7]。以上研究為泵站進水池水流流態優化提供了可靠的理論依據。但受早期測流手段的限制，前人們的研究方法主要集中在物理模型試驗和數值模擬，對前池原型的測流手段較少，有針對性的泵站進水池流態理論研究尚不充分。

近年來，隨著測量技術的發展，測量結果的精度已有了較大幅度的提升，泵站前池流速監測技術已經在工程建設相關領域得到了廣泛應用，具有安全、高效、精度較高等優勢。鑒于此，本次研究基于Realizable k-ε湍流模型對泵站進水池水流流態進行數值計算，并與HXH03-1S型超聲波多普勒流速流向儀在泵站進水池的原型監測結果進行對比，驗證數值模擬結果的準確性。在此基礎上，對泵站進水池水流流動特性進行分析。

1 數學方法和計算模型

1.1 控制方程

1.2 邊界條件設置

①前池斷面尺寸以及泵站設計流量均已知，前池進口斷面處平均流速可求，故選用速度進口邊界條件;②出口壓力和流速均未知，假定出口邊界對上游水流運動無影響，水流視為完全發展流動，設置出口邊界條件為outflow邊界條件;③在近壁面區域，由于流體運動變化比較大，因此選用壁面函數法[8]處理固壁邊界;④泵站進水前池自由表面受外部環境擾動較小，浮動變化可以忽略，同時考慮到網格劃分以及對計算時間的要求，選用剛蓋假定法[9]對自由表面進行模擬計算，將自由表面條件設置為symmetry。

1.3 網格劃分和流場離散化

流動問題數值模擬計算的效率與精度主要取決于網格生成的質量及所采用的算法[10]。由于所選泵站側向進水前池結構形態復雜，因此采用非結構化四面體網格進行劃分，局部位置適當加密。應用Fluent軟件進行渾水兩相均質流模擬，基于有限體積法對計算區域進行離散化，采用二階迎風格式的隱式求解方法以確保計算精度，并選用基于分離求解器的SIMPLEC算法進行流場耦合。

1.4 模型構建

本研究采用的泵站進水池物理模型是按照泵站前池實際尺寸1∶1構建，如圖1所示。圖中吸水管內徑為1.6 m，泵管內流量為2.41 m3/s，喇叭口水深為6 m。管內雷諾數為75 000左右，進水池內雷諾數為10 000左右，Fr和Wb（韋伯系數）為模擬所需確定的相似準則數，其值分別為1.1和840。本研究取前池及進水池作為研究對象，采用非結構化網格進行劃分，在擴散角變化處及水泵進水口附近位置進行網格局部加密，計算區域及局部細節如圖2所示。研究工況是設計工況，其進口取在渠道漸變段與前池連接處，設置為體積進口，流量為8.5 m3/s。出口設置在進水池水泵進水喇叭口處，將自由表面簡化為對稱面并且不考慮其他表面波動。

2 模擬結果驗證

為進一步確認數值計算的準確性，本次研究對泵站前池水流流速采用HXH03-1S型超聲波多普勒流速流向儀進行測量，此型號流速儀采用最先進的電子技術，具有高度靈敏性，尤其對于高含沙水流抗干擾能力強，對弱流和強流都有很好的適用性，在選定的泵站圖紙上建立坐標系，先對泵站前池進行劃分，選取158個測點，測點布置如圖3所示，連接水下探頭沿不同深度對流速進行測量，得到每個測點不同深度的流速數據，再利用Surfer等值線繪制軟件繪制不同斷面處的流速等值線圖。

依據現場測流數據可知，各點流速范圍均在0.00～1.05 m/s之間，運用Surfer軟件分別繪制3個典型區域水深1 m、2 m處流速等值線圖，如圖4所示。圖3 進水前池現場測流點布置

區域Ⅰ受運行機組影響，在各個深度左側流速總體大于右側流速，隨著水深加大，右側低流速區域逐漸增大直到出現靜水區;區域Ⅱ側壁處整體流速較小，且出現較多旋渦回流，整體流速變化較小，旋渦在水深1 m處集中;區域Ⅲ水流入口處流速較大，且受側壁處旋渦的影響，主流有向左偏移的趨勢，隨著流速的增大低速區范圍逐漸變大。由此可見，機組打開的區域水流流速大于機組未打開的區域，且受側壁處旋渦影響主流向左偏移，與現場測流結果趨勢一致。

3 結果與討論

3.1 喇叭管下方速度分布分析

為驗證Realizable k-ε湍流模型在喇叭管周圍流態模擬的準確性，對喇叭口下方進行3種不同工況下的流速測量，圖5按照進水流道內各個方向速度平均值對其3個方向的速度值作無量綱處理，Vx、Vy、Vz分別代表x、y、z方向上的速度，Ly為測點距喇叭口的距離。

圖5（a）為3種工況（工況一、二、三的流量依次增大）條件下基于Realizable k-ε湍流模型與實測數據在x方向的速度分布對比，可以看出，Vx模擬值與實測值變化趨勢一致，隨著流量增大，在近壁處的速度梯度也隨之增大，而在非近壁區域隨工況變化沿x方向速度分布出現一定的波動性，但整體來說3種不同工況下沿x方向的速度分布呈現出相似的趨勢;圖5（b）為y方向的速度分布對比，沿著y方向速度呈現的趨勢為由正到負先增大后減小，其位置分別對應喇叭口測線的兩側，速度最大值出現在Ly的三分之一處，最小值出現在四分之三處;圖5（c）為z方向的速度分布，在此方向3種不同工況下速度分布都呈現出相同的變化趨勢，在壁面兩側因壁面效應出現了速度降低的現象。總的來說，基于Realizable k-ε湍流模型數值模擬與實測數據在3種不同工況條件下3個方向上均呈現相似的趨勢，吻合度較高。

3.2 水流流動特性分析

圖6、圖7分別為平行水流方向和垂直水流方向的平面流速分布。水流在整個前池中均有渦流產生，主要分成3段：靠近進口端，渦流分布在進水口兩邊側壁，渦流旋轉方向主要是切向水流方向翻轉帶動向前，渦帶較長;靠近前池中部的位置，渦帶位于取水前池的中下部，渦帶逐步縮短，表面流速較為均勻;前池尾部的渦帶主要分布在表面，受到側墻反射回來的水流的干擾，渦帶水流方向呈現出周期性。在各個泵取水喇叭口對應的位置上呈現出一定速度的渦，渦的位置偏離取水喇叭口的正前方，沿著水流方向，在兩臺泵取水喇叭口中間位置也出現了渦流，渦流的主要方向是順著水流方向向前，且渦帶較長，沿著水流方向，渦帶的長度在縮短變寬，沿著四周旋轉。沿著水流方向在喇叭口左端渦流方向向上，右端渦流向下，呈現出喇叭口兩側的水流存在搓動現象，與進入喇叭口的水流進行了較明顯的切向旋轉狀態的流動，生成新的渦帶，該類渦帶在喇叭口右前方擴大并拓展到水流的表面，喇叭口周圍產生的渦帶沿著側向前池的寬度方向減小。水體的流動在喇叭口前端產生了不均勻的渦帶，渦帶的旋轉方向不確定，在靠近前池末端的機組管道的前池中相互摻混頂托向上翻滾，在前池前端入口處的機組的水中渦順著入流方向傾斜旋轉，對前池未放置進水管的側墻產生二次沖擊，形成新的不規則的渦流，進而影響取水口取水，在取水口底部交匯出現渦流分布，旋渦產生區域隨著工況變化而變化，也顯示出一定的動態性和不間斷性，影響到水泵前端取水引流管道，致使管道內的流態受到了前池流態的影響。

4 結 論

（1）與HXH03-1S型超聲波多普勒流速流向儀實測數據對比，基于Realizable k-ε湍流模型對泵站進水池水流流態的模擬在數值上存在一定差異，但整體變化趨勢一致，表明該數學模型模擬精度較高，模擬結果可以很好地反映出進水前池的水流流態。

（2）喇叭口兩側出現方向相反的渦流與吸入喇叭口的水流形成了新的渦帶，對前池未放置進水管的側墻產生二次沖擊，在取水口底部形成新的不規則的渦流進而對喇叭口管道內部流態產生影響，表明當進水池結構尺寸設計不當時，側壁撞擊形成的旋渦危害是不可忽視的。

（3）對泵站進水池的數值模擬會略微高估進水池邊壁處的流速，主要原因是進水池兩側水流流速較中間部位小，導致大量泥沙淤積，而模擬過程中未考慮泥沙濃度對水流黏滯性的影響。對于高含沙水流的模擬是復雜的、多維的，高含沙水流顆粒間的相互作用以及水沙兩相的相互作用對進水池流態模擬有一定影響，后續研究將主要集中在這一領域。

參考文獻：

[1] 高傳昌，曾新樂，解克宇，等.泵站進水池超低水位下組合整流方案與驗證[J].農業工程學報，2017，33（23）：101-108.

[2] LI Y， WU Y， MANSA K， et al. The Flow Research in an Open Type Sump by PIV Experiments[C]//ASME Proceedings of ASME Heat Transfer/Fluids Engineering Summer Conference. New York：ASME， 2004：5-12.

[3] SUERICH-GULICK F， GASKIN S J， VILLENEUVE M， et al. Free Surface Intake Vortices： Theoretical Model and Measurements[J].Journal of Hydraulic Research， 2014， 52（4）：502-512.

[4] TEMMERMAN L， HADIABDI M， LESCHZINER M A， et al. A Hybrid Two-Layer URANS-LES Approach for Large Eddy Simulation at High Reynolds Numbers[J].International Journal of Heat and Fluid Flow， 2005， 26（2）：173-190.

[5] JAKIRLI S， KADAVELIL G， KORNHAAS M， et al. Numerical and Physical Aspects in LES and Hybrid LES/RANS of Turbulent Flow Separation in a 3-D Diffuser[J]. International Journal of Heat and Fluid Flow， 2010， 31（5）：820-832.

[6] YU H S， JESSE T. Simulation of Gaseous Pollutant Dispersion Around an Isolated Building Using the k-ω SST （Shear Stress Transport） Turbulence Model[J]. Air Repair， 2017， 67（5）：517-536.

[7] CONSTANTINESCU G S， PATEL V C. Numerical Model for Simulation of Pump-Intake Flow and Vortices[J]. Journal of Hydraulic Engineering， 1998， 24（2）：123-134.

[8] 資丹，王福軍，陶然，等.邊界層網格尺度對泵站流場計算結果影響研究[J].水利學報，2016，47（2）：139-149.

[9] 陸林廣，劉軍，梁金棟，等.大型泵站出水流道三維流動及水力損失數值計算[J].排灌機械，2008，26（3）：51-54.

[10] TRAN T T， KIM D H. A CFD Study into the Influence of Unsteady Aerodynamic Interference on Wind Turbine Surge Motion[J]. Renewable Energy， 2016，90（5）：204-228.

【責任編輯 張 帥】