分-合式微通道流動(dòng)傳熱性能數(shù)值模擬

蔣靜智 劉華峰 安京州

摘 要：針對(duì)矩形微通道進(jìn)出口壓降大、溫度分布不均勻，以及分形微通道受到分形維數(shù)和分支數(shù)限制適用范圍較窄的問(wèn)題，結(jié)合矩形微通道和分形微通道的優(yōu)勢(shì)設(shè)計(jì)一種分-合式微通道散熱器。使用Fluent軟件對(duì)散熱過(guò)程進(jìn)行數(shù)值模擬，研究微通道內(nèi)分支傾斜角度變化對(duì)流動(dòng)和傳熱性能的影響。結(jié)果表明，在100 W/cm²的熱流密度下，Re為970、分支傾斜角度為90°時(shí)，分-合式微通道平均溫度降低了11.9 K，最高溫度降低了14.2 K，Nu增加了85.7%，整體傳熱性能（PEC）也最佳，達(dá)到1.44。分支的引入可以增加微通道內(nèi)部換熱面積，同時(shí)形成新的邊界層，在分支內(nèi)側(cè)產(chǎn)生漩渦，有效提高了微通道散熱器的傳熱性能，為微通道的優(yōu)化設(shè)計(jì)提供了新的理論依據(jù)。

關(guān)鍵詞：傳熱學(xué)；微通道；微型設(shè)備；數(shù)值模擬；傳熱

中圖分類(lèi)號(hào)：TK172

文獻(xiàn)標(biāo)識(shí)碼：A DOI：10.7535/hbkd.2023yx03003

收稿日期：2023-04-27；修回日期：2023-05-05；責(zé)任編輯：馮 民

基金項(xiàng)目：國(guó)家自然科學(xué)基金（51706058）；河北省自然科學(xué)基金（E2019208345）

第一作者簡(jiǎn)介：蔣靜智（1972—），女，河北阜城人，教授，博士，主要從事強(qiáng)化傳熱及儲(chǔ)能方面的研究。E-mail：jjzhi2000@126.com

Numerical simulation of flow and heat transfer performance in branch and combined-branch microchannel

JIANG Jingzhi，LIU Huafeng，AN Jingzhou

（School of Mechanical Engineering， Hebei University of Science and Technology， Shijiazhuang， Hebei 050018， China）

Abstract：The pressure drop at the inlet and outlet of rectangular microchannel is large， and the temperature distribution is uneven; Fractal microchannel is limited by fractal dimension and the number of branches， so its application scope is narrow. A branch and combined-branch microchannel heat sink was designed by combining the advantages of rectangular microchannel and fractal microchannel， and Fluent software was used to numerically simulate its heat dissipation process to study the flow and heat transfer performance when the branch tilt angle inside the microchannel changes. The results show that at the heat flow of 100 W/cm²， when the Re is 970 and the branch tilt angle is 90°， the average temperature of the branch and combined-branch microchannel decreases by 11.9 K， the maximum temperature decreases by 14.2 K， the Nu increases by 85.7%， and the performance evaluation criterion （PEC） is also the best by reaching 1.44. The introduction of branches can increase the heat transfer area inside the microchannel， form a new boundary layer， and generate vortices at the inner side of the branch， effectively improving the heat transfer performance of the microchannel heat sinks， which provides new theoretical basis for optimizing the design of microchannels.

Keywords：heat transfer; microchannel; micro device; numerical simulation; heat transfer

隨著電子工業(yè)技術(shù)的迅速崛起，電子元器件正在向集成化、高頻化方向發(fā)展，電路板熱流密度也在不斷突破新高。當(dāng)電子元器件工作時(shí)，如果熱量無(wú)法被及時(shí)散出，內(nèi)部溫度就會(huì)逐漸上升，從而嚴(yán)重影響自身的性能和可靠性^［1-2］。研究表明：電子元器件的可靠性對(duì)溫度十分敏感，當(dāng)工作溫度達(dá)到70～80 ℃后，每上升1 ℃，可靠性就會(huì)下降5%，溫度過(guò)高已成為現(xiàn)在大多數(shù)設(shè)備失效的主要原因^［3］。微通道散熱器具有體積小、散熱快、反應(yīng)時(shí)間短等優(yōu)勢(shì)，在解決微型設(shè)備器件的高熱負(fù)荷問(wèn)題方面極具應(yīng)用前景^［4］。

近年來(lái)，矩形微通道^［5-8］和分形樹(shù)狀微通道^［9-11］的結(jié)構(gòu)優(yōu)化研究領(lǐng)域出現(xiàn)了很多優(yōu)秀研究成果，為微通道的設(shè)計(jì)提供了良好的理論指導(dǎo)。TUCKERMAN等^［12］首次提出了微通道散熱器，發(fā)現(xiàn)微尺度下流體的散熱效果要比宏觀(guān)尺度好得多。陳孝根等^［13］研究了不同高寬比矩形微通道的散熱性能，結(jié)果表明，微通道的高寬比為2～4時(shí)最好。GUNNASEGARAN等^［14］對(duì)不同截面形狀微通道的傳熱性能進(jìn)行了實(shí)驗(yàn)驗(yàn)證，結(jié)果表明，截面形狀對(duì)微通道的傳熱性能有較大的影響。朱崎峰等^［15］對(duì)不同凹槽形狀微通道進(jìn)行實(shí)驗(yàn)驗(yàn)證，結(jié)果表明，微通道凹槽為扇形且水力直徑較小時(shí)換熱系數(shù)高，凹槽為三角形且水力直徑較大時(shí)換熱系數(shù)高。CHAI等^［16］在傳統(tǒng)平直微通道結(jié)構(gòu)的側(cè)壁上添加扇形肋，結(jié)果表明，這種結(jié)構(gòu)的微通道散熱器綜合性能有了很大提升。PENG等^［17］驗(yàn)證了分形樹(shù)狀微通道的流動(dòng)性能，結(jié)果表明，樹(shù)狀結(jié)構(gòu)可以有效降低流動(dòng)阻力。HUANG等^［18］對(duì)變截面分形樹(shù)狀微通道散熱器的傳熱性能進(jìn)行研究，發(fā)現(xiàn)變截面結(jié)構(gòu)可以進(jìn)一步提高散熱器的傳熱性能。WANG等^［19］使用3D打印技術(shù)制造了Y型和H型樹(shù)狀結(jié)構(gòu)散熱器，發(fā)現(xiàn)H型結(jié)構(gòu)散熱器的傳熱性能較好，Y型結(jié)構(gòu)散熱器的壓降較小。PENCE^［20］分析了分形樹(shù)狀微通道散熱器轉(zhuǎn)角結(jié)構(gòu)和交叉結(jié)構(gòu)對(duì)流動(dòng)和傳熱性能的影響，結(jié)果表明，在轉(zhuǎn)角結(jié)構(gòu)和交叉結(jié)構(gòu)區(qū)域，流體流動(dòng)擾動(dòng)增強(qiáng)并破壞了流動(dòng)邊界層，從而提高微通散熱器的換熱能力。

矩形微通道散熱器結(jié)構(gòu)的改進(jìn)可以提高傳熱性能，但是進(jìn)出口壓降也在增加，并且沿流動(dòng)方向流體的溫度會(huì)逐漸上升。分形樹(shù)狀微通道散熱器在具有更高換熱性能和較為均勻的溫度分布的同時(shí)，具有較低的壓降和泵功消耗，但是分形樹(shù)狀微通道散熱器受到分形維數(shù)和分支數(shù)的限制，其入口段換熱面積小，換熱效率較低，且結(jié)構(gòu)形狀被限制，導(dǎo)致其適用性下降。針對(duì)上述問(wèn)題，設(shè)計(jì)一種分-合式微通道散熱器^［21］，并對(duì)微通道內(nèi)的流體流動(dòng)和傳熱性能進(jìn)行模擬研究。

1 微通道數(shù)值模擬

1.1 物理模型

圖1為矩形微通道散熱器示意圖。該散熱器由40根截面為矩形的微通道組成，熱源基板的尺寸為20 mm ×20 mm，熱流密度為100 W/cm²。微通道的流道高度H為400 μm，寬W為300 μm，微通道上下壁厚δ₁均為50 μm，微通道之間的壁厚δ₂為200 μm。圖2為分-合式微通道散熱器示意圖。微通道內(nèi)分支傾斜角度為θ；分支水力直徑表現(xiàn)為寬度變化，傾斜分支的寬度相等為W₁，通道分流后寬度減小為W₂，合流后寬度增加為W₃，故微通道壁厚略有變化；微組分支與合支的距離為h。引入矩形微通道是為了以矩形微通道作為基準(zhǔn)微通道散熱器，與分-合式微通道散熱器進(jìn)行對(duì)比，由于微通道散熱器中微通道的排列具有周期性，其傳熱性能也呈現(xiàn)周期性分布，為了減少模擬仿真的計(jì)算時(shí)間和計(jì)算成本，以及能夠完整反映分-合式微通道結(jié)構(gòu)的周期性，選取5根微通道作為仿真模型。

MURRAY^［22］通過(guò)研究血管的流阻特性，得出父血管半徑的立方和等于子血管半徑的立方總和，即：

D³=D³₁+D³₂+…，（1）

其中，D為血管的水力直徑。當(dāng)微通道的分支寬度W₁和分流后主流道寬度W₂相等時(shí)，根據(jù)式（1）計(jì)算可得各流道的寬度，分支和合支均勻排布，計(jì)算可得W₁=170 mm，W₂=170 mm，W₃=360 mm，h=115 mm。以分支的傾斜角度為變量，設(shè)計(jì)了5種微通道結(jié)構(gòu)，分支傾斜角度編號(hào)及角度如表1所示。

1.2 數(shù)學(xué)模型及控制方程

為簡(jiǎn)化模型，在模擬過(guò)程中作如下假設(shè)：

1）微通道散熱器的傳熱和流體在流道中的流動(dòng)為穩(wěn)態(tài)過(guò)程，且流體流動(dòng)過(guò)程為充分發(fā)展的層流流動(dòng)（微尺度下Re小于1 000時(shí)，流動(dòng)為層流）；

2）假設(shè)不受熱輻射和重力的影響；

3）流體流動(dòng)為充分發(fā)展的層流流動(dòng)，故不考慮流道內(nèi)軸向傳導(dǎo)和黏性耗散；

4）微通道內(nèi)固體材料為硅，其物性參數(shù)如表2所示，流體材料為去離子水。

根據(jù)上述假設(shè)并依據(jù)傳熱學(xué)^［23］，本模型采用的控制方程如下。

連續(xù)性方程：

動(dòng)量方程：

能量方程：

對(duì)于微通道的固體區(qū)域，能量方程為

式中：i，j=1，2，3；下標(biāo)f表示流體，s表示固體；ρ為密度，kg/m³；u為速度，m/s；p為壓力，Pa；μ為流體動(dòng)力黏度，Pa·s；c_p為定壓比熱容，J/（kg·K）；T為溫度，K；λ為導(dǎo)熱率，W/（m·K）。

模型的邊界條件是：散熱器入口為速度入口邊界條件，流體入口溫度為293 K；出口設(shè)置為壓力出口邊界條件，出口面壓力為1.013×10⁵Pa，即表示壓力為0；流體與固體接觸面為耦合邊界條件，無(wú)速度滑移；散熱器底面為恒定熱流密度邊界，100 W/cm²；其余面為絕熱邊界條件。

流動(dòng)和傳熱的評(píng)價(jià)指標(biāo)如下。

通過(guò)壓降的計(jì)算結(jié)果以及模擬結(jié)果可以求得流體的摩擦阻力系數(shù)f，摩擦阻力系數(shù)的計(jì)算公式如下：

式中：D_c為通道當(dāng)量直徑，m；u_m為平均流速，m/s；L為通道的長(zhǎng)度，m；Δp為通道壓降，Pa。

對(duì)流換熱系數(shù)反映了流體與固體換熱能力的強(qiáng)弱，平均對(duì)流傳熱系數(shù)h的計(jì)算公式為

式中：A為微通道基底面面積，m²；A_b為微通道固液交界面面積，m²；T_f，ave為流體平均溫度，K。

Nu表征了微通道的換熱能力，Nu計(jì)算公式為

為綜合考慮微通道散熱器的傳熱特性與摩擦阻力，采用整體傳熱性能（performance evaluation criterion，PEC）來(lái)驗(yàn)證微通道散熱器的綜合能力，其計(jì)算公式為

式中：Nu₀和f₀為基準(zhǔn)微通道散熱器計(jì)算值；Nu和f為分-合式微通道計(jì)算值。

1.3 網(wǎng)格設(shè)置

采用基于有限元方法進(jìn)行模擬計(jì)算的Fluent軟件進(jìn)行模擬仿真，采用SIMPLEC算法用來(lái)求解連續(xù)性、動(dòng)量和能量方程，當(dāng)能量方程的殘差值小于10^-8，其他方程的殘差值小于10^-5時(shí)，模擬收斂，計(jì)算完成。有限元方法的核心是通過(guò)網(wǎng)格劃分的方式將復(fù)雜的幾何實(shí)體離散成有限個(gè)規(guī)則的幾何單元，網(wǎng)格的質(zhì)量以及數(shù)量直接決定結(jié)果的準(zhǔn)確性。使用微通道分支傾斜角度為45°時(shí)的模型進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證，采用4種不同網(wǎng)格數(shù)目進(jìn)行對(duì)比，驗(yàn)證結(jié)果如表3所示。網(wǎng)格數(shù)目最多時(shí)的值F₀和其他網(wǎng)格數(shù)目所得值F之間的相對(duì)誤差e定義為

如表3所示，隨著網(wǎng)格數(shù)目的不斷增加，計(jì)算結(jié)果之間的誤差也在不斷減小，當(dāng)網(wǎng)格數(shù)目為5 334 608時(shí)，Nu、摩擦阻力系數(shù)和熱阻的誤差分別為0.40%，0.11%和0.08%。為了兼顧模擬計(jì)算的準(zhǔn)確性與計(jì)算機(jī)性能，采用網(wǎng)格數(shù)目為5 334 608時(shí)的網(wǎng)格劃分方式進(jìn)行后續(xù)的模擬計(jì)算。

1.4 數(shù)值計(jì)算結(jié)果有效性驗(yàn)證

STEINKE等^［24］對(duì)微通道中的單相流體摩擦系數(shù)進(jìn)行研究，微通道截面為矩形、流體流動(dòng)為層流時(shí)的對(duì)流傳熱摩擦系數(shù)的理論計(jì)算公式為

f Re=96（1-1.355 3α_c+1.946 7α²_c-1.701 2α³_c+0.956 4α⁴_c-0.253 7α⁵_c），（11）

其中，α_c為通道的寬高比。

圖3為微通道散熱器摩擦阻力系數(shù)對(duì)比曲線(xiàn)。從圖3可以看出，當(dāng)流體的物性參數(shù)為常數(shù)時(shí)，摩擦阻力系數(shù)的模擬值與理論值相差較大；當(dāng)流體為變物性時(shí)，其物性參數(shù)根據(jù)參考文獻(xiàn)［23］選取，從圖3可以看出，流體物性參數(shù)隨溫度變化時(shí)，摩擦系數(shù)的模擬值與理論值很接近，最大差值為2.8%，因此微通道內(nèi)的流體采用變物性。由于變物性模擬結(jié)果與理論公式計(jì)算結(jié)果相差很小，從而驗(yàn)證了模擬結(jié)果的可靠性。

2 結(jié)果與討論

對(duì)6種不同模型進(jìn)行模擬計(jì)算后，從流動(dòng)特性、壓力特性以及溫度特性3方面對(duì)不同結(jié)構(gòu)微通道散熱器進(jìn)行分析，以驗(yàn)證微組的引入對(duì)傳熱性能的影響。

2.1 流動(dòng)特性

當(dāng)Re=322時(shí)，微通道x-y截面的速度云圖如圖4所示（z=0.25 mm，0 mm≤x≤2.5 mm，7.5 mm≤y≤12.5 mm）。

從圖4可以看出，微通道貼近壁面處流體的流速幾乎為0，且越靠近流道中心，其流速越大。矩形微通道由于流道形狀沒(méi)有變化，流動(dòng)狀態(tài)一致，對(duì)流速影響較小。對(duì)于分-合式微通道散熱器，當(dāng)流體沒(méi)有流經(jīng)分支時(shí)，流體在微通道的流動(dòng)狀態(tài)與矩形微通道內(nèi)的流動(dòng)狀態(tài)相似；當(dāng)流體流經(jīng)分支時(shí)，由于局部流道變寬，流體的流速有明顯的下降趨勢(shì)，且在分支處有新的邊界層形成；當(dāng)流體流經(jīng)合支時(shí)，由于局部流道變窄，流體的流速有明顯的上升趨勢(shì)，在合支處也有新的邊界層形成。流體每流經(jīng)一處分支或合支時(shí)，其流動(dòng)狀態(tài)的變化相似。

同時(shí)從圖4可以看出，分支傾斜角度的變化對(duì)微通道分支處的流體流動(dòng)有著較大的影響，對(duì)主流道區(qū)流體的流動(dòng)影響相對(duì)較小。隨著分支傾斜角度的不斷增加，流道內(nèi)的最大流速在不斷變大，并且隨著分支傾斜角度的不斷增加，分支處的流速也出現(xiàn)小幅的上升，并且都形成新的邊界層。當(dāng)分支傾斜角度為90°時(shí)，截面處的最大流速與其他角度時(shí)的最大流速相比最大，為2.24 m/s。

當(dāng)Re=322時(shí)，微通道x-y截面的速度矢量圖如圖5所示（z=0.25 mm，0.25 mm≤x≤0.75 mm，12.9 mm≤y≤13.65 mm）。

從圖5可以看出，矩形微通道的流體速度呈線(xiàn)性分布，幾乎無(wú)擾動(dòng)。分-合式微通道的流體在分支內(nèi)側(cè)產(chǎn)生了漩渦，并且隨著分支傾斜角度的增加，漩渦的大小明顯增加，流體的擾動(dòng)程度不斷增加，這主要是因?yàn)榱黧w流動(dòng)存在慣性，當(dāng)流經(jīng)流道彎曲段時(shí)流體不能突然變向，導(dǎo)致流體在轉(zhuǎn)彎內(nèi)壁面開(kāi)始產(chǎn)生分離現(xiàn)象，彎度越大，分離現(xiàn)象越明顯，甚至出現(xiàn)與主流區(qū)流速相反的回流區(qū)域，即漩渦。當(dāng)分支傾斜角度為90°時(shí)，微通道分支處的漩渦最大。

2.2 壓力特性

不同Re下，微通道總壓降隨分支傾斜角度的變化曲線(xiàn)如圖6所示。從圖6可以看出，隨著Re的增加，不同結(jié)構(gòu)分-合式微通道的總壓降均在不斷增加，且增加幅度也隨著分支傾斜角度的增加而不斷變大。在不同Re下，分-合式微通道的總壓降均高于矩形微通道的總壓降，這主要是因?yàn)樵诜种幃a(chǎn)生了漩渦，引起流體的碰撞與摩擦，使分支處流體的局部阻力損失增加，從而增加了流體流動(dòng)的能量損失。并且微通道的總壓降也隨著分支傾斜角度的增加而增加，這是由于隨著分支傾斜角度的增加，分支處的漩渦不斷變大，導(dǎo)致局部阻力損失增加。當(dāng)Re為970、分支傾斜角度為90°時(shí)，微通道的總壓降最大。

不同Re下，微通道摩擦阻力系數(shù)隨分支傾斜角度的變化曲線(xiàn)如圖7所示。從圖7可以看出，在不同Re下，分-合式微通道的摩擦阻力系數(shù)都大于矩形微通道摩擦阻力系數(shù)。隨著分支傾斜角度的增加，微通道的摩擦阻力系數(shù)也在增加，這主要是因?yàn)樵谖⑼ǖ纼?nèi)分支的存在產(chǎn)生了新的壁面，增加了微通道的流通面積，并且分支處產(chǎn)生漩渦，能量損失增加，使微通道的總壓降增加，摩擦阻力系數(shù)變大，因此隨著分支傾斜角度的增加，摩擦阻力系數(shù)出現(xiàn)逐漸上升的趨勢(shì)。

2.3 溫度特性

當(dāng)Re=322時(shí)，微通道x-y截面的溫度云圖如圖8所示（z=0.25 mm，0 mm≤x≤2.5 mm，7.5 mm≤y≤12.5 mm）。從圖8可以看出，微通道內(nèi)流體沿流動(dòng)方向，其溫度分布都呈現(xiàn)出中心處溫度低、壁面溫度高的趨勢(shì)。矩形微通道內(nèi)由于流速分布一致，沿著流動(dòng)方向，微通道內(nèi)的溫度在不斷上升。分-合式微通道內(nèi)，分支傾斜角度的變化對(duì)微通道內(nèi)流體溫度分布有較大的影響，當(dāng)分支傾斜角度較小時(shí)，分支處的漩渦較小，流體擾動(dòng)程度小，因此分支處的局部換熱系數(shù)較小，分支處的流體溫度較低；當(dāng)分支傾斜角度較大時(shí)，分支處的漩渦變大，流體擾動(dòng)程度增加，局部換熱系數(shù)增加，分支處的流體溫度升高。因此，隨著分支傾斜角度的增加，微通道內(nèi)的溫度分布更加均勻，高溫區(qū)域也在不斷減少。

不同Re下，微通道平均溫度和最高溫度隨分支傾斜角度的變化曲線(xiàn)分別如圖9和圖10所示（橫坐標(biāo)為R處代表矩形微通道散熱器，與下文情況一致）。

從圖9可以看出，隨著Re的增加，各個(gè)結(jié)構(gòu)微通道的平均溫度均不斷降低，但是其變化幅度在不斷減小，這說(shuō)明增加流體的流速可以提高微通道的換熱能力，但當(dāng)流速增加到一定值后，換熱能力的提升效果趨于平緩。在不同Re下，分-合式微通道的平均溫度均低于矩形微通道的平均溫度，這主要是因?yàn)榉种幍牧黧w擾動(dòng)程度大，微通道換熱能力提升。隨著分支傾斜角度的不斷增加，微通道內(nèi)的平均溫度在不斷下降，這是由于隨著分支傾斜角度的增加，分支處的流體擾動(dòng)程度在不斷增加。當(dāng)Re為970、分支傾斜角度為90°時(shí)，微通道的平均溫度最低，為308.9 K，與矩形微通道相比降低了11.9 K。

從圖10可以看出，在不同Re下，分-合式微通道的最高溫度均低于矩形微通道的最高溫度，且隨著分支傾斜角度的增加，微通道內(nèi)的最高溫度在不斷下降。當(dāng)Re為970、分支傾斜角度為90°時(shí)，微通道的最高溫度最低，為316.6 K，與矩形微通道相比降低了14.2 K。

不同Re下，微通道Nu隨分支傾斜角度的變化曲線(xiàn)如圖11所示。從圖11可以看出，隨著Re的增加，各個(gè)結(jié)構(gòu)微通道的Nu均不斷增加，但其增加幅度卻越來(lái)越小。在不同Re下，分-合式微通道的Nu均高于矩形微通道的Nu，并且Nu也隨著分支傾斜角度的增加而不斷增加。當(dāng)Re為970、分支傾斜角度為90°時(shí)，Nu取得最大值，表示其具有最佳的傳熱性能，與矩形微通道相比，在該條件下，Nu最大，為35.8，增加幅度最大，為85.7%。

不同Re下，微通道PEC隨分支傾斜角度的變化曲線(xiàn)如圖12所示。從圖12可以看出，在不同Re下，分-合式微通道的PEC都隨著分支傾斜角度的增加而增加，當(dāng)分支傾斜角度小于45°時(shí)，增加幅度較大，大于45°時(shí)增加幅度較小。分支傾斜角度為45°，Re為646時(shí)，數(shù)據(jù)變化趨勢(shì)明顯不同，分析認(rèn)為，當(dāng)分支傾斜角度為45°時(shí)，Re由484增加到646與Re由646依次增加到970時(shí)相比，摩擦阻力系數(shù)的下降趨勢(shì)明顯要大（見(jiàn)圖7），而Nu下降趨勢(shì)相差不大（見(jiàn)圖11）。PEC與摩擦阻力系數(shù)成反比，所以當(dāng)分支傾斜角度為45°，Re為646時(shí)，PEC增加幅度要大。在本文模擬計(jì)算范圍內(nèi)，當(dāng)Re為970、分支傾斜角度為90°時(shí)，PEC取得最大值，為1.44，表示其具有最佳的綜合傳熱性能。

3 結(jié) 論

通過(guò)對(duì)分-合式微通道散熱器的散熱過(guò)程進(jìn)行模擬計(jì)算，研究了分支傾斜角度變化對(duì)微通道內(nèi)流體流動(dòng)和傳熱性能的影響，主要結(jié)論如下。

1）隨著Re的增加，矩形微通道與分-合式微通道的總壓降不斷增加，但微通道的換熱能力都得到提升，且提升的幅度越來(lái)越小。

2）分支結(jié)構(gòu)打破了微通道的流動(dòng)邊界層，并且流體在分支處產(chǎn)生分離現(xiàn)象，形成漩渦，增加流體擾動(dòng)程度，冷熱流體混合更加充分，高溫區(qū)域不斷減少，微通道內(nèi)的溫度分布更加均勻。

3）隨著分支傾斜角度的增加，微通道的流速不斷增加，分支處的流體漩渦不斷增大，微通道的總壓降和摩擦阻力系數(shù)都呈現(xiàn)不斷增加的趨勢(shì)；最高溫度也在不斷下降，微通道的Nu和PEC也在不斷增加。當(dāng)分-合式微通道的分支傾斜角度為90°時(shí)，微通道的平均溫度為308.9 K，降低了11.9 K；最高溫度為316.6 K，降低了14.2 K；Nu為35.8，增加了85.7%；PEC最佳，為1.44。

研究發(fā)現(xiàn)，分支的引入可以明顯提高微通道散熱器的傳熱性能，并且分支處流體擾動(dòng)程度增加，改善了微通道的溫度分布，但是沿流動(dòng)方向，流體的溫度逐漸升高，換熱能力減少，微通道后端溫度比前端溫度高出很多。為解決這一問(wèn)題，后續(xù)研究可從改進(jìn)微通道內(nèi)分支的排布方式出發(fā)，設(shè)計(jì)更合理的結(jié)構(gòu)，以改善微通道散熱器溫度分布不均的問(wèn)題，為微通道的發(fā)展提供新的思路和方法。

參考文獻(xiàn)/References：

［1］ 寧璐璐，劉清江，楊鳳，等.電子器件冷卻技術(shù)研究進(jìn)展［J］.制冷與空調(diào)，2021，21（12）：1-7.

NING Lulu，LIU Qingjiang，YANG Feng，et al.Research progress of electronic device cooling technology［J］.Refrigeration and Air-conditioning，2021，21（12）：1-7.

［2］ 孫亞松，劉紅敏.微通道強(qiáng)化換熱研究進(jìn)展［J］.應(yīng)用化工，2022，51（4）：1116-1118.

SUN Yasong，LIU Hongmin.Advances in microchannel enhanced heat transfer［J］.Applied Chemical Industry，2022，51（4）：1116-1118.

［3］ PATIL N G， HOTTA T K. A review on cooling of discrete heated modules using liquid jet impingement［J］.Frontiers in Heat and Mass Transfer， 2018，11.DOI： 10.5098/hmt.11.16.

［4］ 趙文博，應(yīng)展烽，宋新陽(yáng)，等.基于熱約束的功率器件散熱器成本最小化方法［J］.科學(xué)技術(shù)與工程，2021，21（30）：12946-12951.

ZHAO Wenbo，YING Zhanfeng，SONG Xinyang，et al.Cost-minimization method for heatsink of power devices based on thermal constraints［J］.Science Technology and Engineering，2021，21（30）：12946-12951.

［5］ DATTA A，SHARMA V，SANYAL D，et al.A conjugate heat transfer analusis of performance for rectangular microchannel with traprzoidal cavities and ribs［J］.International Journal of Thermal Sciences，2019，138：425-446.

［6］ 肖潤(rùn)鋒，侯予，呂坤鵬，等.大深寬比微通道熱沉流動(dòng)換熱特性數(shù)值模擬及優(yōu)化［J］.低溫工程，2018（6）：1-7.

XIAO Runfeng，HOU Yu，LYU Kunpeng，et al.Numerical simulation and optimization of flow and heat transfer characteristics of large-aspect ratio microchannel heat sink［J］.Cryogenics，2018（6）：1-7.

［7］ GHAHREMANNEZHAD A，XU Huijin，ALHUYI NAZARI M，et al.Effect of porous substrates on thermohydraulic performance enhancement of double layer microchannel heat sinks［J］.International Journal of Heat and Mass Transfer，2019，131：52-63.

［8］ 田玉思，焦永剛，孫會(huì)凱，等.通斷型微通道的結(jié)構(gòu)優(yōu)化與流動(dòng)特性研究［J］.機(jī)械工程學(xué)報(bào)，2023，59（4）：274-282.

TIAN Yusi，JIAO Yonggang，SUN Huikai，et al.Study on structural optimization and flow characteristics of on-off microchannels［J］.Journal of Mechanical Engineering，2023，59（4）：274-282.

［9］ BRJAN A，ERRERA M R.Deterministic tree networks for fluid flow：Geometry for minimal flow resistance between a volume and one point［J］.Frantals-Complex Geometry Pattern and Scaling in Nature and Society，2011，5（4）：685-695.

［10］MIGUEI A F.Toward an optimal design principle in symmetric and asymmetric tree flow networks［J］.Journal of Theoretical Biology，2016，389：101-109.

［11］DENG D X， ZENG L， SUN W. A review on flow boiling enhancement and fabrication of enhanced microchannels of microchannel heat sinks［J］.International Journal of Heat and Mass Transfer， 2021， 175.DOI： 10.1016/j.ijheatmasstransfer.2021.121332.

［12］TUCKERMAN D B，PEASE R F W.High-performance heat sinking for VLSI［J］.IEEE Electron Device Letters，1981，2（5）：126-129.

［13］陳孝根，王學(xué)生，陳琴珠，等.微通道冷板換熱及流阻特性計(jì)算與結(jié)構(gòu)優(yōu)化［J］.低溫與超導(dǎo)，2020，48（7）：71-75.

CHEN Xiaogen，WANG Xuesheng，CHEN Qinzhu，et al.Microchannel cold plate heat transfer and flow resistance characteristics calculation and structure optimization［J］.Cryogenics and Superconductivity，2020，48（7）：71-75.

［14］GUNNASEGARAN P，MOHAMMED H A，SHUAIB N H，et al.The effect of geometrical parameters on heat transfer characteristics of microchannels heat sink with different shapes［J］.International Communications in Heat and Mass Transfer，2010，37（8）：1078-1086.

［15］朱崎峰，常昆鵬，薛東東，等.凹槽型微通道內(nèi)液體流動(dòng)換熱特性的實(shí)驗(yàn)研究［J］.河南理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2021，40（3）：100-105.

ZHU Qifeng，CHANG Kunpeng，XUE Dongdong，et al.Experimental study on water flow and heat transfer characteristics in microchannel with grooves in sidewalls［J］.Journal of Henan Polytechnic University（Natural Science），2021，40（3）：100-105.

［16］CHAI Lei，XIA Guodong，WANG Huasheng.Parametric study on thermal and hydraulic characteristics of laminar flow in microchannel heat sink with fan-shaped ribs on sidewalls-part 1： Heat transfer［J］.International Journal of Heat and Mass Transfer，2016，97：1069-1080.

［17］PENG Yi，ZHU Xiugen，CAO Bin，et al.Heat transfer and permeability of the tree-like branching networks［J］.International Journal of Heat and Mass Transfer，2019，129：801-811.

［18］HUANG Pingnan，DONG Guanping，ZHONG Xineng，et al.Numerical investigation of the fluid flow and heat transfer characteristics of tree-shaped microchannel heat sink with variable cross-section［J］.Chemical Engineering and Processing：Process Intensification，2020，147.DOI：10.1016/j.cep.2019.107769.

［19］WANG Gang，GU Yu，ZHAO Luhaibo，et al.Experimental and numerical investigation of fractal-tree-like heat exchanger manufactured by 3D printing［J］.Chemical Engineering Science，2019，195：250-261.

［20］PENCE D.Reduced pumping power and wall temperature in microchannel heat sinks with fractal-like branching channel network［J］.Microscale Thermpphys Engineering，2002，6（4）：319-330.

［21］蔣靜智，劉華峰，安京州.一種分-合式微通道換熱器系統(tǒng)［P］.中國(guó)專(zhuān)利：202123087759.9，2022-05-10.

［22］MURRAY C D.The physiological principle of minimum work（Ⅰ）：The vascular system and the cost of blood volume［J］.Proceedings of the National Academy of Sciences of the United States of America，1926，12（3）：207-214.

［23］楊世銘，陶文銓.傳熱學(xué)［M］.北京：高等教育出版社，2006.

［24］STEINKE M E，KANDLIKAR S G.Single-phase liquid friction factors in microchannels［J］.International Journal of Thermal Sciences，2006，45（11）：1073-1083.