李文 馬驍婧? 徐進(jìn)良2) 王艷 雷俊鵬
1) (華北電力大學(xué), 低品位能源多相流與傳熱北京市重點(diǎn)實(shí)驗(yàn)室, 北京 102206)
2) (華北電力大學(xué), 電站能量傳遞轉(zhuǎn)化與系統(tǒng)教育部重點(diǎn)實(shí)驗(yàn)室, 北京 102206)
液滴在固體表面的潤濕行為是自然界中普遍存在的物理現(xiàn)象, 廣泛應(yīng)用于印刷工藝、涂料技術(shù)、交通運(yùn)輸、航空航天和能源系統(tǒng)等領(lǐng)域[1].眾所周知, 潤濕狀態(tài)(Cassie, Partial Wenzel和Wenzel)是影響冷凝液滴生長和脫離的關(guān)鍵因素[2,3].受自然界中“荷葉效應(yīng)”[4-6]的啟發(fā), 越來越多的學(xué)者關(guān)注超疏水表面的滴狀凝結(jié)換熱及強(qiáng)化機(jī)理.受液滴潤濕模式的影響, 超疏水納米結(jié)構(gòu)表面冷凝存在彈跳、滾動和滑動三種液滴脫落模式[7], 研究液滴在納米結(jié)構(gòu)表面的潤濕行為對減小液滴脫離尺寸, 促進(jìn)冷凝液滴脫落, 進(jìn)而強(qiáng)化滴狀冷凝傳熱有著重要意義.
目前研究的粗糙表面上的潤濕狀態(tài)主要分為Wenzel模型[8]和Cassie模型[9]兩種理論.Wenzel[8]認(rèn)為粗糙表面的存在使得實(shí)際固-液接觸面積大于所觀察到的表面面積, 并假設(shè)液滴完全充滿納米結(jié)構(gòu)并潤濕壁面.潤濕表面的表觀接觸角(θw)和本征接觸角(θe)之間滿足關(guān)系:
其中r為粗糙度因子, 表示液滴與固體表面接觸的實(shí)際面積與投影面積的比值.Cassie和Baxter[9]進(jìn)一步拓展了Wenzel方程, 認(rèn)為液滴在粗糙表面上的接觸是一種復(fù)合接觸, 當(dāng)表面較疏水時(shí), 液滴不能完全潤濕凹槽, 液滴與表面之間有空氣存在,所以表觀上的固-液接觸面實(shí)際上由固-液接觸面和氣-液接觸面共同組成.復(fù)合表面的表觀接觸角θc可以表示為其中f表示與液體接觸的實(shí)際固體表面占比, 稱為相面積分?jǐn)?shù).
液滴在固體表面的潤濕性主要取決于表面的化學(xué)組成[10-15]、微結(jié)構(gòu)形貌[16-18]和幾何結(jié)構(gòu)尺寸[15,19-21].He等[22]采用實(shí)驗(yàn)方法研究了水滴在不同浸潤性和幾何參數(shù)的表面上的潤濕行為, 研究表明, 當(dāng)水滴直接在固體基板沉積時(shí)形成接觸角較大的Cassie潤濕狀態(tài), 當(dāng)水滴以一定高度降落到表面時(shí)形成接觸角較小的Wenzel潤濕狀態(tài).Bico等[23]討論了液滴在納米結(jié)構(gòu)表面的潤濕行為及不同潤濕狀態(tài)轉(zhuǎn)換的可能性.結(jié)果發(fā)現(xiàn)增大柱高和柱寬、減小柱間距均有利于Cassie液滴的形成.Marmur[21]通過計(jì)算潤濕面積定量評價(jià)液滴在粗糙表面的潤濕行為.研究結(jié)果補(bǔ)充了Wenzel[8]和Cassie[9]理論.由于實(shí)驗(yàn)技術(shù)條件和理論模型的限制, 當(dāng)液滴的特征尺度為微納米級別時(shí), 表面積與體積的比值明顯增大, 宏觀尺度下能夠忽略的某些影響因素(分子間的范德瓦耳斯力、液體流動出現(xiàn)的壁面滑移效應(yīng)、傳熱過程出現(xiàn)的溫度階躍等)變得不容忽視.因此, 在微納尺度下, 物理現(xiàn)象所具有的尺度效應(yīng)和界面效應(yīng)逐漸成為探究潤濕和傳熱等宏觀尺度現(xiàn)象影響機(jī)理的重要因素.近年來,分子動力學(xué)(molecular dynamics, MD)模擬技術(shù)快速發(fā)展, 并被廣泛應(yīng)用于探究納米尺度下分子間相互作用及機(jī)理.Shi和Dhir[24]討論了采用MD模擬得到的水滴在鉑表面潤濕的接觸角和氬液滴在鉑表面潤濕的接觸角, 與實(shí)驗(yàn)中水滴在鋁表面潤濕的接觸角的關(guān)系, 發(fā)現(xiàn)上述三種接觸角隨體系溫度的變化規(guī)律一致, 均隨溫度的增加而降低.Li等[19]研究了氬液滴在固體鉑表面上的潤濕行為, 結(jié)果表明, 液滴的接觸角及總蒸發(fā)量隨表面浸潤性的增加而減小.Yong和Zhang[25]研究了改變柱高和相面積分?jǐn)?shù), 汞液滴在固體銅溝槽上的潤濕特性, 并將模擬結(jié)果與Wenzel[8]和Cassie[9]理論預(yù)測結(jié)果對比, 發(fā)現(xiàn)液滴在疏水表面潤濕的接觸角模擬值與理論計(jì)算值符合較好, 然而液滴在親水表面潤濕的接觸角模擬值較理論計(jì)算值偏大, 造成該現(xiàn)象的原因是納米尺度下親水表面高的邊緣密度使得釘扎效應(yīng)明顯.Saha等[18]研究了水滴在不同幾何參數(shù)(柱高、柱寬、柱間距)的石墨表面上的潤濕特性,發(fā)現(xiàn)在最佳柱間距下, 相較水滴在方柱陣列表面上的潤濕, 水滴在圓柱陣列表面的潤濕具有更高的流動性.Chen等[15]研究了不同參數(shù)(柱高、柱間距和水滴尺寸)對Wenzel狀態(tài)向Cassie狀態(tài)自發(fā)轉(zhuǎn)換的影響.結(jié)果表明, 水滴在硅表面上潤濕狀態(tài)的轉(zhuǎn)換分為固-液接觸面積減小和水滴上升兩個(gè)階段.當(dāng)柱高位于21.72—48.87 ?時(shí), 上述潤濕轉(zhuǎn)換過程發(fā)生的概率隨液滴尺寸的減小和柱間距的增大而減小.
目前文獻(xiàn)對液滴潤濕行為的研究主要集中于欄柵形、方柱陣列等規(guī)則形狀的納米結(jié)構(gòu)[15,20], 而低表面張力的流體容易滲入凹槽形成黏滯力較強(qiáng)的潤濕接觸, 導(dǎo)致液滴不易從表面脫離, 給超疏水表面的制備帶來困難, 因此增強(qiáng)能量勢壘實(shí)現(xiàn)復(fù)合接觸成為一個(gè)重大的挑戰(zhàn).可重入結(jié)構(gòu)(如“T”或蘑菇形狀及懸垂形結(jié)構(gòu)等)逐漸成為研究者關(guān)注的焦點(diǎn).Savoy和Escobedo[26]采用MD方法模擬了不同尺寸的液滴在“T”形結(jié)構(gòu)上的潤濕行為, 發(fā)現(xiàn)潤濕與方柱相同高度的該結(jié)構(gòu)表面需要克服更高的能量勢壘, 一定程度上增強(qiáng)了表面疏水性.Wang等[27]通過實(shí)驗(yàn)方法研究了不同表面邊緣角和工質(zhì)對液滴潤濕特性的影響, 結(jié)果表明, 尖銳的表面邊緣帶來的強(qiáng)釘扎作用, 可以顯著改變液滴在表面的接觸角和潤濕狀態(tài).Zhu等[28]采用MD方法研究了不同參數(shù)(固體底角、表面浸潤性)對水滴潤濕狀態(tài)轉(zhuǎn)換的影響, 研究表明, Cassie態(tài)向Wenzel態(tài)轉(zhuǎn)換的臨界本征接觸角對固體底角有很強(qiáng)的依賴性, 而Wenzel態(tài)向Cassie態(tài)轉(zhuǎn)換的臨界本征接觸角與固體底角呈相反的規(guī)律, 并通過模擬宏觀降雨計(jì)算了潤濕轉(zhuǎn)換的自由能壘, 定量評估Cassie狀態(tài)的穩(wěn)定性.本文提出了三種不同形貌(縱截面分別為倒梯形、長方形和正梯形)的納米結(jié)構(gòu)表面,采用MD方法研究納米結(jié)構(gòu)傾角(60°, 90°和120°)及表面浸潤性對液滴潤濕行為的影響, 為提高Cassie狀態(tài)的穩(wěn)定性、制備親水材料的超疏水表面和強(qiáng)化滴狀凝結(jié)換熱提供一定指導(dǎo)意義.本模擬采用開源MD模擬軟件LAMMPS實(shí)現(xiàn), 原子位型可視化采用OVTTO軟件.
本文采用MD方法研究了納米液滴在不同粗糙表面的潤濕行為.圖1(a)為初始狀態(tài)的MD模型.模擬體系的大小為Lx×Ly×Lz= 19.58 nm ×19.58 nm × 17.62 nm,x和y方向采用周期性邊界條件,z方向采用固定邊界條件.整個(gè)體系溫度恒定為100 K (本文計(jì)算側(cè)重于亞臨界壓力下的氣-液相平衡, 該溫度低于氬的臨界溫度150.7 K,100 K對應(yīng)的飽和壓力為0.32377 MPa, 同樣低于氬的臨界壓力4.863 MPa).圖1(b)為物理模型的俯視圖, 其中紅、黃色原子代表固體壁面的Pt原子, 固體基板厚度δ為3σ, 共6層原子.按照面心立方晶格結(jié)構(gòu)排布, 密度為21.45 × 103kg/m3, 晶格常數(shù)為1.15σ, 其(111)晶面與液體原子接觸.納米結(jié)構(gòu)形狀如圖1(c)所示, 對應(yīng)三種傾角φ(60°,90°和120°), 固體原子數(shù)為28806—43958個(gè).藍(lán)色原子代表液滴原子Ar, 因?yàn)锳r屬于單原子分子,氬原子間的Lennard-Jones (LG) 勢函數(shù)具有廣泛的理論和實(shí)驗(yàn)研究基礎(chǔ), 可定性地表征水滴在固體表面的潤濕模擬[22,29].液滴的初始結(jié)構(gòu)為類固體原子排布的面心立方晶格晶格結(jié)構(gòu), 在體系溫度下的飽和密度為1.31 × 103kg/m3, 晶格常數(shù)為1.72σ.初始狀態(tài), 氬原子全部設(shè)置為液體原子, 相空間中氣相氬原子數(shù)目為零.根據(jù)納米尺度與原子尺度量級的比較(D≥ 10σ), 本模擬中液滴直徑選取D=26σ(8.86 nm)[25,27], 原子數(shù)量為7131個(gè), 對應(yīng)的液滴體積為364.2 nm3.
圖1 (a)初始狀態(tài)下的模擬系統(tǒng)圖; (b)系統(tǒng)俯視圖; (c) 納米結(jié)構(gòu)示意圖Fig.1.(a) Diagram of a simulated system in initial state; (b) top view of system; (c) schematic diagram of nanostructure.
前文依據(jù)Wenzel[8]及Cassie和Baxter[9]的理論分別給出了粗糙度因子r及相面積分?jǐn)?shù)f的定義, 并在圖2(a)—(c)中利用不同顏色的線段顯示, 其中Asl代表固-液實(shí)際接觸面積,Alv代表氣-液接觸面積,Apr代表投影面積.下面結(jié)合圖2(d)對相關(guān)參數(shù)作如下設(shè)定, 黑色虛線包圍區(qū)域作為一個(gè)單元, 單個(gè)納米結(jié)構(gòu)頂部寬度設(shè)為w, 即上表面圓直徑, 高度設(shè)為h, 相鄰柱間頂部寬度為s, 納米結(jié)構(gòu)傾角為φ, 粗糙因子r及相面積分?jǐn)?shù)f計(jì)算如下:
圖2 (a) Cassie, (b) Partial Wenzel, (c) Wenzel潤濕狀態(tài)的示意圖; (d)粗糙度因子r及相面積分?jǐn)?shù)f確定的示意圖Fig.2.Sketches for (a) Cassie, (b) Partial Wenzel and (c) Wenzel wetting states; (d) schematic diagram for determining roughness factor r and phase area fraction f.
將各結(jié)構(gòu)的幾何參數(shù)代入(3)式和(4)式, 可得不同φ的粗糙因子r和相面積分?jǐn)?shù)f, 具體見表1.
表1 三種粗糙表面對應(yīng)的結(jié)構(gòu)參數(shù)Table 1.Structural parameters corresponding to three rough surfaces.
不同種原子間的相互作用采用修正后的LJ勢能模型, 表達(dá)式為[30]
式中r為原子間距, 由Lorentz-Berthelot守則[31]可得能量參數(shù)特征尺寸σsl=(σs+σl)/2.調(diào)節(jié)能量系數(shù)α和尺度系數(shù)β可以改變固-液間的相互作用力, 得到不同的浸潤性.其中不同類型的原子對應(yīng)的特征尺寸分別為σs=0.2475 nm,σl= 0.3405 nm; 能量參數(shù)分別為εs=8.35 × 10—20J,εl= 1.67 × 10—21J.模擬過程中的截?cái)喟霃皆O(shè)為3.0σ.對于距離大于3.0σ的原子, 原子間相互作用可忽略不計(jì).
采用Velocity-Verlet算法求解運(yùn)動方程, 時(shí)間步長取Δt= 0.002346τ, 其中為LJ約化時(shí)間常數(shù).模擬過程對系統(tǒng)采用正則系綜(NVT), 維持體系溫度為定值T= 0.827ε/kB,kB為玻爾茲曼常數(shù), 采取Langevin方法控溫.體系運(yùn)行100萬步(對應(yīng)時(shí)間5 ns)后達(dá)到平衡.為保證數(shù)據(jù)的有效性, 平衡后再繼續(xù)運(yùn)行100萬步對相關(guān)結(jié)果進(jìn)行統(tǒng)計(jì)平均.
通過調(diào)節(jié)固-液勢能參數(shù)(α和β)可以改變表面浸潤性, 圖3(a)為體系內(nèi)無量綱勢能隨粒子間距的變化曲線.勢能φ隨原子間距r的增大先減小后增大, 最后趨于穩(wěn)定.勢阱深度代表固體原子對液體原子的束縛程度, 隨固-液勢能參數(shù)的增大而增大.當(dāng)吸收的能量不足以克服勢能壁壘的限制時(shí), 液體原子在平衡位置附近振動, 呈類固體排列.圖3(b)為本征接觸角θe(液滴站立在光滑表面上時(shí), 在氣-液-固三相交點(diǎn)處作氣-液界面切線, 此切線與固-液界線間的夾角)隨固-液作用參數(shù)β的變化曲線.結(jié)果顯示,θe與β呈線性關(guān)系.本文選取了七種不同的浸潤性,α= 0.14,β分別為0.3,0.35, 0.4, 0.45, 0.5, 0.55和0.6, 對應(yīng)的θe分別為157.6°, 143.5°, 129.4°, 114.5°, 101.3°, 84.6°和68.9°.
圖3 (a)無量綱勢能隨粒子間距的變化曲線; (b)本征接觸角隨表面浸潤性的變化曲線Fig.3.(a) Variation curves of dimensionless potential energy with particle spacing; (b) change curves of intrinsic contact angle with surface wettability.
圖4 給出了三種潤濕模式下液滴的密度云圖及接觸角測量方法.不同顏色表示不同的數(shù)密度范圍, 通過擬合連接氣-液界面的等密度點(diǎn), 得到等密度曲線, 圖中藍(lán)色實(shí)線代表液體主體密度與氣體主體密度的平均值, 即ρ*= (ρl+ρv)/2.紅色區(qū)域代表固體基板與納米結(jié)構(gòu), 藍(lán)色虛線代表對氣-液界面進(jìn)行圓形擬合得到的液滴外形輪廓線.不同潤濕模式均選取納米結(jié)構(gòu)頂部與液滴接觸部分作為測量基準(zhǔn)(圖4(d)—(f)中黑色實(shí)線所示), 與液滴外部輪廓線做切線, 與水平面夾角即為液滴在粗糙表面的表觀接觸角θ, 簡化表達(dá)統(tǒng)一寫成接觸角.圖4(e)中標(biāo)示的D代表液滴在粗糙表面的滲透深度, 滲透深度與納米結(jié)構(gòu)高度之比為滲透率P, 表示液滴的潤濕程度.采用文獻(xiàn)[32]中液滴潤濕狀態(tài)的劃分方法, 以潤濕程度P區(qū)分液滴潤濕狀態(tài), 其中Cassie模式,P≤ 20%; Wenzel模式,P= 100%;Partial Wenzel模式, 20% <P< 100%.
圖4 液滴在粗糙表面上三種潤濕模式下的(a)—(c)密度云圖及(d)—(f)接觸角Fig.4.(a)—(c) Density nephograms and (d)—(f) contact angles of droplets in three wetting modes on rough surfaces.
納米液滴在潤濕過程中, 表面張力和分子間作用力是液滴變形與移動的驅(qū)動力, 控制液滴在垂直及水平方向的移動.初始模擬環(huán)境為真空狀態(tài), 由于空間內(nèi)壓差及分子間力的作用, 在液滴演化過程中, 氬原子逐漸離開液滴進(jìn)入氣相區(qū)域發(fā)生擴(kuò)散,隨著氣相原子數(shù)目不斷增多, 氣相空間壓力趨于飽和, 液滴處于穩(wěn)定狀態(tài).圖5(a)—(c)代表不同θe和φ對應(yīng)的液滴潤濕平衡位型圖.當(dāng)θe= 157.6°時(shí), 液滴潤濕過程中固-液接觸面積較小, 導(dǎo)致表面張力在垂直方向的分量遠(yuǎn)大于固-液相互作用力,使得液滴無法穩(wěn)定在固體表面而發(fā)生脫離(Dewetting狀態(tài)).王文輝[33]研究結(jié)果中也提到了類似的現(xiàn)象.
為表示不同φ和θe下液滴在固體表面的潤濕程度, 圖5(d)給出了滲透率(P)隨浸潤性(θe)的變化曲線.隨著θe的減小, 分子間作用力逐漸占主導(dǎo)地位,P由0增至100%, 液滴逐漸向Wenzel模式過渡.隨著φ的減小,P不斷減小, 一定程度上降低了液體原子向納米結(jié)構(gòu)凹槽滲透的能力.
圖5 (a)—(c)當(dāng)φ = 60°, 90°, 120°時(shí), 不同θe的平衡狀態(tài)位型圖; (d) 當(dāng)φ = 60°, 90°, 120°時(shí), 滲透率隨表面浸潤性的變化曲線Fig.5.(a)—(c) Bitmaps of equilibrium states for different θe with φ = 60°, 90°, 120°; (d) permeability curves with surface wettability with φ = 60°, 90°, 120°.
在宏觀實(shí)際中, 固體表面的各種微小結(jié)構(gòu)對潤濕性能都會產(chǎn)生一定影響, 液滴在固體表面的潤濕狀態(tài)和接觸角會隨著表面浸潤性、幾何尺寸和納米結(jié)構(gòu)傾角的改變而發(fā)生變化.如圖6(a)所示, 在相同納米結(jié)構(gòu)傾角(φ)下, 固體表面上的θ隨著θe的增加而增加.而在相同浸潤性(θe)下, 按照潤濕模式和液體進(jìn)入凹槽的數(shù)目Np(Np的統(tǒng)計(jì)方法為: 截取沿z方向7σ的位型圖對應(yīng)的俯視圖), 可將θ隨著φ的變化劃分為I和II兩類.
在I類情況下,θ隨著φ的增加而增加.當(dāng)θe=84.6° (潤濕模式與Np相同)時(shí), 液滴在固體表面呈Wenzel狀態(tài)(圖5(a)—(c)所示), 對應(yīng)的Np= 4(圖6(b)所示), 此時(shí)θ60°<θ90°<θ120°.造成上述現(xiàn)象的原因是, 單個(gè)凹槽的體積隨著φ的增加而減小, 導(dǎo)致納米結(jié)構(gòu)內(nèi)部可容納的總原子數(shù)(Ntotal)隨φ的增加而減小(圖6(c)所示).
在II類情況下,θ隨著φ的增加而減小.當(dāng)θe= 129.4°(潤濕模式不同、Np相同)時(shí), 隨著φ的增加, 潤濕模式依次為Cassie, Partial Wenzel和Partial Wenzel, 對應(yīng)的Np= 1 (圖6(b)), 此時(shí),θ60°>θ90°>θ120°.主要原因?yàn)?φ對液體原子向凹槽內(nèi)部的滲透能力有一定影響, 且Ntotal隨著φ的增加逐漸增加(圖6(c)); 當(dāng)θe= 101.3° (潤濕模式相同、Np不同)時(shí), 液滴在固體表面呈Wenzel狀態(tài)(圖5(a)—(c)), 對應(yīng)的Np,90°=Np,120°= 4和Np,60°= 1 (圖6(c)), 此時(shí),θ60°>θ90°>θ120°.造成上述現(xiàn)象的原因是: 當(dāng)液滴外緣位于納米結(jié)構(gòu)邊緣時(shí), 高的表面邊緣密度導(dǎo)致近壁面處接觸線固定, 限制了液體原子的水平移動.
圖6 (a)當(dāng)φ = 60°, 90°, 120°時(shí), 表觀接觸角隨θe的變化曲線; (b)液體潤濕凹槽數(shù)目隨θe的變化曲線, 其中, φ = 90°和120°時(shí)的曲線完全重合; (c)液體原子滲入到凹槽的數(shù)量隨θe的變化曲線; (d)液滴基底圓半徑隨θe的變化曲線Fig.6.(a) Curves of apparent contact angle with θe when φ = 60°, 90°, 120°; (b) change curves of the number of liquid wetting grooves with θe, where the curves at φ = 90° and 120° are completely coincident; (c) amount of liquid atoms entering the groove varies with θe; (d) change curves of the droplet base circle radius with θe.
此外由圖6(d)可知, 對于相同納米結(jié)構(gòu)傾角(φ)的固體表面, 基底圓半徑B(指納米結(jié)構(gòu)頂部固-液接觸面積對應(yīng)的半徑)隨θe的減小而增大,表明浸潤性越強(qiáng), 液滴在粗糙表面的鋪展程度越大.對于相同浸潤性的固體表面,B隨φ的增大而增大.但φ= 120°與60°時(shí),B存在突降, 如圖6(d)中1和2標(biāo)識所示, 其中標(biāo)識2對應(yīng)的表面浸潤性下, 呈現(xiàn)相反的變化規(guī)律.出現(xiàn)該現(xiàn)象的原因是不同φ對應(yīng)的潤濕狀態(tài)或Np不同.當(dāng)θe= 143.5°時(shí)(圖6(d)中的標(biāo)識1),φ= 60°與90°對應(yīng)的為Cassie狀態(tài), 而φ= 120°對應(yīng)的為Dewetting狀態(tài); 當(dāng)θe= 101.3°時(shí)(圖6(d)中的標(biāo)識2),φ= 90°與120°對應(yīng)的Np= 4, 而φ= 60°對應(yīng)的Np= 1.
3.2.1 理論自由能的計(jì)算
本節(jié)中建立三維空間內(nèi)液滴在粗糙表面的潤濕模型, 假設(shè)液滴分布為球型, 每種狀態(tài)的總能量為液滴和其接觸表面之間所產(chǎn)生的所有能量之和.系統(tǒng)的自由能差可表示為
固體表面上的Young方程為
結(jié)合(6)式和(7)式得到:
(6)式對于微觀系統(tǒng)的能量計(jì)算同樣適用[34,35].值得指出的是, 納米尺度下重力對液滴動力學(xué)行為的影響很小, 所以本文忽略了納米液滴的重力.即公式中的重力項(xiàng)Fg= 0.其中, ΔEint為固體表面上有無液滴時(shí)的固-液-氣界面能差值;Eg為重力勢能;σlv為氣-液表面張力,Alv為氣-液接觸面積,Asl為固-液接觸面積.
對于Cassie模型:
對于Wenzel模型:
這里R為液滴半徑,B為基底圓半徑, 表達(dá)式為
3.2.2 理論預(yù)測與模擬結(jié)果的比較
圖7(a)—(c)為理論自由能的計(jì)算結(jié)果, 且對比MD模擬結(jié)果(圖5), 發(fā)現(xiàn)二者符合較好.當(dāng)θe= 129.4°,φ= 90°時(shí), 理論預(yù)測結(jié)果顯示ΔEW> ΔEC, 基于自由能最低原理判斷液滴呈Cassie潤濕模式, MD模擬結(jié)果同樣得出該表面浸潤性下液滴呈Cassie潤濕模式.
圖7 (a)—(c)不同φ和θe下的理論自由能值; (d)—(f)不同φ和θe下接觸角的變化Fig.7.(a)—(c) Theoretical free energy values under different φ and θe; (d)—(f) change of contact angle under different φ and θe.
由圖7(d)—(f)發(fā)現(xiàn),θ隨θe增加而增加, 變化趨勢與φ無關(guān), 但數(shù)值上存在略微差別.通過對處于C (Cassie)和W (Wenzel)兩種區(qū)域內(nèi)的模擬值(θSIM)與對應(yīng)潤濕模式的預(yù)測值(θC,θw)進(jìn)行比較, 發(fā)現(xiàn)在不同φ值下,θC都近似等于θSIM; 在大多數(shù)情況下θw與θSIM符合良好.但θe= 68.9°時(shí)θSIM>θw, 此時(shí)液滴外緣位于納米結(jié)構(gòu)邊緣(圖5), 高表面邊緣密度使得接觸線固定, 釘扎作用顯著, 但Wenzel與Cassie理論模型未考慮上述釘扎效應(yīng)的影響, 所以理論預(yù)測值偏小.而處于PW (Partial Wenzel)區(qū)域的θSIM與θC,θw的大小關(guān)系為θw<θSIM<θC, 且滲透程度越高,θSIM越接近θw.通過理論公式計(jì)算的C與W潤濕模式的轉(zhuǎn)換接觸角θCW近乎相等(θCW,60°= 123.6°,θCW,90°= 121.4°,θCW,120°=125.8°), 均位于PW區(qū)域.
3.2.3 不同浸潤性及納米結(jié)構(gòu)傾角下的潤濕模式分區(qū)
為了討論不同θe和φ對潤濕模式的影響, 對潤濕模式進(jìn)行統(tǒng)計(jì)(如圖8所示).研究發(fā)現(xiàn), 可將φ-θe相圖劃分為4個(gè)區(qū)域: C, W, PW和DW(Dewetting), 隨著θe的增加, 液滴逐漸實(shí)現(xiàn)W→PW→C→DW的潤濕轉(zhuǎn)換.
圖8 不同φ和θe下的潤濕模式分區(qū)圖Fig.8.Wetting mode partition diagram under different φ and θe.
在60° <θe< 118°時(shí), 液滴在納米結(jié)構(gòu)表面呈現(xiàn)W狀態(tài), 該潤濕狀態(tài)在相圖中占比最高, 約為總區(qū)域面積的1/2, 在該區(qū)域內(nèi)表面浸潤性占主導(dǎo)作用(即固-液間作用力較強(qiáng)).隨著θe的增加(125° <θe< 145°), 液滴在納米結(jié)構(gòu)表面呈現(xiàn)C與PW狀態(tài), 且C與PW區(qū)域分別呈三角形和反梯形分布, 兩種潤濕區(qū)域所占面積近乎相等.其中C區(qū)域面積隨φ的減小而增大, 表明φ減小增強(qiáng)了由銳邊效應(yīng)引起的釘扎作用, 使得壁面外緣可以容納更多的流體原子; 但當(dāng)φ較大時(shí)(φ≥120°)不存在C狀態(tài), 微柱陣列表面的液體原子在毛細(xì)力和固-液相互作用力[27]的共同作用下, 更容易滲透至納米結(jié)構(gòu)內(nèi)部, 不利于形成超疏水表面.在θe> 145°時(shí), 液滴在納米結(jié)構(gòu)表面呈現(xiàn)DW狀態(tài), 表面張力占主導(dǎo)作用(即固-液間作用力較弱), 且發(fā)現(xiàn)在該區(qū)域內(nèi)固體原子數(shù)目隨φ的增加而增加.
采用MD方法研究不同納米結(jié)構(gòu)傾角(φ)及表面浸潤性(θe)下, 結(jié)構(gòu)尺寸與液滴大小相當(dāng)時(shí),納米氬液滴在粗糙鉑表面的潤濕行為.討論表觀接觸角的影響因素及檢驗(yàn)是否符合Wenzel及Cassie理論預(yù)測模型[8,9], 根據(jù)不同結(jié)構(gòu)及浸潤性進(jìn)行潤濕模式分區(qū).具體結(jié)論如下:
1)液滴在不同納米結(jié)構(gòu)傾角(60°, 90°和120°)及浸潤性的表面達(dá)到潤濕平衡后存在四種潤濕狀態(tài), 分別為Dewetting, Cassie, Partial Wenzel和Wenzel, 隨著表面浸潤性的減弱, 潤濕狀態(tài)逐漸由Wenzel向Dewetting過渡;
2)通過比較表觀接觸角模擬值(θSIM)與理論預(yù)測值(θw,θC)的大小關(guān)系, 發(fā)現(xiàn)不同納米結(jié)構(gòu)傾角在Cassie與Wenzel潤濕狀態(tài)的表觀模擬值與預(yù)測值基本符合, 而Partial Wenzel潤濕狀態(tài)由于部分液體原子滲入凹槽內(nèi), 模擬值較預(yù)測值偏小(θSIM<θC), 且滲透程度越高θSIM越接近θw;
3)不同潤濕模式基本滿足自由能最小原理,根據(jù)φ和θe對潤濕模式進(jìn)行分區(qū): 隨著浸潤性θe的減小,φ對最終潤濕狀態(tài)的影響減弱, 尤其當(dāng)θe< 118°時(shí), 液滴在納米結(jié)構(gòu)表面均呈Wenzel狀態(tài); 對于超疏水表面(118° <θe< 145°), 倒梯形(φ= 60°)納米結(jié)構(gòu)有利于保持Cassie狀態(tài), 而正梯形(φ= 120°)納米結(jié)構(gòu)容易形成Partial Wenzel液滴, 其一由于φ減小增強(qiáng)了銳邊效應(yīng)引起的釘扎作用, 相同浸潤性下降低了液體原子向凹槽內(nèi)滲透的能力, 其二是當(dāng)液滴擴(kuò)散至納米結(jié)構(gòu)外緣時(shí),可以容納更多的流體原子.