雙流程電站凝汽器水室的流動分析與改造

邢樂強，張莉，白路平，甘志聯

（上海電力學院，能源與機械工程學院，上海，200090）

雙流程電站凝汽器水室的流動分析與改造

邢樂強，張莉，白路平，甘志聯

（上海電力學院，能源與機械工程學院，上海，200090）

為能夠準確預測凝汽器水側流動特性，應用計算流體力學商用軟件，管束區域采用多孔介質模型，對某電廠雙流程凝汽器水室流場進行了數值分析。結果表明：該電站雙流程凝汽器入口水室結構存在一定的缺陷，導致冷卻水在局部產生了較大的阻力，形成了明顯漩渦。針對凝汽器水室中漩渦產生的原因，對凝汽器水室結構進行了優化改造，并且在此基礎上提出了更多合理的建議。

凝汽器，水室，數值模擬，分析，改造

1　前言

汽輪機系統作為火力發電機組中的重要組成部分，一直都是提高機組效能的主要研究對象。我國火電機組中廣泛采用凝汽式汽輪機，凝汽器工作性能的好壞直接影響到汽輪機系統的性能，而且凝汽器有著較大的節能潛力［1］，所以，長期以來國內外研究者們對凝汽器也傾注著一定的研究精力。但是，多數文獻的研究表明，有關凝汽器的研究主要集中在汽側［2］，即主要研究汽側內部汽氣混合物的流動和換熱特性，而凝汽器水側方面的研究還相對較少［3］，但是凝汽器水側方面的研究與汽側研究具有同樣重要的意義。

凝汽器冷卻水室內的流動是凝汽器水側研究的重點，合理的水室結構能夠消除水室內較大的渦流區域，減小冷卻水阻，可以使冷卻水均勻流入冷卻管內，提高凝汽器冷卻效果［4］，同時，還可以帶動膠球散布進入所有冷卻管，有助于冷卻管的清洗。雖然之前的凝汽器水室已經隨著加工工藝的發展由長方體結構優化為弧形結構，但是這樣的水室流動特性如何？是否還可以在目前制造工藝水平不斷提高的情況下進一步改進？這些都還有待于開展進一步的研究工作。

本文采用數值模擬的研究方法，研究了某電廠雙流程凝汽器水側的流場特性，在得到凝汽器水側流場速度、壓力等參量分布的基礎上，對該凝汽器的水室進行了評價分析，進而對該水室進行了改造嘗試。

2　數值計算方法

2.1物理模型

凝汽器水側的流動包括冷卻水在進水管、水室、冷卻水管、出水管等幾部分區域內的流動。對于雙流程凝汽器，水室又包括進水室、折回水室和出水室。在穩定工況下，凝汽器水側的流動可以看作是三維、定常的湍流流動?？紤]到冷卻水在凝汽器水側的溫升僅有10℃左右，冷卻水的物性變化不大，流動過程中可以不考慮冷卻水能量的變化。

此外，由于凝汽器管束區域的冷卻水管數目較多，而且每根冷卻管的管徑尺寸和凝汽器本體的尺寸相比差異較大，冷卻水在該區域流進冷卻水管，就好像流進了多孔區域。因此，冷卻水流過凝汽器管束區域的流動可以用多孔介質模型表達［5］。

2.2控制方程

2.2.1基本控制方程

根據上面對冷卻水側物理現象的簡化，描述凝汽器水側流動的控制方程包括連續性方程、動量方程和湍流模型方程，該流動過程可統一描述如下［6］

式中，β為蒸汽流經區域的多孔率；變量Ф分別為1，u→，k和ε時，方程分別表示混合物的連續性方程、動量方程、湍流模型。

2.2.2多孔介質模型

多孔介質模型在本文計算中應用于冷卻水在凝汽器管束區冷卻水管內的流動。由于冷卻水一旦進入冷卻水管，就會被限制在管內流動，不會出現在通流截面的二次流動，因此，多孔介質在該區域只需要考慮沿冷卻水管長度方向上阻力的模擬。

根據HEI推薦的凝汽器冷卻管內摩擦損失的計算方法，冷卻水在冷卻管內的摩擦損失Ff可用式 （2）、式 （3）進行計算。

式中，F0為 25.5℃冷卻水在管壁厚為12.4 mm的冷卻水管內流動時的摩擦損失；Vw為凝汽器中冷卻水管內的流速；d2為冷卻水管內徑；Ft和Fw為當冷卻水平均溫度和冷卻水管壁厚度不符合上述要求時所加的修正系數［7］。

此外，冷卻水流進、流出換熱管時會產生局部損失，在數值計算中這部分不能忽略，局部損失采用式 （4）進行計算［8］：

其中，Vw表示冷卻水流速；g為重力加速度；ξ為局部損失系數。冷卻水從水室流入換熱管以及從換熱管流入水室時的局部損失系數近似分別取為0.5和1。

3　研究對象簡介及計算條件設置

3.1研究對象簡介

本文以某單殼體、雙通道、雙流程 （對分式水室，即冷卻水雙進雙出）、表面式凝汽器為研究對象。該機組凝汽器的2個流通通道具有相同的幾何結構，故本文只選取了1個通道作為模型進行計算。圖1給出了凝汽器某一通道流動區域的幾何示意圖。

圖1　凝汽器幾何結構圖

3.2求解設置

采用Gambit軟件對所研究流體區域進行計算網格劃分，計算網格數約為200萬個。

應用FLUENT數值計算軟件進行計算，計算中，對于多孔介質區域阻力的計算啟用了FLUENT軟件的用戶自定義函數 （UDF）功能，編寫、編譯了DEFINE_SOURCE函數用于計算多孔區域動量方程的阻力源項。

基于凝汽器的設計參數，對計算的邊界條件進行了設置：

（1）入口邊界：進水管入口設定為速度入口邊界條件。入口流速的大小按設計工況流量取平均值，方向垂直于進口界面；

（2）出口邊界：出水管出口設定為壓力出口邊界條件；

（3）固體壁面：設定為無滑移的邊界條件，采用標準壁面函數法；

（4）多孔區域：冷卻水在冷卻管內的流動區域，多孔率為0.493。分布阻力的源項由UDF加載到多孔介質模型中。

4　計算結果及分析

采用前面所述的方法，計算得到了凝汽器水側的流動參數分布。由于缺少凝汽器水側詳細的實驗測量結果，本文通過水側壓降的對比進行了數值計算的驗證。通過FLUENT計算得到的冷卻水流程總的壓降大約為64.556 kPa，和設計水阻69.79 Pa相比較，誤差大約為7.5%，符合工程設計的誤差要求，在一定程度上表明了計算結果的合理性。

4.1冷卻水流線分布

圖2給出了該雙流程凝汽器水側冷卻水的流線分布。從圖2可以看出，在進水室下部，當冷卻水流入水室時，隨著流動空間的增大，冷卻水向兩側流動，流線會顯得很雜亂，而在進水室的上部，由于冷卻水在進水室流向的改變，冷卻水不易流到那里，也造成了進水室上部流線一定程度的雜亂。此外，在靠近冷卻水管板處，冷卻水流進冷卻水管的局部阻力也造成了該處的流動雜亂。

圖2　凝汽器水側冷卻水流線分布圖

冷卻水管的管束區內，冷卻水表現出沿軸向的流動平順性，這與冷卻水管將冷卻水限制在管內流動的物理現象相符合，表明了借助多孔介質模型模擬管束區域內流體流動的合理性。在折回水室內，冷卻水的流線相較于進水室平順，出水室中在冷卻水流入出水管前，在出水管的對面出現了一個較小的漩渦，這是由于冷卻水向著出水室出口流動，在另一側發生脫流作用產生的。

從冷卻水流線分布可以預測冷卻水流動整體的平順性，但是只有對冷卻水在各個水室內的速度和壓力分布進行詳細的了解之后才能分析出水室結構的不足之處，從而提出相應的改造方案。

4.2進水室流場分布

圖1中在進水室內沿進水管軸向由下往上依次定義了3個截面來觀察進水室的流場分布。圖3和圖4分別給出這3個截面上的速度和壓力分布。從圖3、圖4可以看出，在進水室的上側壓力較大，冷卻水流速較低，而且距離管板越遠的地方流速越低，這里漩渦的產生是由于冷卻水由y軸的正方向轉向z軸負方向流動所導致的。還可以看到，在進水室下側存在2個明顯的漩渦，原因在于冷卻水流入進水室時，隨著空間的擴大，冷卻水向水室兩側流動，從而形成渦流。在進水室設置有冷卻水管，也就是說這一區域里冷卻水流動的雜亂是由管束布置所決定的，無法通過凝汽器進水室結構的改造來改善這部分的流動狀況。

圖3　進水室內的速度分布

圖4　進口水室內壓力分布

圖5　進口管板截面上的流場分布

圖5所示為凝汽器第一流程管束近管板處流場分布，由于冷卻水的流動方向與幾何建模時z軸的正方向相反，因此，在與進水室相銜接的管板上的冷卻水流速為負。從圖5（a）可以看出，在冷卻水管比較密的區域，流速在z方向的分量大部分在1～2.6 m/s的范圍內，其中，速度范圍在2.4～2.6 m/s的區域主要集中在截面上部。這是因為冷卻水由進水管軸向進入，由于流動的慣性，冷卻水沖到進水室上方的趨勢，因此，位于上側的換熱管內的流量相對于下側的流量要大。此外，在冷卻水管布置比較稀疏的外圍區域，冷卻水流速比較低，大概在1～2 m/s的范圍內，原因是冷卻水流入換熱管時受到周圍管板的阻礙大，降低了冷卻水的流入速度。

管束內冷卻水流場分布的不均勻性，除了有來自進水室幾何結構的影響外，還有換熱管管束布置的影響。所以，凝汽器進水室合理的幾何結構改造可以帶來管束區冷卻水流場分布的改善，但是冷卻水管束的布置方式所造成的影響則無法通過進水室的改進來消除。

4.3折回水室的流場分布

整體上來說，冷卻水在折回水室產生的壓力損失并不多，通過計算得出這部分的壓力損失是964 Pa。從圖6可以看出，在折回水室的中間且遠離管板的區域內水的流速較大，壓力較低，而在上下兩側水的流速較低，壓力較大，這是由于流體在轉彎時與壁面發生了碰撞所造成的，這部分的動能損失是無法避免的。同樣，由于管板中間區域和上、下2個管束區的中央位置沒有布置換熱管，在這些位置也存在低速區。由此折回水室內流場的分布特點是由凝汽器水室的整體布置和管束布置類型所決定的。

圖6　折回水室的流場分布

4.4出水室及出水管流場分布

圖7是出水室和出水管流場分布圖。從圖7可以看出，出水室壓力從管板向出水管側逐漸降低，在遠離管板側出現了高壓低速區域。這是因為冷卻水由沿著z軸正方向流動轉向沿著x軸正方向流動造成的，冷卻水在轉向過程中與壁面發生碰撞會產生動能的損失。冷卻水從空間較大的出水室流入空間較小的出水管，冷卻水流速突然增大，最高可達3.34 m/s，而且剛進入出水管的時候冷卻水主要集中在一側，使得這一側水室壁面和出水管壁會受到很大的沖擊和摩擦。

圖7　出水室及出水管的流場分布

5　水室結構改造方案及效果預測

圖8　水室改造后的計算區域

5.2改造效果預測

本文對改進后的雙流程凝汽器水側流動又進行了數值模擬從而來預測改造方案的效果。從冷卻水流線分布圖 （見圖9）可以看出，凝汽器進水室采用了新的結構后，水室上部的漩渦消失，速度流線更加平順有序，而且在這一區域的冷卻水流速也沒有原來那么低。另外從冷卻水流經各個階段時產生的壓力損失也可以看到改造的效果。

圖9　改造后的冷卻水流線分布

5.1水室結構改造方案

整體上來說，冷卻水在整個流動過程中，除了2個流程的管束區域以外，進水室、出水室這兩部分產生的壓力損失最多，而且從各個水室內的流場分布也可以看出，冷卻水在流動過程中存在一些流動不合理的地方。經過分析，在雙流程凝汽器內冷卻水的流動特性是由多方面因素決定的，包括凝汽器水室的幾何結構以及其他原因。針對水室結構對凝汽器水側流場分布的影響，可以對水室采取相應的改造措施。為了消除進水室上部的渦流，需要對凝汽器進水室結構進行改進。水室改造后的計算區域如圖8所示。改造前在進水室內的壓力損失是5 465.203 Pa，改造之后降為5 261.358 Pa，這表明凝汽器進水室結構改造還是取得了一定的效果。

圖10　改造后進水室流線分布

從圖10還可以看出，在進水室下部依然存在渦流，主要還是因為冷卻水從空間較小的進水管流入空間較大的進水室的原因，如果要消除進水室下部的渦流，就需要將進水室設計為漸擴型的幾何結構，然后同樣對改造后的凝汽器水側進行數值模擬，這也正是凝汽器水室結構優化研究接下來所要做的工作。

6　結論

（1）該電站凝汽器進口水室結構中存在明顯不合理的地方，冷卻水在水室上部和近管板處產生了較大的漩渦，導致冷卻水對管板產生較大的沖擊和摩擦；管板上冷卻水速度分布并不均勻，主要表現為上側有一部分區域的流量比較集中。在后水室內由于冷卻水流動方向發生變化，產生漩渦，但是范圍不大，并不會對循環冷卻水產生太大的阻力。而在出口水室壓力從管板向出水管側逐漸降低，在遠離管板側出現了高壓低速區域；在剛進入出水管時候冷卻水主要集中在一側，使得這一側水室壁面和出水管壁會受到很大的沖擊和摩擦。

（2）根據對凝汽器進水室中漩渦存在原因的分析，為了消除進水室內的漩渦，提出相應的改造方案，然后對改造后的模型重新進行數值模擬，通過與改造之前的計算結果進行對比分析后可以發現，改造后進水室上部的漩渦范圍和程度都有所減小。

（3）雖然在凝汽器進水室的管板附近的漩渦不再那么明顯，但是下部的漩渦依然存在，而且流過管板截面的冷卻水流速分布也依然存在不均勻性。因此，建議以后的改造研究過程中可以考慮將進水室設計為漸擴型的幾何結構，或者在進水室下部靠近進水管的位置處加裝合理的導流裝置，以減小渦流，提高管板截面處冷卻水流速的均勻性。

［1］尤清華．大型燃煤發電機組冷端系統優化設計研究［D］．北京：華北電力大學，2012

［2］曲建麗．電站凝汽器的數值模擬研究與應用［D］．濟南：山東大學，2007

［3］李永華，侯淑娟，孫剛．300 MW機組電站凝汽器水側數值模擬［J］．電站輔機，2009，30（2）：12-15

［4］蔣建飛，黃樹紅，王坤．凝汽器水側流動的三維數值模擬［J］．動力工程，2006，26（2）：249-252

［5］李勇，李立言，曹麗華，等．山谷型雙流程凝汽器管側流場的數值分析［J］．化工機械，2012，39（3）：351-355

［6］汪國山．電站凝汽器熱力性能數值仿真及其應用［M］．北京：中國電力出版社，2010：18-21

［7］張卓澄．大型電站凝汽器［M］．北京：機械工業出版社，1993：54-57

［8］柯葵，朱立明．流體力學與流體機械［M］．上海：同濟大學出版社，2009：85-88

Flow Analysis and Alteration of Water Chamber in a Double-flow Condenser of Power Plant

Xing Leqiang，Zhang Li，Bai Luping，Gan Zhilian

（College of Energy and Mechanical Engineering，Shanghai University of Electric Power，Shanghai，200090）

In order to predict the flow characteristic of condenser water side，the flows field of a double-flow condenser whose pipe bundle has porous medium model is simulated through CFD software．The result shows that the structure of condenser's inlet water chamber has defect which can result obvious resistance in some local parts where large vortex appears．Based on the cause of vortex，some measures have been carried out and some suggestions have been put forward to improve condenser's water chamber.

condenser，water chamber，numerical simulation，analysis，alteration

TK264

A

1674-9987（2015）02-0019-06

10.13808/j.cnki.issn1674-9987.2015.02.004

上海市科委專項科研項目。

邢樂強 （1989-），男，碩士研究生在讀，就讀于上海電力學院能源與機械工程學院。