可逆式水輪機過渡過程水力計算方法研究

黃　祺，聶雄杰，錢　隆，王得蛟

(1. 蘭州理工大學 技術工程學院, 甘肅 蘭州 730050；2.湖南凱利特泵業(yè)有限公司, 湖南 湘潭 411202)

可逆式水輪機過渡過程水力計算方法研究

黃祺1，聶雄杰2，錢隆1，王得蛟1

(1. 蘭州理工大學 技術工程學院, 甘肅 蘭州 730050；2.湖南凱利特泵業(yè)有限公司, 湖南 湘潭 411202)

摘要：在可逆式水輪機復雜的過渡過程中，壓力突變一直沒有準確的算法能夠描述，雖然水擊特征線方程和連續(xù)方程能夠大致計算出過渡過程中的壓力變化，但存在計算量很大的問題。利用FLUENT軟件并改變轉輪的邊界條件，基于雷諾時均方程和RNG k-ε湍流模型, 將壓力和速度耦合，采用半隱式(SIMPLEC)算法模擬分析了可逆機過渡過程中不同工況點的水力變化特點, 得出了可逆式水輪機甩負荷過程中的水力特性,其與理論公式計算結果極為相近，為可逆機過渡過程的水力計算提供了一條新的途徑。

關鍵詞：可逆式水輪機；過渡過程； SIMPLEC算法；水力特性

可逆式水輪機組作為一種行之有效的抽水蓄能裝置，具有電設備起停快速和調(diào)節(jié)靈活的優(yōu)點，能有效地應付負荷的變化[1]。與常規(guī)水電機組相比，其具有工況多樣和變化復雜的特點，所涉及到的可逆機組過渡過程的水力性能預測更是具有非常重要的作用，本文利用理論和數(shù)值模擬計算2種方法分析了水輪機工況甩負荷過程中的水力變化，為可逆式機組其他過渡過程的水力計算提供了一定的理論和試驗依據(jù)。

1數(shù)學模型

1.1水擊過程的運動方程和連續(xù)性方程

水擊過程的運動方程和連續(xù)性方程分別為：

(1)

(2)

式中，a是水擊波速度，單位為m/s；v是水流速度，單位為m/s；D是管道直徑，單位為m；f是摩擦因數(shù)；α是管道中心線和水平線的夾角；H是測壓管壓力;x是水體微元移動距離；g是重力加速度；A是管道斷面面積。

同時：

將式1和式2偏微分方程變?yōu)橐詔、H和Q為未知量的常微分方程：

(3)

(4)

(5)

(6)

計算時引入相應的邊界條件就可以求得過渡過程的瞬態(tài)值。

1.2可逆式水輪機運動方程

可逆式水輪機組的運動方程為：

(7)

式中，M/M0是轉矩相對值，可用m表示；P/P0是功率相對值；Ta是機組時間常數(shù)。

全甩負荷以后，機組的運動方程簡化為：

(8)

1.3水輪機特征方程

水輪機的特性用下述方程來表示：

(9)

Q11=f(n11,a)

(10)

(11)

M11=f(n11,a)

(12)

(13)

Ht=Hi-Hi+1

(14)

式中，n是機組轉速，單位為r/min;n11、Q11、M11是水輪機的單位轉速、流量和力矩；Hi和Hi+1是水輪機機組前、后節(jié)點水壓力。

2理論計算實例

某可逆式水輪機主要參數(shù)為：可逆機工況正常壓力水頭為204 m,設計引水流量為2.35 m3/s,單機容量為4 MW,額定轉速為500 r/min，轉輪直徑為0.45 m。在導葉直線關閉的不同時間內(nèi)，計算時引入相應的邊界條件，利用式3～式6就可以計算出可逆機水輪機工況甩負荷的過程中機組的水力變化[2-4]，計算結果見表1，機組的其他水力過渡工況具有類似的計算過程。

表1　計算可逆機過渡過程數(shù)據(jù)

3FLUENT數(shù)值模擬計算實例

數(shù)值模擬計算域點如圖1所示，計算時在蝸殼內(nèi)部設置9個檢測點位置，如圖2所示，尾水管內(nèi)部設置5個檢測點位置，如圖3所示，邊界條件[5-8]采用雷諾時均方程和RNG k-ε湍流模型，壓力和速度耦合采用SIMPLEC算法。由于可逆機過渡過程需要不同的邊界條件，水輪機甩負荷工況采用壓力進口、壓力出口和轉速漸變的邊界條件。壁面條件滿足無滑移固壁條件，所以采用標準壁面函數(shù)法。計算結果采用最大值算法統(tǒng)計最終數(shù)值，統(tǒng)計結果見表2[9-10]。

圖1　數(shù)值模型計算域　　　　圖2　蝸殼內(nèi)部監(jiān)測點

圖3　尾水管內(nèi)部監(jiān)測點

關閉時間/s機組轉速/r·min-1蝸殼最大壓力/mH2O尾水管最小水壓/mH2O3664247-1.64677234-1.85687.5223-1.56698.5222-1.17706.5219-1.4

4計算結果對比分析

通過對可逆機過渡過程中的理論計算和數(shù)值模擬計算結果(見表1和表2)進行對比分析，從圖4和圖5可以看出，過渡過程中理論計算結果和數(shù)值模擬計算結果較為接近，相對于理論計算而言，數(shù)值模擬計算的讀取數(shù)據(jù)略大，出現(xiàn)這一問題的原因可能與軟件的運行機理有一定關系，所讀取數(shù)據(jù)的精度已經(jīng)能夠達到預測過渡過程壓力突變的目的，可以初步推斷在其他條件下，F(xiàn)LUENT數(shù)值計算過渡過程也有一定的可行性。同時，在導葉直線關閉過程中，機組轉速和關閉時間成明顯的正比關系，而蝸殼中最大壓力變化則與時間成反比關系，相對而言，尾水管中的壓力變化較為平緩，為了平衡蝸殼中最大壓力與機組轉速，本文選定導葉直線關閉時間為5 s，進行了甩負荷的數(shù)值模擬，結果如圖6所示。

圖4　不同機組轉速下理論計算(P1)和　數(shù)值模擬(P2)蝸殼壓力分布

圖5　不同機組轉速下理論計算(P1)和　數(shù)值模擬(P2)尾水管壓力分布

圖6　導葉直線關閉時間為5 s時蝸殼壓力　和機組轉速分布

從圖6可以看出，在導葉直線關閉時間為5 s時，機組的轉速和蝸殼壓力在3 s左右即出現(xiàn)最大值，另外從數(shù)值計算的其他數(shù)據(jù)觀測，導葉關閉時間過快會在引流道內(nèi)出現(xiàn)較大的水擊現(xiàn)象，關于水擊和導葉關閉時間的最優(yōu)調(diào)整還有待進一步分析。

5結語

可逆機甩負荷過程的壓力突變可以用數(shù)值模擬計算的方法取得定量的結果，計算結果較為理想。可逆機甩負荷過程中尾水管的壓力變化波動較小，相對蝸殼劇烈的壓力變化而言可以忽略不計。選定導葉直線關閉時間為5 s進行了甩負荷的數(shù)值模擬，結果顯示在甩負荷后3 s左右達到壓力和轉速的峰值。

參考文獻

[1] 梅祖彥．抽水蓄能發(fā)電技術[M]．北京：機械工業(yè)出版社，2000.

[2] 黃賢榮.水電站過渡過程計算中的若干個問題研究[D].南京：河海大學，2006.

[3] 鮑海艷，楊建東，付亮.水電站功率調(diào)節(jié)過渡過程數(shù)值模擬[J] 武漢大學學報:工學版, 2009,42(1):56-59.

[4] 李建新, 孫美風, 岳曉娜,等.基于MATLAB的水力機組過渡過程研究[J].廣西水利水電，2009(1)：77-80.

[5] 王福軍，黎耀軍，王文娥，等. 水泵CFD應用中的若干問題與思考[J]．排灌機械，2005，23(5)：1-10.

[6] 李琪飛,李仁年,韓偉,等.混流式水輪機引水、導水部件內(nèi)部固液兩相流動的數(shù)值分析[J].蘭州理工大學學報,2008,34(6):47-50.

[7] Martin S C. The hydraulic machinery and systems 20thIAHR symposium[C]. IAHR, Charlotte, North Carolina, USA, 2000.

[8] 李琪飛，李仁年，敏政，等．水輪機PID調(diào)速器最佳參數(shù)整定及尋最優(yōu)計算方法[J].排灌機械，2006，24(4)：33-36.

[9] Hellstroem J G I, Marjavaara B D, Lundstroem T S. Parallel CFD simulations of an original and redesigned hydraulic turbine draft tube[J]. Advances in Engineering Software, 2007,38(5):338-344.

[10] 王瑞金，張凱，王剛.FLUENT技術基礎與應用教程[M]．北京：清華大學出版社，2007.

責任編輯李思文

Calculation Method of Hydraulic Transient for Pump Turbine

HUANG Qi1, NIE Xiongjie2, QIAN Long1，WANG Dejiao1

(1.College of Technology and Engineering, Lanzhou University of Technology, Lanzhou 730050，China;2.Hunan Credo Pump Co., Ltd，Xiangtan 411202,China)

Abstract:There has been no accurate algorithm can describe the pressure jump of the transition process for pump turbine, water hammer characteristic equation and continuous equation to calculate the transition process of roughly pressure changes, but with large amount of calculation. Using FLUENT software and changing the boundary conditions of the runner , based on Reynolds time-averaged equations and RNG k-εturbulence model and use SIMPLEC algorithm for the pressure and velocity coupling. The simulation analysis in the process of pump turbine the hydraulic characteristics of different operation points, concluded that the turbine′s hydraulic characteristics in the load rejection transient, by comparison the result of numerical simulation and theoretical formula is very close, it opened a new way for the transition process of the pump turbine hydraulic calculation.

Key words:pump turbine, transient process, SIMPLEC algorithm, hydraulic characteristics

收稿日期：2014-3-20

作者簡介：黃祺(1984-)，男，碩士，講師，主要從事流體機械等方面的研究。

中圖分類號：TK 05

文獻標志碼：A