張哲，郎元路，陳佳楠，吳巧燕，計宏偉，李星泊，馬妍，陶柳倩，王瑾悅

（天津商業(yè)大學天津市制冷技術重點實驗室，天津 300134）

研究單個固著液滴的表面潤濕行為及原理，不僅對理解液滴內部的熱質傳輸過程至關重要，而且在噴霧冷卻、冷凝換熱、結霜等領域的實驗及工程應用中也很重要。表面潤濕性表現(xiàn)為液相在有機相表面的親水性與疏水性。自然界諸多生物均擁有奇特的疏水性，例如，依附在荷葉上的晶瑩透亮的露珠呈現(xiàn)珍珠狀的球形，它們仿佛有生命力一般在荷盤上滾來滾去，這種現(xiàn)象也被稱為“荷葉效應”（lotus effect）。在工業(yè)領域，由于親水性表面可以有效凝結液滴，致使各種高分子親水性材料被廣泛研制，包括納米復合物（TFN）膜、熱塑性塑料（SFRT）等。有效掌握液滴的潤濕特性對以上工程研究領域具有重要意義。

迄今，已有大量學者利用表面張力梯度、電荷密度梯度、表面粗糙梯度作為驅動力，以驅使液滴定向移動亦或是驅動分子運動，從而分析液滴的相態(tài)變化及水分子跨尺度運動，即液滴的靜態(tài)與動力學研究。潤濕表面的液滴定向移動是表界面物理力學領域的研究重點。Zhang 等通過分子動力學（MD）模擬研究了納米液滴在覆蓋垂直電場的等溫加熱納米結構表面上電潤濕的動力學。Tang等制定了一種結構梯度和化學梯度相結合的研究方法，系統(tǒng)研究了水滴體積和傾斜角度對水滴運輸速度的影響，通過理論模型和機械分析揭示了水滴輸送的過程。部分學者針對液滴撞擊表面后的反彈、振蕩、擴散及收縮過程進行實驗研究。Yuan等研究了不同粗糙梯度下，韋伯數(shù)對液滴反彈及潤濕狀態(tài)的影響。Kong等系統(tǒng)研究了表面潤濕性對過冷液滴碰撞動力學和結冰的影響，建立了預測過冷液滴沖擊不同表面平均凍結擴散率的理論模型。Ahmad等利用產生振動液滴散熱的叉指電極（IDE）結構，研究了不同電壓、頻率和液滴體積下的液滴振蕩動力學和傳熱效率。此外，電荷密度梯度在微觀層面的作用效果顯著，體現(xiàn)為對正離子、負離子及氫鍵的影響。表面張力梯度通常是液滴在表面的靜態(tài)擴散過程的主要驅動力。Chengara等進行了不同表面活性劑的擴散對照實驗，結論有力支持了表面張力梯度驅動擴散的理論。表面潤濕性影響相變過程中的熱質傳遞階段。Kim 等研究了表面潤濕性對沸騰過程從核態(tài)沸騰到臨界熱流的影響，發(fā)現(xiàn)表面潤濕性的變化顯著影響了沸騰表面的液相和氣相行為，如微層和局部干燥點的生長和收縮、潤濕速度。

液滴浸潤過程中幾何參數(shù)（接觸角、尺寸）的研究集中在理論模型與實驗測算方面，通?？紤]重力、流場、磁場調控液滴潤濕過程，而能量勢壘對浸潤性的影響也比較廣泛。溫度梯度的分布同樣影響液滴內部的表面張力梯度及表面能梯度。Bucher等研究了液滴在不同粗糙度形狀表面的接觸問題。Bai 等在不同的表面溫度條件下，研究了不同潤濕性銅表面的水膜爆炸沸騰現(xiàn)象。結果表明，與疏水表面相比，親水表面的水分子具有更高的溫度梯度和初始熱流，更短的沸騰起始時間和更低的Kapitza 阻力。Ye 等研究了在均勻加熱的水平固體襯底上固著液滴的蒸發(fā)動力學，建立了二維液滴高度的演化方程，并分析了蒸發(fā)和液-固、固-氣、液-氣表面張力對溫度的影響。目前，對凍結液滴融化階段的潤濕過程研究較少，液滴融化階段的研究通常設定在恒定底溫下發(fā)生，而工業(yè)應用的液滴熔融發(fā)生在回溫過程中，因此本文將不同凍結條件及不同升溫條件下的液滴融化階段的表面潤濕過程進行對照實驗，觀測溫度對液滴表面潤濕性的影響。

1 實驗

1.1 實驗裝置

本文采用的實驗系統(tǒng)包括觀測裝置和控溫裝置，如圖1 所示。觀測裝置由DSA-30 液滴形狀觀測臺、計算機控制端組成，觀測臺包括液滴觀測腔、升降部件、移動部件高速攝像機及光源。液滴觀測腔材質為透明聚氯乙烯（PVC），可以避免外界環(huán)境因素對液滴相變造成干擾并保證攝像機可以清晰地觀測腔內液滴的變化過程，此外，觀測腔內部放置有水冷頭作為基板的熱源或冷源，水冷頭通過塑膠管與恒溫水浴箱相聯(lián)結。觀測腔放置于升降部件表面，通過移動升降部件、移動部件并在高速攝像機（ALLIED，德國）的監(jiān)測下找尋最佳觀測位置?？販匮b置由恒溫水浴箱（F-12，JULABO，德國）、半導體制冷片及數(shù)據(jù)采集儀（GP-10）組成，其中，恒溫水浴箱的控溫精度為0.1℃，其內部循環(huán)液成分為乙二醇，半導體制冷片則具有70W 的制冷量，制冷功率高達140W。實驗選用的液體成分為去離子水，將其通過微量移液器產生液滴并滴至基板上（圖2）。去離子水具有離子雜質含量、電導率極低的特性，并更接近理想狀態(tài)下純水的熱物性，將去離子水用于實驗中，所獲得的許多參數(shù)會更接近設計或理想數(shù)據(jù)。基板被放置在液滴觀測腔內的水冷頭表面上，水冷頭表面溫度通過恒溫水浴箱進行調控，液滴及水冷頭表面通過極細熱電偶（WRNT-01，T 型，線徑0.08mm）與溫度采集儀（GP-10）聯(lián)結，液滴及水冷頭表面的溫度數(shù)據(jù)由溫度采集儀實時收集，溫度采集儀的通信通道數(shù)為50，功耗較低（≤54V·A），且精度為0.1℃，便于液滴溫度數(shù)據(jù)的精確采集及水冷頭表面的精確控溫。高速攝像機的分辨率為780×580，經由預實驗的結果對比，最終設置了6.5 的變焦倍率，用以達到極佳的實驗拍攝與錄制效果。DSA-30 液滴形狀觀測實驗臺的相關參數(shù)見表1。

圖1 液滴熔融實驗裝置

表1 DSA-30液滴形狀觀測實驗臺組件及相關參數(shù)

圖2 液滴滴定

1.2 實驗流程

采用純鋁片作為液滴基板。為使實驗基板達到較好的傳熱效果，選用的鋁片尺寸為50.0mm×50.0mm×2.0mm。實驗前需要調試觀測腔的位置及角度，然后打開光源并在攝像機錄制畫面中找到合適的觀測位置。將恒溫水浴箱溫度設置為實驗用基底溫度，水浴箱內的載冷劑通過塑膠管向水冷頭內進行充注并形成循環(huán)流動，待水浴箱將載冷劑溫度調控至設定溫度時，水冷頭表面溫度同時達到了水箱設定溫度，即基底溫度。實驗正式開始前需要打開高速攝像機及光源，調試好觀測臺角度，打開恒溫水浴箱，設置基底溫度。恒溫水浴箱通過軟管與水冷頭相連接，水浴箱控制其水槽內的載冷劑溫度，載冷劑流至水冷頭內并改變其表面溫度，再由溫度采集儀監(jiān)測水冷頭表面溫度，以達到溫度控制的效果。溫度采集儀也被稱為無紙記錄儀，可以將介質溫度信號轉化為電動勢信號最終顯示為溫度數(shù)據(jù)。采用的介質為熱電偶，一端與液滴或基底表面相連，另一端則與溫度采集儀相連，其測溫范圍為-200～200℃，精度為0.1℃。

以上為實驗準備階段，正式液滴實驗共分為兩部分：一是不同凍結條件下觀測液滴的凍結及融化過程，共設置了-4℃、-6℃、-8℃、-10℃、-16℃五種基底溫度，將液滴滴至鋁片表面并將鋁片放置在水冷頭表面，同時啟動計時器進行計時，待液滴完全凍結時，將溫度上調至15℃，觀測凍結液滴在基底回溫過程中的融化階段；二是觀測不同升溫條件（基底溫度6℃、11℃、16℃、20℃、24℃）下凍結液滴的融化過程，首先將半導體制冷片溫度設置為定值（-17℃），再將基板放置在其表面上凍結，爾后將凍結液滴移至觀測艙內的水冷頭表面上，放入觀測腔瞬間啟動計時器，開始視頻采集。使用熱電偶測量液滴溫度，將熱電偶探頭伸至液滴中心區(qū)域，該處溫度具有一定的表征作用，測量過程如圖3(a)所示。待固態(tài)冰相完全消失時停止計時，取下液體樣本與試板，將水冷頭表面仔細擦拭、烘干，更換新的試板，重復以上操作。視頻采集由DSA-30 試驗臺配屬的DSA-4 軟件采集。實驗將4 種基板在不同溫度狀態(tài)下分別進行8 組有效實驗。

圖3 液滴溫度與接觸角測量

在圖像處理階段，首先將視頻文件導入Adobe After Effects CS4 軟件，然后將視頻劃分為40 個5s為一段的短視頻，再把每5s 一段的視頻分解為圖像，通過設定1s 為25 幀，使每段視頻至少選定8個有效圖像。對圖像進行標定處理，經由像素比值法計算出液滴的形態(tài)參數(shù)（高度、底面直徑）。然后將圖像導入Image-J 軟件進行接觸角測量，測量標定過程如圖3(b)所示，將每張圖片測量3 組，最后選用平均值。

2 數(shù)據(jù)與分析

2.1 表面潤濕參數(shù)計算

液滴的高度、接觸半徑、接觸角和有效液滴直徑分別用、、及表示，、、由軟件進行圖像像素長度標定得到。使針頭在液滴錄制視頻中顯現(xiàn)，將視頻導入Adobe After Effects CS4軟件后獲得針頭直徑與同幀圖片中液滴的尺寸（接觸直徑、高度、接觸半徑）的像素，并計算二者的像素比，針頭直徑為0.5mm，則液滴尺寸的實際值可表示為針頭直徑與像素比的乘積。圖4顯示了凍結液滴在試驗基板表面融化過程中的形貌特征，示意圖中的液滴形態(tài)一般存在于熔融過程中后階段，持續(xù)時間較長。

圖4 液滴形貌特征

本文所計算的表面潤濕參數(shù)主要為液滴底部潤濕面積與體積，如式(1)、式(2)所示。

式中，為液滴凍結狀態(tài)時與底板的接觸直徑，mm；為液滴高度，mm。

2.2 不確定度測算

為了確保數(shù)據(jù)的精確性，對計算參數(shù)進行了不確定度分析。將測量或計算得到的參數(shù)使用誤差傳遞公式［式(3)］進行計算。

式中，為計算數(shù)據(jù)的總誤差；為計算數(shù)據(jù)與測量數(shù)據(jù)間的函數(shù)關系；x為測量值；δ為該測量值的誤差。實驗中出現(xiàn)的誤差發(fā)生在液滴滴定及圖像處理兩個階段。在液滴滴定過程中，微量移液器量程為2～20μL，其所允許的最大系統(tǒng)誤差小于3%，最大隨機誤差小于2%。

在圖像處理過程中，軟件測量精度為0.3°，為了保證Image-J 軟件的測量液滴接觸角誤差準確性，對每張圖片中的接觸角進行3次測量，待3次測量結果波動小于0.5°時選用其中間值。液滴長度誤差源自圖像處理過程中的像素誤差，經由像素比計算后得出的像素誤差小于0.022mm。潤濕面積與體積誤差可直接利用誤差傳遞公式求得。將測量誤差及表達式代入傳遞公式可得出：值的計算誤差小于3%，值的計算誤差小于5%。因此，本文所選定的表面潤濕參數(shù)值對液滴潤濕特性分析具有較高的參考與分析價值。

3 結果與討論

3.1 液滴相變過程

凍結條件對凍結液滴回溫過程中的融化階段具有顯著影響。將恒定的基底低溫設置為固定冷卻條件，可以觀測液滴的凍結及融化兩個相變過程，從而更加系統(tǒng)地分析，如圖5、圖6 所示。液滴的凝固過程包含冰晶的生長機制，且呈現(xiàn)較強的規(guī)律性。圖5(a)中為液滴的初始過冷狀態(tài)，由于底溫設置較高，使得液滴保持在過冷相的時間較長。當=270s 時，液滴臨界冰核形成，導致宏觀冰晶快速發(fā)生，凍結鋒面首次出現(xiàn)在基部，使得凝固過程得以宏觀觀測。從圖5(b)中可以看出，凍結鋒面由液滴中心延伸至液滴邊界三相界面處，幾乎與底面平行。隨著凍結鋒面向上平移，冰相區(qū)上側液態(tài)層中的冰晶團簇進入粗化階段，晶態(tài)微觀結構的空間分布更加規(guī)則，液滴的透光度隨之降低。圖5(c)～(e)中顯示了鋒面的移動過程，與圖5(a)、(b)相比，凍結鋒面的移動速率減緩。由于液滴與鋁板交界面成為冷源，導致液滴上部與底面基部之間存在溫差，而且冰相區(qū)的驟增使之與上層液態(tài)區(qū)產生了密度差，因此，在溫度梯度與密度梯度的共同作用下，液相熱傳遞過程相對減緩，結晶速率逐漸變小。水的體積膨脹使得冰相體積大于等質量的液態(tài)水體積，導致凍結鋒面不斷上移直至呈尖狀，演變?yōu)閳D5(h)中的類錐形。

圖5 液滴凍結階段的形態(tài)變化圖像（T=-4℃）

圖6 顯示了回溫條件下的凍結液滴的融化過程。實驗發(fā)現(xiàn)，當基底溫度回升至0℃時，凍結液滴便開始融化。此外，與凍結過程不同，融化階段不會釋放大量潛熱。當=413s時，凍結液滴頂部開始融化，液滴表面開始形成液膜并與鋁板表面的微型液珠合并，造成了圖6(b)中表面接觸直徑的擴大?；販刂?426s時，液滴高度發(fā)生驟降，固液相界面出現(xiàn)，接觸直徑進一步擴大，潤濕作用加強。圖6(c)中的液滴內冰相區(qū)透光度較高，說明冰相區(qū)充斥部分液態(tài)水，內部冰晶呈散狀分布，冰相區(qū)實際為固液混合態(tài)，使其后續(xù)的相變過程中出現(xiàn)塌陷現(xiàn)象［圖6(d)］。冰晶的分散分布是由金屬材料高傳導性所致。由于傳熱速率快，導致液膜剛形成時，凍結液滴內部便產生了大量的水，距離基部越遠，熱量傳導變緩，含水量越少，冰晶的密度越高。因此，相界面處的冰晶密度最低，使其受表面張力、對流等因素的影響，呈不規(guī)則變化。當=270s 時，相界面呈弧狀并向一側傾斜，界面最高點處于液滴表面［圖6(e)］。液滴內部液態(tài)層在升溫過程中產生了熱毛細對流，毛細對流貼近液滴表面，令液滴表面附近的冰晶融化最快，也產生了界面最高點。相界面附近的晶粒大小不一，所受浮升力作用不均勻，加之密度驅動力的影響，使得相界面發(fā)生傾斜。熱毛細對流對冰相區(qū)的影響隨著相界面的上移而逐步減小，冰晶主要在浮升力作用下發(fā)生融化，即浮升融化，導致相界面出現(xiàn)明顯的回升過程。當相界面移動至液滴上部時，由于這一區(qū)域冰晶群排布密集，從而使浮升力及密度梯度的驅動性減弱。因此，溫度梯度與晶簇密集度影響該時段的冰晶融化，密集區(qū)的晶簇融化耗時較長，相界面呈現(xiàn)出波浪狀［圖6(f)、(g)］，直至融化結束［圖6(h)］。以上可以看出，當液滴高度發(fā)生驟降后，液滴的接觸直徑增加不明顯，液滴浸潤性減弱。

圖6 液滴融化階段的形態(tài)變化圖像（T=-4℃）

圖7、圖8 展示了基底處于-8℃時的液滴凍結與融化階段。與圖5、圖6相比，液滴的過冷時間、相變時間顯著縮短，而凍結過程、凍結規(guī)律幾近相同。從圖7、圖8 中可以看出，液滴由凍結至融化階段的回溫時間縮短，說明溫差越大，傳熱速率越快。當=244s時，固液相界面與底面平行，頂部冰晶尚未完全融化。相變持續(xù)至258s 時，頂部冰晶塌陷至下側的冰相區(qū)，界面也出現(xiàn)變化。圖8(a)～(d)中的潤濕面積沒有表現(xiàn)明顯的擴散現(xiàn)象，表觀接觸角的變化同樣不明顯，因此，凍結液滴的接觸角大小影響熔融過程中液滴的潤濕過程。圖8(e)、(f)中，冰粒經歷了浮生融化過程，相界面呈弧狀變化，弧線最低點即冰晶群質量最大點。待融化進展至液滴上部時，弧線開始表現(xiàn)出對稱性［圖8(g)］，液滴鋪展過程更加不明顯。在浮生融化階段，重力不再作為液滴內冰晶的主導力，同時對液滴潤濕過程的影響減弱，使得液滴的鋪展過程受限。圖9中展示了液滴相變的兩個階段，經由對比看出，液滴在融化階段的接觸角減小，潤濕性較為凸顯。

圖7 液滴凍結階段的形態(tài)變化圖像（T=-8℃）

圖8 液滴融化階段的形態(tài)變化圖像（T=-8℃）

圖9 液滴相變流程

圖10 所示為=16℃時凍結液滴融化過程。與圖7、圖8相比，圖10中的接觸角較小，液滴高度同樣較低。在凍結液滴融化的初始階段，表面霜層的融化通?？煊谝旱蝺炔康谋诨?，融化形成的液膜在重力的作用下發(fā)生擴散，隨著基部的液態(tài)層形成，液膜開始引導液滴在表面的鋪展運動，如圖10(c)與圖11(c)所示。表面張力在水平方向的分力與垂直方向的分力均阻礙液滴的浸潤過程，金屬與水的黏結力較弱，并未影響液滴的運動。圖10(e)與圖11(e)中的液滴均未發(fā)生明顯的鋪展，表明重力的作用效果減弱，表面張力垂向的分力主導液滴的運動，抑制了液滴的橫向擴散過程。液滴的潤濕過程伴隨著接觸面及高度的變化，液滴的浸潤模式為不完全浸潤，如圖12所示。

圖10 融化階段的形態(tài)變化圖像（T=16℃）

圖11 融化階段的形態(tài)變化圖像（T=20℃）

圖12 液滴潤濕過程中的表面變化

圖13 展示了液滴在融化階段的潤濕過程。圖13(b)～(d)中的液滴發(fā)生了橫向擴散，液滴底面半徑隨之擴大，并伴隨著與周圍液滴的合并過程。此處的摩擦力源于兩相間的黏附力，由于液滴沒有發(fā)生位移，因此表面摩擦力未影響液滴的潤濕過程。在圖13(e)、(f)中，液滴四周呈明顯的環(huán)狀，該環(huán)狀體為液滴的前驅膜，與液滴主體的分界線為液滴足（圖14），液滴的潤濕過程在前驅膜的牽引下進行，前驅膜內的原子團受驅動向四周運動，造成了液滴的擴散，形成了圓環(huán)。此外，圓環(huán)的外邊緣與周側小液珠的合并會擴大前驅膜的面積，同樣促進了液滴的鋪展運動。

圖13 液滴潤濕過程示意圖（俯視）

圖14 潤濕表面微觀示意圖

3.2 表面潤濕參數(shù)

固體表面的潤濕性主要由界面層內部原子團的性質決定，即由固液兩接觸相共同決定。為方便分析融化階段表面潤濕參數(shù)的機理及變化規(guī)律，首先對液滴在固體表面的受力狀況進行解析。接觸角被作為液體浸潤性的表征量，液滴與固體表面的接觸角通常用楊氏方程定義，如式(4)所示，其中，、、分別為固-氣、固-液與液-氣表面張力。

圖15 展示了液滴在兩種接觸角范圍內的潤濕過程中的受力特征，在垂直及朝上的方向，液滴受支撐力、表面張力的作用，在垂直向下的方向，液滴受重力及黏附力的合力作用。此處的表面張力為氣液表面張力，源于液相分子間的內聚力，用以表示在理想情況下單位表面積的自由能。與其他液體相比，液態(tài)水具有更強的內聚力。表面張力大小與指向受表觀接觸角()（圖14）的影響，當小于90°時，表面張力沿拐點切線方向指向液相，當大于90°時，表面張力沿拐點切線方向指向氣相。氣液表面張力與溫度有著密切關系，二者關聯(lián)性蘊含在式(5)中。式中的（N/m）表示參照物標準溫度下的表面張力，為液體臨界溫度，為溫度系數(shù)（1～1.3），、為常數(shù)，與液相本身性質有關。由式(5)可以看出，隨著溫度的增加，表面張力減小，接觸角也隨之減小。溫度與表面張力這一關系在表2中也有所顯現(xiàn)。由圖15 可以看出，微觀前驅膜的徑向移動發(fā)展為液滴的鋪展過程，所以表面張力實際上限制了前驅膜的擴散。

圖15 表面液滴潤濕過程受力分布

由液體和固體之間的分子相互作用產生的黏附力同樣具有限制或促進液滴鋪展過程的作用。固體表面黏附力存在于界面層中，包括黏結力與黏聚力，分別由固液兩相不同分子間的斥力、吸引力構成，黏結力越強，黏聚力越弱，表面潤濕性越強，反之，則表面潤濕性越弱。與黏附力同方向的為重力，在凍結液滴融化階段，物態(tài)轉化使得重力增大，支撐力同時增強，頂部冰晶融化形成的液態(tài)水在重力作用下沿液膜層流向底部液態(tài)層，覆蓋了表面原有的微觀接觸膜，有效促進了潤濕過程，改變了液滴潤濕參數(shù)。同理，黏附力在固體表面性質與液相的共同作用下影響液滴的浸潤性。當液滴處于穩(wěn)態(tài)時，的大小等于重力與黏附力之和。

除了以上幾種常見力，溫度梯度令液滴內產生了熱毛細力。當液滴在融化過程中的內部溫度大于5℃時，液態(tài)區(qū)的不同溫度梯度使得熱毛細對流產生，進而產生了熱毛細力。將熱對流方向作為熱毛細力的作用方向，流動性集中區(qū)域發(fā)生在液滴中心部且向上流動，因此，熱毛細力對液滴的潤濕過程具有明顯的抑制作用。

綜上，凍結液滴相變過程的浸潤性受多種力的影響，令表面潤濕參數(shù)可以較直觀地呈現(xiàn)其規(guī)律性，使得總結并分析潤濕參數(shù)的變化可得出液滴的力學變化。

將液滴在不同凍結條件下的回溫過程與不同升溫條件下的回溫過程進行實驗對比，相變時間與液滴高度的關聯(lián)圖如圖16、圖17 所示。在不同低溫條件下，底板溫度越接近0℃，凍結液滴升溫越快，相變所需時間越短，反之，相變時間增加，時間區(qū)間為50～90s。該實驗結果受初始溫度及水箱升溫速率的影響，初始溫度越低，水箱升至0℃的時間越長，由于升溫速率恒定，導致相變時間表現(xiàn)出一定的規(guī)律。在不同高溫條件下，基底溫度越高，凍結液滴熔化越快，且時間縮短跨度較為明顯，最低反應時間為15s（6℃），最高為104s（24℃）。顯然，基底溫度對這一階段實驗具有主要影響，基底溫度越高，熱傳遞速率越快，凍結液滴升溫越快，冰相達到融點耗費的時間越短，反之則越長。

圖16 不同基底溫度下的相變時間規(guī)律

圖17 液滴高度變化

兩種條件下的凍結液滴高度的形變規(guī)律具有一定的相似性，這種規(guī)律在圖17 中有所展示。水的物態(tài)轉換伴隨著體積的變化，凍結過程中，液滴的體積增大，頂部產生凸起的尖端，融化過程中，液滴體積縮小，高度變化首先由頂端的融化引發(fā)。從圖17 中可以看出，液滴的頂部首先緩慢融化，高度緩慢降低，液膜開始包裹液滴周身，液滴基部液態(tài)層形成后與液膜匯聚，在重力的作用下，液膜底部開始向四周鋪展。當液滴頂部發(fā)生冰晶塌陷時，液滴高度驟降，基部液態(tài)層體積占比較大，重力作用增強促使前驅膜擴散，液滴潤濕過程進一步發(fā)展。待高度驟降后，液滴高度開始緩慢降低，液滴接觸面緩慢擴大，這是由于表面張力的抑制作用逐漸增強，在表面相的分子會在表面張力作用下更有序地排列，而體相的液體分子是無序排列的，潤濕進程的驅動力隨之減弱，直至高度不再改變。

接觸角的變化規(guī)律能直觀地展示液滴的潤濕特征。從圖18中可以看出，接觸角處于63°～88°，且始終呈下降趨勢，接觸角越小，底面潤濕面積越大，表明潤濕性越強。由于接觸角的大小受每次在常溫基板上滴定產生的穩(wěn)態(tài)液滴的靜態(tài)接觸角影響，而每次產生的靜態(tài)接觸角具有較強的隨機性，但一般在10°以內波動，在滴定完成后將承載液滴的常溫基板放置在低溫表面上進行凍結，所以接觸角的大小整體分布與低溫工況（-4℃、-6℃、-8℃、-10℃、-16℃）無關。此處將接觸角下降過程分為兩個階段，在初始階段，接觸角的下降趨勢較強，該階段對應凍結液滴的初始融化過程，凍結液滴在升溫過程中發(fā)生融化并展開液滴鋪展運動，導致初始接觸角變化。液滴頂部輪廓對接觸角測量影響較大，當＜90°時，表面張力的水平方向分力與液滴的鋪展方向相反，表現(xiàn)為抑制性，重力則促進液滴的鋪展過程，頂端冰晶融化形成的液態(tài)水在重力的主導作用下沿壁面擴散，使得液滴鋪展性明顯。而隨著表面張力的水平分力不斷增強，液滴潤濕過程受阻，因此，在第二階段的潤濕過程，接觸角下降趨勢平緩，后期的表面張力為主導力，液滴不再發(fā)生擴散。

圖18 接觸角動態(tài)變化

圖19 中展示了兩種條件下液滴的潤濕面積變化過程?？梢钥闯?，兩種條件下液滴的潤濕性具有明顯的不同。圖19(a)中將數(shù)據(jù)進行了擬合處理，液滴底面積呈遞增趨勢，整體數(shù)值處于2.046～4.47之間，液滴潤濕性較強，而圖19(b)中的面積只在初始時發(fā)生上升，后期幾乎沒有明顯變化。兩部分實驗均伴隨著溫度的升高，表面張力都發(fā)生了減弱，而接觸角始終小于90°，表面張力的作用性相同，因此，表面張力不是造成這一現(xiàn)象的原因。在凍結后的回溫過程中，液滴溫度由零下逐步升至零上，溫差變化具有規(guī)律性，這一點可由圖16 中的時間變化規(guī)律側面證實。由于溫度的升高，固體表面有序單層水的氫鍵結構被破壞，水層內水分子間氫鍵壽命變短，進而導致單層水內的水分子與其上方水滴形成更多氫鍵，從而使得單層水更親水，表面接觸角變小。將低溫液滴移至不同高溫面上進行升溫時，液滴的溫度變化趨勢加快，不同區(qū)域的溫度梯度相差較大，液滴幾乎沒有鋪展運動。表明溫度梯度對液滴潤濕性具有重要影響，甚至抑制液滴的整個融化階段的潤濕過程。液滴體積通常伴隨潤濕過程而逐漸減小，如圖20 所示。液滴體積受接觸面積與高度的影響，在融化初始階段高度主要影響液滴體積變化，使得體積趨勢與高度變化趨勢相似，初始值區(qū)間為1.337～1.815，1.097～1.581。這一特征在不同升溫條件下更加明顯，由于潤濕面積幾乎不變，液滴體積便只與高度相關。從圖20(a)中可以看出，基底溫度為-4℃的條件下，液滴體積降低最快，溫度最快達到穩(wěn)定，這是由于基底溫度較為接近0℃。當基底溫度在-16℃條件下回溫時，液滴的體積在降低后出現(xiàn)微弱升高的現(xiàn)象。溫度越低，基板表面的霜層覆蓋越密集，在回溫過程中，霜晶融化形成的液珠與液滴發(fā)生合并，隨著溫度的升高，融合現(xiàn)象進一步發(fā)生，潤濕面開始擴大，最終導致體積的緩慢增加。升溫條件下，由于溫度的驟然升高，基板表面的霜晶熔化后不會在板面長期留存，不會發(fā)生液滴合并現(xiàn)象，而且隨著基底溫度的升高，液滴三相線處的分子甚至在獲得動能后脫離液滴。

圖19 液滴潤濕面積變化

圖20 液滴體積變化

上文提到熱毛細力同樣影響液滴的潤濕過程。液態(tài)層加熱方向的溫度梯度產生了表面張力梯度，表面張力梯度驅動液體由低表面張力區(qū)域向高表面張力區(qū)域流動，稱為Marangoni 效應（熱毛細效應），因此對液滴溫度梯度與熱毛細效應展開研究。圖21 中展現(xiàn)了不同升溫條件下液滴內部溫度變化規(guī)律。從圖中看出，在融化過程中，液滴內部溫度在0℃時作短暫停留，在40s 時溫度便不再改變，但部分液滴的融化過程依然進行。液滴的熔點為0℃，冰晶溫度達到熔點時開始持續(xù)吸熱，溫度上升標志著物態(tài)轉換的開始，冰相逐漸轉化為液態(tài)水?；诇囟仍礁?，液滴的升溫速率越大，融化進程越快，使得冰相達到完全融化狀態(tài)之前，液滴溫度已經穩(wěn)定。

圖21 液滴內部溫度變化（升溫條件）

在表面張力梯度的驅動作用下，形成了克服重力和液體內部浮力的熱毛細流動，流體流動是由熱毛細作用引起的。對于液體來說，Marangoni效應表現(xiàn)為高表面張力區(qū)域對低表面張力區(qū)域的拉力作用，使得液體流向高表面張力的區(qū)域。此處引入外凸液面的數(shù)計算公式，用以表征液滴內Marangoni流動表面張力與黏性力之比，如式(6)所示。式中的為液滴的氣液表面的溫度梯度，為液滴自由表面的弧長，計算過程如式(7)、式(8)所示。

式中，為水的表面張力（表2）；為液滴的中心溫度，℃；為三相點處的溫度，此處采用底板表面溫度，與中心溫度的差值Δ見表3；為熱導率，W/cm；為水的密度，g/cm；c為比熱容，kJ/(kg·℃)；為黏度，N·s/m；熱物性參數(shù)值在表4中列出。

表2 去離子水的表面張力

表3 液滴溫度差

表4 去離子水的熱物性參數(shù)

當基底溫度處于10℃左右時，液滴中心溫度與底面溫度幾乎相同。隨著基底溫度的升高，液滴中心溫度與底面溫度的差值逐漸增大，|Δ|也隨之增大，見表3。

在液滴熱量傳導過程中，基底面將熱量傳至液滴底部并融化為液態(tài)層，隨著冰相的不斷融化與熱量的持續(xù)傳遞，液態(tài)層開始產生溫差。由于表面張力對溫度的依賴性，產生了表面張力梯度，進而造成了界面切應力的不平衡，流體開始由下向上移動，并沿氣液表面下降，然后在底面高溫處再次上升，并在腔內形成一個逆時針旋轉的渦流，液滴的形貌決定了這種渦流主要沿液滴表面排布。液滴內存在一些具有溫度梯度的微尺度區(qū)域，這些區(qū)域同樣形成了很多的熱毛細小胞流，進而構成某個微區(qū)的循環(huán)流動，影響主渦流的作用效果。圖22 中展示了不同升溫條件的溫度差與數(shù)的關聯(lián)性。由圖中可以看出，基底溫度越高，溫差越大，數(shù)越大，液滴內流體的熱毛細力越強。液滴內部各熱毛細環(huán)狀流在液滴底面邊緣區(qū)域匯聚上升流動，重力會阻礙其流動，熱毛細力越強，重力對液滴運動的影響便越弱，流體流動強度增強。基底溫度越高，熱流密度越大，熱傳導過程得以加強，令Marangoni 對流越強烈，導致液滴在融化階段幾乎沒發(fā)生潤濕運動。

圖22 Ma數(shù)與溫差變化規(guī)律（升溫條件）

4 結論

通過將鋁板表面上不同凍結條件（基底溫度-4℃、-6℃、-8℃、-10℃及-16℃）與不同升溫條件（基底溫度6℃、11℃、16℃、20℃、24℃）下的液滴相變過程與融化階段的潤濕特性進行了可視化實驗，整理了不同條件下液滴相變規(guī)律，并結合表面潤濕參數(shù)對液滴的表面潤濕過程進行了力學分析，研究結果如下。

（1）在融化階段，液滴的潤濕性主要受重力、表面張力、熱毛細力的影響，重力對液滴的橫向擴散促進作用、表面張力的作用與接觸角相關，本文表現(xiàn)為抑制性，熱毛細力受底板溫度影響，同樣具有抑制液滴潤濕過程的作用。

（2）不同凍結條件下，凍結液滴的潤濕過程主要發(fā)生在融化初始階段，重力促進液滴的潤濕過程，液滴接觸角處于63°～88°之間，接觸角越小，潤濕面積越大。在潤濕后階段，接觸角減小，重力的作用減弱，表面張力的作用增強，液滴的擴散進程受阻，體積的降低趨勢也變緩，隨著融化的進行，溫度保持穩(wěn)定，表面張力不變，液滴保持靜止。

（3）不同升溫條件下，凍結液滴的潤濕過程幾乎沒有發(fā)生，熱毛細力與表面張力在潤濕過程中占據(jù)主導性，隨著基底溫度的升高，液滴內部與三相線溫差逐漸增大，數(shù)呈增加的趨勢，數(shù)值由1802 增至22876，熱毛細力始終抑制液滴的運動，使得液滴沒有明顯的移動過程。

（4）凍結條件與升溫條件下的回溫過程對比發(fā)現(xiàn)，凍結液滴高度變化規(guī)律相同，隨著融化的進行，液滴高度驟降，然后緩慢降低。升溫條件下的熱流密度較高，液滴內Marangoni 對流強烈，熱毛細力的抑制性作用增強，液滴浸潤過程不明顯。

通過實驗對比，升溫條件下的實驗相比實踐生產具有一定參考性，相變時間區(qū)域差值相對較小，且可用于液滴內流動傳熱機理的研究。然而，實際生產中的霜層融化階段具有明顯的潤濕特征，凍結條件下的潤濕性分析用液滴可以作為研究單元，升溫條件下的液滴未表現(xiàn)出明顯的潤濕性，不能作為底面加溫熱除冰技術的研究單元。