基于擴(kuò)散界面法的熱管毛細(xì)芯氣液界面蒸發(fā)過(guò)程模擬研究

常 建，蔡杰進(jìn)，譚 冰

(華南理工大學(xué) 電力學(xué)院，廣東 廣州 510640)

空間熱管堆采用熱管作為堆芯向外界熱傳導(dǎo)的媒介，由于缺少重力的作用，毛細(xì)芯對(duì)冷卻劑的泵送就更重要[1-2]。為更有效地進(jìn)行傳熱，對(duì)毛細(xì)芯的相變研究越來(lái)越有必要。由于毛細(xì)壓力的限制，在毛細(xì)芯中所能提供的最大供液流量限制了熱管所能承受的最大熱輸入，過(guò)高則會(huì)導(dǎo)致毛細(xì)芯的燒干。為改善熱管性能，毛細(xì)結(jié)構(gòu)需優(yōu)化其泵送能力和滲透率[3-4]。在大多數(shù)情況下，毛細(xì)芯結(jié)構(gòu)被認(rèn)為是一個(gè)連續(xù)的多孔介質(zhì)，但越來(lái)越多的研究發(fā)現(xiàn)薄膜蒸發(fā)和局部孔隙的馬拉格尼對(duì)流等微觀尺度現(xiàn)象的影響在相變過(guò)程中起重要作用。一般來(lái)說(shuō)，薄膜蒸發(fā)發(fā)生在幾個(gè)微米長(zhǎng)的區(qū)域附近的固-液-汽交界處，長(zhǎng)期以來(lái)一直被認(rèn)為是在這樣的系統(tǒng)中傳熱的主要模式。薄膜傳熱的高效率是由大的分離壓力梯度以及極低的熱阻共同影響的[5-6]。三相接觸點(diǎn)附近的強(qiáng)烈蒸發(fā)導(dǎo)致了沿半月板界面方向的溫度梯度，進(jìn)而導(dǎo)致表面張力梯度，從而產(chǎn)生熱毛細(xì)對(duì)流[7-8]。

本文通過(guò)對(duì)單個(gè)毛細(xì)芯孔隙內(nèi)的界面蒸發(fā)過(guò)程進(jìn)行模擬，以對(duì)整個(gè)毛細(xì)芯的相變提供更準(zhǔn)確的模型。首先在相場(chǎng)法的基礎(chǔ)上，建立一個(gè)從毛細(xì)彎液面蒸發(fā)到空氣中的三維模型，考慮界面處的馬拉格尼效應(yīng)、熱浮力效應(yīng)、界面處的蒸發(fā)相變過(guò)程以及蒸汽擴(kuò)散產(chǎn)生的反作用力。然后將模擬結(jié)果與實(shí)驗(yàn)進(jìn)行對(duì)比驗(yàn)證。最后對(duì)界面處的蒸發(fā)進(jìn)行敏感性分析，討論壁面過(guò)熱度、接觸角和管徑對(duì)相變過(guò)程的影響。

1 薄膜蒸發(fā)傳熱傳質(zhì)模型

擴(kuò)散界面法不僅是一個(gè)簡(jiǎn)單的數(shù)學(xué)計(jì)算，而且包含著一定的物理意義。它是通過(guò)相場(chǎng)變量φ來(lái)得到界面層的信息，而不是去直接追蹤兩種流體界面的變化。它把表面張力等效為場(chǎng)變量的梯度與化學(xué)勢(shì)的乘積，并將它作為一個(gè)體積力加入到Navier-Stokes方程中。相場(chǎng)法不但可計(jì)算流體界面的對(duì)流，而且還保證系統(tǒng)總能量合理減少[9]。

擴(kuò)散界面法認(rèn)為界面是一個(gè)薄的、物理性質(zhì)連續(xù)變化的區(qū)域，且兩相的物性參數(shù)也在界面連續(xù)變化。本文用相場(chǎng)變量φ區(qū)分兩相，類似于水平集方法。計(jì)算時(shí)，第1相的值取-1，第2相的值取1，界面處的值從-1到1變化[10]。計(jì)算區(qū)域內(nèi)每相所占的份額V可寫成相場(chǎng)變量的函數(shù)：

(1)

式中，下標(biāo)v、l分別表示氣體和液體。

在模型計(jì)算中涉及的物性參數(shù)，包括導(dǎo)熱系數(shù)、黏度、密度以及定壓熱容等，可定義為：

B=BlVf,l+BvVf,v

(2)

式中，B為相界面對(duì)應(yīng)的物性參數(shù)。

λ為混合能量密度，其與表面張力和擴(kuò)散界面厚度關(guān)系如下所示：

(3)

式中，εpf為界面厚度參數(shù)，取決于界面處網(wǎng)格的大小。為提高界面變化的準(zhǔn)確度，εpf應(yīng)盡可能小，以保證相場(chǎng)變量φ及兩相物理性質(zhì)在相界面的均勻變化，但同時(shí)為保證數(shù)值計(jì)算的收斂性，相界面的厚度也要與在相界面處的網(wǎng)格密度相匹配。通常取界面處最大網(wǎng)格的1/2作為εpf的值，具有較好的模擬結(jié)果。

(4)

式中：ρ為密度；n為界面法向量。對(duì)n向量求偏導(dǎo)，并整理可得：

(5)

(6)

式中：u為速度場(chǎng)；ξ為相界面平滑函數(shù)：

(7)

在核態(tài)沸騰的過(guò)程中，通過(guò)相界面所發(fā)生的傳熱傳質(zhì)過(guò)程是最為重要和特殊的過(guò)程，在相變過(guò)程中，熱流密度和汽化潛熱對(duì)相變過(guò)程的影響十分重要。相變過(guò)程的質(zhì)量轉(zhuǎn)換率可寫為：

(8)

式中：q為熱流密度；ΔHvl為汽化潛熱；Ml為液體摩爾質(zhì)量。在相界面處，考慮熱流僅通過(guò)導(dǎo)熱方式傳導(dǎo)。將熱流密度通過(guò)傅里葉導(dǎo)熱定律展開(kāi)，式(8)變?yōu)椋?/p>

(9)

式中：T為溫度；k為導(dǎo)熱系數(shù)。上述推導(dǎo)中，主要考慮了流體蒸發(fā)造成的界面移動(dòng)，由于有出口存在，蒸汽的產(chǎn)生未對(duì)界面的移動(dòng)產(chǎn)生直接影響，但蒸發(fā)時(shí)相變帶來(lái)的密度變化大，在微小的空間內(nèi)，汽體從出口擴(kuò)散時(shí)對(duì)界面的反作用力不容忽視。本文考慮由于相變導(dǎo)致的氣液移動(dòng)速度，汽體對(duì)界面的反作用力fvl可表示為：

(10)

在確定了相變質(zhì)量通量后，可對(duì)其他控制方程同樣通過(guò)添加源項(xiàng)和界面平滑參數(shù)的方式，修正在相界面處的表達(dá)式。

相場(chǎng)控制方程(C-H方程)：

(11)

(12)

式中：t為時(shí)間；γ為遷移率。由于各方程所引入的參數(shù)較多，有必要對(duì)各參數(shù)進(jìn)行更詳細(xì)的解釋，以幫助理解方程意義。在對(duì)氣泡運(yùn)動(dòng)的研究中，取常溫下水的表面張力σ為0.076 N/m。

(13)

在微尺度傳熱下，這些溫差變化會(huì)大到足以對(duì)流場(chǎng)造成實(shí)質(zhì)性的影響。此外，由于流體會(huì)傳遞熱量，因此溫度場(chǎng)也會(huì)受到流場(chǎng)變化的影響。

動(dòng)量方程(N-S方程)：

(14)

其中：

(15)

式中：Fg為重力項(xiàng)；I為應(yīng)力張量。

能量方程：

(16)

式中，cp為比定壓熱容。

2 計(jì)算模型與網(wǎng)格驗(yàn)證

通道內(nèi)最小的網(wǎng)格單元邊長(zhǎng)為10 μm，其長(zhǎng)度遠(yuǎn)大于滑移長(zhǎng)度，計(jì)算中視為不可滑移邊界條件。首先，壁面處沒(méi)有熱通量，流體在毛細(xì)力的作用下上升，彎液面在重力的作用下重構(gòu)，并達(dá)到穩(wěn)態(tài)，然后添加壁面初始溫度，液體在彎液面通過(guò)蒸發(fā)被轉(zhuǎn)化為蒸汽并通過(guò)管中擴(kuò)散，同時(shí)下方的入口處由于毛細(xì)作用對(duì)流體進(jìn)行補(bǔ)充，彎液面隨著蒸發(fā)產(chǎn)生波動(dòng)。表1列出了計(jì)算中所需的物性參數(shù)。

表1 數(shù)值模擬計(jì)算采用物性參數(shù)Table 1 Physical property parameter used in numerical simulation calculation

圖1為微通道彎液面蒸發(fā)計(jì)算模型及自適應(yīng)網(wǎng)格劃分。數(shù)值模型為采用軸對(duì)稱的簡(jiǎn)化的二維通道(圖1a)，藍(lán)色部分為液體，灰色部分為氣體，初始彎液面構(gòu)造為圓弧，在相界面初始化過(guò)程中，會(huì)根據(jù)重力和接觸角對(duì)界面進(jìn)行重構(gòu)，入口為壓力入口，其值為液面產(chǎn)生的壓力，出口壓力為0。本研究的重點(diǎn)是彎液面處的蒸發(fā)速率，由于彎液面會(huì)隨著蒸發(fā)不斷下移，故本文在模擬該模型時(shí)，力圖在界面梯度較大的區(qū)域中創(chuàng)建密集網(wǎng)格。但有時(shí)很難預(yù)測(cè)梯度驟變出現(xiàn)的位置。還有一些瞬態(tài)情況，驟變的梯度會(huì)移動(dòng)。這可通過(guò)在出現(xiàn)梯度驟變的所有區(qū)域創(chuàng)建細(xì)化網(wǎng)格來(lái)解決，但通常計(jì)算成本很高，較好的解決方法是采用自適應(yīng)網(wǎng)格劃分?；A(chǔ)粗化網(wǎng)格用于在一定時(shí)間間隔內(nèi)推進(jìn)求解，隨后使用求解結(jié)果基于某個(gè)指示函數(shù)來(lái)細(xì)化網(wǎng)格。最后，使用自適應(yīng)網(wǎng)格再次模擬時(shí)間間隔。

本文采用了三角形網(wǎng)格進(jìn)行區(qū)域劃分。由于這是一個(gè)兩相模型，且主要考慮界面處的蒸發(fā)，因此自適應(yīng)基于進(jìn)行相場(chǎng)變量的梯度。圖1d中黑色線為初始界面，紅色線為移動(dòng)后的界面，顯示了自適應(yīng)時(shí)間間隔結(jié)束時(shí)的解。由圖1可看出，隨著界面的移動(dòng)，對(duì)界面區(qū)域，即相場(chǎng)變量的梯度變化大的范圍進(jìn)行了網(wǎng)格細(xì)化，達(dá)到了所需的計(jì)算精度。

本文使用水作為工作液進(jìn)行了網(wǎng)格無(wú)關(guān)性研究，如表2所列，生成了4組網(wǎng)格，其對(duì)計(jì)算結(jié)果的影響列于表3。經(jīng)過(guò)對(duì)比發(fā)現(xiàn)，網(wǎng)格2相比網(wǎng)格1網(wǎng)格數(shù)量減少了1/2，計(jì)算結(jié)果的誤差在1%左右，滿足網(wǎng)格獨(dú)立性的驗(yàn)證，對(duì)下文的計(jì)算中均采用網(wǎng)格3作為網(wǎng)格劃分的密度。

a——微通道彎液面蒸發(fā)計(jì)算模型示意圖；b——微通道兩相及彎液面處的基礎(chǔ)網(wǎng)格； c——根據(jù)相界面的移動(dòng)和梯度變化后經(jīng)過(guò)細(xì)化的界面附近網(wǎng)格；d——根據(jù)加密后的網(wǎng)格計(jì)算得到的移動(dòng)過(guò)后兩相界面圖1 微通道彎液面蒸發(fā)計(jì)算模型及自適應(yīng)網(wǎng)格劃分Fig.1 Microchannel meniscus evaporation calculation model and adaptive meshing

表2 網(wǎng)格獨(dú)立研究中使用的不同網(wǎng)格Table 2 Different grids used in grid independent research

表3 網(wǎng)格獨(dú)立驗(yàn)證結(jié)果Table 3 Grid independent verification result

3 計(jì)算結(jié)果驗(yàn)證與討論

Buffone等[11]通過(guò)粒子成像系統(tǒng)拍攝不同管徑的石英毛細(xì)孔隙內(nèi)的流場(chǎng)和溫度場(chǎng)分布，實(shí)驗(yàn)工質(zhì)為乙醇和丙酮，采用ITO鍍膜進(jìn)行電加熱，得到了在彎液面附近看到環(huán)形渦旋流場(chǎng)。本文采用實(shí)驗(yàn)中的相似的物性和條件，模擬了豎直放置的微通道彎液面蒸發(fā)，結(jié)果如圖2所示。模擬所得彎液面附近的流場(chǎng)與Buffone等實(shí)驗(yàn)結(jié)果進(jìn)行了對(duì)比，驗(yàn)證了數(shù)值模型的正確性。在彎液面附近液面中心的流體沿著相界面向三相接觸線附近流動(dòng)，推動(dòng)壁面附近的流體回流，最終在入口補(bǔ)償流的作用下，發(fā)生偏轉(zhuǎn)，補(bǔ)充位于彎液面中心的流體。這種現(xiàn)象稱之為馬拉格尼流動(dòng)。液-汽界面處的馬拉格尼對(duì)流驅(qū)動(dòng)液體從管中心線的較熱區(qū)域流向壁的較冷區(qū)域，加速了熱對(duì)流，促進(jìn)了蒸發(fā)。入口處的液體流入速度約為1 mm/s，渦旋附近的流速遠(yuǎn)大于入口的流速，因此界面附近的流動(dòng)主要是由溫度梯度導(dǎo)致的馬拉格尼流引起的。剛流入時(shí)，流體流線平行于管軸，直到接近彎液面在渦旋的影響下發(fā)生偏轉(zhuǎn)，最后從彎液面的中心流向接觸線附近。從Buffone等實(shí)驗(yàn)結(jié)果也可清晰看到兩個(gè)渦旋產(chǎn)生在界面附近，同時(shí)補(bǔ)給的流體明顯發(fā)生了偏轉(zhuǎn)，與模擬所得的流場(chǎng)圖幾乎一致。該結(jié)果驗(yàn)證了模型的正確性，表明建立的基于相場(chǎng)法的界面蒸發(fā)模型可用于蒸發(fā)速率的計(jì)算。

圖2 豎直放置的微通道中流場(chǎng)圖Fig.2 Flow field diagram in vertically placed microchannel

不同壁面過(guò)熱度下彎液面蒸發(fā)時(shí)的溫度場(chǎng)如圖3所示?？煽闯龀跏紩r(shí)壁面附近溫度分層分布，因?yàn)檎舭l(fā)的作用，界面處于飽和溫度，相對(duì)于過(guò)熱的液體，界面溫度最低，且從界面向下溫度逐漸增高。隨著液體在馬拉格尼流渦旋作用下，界面處的高溫流體被帶入管道的中心，最終在彎液面附近產(chǎn)生較大的溫度梯度，且越靠近壁面，溫度梯度越大，說(shuō)明壁面附近的蒸發(fā)越強(qiáng)烈[12-14]。對(duì)比圖2中的流場(chǎng)圖可看出，溫度場(chǎng)的分布和流場(chǎng)的分布具有一致性。不同于毛細(xì)孔隙在層流下的壁面分層傳熱，馬拉格尼渦旋的作用下，壁面處的溫度可更快地傳遞到蒸發(fā)界面，增強(qiáng)了換熱。對(duì)于不同的過(guò)熱度，具有相似的溫度場(chǎng)分布，過(guò)熱度的影響將在接下來(lái)的結(jié)果討論中分析。

圖3 不同壁面過(guò)熱度下彎液面蒸發(fā)時(shí)的溫度場(chǎng)Fig.3 Temperature field of meniscus evaporation under different wall superheat degrees

圖4為不同尺寸的微通道水平放置時(shí)在重力作用下渦旋發(fā)生偏移的模擬結(jié)果與Hemanth等的實(shí)驗(yàn)結(jié)果對(duì)比。在模擬中，毛細(xì)孔隙的軸線是水平的，因此重力效應(yīng)會(huì)影響彎液面附近的流動(dòng)方式[15]。圖中觀察到的不對(duì)稱流動(dòng)是浮力和馬拉格尼效應(yīng)的綜合結(jié)果。如圖中實(shí)驗(yàn)結(jié)果所示，圖中每個(gè)部分均疊加有多個(gè)圖像，以提高對(duì)比度。通過(guò)沿垂直中心平面獲得的長(zhǎng)曝光顆粒條紋圖像表征了浮力和熱毛細(xì)作用對(duì)流動(dòng)的相對(duì)影響，渦旋的偏移與模擬結(jié)果一致。圖中的箭頭指向彎液面上停滯點(diǎn)的位置，模擬中該位置處的界面由于補(bǔ)償流的作用溫度相對(duì)最低，熱毛細(xì)作用最小，與兩側(cè)所產(chǎn)生的切應(yīng)力處于平衡狀態(tài)，當(dāng)在界面處發(fā)生蒸發(fā)時(shí)，它使彎液面區(qū)域液體較離開(kāi)彎液面的回流液體相對(duì)溫度低。溫度的這種差異同時(shí)會(huì)產(chǎn)生浮力效應(yīng)，其強(qiáng)度隨管徑的增加而增加。對(duì)于研究的最小直徑為75 μm的毛細(xì)孔隙，浮力的影響可忽略不計(jì)，并且觀察到與豎直的毛細(xì)孔隙相同的對(duì)稱流型(圖4a)。對(duì)于200 μm的毛細(xì)孔隙，渦旋的不對(duì)稱性開(kāi)始產(chǎn)生，發(fā)現(xiàn)逆時(shí)針底部渦旋大于順時(shí)針上部渦旋，這表明存在不對(duì)稱的環(huán)形渦旋，前面提到的停滯點(diǎn)不再位于液面的中心，開(kāi)始向上偏移。浮力效應(yīng)隨著毛細(xì)孔隙直徑的增加而增大，如圖4c、d所示，逆時(shí)針渦旋的大小隨著毛細(xì)孔隙直徑的增加而變大，停滯點(diǎn)沿彎液面向上移動(dòng)，浮力增加?？煽闯觯摿鲃?dòng)被分為兩個(gè)區(qū)域：由表面張力引起運(yùn)動(dòng)的上部區(qū)域和由密度差驅(qū)動(dòng)的下部區(qū)域。同時(shí)在模擬中還可看到，直徑的增大導(dǎo)致重力在界面的影響也越明顯，圖4a、b中，上下兩側(cè)的三相接觸點(diǎn)大致在同一垂直平面上，圖4c、d中，則由于重力的作用，下方界面向前伸展，在加熱時(shí)，這種不平衡增大了下方的蒸發(fā)界面，加劇了渦旋的偏移。因此，浮力驅(qū)動(dòng)渦旋的尺寸隨毛細(xì)直徑的增加而增大。在界面處產(chǎn)生的相對(duì)溫差驅(qū)動(dòng)了浮力引起的流動(dòng)，其大小取決于蒸發(fā)的質(zhì)量通量。

圖5示出了不同直徑毛細(xì)孔隙彎液面處的蒸發(fā)速率。x0為管徑，x為界面與對(duì)稱軸的軸向距離。從圖中可看出，對(duì)稱軸附近的蒸發(fā)速率均較小，且在約60%的長(zhǎng)度區(qū)域內(nèi)，蒸發(fā)速率均變化不大，在靠近壁面的區(qū)域內(nèi)，蒸發(fā)速率快速增長(zhǎng)，在接近壁面處，不同管徑的界面具有相同的蒸發(fā)速率。在對(duì)稱軸附近，直徑較小的彎液面蒸發(fā)速率較高，其值最多僅有壁面附近蒸發(fā)率的1/10，這是由于較小的彎液面從壁面處由渦旋帶來(lái)的熱量更多，其溫度梯度更大，因而小管徑的彎液面局部蒸發(fā)速率更大。從圖5還可看出，蒸發(fā)主要集中在壁面附近，壁面附近的蒸發(fā)速率與管徑尺寸無(wú)關(guān)。

圖5 不同直徑毛細(xì)孔隙界面處的蒸發(fā)速率Fig.5 Evaporation rate at interface of different diameter capillarie interstices

圖6示出管徑為0.4 mm時(shí)不同過(guò)熱度ΔT和接觸角下毛細(xì)孔隙界面處的蒸發(fā)速率。從圖6a可看出，蒸發(fā)速率是過(guò)熱度的函數(shù)，呈正比關(guān)系，過(guò)熱度的增加并未改變彎液面的相對(duì)蒸發(fā)速率，即過(guò)熱度對(duì)渦旋的位置未產(chǎn)生影響，壁面處的蒸發(fā)速率遠(yuǎn)大于彎液面中心停滯點(diǎn)。圖3中的溫度場(chǎng)也同樣證明了這一結(jié)論，壁面處的界面溫度幾乎恒定，且溫度梯度變化大，蒸發(fā)劇烈；中心處由于補(bǔ)償流的作用，溫度梯度小，蒸發(fā)慢。

微通道尺寸：a——75 μm；b——200 μm； c——400 μm；d——762 μm圖4 微通道水平放置時(shí)在重力作用下 渦旋發(fā)生偏移的結(jié)果Fig.4 Vortex shift results of microchannels placed horizontally under action of gravity

壁面附近的蒸發(fā)速度是壁面上液體厚度的函數(shù)，而這個(gè)厚度是由界面的接觸角決定的，為接觸角的函數(shù)。從圖6b可看出壁面處液膜厚度與蒸發(fā)速率的關(guān)系，其值相對(duì)于接觸角的增加，變化不大。在接近壁面處，接觸角減小導(dǎo)致更好的潤(rùn)濕性，液體更有利于向壁面流動(dòng)，補(bǔ)償壁面處的蒸發(fā)，且液膜的厚度更小，蒸發(fā)相對(duì)更劇烈，考慮到軸向液膜的長(zhǎng)度，則接觸角小的液膜蒸發(fā)速率遠(yuǎn)大于大接觸角的情況。本文潤(rùn)濕性在蒸發(fā)過(guò)程中更主要的是提供更薄的液膜，在厚度一致的情況下，蒸發(fā)的速率基本相同，接觸角小意味著有更長(zhǎng)的薄膜蒸發(fā)區(qū)域，因此在整個(gè)蒸發(fā)表面，有更大的蒸發(fā)速率。

圖6 管徑為0.4 mm時(shí)不同過(guò)熱度和接觸角下毛細(xì)孔隙界面處的蒸發(fā)速率Fig.6 Evaporation rates at capillary interstice interface under different superheat degrees and contact angles with pipe diameter of 0.4 mm

4 半月板上蒸汽擴(kuò)散通量的不均勻性

通過(guò)圖6可知界面處某點(diǎn)的蒸發(fā)速度是液面厚度和溫度梯度的函數(shù)，對(duì)于整個(gè)界面的蒸發(fā)速率，需對(duì)界面作積分以體現(xiàn)在軸向上的液面分布。毛細(xì)孔隙界面處蒸發(fā)模型示意圖如圖7所示，d為半徑，r為界面上點(diǎn)到對(duì)稱軸距離，r0為彎液面的等效半徑。這里假設(shè)重力對(duì)液面的影響簡(jiǎn)化為圓弧，設(shè)彎液面最低點(diǎn)的高度為0，則彎液面的某點(diǎn)處的位置表示為：

圖7 毛細(xì)孔隙界面處蒸發(fā)模型示意圖Fig.7 Diagram of evaporation model at interface of capillary interstice

r0=d/cosθ

(17)

(18)

實(shí)際中，熱管毛細(xì)芯的流體均是選擇潤(rùn)濕性較好的流體，以達(dá)到流體的泵送目的，其接觸角大多在60°以內(nèi)，對(duì)圖6b中60°接觸角進(jìn)行膜厚與蒸發(fā)速率的曲線擬合得到：

m″=22.5+(2.1er/0.11+1.3×10-6er/0.023)/cosθ

(19)

式中，m″為蒸發(fā)速率，g/(m2·s)。

5 結(jié)論

本文通過(guò)對(duì)單個(gè)毛細(xì)芯微孔內(nèi)的毛細(xì)蒸發(fā)過(guò)程進(jìn)行建模，考慮了馬拉格尼效應(yīng)引起的界面剪切作用和液體區(qū)域的熱浮力效應(yīng)，研究了微管內(nèi)部的液體在彎液面的傳熱和傳質(zhì)過(guò)程，具體結(jié)論如下。

1) 在彎液面和管壁的三相接觸處較在彎液面中心處的蒸發(fā)更強(qiáng)。從彎液面中心到壁的溫度梯度導(dǎo)致從管的中心朝向壁的馬拉格尼流。

2) 重力的作用導(dǎo)致液體中的渦旋不對(duì)稱，這種不對(duì)稱對(duì)小管徑影響不大，在大于0.4 mm的情況下更顯著。

3) 液面中心處的蒸發(fā)速率均較小，在靠近壁面的區(qū)域內(nèi)，蒸發(fā)速率快速增長(zhǎng)，在接近壁面處，不同管徑的界面具有相同的蒸發(fā)速率。

4) 壁面附近的蒸發(fā)速度是壁面上液體厚度和過(guò)熱度的函數(shù)，而該厚度是由界面的接觸角決定的，為接觸角的函數(shù)。

5) 隨著孔隙在毛細(xì)芯中的液體有更多的比表面積，文中可看出接觸線附近蒸發(fā)速率最大，因此孔隙減小，蒸發(fā)速率增大。但較小的孔隙會(huì)導(dǎo)致較差的滲透率，發(fā)生燒干，需綜合考慮孔隙減小對(duì)蒸發(fā)的影響。