馮振飛,羅小平,郭 峰,李海燕,王兆濤

(1.華南理工大學(xué)機(jī)械與汽車工程學(xué)院,廣東廣州510641;2.廣西大學(xué)廣西石化資源加工及過程強(qiáng)化技術(shù)重點(diǎn)實(shí)驗(yàn)室,廣西南寧530004)

微柱體對(duì)微通道熱沉綜合性能影響的數(shù)值分析

馮振飛1,2,羅小平1,郭 峰1,李海燕1,王兆濤1

(1.華南理工大學(xué)機(jī)械與汽車工程學(xué)院,廣東廣州510641;2.廣西大學(xué)廣西石化資源加工及過程強(qiáng)化技術(shù)重點(diǎn)實(shí)驗(yàn)室,廣西南寧530004)

為強(qiáng)化微通道熱沉的傳熱性能,設(shè)計(jì)一種內(nèi)置微柱體的微通道熱沉,并采用數(shù)值方法研究微柱體對(duì)微通道熱沉內(nèi)流體流動(dòng)、傳熱及綜合性能的影響。分析了進(jìn)口雷諾數(shù)、微柱體的錯(cuò)位量對(duì)內(nèi)置微柱體微通道熱沉(微柱通道)的壓降、熱阻和努塞爾數(shù)的影響,并與光滑微通道熱沉(光滑通道)進(jìn)行對(duì)比。采用熱阻與泵功的關(guān)系、熵產(chǎn)原則及性能評(píng)價(jià)準(zhǔn)則對(duì)微通道熱沉的綜合性能進(jìn)行評(píng)價(jià)。結(jié)果表明,微柱通道壓降和努塞爾數(shù)隨雷諾數(shù)增大而增大,熱阻反而減小;在研究的雷諾數(shù)范圍內(nèi),微柱通道壓降比光滑通道的平均高84.3%,熱阻平均低27.8%,而努塞爾數(shù)平均高54.5%;有錯(cuò)位量的微柱通道熱阻比無錯(cuò)位量的平均低8.9%,而努塞爾數(shù)平均高12.6%;微柱通道綜合性能優(yōu)于光滑通道,且有錯(cuò)位量的微柱通道更優(yōu)。

微通道;微柱體;數(shù)值模擬;流動(dòng);傳熱

工業(yè)中諸多領(lǐng)域的突破性進(jìn)展是由微型化技術(shù)和先進(jìn)的冷卻技術(shù)推動(dòng)的。在現(xiàn)代高科技中起到重要作用的微芯片是典型的代表。由于其高度集成化和微型化,使得工作時(shí)熱流密度急劇增加,穩(wěn)定性不斷下降。為保證高熱流密度微芯片正常工作,迫切需要高效的散熱技術(shù)和微型化的散熱裝置[1]。微通道熱沉因具有結(jié)構(gòu)緊湊、散熱效率高、功耗低及使用的冷卻液少等優(yōu)點(diǎn)[2-3],自其被Tuckerman等[4]提出以來,學(xué)者們不斷研究以進(jìn)一步強(qiáng)化其傳熱性能。強(qiáng)化微通道熱沉傳熱有很多種方法,如使用強(qiáng)化的工質(zhì)、改進(jìn)微通道的結(jié)構(gòu)形狀等。由于這些方法不需要增加額外的動(dòng)力(除進(jìn)出口泵功外),因此成為學(xué)者們研究的熱點(diǎn)[5]。納米流體屬于強(qiáng)化工質(zhì)中的一種,其較基液具有高導(dǎo)熱率。Azizi等[6]對(duì)圓柱型微通道內(nèi)Cu-水納米流體的對(duì)流傳熱研究發(fā)現(xiàn),納米流體可強(qiáng)化微通道傳熱,強(qiáng)化效果最大可達(dá)23%。然而,Rimbault等[7]對(duì)微通道熱沉內(nèi)Cu-水納米流體的流動(dòng)與傳熱的試驗(yàn)研究發(fā)現(xiàn),低濃度納米流體的傳熱效果稍有增強(qiáng),但是壓降也提高。此外,由于納米流體易發(fā)生團(tuán)聚使納米粒子沉積,致使納米流體在工作時(shí)穩(wěn)定性降低,還有研究者[8]發(fā)現(xiàn),從熱力的觀點(diǎn)來看,使用強(qiáng)化微通道表面(如正弦波形通道)的方法來增強(qiáng)傳熱的效果比使用強(qiáng)化工質(zhì)(如納米流體)的方法更有效,且穩(wěn)定性和可靠性更高。夏國(guó)棟等[9]研究了硅基長(zhǎng)菱形微針肋熱沉的流動(dòng)和傳熱特性。Chai等[10]研究了微通道兩側(cè)壁面帶有不同形狀肋片的層流流動(dòng)與傳熱特性。Yadav等[11]研究了微通道內(nèi)置微圓柱肋的流動(dòng)和傳熱特性。Chen等[12]研究了帶有縱向渦發(fā)生器微通道內(nèi)流體的流動(dòng)和傳熱。筆者采用數(shù)值方法研究微柱體對(duì)微通道熱沉內(nèi)流體流動(dòng)、傳熱和綜合性能的影響,分析不同錯(cuò)位量、不同雷諾數(shù)條件下,微柱體對(duì)微通道熱沉的流動(dòng)壓降、熱阻和努塞爾數(shù)的影響,并從熱阻與泵功的關(guān)系、熵產(chǎn)原則及性能評(píng)價(jià)準(zhǔn)則等方面對(duì)微通道熱沉綜合性能進(jìn)行評(píng)價(jià)。

1 數(shù)值模擬

1.1 物理模型

本研究的基本微通道熱沉結(jié)構(gòu)與文獻(xiàn)[13]中的一致。銅基熱沉共蝕刻有21條微槽,槽寬W=231 μm,槽高H=713 μm,槽長(zhǎng)L=44.764 mm,并蓋上PC(聚碳酸酯)板形成微通道結(jié)構(gòu)。基于熱沉里微通道結(jié)構(gòu)的周期性和計(jì)算機(jī)的運(yùn)算能力,取一微通道單元進(jìn)行模擬,并考慮鄰近通道固體表面的影響。通道內(nèi)設(shè)有10個(gè)微柱體,即為微柱通道(MCC),見圖1(a)。模型的具體幾何參數(shù)見圖1?；诓煌奈⒅w錯(cuò)位量S,可得3種不同規(guī)格微柱通道MCC-0、MCC-1和MCC-2,對(duì)應(yīng)的S值分別為0、60和120 μm。為對(duì)比有無微柱體的微通道熱沉內(nèi)流體的流動(dòng)、傳熱和綜合性能,還模擬了光滑矩形微通道(PMC)。

圖1 微柱通道結(jié)構(gòu)示意圖Fig.1 Structure of micro-channel with micro-cylinders

1.2 數(shù)值方法

微通道內(nèi)的工質(zhì)為去離子水。為簡(jiǎn)化計(jì)算,假設(shè):(1)流體流動(dòng)為單相連續(xù)不可壓縮的穩(wěn)態(tài)層流;(2)模型的熱傳遞為穩(wěn)態(tài),且忽略熱輻射影響;(3)除流體的黏度隨溫度呈分段線性變化外,固體和流體的物性參數(shù)均為常數(shù);(4)忽略體積力。

基于這些假設(shè),可得到流體域的連續(xù)方程、動(dòng)量方程和能量方程分別為

對(duì)于固體域,U=0。固體域的傳熱分析僅需要能量方程

式中,下標(biāo)f和s分別表示流體和固體;U為流體速度矢量,m/s;p為壓力,Pa;T為溫度,K;μ、ρ、cp、λf分別為流體相應(yīng)溫度的黏度、密度、比熱和熱導(dǎo)率,Pa·s、kg/m3、J/(kg·K)和W/(m·K);λs為固體熱導(dǎo)率,W/(m·K)。

模型采用均勻速度入口邊界條件,給定進(jìn)口速度uin分別為1.33、1.99、2.6、3.32、3.98 m/s,入口溫度Tin為288.1 K;通道出口采用壓力出口條件,相對(duì)壓力pout=0 Pa;通道底面恒熱流加熱,熱流密度q為2 MW/m2;兩側(cè)面設(shè)為周期性邊界條件,固液接觸面設(shè)為固液交界面邊界條件,其余壁面絕熱。采用有限體積法離散控制方程組,并由CFD軟件進(jìn)行求解,收斂殘差為10-6。

1.3 網(wǎng)格劃分

鑒于微柱通道幾何結(jié)構(gòu)的復(fù)雜性,采用六面體和四面體混合的網(wǎng)格模式對(duì)計(jì)算域進(jìn)行劃分。在流體域的壁面,特別是固體域和流體域交界處的壁面,進(jìn)行網(wǎng)格細(xì)化處理。為保證計(jì)算精度及節(jié)省計(jì)算時(shí)間和內(nèi)存,對(duì)每種規(guī)格的微通道熱沉單元模型進(jìn)行網(wǎng)格獨(dú)立性驗(yàn)證。各模型的最終網(wǎng)格數(shù)見表1。

表1 各計(jì)算模型的網(wǎng)格單元數(shù)Table 1 Numbers of grid elements for each model

2 結(jié)果及其討論

2.1 模型有效性驗(yàn)證

為驗(yàn)證數(shù)值模型的有效性,將光滑矩形微通道的數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果[13]進(jìn)行對(duì)比。圖2為光滑矩形微通道進(jìn)出口的壓降及溫差的數(shù)值結(jié)果與試驗(yàn)結(jié)果對(duì)比。壓降和溫差計(jì)算式分別為

式中,Δp為進(jìn)出口壓降,Pa;ΔT為進(jìn)出口溫差,K;pin為進(jìn)口壓力,Pa;Tout為出口溫度,K;Ain、Aout分別為進(jìn)出口面積,m2;p(x,y,z)、u(x,y,z)、Tf(x,y,z)分別為流體域坐標(biāo)(x,y,z)處微元體的壓力、速度和溫度,Pa、m/s和K。

圖2中Re為雷諾數(shù),計(jì)算式為

式中,dh為通道的水力直徑,m;μ為流體進(jìn)口溫度對(duì)應(yīng)的黏度,Pa·s。

圖2 數(shù)值結(jié)果與試驗(yàn)結(jié)果對(duì)比Fig.2 Comparison of numerical results with experimental results

由圖2可見,進(jìn)出口壓降和溫差的數(shù)值結(jié)果與試驗(yàn)結(jié)果比較吻合,最大誤差分別為14.8%和2.8%,考慮到試驗(yàn)測(cè)量具有一定的誤差,因此本數(shù)值方法計(jì)算的結(jié)果是可靠的,進(jìn)而驗(yàn)證了本數(shù)值方法的有效性,可以用于微柱通道熱沉模型的計(jì)算。由圖2還可見,進(jìn)出口溫差隨著雷諾數(shù)的增加逐漸減小,且減小的趨勢(shì)變緩。這是因?yàn)樵跓崃髅芏纫欢ǖ臈l件下,由能量平衡得到的進(jìn)出口溫差與流體的流速成反比關(guān)系。進(jìn)出口壓降隨雷諾數(shù)的增加逐漸增大,但增速較平緩。這是因?yàn)榱黧w的黏度隨溫度變化所致。

2.2 壓降特性

流體在通道內(nèi)流動(dòng)時(shí),因流體的黏性而產(chǎn)生流動(dòng)阻力,表現(xiàn)為流體壓力沿程逐漸降低。圖3為3種不同規(guī)格的微柱通道及光滑通道進(jìn)出口壓降隨雷諾數(shù)變化關(guān)系。由圖3可見,在相同雷諾數(shù)下,微柱通道(MCC-0、MCC-1、MCC-2)壓降明顯高于光滑通道(PMC)。在研究的雷諾數(shù)范圍內(nèi),微柱通道壓降比光滑通道平均高84.3%,最高達(dá)112.3%。這是由于微柱體的存在使得通道的流通面積突縮突放,流阻增大,且微柱體后面產(chǎn)生尾渦區(qū),尾渦區(qū)內(nèi)的旋渦消耗有用的機(jī)械能,亦增大流阻。由圖3還可見,隨著雷諾數(shù)的增大,微柱通道與光滑通道壓降的差異逐漸變大,這是由于增大雷諾數(shù)使得尾渦區(qū)內(nèi)的旋渦加劇,從而增大了流阻。微柱通道中MCC-1的壓降最大,MCC-2次之,MCC-0最小。表明微柱通道壓降隨錯(cuò)位量增大先增大而減小。這是由于錯(cuò)位量影響了流體流動(dòng)的路徑及混合程度。錯(cuò)位量為0時(shí),對(duì)應(yīng)為MCC-0通道,沿流動(dòng)方向看,相鄰兩微柱體重疊,流體可同時(shí)從微柱體兩側(cè)繞流。而隨錯(cuò)位量增大,對(duì)應(yīng)為MCC-1通道,沿流動(dòng)方向看,相鄰兩微柱體部分重疊,流體會(huì)選擇流通面積較大一側(cè),即阻力較小一側(cè)流過,使得流體流動(dòng)的路徑呈S字型,增長(zhǎng)流動(dòng)路徑,促進(jìn)流體混合,進(jìn)而增大了流阻。隨錯(cuò)位量進(jìn)一步增大,對(duì)應(yīng)為MCC-2通道,沿流動(dòng)方向看,相鄰兩微柱體有間隙,大部分流體可在通道中間區(qū)域流動(dòng),流動(dòng)路徑減小,流阻下降。

圖3 進(jìn)出口壓降隨雷諾數(shù)變化Fig.3 Variation of pressure drop of inlet and outlet with Reynolds number

2.3 傳熱特性

電子器件散熱必須均勻,以免局部過熱或過冷而產(chǎn)生熱應(yīng)力降低電子器件壽命。對(duì)于使用微通道熱沉散熱的電子器件,要想提高壽命,必須提高微通道熱沉底面溫度的均勻性,這可通過降低微通道熱沉的熱阻實(shí)現(xiàn)。微通道熱沉對(duì)流傳熱過程中總熱阻定義為

式中,RT為總熱阻,K/W;Aw為熱沉底面面積,m2;Tw,max為底面最高溫度,K。

圖4為3種不同規(guī)格的微柱通道及光滑通道熱阻隨雷諾數(shù)變化的關(guān)系。由圖4可見,熱阻隨雷諾數(shù)增加逐漸下降,低雷諾數(shù)時(shí)下降趨勢(shì)稍快,高雷諾數(shù)時(shí)下降趨勢(shì)略緩。在相同雷諾數(shù)下,微柱通道熱阻明顯低于光滑通道,有錯(cuò)位量的微柱通道(MCC-1、MCC-2)熱阻略低于無錯(cuò)位量微柱通道(MCC-0)。在研究的雷諾數(shù)范圍內(nèi),微柱通道熱阻比光滑通道平均低27.8%,有錯(cuò)位量的微柱通道熱阻比無錯(cuò)位量平均低8.9%。表明微柱通道熱沉底面溫度的均勻性優(yōu)于光滑通道熱沉,且有錯(cuò)位量的微柱通道熱沉更優(yōu)。

圖4 熱阻隨雷諾數(shù)的變化Fig.4 Variation of thermal resistance with Reynolds number

微通道的傳熱特性可用努塞爾數(shù)Nu進(jìn)行評(píng)估。Nu定義為

式中,Aif為流固接觸面積,m2;Tw為底面平均溫度,K;Ts(x,y,0)為固體域坐標(biāo)(x,y,0)處微元體的溫度,K。

圖5為3種規(guī)格微柱通道及光滑通道努塞爾數(shù)隨雷諾數(shù)變化的關(guān)系。由圖5可見,努塞爾數(shù)隨雷諾數(shù)增大而增大。在相同雷諾數(shù)下,微柱通道的努塞爾數(shù)明顯高于光滑通道,表明微柱通道的傳熱性能優(yōu)于光滑通道。這是由于微柱體的存在加劇了通道內(nèi)冷熱流體的混合,重新發(fā)展了微柱體附近的邊界層,擴(kuò)展了傳熱面積。在研究的雷諾數(shù)范圍內(nèi),微柱通道努塞爾數(shù)比光滑通道平均高54.5%,有錯(cuò)位量的微柱通道努塞爾數(shù)比無錯(cuò)位量平均高12.6%。由圖5還可見,相同雷諾數(shù)下,微柱通道中MCC-1的努塞爾數(shù)最大,MCC-2次之,MCC-0最小。MCC-1的流動(dòng)路徑近似S形,這可促進(jìn)流體混合,增長(zhǎng)傳熱時(shí)間,進(jìn)而強(qiáng)化了傳熱。MCC-2和MCC-1通道在低雷諾數(shù)時(shí)努塞爾數(shù)相近,且比MCC-0高許多,表明錯(cuò)位微柱體在低雷諾數(shù)時(shí)強(qiáng)化傳熱效果更為顯著。

圖5 努塞爾數(shù)隨雷諾數(shù)的變化Fig.5 Variation of Nusselt number with Reynolds number

2.4 綜合性能評(píng)價(jià)

在實(shí)際工程設(shè)計(jì)中,低阻高效傳熱的微通道熱沉是工程師所追求的。微柱體強(qiáng)化微通道傳熱的同時(shí)也增大了流阻。須對(duì)微柱通道的綜合性能進(jìn)行評(píng)價(jià),以對(duì)其性能的優(yōu)劣進(jìn)行判斷。由于熱阻與泵功的關(guān)系直接反映微通道熱沉運(yùn)行的經(jīng)濟(jì)性,因此采用熱阻與泵功的關(guān)系評(píng)價(jià)微通道熱沉的性能[14]。泵功的定義為

式中,P為泵功,W;qv為體積流量,m3/s。

圖6為3種不同規(guī)格的微柱通道及光滑通道熱阻隨泵功變化的關(guān)系。由圖6可見,在相同泵功下,PMC通道熱阻最大,MCC-0、MCC-1和MCC-2通道熱阻相近,MCC-1和MCC-2通道略低于MCC-0,表明微柱通道熱阻明顯比光滑通道低,有錯(cuò)位量的微柱通道熱阻略低于無錯(cuò)位量的微柱通道。這也表明微柱通道熱沉的經(jīng)濟(jì)性明顯優(yōu)于光滑通道,且有錯(cuò)位量的更優(yōu)些。

微通道熱沉的對(duì)流換熱過程包含熱量傳遞過程和流體流動(dòng)過程,這些不可逆過程都會(huì)引起系統(tǒng)的熵產(chǎn)。熵產(chǎn)越小,系統(tǒng)的有用能損失就越少,即能源的有效利用程度越高。有學(xué)者采用熵產(chǎn)原則評(píng)價(jià)微通道熱沉的性能[15-17]。根據(jù)Bejan[18]的熵產(chǎn)模型可得,因傳熱不可逆引起的局部體積熵產(chǎn)為

式中,ST為傳熱局部體積熵產(chǎn),W/(m3·K)。

因流動(dòng)不可逆引起的局部體積熵產(chǎn)為

式中,SF為流動(dòng)局部體積熵產(chǎn),W/(m3·K);u、v、w為速度矢量U在x、y、z方向的分量,m/s。

圖6 熱阻隨泵功的變化Fig.6 Variation of thermal resistance with pumping power

為了便于比較,定義無量綱總熵產(chǎn)率為[15]

式中,S*為總熵產(chǎn)率,無量綱;V為流體的體積,m3;M為質(zhì)量流量,kg/s。

圖7為3種不同規(guī)格的微柱通道及光滑通道無量綱總熵產(chǎn)率隨雷諾數(shù)的變化關(guān)系。

圖7 總熵產(chǎn)率隨雷諾數(shù)的變化Fig.7 Variation of total entropy generation rate with Reynolds number

由圖7可見,無量綱總熵產(chǎn)率隨雷諾數(shù)增大逐漸減小,表明在研究的雷諾數(shù)范圍內(nèi),因傳熱不可逆引起的熵產(chǎn)占主導(dǎo)地位。這是由于微通道熱沉的壓降和熱阻隨雷諾數(shù)增大分別增大和減小,即流動(dòng)熵產(chǎn)和傳熱熵產(chǎn)隨雷諾數(shù)增大分別增大和減小,只有傳熱熵產(chǎn)占總熵產(chǎn)的比重較大才會(huì)使得總熵產(chǎn)率隨雷諾數(shù)增大逐漸減小。由圖7還可見,在相同雷諾數(shù)下,微柱通道總熵產(chǎn)率明顯比光滑通道低,表明微柱體能改善微通道熱沉的傳熱性能,從而減少在傳熱過程中的有用能損失,而且有錯(cuò)位量的微柱通道的有用能損失低于無錯(cuò)位量的微柱通道,低雷諾數(shù)時(shí)尤為顯著。這也表明了有錯(cuò)位量的微柱通道的傳熱性能優(yōu)于無錯(cuò)位量的。

性能評(píng)價(jià)準(zhǔn)則被廣泛地應(yīng)用于評(píng)價(jià)微通道熱沉的綜合性能[15,17,19]。為此引入綜合評(píng)價(jià)因子,定義為

式中,η為綜合評(píng)價(jià)因子;下標(biāo)0代表光滑微通道,即為PMC通道。

圖8 綜合評(píng)價(jià)因子隨雷諾數(shù)的變化Fig.8 Variation of overall performance factor with Reynolds number

圖8為3種規(guī)格微柱通道綜合評(píng)價(jià)因子隨雷諾數(shù)變化的關(guān)系。由圖8可見,綜合評(píng)價(jià)因子均大于1,最大值達(dá)1.31,表明微柱體對(duì)微通道熱沉性能的提升具有重要性。在整個(gè)雷諾數(shù)的研究范圍內(nèi),MCC-1和MCC-2通道的綜合評(píng)價(jià)因子均大于MCC-0,低雷諾數(shù)時(shí)更顯著。表明有錯(cuò)位量的微柱通道的綜合性能優(yōu)于無錯(cuò)位量的微柱通道,尤其在低雷諾數(shù)時(shí)。由圖8還可見,隨雷諾數(shù)的增加,MCC-0通道的綜合評(píng)價(jià)因子逐漸增加,MCC-1和MCC-2通道的綜合評(píng)價(jià)因子先緩慢增大,后緩慢減小。表明無錯(cuò)位量的微柱通道MCC-0在整個(gè)雷諾數(shù)的研究范圍內(nèi),Nu增量占主導(dǎo),而對(duì)于有錯(cuò)位量的微柱通道MCC-1和MCC-2,低雷諾數(shù)時(shí),Nu增量占主導(dǎo),高雷諾數(shù)時(shí)Nu增量受壓降增量抑制。由此可見,不能一味地使用提高流速的方法強(qiáng)化傳熱。低雷諾數(shù)時(shí)MCC-2通道綜合評(píng)價(jià)因子略高于MCC-1,高雷諾數(shù)時(shí)反而略低,表明高錯(cuò)位量的微柱通道在低雷諾數(shù)時(shí)綜合性能更優(yōu)。

3 結(jié) 論

(1)微柱通道壓降隨雷諾數(shù)增加而增大;同雷諾數(shù)下,微柱通道壓降明顯高于光滑通道,且微柱通道壓降隨錯(cuò)位量增大先增后減。在研究的雷諾數(shù)范圍內(nèi),微柱通道壓降比光滑通道平均高84.3%。

(2)微柱通道熱阻和努塞爾數(shù)隨雷諾數(shù)增加分別減小和增大;在研究的雷諾數(shù)范圍內(nèi),微柱通道熱阻比光滑通道平均低27.8%,而努塞爾數(shù)平均高54.5%;有錯(cuò)位量的微柱通道熱阻比無錯(cuò)位量平均低8.9%,而努塞爾數(shù)平均高12.6%。

(3)采用熱阻與泵功的關(guān)系、熵產(chǎn)原則及性能評(píng)價(jià)準(zhǔn)則對(duì)微通道熱沉的綜合性能進(jìn)行評(píng)價(jià),微柱通道熱沉綜合性能優(yōu)于光滑通道熱沉,且有錯(cuò)位量的微柱通道熱沉更優(yōu)。

[1] 劉用鹿,羅小兵,劉偉.基于擾流的蜂窩層疊微通道熱沉散熱試驗(yàn)研究[J].半導(dǎo)體光電,2010,31(4):563-566.LIU Yonglu,LUO Xiaobing,LIU Wei.Experimental study on the multilayered honeycomb microchannel heat sink design based on the flow disruption[J].Semiconductor Optoelectronics,2010,31(4):563-566.

[2] 王麗鳳,邵寶東,程赫明,等.矩形微通道熱沉內(nèi)單相穩(wěn)態(tài)層流流體的流動(dòng)與傳熱分析[J].應(yīng)用數(shù)學(xué)和力學(xué),2013,34(8):855-862.WANG Lifeng,SHAO Baodong,CHENG Heming,et al.Fluid flow and heat transfer analysis of single-phase steady laminar flow in rectangular micro-channel heatsinks[J].Applied Mathematics and Mechanics,2013,34(8):855-862.

[3] 馮振飛,何榮偉,朱禮,等.周向平行細(xì)通道夾套的換熱特性[J].過程工程學(xué)報(bào),2015,15(6):901-908.FENG Zhenfei,HE Rongwei,ZHU Li,et al.Heat transfer characteristics of a jacket with circumferential parallel minichannels[J].The Chinese Journal of Process Engineering,2015,15(6):901-908.

[4] TUCKERMAN D B,PEASE R E W.High-performance heat sinking for VLSL[J].IEEE Electronic Device Letters,1981,2(5):126-129.

[5] 夏國(guó)棟,翟玉玲,崔珍珍,等.帶有內(nèi)肋的凹穴型微通道熱沉的熵產(chǎn)及傳熱特性分析[J].中國(guó)科學(xué):技術(shù)科學(xué),2013,43(8):860-866.XIA Guodong,ZHAI Yuling,CUI Zhenzhen,et al.Characteristics of entropy generation and heat transfer in a microchannel with fan-shaped reentrant cavities and inter-nal ribs[J].Science China Technological Sciences,2013,43(8):860-866.

[6] AZIZI Z,ALAMDARI A,MALAYERI M R.Convective heat transfer of Cu—water nanofluid in a cylindrical microchannel heat sink[J].Energy Conversion and Management,2015,101:515-524.

[7] RIMBAULT B,NGUYEN C T,GALANIS N.Experimental investigation of CuO—water nanofluid flow and heat transfer inside a microchannel heat sink[J].International Journal of Thermal Sciences,2014,84:275-292.

[8] KHOSHVAGHT-ALIABADI M,SAHAMIYAN M,HESAMPOUR M,et al.Experimental study on cooling performance of sinusoidal—wavy minichannel heat sink[J].Applied Thermal Engineering,2016,92:50-61.

[9] 夏國(guó)棟,崔珍珍,翟玉玲,等.長(zhǎng)菱形微針肋熱沉的流動(dòng)與換熱特性[J].中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,38(2):130-134.XIA Guodong,CUI Zhenzhen,ZHAI Yuling,et al.Flow and heat transfer characteristics for long-diamond shaped micro pin fin[J].Journal of China University of Petroleum(Edition of Natural Science),2014,38(2):130-134.

[10] CHAI L,XIA G D,WANG H S.Numerical study of laminar flow and heat transfer in microchannel heat sink with offset ribs on sidewalls[J].Applied Thermal Engineering,2016,92:32-41.

[11] YADAV V,BAGHEL K,KUMAR R,et al.Numerical investigation of heat transfer in extended surface microchannels[J].International Journal of Heat and Mass Transfer,2016,93:612-622.

[12] CHEN C,TENG J T,CHENG C H,et al.A study on fluid flow and heat transfer in rectangular microchannels with various longitudinal vortex generators[J].International Journal of Heat and Mass Transfer,2014,69:203-214.

[13] QU W,MUDAWAR I.Experimental and numerical study of pressure drop and heat transfer in a single-phase microchannel heat sink[J].International Journal of Heat and Mass Transfer,2002,45(12):2549-2565.

[14] XIE G N,SHEN H,WANG C C.Parametric study on thermal performance of microchannel heat sinks with internal vertical Y-shaped bifurcations[J].International Journal of Heat and Mass Transfer,2015,90:948-958.

[15] WANG Y S,ZHOU B,LIU Z C,et al.Numerical study and performance analyses of the mini-channel with discrete double-inclined ribs[J].International Journal of Heat and Mass Transfer,2014,78:498-505.

[16] SHI Z Y,DONG T.Entropy generation and optimization of laminar convective heat transfer and fluid flow in a microchannel with staggered arrays of pin fin structure with tip clearance[J].Energy Conversion and Management,2015,94:493-504.

[17] CHAI L,XIA G D,WANG H S.Parametric study on thermal and hydraulic characteristics of laminar flow in microchannel heat sink with fan-shaped ribs on sidewalls—part 3:performance evaluation[J].International Journal of Heat and Mass Transfer,2016,97:1091-1101.

[18] BEJAN A.Entropy generation minimization[M].New York:CRC Press,1996.

[19] LI Y L,ZHANG F L,SUNDEN B,et al.Laminar thermal performance of microchannel heat sinks with constructal vertical Y-shaped bifurcation plates[J].Applied Thermal Engineering,2014,73(1):185-195.

(編輯 沈玉英)

Effects of micro-cylinders on overall performance of micro-channel heat sinks:numerical simulation

FENG Zhenfei1,2,LUO Xiaoping1,GUO Feng1,LI Haiyan1,WANG Zhaotao1
(1.School of Mechanical and Automotive Engineering,South China University of Technology,Guangzhou 510641,China;2.Guangxi Key Laboratory of Petrochemical Resource Processing and Process Intensification Technology,Guangxi University,Nanning 530004,China)

In order to enhance the heat transfer performance of micro-channel heat sinks,a micro-channel heat sink with micro-cylinders is presented and the effects of micro-cylinders on the fluid flow,heat transfer and overall performance in the micro-channel heat sinks are studied numerically.The effects of the inlet Reynolds number and micro-cylinder offset value on the pressure drop,thermal resistance and Nusselt number are analyzed,and these results are compared with those of the plain micro-channel heat sink.Then the overall performance on the micro-channel heat sink is evaluated by using the thermal resistance-pumping power relation,entropy generation principle and performance evaluation criterion.It is found that the pressure drop and Nusselt number increase and the thermal resistance decreases with the increase of the Reynolds number.For the studied Reynolds number range,there are the 84.3%average increase in the pressure drop,a 27.8%average decrease in the thermal resistance and a 54.5%average increase in the Nusselt number for the micro-channel heat sink with micro-cylinders,in comparison with the plain micro-channel heat sink.There are a 8.9%average decrease in the thermal resistance and a 12.6%average increase in the Nusselt number for micro-channel heat sinks with the micro-cylinder offsetvalue being nonzero,compared with that with the micro-cylinder offset value being zero.The overall performance in the micro-channel with micro-cylinders is much better than that of the plain micro-channel.Furthermore,the overall performance in the micro-channel with micro-cylinder offset value being nonzero is much better.

micro-channel;micro-cylinder;numerical simulation;flow;heat transfer

TK 124

:A

馮振飛,羅小平,郭峰,等.微柱體對(duì)微通道熱沉綜合性能影響的數(shù)值分析[J].中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2017,41(3):122-128.

FENG Zhenfei,LUO Xiaoping,GUO Feng,et al.Effects of micro-cylinders on overall performance of micro-channel heat sinks:numerical simulation[J].Journal of China University of Petroleum(Edition of Natural Science),2017,41(3):122-128.

1673-5005(2017)03-0122-07doi:10.3969/j.issn.1673-5005.2017.03.015

2016-06-17

國(guó)家自然科學(xué)基金項(xiàng)目(21276090);廣西自然科學(xué)基金項(xiàng)目(2014GXNSFBA118051);廣西石化資源加工及過程強(qiáng)化技術(shù)重點(diǎn)實(shí)驗(yàn)室主任基金項(xiàng)目(2015Z012)

馮振飛(1980-),男,講師,博士研究生,研究方向?yàn)槲⒓{米尺度傳熱。E-mail:zffeng@gxu.edu.cn。

羅小平(1967-),男,教授,博士,博士生導(dǎo)師,研究方向?yàn)閺?qiáng)化傳熱與節(jié)能。E-mail:mmxpluo@scut.edu.cn。