雙開縫翅片空氣側(cè)換熱和流動特性的數(shù)值模擬

魏 雙 鄭傳祥 何建龍 馮苗根

(1.浙江大學(xué)化工機械研究所；2.杭州杭氧換熱設(shè)備有限公司)

雙開縫翅片空氣側(cè)換熱和流動特性的數(shù)值模擬

魏 雙*1鄭傳祥1何建龍2馮苗根2

(1.浙江大學(xué)化工機械研究所；2.杭州杭氧換熱設(shè)備有限公司)

利用數(shù)值耦合傳熱計算方法研究了四排管雙開縫翅片管式換熱器的傳熱特性和阻力特性，獲得翅片附近空氣的速度場和溫度場，分析了開縫對翅片管換熱器換熱性能的影響。并針對翅片間距和開縫高度對翅片管換熱性能的影響進行了數(shù)值模擬分析，得到了開縫高度和翅片間距針對換熱系數(shù)的最佳組合。

換熱器 換熱性能 雙開縫翅片 翅片間距 開縫高度

翅片開縫是一種有效的強化空氣側(cè)傳熱的方法。國內(nèi)外對開縫翅片傳熱與流動阻力進行了大量的實驗和數(shù)值仿真研究。Jeom-Yul Yun和Kwan-Soo Lee針對開縫翅片的翅片間距、縫高、縫長及開縫數(shù)目等7個結(jié)構(gòu)參數(shù)進行了三倍放大模型的正交實驗，研究表明，翅片間距對換熱性能的影響最大[3]。Hiroaki K等對X形排列的開縫翅片進行了實驗比較，表明這種翅片具有較好的換熱性能[4]。Yonghan Kim和Yongchan Kim通過實驗研究得知平直翅片的翅片間距從15.0mm減少到7.5mm時換熱性能降低[5]。由于開縫翅片流動傳熱的復(fù)雜性，其實驗研究具有很大的局限性，因此與開縫翅片結(jié)構(gòu)尺寸參數(shù)有關(guān)的數(shù)值模擬研究仍處于發(fā)展過程中[6]，如屈治國等對3種開縫翅片進行了數(shù)值模擬，并應(yīng)用場協(xié)同原理解釋后開縫翅片比前開縫翅片優(yōu)越的原因[7]。

筆者利用fluent軟件對四排管雙開縫翅片管式換熱器的換熱特性和阻力特性進行數(shù)值模擬研究，比較這種翅片管換熱器相較于平直翅片管換熱器的優(yōu)勢，分析雙開縫翅片的翅片間距和開縫高度對換熱器換熱特性和阻力特性的影響，為進一步的工程實踐應(yīng)用提供參考。

1 數(shù)值模擬

1.1 物理問題描述和控制方程

筆者的模擬對象為翅片管，管內(nèi)的冷卻水與翅片間的熱空氣進行強制對流換熱，翅片和管子之間導(dǎo)熱。根據(jù)換熱器結(jié)構(gòu)的對稱性和周期性，選取兩翅片之間的空間進行模擬，考慮到回流的問題，在流動方向兩端都做相應(yīng)的延長，三維計算區(qū)域選取如圖1虛線所示部分[8]。由于模擬空氣的馬赫數(shù)小于0.3，因此，模型確定為三維不可壓縮、穩(wěn)態(tài)、常物性的紊流控制方程。數(shù)值模擬參數(shù)如下：

翅片間距 1.8、2.0、2.2、2.5、3.0mm

管子外徑 19.5mm

橫向/縱向管間距 25mm

開縫寬度 2mm

開縫高度 1.25、1.50mm

管排數(shù) 4

牛羊用顆粒飼料主要是精料補充料顆粒飼料。與豬雞飼料構(gòu)成不同的是，牛羊精料補充料中常含有較高比例的高纖維原料，如DDGS、麥麩、米糠、大豆皮、苜蓿草粉、甜菜粕等，這些原料吸水性較差，調(diào)質(zhì)目標(biāo)水分值低，因而需要采用較高的蒸汽壓力（0.25～0.4 MPa），較小的蒸汽量，較長的調(diào)質(zhì)時間?？梢圆捎谜{(diào)質(zhì)器（80～85 ℃，20～30 s）+保持器（45～90 s）+制粒的一次制粒工藝。

翅片厚度 0.2mm

圖1 雙開縫翅片管三維計算區(qū)域示意圖

1.2 網(wǎng)格劃分和邊界條件設(shè)置

為了在保證一定計算精度的同時節(jié)省運算時間，在入口段和出口段采用結(jié)構(gòu)化網(wǎng)格，在復(fù)雜的翅片附近采用六面體非結(jié)構(gòu)網(wǎng)格，在管壁處對網(wǎng)格進行加密，翅片表面劃邊界層(圖2)。在檢查網(wǎng)格獨立性后采用網(wǎng)格數(shù)為40～50萬。

圖2 計算模型網(wǎng)格

假設(shè)管壁溫度為常數(shù)288K，翅片溫度由耦合計算得到。耦合傳熱計算設(shè)置管子的導(dǎo)熱系數(shù)為無窮大，以保證管壁溫度為常數(shù)。流體入口溫度為373K，進口風(fēng)速V為1～10m/s。計算區(qū)域邊界條件如圖3所示。

圖3 計算區(qū)域邊界條件

1.3 計算變量

2 模擬結(jié)果與分析

2.1 開縫對翅片管的影響

針對平直翅片和雙開縫翅片，通過數(shù)值模擬得到兩者在進口風(fēng)速V=10m/s時翅片表面附近空氣的速度場和溫度場(圖4、5)。從圖4可以明顯看出，雙開縫翅片相比于平直翅片，管后的間斷縫使尾流區(qū)不同程度地遭到了破壞，空氣的主流通道得到擴展。圖5顯示雙開縫翅片表面高溫區(qū)域的面積比平直翅片表面高溫區(qū)域的面積大，而且前者等溫線密集。這是由于主流通道的擴展使得熱空氣把更多熱量傳遞給翅片表面，使翅片表面溫度更高，從而使翅片管獲取更好的換熱效果。翅片上的開縫結(jié)構(gòu)不僅能加強空氣擾動，擴展空氣的主流通道，還能使翅片空氣的邊界層分層發(fā)展，有效減薄了邊界層厚度，從而起到強化傳熱的作用。

a. 平直翅片 b. 雙開縫翅片

a. 平直翅片 b. 雙開縫翅片

圖6為平直翅片和雙開縫翅片的換熱系數(shù)、阻力系數(shù)曲線?？梢钥闯觯p開縫翅片的換熱系數(shù)比平直翅片高25%，但開縫會導(dǎo)致流動阻力迅速增加，而空氣側(cè)阻力系數(shù)的增加程度低于23%。綜合考慮對換熱系數(shù)和阻力系數(shù)的影響可知，翅片開縫提高了換熱器的換熱性能。

a. 換熱系數(shù)

b. 阻力系數(shù)

2.2 翅片間距和開縫高度對翅片管換熱器性能的影響

當(dāng)換熱器的空氣進口速度較小時，邊界層較厚，翅片間距和開縫高度對換熱性能的影響明顯。平直翅片管翅片間距減小，換熱效率升高；而雙開縫翅片管翅片和縫條的相對高度也對空氣流動邊界層的破壞程度有影響，從而影響整體的換熱效率。

筆者在進口流速V為2～4m/s范圍內(nèi)，對開縫高度為1.50mm，翅片間距分別為2.0、2.2、2.5、3.0mm和開縫高度為1.25mm，翅片間距分別為1.8、2.0、2.2、2.5mm的雙開縫翅片換熱性能進行數(shù)值模擬，計算結(jié)果如圖7所示。

a. 開縫高度1.50mm

b. 開縫高度1.25mm

由圖7可知，除個別特例以外，換熱系數(shù)大體上隨翅片間距的減小而提高。在開縫高度為1.50mm時，翅片間距為2.2mm的換熱系數(shù)最大；開縫高度為1.25mm時，翅片間距為2.0mm的換熱系數(shù)最大。顯然，對于雙開縫翅片來說，最佳換熱系數(shù)的翅片間距與開縫高度有關(guān)，一般取高于開縫高度0.7mm左右。由圖8可看出，隨著翅片間距的減小，翅片的阻力系數(shù)增大。因此，兼顧換熱系數(shù)和阻力系數(shù)，雙開縫翅片管間距的取值應(yīng)考慮開縫高度的影響。

a. 開縫高度1.50mm

b. 開縫高度1.25mm

圖9為進口流速V=3m/s，翅片開縫高度分別為1.25、1.50mm時，換熱系數(shù)和阻力系數(shù)隨翅片間距的變化曲線。從圖9a可以看出，換熱系數(shù)隨翅片間距的變化趨勢受開縫高度的影響；從圖9b可看出，開縫高度越大翅片的阻力系數(shù)越大。

a. 換熱系數(shù)

b. 阻力系數(shù)

3 結(jié)束語

筆者通過對平直翅片和雙開縫翅片運用fluent進行耦合傳熱數(shù)值模擬，得出了兩種翅片附近空氣的速度場和溫度場分布情況；同時針對翅片間距和開縫高度這兩個結(jié)構(gòu)參數(shù)對翅片管換熱性能影響進行了數(shù)值模擬分析，由此可知：翅片開縫能提高管后空氣的流速，使空氣的主流通道得到擴展，在相同的雷諾數(shù)下，雙開縫翅片的換熱系數(shù)比平直翅片高25%，空氣側(cè)阻力系數(shù)的增加低于23%?？偠灾_縫能提高翅片的綜合換熱性能；對于雙開縫翅片，增加翅片間距，阻力系數(shù)會減小，換熱系數(shù)大體呈遞減趨勢，但當(dāng)翅片間距取大于開縫高度0.7mm時,換熱系數(shù)能達到最高；開縫高度對換熱系數(shù)的影響隨翅片間距的不同而不一樣，但開縫高度增大，阻力系數(shù)會增大。

[1] 崔海亭，彭培英.強化傳熱新技術(shù)及其應(yīng)用[M].北京：化學(xué)工業(yè)出版社，2002.

[2] Mattarolo L.Recent Development in Finned Tube Heat Exchangers: Theoretical and Practical Aspect：DTI，Energy Technology，Po Box 141，DK-2630 Taastrup，Denmark[J].International Journal of Refrigeration，1994，17(5)：359～362.

[3] Jeom-Yul Yun，Kwan-Soo Lee.Influence of Design Parameters on the Heat Transfer and Flow Friction Characteristics of the Heat Exchanger with Slit Fins[J].International Journal of Heat and Mass Transfer，2000，43(14)：2529～2539.

[4] Hiroaki K，Shinicki I，Osamu A，et al.High-efficiency Heat Exchanger[J].National Technical Report，1989，35(6)： 653～ 661.

[5] Yonghan Kim，Yongchan Kim.Heat Transfer Characteristics of Flat Plate Finned-tube Heat Exchangers with Large Fin Pitch[J].International Journal of Refrigeration，2005，28(3)：851～858.

[6] Shah R K.Progress in the Numerical Analysis of Compact Heat Exchanger Surfaces[J].Advances in Heat Transfer，2001，40(1)：579～ 599.

[7] 屈治國，何雅玲，陶文銓.平直開縫翅片傳熱特性的三維數(shù)值模擬及場協(xié)同原理分析[J].工程熱物理學(xué)報，2003，24(5)：825～827.

[8] 陶文銓.數(shù)值傳熱學(xué)[M].西安：西安交通大學(xué)出版社，2001.

NumericalSimulationofHeatTransferandFluidFlowCharacteristicsatAir-sideofDouble-slottedFinSurface

WEI Shuang1, ZHENG Chuan-xiang1, HE Jian-long2, FENG Miao-gen2

(1.InstituteofChemicalMachinery,ZhejiangUniversity,Hangzhou310027,China;2.HangzhouHangyangHeatExchangeEquipmentCo.,Ltd.,Hangzhou310013,China)

A numerical investigation of air-side heat transfer and fluid flow characteristics of a four-row finned tube heat exchanger with double-slotted fin surface was implemented to get air velocity and temperature field nearby the fins and to analyze the fin effect on the heat transfer performance of finned tube heat exchanger. Simulating and analyzing both fin pitch and slit height’s influence on the heat exchanger performance shows that the optimal combination of the fin pitch and slit height can be reached to improve heat transfer coefficient.

heat exchanger, heat exchange performance, slotted fin, fin pitch, slit height

TQ051.5

A

0254-6094(2016)01-0072-05

*魏 雙，女，1990年2月生，碩士研究生。浙江省杭州市，310027。

2015-01-29，

2016-01-06)