板翅式換熱器中新型微凸翅片結(jié)構(gòu)設(shè)計及性能分析

彭 翔 李登洪 李吉泉 姜少飛

1.浙江工業(yè)大學(xué)機(jī)械工程學(xué)院，杭州，3100232.浙江大學(xué)流體動力與機(jī)電系統(tǒng)國家重點(diǎn)實驗室，杭州，310027

0 引言

板翅式換熱器具有比表面積大、質(zhì)量小、溫差小等優(yōu)點(diǎn)，已在大型空分裝備、航空航天、核工程等行業(yè)得到了廣泛應(yīng)用[1-2]。翅片是板翅式換熱器中進(jìn)行冷熱流體熱量交換的關(guān)鍵部件。伴隨著板翅式換熱器超大型化和低能耗化的發(fā)展趨勢，挖掘翅片區(qū)域的換熱潛能、提高換熱效率成為了板翅式換熱器發(fā)展的關(guān)鍵。

典型翅片結(jié)構(gòu)包括平直翅片、波紋翅片、鋸齒翅片等。翅片結(jié)構(gòu)參數(shù)與出入口工況參數(shù)協(xié)調(diào)作用，共同影響板翅式換熱器的熱功率、壓降、換熱效率等綜合性能指標(biāo)。國內(nèi)外通過計算流體力學(xué)(computational fluid dynamics，CFD)仿真分析[3-4]、風(fēng)洞測試與粒子圖像測速(particle image velocimetry，PIV)等實驗手段[5]、理論模型計算[6]等方式，獲取了翅片結(jié)構(gòu)參數(shù)與換熱性能間的耦合關(guān)聯(lián)關(guān)系，并進(jìn)行了翅片結(jié)構(gòu)參數(shù)的優(yōu)化設(shè)計。SAHA等[7]研究了換熱器中局部流動與換熱強(qiáng)化之間的關(guān)系，利用場協(xié)同原理確定了換熱器內(nèi)部局部傳熱增強(qiáng)機(jī)理。王威[8]通過風(fēng)洞測試研究了三角波紋翅片和雙尺度鋸齒翅片的流動與傳熱性能，研究發(fā)現(xiàn)，兩種結(jié)構(gòu)都能在保證壓降增高幅度很小的情況下實現(xiàn)傳熱強(qiáng)化。

在常規(guī)翅片結(jié)構(gòu)上進(jìn)行結(jié)構(gòu)參數(shù)優(yōu)化，雖然能夠改善翅片的換熱性能，但受限于翅片結(jié)構(gòu)本身的局限性，換熱潛能難以完全挖掘，因此，需要對翅片結(jié)構(gòu)進(jìn)行改進(jìn)，使其結(jié)構(gòu)更有利于破壞流道內(nèi)流體的熱邊界層，并增強(qiáng)流體的湍流性能?，F(xiàn)有的改進(jìn)方法主要分為以下兩種。一種是將不同傳統(tǒng)翅片的優(yōu)勢結(jié)合在一起形成組合翅片結(jié)構(gòu)，如仿生導(dǎo)流增效翅片結(jié)構(gòu)[9]、帶凹坑和擾流柱的翅片結(jié)構(gòu)[10]、拓?fù)鋬?yōu)化確定的非常規(guī)翅片結(jié)構(gòu)[11-12]等。IBRAHIM等[13]分析了穿孔形狀對穿孔翅片傳熱性能的影響，研究發(fā)現(xiàn)，圓形穿孔的傳熱性能優(yōu)于矩形穿孔和三角穿孔的傳熱性能。劉景成等[14]提出了一種仿生導(dǎo)流增效翅片結(jié)構(gòu)，通過促使相鄰流道內(nèi)的流體相互流通來加強(qiáng)流道內(nèi)的擾流，從而提高了換熱器的換熱性能。另一種是在流道內(nèi)增設(shè)渦流發(fā)生輔助裝置，使流道內(nèi)產(chǎn)生渦流，讓流體混合更均勻，如三角形渦發(fā)生器[15]、螺旋曲面渦旋發(fā)生器[16]、翼片式渦流發(fā)生器[17]等。WANG等[18]提出了一種波浪肋結(jié)構(gòu)，提高了渦輪葉片中冷卻流道的流動特性和傳熱性能。喻茹[19]通過在流道內(nèi)增設(shè)V形導(dǎo)流結(jié)構(gòu)，改善了鋸齒錯列翅片換熱器的傳熱效率。

現(xiàn)有的組合翅片和渦流輔助裝置都在一定程度上增大了翅片尺寸，從而導(dǎo)致流阻增大、壓降升高，并且結(jié)構(gòu)優(yōu)化過程均較為復(fù)雜。此外，擾流柱、渦流發(fā)生裝置等新增結(jié)構(gòu)也使換熱器成本顯著增加。本文綜合兩種翅片結(jié)構(gòu)改進(jìn)方式的優(yōu)點(diǎn)，在平直翅片基礎(chǔ)上提出了一種新型微凸翅片結(jié)構(gòu)，通過微凸體以較小的流阻增高幅度來改善翅片流層內(nèi)部的橫向流動，從而提高板翅式換熱器的整體換熱性能。

1 新型微凸翅片結(jié)構(gòu)

1.1 微凸翅片結(jié)構(gòu)模型構(gòu)建

板翅式換熱器常用的翅片形式是平直翅片，翅片內(nèi)部流阻小，但翅片區(qū)域流道形式不變、缺乏橫向流動，從而導(dǎo)致流動不均勻性高、換熱效率低。借鑒組合翅片的換熱優(yōu)勢和平直翅片的低阻特性，在平直翅片的壁面上開縫，沖壓出犁型微凸，形成一種新型微凸翅片結(jié)構(gòu)，如圖1a所示，在盡量保持平直翅片的低阻流動特性的同時，依靠微凸體形成的局部橫向流動來改善換熱性能。翅片壁面上的微凸結(jié)構(gòu)朝向兩邊流道呈交叉狀周期分布，當(dāng)流體流過微凸段時，一部分沿著微凸體的壁面流動，一部分穿過微凸孔與相鄰流道的流體混合，如圖1b所示，起到了流動邊界層破壞及流體橫向混合的效果，使換熱器單流道內(nèi)邊緣區(qū)域和中心區(qū)域的溫度差減小，橫向方向上各流道間的溫度分布也更加均勻。微凸翅片的主要結(jié)構(gòu)參數(shù)如圖2所示，包括翅片間距fc、翅片高度hc、翅片厚度δc、微凸高度hw、微凸開口長度Bw、微凸夾角αw、微凸厚度δw。

(a) 結(jié)構(gòu)示意圖 (b) 內(nèi)部流動示意圖圖1 新型微凸翅片F(xiàn)ig.1 New micro convex fin

圖2 微凸翅片主要結(jié)構(gòu)參數(shù)Fig.2 Main structure parameters of micro convex fin

1.2 微凸翅片結(jié)構(gòu)性能分析

在FLUENT軟件中對新型微凸翅片的換熱性能進(jìn)行仿真分析。首先建立圖3所示的簡化仿真模型，包括入口面、出口面、周期對稱面、上下隔板面等6個邊界面，翅片部分為固體計算域，邊界面與翅片圍成的空間為流體計算域。FLUENT仿真分析使用的計算模型、出入口條件、收斂條件如下。

圖3 仿真模型Fig.3 Simulation model

(1)計算模型選用RNGk-ε湍流模型，采用分離變量法隱式求解，采用二階迎風(fēng)格式求解動量、能量以及湍流參量，壓力與速度耦合選用SIMPLE算法。

(2)只進(jìn)行單層翅片的流動仿真分析，翅片材料為鋁，翅片中流體為理想不可壓縮空氣。

(3)壁面條件：上下隔板面設(shè)置為恒溫300 K的壁面，翅片為耦合壁面，左右側(cè)面設(shè)置為周期對稱面。

(4)出入口條件：入口為14 m/s的速度入口，湍流耗散率為5.8%，水力直徑為3.3 mm，溫度為273 K；出口為壓力出口、背壓101 325 Pa。

(5)收斂條件：能量方程的殘差絕對值達(dá)到10-4，其他方程組的殘差絕對值達(dá)到10-3，進(jìn)出口流體質(zhì)量守恒，系統(tǒng)內(nèi)能量守恒，且出口處流體壓力達(dá)到穩(wěn)定。

在FLUENT中進(jìn)行仿真分析，得到了新型微凸翅片區(qū)域的溫度分布及整體熱功率，以及沿流體流動方向上不同截面位置處的溫度分布。同時還進(jìn)行了相同尺寸下平直翅片的換熱性能仿真分析，距入口面距離為Z的不同截面處的溫度對比如圖4所示。從圖4中可以看出，平直翅片流道內(nèi)存在邊界規(guī)則的熱邊界層，靠近翅片壁面處的溫度顯著高于流道中心溫度，并且隨著流體的流動，熱邊界層厚度不斷增大。對于本文所提出的新型微凸翅片結(jié)構(gòu)，微凸結(jié)構(gòu)的擾流作用使得熱邊界層得到了破壞，且溫度場分布也得到了有效的改善。由整體熱功率的對比結(jié)果可知，微凸翅片的熱功率為7.672 W，與平直翅片的熱功率(6.477 W)相比，提高了18.45%。為提高微凸翅片的換熱性能，需要進(jìn)行單個微凸體結(jié)構(gòu)參數(shù)、翅片上微凸體分布的優(yōu)化設(shè)計，以獲取最優(yōu)的微凸翅片結(jié)構(gòu)參數(shù)。

圖4 平直翅片和微凸翅片在不同截面位置處的溫度分布云圖Fig.4 Temperature contours of straight fin and micro convex fin at different positions of cross-sections

2 翅片微凸體參數(shù)優(yōu)化設(shè)計

為了與平直翅片進(jìn)行性能對比，設(shè)置本文所提出的新型微凸翅片結(jié)構(gòu)的參數(shù)(翅片間距fc、翅片高度hc、翅片厚度δc)與常用平直翅片尺寸一致，分別為fc=2.8 mm、hc=4.7 mm、δc=0.2 mm。取翅片區(qū)域沿流體流動方向的總長度為L=70 mm，進(jìn)行微凸體結(jié)構(gòu)參數(shù)hw、Bw、αw的優(yōu)化設(shè)計。

考慮到軋制過程中翅片變薄這一特點(diǎn)，微凸高度hw不宜過高，但如果hw過低，所起到的擾流作用以及流道間的流通性都將下降，因此取微凸高度hw∈[0.8 mm，1.2 mm]；為簡化仿真模型，取微凸厚度δw為定值(即δw=0.1 mm)。當(dāng)微凸高度δw確定后，微凸夾角αw的大小將決定單個微凸的長度，αw越小，單個微凸在流體流動方向上的長度則越大，這將導(dǎo)致流體流過微凸的時間延長，流動變化平緩，從而會削弱擾流效果，反之αw越大，微凸長度越小，雖然能在一定程度上強(qiáng)化擾動，但勢必會造成壓降的大幅增高，綜合分析后，取微凸夾角αw∈[15°，35°]。微凸開口長度Bw受到翅片高度和微凸高度的制約，Bw不能大于4.3 mm，且考慮到后續(xù)網(wǎng)格劃分時網(wǎng)格質(zhì)量的問題，在開口兩端需各留出0.3 mm的間隙，但微凸開口長度過小、微凸高度較大又與實際制造不相符合，所以取單層微凸開口長度Bw∈[2.5 mm，3.7 mm]。

結(jié)構(gòu)參數(shù)hw、Bw、αw分別取最大值、最小值、中值三個水平，按照L9(33)正交表生成初步采樣點(diǎn)，并隨機(jī)添加另外9組采樣點(diǎn)形成了采樣方案，同時還添加了4組參考組采樣點(diǎn)，與前面18組采樣點(diǎn)中任意3組構(gòu)成誤差分析樣本。根據(jù)各個采樣點(diǎn)的結(jié)構(gòu)參數(shù)值，在FLUENT中構(gòu)建幾何模型、劃分網(wǎng)格、進(jìn)行仿真計算，讀取翅片部分的熱功率R，以及讀取入口和出口面的壓力以求得壓降值Δp，具體結(jié)果如表1所示。

表1 采樣點(diǎn)處仿真分析結(jié)果

根據(jù)采樣點(diǎn)仿真結(jié)果得到的熱功率R和壓降Δp數(shù)據(jù)，利用二階響應(yīng)面法建立了微凸結(jié)構(gòu)參數(shù)與R、Δp的近似模型，其計算表達(dá)式分別如下：

(1)

Δp=-9.2873αw-359.1559hw-47.3048Bw+

8.3974αwhw+1.9107αwBw+117.6960hwBw-

(2)

通過各試驗點(diǎn)數(shù)據(jù)的交叉驗證，可得到基于近似模型的預(yù)測值與真實值的相關(guān)系數(shù)分別為0.9938(熱功率R)和0.9979(壓降Δp)?；跇?gòu)建的近似模型，以熱功率R最高為優(yōu)化目標(biāo)、壓降Δp不超過550 Pa為約束，并考慮加工制造限制(微凸角度αw取整數(shù)，微凸高度hw和微凸開口長度Bw保留小數(shù)點(diǎn)后一位小數(shù))，在Isight軟件中采用序列二次規(guī)劃法進(jìn)行了微凸體結(jié)構(gòu)參數(shù)的優(yōu)化設(shè)計，得到的最優(yōu)微凸體結(jié)構(gòu)參數(shù)及其仿真結(jié)果如表2所示。

表2 最優(yōu)微凸體結(jié)構(gòu)參數(shù)值及仿真結(jié)果

3 含微凸體的新型翅片結(jié)構(gòu)參數(shù)優(yōu)化

3.1 排列間距

翅片上微凸分布的疏密程度直接影響翅片的傳熱效率。為保證相鄰流道的流量相近，需使不同朝向開口的微凸數(shù)量相等。微凸排列間距通過微凸組數(shù)g確定，g=1表示分別向左和向右開口所組成的一組微凸。翅片區(qū)域總長度為70 mm，微凸組數(shù)g會受到翅片長度的限制，g過小時，微凸數(shù)量很少，翅片區(qū)域內(nèi)部流體的橫向流動少，起到的強(qiáng)化傳熱效果差；g過大時，會造成微凸之間的間距很小，且結(jié)構(gòu)制造困難、壓降損失也會很大。綜上，取微凸組數(shù)g分別為3、4、5、6，對應(yīng)的微凸間距fw依次為10.5 mm、8.0 mm、6.7 mm、5.7 mm。不同微凸組數(shù)下的微凸分布情況如圖5所示。

圖5 不同微凸組數(shù)下翅片的結(jié)構(gòu)Fig.5 Fin structures under different number of micro convex

不同排列間距下的仿真結(jié)果如表3所示。當(dāng)微凸組數(shù)g從3增加到6時，R的增長值依次為0.26 W、0.22 W、0.15 W，增長幅度逐漸減小，主要原因是隨著流體在流道內(nèi)流動，流體靜態(tài)溫度T在不斷地升高，在翅片前半段，翅片與流體的溫差較大，傳熱率高，翅片后半段二者溫差較小，傳熱率低。從圖5和圖6中可以看出，當(dāng)微凸數(shù)量越多時，第一個微凸的位置距離入口處越近、微凸間距越小，流道內(nèi)流體升溫越快，高傳熱率流段逐漸變短、低傳熱率流段逐漸變長，并且從Z=70 mm位置處可以觀察到流體出口溫度升高的空間越來越小，從而出現(xiàn)熱功率R增幅逐漸減小的現(xiàn)象。Δp的增長值依次為53.72 Pa、55.97 Pa、58.95 Pa，增長幅度均在緩慢增大，從圖7中可以看出，只有在經(jīng)過每個微凸時，壓降才會發(fā)生明顯變化，這就意味著增加微凸數(shù)量必然會導(dǎo)致壓降的持續(xù)增高。綜上可知，基于壓降不超過550 Pa的約束，選擇g=4、fw=8.0 mm作為最佳排列間距，實際應(yīng)用中可以根據(jù)具體的壓降條件來選擇微凸組數(shù)，從而確定出相應(yīng)的微凸排列間距。

表3 微凸翅片不同排列間距下的仿真結(jié)果

圖6 不同微凸組數(shù)下流道中心線上的溫度分布Fig.6 Temperature distribution on centerline of fin channel with different number of micro convex

圖7 不同微凸組數(shù)下流道中心線上的壓力分布Fig.7 Pressure distribution on centerline of fin channel with different number of micro convex

3.2 排列方式

微凸凸起朝向的不同會對流道內(nèi)流體的流動狀態(tài)造成不同影響，為此，本文建立了交叉排列、對稱排列、雙交叉排列這三種主要的微凸排列方式結(jié)構(gòu)，如圖8所示。交叉排列每組微凸(4個)中處于流道同一縱向位置的兩個微凸的朝向相同并與另外兩個微凸的朝向相反，對稱排列翅片上的微凸關(guān)于流道中心面呈對稱分布且單側(cè)翅片上相鄰兩個微凸的朝向相反，雙交叉排列每組微凸的朝向相同但相鄰微凸組的朝向相反。

(a) 交叉排列

每種微凸排列方式下翅片區(qū)域熱功率和壓降的仿真結(jié)果如表4所示。相較于交叉排列方式，對稱排列方式的熱功率有所提高，其壓降也大幅增大，雙交叉排列方式的熱功率最低，其壓降卻大于交叉排列方式的壓降，因此可以判斷雙交叉排列方式最不可取，對壓降要求寬松的工況可以考慮對稱排列方式，本文以壓降不超過550 Pa為約束，因此選用交叉排列方式。

表4 三種微凸排列方式下的仿真結(jié)果

4 最優(yōu)微凸翅片結(jié)構(gòu)的性能驗證

通過微凸體參數(shù)優(yōu)化和微凸排列優(yōu)化，可確定的微凸翅片結(jié)構(gòu)參數(shù)為：微凸高度hw=1.2 mm、微凸夾角αw=17°、微凸開口長度Bw=3.6 mm、微凸組數(shù)g=4、微凸間距fw=8.0 mm，微凸排列方式為交叉排列。本文從溫度場分析、壓力場分析、與平直翅片換熱性能對比三個角度來驗證新型微凸翅片的有效性。

4.1 溫度場分析

圖9為流道水平中心面和垂直中心面上流體溫度分布云圖。由水平中心面溫度云圖(圖9a)可以直觀地看出流道內(nèi)流體溫度的變化過程，每經(jīng)過一組微凸，流體溫度都有較為明顯的升高，從顏色區(qū)域的分布可以看出前半段顏色變化顯著(即溫度變化幅度大)， 后半段溫度變化減緩(即溫度變化幅度減小)。由垂直中心面溫度云圖(圖9b)可以看出，微凸對流道中間部位產(chǎn)生了擾動作用，使位于流道中間部位的流體熱邊界層不斷被破壞，流體混合更為均勻，從而起到強(qiáng)化傳熱的作用。

(a) 流道水平中心面

圖10給出了流體流動方向上不同截面位置處的溫度分布云圖，其中Z=0表示入口截面，Z=16 mm、32 mm、48 mm、64 mm分別表示流體流過第1、2、3、4組微凸位置邊緣處的截面位置，Z=70 mm表示出口截面，從Z=0到Z=16 mm，流體的熱邊界層被微凸破壞，溫度分布情況發(fā)生了明顯變化；從Z=16 mm到Z=64 mm，流體在微凸后方形成了一個相對的高溫區(qū)域，并且這個區(qū)域在逐漸增大，這是因為流過微凸的流體流量在增加，流動混合更加均勻，傳熱效率更高；由出口截面位置處的溫度云圖可以看出，流體在流出流道時仍未形成明顯的熱邊界層，這表明微凸結(jié)構(gòu)有利于改善流道內(nèi)的溫度場分布。

(a) Z=0 (b) Z=16 mm (c) Z=32 mm

4.2 壓力場分析

圖11為翅片表面壓力分布云圖，圖中紅色圈出部分為微凸位置，流體流至微凸頂部時壓力最大，流過微凸后，壓力又明顯降低，流道內(nèi)壓力是以微凸位置為基準(zhǔn)呈階梯狀降低。這是因為每個微凸頂部都是流體流動空間發(fā)生突變的臨界位置，沿著流動方向，流動空間先減小后增大，空氣在這個過程中產(chǎn)生一個先壓縮后膨脹的變化，在壓縮階段，空氣對翅片的壓力逐漸升高，在膨脹階段，空氣對翅片的壓力瞬間降低，空氣每流過一個微凸時都會重復(fù)上述過程，這也解釋了空氣流動壓降隨微凸間距的減小而不斷增大的原因。

圖11 微凸翅片表面壓力分布云圖Fig.11 Pressure contours of micro convex fin

4.3 與平直翅片的性能對比

本文所提出的新型微凸翅片是在常用的平直翅片基礎(chǔ)上改進(jìn)得到的。為了驗證該新型翅片結(jié)構(gòu)的有效性，構(gòu)建了與微凸翅片結(jié)構(gòu)參數(shù)一樣的平直翅片結(jié)構(gòu)，在所有仿真參數(shù)設(shè)置一樣的情況下，對兩種翅片結(jié)構(gòu)的出口溫度進(jìn)行了對比分析。

圖12為仿真后平直翅片和微凸翅片出口面上的溫度分布圖，與平直翅片相比，微凸翅片出口面上最低溫度基線大幅度上升，溫度極差減??；從線條密集程度上可以看出，微凸翅片流道中間部位的溫度分布得到了極大的改善。

(a) 平直翅片

圖13為平直翅片和微凸翅片的出口面溫度定量分析統(tǒng)計圖，可以看出，兩種翅片出口面的溫度都主要集中在290 K以上的區(qū)域，并且微凸翅片幾乎沒有溫度低于288 K的區(qū)域；微凸翅片出口溫度在292 K以上占比為87.04%，遠(yuǎn)高于平直翅片的68.24%。

圖13 平直翅片和微凸翅片出口面上不同溫度區(qū)間比例對比Fig.13 Comparison of ratio of different temperature ranges on outlet surface of straight fin and micro convex fin

依據(jù)圖13中兩種翅片的溫度區(qū)間分布概率，可計算得出平直翅片和微凸翅片的出口面平均溫度分別為293.60 K和295.18 K。由CFD仿真結(jié)果計算可得，平直翅片中翅片熱功率為6.477 W，壓降為290.62 Pa；微凸翅片對應(yīng)的熱功率和壓降依次為8.282 W、 544.35 Pa。

對比結(jié)果顯示，在以一定的壓力損失為代價后，微凸翅片的傳熱性能得到了提高，從而驗證了其對板翅式換熱器性能改善的有效性。

5 結(jié)語

本文在平直翅片結(jié)構(gòu)基礎(chǔ)上提出了一種新型微凸翅片結(jié)構(gòu)。進(jìn)行了微凸體參數(shù)和翅片結(jié)構(gòu)參數(shù)的優(yōu)化設(shè)計，使用序列二次規(guī)劃法，以熱功率最高為優(yōu)化目標(biāo)、壓降不超過550 Pa為約束，確定了微凸翅片的最佳結(jié)構(gòu)參數(shù)。通過對微凸翅片流道內(nèi)溫度場和壓力場進(jìn)行分析，并與相同結(jié)構(gòu)的平直翅片的傳熱性能進(jìn)行對比，驗證了微凸結(jié)構(gòu)對改善板翅式換熱器性能的有效性。

本文著重從仿真分析與結(jié)構(gòu)優(yōu)化的角度進(jìn)行了新型微凸翅片結(jié)構(gòu)的參數(shù)設(shè)計，后續(xù)將開展實物制備及實驗，以進(jìn)一步驗證提出的新型翅片結(jié)構(gòu)的有效性。