縱向渦發(fā)生器攻角對(duì)翅片傳熱性能的影響

王成剛，劉俊，劉慧，張博

武漢工程大學(xué)機(jī)電工程學(xué)院，湖北 武漢 430205

縱向渦發(fā)生器攻角對(duì)翅片傳熱性能的影響

王成剛，劉俊，劉慧，張博

武漢工程大學(xué)機(jī)電工程學(xué)院，湖北 武漢 430205

為了提高翅片式換熱器的強(qiáng)化傳熱性能，對(duì)不同攻角的縱向渦發(fā)生器的翅片區(qū)域進(jìn)行了研究.采用流體仿真軟件FLUENT對(duì)發(fā)生器的翅片區(qū)域建立六面體網(wǎng)格模型，對(duì)不同攻角的縱向渦發(fā)生器的努塞爾數(shù)、阻力因子、綜合性能分別進(jìn)行數(shù)值分析和對(duì)比.結(jié)果顯示：隨著縱向渦發(fā)生器的攻角增加，縱向渦翅片的努塞爾數(shù)增強(qiáng)的越來(lái)越明顯，其中縱向渦發(fā)生器攻角為45°的翅片努塞爾數(shù)最大；同時(shí)隨著縱向渦發(fā)生器的攻角增加，阻力因子也隨之增加；帶縱向渦發(fā)生器的翅片的傳熱效果強(qiáng)于不帶縱向渦發(fā)生器的翅片；通過(guò)比較綜合評(píng)價(jià)因子，攻角為45°的縱向渦翅片在雷諾數(shù)為2000～6000內(nèi)，綜合性能最好.

翅片換熱器；強(qiáng)化傳熱；攻角；縱向渦

0 引言

為了提高翅片式換熱器的強(qiáng)化傳熱性能，將研究一種翅片擾流裝置：矩形縱向渦發(fā)生器.由于帶縱向渦發(fā)生器的幾何結(jié)構(gòu)復(fù)雜，影響縱向渦發(fā)生器翅片傳熱特性的因素比較多，做實(shí)驗(yàn)來(lái)研究不僅成本高、耗時(shí)長(zhǎng).采用數(shù)值模擬可以簡(jiǎn)化這一過(guò)程，而且現(xiàn)在的數(shù)值計(jì)算技術(shù)已經(jīng)可以很高效地模擬翅片中流場(chǎng).目前對(duì)縱向渦發(fā)生器的研究已經(jīng)普遍，但是對(duì)于其參數(shù)的影響卻說(shuō)法不一.所以本文采用矩形縱向渦發(fā)生器在一種管翅式換熱器中進(jìn)行數(shù)值模擬.得出矩形縱向渦發(fā)生器在這種換熱器中的最佳攻角參數(shù).

1 數(shù)值計(jì)算模型

1.1 物理模型

針對(duì)4排管換熱器作為研究對(duì)象，單個(gè)翅片基本結(jié)構(gòu)呈正方形結(jié)構(gòu)，換熱管錯(cuò)開(kāi)排列.如圖1所示翅片寬度為140 mm，間距2.2 mm，換熱管直徑為16.3mm，其它參數(shù)如表1所示.因?yàn)槌崞w結(jié)構(gòu)為循環(huán)對(duì)稱結(jié)構(gòu)，所以只取其一部分進(jìn)行研究［1－2］.主要采取控制變量的方法進(jìn)行數(shù)值研究.通過(guò)改變矩形縱向渦的發(fā)生器的攻角來(lái)研究矩形縱向渦發(fā)生器對(duì)翅片式換熱器傳熱特性的影響.

圖1 模型圖Fig.1model figure

表1 翅片換熱器參數(shù)Table.1 The parameters of the Finned heat exchanger

1.2 計(jì)算域控制方程

模型的假設(shè)條件：流體不可壓縮、穩(wěn)態(tài)流動(dòng)，忽略流體粘性耗散、忽略重力和流體升浮力，換熱管溫度恒定為常數(shù)，且不考慮接觸熱阻和輻射放熱.其控制方程如下：

連續(xù)性方程為

動(dòng)量方程為

其中，μ為流體動(dòng)力粘度.

能量守恒方程為

式（3）中a為流體熱擴(kuò)散率.固體導(dǎo)熱方程為

對(duì)流方程為

式（5）中，i，j＝（x，y，z），ρ為氣體密度，λ1，λ2分別為管壁以及氣體導(dǎo)熱系數(shù)；T為溫度，K；p為壓力，Pa.

設(shè)定通過(guò)管翅式換熱器的流體速度為2~9 m／s，流動(dòng)大部分為湍流.所以在FLUENT中選擇k-ε模型計(jì)算.

標(biāo)準(zhǔn)k-ε模型方程

式（7）中：k是湍動(dòng)能，ε為湍流耗散率，Gk是由平均速度梯度引起的湍動(dòng)能k的產(chǎn)生項(xiàng)，Gb升浮力引起的湍動(dòng)能k的產(chǎn)生項(xiàng)，ui為液相速度分量，μt為湍流粘性系數(shù)，σk和σε分別是與湍動(dòng)能和耗散率相關(guān)的普朗特準(zhǔn)數(shù)，σk＝1.0，σε=1.3；Cμ，C1ε，C2ε為經(jīng)驗(yàn)常數(shù)，在FLUENT中通常取0.09、1.44和1.92［3］.

1.3 邊界條件

流體為理想空氣，且忽略重力和升浮力.計(jì)算域入口設(shè)為速度入口（velocity inlet），翅片換熱器中氣體流速一般很低，所以取入口為1～2 m／s，且溫度設(shè)為293 K.上下表面以及延長(zhǎng)段上下表面都設(shè)為周期邊界（periodic）.翅片的厚度相當(dāng)于間距的大約1／10，不考慮其溫度在翅片厚度方向上的變化，設(shè)為無(wú)滑移絕熱壁面（wall）.換熱管管壁設(shè)為無(wú)滑移壁面（wall），且溫度設(shè)為300 K；翅片兩側(cè)設(shè)為對(duì)稱邊界條件（symmetry）［4－6］.

1.4 劃分模型網(wǎng)格

采用Spaceclaim建立模型，并使用ANSYSICEM劃分模型網(wǎng)格.計(jì)算結(jié)果的好壞很大程度上取決于網(wǎng)格質(zhì)量的優(yōu)劣，而ICEM網(wǎng)格劃分工具對(duì)模型有更好的結(jié)構(gòu)適應(yīng)性，能夠?yàn)榉墙Y(jié)構(gòu)體模型劃分結(jié)構(gòu)型網(wǎng)格.這樣不僅節(jié)約計(jì)算資源也很大地減少計(jì)算所用的時(shí)間.模型都是采用結(jié)構(gòu)六面體網(wǎng)格，如圖2所示.采用局部加密的措施保證計(jì)算的可靠性.

圖2 縱向渦區(qū)域Fig.2Area of longitudinal vortex

縱向渦區(qū)域?yàn)橹饕芯繉?duì)象，所以將其所在區(qū)域網(wǎng)格加密.網(wǎng)格從下翅片向上增長(zhǎng)率為1.2，且設(shè)定邊界層網(wǎng)格，因?yàn)樵诹魉贋?.8 m／s時(shí)，流體將從層流轉(zhuǎn)為湍流.且根據(jù)估算第一層網(wǎng)格高度為0.01 mm.總網(wǎng)格數(shù)量為10萬(wàn)左右.

數(shù)據(jù)整理與結(jié)果分析：

總換熱量：q=qmCp(T0-Ti)

總換熱系數(shù)：h=q/(AΔT)

對(duì)數(shù)平均溫度：ΔT=(Tmax-Tmin)/ln(Tmax/Tmin)

雷諾數(shù)：Re=umaxd3/v

努塞爾數(shù)：Nu=hd3/λ

其中，Ti、T0分別為進(jìn)出口流體的平均溫度；q為入口流體的質(zhì)量流量；Tmax=max(Twi-Ti,Two-To)，Tmin=min(Twi-Ti,Two-To)分別為進(jìn)出口處的最大和最小溫度差值；umax為最小截面流速；d3為特征長(zhǎng)度d3=d2+2δ，d2為換熱管外徑，δ為翅片厚度；L為翅片長(zhǎng)度；ΔP為壓降.

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

圖3為Nu不同攻角縱向渦翅片下比較.帶縱向渦發(fā)生器的翅片改變縱向渦的攻角，其傳熱性能也有所不同.圖3中帶wu標(biāo)簽的是光翅片.可以看出隨著雷諾數(shù)增加，傳熱系數(shù)不斷增強(qiáng).其中縱向渦攻角為45°的換熱性能最優(yōu).而光翅片的換熱性能最差.在不同的雷諾數(shù)下，攻角為15°的縱向渦翅片的Nu比光翅片平均提高8%左右，根據(jù)上面的原理分析和圖的流線布局來(lái)分析：這是因?yàn)榭v向渦發(fā)生器的夾角過(guò)小，導(dǎo)致垂直縱向渦發(fā)生器的速度太小，形成的縱向渦強(qiáng)度不夠，并沒(méi)有很明顯的增強(qiáng)傳熱.隨著攻角的增加，縱向渦翅片的Nu增強(qiáng)的越來(lái)越明顯，但這種趨勢(shì)在低雷諾數(shù)下體現(xiàn)的并不明顯，在雷諾數(shù)低于2 000的情況下，由于該情況屬于層流，湍流度不強(qiáng)，導(dǎo)致這幾種不同攻角的Nu都在4～6之間，這與何雅玲和楚攀等［7］人的實(shí)驗(yàn)結(jié)果一致.

圖3 不同攻角縱向渦翅片Nu比較Fig.3Comparison of Nu with different attack Angle longitudinal vortex fin

圖4為六種不同翅片阻力因子f的比較圖，從圖中可以看出，縱向渦迎流攻角為75°的縱向渦翅片阻力因子f最大.不帶縱向渦的翅片阻力因子f最小.原因是縱向渦發(fā)生器的形狀阻力，使得整個(gè)翅片的阻力因子f變大.且隨著攻角的增加阻力因子f逐漸增加.當(dāng)縱向渦發(fā)生器迎流攻角為75°時(shí)，該發(fā)生器幾乎與來(lái)流方向垂直，所以其阻力因子f最大.而光翅片沒(méi)有增加縱向渦，所以其阻力因子f最小.從圖中可以看出在攻角從45°增加到60°之間，阻力因子差距比其他幾組之間的差距要大的多，說(shuō)明在45°和60°之間流體的湍流度有較大的增加.

圖4 不同攻角縱向渦翅片f比較Fig.4Comparison of f with different attack Angle longitudinal vortex fin

流速增加，雷諾數(shù)也同時(shí)增大，流體的湍流程度增加，邊界層與主流中心的混合更加充分，傳熱增強(qiáng)，但也導(dǎo)致阻力因子增加.考慮到綜合效能，所以考慮Nu/f1/3與雷諾數(shù)之間的關(guān)系.從圖5中可以看出，45°攻角的縱向渦翅片在雷諾數(shù)為2 000～6 000綜合性能最好，在6 000～8 000之間略有下降，其中光翅片的綜合性能最差，有縱向渦的翅片綜合性能要強(qiáng)于光翅片.這是由于引入縱向渦使得翅片間流體的湍流度增加，主要熱阻區(qū)域的熱邊界層與主流區(qū)之間形成擾動(dòng)，破壞熱邊界層厚度，增強(qiáng)傳熱.從圖5中也可以看出，在高雷諾數(shù)下（Re＞6 500以上區(qū)域）15°、30°和45°翅片的綜合效能隨著雷諾數(shù)的增加會(huì)明顯下降.說(shuō)明高雷諾數(shù)下縱向渦發(fā)生器的影響受到限制，甚至?xí)蓴_傳熱.

圖5 不同攻角縱向渦翅片Nu/f1/3與Re的關(guān)系Fig.5Nu/f1/3and Re in different attack Anglelongitudinal vortex fin

3 結(jié)語(yǔ)

通過(guò)對(duì)不帶縱向渦的翅片和攻角分別是15°、30°、45°、60°、75°的縱向渦發(fā)生器的傳熱效果進(jìn)行數(shù)值分析和對(duì)比顯示：

（1）帶縱向渦發(fā)生器的翅片的傳熱效果強(qiáng)于不帶縱向渦發(fā)生器的翅片.

（2）通過(guò)對(duì)不帶縱向渦的翅片和攻角分別是15°、30°、45°、60°、75°的縱向渦發(fā)生器的傳熱效果進(jìn)行數(shù)值分析和對(duì)比，攻角為45°的縱向渦翅片具有更好的傳熱效果.

（3）隨著縱向渦發(fā)生器的攻角增加，縱向渦翅片的努塞爾數(shù)增強(qiáng)的越來(lái)越明顯，其中縱向渦發(fā)生器攻角為45°的翅片努塞爾數(shù)最大.

（4）隨著縱向渦攻角的增加，阻力因子逐漸增強(qiáng)，在45°～60°翅片之間出現(xiàn)了較大的差距，45°之后隨著角度增加阻力增加程度很大.

（4）六種翅片的綜合效能Nu/f1/3比較，45°翅片在雷諾數(shù)2 000～6 000之間具有最好的傳熱效能.

（5）在高雷諾數(shù)下（Re＞6 500以上區(qū)域）15°、30°和45°翅片的綜合效能隨著雷諾數(shù)的增加會(huì)明顯下降.說(shuō)明高雷諾數(shù)下縱向渦發(fā)生器的影響受到限制，甚至?xí)蓴_傳熱.

致謝

感謝國(guó)家自然科學(xué)基金委員會(huì)和武漢工程大學(xué)研究生創(chuàng)新基金的資助！

［1］王成剛，鄭曉敏，王小雨，等.圓弧型組合開(kāi)縫翅片三維數(shù)值模擬［J］.化工機(jī)械，2012（1）：62－65.

WANG Cheng－gang，ZHENG Xiao-min，WANG Xiaoyu，et al.3D Numerical simulation on arc-type combined slotted fins［J］.Chemical Engineering＆Machinery，2012（1）：62－65.（in chinese）

［2］王成剛，鄭曉敏，謝小恒，等.4種圓弧型開(kāi)縫翅片特性數(shù)值模擬［J］.石油化工設(shè)備，2012.4（3）：29－33

WANG Cheng-gang，ZHENG Xiao-min，XIE Xiaoheng et al.Numerical simulation on characteristics of four Arc-type slotted fin［J］.Petro-Chemical Equipment，2012，4（3）：29－33（in chinese）

［3］Wang Chenggang，Zheng Xiaomin，Wang Xiaoyu，et al.Numerical simulation for heat transfer performance off four Arc-type slotted fin［C］//Proceedings 2011 International Conference on Materials for Renewable Energy＆Environment，2011：1197－1201.

［4］陽(yáng)祥，王良璧.縱向渦強(qiáng)度衰減及其干涉的數(shù)值分析［J］.工程熱物理學(xué)報(bào)，2006，27（5）：845－846.

YANG Xiang,WANG Liang-Bi.Numerrical analysis of intensity decay and interaction of longitudinal vortices［J］.JournalofEngineeringThermophysics，2006.27（5）：845－846.（in chinese）

［5］葉秋玲，周國(guó)兵，程金明，等.矩形通道中不同渦流發(fā)生器對(duì)換熱和壓降的影響［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2010，30（11）：86－91.

YE Qiu-ling，ZHOU，Guo-bing，CHENG Jin-ming，et al.Influence of different vortex generators on heat transfer enhancement and pressure drop characteristics in a rectangular channel［J］.Proceedings of the CSEE，2010，30（11）：86－91.（in chinese）

［6］李惠珍，屈治國(guó)，程永攀，等.開(kāi)縫翅片流動(dòng)和傳熱性能的實(shí)驗(yàn)研究及數(shù)值模擬［J］.西安交通大學(xué)學(xué)報(bào)，2005，39（3）：229－232.

LI Hui-zhen，Qu Zhi-guo，Cheng Yong-pan，et al.Experimental and numerical study on heat transfer and fluid Flow characteristics of slotted fin2 and 2tube heat transfer surfaces［J］.Journal of Xi’an Jiaotong University，2005，39（3）：229－232.（in chinese）

［7］何雅玲，楚攀，謝濤.縱向渦發(fā)生器在管翅式換熱器中的應(yīng)用與優(yōu)化［J］.化工學(xué)報(bào).2012，3（3）：746－760.

HE Ya-ling，CHU Pan，XIE Tao.Applition and optimization of fin-and-tube heat exchangers with longitudinal vortex generators［J］.CIESC Journal，2012，3（3）：746－760.（in chinese）

Influence of attack angle of longitudinal vortex generator on fin heat transfer performance

WANG Cheng－gang，LIU Jun，LIU Hui，ZHANG Bo
School of Mechanical and Electrical，Wuhan Institute of Technology，Wuhan 430205，China

To improve the heat transfer performance of fin heat exchanger，the fin areas of longitudinal vortex generator with different attack angles were studied.Hexahedral mesh model of fin area of longitudinal vortex generator was established by fluid simulation software FLUENT，Nusselt number，Resistance Factor and comprehensive performance of longitudinal vortex generator with different attack angles were numerically analyzed and compared.The results show that the Nusselt number of longitudinal vortex fin increases with the attack angles of longitudinal vortex generator increasing，and it is biggest at 45°attack angle；the friction factor increases with the attack angles of longitudinal vortex generators increasing；the heat transfer effect of fin with longitudinal vortex generator is better than that without longitudinal vortex generator；the comprehensive performance of longitudinal vortex fins is best when the attack angle is 45°in 2000－6000 Reynolds number.

finned heat exchanger；heat transfer；attack angle；longitudinal vortex

TQ051.5

A

10.3969/j.issn.1674-2869.2015.01.012

本文編輯：陳小平

1674－2869（2015）01－0054－04

2014-09-08

國(guó)家自然科學(xué)基金（50976080）；武漢工程大學(xué)科學(xué)研究基金

王成剛（1974－），男，湖南祁東人，副教授，博士.研究方向：應(yīng)力分析、計(jì)算機(jī)測(cè)控技術(shù)、PLC控制、新型化工設(shè)備的研究和開(kāi)發(fā).