大套一站軸流泵裝置事故飛逸特性研究

徐建葉，張海勝，劉躍飛，周大慶，鄭 源

(1.鹽城市通榆河樞紐工程管理處，江蘇 鹽城 224511；2.河海大學 能源與電氣學院，南京 211100)

當軸流泵裝置因突然斷電或者誤操作等原因造成事故停機，且泵出水管道上的閥門或者閘門等斷流設施失靈，則軸流泵裝置將由水泵工況轉至水輪機工況，當倒流流速達到最大且持續運行時，軸流泵裝置將處于事故飛逸狀態[1]。較大的飛逸轉速會對軸流泵裝置帶來損害。因而針對泵飛逸狀態的研究顯得尤為重要。

目前針對軸流泵飛逸特性的研究主要基于模型試驗[2-5]，通過對水泵模型進行飛逸特性試驗，獲得模型泵的單位飛逸轉速，進而換算出原型泵的飛逸轉速。隨著計算流體力學(CFD)在流體機械穩態數值模擬中的廣泛運用[6,7]，應用CFD技術解決水利工程問題愈發成熟。國內已有學者嘗試利用CFD方法對泵的飛逸特性進行研究[8-12]，而其中，對包含運動邊界的非定常流動進行數值模擬是三維瞬態研究的難點。

本文基于有限體積法的動網格技術對大套一站軸流泵裝置的飛逸狀態進行了三維瞬態數值模擬研究，監測了不同凈揚程下葉片壓力、軸向力等參數變化。同時，與模型試驗結果比較，顯示了很好的可靠性。

1 計算對象

1.1 模型參數

本文計算模型基于大套一站立式軸流泵裝置，停機方式為真空破壞閥開啟引入空氣進入虹吸式出水流道斷流，懸掛式電動機直接傳動，其具體參數如表1所示。模型包含進水池、肘型進水流道、葉輪區、導葉區、虹吸式出水流道、出水池等部件，結構如圖1所示。數值計算模型與軸流泵裝置實際尺寸比例為1∶1，葉片安裝角度0度。

表1 具有虹吸式出水流道軸流泵裝置參數Tab.1 Parameters of axial pump with siphon oulet

圖1 軸流泵裝置幾何模型Fig.1 Geometrical model of axial pump

1.2 網格劃分

采用非結構化網格來劃分進水流道、葉輪、導葉及出水流道；采用結構化網格來劃分進水池、出水池；由于葉輪區和導葉區流態復雜，對其進行網格加密。經網格無關性驗證計算，發現網格超一定數量后對裝置性能影響很小，最終選擇方案3來劃分計算模型，網格總數為245 萬個，不同網格劃分方案如表2所示。

表2 不同網格劃分方案Tab.2 Results of different meshing programs

2 數值計算方法

2.1 控制方程與湍流模型

本次數值模擬屬于三維非定常不可壓縮湍流流動，對非穩態Navier-Stokes方程采用時間平均法，得到時均形式的控制方程[13]。

連續方程：

(1)

動量方程：

(2)

式中：ui，uj表示流體速度分量。

本文在對軸流泵裝置飛逸狀態進行數值模擬時，轉速的數值變化較大，過程中葉輪及導葉流體區域的流態較為復雜。基于此，本文采用基于k-ε方程改進而來的Realizablek-ε兩方程模型為本文三維數值模擬的湍流模型[14]。

2.2 葉輪及動網格控制

本文利用動網格技術與UDF自定義的方法實現真空破壞閥未開啟出水流道內沒有空氣進入條件下葉輪飛逸過程的轉速變化。為使網格能夠適應運動邊界的移動和幾何形狀的變化，必須要對計算網格進行修正。邊界移動的任意控制體積V一般標量 的守恒方程為[15]：

(3)

同時利用UDF自定義技術通過葉輪力矩平衡方程控制轉速變化，方程如下[16]：

(4)

式中：J為泵裝置轉動慣量，kg/m2；ω為飛逸過程中葉輪角速度，rad/s；M0為電機電磁力矩，N·m，在飛逸過程中，電動機斷電，M0=0；M1為水泵的水力矩，N·m，其由UDF功能實時讀取葉片上轉矩得到；M2為軸承摩擦力矩，N·m，僅考慮推力軸承摩擦力矩，省略徑向摩擦力矩；M3為電機風損力矩，N·m，較小省略。

2.3 離散格式及定解條件

離散格式：本次數值計算利用Fluent 6.3軟件完成，用有限體積法對上述數學模型進行離散，壓力項采用PRESTO格式，體積分數項采用Geo-Reconstruct格式，湍動能和對流項采用一階迎風格式，采用適合瞬態計算的PISO算法對流場速度壓力進行求解，數值計算迭代時間步長為0.002 s，初始時間為0 s，總計算時長60 s。

定解條件：進水池水面采用壓力進口條件，壓力值由進水池水位確定；出水池水面采用壓力出口條件，壓力值由出水池水位決定；

初始時刻，水流流速為0，葉輪轉速為0。轉速由葉輪力矩方程控制，直至到達飛逸轉速并保持穩定。

3 數值計算結果

圖2為軸流泵裝置達到飛逸狀態時流到子午剖面的速度矢量圖。由圖2可見裝置達到飛逸狀態時，裝置流道內水流為反向流動，水流從出水池通過出水流道、導葉、葉輪、進水流道等部件倒流至進水池且流線平順。流道內壓力分布圖顯示，隨水平高度的上升，水流壓力值逐漸縮小，飛逸狀態下虹吸式出水流道虹頂處處于負壓狀態。

圖2 事故飛逸狀態速度矢量圖Fig.2 Velocity vectors in runaway state

圖3為軸流泵裝置從靜止到飛逸狀態的葉輪轉速變化過程曲線，分別計算了0.5、1.0、2.0、3.0、4.0、5.4、6.8 m總共7個不同凈揚程工況。由圖可見，60 s時，各個凈揚程下葉輪轉速均已達到反向最大，且保持穩定，此時軸流泵裝置已進入飛逸狀態。隨著裝置凈揚程的升高，回流水流獲得的能量及其回流速度隨之增大，從而飛逸狀態下葉輪的飛逸轉速隨凈揚程升高而增大。6.8 m時飛逸轉速為348.8 r/min，為額定轉速214.3 r/min的1.62倍，此時機組運行會產生安全隱患，因而在高揚程下運行時要采取預防措施，防止進入飛逸狀態。

圖3 不同凈揚程飛逸轉速變化曲線Fig.3 Changing curves of runaway speed of various head

飛逸狀態時，軸向力的變化是十分關鍵的參數，本次數值模擬中不同凈揚程，飛逸狀態下的葉輪軸向力如表3所示，葉輪軸向力的方向為垂直向下。如表3所示，隨凈揚程增加，軸向力逐步增大，這是因為裝置處于飛逸工況與水輪機工況相似，較高的揚程導致了葉輪葉片上較高的轉矩和軸向力。數值模擬中最高凈揚程6.8 m下軸向力46.48 kN為轉動部件重量83 kN的0.56倍。

表3 不同凈揚程飛逸狀態葉輪軸向力Tab.3 Axial force in runaway state of various head

圖4 不同凈揚程飛逸狀態下葉片壓力分布Fig.4 Pressure distribution on blades in runaway state of various head

不同揚程下葉片表面壓力分布如圖4所示，左側為葉片壓力面，右側為葉片吸力面。0.5 m揚程下，由于揚程數值較小，葉片壓力面和吸力面的壓力較小；3.0 m揚程時，可明顯發現，葉片壓力面壓力從進水邊至出水邊逐漸增大，吸力面壓力從進水邊至出水邊逐漸減小，壓力梯度明顯；6.8 m揚程下，葉片壓力分布的變化規律和3.0 m揚程時相似，但局部高壓及低壓區域均增大；高壓及低壓區域的分布顯示，吸力面進水邊為水流撞擊，壓力面進水邊為水流脫流；隨著揚程的升高，撞擊及脫流現象逐漸加重，葉片表面進出水邊的壓力差亦逐漸增大。

4 計算結果驗證

圖5 為三維數值模擬與模型試驗在0度葉片安放角下不同揚程所對應的飛逸轉速的對比圖，飛逸轉速曲線中數值模擬部分由圖3得到。由圖可見數值模擬與模型試驗的飛逸轉速延揚程變化趨勢是一致的，數值相差較小(最大誤差為5%)。數值模擬中的飛逸轉速略高于模型試驗值，一方面是由于本文中對葉輪的力矩計算中僅考慮了軸向推力軸承的摩擦力矩，而徑向軸承摩擦力矩及轉子風阻力矩均近似省略，另一方面，由于模型試驗中原型泵飛逸轉速曲線是通過單位轉速換算得到的，與實際數值有一定的偏差。因而，利用本文所用模型及數值計算方法對軸流泵裝置的起動及停及過渡過程進行數值模擬是可行的。

圖5 數值模擬與試驗數據對比圖Fig.5 Comparison of experiment and numerical simulation

5 結 語

(1)利用動網格及UDF自定義技術可對軸流泵裝置飛逸工況進行三維瞬態數值模擬，可直觀獲得流道內的動態特性。

(2)隨裝置凈揚程的升高，其對應飛逸工況下，飛逸轉速，葉片軸向力均隨之增大，且高揚程下飛逸轉速較大，影響機組運行安全。

(3)數值模擬結果與模型試驗比較顯示：數值模擬下飛逸轉速與模型試驗結果吻合較好，其結果具有很好的準確性。

□

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

[2] 關醒凡, 聶書彬, 黃道見,等. 軸流泵飛逸轉速特性試驗研究[J]. 水泵技術, 2002,(5):14-15,28.

[3] 劉潤根, 馬曉忠, 詹 磊. 黃家壩30°斜式軸流泵裝置模型試驗研究[J]. 中國農村水利水電, 2016,(2):109-111.

[4] ZHANG De-sheng, SHI Wei-dong, CHEN Bin, et al. Unsteady flow analysis and experimental investigation of axial-flow pump[J]. Journal of Hydrodynamics, 2010,22(1):35-43.

[5] 沙 毅, 宋德玉, 段福斌, 等. 軸流泵變轉速性能試驗及內部流場數值計算[J]. 機械工程學報, 2012,(6):187-192.

[6] SHI Wei-dong, ZHANG De-sheng, Guan Xing fan, et al. Numerical and experimental investigation on high-efficiency axial-flow pump[J]. Chinese Journal of Mechanical Engineering, 2010,23(1):38-44.

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

[8] Liu S, Zhou D, Liu D, et al． Runaway transient simulation of a model Kaplan turbine[C]//Proceedings of the IOP Conference Series: Earth and Environmental Science, 2010.

[9] 周大慶, 鐘淋涓, 鄭 源, 等. 軸流泵裝置模型斷電飛逸過程三維湍流數值模擬[J]. 排灌機械工程學報, 2012,(4):401-406.

[10] 金國棟, 潘志軍, 孟金波, 等. 斜式軸流泵裝置模型的飛逸特性研究[J]. 水動力學研究與進展A輯, 2013,(5):591-596.

[11] 周大慶, 張任田, 屈 波, 等. 大型立式軸流泵站停泵過渡過程研究[J]. 河海大學學報(自然科學版), 2006,34(3):272-275.

[12] ZHANG Fan, YUAN Shou-qi, FU Qiang, et al. Investigation on transient flow of a centrifugal charging pump in the process of high pressure safety injection[J]. Nuclear Engineering and Design, 2015,293:119-126.

[13] H K Versteeg, W Malalasekera. An introduction to computational fluid dynamics: the finite volume method[M]. New York: Wiley,1995.

[14] 王福軍. 計算流體動力學分析——CFD 軟件原理與應用[M]. 北京：清華大學出版社,2004:118-119.

[15] 黃 思, 郭 京, 張建可,等. 運用三維動網格技術模擬計算旋噴泵的非定常流動[J]. 中國農村水利水電, 2014,(4):112-115.

[16] 陳松山, 蔣紅梅, 周正富, 等. 大型貫流泵站機組啟動過渡過程仿真計算[J]. 揚州大學學報(自然科學版), 2009,12(3):74-78.