基于CFD的離心式高效噴濺裝置水力特性

陳倩 鄭源 章國芳 栗文靈 王秋景

摘要：為提高離心式高效噴濺裝置的布水均勻性，采用歐拉（Eulerian）多相流模型和重整規劃群（RNG）k-ε湍流模型對離心式高效噴濺裝置的水力特性進行了三維數值模擬，將其模擬的濺灑情況與試驗情況進行了對比，驗證了模型的正確性。并用此方法研究了離心式噴濺裝置在3種水頭和4種轉速下的水力特性，結果表明其在較廣的水頭和轉速范圍下均能實現均勻布水，且一定范圍內提高水頭和降低轉速可以提高布水均勻性，提出離心式高效噴濺裝置的建議工作轉速為60～120 r/min。

關鍵詞：冷卻塔；噴濺裝置；數值模擬；試驗；水力性能；均勻布水

中圖分類號：TK050 文獻標志碼：A 文章編號：1672-1683（2016）04-0158-05

Abstract：In order to improve the uniformity of water distribution of efficient centrifugal spray nozzle，the Eulerian model，one of multiphase flow，and the RNG turbulence model were chosen in three-dimensional numerical simulation to investigate the hydraulic characteristics of efficient centrifugal spray nozzle.Comparison of numerical simulation results with experimental results verified the correctness of the model.The hydraulic characteristics of efficient centrifugal spray nozzle were comparatively analyzed under three head conditions and four speed conditions with this numerical simulation method.The results showed that it could achieve uniform water distribution in a wide head range and a wide speed range.The uniformity of water distribution of efficient centrifugal spray nozzle could be improved by raising the head and reducing the speed in a certain range.Accordingly，the recommended working speed range of the efficient centrifugal spray nozzle was 60～120 r/min.

Key words：cooling tower；spray nozzle；numerical simulation；experiment；hydraulic characteristics；uniform water distribution

自然通風逆流濕式冷卻塔廣泛用作火電、核電的冷端設備，其換熱效率影響著整個電站的效益[1]。50 MW的機組，出塔水溫升高1 ℃，煤耗率增加1.52 g/（kW·h）[2]。噴濺裝置是安裝在冷卻塔內，噴淋濺散水滴的，對其進行水力分析及結構優化可有效改善布水均勻性，提高冷卻塔換熱效率[3]。華能巢湖電廠1號機組冷卻塔進行噴濺裝置改造后，出塔水溫降低了1.51 ℃，每年減少標煤消耗

3 750 t以上[4]。針對使用較廣的旋轉型噴濺裝置[5-7]，韓永春[8]提出了多層旋轉噴濺裝置，劉劍卿[9]提出了離心式高效噴濺裝置，本質上都是為使布水更均勻而進行的結構上的改進。目前，研究噴濺裝置水力性能多采用實驗方法[10-13]，但該方法需要實體的噴濺裝置，周期長，經費高。近年來，隨著CFD技術的發展，許多水力機械問題都采用CFD進行研究，并得到了較好的效果[14-17]，亦有采用CFD研究冷卻塔內布水特性對整塔效率的影響[1，18]，但用其對冷卻塔內具體結構進行研究的較少。鑒此，本文采用CFD軟件Fluent對離心式高效噴濺裝置進行了數值模擬，分別對其在3種水頭，4種轉速下的水力特性進行了研究，分析其布水性能。

1 離心式高效噴濺裝置

離心式高效噴濺裝置見圖1，其具有連接管、支撐架和濺水碟等結構，連接管上設有套管，支撐架上連有轉輪，轉輪有40個葉片。離心式高效噴濺裝置配有一系列不同口徑的套管，以滿足不同工況下的流量要求。

2 計算模型

2.1 模型的建立與邊界條件

按1∶1的實際尺寸建立離心式高效噴濺裝置模型，為避免邊界條件對計算結果的影響，建立如圖2 所示圓柱形計算區域，圓柱底面半徑5 m，高度為3 m。

考慮速度進口收斂性好，套管進口處設為速度進口邊界條件，速度大小與水頭匹配；計算區域進、出口分別設為壓強進、出口邊界條件，取各處環境值；并通過設置體積率來界定各邊界面上所流過的流體種類；對于旋轉區域轉輪部分采用多重參考系（MRF）處理，其中轉速與水管壓力的對應關系由模型試驗測得；其余管壁、計算域邊界等采用無滑移固體墻壁邊界條件；速度分布采用標準壁面函數處理。

2.2 多相流模型和控制方程

為了跟蹤不同時刻水流的運動軌跡和分布情況，選取非穩態求解器進行求解，時間步長0.002 s。將空氣相設為主相，液態水相設為第二相。考慮到液態水相在噴濺裝置工作過程中破碎成液滴體積分數大于10%，采用歐拉（Eulerian）模型。該模型把氣體和水滴看成兩種流體，在空間的各節點上都有這兩種流體各自不同的速度和密度，這些流體存在于同一空間中，并且發生相互滲透，但各自有著不同的體積分數，兩相體積分數和為1。

考慮到噴濺裝置在開敞式的空間工作，不考慮空氣和水流的壓縮性；考慮到在冷卻塔中噴淋區的水流溫降只有整體溫降的4%左右[19]，而冷卻塔里的溫降95%～97%左右來自于蒸發[7]，可得循環水在噴淋區的蒸發量為整體蒸發量的4%左右，可認為蒸發對噴濺裝置的水力特性影響很小，本文主要研究噴濺裝置的水力特性，所以不考慮蒸發等能量影響；歐拉模型的連續性方程和動量方程[20]為

2.3 網格劃分與求解算法

考慮到離心式高效噴濺裝置結構復雜，尤其是轉輪葉片等結構細密，選用了適應性較強的四面體非結構化網格，并在轉輪葉片處進行了網格加密。選用了網格單元為283萬、348萬、428萬三種網格劃分方案分別在1.6 m水頭下模擬計算，噴濺半徑分別為2.6 m，2.8 m和2.85 m，結果表明網格達到348萬時，網格再增加對結果影響很小，考慮到網格越多，計算成本越大，因此采用348萬網格單元。

由于水流從套管中流出，再撞擊濺水碟，屬于射流撞擊，選取重整規劃群（RNG）[WTB1X]k[WTBX]-ε湍流模型，其在ε方程中加了一個條件，有效地提高了精度[21]。壓力與速度的耦合方式選擇相耦合的SIMPLE算法，并保持其亞松弛因子為默認值，體積率的離散方式設置為QUICK，其余各項皆為二階精度。計算時，以殘差收斂到1e-5數量級和滿足質量守恒定律作為判斷計算收斂的依據。

3 模擬計算結果

3.1 模型驗證

為驗證模擬結果的可靠性，現將離心式噴濺裝置配以口徑32 mm短管，在0.8水頭工況下的模擬計算結果與實際工作情況進行對比，見圖3。從圖中可看出，模擬結果中的流線圖呈傘狀，與實際工作時水流沖擊濺水碟呈傘狀噴灑吻合；模擬結果中豎直方向離濺水碟0.8 m的水平面的水量分布云圖為環狀，中央為無水區，與實際工作時地面上的水量情況吻合；說明CFD數值模擬計算方法具有較好的準確性，可以借此研究噴濺裝置的水力特性。

3.2 離心式高效噴濺裝置水力特性分析

圖4 為離心式高效噴濺裝置配以口徑32 mm的套管在水頭分別為0.8 m、1.2 m、和1.6 m下工作時距離濺水碟0.8 m水平面上水相分布云圖，從圖中可看出：（1）離心式高效噴濺裝置的有效噴濺范圍均為近似環形，受支撐架結構的影響，其對應的位置上有缺口，即支撐架結構對應位置噴濺半徑較小，在有效噴濺范圍內水相體積分數較均勻，說明離心式高效噴濺裝置在水頭為0.8～1.6 m的范圍內能有效實現均勻噴灑。（2）隨著水頭增大，噴濺范圍內中空區域增大，噴濺半徑增大，即噴濺的整體范圍增大，水相體積分數普遍減小，即水相體積分數的最大值和平均值均減小。

圖5 為離心式高效噴濺裝置配以口徑32 mm的套管在不同水頭下工作時距離濺水碟0.8 m水平面上水相沿徑向分布圖，由圖可看出：（1）隨著水頭的增大，離心式高效噴濺裝置的有效噴濺范圍沿著徑向整體外移，即噴濺半徑和中空半徑都增大，分別從1.8 m增大到2.8 m和從0.8 m增大到1.4 m，且噴濺半徑與中空半徑的差值在增大，從1 m增大到1.4 m，即有效噴濺范圍增大。（2）隨著水頭的增大，離心式高效噴濺裝置的有效噴濺范圍內的水相體積分數最大值在減小，從0.1252%減小到0.08686%，即噴灑的更均勻。

表1為離心式高效噴濺裝置配以口徑32 mm的套管在不同水頭H下工作時的水力特性參數，由表可得出隨著H增大，噴濺裝置的流量系數增加，流量Q增大，噴濺半徑R增大，且增加的速度減緩，有效噴濺面積S增大，Q/S減小。根據管嘴的流量系數計算公式可知，噴濺裝置的流量系數增加，說明裝置套管處的局部水頭損失減小，即在一定范圍內適當提高裝置工作的H可有效提高能源利用率；Q增大顯然是由于H增大引起的套管內的流速增大；R增大是由于管內的流速增大，水流撞擊濺水碟獲得的徑向速度增大，這也導致了噴濺面積和中空面積都增大，但噴濺面積增大幅度比中空面積大，導致S增大，根據R與水流濺灑厚道額徑向速度正相關，可以將套管中出來的水流分股，部分直接撞擊濺水碟，部分降低速度后撞擊濺水碟，以減小中空面積，這也是目前應用較多的多級噴濺裝置的理論原理；Q/S減小說明單位面積上濺灑的水量較低，易實現更均勻布水。

4 轉輪轉速影響

由于國內水質的日益惡化，噴濺裝置工作的環境也日益惡劣，水硬度增加常會導致在噴濺裝置的濺水碟和轉輪葉片，甚至轉輪與支撐架連接的軸孔內結垢，進而降低轉輪轉速，因此研究噴濺裝置在不同轉速下的水力特性很有意義。

圖6為離心式高效噴濺裝置配以口徑32 mm套管在水頭1.6 m下，分別以60、120、180和240 r/min轉速工作時距離濺水碟0.8 m水平面上水相沿徑向分布圖，由圖可看出：（1）在各轉速下，水相沿徑向分布曲線趨勢一致，在徑向距離1.9 m左右出現水相體及分數最大值，即離心式噴濺裝置在各轉速下都能實現均勻噴灑。（2）隨著轉速增大，噴濺半徑減小，但轉速增加到180 r/min時，轉速再增大，噴濺半徑變化不大。（3）隨著轉速增大，中空范圍變化不大，結合噴濺半徑考慮，則噴濺有效范圍減小。（4）隨著轉速的增大，水相體積分數的最大值相差不大，相差最大是轉速60 r/min和120 r/min，為0.0046%。

表2為離心式高效噴濺裝置配以口徑32 mm套管在水頭1.6 m下以不同轉速工作的水力特性參數，由表可得：（1）隨著轉速的提高，噴濺半經減小， 180 r/min工況下其比60 r/min工況下小0.3 m，這是由于離心式噴濺裝置是自流旋轉，即轉輪旋轉的能量來源于水流的動能，當轉速增大時，消耗水流的動能增大，水流自身的動能降低，速度降低，同時轉速增大，使得周向速度增大，徑向速度減小，所以噴濺半徑就降低了；到180 r/min 時轉速再增加，噴濺半徑幾乎不變，是因為240 r/min葉片的扭矩值為負值，即在實際的自流旋轉時轉速達不到240 r/min，轉輪的轉速最大極限值要略低于240 r/min。（2）隨著轉速增大，有效面積減小，180 r/min工況下其比60 r/min工況下減小明顯，可見轉速在60～120 r/min內是比較合理的，即在選擇材料和加工轉輪與支撐架的軸連接部件時，無需盡可能減小摩擦，只要保持摩擦在合適的范圍內即可。

綜上：不考慮轉速對水相濺灑后的水滴粒徑大小的影響，離心式噴濺裝置的轉輪轉速對其水力特性影響較小，因此離心式噴濺裝置在長期工作后不會因結垢問題引起的轉速下降導致嚴重不利影響，為獲得較優的噴灑效果，建議在60～120 r/min內工作。

5 結論

（1）離心式高效噴濺裝置可以在較廣的水頭下實現均勻布水，且一定范圍內，提高水頭可以提高單個噴濺裝置的布水均勻性。

（2）離心式高效噴濺裝置水力特性受轉速影響較小，其長時間工作后在結垢導致轉輪轉速下降的情況下仍能實現均勻布水，且轉速在60～120 r/min時，能獲得較優的水力特性，實現較均勻布水。

（3）由于本文研究中試驗測量的數據不夠，需進一步完善實驗臺，獲得更多數據。

參考文獻（References）：

[1] 金臺，張力，唐磊，等.自然通風濕式冷卻塔配水優化的三維數值研究[J].中國電機工程學報，2012，32（2）：9-15.（JIN Tai，ZHANG Li，TANG Lei，et al.Three-dimensional numerical study on water-distribution optimization in a natural draft wet cooling tower[J].Proceeding of the CSEE，2012，32（2）：9-15.（in Chinese））

[2] 李秀云，林萬超，嚴俊杰，等.冷卻塔的節能潛力分析[J].中國電力，1997，30（10）：34-36.（LI Xiu-yun，LIN Wan-chao，YAN Jun-jie，et al.To analyze energy-saving potential of cooling tower[J].Electric Power，1997，30（10）：34-36.（in Chinese））

[3] 馬文忠，王國春.改善水塔的布水均勻性提高水塔的換熱效率[A].2012（第二屆）全國發電企業設備優化檢修技術研討會論文集[C].北京，2012.（MA Wen-zhong，WANG Guo-chun.Improvement of the heat transfer efficiency and water distribution uniformity in water tower[A].The 2th national Symposium on equipment optimization and maintenance technology of power generation enterprise [C].Beijing，2012.（in Chinese））

[4] 張世宏，孫立春，樊巖.華能巢湖電廠1#機組冷卻塔噴濺裝置節能改造[J].節能，2013（10）：74-76.（ZHANG Shi-hong，SUN Li-chun，FAN Yan.Huaneng Chaohu power plant 1# unit energy-saving cooling tower spray device[J].Energy Conservation，2013（10）：74-76.（in Chinese））

[5] 袁小偉，劉曉朵.冷卻塔噴嘴性能測試[J].山東工業技術，2015，（06）：54-55.（YUAN Xiao-wei，LIU Xiao-duo.Cooling tower nozzle performance test[J].Shandong Industrial Technology，2015，（06）：54-55.（in Chinese））

[6] 趙振國.冷卻塔[M].北京：中國水利水電出版社，1996.（ZHAO Zhen-guo.Cooling Tower[M].Beijing：China Water & Power Press，1996.（in Chinese））

[7] 史佑吉.冷卻塔運行與試驗[M].北京：水利電力出版社，1990.（SHI You-ji.Operation and test of cooling tower[M].Beijing：Water Conservancy and Electric Power Press，1990.（in Chinese））

[8] 韓永春.多層旋轉噴濺裝置[P].中國：200920010212.2，2009-10-28.（HAN Yong-chun.The multi-layer rotary spray nozzle[P].China：200920010212.2，2009-10-28.）

[9] 劉劍卿，于秋江.離心式高效噴濺裝置[P].中國：200920302958.0，2010-05-19.（LIU Jian-qin，YU Qiu-jiang.efficient centrifugal spray nozzle[P].China：200920302958.0，2010-05-19.）

[10] Thirapong Muangnoi，Wanchai Asvapoositkul，Pibool Hungspreugs.Performance characteristics of a downward spray water-jet cooling tower[J].Applied Thermal Engineering，2014，69：165-176.

[11] P Panwaree，W Asvapoositkul.Droplets Analysis Of Cooling Tower Spray Nozzle by Photography[A].The 11th Asian Symposium on Visualization[C].Japan，2011.

[12] Gao M，Sun F-z，Turan A.Experimental study regarding the evolution of temperature profiles inside wet cooling tower under crosswind conditions[J].International Journal of Thermal Sciences.2014，86：284-91.

[13] Chen Y，Sun F，Wang H，et al.Experimental research of the cross walls effect on the thermal performance of wet cooling towers under crosswind conditions[J].Applied Thermal Engineering.2011，31（17-18）：4007-13.

[14] 熊妍，屈波，霍志紅，等.冷卻塔專用超低比轉速水輪機的設計與數值模擬[J].南水北調與水利科技，2014，12（3）：112-115.（XIONG Yan，QU Bo，HUO Zhi-hong，et al.Design and numerical simulation of turbine with super low specific speed used in cooling tower[J].South-to-North Water Transfers and Water Science &Technology，2014，12（3）：112-115.（in Chinese））

[15] 羅欣，鄭源，馮俊.軸流泵內部流場的數值模擬[J].水電能源科學，2012，30（12）：120-122，145.（LUO Xin，ZHENG Yuan，FENG Jun.Numerical simulation of interior flow filed of axial flow pump[J].Water Resources and Power，2012，30（12）：120-122，145.（in Chinese））

[16] 牟童，鄭源，鄭建坤，等.T型攪拌器數值模擬研究[J].南水北調與水利科技，2014，12（3）：89-94.（MOU Tong，ZHENG Yuan，ZHENG Jian-kun，et al Numerical simulation of T-shaped stirrer[J].Sou-to-North Water Transfers and Water Science & Technology，2014，12（3）：89-94.（in Chinese））

[17]張強，鄭源，陳會向.基于CFD 的潮流能水輪機數值模擬研究[J].南水北調與水利科技，2015，13（3）：431-434.（ZHANG Qiang，ZHEN Yuan，CHEN Hui-xiang.Numerical simulation of tidal current energy turbine based on CFD[J].Sou-to-North Water Transfers and Water Science&Technology，2015，13（3）：431-434.（in Chinese））

[18] 鄭水華.超大型冷卻塔內氣液兩相流動和傳熱傳質過程的數值模擬研究[D].杭州：浙江大學，2012.（ZHENG Shui-hua.Numerical simulation of gas-liquid flow，heat and mass transfer process in super-large cooling tower[D].Hangzhou：Zhejiang University，2012.（in Chinese））

[19] Rafat Al-Waked，Masud Behnia.CFD simulation of wet cooling towers[J].Applied Thermal Engineer，2006，26：382-395.

[20] FLUENT.User's guide[M].Lebanon：Fluent Inc.，2003.

[21] 李鵬飛，徐敏義，王飛飛.精通CFD工程案例與案例實戰[M].北京：人民郵電出版社，2011（10）：122-131.（LI Peng-fei，XU Min-yi，WANG Fei-fei.Master CFD engineering case and case study[M].Beijing：People's Posts and Telecommunications Press，2011（10）：122-131.（in Chinese））