金曉宇，常鵬程，胡文竹，曹立力，張一祁

(1.揚州大學 水利科學與工程學院，江蘇 揚州225009； 2.揚州市勘測設計研究院有限公司，江蘇 揚州225009)

0 引 言

立式軸流泵具有低揚程、高流量、運行穩(wěn)定、易于安裝檢修等特點[1]，被廣泛應用于防洪排澇、農(nóng)田灌溉、田間排水、跨區(qū)域調(diào)水等多個領域[2-4]，在國民經(jīng)濟和發(fā)展中具有重要的地位和作用。當前，國內(nèi)外學者對立式軸流泵裝置已開展較多的研究。戴景等[5]采用Fluent16.0與UDF二次開發(fā)程序相結(jié)合的數(shù)值模擬方法研究了立式軸流泵裝置的飛逸過渡過程，獲得機組在飛逸過渡過程中流量等特性參數(shù)的變化以及肘形流道內(nèi)流線的演化特性。朱曉明等[6]采用Standardk-ε、RNGk-ε和Realizablek-ε湍流模型的RANS方法，研究了不同流量下肘形進水流道內(nèi)的壁面壓強，得出流量和壓差的關系。王麗慧等[7]采用基于三維湍流流動雷諾時均N-S方程和RNGk-ε湍流模型，研究了不同出水流道組成的立式軸流泵裝置模型內(nèi)部流動，得出泵裝置肘形進水流道內(nèi)的流場分布情況。本文采用數(shù)值模擬的方法對立式軸流泵裝置全流道開展數(shù)值模擬分析工作，分析不同工況時立式軸流泵裝置全流道的壓力和速度分布規(guī)律，探討不同工況時各過流結(jié)構(gòu)的能量變化。

1 數(shù)值計算模型

1.1 三維模型

立式軸流泵裝置包括4個過流部件：肘形進水流道、葉輪、導葉體、虹吸式出水流道。葉輪的葉片數(shù)為4，葉片安放角為0°，葉輪名義直徑為120 mm，輪轂比為0.40，額定轉(zhuǎn)速為2 200 r/min；導葉體的葉片數(shù)為7 。立式軸流泵裝置的三維模型見圖1。

圖1 立式軸流泵裝置

1.2 控制方程與湍流模型

數(shù)值計算采用雷諾時均Navier-Stokes(RANS)控制方程，并將其基于有限元的有限體積法(CV-FEM)進行離散化，離散方程采用全隱式耦合代數(shù)多重網(wǎng)格方法進行求解。低揚程立式軸流泵屬于不可壓縮三維湍流流動，而適用于此流動的湍流模型有多種，其中RNGk-ε湍流模型具有適于計算流線曲率大、有旋轉(zhuǎn)和分離流動的優(yōu)點[8-9]，能有效反映葉輪旋轉(zhuǎn)產(chǎn)生的動靜干涉問題，已在水泵及泵裝置三維湍流流動數(shù)值模擬中得到較多應用[10-12]。

1.3 邊界條件設置

選用的立式軸流泵裝置包括4個過流部件：肘形進水流道、葉輪、導葉體、虹吸式出水流道。進口邊界設置于肘形進水流道進口延伸段處，采用速度進口；出口邊界設置于出水流道出口延伸段，采用壓力出口條件；壁面采用無滑移邊界條件；近壁區(qū)采用可伸縮壁面函數(shù)(Scalable Wall Function)處理；各物理量的殘差收斂精度設置為1.0×10-5，并通過設置監(jiān)測點監(jiān)測泵裝置揚程的變化情況。進水延伸段和出水延伸段的自由液面采用對稱平面處理，葉輪與導葉體、肘形進水流道的交界面采用動-靜交界面，其余區(qū)域的交界面選用靜-靜交界面處理。

1.4 網(wǎng)格無關性分析

本文采用ICEM CFD軟件對肘形進水流道和虹吸式出水流道進行結(jié)構(gòu)化網(wǎng)格劃分，通過Ansys TurboGrid軟件對葉輪和導葉體進行結(jié)構(gòu)化網(wǎng)格劃分。為了提升計算精度且兼顧計算效率，需要進行網(wǎng)格無關性分析，結(jié)果見表1。當網(wǎng)格數(shù)量增加到510萬以后，泵裝置的效率絕對誤差不超過0.2%，滿足精度要求，故選擇510萬網(wǎng)格作為計算網(wǎng)格數(shù)。

表1 不同網(wǎng)格數(shù)量泵裝置效率表

1.5 立式軸流泵能量性能預測

基于CFD數(shù)值計算結(jié)果，對立式軸流泵裝置的能量性能進行預測，得出泵裝置能量性能曲線(圖2)。當Q/Qbep=1時，泵裝置效率達到最大為74.71%，此時揚程為1.51 m。

圖2 泵裝置能量性能曲線

2 內(nèi)流特性結(jié)果及分析

為分析不同工況下立式軸流泵裝置的內(nèi)流特性，選取3個特征工況(小流量工況0.3Qbep、最優(yōu)流量工況1.0Qbep和大流量工況1.2Qbep)進行泵裝置內(nèi)流分析。

2.1 泵裝置全流道流線圖分析

由泵裝置全流道流線圖(圖3)可知，肘形進水流道內(nèi)部流態(tài)較好，流線平順，但在小流量工況下，肘形進水流道出口處流態(tài)較為混亂，存在大尺度螺旋。0.3Qbep工況下，水流經(jīng)過葉輪和導葉體后，由于流體流速較小，導葉體對水流的速度環(huán)量回收效果較差，導致在彎管段存在較為明顯的螺旋水流，在出水流道上升段水流旋轉(zhuǎn)減弱，流態(tài)較為平順；1.0Qbep工況下，導葉體對葉輪出口環(huán)量的回收效果較為明顯，但在彎管段仍然存在螺旋水流，隨著水流進入上升段，流態(tài)逐漸平順；1.2Qbep工況下，導葉體對葉輪出口環(huán)量的回收效果非常明顯，螺旋流進一步減弱，流出水流道中流線更為平順。

圖3 泵裝置全流道流線圖

2.2 泵裝置縱斷面壓力云圖分析

由泵裝置縱斷面壓力云圖(圖4)可知，在肘形進水流道段，縱斷面壓力分布規(guī)律基本一致，且在肘形進水流道出口段、葉輪域進口前存在一定范圍的低壓區(qū)。在0.3Qbep工況下，由于流態(tài)較差，壓力分布變化幅度較大，虹吸式出水流道壓力值高于肘形進水流道37 500 Pa左右；在1.0Qbep工況下，壓力分布變化幅度進一步減小，虹吸式出水流道壓力值高于肘形進水流道15 000 Pa左右，在葉輪域前側(cè)的低壓區(qū)沿著肘形進水流道內(nèi)側(cè)進一步擴大；1.2Qbep工況下，壓力分布變化幅度進一步減小，虹吸式出水流道與肘形進水流道無明顯壓力差距，彎管彎折處與出水流道出口下側(cè)有一定的高壓區(qū)，葉輪域前側(cè)低壓區(qū)進一步擴大，分布與1.0Qbep工況下類似。

2.3 泵裝置速度云圖分析

由泵裝置速度云圖(圖5)可知，3個典型工況下，縱斷面速度分布規(guī)律類似，主要差異在于肘形進水流道出口、葉輪域、導葉體域、彎管段和虹吸式出水流道上升段。隨著流速的增大，肘形進水流道出口處流速不斷擴大；葉輪域靠近輪緣處流速最大；在彎管內(nèi)部彎折處出現(xiàn)局部高速流動區(qū)域，且隨著流速的增大擴散，在1.2Qbep工況延伸到虹吸式出水流道上升段。

圖5 泵裝置縱斷面速度云圖

不同工況時，泵裝置各過流部件能量變化時均值見表2。隨著流量的增加，進水流道水力損失逐漸增大，水力損失從約占葉輪揚程的0.12%上升到3.29%；導葉體內(nèi)水力損失隨著流量的增大先減小后上升；彎管段水力損失隨著流量的增大逐漸減??；虹吸式出水流道水力損失隨著流量的增大逐漸上升，且從約占葉輪揚程的0.6%上升到13.9%。

表2 不同工況各過流部件能量變化 /m

3 結(jié) 論

1) 肘形進水流道內(nèi)部流態(tài)較好，隨著流量的增加，肘形進水流道出口處流態(tài)愈發(fā)平順，大尺度螺旋消失，流線平順導葉體對葉輪出口環(huán)量的回收效果非常明顯，螺旋流進一步減弱，出水流道中流線更為平順。

2) 肘形進水流道段，縱斷面壓力分布規(guī)律基本一致，且在肘形進水流道出口段、葉輪域進口前存在一定范圍的低壓區(qū)。隨著流量的增加，壓力分布變化幅度越來越??；至1.2Qbep工況下，壓力分布變化幅度進一步減小，虹吸式出水流道與肘形進水流道無明顯壓力差距，彎管彎折處與出水流道出口下側(cè)有一定的高壓區(qū)，葉輪域前側(cè)低壓區(qū)進一步擴大。

3) 3個典型工況下，縱斷面速度分布規(guī)律類似。隨著流速的增大，肘形進水流道出口處流速不斷擴大；葉輪域靠近輪緣處流速最大；在彎管內(nèi)部彎折處出現(xiàn)局部高速流動區(qū)域，且隨著流速的增大擴散，在1.2Qbep工況延伸到虹吸式出水流道上升段。

4) 隨著流量的增加，進水流道水力損失逐漸增大，導葉體內(nèi)水力損失隨著流量的增大先減小后上升；彎管段水力損失隨著流量的增大逐漸減?。缓缥匠鏊鞯浪p失隨著流量的增大逐漸上升。