電站鍋爐調(diào)節(jié)閥管配系統(tǒng)流場特性的數(shù)值分析

劉曉敘，劉 佳，何慶中，郭 斌

四川理工學院機械工程學院，四川自貢 643000

【環(huán)境與能源 / Environment and Energy】

電站鍋爐調(diào)節(jié)閥管配系統(tǒng)流場特性的數(shù)值分析

劉曉敘，劉 佳，何慶中，郭 斌

四川理工學院機械工程學院，四川自貢 643000

分析某火電機組儲水罐調(diào)節(jié)閥在運行中出現(xiàn)的閥后配管端蓋損壞問題，建立該調(diào)節(jié)閥和配管的計算流體動力學模型，利用CFD軟件CFX對其進行兩相流流場特性分析．結(jié)果表明，調(diào)節(jié)閥在節(jié)流口處發(fā)生了強烈的汽化現(xiàn)象，導致配管底板處沖蝕力較大，汽蝕沖擊現(xiàn)象嚴重，分析結(jié)果與實際破壞情況相符．結(jié)合分析改進了調(diào)節(jié)閥配管系統(tǒng)結(jié)構(gòu)，并對改進后的結(jié)構(gòu)進行流場數(shù)值模擬，結(jié)果表明，在配管底板處流場速度降低近70%，汽相體積減少90%以上，汽蝕沖擊現(xiàn)象減弱，分析結(jié)果可供該類閥門和配管設計參考．

計算流體力學；調(diào)節(jié)閥；兩相流；CFX軟件；汽蝕；結(jié)構(gòu)優(yōu)化

超臨界火電機組儲水罐調(diào)節(jié)閥主要配套過熱系統(tǒng)使用，用于控制儲水罐內(nèi)的壓力和水位．在使用過程中，高溫高壓介質(zhì)經(jīng)調(diào)節(jié)閥減壓后，產(chǎn)生強烈的汽化現(xiàn)象，并伴隨劇烈的振動和噪聲，給閥門和管配系統(tǒng)造成嚴重的沖擊，影響整個閥門及管配系統(tǒng)的使用性能和壽命，甚至危害整個鍋爐系統(tǒng)的運行安全．因此，有必要對該類閥門及管配系統(tǒng)進行流場特性研究和結(jié)構(gòu)優(yōu)化．

目前，對相關(guān)流場的研究主要集中在離心泵[1-2]、軸流泵[3]、常規(guī)調(diào)節(jié)閥[4]、球閥[5-6]以及相關(guān)水利機械上，而對高溫高壓條件下的調(diào)節(jié)閥[7-8]管配系統(tǒng)流場分析還較少．因此，本研究基于重整化群 (renormalization group，RNG)k-ε湍流模型和混合物均相流空化模型，對某火電機組儲水罐調(diào)節(jié)閥進行流場特性模擬．給出了設計工況下閥門和管配的流場分布，結(jié)合流場分析時發(fā)現(xiàn)的問題，提出相應改進方案，并對改進后的結(jié)構(gòu)進行數(shù)值模擬，研究結(jié)果對閥門和管配結(jié)構(gòu)優(yōu)化設計、提高可靠性及其國產(chǎn)化有積極意義．

1 破壞原因分析及模型建立

閥門流動介質(zhì)為高溫高壓水，溫度達到311 ℃，壓力為10.2 MPa，此時，介質(zhì)處于壓力和溫度動平衡的臨界狀態(tài)，細微的參數(shù)變化就會引起介質(zhì)屬性的轉(zhuǎn)變．根據(jù)閥門運行工況要求，調(diào)節(jié)閥主要技術(shù)參數(shù)設計如下：壓力為30.5 MPa、溫度375 ℃、運行條件為直流升負荷、開度21%、進口壓力10.2 MPa、出口壓力4.78 MPa、閥門入口溫度311 ℃．

圖1 儲水罐調(diào)節(jié)閥Fig.1 Regulating valve of water storage

圖1為該調(diào)節(jié)閥實物圖．在電站經(jīng)3次啟動試運行后，發(fā)現(xiàn)閥后直管末端直角轉(zhuǎn)彎處端蓋被擊穿．此時閥門開度為21%，處于升負荷至直流升負荷階段．在機組啟動過程進入直流運行前，閥門入口介質(zhì)(來自儲水罐)為飽和水，通過調(diào)節(jié)閥節(jié)流后的介質(zhì)為汽水兩相混合物．若閥后管路布置不當，閥門出口汽水兩相介質(zhì)在閥后直管段未形成穩(wěn)定流態(tài)前進入90°直角轉(zhuǎn)彎處，流態(tài)受管路原因急劇改變將在端蓋底部形成湍流強度較大的漩渦，介質(zhì)高頻沖擊將造成端蓋的破壞．

本研究根據(jù)某企業(yè)在實際運行中發(fā)生端蓋破壞的調(diào)節(jié)閥尺寸，建立流場分析模型，如圖2．

圖2 調(diào)節(jié)閥結(jié)構(gòu)簡圖Fig.2 Structure sketch of regulating valve

利用HyperMesh對調(diào)節(jié)閥流道模型進行網(wǎng)格劃分，以非結(jié)構(gòu)化網(wǎng)格為主．由于閥芯為關(guān)鍵部位，其周邊流場較為復雜，劃分網(wǎng)格時對閥芯進行了局部細化并考慮了各塊間連接區(qū)域網(wǎng)格具有的連續(xù)性．最終得到網(wǎng)格數(shù)為2 941 798，節(jié)點數(shù)為580 096．經(jīng)過網(wǎng)格無關(guān)性檢查，在此情況下，網(wǎng)格數(shù)量增加并不會進一步提高計算結(jié)果的準確性，故網(wǎng)格劃分是可信的．流道網(wǎng)格劃分如圖3．

圖3 流道網(wǎng)格劃分示意圖Fig.3 Grid of flow field model

2 空化湍流控制理論

本研究采用RNG k-ε湍流模型，以及基于Reyleigh-Plesset方程的混合物均相流空化模型來求解調(diào)節(jié)閥復雜的空化流動．對于調(diào)節(jié)閥流道內(nèi)的空化流動，可認為是空泡相和水流相達到了熱平衡和動力平衡[9]，將空泡相和水流相當成統(tǒng)一的流體進行分析，因此采用混合物均相流模型對調(diào)節(jié)閥的空化流動進行流場分析．

2.1 連續(xù)性方程

(1)

2.2 動量方程

(2)

2.1 Reyleigh-Plesset方程

(3)

考慮到流場處于低頻狀態(tài)和忽略表面張力的情況下，可將式(3)簡化為

其中，RB為氣泡半徑； pV為局部溫度下的汽化壓力； ρf為液體密度； σ為表面張力．

3 流場計算結(jié)果及分析

將網(wǎng)格模型導入CFX中，按實際使用工況參數(shù)設置邊界條件，進口壓力10.2 MPa，出口壓力4.78 MPa．給定流體介質(zhì)溫度311 ℃，飽和蒸汽壓約為9.9 MPa．此時，流動介質(zhì)屬于低于350 ℃的飽和水，此溫度下介質(zhì)處于臨界飽和線附近，當溫度不變時，微小的壓力波動就會引起介質(zhì)狀態(tài)的改變．閥門數(shù)值模擬選用高階求解模式，物理時間尺度設置為0.1 s，收斂最大殘差值為1.0×10-4．計算收斂圖顯示收斂效果很好，說明所建流道模型合理．

3.1 流場分析結(jié)果

根據(jù)流場數(shù)值計算結(jié)果顯示，在4 s后介質(zhì)流動狀態(tài)趨于穩(wěn)定，此時的流場分布狀況如圖4和圖5．圖4和圖5為截取調(diào)節(jié)閥中心對稱面流場分布情況．從圖4速度云圖可知，介質(zhì)經(jīng)過調(diào)節(jié)閥閥芯喉口時速度急劇增大，最高可達149 m/s，而高速介質(zhì)隨著管道流動，能量有所衰減，在調(diào)節(jié)閥直角轉(zhuǎn)彎處速度達到60 m/s左右，同時可見在端蓋底板附近形成了一定湍流強度的急劇漩渦，頻繁撞擊端蓋底板．從圖5汽相體積云圖可見，介質(zhì)經(jīng)過喉口時，發(fā)生大量汽化且最先發(fā)生在壁面附近．隨著管道的流動，能量有所衰減，部分介質(zhì)壓力恢復至空化壓力之上，使得這部分氣泡壓縮潰滅，在氣泡潰滅過程中，散發(fā)出巨大的能量，該能量會在短時間內(nèi)反復沖擊閥體壁面．尤其是在直管末端，由于底板的阻擋，氣泡沖擊底板而破裂，因此，圖5中顯示在底板附近汽相介質(zhì)含量很低．

圖4 調(diào)節(jié)閥速度云圖Fig.4 Contours of the velocity

圖5 調(diào)節(jié)閥汽相體積分布Fig.5 Contours of the vapor volume fraction

為更加直觀表現(xiàn)調(diào)節(jié)閥的調(diào)節(jié)過程，本研究提取了喉口和底板處樣點的速度和汽相體積分數(shù)變化數(shù)據(jù)并繪制成曲線，如圖6和圖7．

圖6 調(diào)節(jié)閥樣點速度曲線Fig.6 Sample velocity curve of the valve

圖7 調(diào)節(jié)閥樣點汽相體積曲線Fig.7 Sample vapor volume fraction curve of the valve as a function of time

結(jié)合圖6和圖7曲線可看出，喉口附近介質(zhì)汽化體積分數(shù)達到90%以上，減壓效果明顯，同時，壓力內(nèi)能轉(zhuǎn)化為動能使得速度急劇上升．而在直管末端底板附近，介質(zhì)流動狀態(tài)極不穩(wěn)定，波動的速度頻繁對底板形成沖擊，這也驗證了實際使用過程中底板的失效是由于介質(zhì)高速沖擊造成的．

3.2 閥后直管段結(jié)構(gòu)改進及分析結(jié)果

綜上分析可知，在閥后直管末端底板處存在較強的汽蝕沖擊現(xiàn)象，其原因在于介質(zhì)流過喉口時劇烈的壓降導致速度急劇增加，介質(zhì)汽化嚴重，由于管道長度不夠，介質(zhì)速度衰減不充分，高速運動的介質(zhì)對底板造成強烈的沖擊，最終導致底板被擊穿．結(jié)合分析結(jié)果，為減少汽蝕的影響，在滿足工程整體布局的前提下，增加閥后配管的長度并進行模擬計算．設置相同的邊界條件和介質(zhì)參數(shù)，分別增長管配200～1 000 mm等5種方案的模擬計算．下面選取增長1 000 mm的模擬結(jié)果，如圖8、圖9和圖10．

圖8 改進后調(diào)節(jié)閥速度云圖Fig.8 Contours of the velocity after optimizing

圖9 改進后調(diào)節(jié)閥汽相體積分布Fig.9 Vapor volume fraction distribution of the valve after optimizing

圖10 改進后速度與汽相體積對比Fig.10 Curves of velocity and vapor volume fraction after optimizing

從圖8可見，介質(zhì)經(jīng)過喉口后速度仍然較高，接近132 m/s，而隨著配管的增長，速度逐漸衰減，在底板附近降至20 m/s以下，符合常規(guī)閥門流速的要求．從圖9可知，隨著管道的延長，底板附近汽相體積幾乎為0．結(jié)合圖10延長后底板樣點速度與汽相體積曲線圖可以看出，在1.5 s后，底板附近介質(zhì)流速降低，汽相體積幾乎為0，說明增長管道后底板附近流場改善比較明顯，該區(qū)域處于相對安全的狀態(tài)．

為深入探究閥后管配長度對整體流場的影響，分別模擬不同管配長度下流場分布情況，并提取了相應管配長度下底板附近樣點的流場量變化，繪制成速度和汽相體積變化對比曲線，如圖11和圖12．

圖11 不同增長量下底板處速度變化Fig.11 Velocity distribution at the baseboard under different lengthen pipelines

圖12 不同增長量下底板處汽相體積變化Fig.12 Vapor volume fraction distribution at the baseboard under different lengthen pipelines

從圖11可見，增加管配流道長度后，流場速度下降較為明顯，在加長至800 mm后，速度降低至18 m/s以下，與圖6相比速度降低近70%，隨著時間推移還有進一步降低的趨勢，表明加長管配管道可以有效改善底板處的流場分布情況．

此外，從圖12可見，增加管配管道長度后，底板附近的汽相體積占比減少，與圖7相比汽相體積減少近95%．在管道加長800 mm后，該區(qū)域的汽相體積占比幾乎為0，說明在該種情況下底板附近流場幾乎為液相，即流場分布較為均勻，表明閥門結(jié)構(gòu)滿足了運行工況的要求．

4 結(jié) 論

本研究結(jié)合某超臨界火電機組調(diào)節(jié)閥在實際運行中出現(xiàn)的問題，基于CFX空化湍流模型，對其進行了流場分析，得到了該調(diào)節(jié)閥在設計工況下的流場規(guī)律，并通過對其進行結(jié)構(gòu)改進，得出以下結(jié)論：

1)通過對調(diào)節(jié)閥和管配系統(tǒng)的流場特性分析，其結(jié)果和實際閥后直角轉(zhuǎn)彎處端蓋破壞情況比較吻合，表明使用CFX分析調(diào)節(jié)閥流場的方法是合理的．

2)通過結(jié)構(gòu)改進后，底板附近流場速度降低近70%，氣泡潰滅沖蝕現(xiàn)象明顯減弱，汽相體積減少近90%，流動相對穩(wěn)定，有效改善了底板處汽蝕問題．

3)管配管道的長度與流場分布的均勻性有較大的關(guān)系，管道越長，流場相對越穩(wěn)定，體現(xiàn)了多相流流動的松弛性．

4)改進后的結(jié)構(gòu)能有效改善配管銜接處的流場特性，降低介質(zhì)對配管的沖擊，提高配管的使用壽命．對今后同類閥門配管布置具有重要的參考價值．

/ References：

[1] 劉厚林,劉東喜,王 勇,等.三種空化模型在離心泵空化流計算中的應用評價[J].農(nóng)業(yè)工程學報,2012，28(16):54-59. Liu Houlin, Liu Dongxi, Wang Yong, et al. Applicative evaluation of three cavitating models on cavitating flow calculation in centrifugal pump[J]. Transactions of the Chinese Society of Agricultural Engineering, 2012，28(16):54-59.(in Chinese)

[2] 王秀禮，袁壽其，朱榮生,等.離心泵汽蝕穩(wěn)定流動特性數(shù)值模擬[J].農(nóng)業(yè)機械學報，2012,43(3)：67-72. Wang Xiuli, Yuan Shouqi, Zhu Rongsheng, et al. Numerical simulation on cavitation unsteady characteristics in centrifugal pump[J]. Transactions of the Chinese Society for Agricultural Machinery,2012,43(3):67-72.(in Chinese)

[3] 楊敏官，姬 凱，李 忠.軸流泵葉輪內(nèi)空化流動的數(shù)值計算[J].農(nóng)業(yè)機械學報，2010,41:10-14. Yang Minguan, Ji Kai, Li Zhong. Numerical calculation of cavitating flow in impeller of axial-flow pump[J]. Transactions of the Chinese Society for Agricultural Machinery, 2010, 41: 10-14.(in Chinese)

[4] 吳國洋，何澤銀，李國云，等. 調(diào)壓閥內(nèi)流道流場分析及閥芯結(jié)構(gòu)改進[J].蘭州理工大學學報,2013，39(3):38-42. Wu Guoyang, He Zeyin, Li Guoyun, et al. Flow field analysis of internal flow passage in presure-regulating valve and spool structure improvement[J]. Journal of Lanzhou University of Technology, 2013, 39(3):38-42.(in Chinese)

[5] 張生昌，張玉林，方志明，等. 基于CFD 的球閥三維流場數(shù)值模擬[J].排灌機械工程學報，2013,31(2)：157-161. Zhang Shengchang, Zhang Yulin, Fang Zhiming,et al.3D flow numerical simulation of ball valve based on CFD[J]. Journal of Drainage and Irrigation Machinery Engineering,2013,31(2)：157-161.(in Chinese)

[6] 屈 鐸，彭利坤，林俊興，等. 基于CFD的球閥閥腔流場的數(shù)值仿真[J].液壓與氣動，2014(11)：112-115. Qu Duo, Peng Likun, Lin Junxing, et al. Numerical simulation of flow field of the ball valve cavity based on CFD[J]. Chinese Hydraulics & Pneumatics, 2014(11)：112-115.(in Chinese)

[7] 何慶中,郭 斌,陳雪峰,等.超超臨界電站鍋爐調(diào)節(jié)閥及管配系統(tǒng)空蝕特性研究[J].制造業(yè)自動化，2015,8(15):92-95. He Qingzhong, Guo Bin, Chen Xuefeng, et al. Cavitation research of regulating valve and pipe system used in ultra-supercritical Boilers[J]. Manufacturing Automation, 2015, 8(15): 92-95.(in Chinese)

[8] 郭 斌，何慶中，王渝浩，等.超超臨界火電機組高溫減壓調(diào)節(jié)閥空蝕問題分析及結(jié)構(gòu)優(yōu)化[J].熱力發(fā)電，2015,10(10):82-85. Guo Bin, He Qingzhong, Wang Yuhao, et al. Cavitation analysis and structural improvement of regulating?valve 361 used in the system of thermal ultra supercritical[J]. Thermal Power Generation, 2015,10(10):82-85.(in Chinese)

[9] Wang Chichuan, Tseng C Y, Chen I Y. A new correlation and the review of two-phase flow pressure change across sudden expansion in small channels[J].International Journal of Heat and Mass Transfer， 2010, 53(19/20)：4287-4295．

[10] 洪 鋒，袁建平，周幫倫.空泡半徑非線性變化的空化模型及應用[J].華中科技大學學報自然科學版,2015,43(10)：15-20. Hong Feng, Yuan Jianping, Zhou Banglun. Cavitation model considering non-linear variation of bubble radius and its application[J]. Journal of Huazhong University of Science and Technology Nature Science Edition, 2015,43(10)：15-20.(in Chinese)

[11] Chern M J, Hsu P H, Cheng Y J, et al. Numerical study on cavitation occurrence in globe valve[J]. Journal of Energy Engineering, 2013,139(3):25-34.

[12] Kwon W C，Kim G R，Park S C，et al. Design of a tortuous path trim for a high-pressure turbine bypass valve[J]. Proceedings of the Institution of Mechanical Engineers，Part E: Journal of Process Mechanical Engineers, 2010,224(E2):149-153.

[13] 李樹勛，胡建華，李連翠，等.超(超)臨界蒸汽疏水閥內(nèi)部空化模擬研究[J].四川大學學報工程科學版，2013,45(5):145-153. Li Shuxun, Hu Jianhua, Li Liancui, et al. Numerical study of cavitation flows in ultra-supercritical steam trap[J]. Journal of Sichuan University Engineering Science Edition, 2013, 45(5):145-153.(in Chinese)

【中文責編：晨 兮；英文責編：之 聿】

Simulation of the flow field characteristics of the boiler control valve pipe system of a power station

Liu Xiaoxu?, Liu Jia, He Qingzhong, and Guo Bin

Institute of Mechanical Engineering, Sichuan University of Science & Engineering, Zigong 643000, Sichuan Province, P.R.China

In order to analyze the problem that the piping end cover behind the regulator valve used in the water storage tank regulator of one thermal power unit is damaged during operating, the computational fluid dynamics (CFD) model for the valve and the piping is established, and the characteristic of two-phase flow field is studied by using CFD software CFX. The results show that cavitation is serious near the gap of the stem, which causes greater erosion force near the baseboard, and the computational results are consistent with that of the actual damage case. The piping structure of the regulator valve is improved according to the results, and numerical simulation of the flow field is done for the improved piping structure. The simulation results of the improved piping structure show that the flow speed is reduced by nearly 70%, and the bubble collapse is reduced by over 90%, the cavitation phenomenon is suppressed. The results can be used as the reference for the structural design and optimization of this kind of valves and piping.

computational fluid dynamics; regulating valve; two-phase flow; CFX software; cavitation erosion; structure optimization

Received：2016-03-09；Revised：2016-11-10；Accepted：2016-11-25

Foundation：National Science and Technology of National Energy Administration(2011HVCHT095)

? Corresponding author：Professor Liu Xiaoxu. E-mail: xiaoxu_liu@163.com

：Liu Xiaoxu, Liu Jia, He Qingzhong, et al. Simulation of the flow field characteristics of the boiler control valve pipe system of a power station[J]. Journal of Shenzhen University Science and Engineering, 2017, 34(1): 105-110.(in Chinese)

TH 138.52

A

10.3724/SP.J.1249.2017.01105

國家能源局國能科技基金資助項目(2011HVCHT095)；自貢市科技創(chuàng)新苗子工程資助項目(2015CXM03)

劉曉敘(1957—)，男，四川理工學院教授. 研究方向：機械設計制造及其自動化. E-mail:xiaoxu_liu@.163.com

引 文：劉曉敘，劉 佳，何慶中，等. 電站鍋爐調(diào)節(jié)閥管配系統(tǒng)流場特性的數(shù)值分析[J]. 深圳大學學報理工版，2017，34(1)：105-110.