林昌盛，朱張磊，印萬(wàn)忠1,，姜又華

（1.福州大學(xué)，福州 350108；2.廣東以色列理工學(xué)院，廣東 汕頭 515063；3.東北大學(xué)，沈陽(yáng) 110819）

潤(rùn)濕是最基本的物理現(xiàn)象之一[1]，在傳熱傳質(zhì)、微流體及礦物浮選等領(lǐng)域有著十分重要的應(yīng)用價(jià)值[2-7]。液滴在固體表面的釘扎是潤(rùn)濕現(xiàn)象中最重要的特征之一，而釘扎程度一般由接觸角表征[8]。接觸角是氣液界面與固液界面的夾角[9]，包括前進(jìn)接觸角θa與后退接觸角θr，即分別為接觸線向前擴(kuò)張與向后收縮之前的接觸角，兩者之差被稱為滯后接觸角（也即接觸角遲滯現(xiàn)象），是表征液滴釘扎程度的重要參數(shù)[10]。目前，關(guān)于液滴前進(jìn)與后退接觸角的工作大多在二維平面，例如研究基底微結(jié)構(gòu)對(duì)液滴釘扎特征的影響機(jī)理等[11-13]。

然而，許多涉及到液滴釘扎的應(yīng)用都離不開球面，例如浮選中礦物顆粒的表面。目前已有一些針對(duì)液滴在球面（或者曲面）潤(rùn)濕現(xiàn)象的研究[14-18]。Jasper等[14]采用廣義變分法建立預(yù)測(cè)納米液滴在曲面的接觸角模型，并利用壬烷、十二烷、富勒烯等有機(jī)物作為液滴，成功地證實(shí)接觸角模型的適用性。Viswanadam等[16]將接觸角與曲率半徑和液滴體積相結(jié)合，得出相應(yīng)的模型，并利用試驗(yàn)數(shù)據(jù)來(lái)驗(yàn)證該模型。Wu 等[17]建立了預(yù)測(cè)液滴在曲面的宏觀接觸角模型，并考察了曲率半徑以及固體表面的潤(rùn)濕性對(duì)宏觀接觸角的影響程度。然而，研究液滴在球面上前進(jìn)與后退時(shí)的釘扎特征的工作較少。對(duì)于在球面上的液滴，其接觸角較難測(cè)得（因?yàn)檫@需要明確氣液界面切線與固液界面切線的夾角，其中固液界面的切線極難準(zhǔn)確地觀測(cè)，如圖1b 所示），而在三相接觸線處的氣液界面切線與三相接觸線所在平面形成的宏觀角（如圖1a 所示的θH）較易得到相對(duì)準(zhǔn)確的值。更為重要的是，液滴的受力特征，如氣液界面張力的方向，由宏觀角θH決定[18]。因此，應(yīng)利用液滴在球面的宏觀前進(jìn)與后退角來(lái)研究液滴在曲面上的釘扎特征，但目前關(guān)于曲率半徑對(duì)液滴宏觀前進(jìn)角θHa與后退角θHr的影響機(jī)制還未明確。

圖1 液滴在球面的形貌示意圖及液滴在球面上的實(shí)驗(yàn)測(cè)量圖Fig.1 a) Schematic diagram of a droplet on a spherical surface, b)Experimental diagram of a droplet wetting a spherical surface

本工作基于液滴在球面釘扎時(shí)的吉布斯自由能與體積公式，得到吉布斯自由能、液滴體積、球面曲率半徑、接觸角和宏觀角之間的關(guān)系。其次，把介觀接觸線釘扎的能量引入體系吉布斯自由能，并依據(jù)液滴在吉布斯自由能最小時(shí)呈現(xiàn)出宏觀前進(jìn)角θHa與后退角θHr的現(xiàn)象，利用液滴的本征前進(jìn)接觸角θa與后退接觸角θr預(yù)測(cè)液滴在球面的θHa與θHr，并討論曲率半徑對(duì)宏觀角的影響機(jī)制。最后，利用測(cè)得的宏觀角數(shù)據(jù)驗(yàn)證所預(yù)測(cè)的模型。

1 試驗(yàn)

如圖1 所示，試驗(yàn)通過座滴法測(cè)量去離子水液滴在球面上的宏觀前進(jìn)角θHa和后退角θHr。球面是由PET（聚對(duì)苯二甲酸乙二醇酯）膜包覆的玻璃珠（半徑分別為2、4、6、8 mm）。首先，通過微型注射器在球面上生成一個(gè)體積約為3.0 μL（VL）的液滴；然后，利用注射泵以0.04 μL/s 的速度向液滴內(nèi)部添加或吸取去離子水，當(dāng)三相接觸線出現(xiàn)明顯移動(dòng)之前所觀測(cè)到的宏觀角θH即為宏觀前進(jìn)θHa或后退角θHr。本研究所得的數(shù)據(jù)是至少5 次試驗(yàn)的平均值。

2 結(jié)果與分析

2.1 試驗(yàn)結(jié)果

如表1 所示，測(cè)得液滴在二維PET 平面上的本征前進(jìn)角θa與后退接觸角θr分別為(91±2)°與(55±3)°。試驗(yàn)測(cè)得的液滴在不同曲率半徑球面上的宏觀前進(jìn)角θHa和后退角θHr均隨著曲率半徑RS的增大而減小，并分別趨近于θa與θr。

表1 液滴在不同曲率半徑球面的宏觀前進(jìn)角和后退角Tab.1 Values of droplet apparent angles on spherical surfaces with different curvatures

2.2 結(jié)果分析

圖2 基底表面上液滴的前進(jìn)接觸角與后退接觸角示意圖Fig.2 Schematic diagram of the advancing contact angle (a) and the receding contact angle of a droplet on a realistic flat surface (b)

式中：RS、VL和θ（θa或θr）為已知量，G*、α和θH（θHa或θHr）為未知變量。依據(jù)公式（7）和公式（9）可知，存在無(wú)數(shù)組α和θH的組合及相對(duì)應(yīng)的G*值，且其中必定出現(xiàn)使G*值最小的α和θH組合。

圖3 球缺液滴示意圖（其中黑色曲面與紅色曲面之間的體積為VL，黃色平面與紅色曲面之間的球缺體積為V，VL+V 即黑色曲面與黃色平面之間的球缺體積）Fig.3 Schematic diagram of a spherical cap, where VL represents the volume between the black and red curved surfaces, V represents the volume between the yellow plane and the red curved surface, and VL+ V represents the volume between the black curved surface and the yellow plane

如圖4 所示，以滿足公式（7）與公式（9）為前提，在不同曲率半徑條件下，G*值隨著θH?α（即液滴在球面上的接觸角）變化的情況。從圖4 中可知，不論曲率半徑為何值，G*值最小時(shí)所對(duì)應(yīng)的θHa?α（或θHr?α）值與液滴在平面上的本征前進(jìn)接觸角θa（或后退接觸角θr）的大小相等。這表明在考慮介觀接觸線釘扎最大能量的前提下，當(dāng)體系吉布斯自由能最小時(shí)，液滴在球面上的接觸線處于即將前進(jìn)或后退的狀態(tài)，此時(shí)相對(duì)應(yīng)的宏觀角即為宏觀前進(jìn)角θHa或后退角θHr。體系吉布斯自由能最小意味著體系處于最穩(wěn)態(tài)[19]，因此若將接觸線的最大釘扎能量納入體系吉布斯自由能時(shí)，液滴的前進(jìn)或后退接觸角出現(xiàn)于體系最穩(wěn)狀態(tài)下。此外，圖4 還具有更為重要的實(shí)際意義，即基于公式（7）與公式（9）并使G*最小時(shí)，可直接預(yù)測(cè)液滴在球面上的宏觀前進(jìn)角θHa與后退角θHr。

圖4 在不同曲率半徑條件下G*隨θHr?α 變化的情況和G*隨θHa?α 變化的情況Fig.4 Diagrams of the value of G* with respect to θHr?α (a) and the value of G* with respect to θHa?α for surfaces with different curvatures (b)

圖5 顯示了在滿足公式（7）與公式（9）的條件下，當(dāng)體系吉布斯自由能最小時(shí)，不同曲率半徑的球面基底上液滴宏觀前進(jìn)角與后退角的預(yù)測(cè)情況。從圖5 中可知，隨著曲率半徑的增大，G*的最小值減小，且G*值最小時(shí)所對(duì)應(yīng)的θHa或θHr也減小。此外，依據(jù)上述推導(dǎo)的球面體系吉布斯自由能公式，同理可得，液滴在二維平面上，其前進(jìn)或后退接觸線處于最大釘扎狀態(tài)時(shí)，無(wú)量綱化的體系自由能公式為：

其中θH為液滴在平面上的宏觀角，θ為測(cè)得的液滴在基底上的本征前進(jìn)或后退接觸角[8,20]。因此，如圖5 所示，當(dāng)曲率半徑為無(wú)窮大時(shí)，即球面成為二維平面時(shí)，液滴的宏觀前進(jìn)角或后退角的大小等于本征前進(jìn)接觸角θa（91°）或后退接觸角θr（55°）。

圖5 在不同曲率半徑條件下G*隨θHr 和θHa 變化的情況（其中箭頭代表G*的最小值）Fig.5 Diagrams of the value of G* with respect to θHr (a) and the value of G* with respect to θHa (b) for surfaces with different curvatures, where the arrow indicates the minimum of G*

圖6 展示了在不同曲率半徑的球面基底上，預(yù)測(cè)所得宏觀角（虛線）與試驗(yàn)測(cè)得數(shù)據(jù)（已在表1 列出）的對(duì)比?？芍?，試驗(yàn)值與預(yù)測(cè)值（宏觀前進(jìn)角和后退角）隨著曲率半徑變化的趨勢(shì)一致，但宏觀后退角的試驗(yàn)值普遍略小于預(yù)測(cè)值，這可能是因?yàn)樵谠囼?yàn)中注射器的存在，使得液滴的形貌偏離球缺狀?？偠灾?，試驗(yàn)測(cè)得的液滴在球面上的宏觀前進(jìn)角和后退角與預(yù)測(cè)的數(shù)據(jù)大體上相吻合，這表明本研究所建立的宏觀角預(yù)測(cè)模型的可行性。

圖6 在不同曲率半徑的球面基底上預(yù)測(cè)和試驗(yàn)所得的θHa和θHr 的對(duì)比Fig.6 The comparison of predicted and measured θHa and θHr on spherical surfaces with different curvatures

3 結(jié)語(yǔ)

與以往大多在二維基底的研究不同，本工作瞄準(zhǔn)球面上液滴釘扎的問題，建立了液滴在球面上的宏觀前進(jìn)角與后退角預(yù)測(cè)模型，并討論曲率半徑對(duì)體系吉布斯自由能以及宏觀角的影響機(jī)制。由于不存在絕對(duì)光滑的表面，故將介觀接觸線的釘扎能量納入即將移動(dòng)液滴的吉布斯自由能內(nèi)，且接觸線最大釘扎所產(chǎn)生的能量被假想為勢(shì)能。基于液滴在球面上的吉布斯自由能與體積表達(dá)式，即公式（7）和公式（9），得出球面曲率半徑、液滴體積、本征接觸角與宏觀角之間的關(guān)系式。依據(jù)吉布斯自由能最小時(shí)液滴呈現(xiàn)出宏觀前進(jìn)角與后退角的發(fā)現(xiàn)，利用已知的本征前進(jìn)與后退接觸角預(yù)測(cè)液滴在球面上所對(duì)應(yīng)的宏觀前進(jìn)角和后退角，且所預(yù)測(cè)的結(jié)果被試驗(yàn)測(cè)得的數(shù)據(jù)證實(shí)。所得主要結(jié)論如下：

1）隨著球面曲率半徑的增大，體系吉布斯自由能減小，且最后趨近于液滴在二維平面上的自由能。這表明越彎曲的表面，越易于讓液滴處于不穩(wěn)定狀態(tài)。

2）液滴的宏觀前進(jìn)角θHa和后退角θHr隨曲率半徑的增大而減小。這意味著表面曲率半徑的減小，可增大液滴的宏觀角。

3）在一定誤差范圍內(nèi)，當(dāng)表面的本征前進(jìn)與后退接觸角、液滴體積、曲率半徑已知時(shí)，液滴在球面上的宏觀前進(jìn)角與后退角可以被成功地預(yù)測(cè)。

本研究建立了液滴在球面的宏觀前進(jìn)角與后退角預(yù)測(cè)模型，這在一定程度上豐富了潤(rùn)濕理論的內(nèi)容，但其中仍存在一些不足之處，例如試驗(yàn)所采用的液滴僅僅是純?nèi)軇ㄈルx子水），固體曲面僅包覆親水性材料膜（PET 膜），且固體曲面是對(duì)稱的。但在實(shí)際應(yīng)用中，例如浮選領(lǐng)域，液滴一般是混合表面活性劑或是含有微納米顆粒的混合表面活性劑，而固體表面則可能具有疏水性，且存在微納結(jié)構(gòu)或不均勻的化學(xué)基團(tuán)[21-28]。因此，需要進(jìn)一步研究更復(fù)雜的體系，創(chuàng)建出適合實(shí)際應(yīng)用的理論公式或模型，以期更好地為工業(yè)應(yīng)用提供指導(dǎo)。