具有不同孔隙率多孔介質內(nèi)的蒸發(fā)特性

李鴻如，陳巖，張勁草，辛公明

?

具有不同孔隙率多孔介質內(nèi)的蒸發(fā)特性

李鴻如，陳巖，張勁草，辛公明

（山東大學能源與動力工程學院，山東濟南 250061）

制作了具有3種不同孔隙率（0.85、0.75和0.60）的多孔鎳粉樣本，并進行了毛細抽吸實驗和蒸發(fā)冷卻實驗（包括穩(wěn)態(tài)冷卻實驗，初始狀態(tài)為干態(tài)和濕態(tài)的冷卻實驗），以探究多孔介質抽吸和蒸發(fā)過程規(guī)律。毛細抽吸實驗發(fā)現(xiàn)孔隙率大的鎳粉樣本抽吸速率更快，且根據(jù)毛細單管抽吸實驗和毛細抽吸理論分析結果表明：孔隙率大的鎳粉樣本具有更大孔徑是其抽吸速率更快的根本原因。穩(wěn)態(tài)冷卻實驗結果表明孔隙率0.85的鎳粉樣本具有最快的蒸發(fā)速率，0.60鎳粉樣本蒸發(fā)速率最慢，這是受到多孔樣本抽吸速率主導的結果。在初始干態(tài)冷卻實驗中，孔隙率為0.85的鎳粉樣本具有最大的瞬態(tài)冷卻熱通量，但要達到最大冷卻溫度幅度須取得瞬態(tài)冷卻熱流與冷卻時間的平衡，因此孔隙率0.75樣本冷卻溫度幅度最大。在初始濕態(tài)實驗中，孔隙率為0.85的鎳粉樣本具有最高的燒干溫度和最長的燒干時間，表現(xiàn)出最強的抗燒干性能，但是孔隙率0.75樣本劇烈蒸發(fā)過熱度最小且劇烈蒸發(fā)時間最長，最能有效抑制樣本熱端溫度升高。

多孔介質；微通道；孔隙率；蒸發(fā)；燒干；表面張力

引 言

多孔介質蒸發(fā)問題備受關注，許多研究都表明多孔表面蒸發(fā)傳熱系數(shù)高出光滑表面數(shù)倍，具有明顯優(yōu)勢[1-3]。很多學者圍繞多孔介質參數(shù)對其流動及蒸發(fā)影響開展了數(shù)值模擬研究[4]或實驗研究。李金旺等[5]研究了毛細壓力與孔隙率對多孔介質抽吸性能的影響，表明毛細壓力和滲透率共同主導抽吸性能。員盼峰等[6]研究了毛細壓力和滲透率在多孔介質蒸發(fā)中的作用，表明它們共同影響多孔介質的回液能力，對防止多孔介質局部燒干現(xiàn)象具有關鍵作用。張曉杰等[7]研究了多孔介質的初始顆粒粒徑對蒸發(fā)效率的影響，結果表明不同顆粒粒徑所對應的傳熱系數(shù)隨熱通量變化的規(guī)律并不線性單一。施娟[8]的研究表明多孔介質孔隙直徑與數(shù)量對其內(nèi)部流體流動和氣泡生長及逃逸均具有多重復合作用，綜合影響蒸發(fā)效率。曹小林等[9]對雙多孔介質蒸發(fā)實驗的結果表明存在最佳孔徑使其傳熱系數(shù)和臨界熱通量最大。Kiseev等[10]提出了多孔介質最優(yōu)孔徑搭配的計算方法，并通過實驗證明了孔參數(shù)優(yōu)化搭配可使多孔介質相變傳熱能力達到最佳。

多孔介質蒸發(fā)有很多應用領域，環(huán)路熱管毛細芯是其典型應用之一[11-12]。利用相變換熱和多孔芯的毛細作用力驅動工質循環(huán)流動完成熱量傳遞是環(huán)路熱管毛細芯的主要特點，而孔隙率作為毛細芯關鍵參數(shù)影響很大。Wang等[13]通過調節(jié)填粉密度對毛細芯孔隙率及孔徑等進行調控并實驗驗證了孔隙率59%和滲透率1.25×10?13搭配的毛細芯可使其熱導率、熱容和熱擴散率最佳。王野等[14]實驗發(fā)現(xiàn)具有不同孔隙率的3層毛細芯可以有效降低環(huán)路熱管蒸發(fā)器溫度，減少熱漏，性能優(yōu)于單孔隙率毛細芯。Ling等[15]研究發(fā)現(xiàn)使用60%孔隙率纖維銅板在蒸發(fā)器側和70%孔隙率的纖維銅板在補償器側是室溫范圍內(nèi)最優(yōu)孔隙率搭配，可以得到最低的環(huán)路熱管運行溫度和系統(tǒng)熱阻。Franchi等[16]研究指出不同形狀及大小孔的組合并存是保持毛細芯具有大毛細抽吸力且維持較大滲透率的關鍵，對提高環(huán)路熱管整體性能很有必要。Zan等[17]研究了具有不同孔隙率、孔徑和滲透率毛細芯的抽吸性能和傳熱性能，結果表明存在最優(yōu)參數(shù)組合可使環(huán)路熱管在500 W功率下低溫運行且保持較低熱阻值。徐計元等[18]研究通過添加不同造孔劑得到毛細芯，實驗表明孔隙率越大，孔徑越小且分布越集中的毛細芯抽吸性能越優(yōu)異，所以認為對于提高環(huán)路熱管運行性能越有利。

多孔介質蒸發(fā)問題非常復雜，對其有影響的參數(shù)變量較多，難以得到統(tǒng)一結論；且參數(shù)變化研究尚多集中于對多孔介質及其應用元件的整體性能，難以總結獲得單一變量的影響規(guī)律。

本文主要研究具有不同孔隙率的多孔介質蒸發(fā)換熱特征規(guī)律，通過實驗觀察其蒸發(fā)現(xiàn)象，探究單一參數(shù)孔隙率變化時多孔介質蒸發(fā)換熱的特性，以期為優(yōu)化基于多孔介質換熱過程和系統(tǒng)應用提供實驗支持。

1 實驗思路和方法

1.1 實驗思路

制備3種具有不同孔隙率的非燒結鎳粉樣本，分別進行毛細抽吸實驗和冷卻實驗。毛細抽吸實驗的目的是探究兩種孔隙率的鎳粉樣本不同的抽吸性能，冷卻實驗的目的則是探究在局部加熱情況下不同孔隙率樣本利用毛細力驅動工質流動并對熱端的冷卻特性。實驗所用鎳粉為加拿大Inco公司生產(chǎn)的T-255型粒徑2.2～2.8 μm和T-123型粒徑3～7 μm的高純鎳粉。采用非燒結松散粉末有利于精確控制樣本孔隙率，避免了燒結粉末后體積收縮造成的孔隙率不可控誤差?？刂瓶紫堵实姆椒ㄊ窃谙嗤?guī)格的槽道內(nèi)填充不同質量的鎳粉，并通過施加適當壓力得到體積相同的樣本。兩種樣本設計尺寸均為××= 40 mm × 4 mm × 2 mm（圖1）。每次樣本制作完后使用電子天平和毫米刻度尺測其實際質量與外觀尺寸。鎳金屬密度取8.9 g·cm-3。根據(jù)體積和質量計算出樣本孔隙率分別為0.85（以下簡稱樣本A）、0.75（以下簡稱樣本B）和0.60（以下簡稱樣本C）。所制備的樣本為長且薄的多孔樣本，可以在抽吸實驗中清晰觀察到液體浸潤界面的移動，從而直接觀察到樣本的濕潤或干燥狀態(tài)。

1.2 實驗方法

1.2.1 毛細抽吸實驗 實驗工質選用無水乙醇。抽吸實驗中通過升降臺將裝滿乙醇工質的培養(yǎng)皿從下方緩慢靠近垂直放置的待測樣本。工質與樣本接觸時刻即抽吸開始時刻。整個抽吸過程利用高速攝像機拍攝記錄全程。根據(jù)高速攝像記錄數(shù)據(jù)統(tǒng)計獲得抽吸距離與時間的關系。實驗裝置如圖2所示。

1.2.2 粉末冷卻實驗 冷卻實驗分為3種工況：穩(wěn)態(tài)冷卻實驗、粉末初始狀態(tài)為干態(tài)的冷卻實驗和初始狀態(tài)為濕態(tài)的冷卻實驗（以下分別簡稱為干態(tài)冷卻實驗和濕態(tài)冷卻實驗）。穩(wěn)態(tài)實驗目的是測試3種樣本的蒸發(fā)速率，干態(tài)實驗的目的是測試3種樣本在較高溫度下通過蒸發(fā)作用瞬態(tài)降溫的能力，濕態(tài)實驗的目的是測試它們在供液不足、持續(xù)加熱條件下的抗燒干能力。

在3組實驗中，鎳粉樣本全部垂直放置。樣本在底部通過與截面積相同的銅加熱片對接實現(xiàn)加熱。加熱功率通過直流電源和電加熱片實現(xiàn)調控。銅加熱片由上至下依次布置1、2和3共3根T形熱電偶。其中1距離鎳粉樣本2 mm，近似看作樣本熱端溫度。加熱片與鎳粉樣本之間涂抹導熱硅酯減小接觸熱阻。加熱裝置包裹絕熱材料以減少其散熱損失，實驗臺旁邊布置一根熱電偶記錄環(huán)境溫度，外圍則罩有透明有機玻璃防護罩以減少外界影響和維持實驗環(huán)境穩(wěn)定。實驗所測溫度由Fluke采集和記錄，如圖3所示。

穩(wěn)態(tài)冷卻實驗中，加熱功率3.5 W，首先對鎳粉樣本加熱。加熱片1溫度持續(xù)上升，達到穩(wěn)定后暫停加熱，待1降到110℃后重新開始加熱。當1達到110℃時從樣本上方用注射器對其加入無水乙醇。之后乙醇被抽吸到樣本底端并通過蒸發(fā)作用開始對熱端降溫。實驗中保持乙醇工質足量供應，直到1下降達到穩(wěn)定為止。達到穩(wěn)定蒸發(fā)后，通過測量單位時間內(nèi)所蒸發(fā)乙醇體積計算工質穩(wěn)定蒸發(fā)速率

干態(tài)冷卻實驗操作步驟與穩(wěn)態(tài)冷卻實驗基本相同，唯一不同之處在待1達到110℃后只對樣本滴加定容積0.3 ml無水乙醇。之后乙醇被抽吸到樣本底端并通過蒸發(fā)作用開始對熱端進行降溫，實驗可以觀察到被浸潤的樣本顯現(xiàn)出明顯不同于干燥樣本的深灰色。加熱持續(xù)進行直到所有乙醇通過表面被蒸干，即樣本表面再次變成干態(tài)為止，此時觀察到樣本顏色重新變回淺灰色。乙醇被蒸干時1溫度已經(jīng)接近120℃，然后暫停加熱使1降溫至110℃之下，重復此過程3次。

濕態(tài)冷卻實驗中，鎳粉樣本首先自上向下抽吸0.3 ml乙醇工質使之完全浸潤，觀察到浸潤后的樣本表面呈深灰色。然后對樣本施加3.5 W加熱功率，直至樣本表面變回淺灰色，所含乙醇被完全蒸干為止。重新制備相同參數(shù)樣本后重復此過程。

實驗所用高速攝像設置為每秒拍攝24幀圖像，即圖像捕捉最小時間間隔0.04 s，實驗中時間讀取最小間隔為0.5 s，最大相對誤差8%；實驗用電子天平最小稱量分度0.01 g，最小稱量質量3.4 g，最大相對誤差3.4%；T形熱電偶最大測量誤差為±0.5℃，最小測量溫度為20℃，最大相對誤差2.5%；流動距離最大測量誤差為0.25 mm，由于工質流動前5 mm受入口效應影響明顯不予分析，且被舍棄流動部分數(shù)據(jù)不影響實驗結論的正確性。數(shù)據(jù)從5 mm后開始分析，最大相對誤差5%；被舍棄的前5 mm流動所用時間均小于3.5 s，相對于最小流動全程時間114 s，相對誤差為3%。為獲得有效實驗數(shù)據(jù)，每組實驗都進行了獨立重復實驗以保證結果的重復性和準確性。

2 實驗結果與討論

2.1 抽吸實驗結果

圖4為工質被樣本抽吸過程樣本表面被浸潤所產(chǎn)生的液面前鋒。圖5所示為通過高速攝像所觀測記錄和統(tǒng)計得到的兩類樣本的曲線。從圖中可以看出，樣本A、B、C在相同時間內(nèi)抽吸距離依次減小，即抽吸速度依次降低，達到最終抽吸距離的時間也依次變長。

式(2)表明，微通道內(nèi)的毛細流動距離與通道半徑的1/2次方呈正比關系，即通道半徑越大，工質流動速度越快[30]。第4階段毛細流動由Zhmud等[31]給出了近似顯式解，同樣表明在未迫近毛細抽吸平衡長度的區(qū)間內(nèi)（即本實驗范圍內(nèi)），流動距離（速度）與通道半徑是正相關的。

也即以第3、4階段主導的毛細流動在未迫近平衡高度前，流動距離總與通道半徑呈正相關。為驗證這個結論，進行了單管抽吸實驗，采用長90 mm、內(nèi)徑分別為0.1和0.2 mm的石英玻璃管，測試其在豎直放置時的抽吸水工質特性。

圖6給出了內(nèi)徑為0.1和0.2 mm的單管豎直抽吸距離與時間的關系?？梢钥闯鲈诹鲃泳嚯x未迫近抽吸平衡高度之前（約70 mm前），0.2 mm單管的抽吸速度從一開始就高于0.1 mm單管的抽吸速度，而且差距逐漸拉大。

對比圖5、圖6可以看出，大孔隙率的樣本相比于小孔隙率樣本，可對應于大直徑的0.2 mm單管對比于0.1 mm的單管，都具有更快的抽吸速度和更多的抽吸量。這說明大孔隙率樣本內(nèi)部孔隙空間相對更大，比小孔隙率樣本具有更大的平均當量孔徑，因此大孔隙率樣本可以更快抽吸工質，且抽吸量更多。

2.2 穩(wěn)態(tài)冷卻實驗結果

圖7為樣本A在穩(wěn)態(tài)冷卻實驗中的溫度變化趨勢。由圖可知乙醇工質對樣本底端有明顯降溫作用，至降溫結束達到穩(wěn)態(tài)蒸發(fā)后，平均降溫幅度為19.0℃。樣本B，C溫度變化趨勢與樣本A相同，降溫幅度分別為16.8和7.0℃。可以看出3種樣本降溫幅度隨孔隙率降低而降低。經(jīng)式(1)計算，樣本A、B、C的穩(wěn)態(tài)蒸發(fā)熱通量分別為60、56和41 kW?m-2（穩(wěn)態(tài)蒸發(fā)溫度分別為97、101和109℃），呈現(xiàn)出隨孔隙率減小依次遞減的趨勢。這與3種樣本在抽吸實驗中表現(xiàn)出的抽吸速率規(guī)律一致。這表明本實驗中3種多孔樣本穩(wěn)態(tài)蒸發(fā)速率受限于抽吸工質的速率，孔隙率大的多孔材料因具有更快從冷端補充工質的速率而具有更強的蒸發(fā)冷卻能力。

2.3 干態(tài)冷卻實驗結果

圖8為樣本A的干態(tài)冷卻實驗得到1、2和3溫度隨時間的變化曲線。

可以看出加熱開始后1、2和3迅速上升，至1達到約120℃以后暫停加熱使1降溫至110℃后再恢復加熱并對鎳粉樣本滴加乙醇。滴加乙醇時刻即是干態(tài)冷卻實驗起點。由于乙醇工質被抽吸到底部需要一定時間，所以待乙醇被樣本抽吸至底端，1已上升至110℃以上（此時溫度達到最高，記作h，對應時間記為h），然后在乙醇的蒸發(fā)作用下1開始迅速下降。1下降到最低溫度（記作T）以后，溫度又重新上升。待1上升到一定值后觀察到樣本底端開始出現(xiàn)燒干現(xiàn)象，即出現(xiàn)明顯的干態(tài)和濕態(tài)粉末的分界線（記此時的溫度為臨界溫度c，時間是臨界時間c）。然后干濕態(tài)界限逐漸上移，1溫度也不斷升高。直到界限上移到樣本頂端，樣本被完全燒干，認為一次重復實驗完成。暫停加熱待1下降到110℃后再重新開始加熱進行第2次重復實驗。每個樣本如此重復冷卻實驗至少3次。

圖9和圖10為樣本B、C的冷卻實驗曲線，變化趨勢與樣本A相似。圖11為根據(jù)一維傅里葉導熱定律=?d/d計算得到的流過銅加熱片的熱通量。計算中取108.9 W?m-1?K-1，d/d采用1、2和3之間平均溫度梯度。表1給出了樣本A、B、C的重復實驗統(tǒng)計結果，每個樣本的冷卻實驗至少重復3次，表中數(shù)據(jù)為多次實驗結果的平均數(shù)據(jù)。

表1 樣本A,B,C干態(tài)冷卻實驗結果匯總

?HF是銅加熱片上由乙醇工質蒸發(fā)引起的最大熱通量差（開始滴加乙醇工質時的熱通量與最大熱通量之間的差值）。由表1可知熱通量增長幅度隨孔隙率減小而減小，表明大孔隙率樣本由于抽吸工質速率快而具有更強的瞬態(tài)降溫能力，即大孔隙率樣本具有更強的抽吸能力可以將冷量更快輸運到加熱端。但是因為滴加乙醇工質總量一定，樣本A也因此具有最短的降溫時間h-t。表1顯示樣本A雖然瞬態(tài)降溫能力最強，但是由于其過短的降溫時間導致其降溫幅度h?T低于樣本B。樣本C雖然具有最長的降溫時間，但是受限于其過低的降溫能力，因此降溫幅度h?T最低。樣本B在樣本A、C之間取得平衡，因此降溫幅度最大。

表1表明臨界時間c隨孔隙率減小而增加。這是因為大孔隙率樣本在瞬態(tài)冷卻過程中蒸發(fā)了更大量的乙醇工質，且大的孔隙結構有利于蒸氣的逸出到環(huán)境中，因此會更早出現(xiàn)燒干。

綜上所述，多孔材料瞬態(tài)冷卻能力隨孔隙率增加而增大，但相應冷卻時間卻縮短，被燒干時間也會變短。要達到最大的降溫幅度，需取得冷卻能力與冷卻時間的平衡。本實驗中孔隙率0.75樣本冷卻幅度最大，冷卻效果最佳。

2.4 濕態(tài)冷卻實驗結果

圖12為樣本A、B、C在濕態(tài)冷卻實驗中1的變化趨勢。隨著加熱進行，乙醇工質通過樣本外表面持續(xù)發(fā)生自然蒸發(fā)，但是自然蒸發(fā)所帶走熱量小于輸出熱量，所以3種樣本1溫度均持續(xù)上升。樣本B在1為85℃左右出現(xiàn)一段停滯或緩慢上升階段，樣本A在95℃左右出現(xiàn)停滯或緩慢上升階段，樣本C全程未出現(xiàn)此階段。停滯或緩慢上升階段是由于樣本中的乙醇工質超過沸點、達到一定過熱度并發(fā)生了劇烈蒸發(fā)使加熱端達到熱平衡所致。由圖可知樣本B達到劇烈蒸發(fā)的過熱度最小，且劇烈蒸發(fā)維持時間最長。因此整體運行溫度最低。樣本A達到劇烈蒸發(fā)時過熱度較高，且維持劇烈蒸發(fā)時間較短。樣本C未出現(xiàn)明顯劇烈蒸發(fā)現(xiàn)象。之后1又開始快速上升，樣本底端被觀察到出現(xiàn)干態(tài)粉末時溫度和時間分別記為臨界溫度c、臨界時間c。

加熱持續(xù)進行直到樣本被完全變成干態(tài)，實驗完成。由圖可知由于3種樣本蒸發(fā)條件非常接近，所以可以得到非常相近的終止溫度。

3種樣本的濕態(tài)冷卻實驗各進行3次，對每次實驗結果進行統(tǒng)計，見表2。

表2 樣本A,B,C濕態(tài)冷卻實驗結果匯總

從表中可以看出樣本A的c和c都高于樣本B和C，這是因為樣本A具有最快的抽吸速率，可以快速補充液態(tài)乙醇工質到加熱端，更加長久地維持樣本自身濕潤，推遲燒干極限的到來。樣本C因抽吸速率最慢，補充乙醇工質到加熱端的能力不足，因此最容易被燒干，且燒干溫度較低。樣本B因具有最低的劇烈蒸發(fā)過熱度和最長的劇烈蒸發(fā)時間，因此乙醇蒸發(fā)速率最快，蒸發(fā)量最大。又因為相比于樣本A具有較小的抽吸速率，無法充分補充乙醇工質到加熱端，所以燒干現(xiàn)象出現(xiàn)更早。

綜上所述，樣本A具有最高的臨界溫度和臨界時間，表現(xiàn)出最強的抗燒干能力。樣本B雖然較早地出現(xiàn)燒干現(xiàn)象，但是由于其具有最低劇烈蒸發(fā)過熱度與最長劇烈蒸發(fā)時間，在濕態(tài)實驗中整體運行溫度是最低的。

3 結 論

通過抽吸實驗和冷卻實驗研究了具有不同孔隙率的多孔樣本的蒸發(fā)換熱特性?；趯嶒?，得出主要結論如下。

（1）大孔隙率多孔材料抽吸工質速率更快，抽吸量更大。這是由于大孔隙率材料具有大的孔徑，更有利于提高工質抽吸速率。

（2）多孔材料穩(wěn)態(tài)蒸發(fā)速率隨孔隙率增大而增大，這是由多孔材料的抽吸速率所主導的。

（3）多孔材料在干態(tài)冷卻實驗中的瞬態(tài)冷卻能力隨孔隙率增大而增大，但相應冷卻時間卻隨孔隙率增大而減小。要達到最大的冷卻幅度，須取得冷卻能力和冷卻時間的平衡。本實驗表明0.75孔隙率多孔樣本冷卻幅度最大，冷卻效果最佳。

（4）多孔材料在濕態(tài)冷卻實驗中的燒干時間隨孔隙率增大而增大，即大孔隙率多孔樣本具有更強的推遲燒干現(xiàn)象到來的能力。但是本實驗中0.75孔隙率多孔材料具有最低的劇烈蒸發(fā)過熱度和最長的劇烈蒸發(fā)時間，整體運行溫度最低，表現(xiàn)出最好的抑制熱端溫度上升的效果。其相關機理與調控還需更多實驗研究。

符 號 說 明

cp——比定壓熱容，kJ?kg-1?K-1 g——重力加速度，m?s-2 h——流動距離，mm q——熱通量，W?m?2 R——孔半徑或管半徑，mm Tc——臨界溫度，℃ Th, Tl——最高溫度和最低溫度，℃ t——時間，s tc——臨界時間，s th, tl——最高溫度和最低溫度的對應時間點，s γ——汽化潛熱，kJ?kg-1 θ——接觸角，(°) λ——熱導率，W?m-1?K-1 μ——動力黏度，Pa?s ρ——密度，kg?m-3 σ——表面張力系數(shù)，N?m-1 下角標 c——臨界點 h——最高點 l——最低點

References

[1] 扈本學, 徐彥峰, 蘇光輝, 等. 水平窄縫內(nèi)多孔介質傳熱特性研究[J]. 原子能科學技術, 2009, 43(5): 431-434. HU B X, XU Y F, SU G H,. Characteristics of heat transfer for narrow gap on horizontal heated surface with porous media[J]. Atomic Energy Science and Technology, 2009, 43(5): 431-434.

[2] 左少華, 趙曉玥, 王杰陽, 等. 鋁基Al2O3納米多孔表面大容積池沸騰實驗[J]. 化工進展, 2015, 34(5): 1254-1258. ZUO S H, ZHAO X Y, WANG J Y,. Experimental study on pool boiling of aluminum base Al2O3nano-porous surface[J]. Chemical Industry and Engineering Progress, 2015, 34(5): 1254-1258.

[3] 劉阿龍, 徐宏, 王學生, 等. 復合粉末多孔表面管的沸騰傳熱[J]. 化工學報, 2006, 57(4): 726-730. LIU A L, XU H, WANG X S,. Boiling heat transfer on composite powder porous surface tubes[J]. Journal of Chemical Industry and Engineering (China), 2006, 57(4): 726-730.

[4] 馬強, 陳俊, 陳振乾. 分形多孔介質傳熱傳質過程的格子Boltzmann模擬[J]. 化工學報, 2014, 65(s1): 180-187. MA Q, CHEN J, CHEN Z Q. Lattice Boltzmann simulation for heat and mass transfer in fractal porous media[J]. CIESC Journal, 2014, 65(s1): 180-187.

[5] 李金旺, 鄒勇, 程林. 多孔結構毛細抽吸模型及實驗驗證[J]. 中國電機工程學報, 2011, 31(11): 57-61. LI J W, ZOU Y, CHENG L. Development and experimental verification of capillary pumping model for porous structure[J]. Proceedings of the CSEE, 2011, 31(11): 57-61.

[6] 員盼鋒, 史琳, 晉中華, 等. 薄層多孔層內(nèi)池沸騰液體回流機理分析[J]. 發(fā)電設備, 2013, 27(5): 307-310. YUAN P F, SHI L, JIN Z H,. Mechanism analysis of pool-boiling liquid backflow in thin porous coatings[J]. Power Equipment, 2013, 27(5): 307-310.

[7] 張曉杰, 朱彥雷, 白博峰, 等. 多孔球層內(nèi)沸騰現(xiàn)象與傳熱特性研究[C]// 2007多相流學術會議, 2007. ZHANG X J, ZHU Y L, BAI B F,. Boiling observation and heat transfer characteristics analysis in porous spherical layer[C]// 2007 Multiphase Flow Conference, 2007.

[8] 施娟. 泡沫金屬強化沸騰傳熱過程的研究[D]. 南京: 東南大學, 2015. SHI J. Study on boiling heat transfer enhancement with metal foams[D]. Nanjing: Southeast University, 2015.

[9] 曹小林, 鄭平, 晏剛, 等. 燒結雙多孔介質中沸騰換熱的實驗研究[J]. 宇航學報, 2004, 25(6): 690-693. CAO X L, ZHENG P, YAN G,. An experimental study on evaporative heat transfer in bi-dispersed wick structures[J]. Journal of Astronautics, 2004, 25(6): 690-693.

[10] KISEEV V M, VLASSOV V V, MURAOKA I. Optimization of capillary structures for inverted meniscus evaporators of loop heat pipes and heat switches[J]. International Journal of Heat & Mass Transfer, 2010, 53(9/10): 2143-2148.

[11] MAYDANIK Y F. Loop heat pipes[J]. Applied Thermal Engineering, 2005, 25(5/6): 635-657.

[12] LAUNAY S, SARTRE V, BONJOUR J. Parametric analysis of loop heat pipe operation: a literature review[J]. International Journal of Thermal Sciences, 2007, 46(7): 621-636.

[13] WANG D Z, WANG X Y, ZHOU P,. Influence of packing density on performance of porous wick for LHP[J]. Powder Technology, 2014, 258(258): 6-10

[14] 王野, 紀獻兵, 鄭曉歡, 等. 多尺度復合毛細芯環(huán)路熱管的傳熱特性[J]. 化工學報, 2015, 66(6): 2055-2061. WANG Y, JI X B, ZHENG X H,. Heat transfer characteristics of loop heat pipe with modulated composite porous wick[J]. CIESC Journal, 2015, 66(6): 2055-2061.

[15] LING W S, ZHOU W, LIU R L,. Thermal performance of loop heat pipe with porous copper fiber sintered sheet as wick structure[J]. Applied Thermal Engineering, 2016, 108: 251-260.

[16] FRANCHI G, HUANG X. Development of composite wicks for heat pipe performance enhancement[J]. Heat Transfer Engineering, 2008, 29(10): 873-884.

[17] ZAN K J, ZAN C J, CHEN Y M,. Analysis of the parameters of the sintered loop heat pipe[J]. Heat Transfer — Asian Research, 2004, 33(8): 515-526.

[18] 徐計元, 鄒勇, 程林, 等. 燒結鎳毛細芯的孔參數(shù)控制及其對抽吸性能的影響[J]. 化工學報, 2012, 63(2): 463-469. XU J Y, ZOU Y, CHENG L,. Control of pore parameters and influence on capillary pumping performance of sintered nickel wicks[J]. CIESC Journal, 2012, 63(2): 463-469.

[19] WASHBURN E W. The dynamics of capillary flow[J]. Physical Review, 1921, 17(3): 273-283.

[20] CHEBBI R. Dynamics of liquid penetration into capillary tubes[J]. Journal of Colloid & Interface Science, 2007, 315(1): 255-260.

[21] EXTRAND C W. Forces, pressures and energies associated with liquid rising in nonuniform capillary tubes[J]. Journal of Colloid & Interface Science, 2015, 450: 135-140. .

[22] MARTIC G, GENTNER F, SEVENO D,. A molecular dynamics simulation of capillary imbibition[J]. Langmuir, 2002, 18(21): 7971-7976.

[23] FRIES N, DREYER M. The transition from inertial to viscous flow in capillary rise[J]. Journal of Colloid & Interface Science, 2008, 327(1): 125-128.

[24] FRIES N, DREYER M. An analytic solution of capillary rise restrained by gravity[J]. Journal of Colloid & Interface Science, 2008, 320(1): 259-63

[25] QUERE D. Inertial capillarity[J]. Europhysics Letters, 1997, 39(5): 533-538.

[26] ANDRUKH T, MONAENKOVA D, RUBIN B,. Meniscus formation in a capillary and the role of contact line friction[J]. Soft Matter, 2013, 10(4): 609-615.

[27] BOSANQUET C H. LV. On the flow of liquids into capillary tubes[J]. Philosophical Magazine, 1923, 45(267): 525-531.

[28] ICHIKAWA N, SATODA Y. Interface dynamics of capillary flow in a tube under negligible gravity condition[J]. Journal of Colloid & Interface Science, 1994, 162(2): 350-355.

[29] LUCAS R. Rate of capillary ascension of liquids[J]. Kolloid- Zeitschrift, 1918, 23: 15-22

[30] TAS N R, HANEVELD J, JANSEN H V,. Capillary filling speed of water in nanochannels[J]. Applied Physics Letters, 2004, 85(15): 3274-3276.

[31] ZHMUD B V, TIBERG F, HALLSTENSSON K. Dynamics of capillary rise[J]. Journal of Colloid & Interface Science, 2000, 228(2): 263.

Evaporation in porous media with different porosity

LI Hongru, CHEN Yan, ZHANG Jingcao, XIN Gongming

(School of Energy and Power Engineering, Shandong University, Jinan 250061, Shandong, China)

To investigate the evaporation characteristics of porous media with different porosity, 2 pressed slice samples composed of nickel powder are manufactured with porosities 0.85, 0.75 and 0.60 (named sample A,B and C). The slice samples are long and thin that the wetting frontier movement can be clearly visible during suction. The samples were firstly used in suction experiments in which ethanol penetration was recorded by a high-speed camera. The results show that the sample with larger porosity has higher suction rate. It is because that porous media with larger porosity have larger pore sizes, which can promote the suction rate according to theoretical analysis. Next the samples were heated with constant heat load to test their cooling ability and research their heat transfer efficiency by evaporation. Heating experiments were separated into 3 parts: steady state-cooling, dry-cooling and wet-cooling. Steady state-cooling results show that larger porosity sample has higher evaporation rate mainly due to its higher suction ability. Dry-cooling is aimed to compare samples’ transient cooling ability. The results indicate that the highest heat flux caused by evaporation belongs to sample A while its cooling time is the shortest. Although sample C keeps cooling for the longest, its limited evaporation heat flux constrains its cooling degree. Sample B, with both modest evaporation heat flux and cooling time decreases the heater temperature most. In wet-cooling test (on purpose of testing samples’ performance in dry-out condition), sample A keeps itself wet for the longest time and has the highest dry-out temperature, representing the best ability to avoid dry-out phenomenon. Sample B shows relative low superheat degree and the longest time in intensive evaporation, resulting in its overall working temperature being the lowest.

porous media; microchannels; porosity; evaporation; dryout; interfacial tension

10.11949/j.issn.0438-1157.20170134

TQ 028.8

A

0438—1157（2017）09—3380—08

2017-02-13收到初稿，2017-05-09收到修改稿。

辛公明。

李鴻如（1992—），男，碩士研究生。

國家自然科學基金項目（51641607）。

2017-02-13.

Prof. XIN Gongming, xingm@sdu.edu.cn

supported by the National Natural Science Foundation of China (51641607).