陳鈺龍 蘇開華 吳磊

摘要：矩形微通道已廣泛用于均熱板中，并且已有大量研究對(duì)微通道進(jìn)行優(yōu)化。為了對(duì)矩形微通道的尺寸進(jìn)行優(yōu)化，選擇微通道的溫度分布標(biāo)準(zhǔn)偏差作為標(biāo)準(zhǔn)。首先，根據(jù) Navier –Stokes方程推導(dǎo)了矩形微通道的平均流速公式，然后，又推導(dǎo)了微通道的壁面溫度分布公式，最后采用多項(xiàng)式擬合的方式對(duì)平均流速公式進(jìn)行簡(jiǎn)化，得出微通道的溫度分布標(biāo)準(zhǔn)偏差。數(shù)值研究表明，寬長比越大，溫度分布越好。僅當(dāng)寬度與深度的比率滿足接觸角的函數(shù)時(shí)，可以獲得在微通道的縱向上的最佳溫度分布。為均熱板矩形微通道的尺寸優(yōu)化提供了一定的理論計(jì)算依據(jù)，提出了以微通道的溫度分布標(biāo)準(zhǔn)偏差最小作為優(yōu)化的方向。

關(guān)鍵詞：矩形微通道;平均流速;溫度分布

中圖分類號(hào)：TK124文獻(xiàn)標(biāo)志碼：A文章編號(hào)：1009-9492（2021）11-0109-04

Optimization of Rectangular Microchannel of Vapor Chamber Based on Optimal Temperature Distribution

Chen Yulong1，Su Kaihua2， Wu Lei2

（1. Zhongshan Sansheng IOT Technology Co.， Ltd.， Zhongshan， Guangdong 528437， China;2. Department of Equipment Manufacturing， Zhongshan Torch Polytechnic， Zhongshan， Guangdong 528437， China）

Abstract： The rectangular microchannel has been extensively used in flat heat pipe and lots of study have been conducted to optimize thechannel on the maximum thermal conductivity in a microchannel. To optimize the size of the rectangular microchannel， the standard deviationof the temperature distribution of the microchannel is selected as the standard. First， the mean flow velocity was established by solving theNavier –Stokes equation. Then， the wall temperature distribution formula of the microchannel was derived. Finally， a polynomial fitting methodwas used to simplify the average flow rate formula， and obtain the standard deviation of the temperature distribution of the microchannel. Thenumerical investigation shows that the larger the ratio of width to length， the better the temperature distribution. Only when the ratio of width todepth satisfies a function of the contact angle， the best temperature distribution in longitudinal direction of the microchannel can be got. Theresearch provides a certain theoretical basis for the optimization of the size of the rectangular microchannel， and proposes the minimumstandard deviation of the temperature distribution of the microchannel as the optimization direction.

Key words： rectangular microchannel; mean flow velocity; temperature distribution

0 引言

電子元器件在使用過程中不可避免地都會(huì)出現(xiàn)發(fā)熱，因此需要對(duì)其進(jìn)行散熱。而均熱板具有高導(dǎo)熱性的能力，特別適合用于冷卻高功耗電子元件[1]，因此，對(duì)均熱板散熱性能研究具有重要意義。已經(jīng)有大量的文獻(xiàn)對(duì)均熱板中的流動(dòng)和傳熱特性進(jìn)行研究，姜勇[2]對(duì)兩款內(nèi)部毛細(xì)芯結(jié)構(gòu)不同的均熱板進(jìn)行數(shù)值分析，認(rèn)為內(nèi)部具有燒結(jié)銅粉柱體的均熱板熱性能更加優(yōu)良;趙亮[3]測(cè)試了均熱板和具有相同尺寸電子設(shè)備模塊殼體在不同工況下的散熱性能，發(fā)現(xiàn)均熱板的均溫效果優(yōu)于同尺寸常規(guī)鋁材模塊殼體;基于植物葉片蒸騰作用帶走大量熱的特性，劉旺玉[4]提出了新型仿葉脈分形結(jié)構(gòu)用以優(yōu)化均熱板的熱結(jié)構(gòu);曹紅[5]對(duì)功放腔體結(jié)構(gòu)采用了一體化均熱板的形式，解決了毫米波功率放大器芯片的散熱問題，通過熱仿真和熱測(cè)試表明一體化均熱板可有效降低功放芯片的工作溫度。然而，到目前為止，研究文獻(xiàn)幾乎都是關(guān)注均熱板與其他類型散熱器件的對(duì)比研究，而作為均熱板的內(nèi)部微通道的研究卻很少有關(guān)注。本文關(guān)注溫度在微槽道中的均衡性，通過基于溫度分布均勻性這一標(biāo)準(zhǔn)對(duì)矩形微通道的尺寸進(jìn)行優(yōu)化，為均熱板的優(yōu)化提供一定的理論指導(dǎo)。

1 分析模型及矩形微通道平均流速

本文討論的模型是一種毛細(xì)管通道，具有矩形橫截面，并分布在葉脈狀的分形結(jié)構(gòu)中，如圖1所示。假定微通道在縱向方向上長度為 l ，寬度為2W ，深度（高度）為 H （圖1（b））。其他3個(gè)側(cè)面由銅壁制成。

為了使計(jì)算更加方便，開放的微通道等效于圖2中封閉的矩形通道的下半部分，因?yàn)殚_放通道表面中的流體和蒸氣之間的內(nèi)部黏性摩擦近似為0。對(duì)于封閉的微通道，由于其對(duì)稱性，中間層界面的流速均相等。根據(jù)牛頓內(nèi)摩擦定律，中間層流體之間的內(nèi)黏性摩擦也為0。

毛細(xì)矩形槽道內(nèi)的流動(dòng)控制方程可以由 Navi- er-Stokes 方程得到，Navier-Stokes方程在本笛卡爾坐標(biāo)系下的表達(dá)形式如式（1） 所示。其中， u=（u， v， w ）為速度矢量，μ為液體粘度，p 為毛細(xì)壓。

本文中，毛細(xì)矩形槽道為充分發(fā)展好的流動(dòng)滿足以下條件：（1） 流動(dòng)是穩(wěn)定的，屬于定常流動(dòng)，因而各參量不隨時(shí)間變化，即 =0;（ 2）流動(dòng)速度在 y、z 軸方向的速度分量為0，只有 x 方向的速度，即p=0，w=0;（3） 流體為不可壓縮流體，?.u =0; （4） 毛細(xì)壓力在?p ?p對(duì)于毛細(xì)壓力，重力影響很小，忽略重力因素，即g=0。

式（ 1）可簡(jiǎn)化為：

邊界條件如下所示：

式中：b1、b2為速度滑移長度。

在水的自由界面處，剪切應(yīng)力被忽略，因此：

根據(jù)牛頓內(nèi)摩擦定律，得到：

為了計(jì)算公式（2），應(yīng)有：

結(jié)合式（3）、（4）、（6），得：

其中，βn 滿足：

結(jié)合（2）、（3）、（4）、（7） 得：

對(duì)于式（7），滑移長度 b 的數(shù)量級(jí)為10-3 mm[6]， βn 可以表示為：

然后，可以通過以下公式計(jì)算通過開放式微通道的總流量（定義為每單位面積的流量）：

結(jié)合式（8）（9），流量密度可以表達(dá)為：

Wang等[7]根據(jù)能量原理計(jì)算了入口和出口之間的毛細(xì)管壓力差。在當(dāng)前建立的坐標(biāo)系中，毛細(xì)壓力可以定義為：

式中：σ為表面張力;θ為接觸角常數(shù)。而

則開放式微通道中的平均流速可以寫為：

結(jié)合式（10）～ （13），得平均流速表達(dá)式為：

2 微通道的壁溫

由于矩形通道中的溫度分布非常復(fù)雜，因此可以假定通道橫截面中的壁溫是均勻的。 Wang[8]， Lee[9]和Aghanajafi[10]等許多研究通過水力直徑研究了傳熱特性，即圓形微通道等效于矩形微通道。水力直徑可以定義為：

與流場(chǎng)一樣，溫度場(chǎng)也可以分為入口段和充分發(fā)展段。在二維坐標(biāo)系中，如圖3所示，x 表示流速的方向。入口段溫度場(chǎng)的長度可以表示為：

其中 Pe 是 Berkeley數(shù)，可以定義為：

式中：α為熱擴(kuò)散系數(shù)。

結(jié)合式（16）、（17），入口段溫度場(chǎng)的長度可以寫為：

其中α=1.11×10-4 m2/s[11]。H 和 W 的大小為10-1 mm ，uB的大小為10-1 mm/s ，因此 LT 的大小為10-5 mm 級(jí)別。因此，與微通道的長度相比，它非常小，可以忽略不計(jì)。在這項(xiàng)研究中，只考慮了充分發(fā)展的溫度場(chǎng)。

2.1 x 方向流體溫度分布

對(duì)于恒定熱流，充分發(fā)展的溫度場(chǎng)是相似的，滿足：

式中： Tw 為壁溫; TB 為流體的平均溫度。

流體能量方程可以寫成：

對(duì)于充分發(fā)展的溫度場(chǎng)：

用式（22）代替式（21）中的v ，得：

其中導(dǎo)熱系數(shù) k 可以定義為：

由于通道是對(duì)稱的，對(duì)式（23）在 y （0， d/2）內(nèi)進(jìn)行積分，并結(jié)合式（19），可以得到以下方程：

式中：q 為熱通量。

2.2 x 方向的壁面溫度分布

在恒熱流條件下，TB 和 Tw 滿足方程：

其中 h 是傳熱系數(shù)，可以定義為：

式中：Nu 為 Nusselt常數(shù)，假定 Nu 不會(huì)隨x 改變。

結(jié)合式（25）～ （27）， x 方向的壁溫分布可寫為：

2.3 長度方向溫度分布均勻性的標(biāo)準(zhǔn)

長度方向的平均壁面溫度為：

壁面溫度在長度方向的方差為：

其中標(biāo)準(zhǔn)偏差 s 被認(rèn)為是長度方向溫度分布均勻性的標(biāo)準(zhǔn)，可以通過式（28）～ （30）得：

其中 q 為負(fù)，因?yàn)槲⑼ǖ罏樯帷?/p>

3 結(jié)果分析

如前所述，滑移長度b 的數(shù)量級(jí)為10-3 mm ，幾乎可以忽略。因此，不考慮滑移長度時(shí)，式（14）可以簡(jiǎn)化為：

式中：λ=H/W。

定義：

經(jīng)過數(shù)值擬合后， f （λ）為關(guān)于λ的多項(xiàng)式，當(dāng) 0.1<λ<1時(shí)，表達(dá)式如下：

結(jié)合式（32）～ （34），平均流速的簡(jiǎn)化形式可以表示為：

由于標(biāo)準(zhǔn)偏差 s 是微通道尺寸和uB的函數(shù)，而uB是微通道尺寸和接觸角的函數(shù)，因此，可以推斷出標(biāo)準(zhǔn)偏差 s 是微通道尺寸和接觸角的函數(shù)?？梢赃M(jìn)一步推斷出，在一定的微通道尺寸下可以獲得最佳的溫度分布。結(jié)合公式（31）、（35），標(biāo)準(zhǔn)偏差可以表達(dá)為：

為了獲得標(biāo)準(zhǔn)偏差 s 的最小值，假設(shè) l / H 是一個(gè)確定值，并定義：

對(duì)于銅的矩形微通道，有 cosθ=0.978，圖 4所示為此時(shí) g （λ）隨深度與寬度的一半之比λ的變化曲線。由圖可知， g （λ）首先在一定范圍內(nèi)隨λ而減小，然后增大。因此，當(dāng) l/H 為常數(shù)，且滿足0.75

式（38）表明，深度與寬度的一半之比λ對(duì)溫度分布均勻性有很大影響。盡管它是根據(jù)銅通道的實(shí)驗(yàn)結(jié)果得出的，但仍然適用于其他材料。對(duì)于銅矩形微通道， cosθ=0.978，λ=0.6676，即長寬比為0.3338時(shí)，可以獲得最佳的溫度分布。

4 結(jié)束語

用數(shù)值研究的方法來優(yōu)化微通道的長度方向上的溫度分布均勻性。通過穩(wěn)定的不可壓縮的 Navier-Stokes方程，并經(jīng)過多項(xiàng)式函數(shù)的擬合簡(jiǎn)化平均速度模型。在恒定熱通量的條件下獲得了微通道長度向上溫度分布的函數(shù)，得出溫度的標(biāo)準(zhǔn)偏差。經(jīng)過數(shù)值研究，發(fā)現(xiàn)寬度與長度之比越大，溫度分布均勻性就會(huì)越好，深度與寬度的一半之比滿足接觸角的函數(shù)時(shí)，可以得到最佳的溫度分布。

參考文獻(xiàn)：

[1] HONG X， ALI A， BHATIA R. The use of heat pipes in personal computers[C]//Proceedings of the Thermal and Thermomechani- cal Phenomena in Electronic Systems， 1998 ITHERM The Sixth Intersociety Conference， 1998.

[2]姜勇， 李驥.均熱板散熱器的數(shù)值分析與結(jié)構(gòu)優(yōu)化[J].中國科學(xué)院研究生院學(xué)報(bào)， 2012， 29（2）：169-174.

[3]趙亮， 田灃， 楊龍.均熱板散熱性能實(shí)驗(yàn)研究[J].機(jī)械工程師， 2016（2）：23-25.

[4]劉旺玉， 王力， 羅遠(yuǎn)強(qiáng).仿葉脈分形結(jié)構(gòu)在均熱板蒸發(fā)端的實(shí)驗(yàn)研究[J].華南理工大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2017， 45（1）：118-122.

[5]曹紅.一體化均熱板在某毫米波功率放大器熱設(shè)計(jì)中的應(yīng)用[J].機(jī)械與電子， 2015（2）：38-40.

[6] ROY P， ANAND N K， BANERJEE D. Liquid slip and heat trans- fer in rotating rectangular microchannels [J]. International Jour- nal of Heat and Mass Transfer， 2013， 62（7）：184-199.

[7] WANG H J， TSAI H C， CHEN H K， et al. Capillarity of Rectan- gular Micro Grooves and Their Application to Heat Pipes [J]. TamkangJournalofScienceandEngineering， 2005， 8（3）：249-255.

[8] WANG F， ZHANG J， WANG S. Investigation on flow and heat transfer characteristics in rectangular channel with drop-shaped pin fins [J]. Propulsion and Power Research， 2012（1）：64-70.

[9] LEE P， GARIMELLA S V， LIU D. Investigation of heat transfer in rectangular microchannels [J]. International Journal of Heat and Mass Transfer， 2005， 48（9）：1688-1704.

[10] AGHANAJAFI C， VANDADI V， SHAHNAZARI M R. Investi- gation of Convection and Radiation Heat Transfer in Rhombus Microchannels [J]. International Journal of Research and Re- views in Applied Sciences， 2010， 3（2）：167-176.

[11] CASALEGNO V， VAVASSORI P， VALLE M， et al. Measure- mentof thermal propertiesof aceramic/metal joint by laser flash method [J]. Journal of Nuclear Materials， 2010， 407（2）：83-87.

第一作者簡(jiǎn)介：陳鈺龍（1974-），男，陜西扶風(fēng)人，碩士，工程師，研究領(lǐng)域?yàn)槲锫?lián)網(wǎng)技術(shù)。

（編輯：王智圣）