賀曉怡，陰繼翔，王 濤，馬素霞

?

擾流孔及翅片間距對直接空冷蛇形翅片管流動與換熱的影響

賀曉怡，陰繼翔，王 濤，馬素霞

（太原理工大學(xué)電氣與動力工程學(xué)院，山西 太原 030024）

在傳統(tǒng)直接空冷單排蛇形翅片扁平管上增設(shè)擾流孔以強化換熱。本文應(yīng)用Fluent軟件對 有、無擾流孔的2種蛇形翅片扁平管（蛇形翅片管）模型進行三維數(shù)值模擬，研究了1.4、2.2、3.0、3.8 m/s 4種迎面風(fēng)速下不同翅片間距對有、無擾流孔的2種型式翅片管換熱、阻力以及綜合性能的影響。結(jié)果表明：與未開孔的蛇形翅片管相比，增設(shè)擾流孔的蛇形翅片管在各種翅片間距、迎面風(fēng)速下使空氣側(cè)的換熱系數(shù)提高1%～18%，且流動阻力略微減小；在較小的迎面風(fēng)速范圍內(nèi)，增設(shè)擾流孔可以使相同工況下的最佳翅片間距減小。

直接空冷；蛇形翅片管；擾流孔；翅片間距；換熱系數(shù)；流動阻力；數(shù)值模擬

直接空冷系統(tǒng)廣泛應(yīng)用于我國富煤缺水地區(qū)的火力發(fā)電機組中[1-3]，其核心部件是翅片管束結(jié)構(gòu)的空冷凝汽器。20世紀70年代至今，直接空冷凝汽器所采用的翅片管經(jīng)歷了4排等管徑連續(xù)翅片結(jié)構(gòu)、雙排橢圓管矩形翅片結(jié)構(gòu)、3排橢圓管橢圓形翅片結(jié)構(gòu)及單排扁管蛇形翅片結(jié)構(gòu)4個階段[4]。其中，單排蛇形翅片扁平管因其具有換熱面積大、流動阻力小、防凍性能好、傳熱系數(shù)高等諸多優(yōu)點而成為直接空冷換熱元件的主要結(jié)構(gòu)形式[5-6]。

針對直接空冷單排蛇形翅片扁平管空氣側(cè)的換熱與流動性能，已有很多學(xué)者進行了研究。文獻[7-8]分別通過數(shù)值模擬和風(fēng)洞實驗的方法研究了迎面風(fēng)速對蛇形翅片扁平管（蛇形翅片管）換熱系數(shù)的影響，結(jié)果表明氣側(cè)換熱系數(shù)隨迎面風(fēng)速的增大而逐漸增大。申微波[9]研究了蛇形翅片的厚度對翅片管換熱系數(shù)與流動阻力的影響，發(fā)現(xiàn)不同迎面風(fēng)速下翅片管的對流換熱系數(shù)與流動阻力均隨翅片厚度的增加而增大。馬榮榮[10]研究發(fā)現(xiàn)，蛇形翅片管的換熱系數(shù)與流動阻力隨著翅片間距的的增大而減小。侯雪梅[11]模擬研究了蛇形翅片管擾流孔距對換熱系數(shù)的影響，發(fā)現(xiàn)在一定范圍內(nèi)，換熱系數(shù)隨擾流孔間距的減小而增大，但在超過一定界限后，換熱系數(shù)無明顯提高。文獻[12-13]研究發(fā)現(xiàn)，在矩形翅片上加設(shè)擾流孔后翅片管的換熱性能得到提高。

綜上，開孔翅片可以提高翅片管的換熱性能，但針對直接空冷蛇形翅片管的研究主要集中在迎面風(fēng)速與蛇形翅片尺寸上，對于蛇形翅片上增設(shè)擾流孔強化換熱的文獻較少。本文建立了相同結(jié)構(gòu)尺寸下有、無擾流孔的2種翅片管數(shù)值模型，通過應(yīng)用Fluent軟件研究了1.55、1.85、2.15、2.45、2.75、3.05、3.35 mm 7個翅片間距的蛇形翅片管在4種迎面風(fēng)速下的換熱、流動阻力以及綜合性能，相關(guān)結(jié)論可為蛇形翅片管的研究及選擇提供參考依據(jù)。

1 幾何模型及區(qū)域網(wǎng)格化

圖1為實際工程中采用的蛇形翅片管。為便于建模，將蛇形翅片簡化為矩形翅片。由于蛇形翅片的幾何形狀與流動換熱皆具有對稱性，選擇翅片單元（一個“Z”形的翅片與翅片兩側(cè)的各一半流動空間）為研究對象，計算區(qū)域如圖2所示。

為防止入口效應(yīng)和出口回流的影響，分別在進、出口區(qū)域設(shè)置長度為50 mm和400 mm的延長區(qū)。扁平管采用鋼管，管外釬焊鋁合金蛇形翅片。扁平管長軸為=219 mm，短軸為=19 mm，壁厚為=1.5 mm；翅片的長度為=200 mm，高度為=19 mm，厚度為=0.3 mm，間距為；位于翅高中心處的擾流孔個數(shù)為=15，孔徑為=3 mm。

使用Gambit軟件生成計算域網(wǎng)格。為了便于劃分網(wǎng)格并獲得高質(zhì)量的網(wǎng)格，對計算域進行了分塊處理。進口區(qū)域與出口區(qū)域結(jié)構(gòu)簡單，使用結(jié)構(gòu)化網(wǎng)格。翅片換熱區(qū)域速度變化劇烈，采用給壁面附近和開孔處及其周圍加密的方法。開孔處使用非結(jié)構(gòu)化網(wǎng)格，其余部分采用結(jié)構(gòu)化、非均勻的六面體網(wǎng)格。翅片的局部區(qū)域網(wǎng)格結(jié)構(gòu)如圖3所示。

2 控制方程與邊界條件

本文的模擬計算基于以下假設(shè)：1）流體為不可壓縮的常物性空氣；2）不考慮翅片與基管間的接觸熱阻；3）不考慮翅片管的自然對流與輻射換熱；4）扁平管內(nèi)壁面溫度恒定；5）流體在管外翅片間的流動為穩(wěn)態(tài)。

控制方程[14]如下：

動量方程：

式中：u為流體在X方向上的速度；、=1,2,3，分別代表、、方向；為流體密度；為流體微元上的靜壓力。

流體進口面為速度進口，界面上速度及溫度均勻。出口為壓力邊界條件，無回流發(fā)生。扁平管內(nèi)壁面設(shè)為等溫壁面，外壁面與翅片表面均為耦合面，其他表面設(shè)置為對稱邊界條件。

選擇適用于求解不可壓縮流體流動計算的壓力求解器。其中，速度與壓力耦合采用Simple算法，對流項的離散采用Power Law格式。

3 模型驗證

3.1 網(wǎng)格無關(guān)性驗證

數(shù)值模擬作為一種近似求解方法，其網(wǎng)格數(shù)目直接影響計算結(jié)果精度。為了選取合適的網(wǎng)格 數(shù)目，須先進行網(wǎng)格無關(guān)性檢驗。以翅片間距 2.75 mm、迎面風(fēng)速2.2 m/s工況下的開孔蛇形翅片管為例，采用6種網(wǎng)格數(shù)目對開孔蛇形翅片管的數(shù)值模型進行計算，結(jié)果如圖4所示。由圖4可知，網(wǎng)格數(shù)為1 186 752與1 325 454時換熱系數(shù)的相對誤差為0.1%，這表明選取網(wǎng)格數(shù)約為118萬的模型進行后續(xù)模擬計算較為合理。

3.2 算法驗證

為了驗證本文算法及模型的合理性與有效性，使用文獻[14]中扁平管蛇形翅片模型進行模擬計算，得到傳熱系數(shù)與流動阻力的變化關(guān)系。圖5為本文模擬結(jié)果與文獻[14]的比較。由圖5可知，模擬值與文獻[14]數(shù)據(jù)相比誤差較小，傳熱系數(shù)與流動阻力的最大誤差為2.94%和7.80%。該誤差在可接受范圍內(nèi)，驗證了本文算法及模型的正確性。

4 計算結(jié)果分析

本文數(shù)據(jù)處理的相關(guān)參數(shù)定義[15]如下：

式中，m為空氣定性溫度，in、out分別為計算域進、出口空氣溫度，Δ、e、、m、分別為進出口壓差、翅片管當(dāng)量直徑、空氣密度、翅片管最小截面處流速、翅片沿空氣流動方向的長度（取0.2 m），、、分別為努塞爾數(shù)、雷諾數(shù)和普朗特數(shù)。

4.1 速度場

圖6為迎面風(fēng)速=1.4 m/s時有、無擾流孔的翅片管翅片間距中心截面的速度分布?？諝庥勺蠖诉M入，流動至扁管時流速逐漸增大。當(dāng)空氣流經(jīng)翅片通道時其流速達到最大值。翅片間空氣的平均流速為進口迎面風(fēng)速的1.6~1.8倍。隨著翅片間距的增加，流速提高的倍數(shù)越小。對比圖6b)、d)可見，在翅片間距較大時翅片間空氣流速更小。產(chǎn)生此現(xiàn)象的原因在于，與進口處相比，翅片區(qū)域的流動截面積減小。隨著翅片間距的增大，翅片厚度不變，翅片厚度截面在整個計算域截面上所占比例逐漸變小，而流通截面面積所占比例越來越大，速度提高的比例也因此越來越小。

從圖6中還可以看出：空氣在進入翅片通道前，有、無擾流孔2種翅片通道內(nèi)的流動狀況基本相同；當(dāng)空氣流經(jīng)翅片開孔處，流動空間瞬間增大，流速明顯減小，空氣在流動過程中受到擾動，且阻礙了邊界層的持續(xù)發(fā)展，不過擾流孔對邊界層的影響僅局限于翅片高度方向上的中心部分；在擾流孔翅片間的流動中，翅片根部與頂部的空氣流速未受到影響；在同一迎面風(fēng)速下，翅片間距增大，擾流孔對翅片間空氣流動的影響減弱；在翅片間距為3.35 mm的流道中剖面的速度分布基本沒有受到孔的擾動，因此，隨著翅片間距的增加，擾動減小，換熱強化效果減弱。

圖7為翅片間距=1.55、2.45 mm時有、無擾流孔翅片管在翅片高度9.5 mm處橫截面的速度場。由圖7可見：邊界層變化明顯，平直蛇形翅片管的流道內(nèi)速度分布相同、流動均勻；由于擾流孔的存在，翅片截面被間斷，在間斷處，空氣流動明顯受到擾流孔的影響，邊界層經(jīng)歷了不斷地斷裂和形成，從而增強了換熱。

圖7 不同翅片間距時翅片高度9.5 mm處橫截面速度分布

Fig.7 The velocity distribution in cross section with fin height of 9.5 mm under conditions with different fin pitches

4.2 換熱性能

不同迎面風(fēng)速下?lián)Q熱系數(shù)隨翅片間距的變化關(guān)系如圖8所示。

Fig.8 Changes of the heat transfer coefficient with fin pitch

由圖8可見：在不同迎面風(fēng)速下，翅片管換熱系數(shù)的變化趨勢一致；翅片管的換熱系數(shù)隨著翅片間距的增加、迎面風(fēng)速的減小而減??；在所計算的迎面風(fēng)速與翅片間距下，帶擾流孔的翅片管的換熱系數(shù)均大于無擾流孔的翅片管，換熱系數(shù)提升1%~18%；擾流孔對換熱系數(shù)的提高程度與翅片間距、迎面風(fēng)速有關(guān)，在較小的迎面風(fēng)速與翅片間距范圍內(nèi)，增設(shè)擾流孔對換熱系數(shù)的提高更為顯著。

4.3 流動阻力

4種迎面風(fēng)速下有、無擾流孔2種類型翅片管的阻力系數(shù)隨翅片間距的變化關(guān)系如圖9所示。由圖9可見，有、無擾流孔2種翅片管的阻力系數(shù)隨著迎面風(fēng)速與翅片間距的增大而減小，減小的趨勢逐漸變緩。對比各工況下2種翅片管的流動性能發(fā)現(xiàn)，擾流孔對阻力系數(shù)的影響隨著翅片間距的增大逐漸減小，直至影響基本消失。

4.4 綜合性能

換熱能力的提高通常是以阻力損失的增大為代價，為此需要對換熱技術(shù)進行綜合性能評估。本文以/[9]作為綜合性能的評價指標，考察翅片管的綜合性能。該值越大，在相同泵功下完成同樣的換熱量所需的迎風(fēng)面積越小。圖10為有、無擾流孔翅片管在不同迎面風(fēng)速下/隨翅片間距的變化。由圖10可見：除小迎面風(fēng)速（=1.4 m/s）外，其他迎面風(fēng)速下隨著翅片間距的增大，/呈先增大后減小趨勢，且翅片間距較小時更為明顯；迎面風(fēng)速越大，/變化趨勢變緩，最佳翅片間距隨著迎面風(fēng)速的增加而減小。對于2種蛇形翅片管，不同迎面風(fēng)速下的最佳綜合換熱結(jié)構(gòu)也不盡相同。按迎面風(fēng)速由小至大排序，無孔翅片管的最佳翅片間距依次為3.35、3.05、2.75、2.45 mm，開孔翅片管的最佳翅片間距依次為3.05、2.75、2.75、2.45 mm。這表明，在迎面風(fēng)速較小時增設(shè)擾流孔可以減小翅片間距，從而使得翅片換熱結(jié)構(gòu)更為緊湊。

5 結(jié) 論

1）翅片間空氣流速為進口迎面風(fēng)速的1.6~ 1.8倍。翅片間距越小，翅片區(qū)域的風(fēng)速提高的倍數(shù)越大，這是導(dǎo)致翅片間距較小時取得較高換熱系數(shù)的原因之一。

2）增設(shè)擾流孔后翅片管的換熱系數(shù)提高了1%~18%，且提高的程度隨著翅片間距與迎面風(fēng)速的增大而逐漸減弱；擾流孔對翅片管阻力性能的影響甚微?？梢姡瑢τ谳^小的迎面風(fēng)速和翅片間距，可以在不增加能耗的條件下改善翅片管換熱性能。

3）隨著翅片間距的增加，4種翅片管的綜合性能系數(shù)先增大后減?。?.4 m/s除外）。當(dāng)迎面風(fēng)速較小時，增設(shè)擾流孔使得最佳翅片間距減小，單位容積內(nèi)可以布置較多翅片使得換熱器結(jié)構(gòu)更為緊湊。

［1］ 溫高. 發(fā)電廠空冷技術(shù)[M]. 北京: 中國電力出版社, 2008: 24-26.

WEN Gao. Power plants air-cooling technology[M]. Beijing: China Electric Power Press, 2008: 24-26.

［2］ 辛世崗. 火力發(fā)電機組直接空冷系統(tǒng)研究[J]. 能源與節(jié)能, 2012(8): 119-120.

XIN Shigang. Research on the direct air-cooled condenser system of power plants[J]. Energy and Energy Conservation, 2012(8): 119-120.

［3］ 張曉魯, 汪建平, 孫銳, 等. 火電機組直接空冷系統(tǒng)優(yōu)化設(shè)計方法研究[J]. 中國電機工程學(xué)報, 2011, 31(11): 1-5.

ZHANG Xiaolu, WANG Jianping, SUN Rui, et al. Optimal design method of direct air-cooling system for thermal power plants[J]. Proceedings of the CSEE, 2011, 31(11): 1-5.

［4］ 劉煒, 趙蘭萍, 楊志剛. 空冷凝汽器蛇形翅片扁平管空氣側(cè)流動與換熱特性研究[J]. 發(fā)電與空調(diào), 2014(5): 67-71.

LIU Wei, ZHAO Lanping, YANG Zhigang. Study on the air-side flow and heat transfer characteristic of wave-finned flat tube in air-cooling condensers[J]. Power Generation & Air Condition, 2014(5): 67-71.

［5］ KR?GER D G. Air-cooled heat exchangers and cooling towers[J]. Thermal Power Plants, 2005, 3(6): 96.

［6］ 喬英杰, 楊建國, 張兆營, 等. 直接空冷凝汽器單排翅片管外部流動及換熱特性分析[J]. 電站系統(tǒng)工程, 2010, 26(5): 45-48.

QIAO Yingjie, YANG Jianguo, ZHANG Zhaoying, et al. Analysis of extratubal flow and heat transfer performance for finned-tube of direct dry air-cooled condenser[J]. Power System Engineering, 2010, 26(5): 45-48.

［7］ 代善良. 迎面風(fēng)速對直接空冷凝汽器流動換熱特性的影響[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2013:10-11.

DAI Shanliang. Influence of head wind velocity on flow and heat transfer performance of direct air-cooled condenser[D]. Harbin: Harbin Institute of Technology, 2013:10-11.

［8］ 胡漢波. 直接空冷式凝汽器翅片散熱器流動傳熱性能及單元流場特性研究[D]. 重慶: 重慶大學(xué), 2006: 15-24.

HU Hanbo. Research on flow and heat transfer of finned heat exchanger of direct-cooled condenser and flow characteristic of condenser unit[D]. Chongqing: Chongqing University, 2006:15-24.

［9］ 申微波. 直接空冷器空氣側(cè)傳熱和流動的數(shù)值模擬[D]. 武漢: 華中科技大學(xué), 2008: 15-32.

SHEN Weibo. Numerical simulation for heat transfer and flow characteristics in air-side of direct air coolers[D]. Wuhan: Huazhong University of Science and Technology, 2008: 15-32.

［10］ 馬榮榮. 直接空冷凝汽器翅片管空氣流動換熱特性的數(shù)值研究[D]. 保定: 華北電力大學(xué), 2011: 27-33.

MA Rongrong. Numerical study on air flow and heat transfer characteristics for the direct air-cooled condneser fin[D]. Baoding: North China Electric Power University, 2011: 27-33.

［11］ 侯雪梅. 火電廠直接空冷凝汽器蛇形翅片管流動換熱性能數(shù)值研究[D]. 濟南: 山東大學(xué), 2011: 55-74.

HOU Xuemei. Numerical study of finned serpentine tubes flow and heat transfer performance in direct air cooling power plant[D]. Ji’nan: Shandong University, 2011: 55-74.

［12］ 王厚華, 吳偉偉, 李大偉, 等. 穿孔型翅片幾何結(jié)構(gòu)對換熱性能的影響[J]. 同濟大學(xué)學(xué)報(自然科學(xué)版), 2014, 42(12): 1912-1919.

WANG Houhua, WU Weiwei, LI Dawei, et al. Refrigeration performance analysis of fin-tube heat exchangers by effect of perforated fin geometry[J]. Journal of Tongji University (Natural Science), 2014, 42(12): 1912-1919.

［13］ 趙蘭萍, 楊志剛. 管間距對矩形翅片橢圓管換熱管束性能的影響[J]. 同濟大學(xué)學(xué)報(自然科學(xué)版), 2016, 44(1): 150-154.

ZHAO Lanping, YANG Zhigang. Effect of tube pitches on performance of rectangular finned elliptical tube bundled[J]. Journal of Tongji University (Natural Science) 2016, 44(1): 150-154.

［14］ 程遠達. 直接空冷式翅片管及空冷凝汽器單元流動與換熱的數(shù)值研究[D]. 北京: 北京交通大學(xué), 2008: 19-20.

CHENG Yuanda. Numerical investigation on the heat transfer and fluid flow characteristic of finned-tube heat exchanger and direct air cooled condenser unit[D]. Beijing: Beijing Jiaotong University, 2008:19-20.

［15］ MEYER C J, KR?GER D G. Air-cooled heat exchanger inlet flow losses[J]. Applied Thermal Engineering, 2001, 21(7): 771-786.

Influence of interrupted holes and fin pitch on flow and heat transfer of wavy-finned tubes in direct air-cooled condensers

HE Xiaoyi, YIN Jixiang, WANG Tao, MA Suxia

(School of Electrical and Power Engineering, Taiyuan University of Technology, Taiyuan 030024, China)

Interrupted holes are added to the conventional wave-fins of direct air-cooled single-row tube to enhance heat transfer. The three-dimensional numerical simulation on two kinds of finned tube models with and without interrupted holes was performed by Fluent software. The effects of fin pitch on the heat transfer, flow resistance and overall performance of the two finned tubes at four facing velocities (1.4 m/s, 2.2 m/s, 3.0 m/s, 3.8 m/s) were studied. The results show that, compared with the wavy-finned tubes with no holes, the wavy-finned tubes with interrupted holes increase the air-side heat transfer coefficient by 1%~18% at various fin pitches and facing velocities, and the flow resistance slightly reduces. In the smaller facing velocity range, adding interrupted holes can reduce the optimal fin pitch under the same working conditions.

direct air cooling, wavy-finned tube, interrupted hole, fin pitch, heat transfer coefficient, flow resistance, numerical simulation

Science and Technology Major Special Project of Shanxi Province (MD2016-02); National Natural Science Foundation of China(51476108)

賀曉怡(1991—)，女，碩士研究生，主要研究方向為對流換熱強化技術(shù)，1522066082@qq.com。

K26

A

10.19666/j.rlfd.201805152

賀曉怡, 陰繼翔, 王濤, 等. 擾流孔及翅片間距對直接空冷蛇形翅片管流動與換熱的影響[J]. 熱力發(fā)電, 2019, 48(2): 72-77. HE Xiaoyi, YIN Jixiang, WANG Tao, et al. Influence of interrupted holes and fin pitch on flow and heat transfer of wavy-finned tubes in direct air-cooled condensers[J]. Thermal Power Generation, 2019, 48(2): 72-77.

2018-05-13

山西省科技重大專項(MD2016-02)；國家自然科學(xué)基金項目(51476108)

陰繼翔(1964—)，女，副教授，碩士生導(dǎo)師，yinjixiang@tyut.edu.cn。

（責(zé)任編輯 劉永強）