張福建 陳悅 高翔 劉珍 張忠強(qiáng)3)?

1) (江蘇大學(xué)智能柔性機(jī)械電子研究院, 鎮(zhèn)江 212013)

2) (江蘇科技大學(xué)船舶與海洋工程學(xué)院, 鎮(zhèn)江 212003)

3) (常州大學(xué), 江蘇省光伏科學(xué)與工程協(xié)同創(chuàng)新中心, 常州 213164)

(2021 年5 月13 日收到; 2021 年5 月28 日收到修改稿)

1 引 言

液滴定向自驅(qū)動(dòng)在自然界中廣泛存在, 例如沙漠甲蟲利用背部的凸起結(jié)構(gòu)收集液滴[1], 豬籠草吻部特殊的溝槽結(jié)構(gòu)實(shí)現(xiàn)水的單向運(yùn)輸[2]等. 明確液滴自驅(qū)動(dòng)的特性和機(jī)理對(duì)于實(shí)現(xiàn)水收集[3,4]、固體表面自清潔、微流控芯片[5]等應(yīng)用具有重要價(jià)值.目前, 液滴驅(qū)動(dòng)主要分為液滴自驅(qū)動(dòng)和外場(chǎng)作用下的驅(qū)動(dòng). 外場(chǎng)作用下的驅(qū)動(dòng)具有需要外場(chǎng)輸入、控制復(fù)雜等特點(diǎn)[6-8], 在宏觀尺度存在能量消耗, 在微納尺度難以精確添加外場(chǎng). 而液滴自驅(qū)動(dòng)則是依靠基底的形貌梯度[9,10]、表面潤(rùn)濕梯度[11]、剛度梯度[12]等來驅(qū)動(dòng)液滴, 由基底自身的潤(rùn)濕梯度決定液滴自驅(qū)動(dòng)的方向, 無需外場(chǎng)輸入, 實(shí)現(xiàn)了簡(jiǎn)單高效的液滴驅(qū)動(dòng). 例如, 梯度密度排列的柱狀結(jié)構(gòu)實(shí)現(xiàn)水滴從柱稀疏端向密集端移動(dòng)[13], 錐形曲率結(jié)構(gòu)實(shí)現(xiàn)水滴從圓錐尖端向尾端的自主運(yùn)輸[14,15]. 然而, 基底表面缺陷產(chǎn)生的液滴移動(dòng)接觸線釘扎通常會(huì)影響液滴的運(yùn)動(dòng)[16-18]. 基于豬籠草吻部和U 型島陣列的非對(duì)稱結(jié)構(gòu)可以利用釘扎現(xiàn)象防止液體回流[19,20], 從而實(shí)現(xiàn)連續(xù)液體的定向鋪展, 但是無法驅(qū)動(dòng)單個(gè)液滴, 而基于萊頓弗羅斯特效應(yīng)克服液滴接觸線釘扎以驅(qū)動(dòng)液滴的方法又表現(xiàn)出控制困難的特點(diǎn)[21,22], 且萊頓弗羅斯特效應(yīng)需要高溫, 無法適用于生物微流控芯片、水收集等領(lǐng)域. 因此,在常溫狀態(tài)下克服釘扎阻力實(shí)現(xiàn)液滴快速自驅(qū)動(dòng)顯得尤為重要.

單層石墨烯作為一種二維熱點(diǎn)材料, 近年來已受到廣泛關(guān)注. 已有研究表明, 基于單層石墨烯的潤(rùn)濕透明性和低流-固界面摩擦[23,24], 單層石墨烯覆層在缺陷基底上可有效降低液滴進(jìn)入和脫離缺陷時(shí)的釘扎阻力[25]. 此外, 水滴在單層石墨烯上的快速擴(kuò)散已被報(bào)道[26], 進(jìn)一步地, 液滴在單層石墨烯覆層的柱狀和凹槽銅基底上實(shí)現(xiàn)單向快速自驅(qū)動(dòng)[27,28],表明單層石墨烯在液滴自驅(qū)動(dòng)領(lǐng)域的巨大潛力.

本文提出了一種楔形銅基底-單層石墨烯覆層結(jié)構(gòu)(monolayer graphene-covered wedge-shaped copper substrate, GWCS), 用以實(shí)現(xiàn)液滴的單向超快速自驅(qū)動(dòng). 采用分子動(dòng)力學(xué)方法研究了液滴在GWCS 上的自驅(qū)動(dòng)行為, 探索了基底表面潤(rùn)濕梯度和Young-Laplace 壓力梯度驅(qū)動(dòng)液滴的機(jī)理,明確了液滴最大位移與基底潤(rùn)濕性、基底角度和液滴表面張力之間的關(guān)系. 進(jìn)一步分析了楔形基底寬度非連續(xù)線性增加對(duì)液滴自驅(qū)動(dòng)的作用機(jī)制, 并確定了基底角度和液滴尺寸對(duì)液滴自驅(qū)動(dòng)的影響規(guī)律. 最后, 研究了楔形基底兩側(cè)凹槽深度對(duì)液滴自驅(qū)動(dòng)機(jī)理和液滴輸運(yùn)效率的影響機(jī)制. 研究結(jié)果對(duì)實(shí)現(xiàn)液滴自驅(qū)動(dòng)的單層石墨烯覆層的功能織構(gòu)表面設(shè)計(jì)具有理論指導(dǎo)意義.

2 模型和方法

分子動(dòng)力學(xué)模型由帶有楔形結(jié)構(gòu)的銅基底、單層石墨烯和納米水液滴組成, 銅原子尺度σCu-Cu=0.2314 nm[29], 碳原子尺度σC-C= 0.340 nm[30],單個(gè)水分子直徑約為0.385 nm, 模擬盒子大小為60 nm × 33 nm × 16 nm. 銅基底楔形結(jié)構(gòu)的楔角為α, 長(zhǎng)度為34 nm, 楔形銅基底左側(cè)延伸出長(zhǎng)約4 nm、寬約1.3 nm 的長(zhǎng)方體結(jié)構(gòu), 用于在放置液滴時(shí)起到支撐作用, 為了使楔形銅基底兩側(cè)的石墨烯不被吸附進(jìn)楔形基底兩側(cè)的凹槽, 將凹槽深度設(shè)置為2 nm，如圖1 所示. 將單層石墨烯覆蓋在帶有楔形結(jié)構(gòu)的銅基底上0.2857 nm (C-Cu 原子間的平衡距離)處. 將銅基底楔角底部的底板、楔角左側(cè)的長(zhǎng)方體結(jié)構(gòu)和x方向上的單層石墨烯兩端各1 nm 設(shè)為剛性, 弛豫0.5 ns, 使得單層石墨烯吸附在銅基底上. 水納米液滴的直徑為d, 水分子采用TIP4P-EW模型[31], 將水納米液滴弛豫0.5 ns, 使其達(dá)到能量最小化的平衡狀態(tài), 然后將液滴的質(zhì)心放置在銅基底楔角結(jié)構(gòu)最左側(cè)的正上方, 液滴底部距離單層石墨烯0.328 nm (C-O 原子間的平衡距離)[27].

圖1 水納米液滴在GWCS 上定向運(yùn)輸模型 (a)楔形銅基底模型的俯視圖和A-A 方向剖視圖, 楔形銅基底包括中間的楔形基底、楔形基底兩側(cè)的凹槽和凹槽底部底板; (b)水納米液滴在GWCS 上的俯視圖和正視圖. 粉色是銅基底, 黑色是單層石墨烯, 紅白相間的球是水納米液滴Fig. 1. Uni-directional self-driving model of a water nanodroplet on GWCS: (a) Top view and A-A sectional view of wedge-shaped copper substrate, the wedge-shaped copper substrate includes wedge-shaped substrate in the middle, grooves on both sides of the wedge-shaped substrate and a copper plate at the bottom of the grooves; (b) top view and front view of water nanodroplet on GWCS. The pink is Cu substrate, the black is monolayer graphene, and the red and white sphere is water nanodroplet.

模擬采用開源軟件LAMMPS, 使用正則系宗將模擬溫度控制在298 K, 時(shí)間步長(zhǎng)為L(zhǎng)AMMPS模擬的最小時(shí)間步長(zhǎng)1 fs. 碳原子之間的相互作用由第二代REBO 勢(shì)函數(shù)來描述, 銅原子之間的相互作用采用EAM 來描述. 水分子、銅、石墨烯之間的相互作用由Lennard-Jones (LJ)勢(shì)函數(shù)來描述, 不同原子間相互作用勢(shì)函數(shù)參數(shù)由Lorentz-Berthelot混合準(zhǔn)則得到, 勢(shì)函數(shù)參數(shù)如表1 所列[29,30,32]. LJ相互作用和庫侖相互作用的截?cái)喟霃骄O(shè)置為12 ?,長(zhǎng)程靜電相互作用采用Particle-Particle Particle-Mesh 算法計(jì)算. 氫原子和其他原子之間的平衡距離和勢(shì)阱深度設(shè)置為0[31].

表1 不同原子之間的LJ 勢(shì)能參數(shù)Table 1. Parameters of LJ potential between different atoms.

3 結(jié)果和討論

首先研究了直徑為8 nm 的水納米液滴(8960 個(gè)水分子)在楔角α= 40°的GWCS 上的自驅(qū)動(dòng)行為. 液滴的速度、位移隨時(shí)間的變化關(guān)系如圖2 所示. 在無任何外場(chǎng)的作用下, 液滴沿著楔角寬度增加的方向運(yùn)動(dòng), 液滴的運(yùn)動(dòng)速度先隨著時(shí)間的增加而增大, 在t= 0.358 ns 時(shí), 速度達(dá)到最大值56.2 m/s (202.32 km/h). 而后液滴速度快速減小, 并在t≈ 0.68 ns 時(shí)減小到0, 此時(shí)液滴的最大位移為21.1 nm, 液滴質(zhì)心所在位置的楔形基底寬度為16.2 nm. 最后, 當(dāng)液滴位移達(dá)到最大值時(shí), 液滴完全處于楔形基底上面, 由于液滴在前進(jìn)方向上的彈性形變, 液滴發(fā)生原地振蕩, 位移略有回落.以上結(jié)果表明GWCS 可以實(shí)現(xiàn)液滴從楔形基底尖端向?qū)挾说亩ㄏ蜃则?qū)動(dòng).

圖2 液滴自驅(qū)動(dòng)的速度、位移與時(shí)間的關(guān)系, 液滴直徑為8 nm, 楔形基底角度為40°, 插圖分別為液滴在0.35 和1.00 ns 時(shí)的運(yùn)動(dòng)狀態(tài), 單層石墨烯模型設(shè)置為線條結(jié)構(gòu)且未顯示銅基底底板Fig. 2. Velocity and displacement of the water droplet selfdriving on GWCS. The diameter of the water nanodroplet is 8 nm, and the wedge angle of GWCS is 40°, the illustrations show the motion of the droplet at 0.35 and 1.00 ns, respectively. The monolayer graphene model is set to linear structure and the bottom plate of copper substrate does not show.

為了揭示水納米液滴在GWCS 上的自驅(qū)動(dòng)行為, 計(jì)算了水納米液滴在前進(jìn)方向上的受力F和液滴與GWCS 之間的相互作用能隨時(shí)間的變化關(guān)系, 如圖3 所示. 液滴在前進(jìn)方向上的受力F決定了液滴的速度變化, 力F從正數(shù)開始, 先增大后減小, 當(dāng)F減小到0 時(shí), 液滴速度達(dá)到最大(圖3(a)).然后, 力F繼續(xù)減小, 導(dǎo)致液滴速度快速下降. 最后, 隨著液滴速度的減小, 力F逐漸增加至0 附近,液滴速度最終減小到0. 為了進(jìn)一步探索液滴的受力來源, 計(jì)算了液滴與GWCS 之間的相互作用能隨時(shí)間的變化關(guān)系, 如圖3(b)所示. 在液滴從楔形基底的尖端向?qū)挾诉\(yùn)動(dòng)的過程中, 液滴與銅基底之間的相互作用能Ew-Cu逐漸減小, 表明液滴可以透過單層石墨烯與銅基底發(fā)生相互作用, 印證了單層石墨烯的潤(rùn)濕透明性. 然而, 液滴與銅基底之間的相互作用能Ew-Cu的變化量遠(yuǎn)小于液滴與單層石墨烯之間的相互作用能變化量, 表明單層石墨烯在液滴自驅(qū)動(dòng)過程中占主導(dǎo)作用. 同時(shí), 隨著液滴位移的增大, 液滴所在位置的楔形基底寬度逐漸增大, 導(dǎo)致液滴與GWCS 之間的相互作用能Ew-GCu逐漸減小, 減小的流-固界面相互作用能為液滴的自驅(qū)動(dòng)提供能量. 當(dāng)液滴達(dá)到最大位移后, 液滴完全處于楔形基底上, 液滴與基底之間的相互作用能基本不再變化. 此外, 由于液體的表面張力, 液滴表面積的變化會(huì)改變液滴的勢(shì)能. 在液滴自驅(qū)動(dòng)過程中, 液滴的勢(shì)能呈現(xiàn)出增加的趨勢(shì), 當(dāng)液滴達(dá)到最大位移后, 液滴原地振蕩, 勢(shì)能在一定范圍內(nèi)波動(dòng), 趨勢(shì)上不再增加. 考慮到液滴總是往能量降低的方向運(yùn)動(dòng), 而在液滴達(dá)到最大位移之前, 流-固界面相互作用能的減小和液滴勢(shì)能的增加表明流-固界面相互作用能和液滴的勢(shì)能之間存在競(jìng)爭(zhēng)關(guān)系,即減小的流-固界面能有部分轉(zhuǎn)化為液滴的勢(shì)能,并非完全作用于液滴驅(qū)動(dòng).

圖3 (a) 液滴前進(jìn)方向所受的合力F 和速度隨時(shí)間的變化關(guān)系; (b) 液滴的勢(shì)能、液滴與GWCS 之間的相互作用能隨時(shí)間的變化關(guān)系. Ew-Cu, Ew-G 和Ew-GCu 分別為液滴與銅基底、石墨烯、GWCS 之間的相互作用能, Ew 為液滴的勢(shì)能Fig. 3. (a) Resultant force F and velocity of the droplet in +x direction as a function of time; (b) potential energy of water droplet and the interfacial interaction energy between water droplet and GWCS versus time. Ew-Cu is the interaction energy between water droplet and copper substrate, Ew-G is the interaction energy between water droplet and graphene, Ew-GCu is the interaction energy between water droplet and GWCS, Ew is the droplet potential energy.

從理論上分析, 液滴自驅(qū)動(dòng)過程中受到的驅(qū)動(dòng)力可以分為兩個(gè)階段: 表面潤(rùn)濕梯度驅(qū)動(dòng)和Young-Laplace 壓力差驅(qū)動(dòng). 當(dāng)液滴在GWCS 上的接觸面包括楔形基底和溝槽正上方的單層石墨烯時(shí)(圖4(a)), 水納米液滴在單層石墨烯、單層石墨烯覆層的銅基底上的接觸角分別為57°和29°, 所以單層石墨烯的親水性小于單層石墨烯覆層的銅基底, 表面潤(rùn)濕梯度占主導(dǎo)作用, 將液滴從溝槽正上方的單層石墨烯表面推到單層石墨烯覆層的楔形基底上, 推力F0垂直于楔形基底的邊線, 液滴從楔形基底的尖端向?qū)挾诉\(yùn)動(dòng). 此時(shí), 液滴前進(jìn)方向上受到的驅(qū)動(dòng)力Fdri可以表示為

圖4 液滴自驅(qū)動(dòng)過程中在GWCS 上的兩種接觸狀態(tài)的俯視圖和正視圖(未顯示銅基底底板, 且俯視圖中的單層石墨烯為線條結(jié)構(gòu)) (a) 液滴在GWCS 上的接觸面包括楔形基底和溝槽正上方的單層石墨烯; (b) 液滴完全處于楔形基底上Fig. 4. Top and front view of two contact states of the water droplet on GWCS in the self-driving process (the bottom plate of copper substrate is not shown, and the monolayer graphene in the top view is linear structure): (a) Contact area of the droplet on the substrate includes the wedge-shaped substrate and the graphene right above the groove; (b) droplet is completely on the wedgeshaped substrate.

式中,l為液滴與楔形基底邊線的接觸線長(zhǎng)度,θ1和θ2分別是液滴在單層石墨烯覆層的銅基底和單層石墨烯上的接觸角,γ為水的表面張力, 如圖5(a)所示. 液滴與楔形基底邊線的接觸線長(zhǎng)度l可以進(jìn)一步表示為

式中R為液滴在GWCS 上接觸線的輪廓半徑,x0為液滴所在位置到楔形基底虛擬尖端的距離.

隨著液滴位移的增加, 當(dāng)液滴開始完全處于楔形基底上時(shí), 液滴在基底上的潤(rùn)濕形態(tài)如圖4(b)所示, 液滴內(nèi)部的Young-Laplace 壓力差驅(qū)動(dòng)占主導(dǎo)作用. 液滴在楔形基底上運(yùn)動(dòng)時(shí)的Young-Laplace驅(qū)動(dòng)力Fdri可以表示為[5]

式中, ΔP為液滴內(nèi)部的Young-Laplace 壓力差,Ra和Rr分別為液滴在楔形基底上的前、后接觸線的輪廓半徑.

當(dāng)液滴完全處于楔形基底上, 且液滴完成自驅(qū)動(dòng)時(shí), 液滴在基底上的位置如圖5(b)所示, 液滴在基底上的最大位移xmax可以表示為

圖5 (a) 液滴在楔形基底上的驅(qū)動(dòng)力示意圖; (b) 液滴達(dá)到最大位移時(shí)所在基底上的位置Fig. 5. (a) Driving force diagram of the droplet on GWCS; (b) position of the droplet on the substrate when the droplet reaches the maximum displacement.

根據(jù)(1)式—(4)式可以得出在本模型中影響液滴自驅(qū)動(dòng)行為的參數(shù)有楔角和液滴直徑. 為研究銅基底楔角對(duì)液滴驅(qū)動(dòng)的影響, 楔角α選取40°,30°, 24°和20°, 分別計(jì)算了8 nm 液滴在不同楔角的GWCS 上的自驅(qū)動(dòng)速度、位移隨時(shí)間的變化關(guān)系, 結(jié)果如圖6 所示. 可以看出, 在加速階段, 液滴在GWCS 上的加速度隨著楔角α的增大而增大,且液滴達(dá)到最大速度所用的時(shí)間隨著楔角角度的增大而減小(圖6(a)). 這是因?yàn)楸砻鏉?rùn)濕梯度的方向垂直于楔形基底的邊線, 楔形基底角度越大, 表面潤(rùn)濕梯度產(chǎn)生的液滴推動(dòng)力F0在液滴前進(jìn)方向上的分力越大, 即液滴在較大角度的楔形基底上擁有更大的驅(qū)動(dòng)力Fdri. 然而, 液滴的最大位移隨著基底楔角的減小而增大, 但楔角越大, 液滴越快達(dá)到最大位移(圖6(b)). 同時(shí), 隨著楔角的減小, 液滴的速度產(chǎn)生了較大的波動(dòng), 這是因?yàn)橐旱卧谛〗嵌刃ㄐ位咨鲜艿降尿?qū)動(dòng)力更小, 導(dǎo)致自驅(qū)動(dòng)行為更易受到原子熱振蕩的影響. 此外, 楔形銅基底的寬度是每隔一段距離k增加兩層銅原子(圖7), 而不是連續(xù)的線性增加, 使得GWCS 的潤(rùn)濕梯度并非連續(xù), 造成液滴在單層石墨烯覆層的楔形銅基底上的前進(jìn)方向驅(qū)動(dòng)力會(huì)隨著楔形基底寬度的增加而產(chǎn)生跨越式改變. 受液滴尺寸的影響, 間隔距離k越大, 液滴前進(jìn)方向上受到的驅(qū)動(dòng)力波動(dòng)越大,這是造成液滴在小角度楔形基底上速度具有較大波動(dòng)的另一重要原因. 對(duì)于本文, 當(dāng)楔形基底的楔角α為20°時(shí), 間隔距離k為 1.1 nm, 相對(duì)于液滴直徑8 nm 來說已經(jīng)是一個(gè)不可忽略的重要因素,表現(xiàn)在數(shù)據(jù)上則是由于液滴前進(jìn)方向驅(qū)動(dòng)力不穩(wěn)定導(dǎo)致液滴在20°楔角的基底上的速度相比于40°楔角的基底具有較大的速度波動(dòng), 在液滴自驅(qū)動(dòng)的末端尤為明顯(圖6(a)).

圖6 楔形基底角度對(duì)液滴自驅(qū)動(dòng)的影響 (a) 液滴運(yùn)動(dòng)速度; (b) 液滴的位移隨時(shí)間的變化關(guān)系Fig. 6. Influence of the wedge angle of GWCS on the droplet self-driving behaviors: (a) Velocity and (b) displacement of the droplet on GWCS with different wedge angles versus time.

圖7 不同角度的楔形銅基底, x 方向每間隔距離k, y 方向增加兩層銅原子寬度Fig. 7. For the GWCS with different wedge angles, when the interval distance is k in x direction, the width of the wedge-shaped copper substrate increases two layers of copper atoms in y direction.

為了進(jìn)一步明確楔形基底角度對(duì)液滴自驅(qū)動(dòng)的影響, 繪制了液滴自驅(qū)動(dòng)的最終位移與所在位置楔形基底寬度的理論值與模擬值對(duì)比圖, 如圖8 所示. 理論值分別由(4)式和(5)式計(jì)算得到. 受間隔距離k的影響, 液滴需要跨過間隔距離k才能達(dá)到下一個(gè)潤(rùn)濕梯度, 導(dǎo)致液滴在小角度楔形基底上的最終位移和所在位置基底寬度的模擬值小于理論值. 隨著楔形基底角度的增大, 間隔距離k逐漸減小, 模擬值與理論值的差距逐漸縮小, 當(dāng)楔形基底角度α為40°時(shí), 位移的理論值和模擬值相近.同時(shí), 考慮到液滴運(yùn)動(dòng)速度較大, 在無潤(rùn)濕梯度的基底上仍能減速運(yùn)動(dòng)一段距離, 導(dǎo)致液滴在40°楔形基底上的最終位移和所在位置的基底寬度的模擬值略大于理論值.

圖8 液滴在不同楔角GWCS 上的自驅(qū)動(dòng)最終位移、所在位置楔形基底寬度的模擬值與理論值對(duì)比Fig. 8. Comparison of the simulation and theoretical results of the final self-driving displacement and the width of wedge-shaped substrate with different wedge angles respectively.

此外, 研究了液滴尺寸對(duì)液滴在GWCS 上自驅(qū)動(dòng)行為的影響, 計(jì)算了不同直徑的納米水液滴(6, 7, 8 和9 nm, 水分子數(shù)分別為3780, 6000, 8960和12750)在楔角α為40°的GWCS 上的速度、位移隨時(shí)間的變化關(guān)系, 如圖9 所示. 結(jié)果表明, 較小的液滴在加速階段擁有較大的加速度, 其中6 nm 液滴的最大速度可以達(dá)到73.8 nm/s, 但是較大液滴的運(yùn)動(dòng)速度更為平穩(wěn)(圖9(a)). 然而,6 和7 nm 直徑的液滴尺寸較小, 易受到原子熱振蕩的影響, 且兩者體積差別小于7, 8 和9 nm 直徑液滴之間的差距, 導(dǎo)致6 nm 直徑液滴的速度與7 nm 液滴的速度和最終位移相近. 隨著液滴直徑增加到8 和9 nm, 液滴的最終位移表現(xiàn)出隨液滴直徑的增加而增加的明顯趨勢(shì)(圖9(b)). 為了研究高速運(yùn)動(dòng)的液滴在無織構(gòu)基底(驅(qū)動(dòng)力消失后液滴所在位置的基底)上的減速運(yùn)動(dòng)對(duì)液滴最終位移的影響, 繪制了液滴最終位移的模擬值與理論值((4)式)隨液滴直徑的變化關(guān)系和模擬值與理論值的差值相對(duì)于液滴直徑的比率圖, 如圖9(c)所示.直徑在6—9 nm 范圍內(nèi)的液滴, 最終位移的模擬值與理論值的差值隨著液滴直徑的增大而減小, 表明6 nm 直徑的液滴在無織構(gòu)基底上減速運(yùn)動(dòng)的距離要大于7 nm 液滴, 這也是造成6 nm 液滴和7 nm 液滴最終位移相近的原因之一. 此外, 液滴最終位移的模擬值xsim與理論值xmax的差值相對(duì)于液滴直徑d的比率 ratio=(xsim?xmax)/d隨液滴直徑的增大而減小, 表明高速運(yùn)動(dòng)的小液滴在無織構(gòu)基底上更易產(chǎn)生滑動(dòng), 而大液滴能夠抑制這種滑動(dòng)使得液滴最終位移的模擬值與理論值更加契合.

圖9 液滴直徑在6-9 nm 范圍內(nèi)時(shí)液滴尺寸對(duì)液滴自驅(qū)動(dòng)的影響 (a) 液滴的速度和(b)液滴的位移隨時(shí)間的變化關(guān)系, 插圖為6 nm 液滴的最終位置; (c) 液滴最大位移的模擬值與理論值隨液滴直徑的變化關(guān)系和模擬值與理論值的差值相對(duì)于液滴直徑的比率Fig. 9. Influence of water droplet size on the droplet self-driving behaviors when the droplet diameter in the range of 6-9 nm:(a) Velocity and (b) displacement versus time, the illustration is the final position of the 6 nm droplet; (c) simulated and theoretical values of the droplet displacement versus droplet diameter, and the ratio of the difference between the simulated value and the theoretical value to the droplet diameter.

最后, 討論楔形銅結(jié)構(gòu)兩側(cè)的凹槽深度對(duì)液滴自驅(qū)動(dòng)行為的影響. 分別計(jì)算了凹槽深度為0.542 nm(3 層銅原子)和0.904 nm(5 層銅原子)情況下的8 nm 直徑液滴自驅(qū)動(dòng)效果. 液滴的速度、位移隨時(shí)間的變化關(guān)系如圖10 所示. 結(jié)果表明, 在凹槽深度為0.542 和0.904 nm 情況下, 液滴整體速度小于凹槽為2 nm 時(shí)的液滴速度(圖10(a)), 且在凹槽深度為0.542 和0.904 nm 情況下液滴的最大位移約為8 nm, 遠(yuǎn)小于凹槽深度為2 nm 時(shí)21.1 nm的最大位移(圖10(b)). 這是因?yàn)? 當(dāng)凹槽深度為0.504 和0.904 nm 時(shí), 單層石墨烯被吸附進(jìn)凹槽內(nèi)(圖10(c)), 造成基底的表面潤(rùn)濕梯度消失, 不會(huì)產(chǎn)生凹槽深度為2 nm 時(shí)凹槽正上方的單層石墨烯將液滴推到楔形基底上的行為. 此外, 當(dāng)單層石墨烯被吸附進(jìn)凹槽時(shí), 液滴兩側(cè)會(huì)鋪展到凹槽內(nèi), 液滴僅依靠楔形基底左側(cè)的結(jié)構(gòu)梯度和楔形基底左端的長(zhǎng)方體結(jié)構(gòu)兩側(cè)與石墨烯之間的吸附間隙產(chǎn)生的微小潤(rùn)濕梯度驅(qū)動(dòng)液滴. 因此, 凹槽深度決定了單層石墨烯是否會(huì)被吸附進(jìn)凹槽內(nèi), 進(jìn)而決定了液滴的驅(qū)動(dòng)力來源和自驅(qū)動(dòng)效率. 在單層石墨烯不被吸附進(jìn)凹槽時(shí), 液滴擁有最大的自驅(qū)動(dòng)速度和最遠(yuǎn)的輸運(yùn)距離.

圖10 楔形銅結(jié)構(gòu)兩側(cè)凹槽深度對(duì)液滴自驅(qū)動(dòng)的影響 (a) 液滴的速度和(b)液滴的位移隨時(shí)間的變化關(guān)系; (c) 8 nm 直徑液滴在凹槽深度為0.904 nm 的GWCS 上最終運(yùn)動(dòng)狀態(tài)Fig. 10. Influence of groove depth on both sides of wedge-shaped copper structure on the droplet self-driving: (a) Velocity and (b) displacement of droplet versus time; (c) final motion state of 8 nm diameter droplet on GWCS with 0.904 nm groove depth.

4 結(jié) 論

本文采用分子動(dòng)力學(xué)方法研究了液滴在GWCS上的自驅(qū)動(dòng)行為, 計(jì)算了楔形基底角度和液滴尺寸對(duì)液滴自驅(qū)動(dòng)的影響. 結(jié)果表明, 液滴受基底表面潤(rùn)濕梯度和Young-Laplace 壓力梯度驅(qū)動(dòng), 實(shí)現(xiàn)了從楔形基底尖端向?qū)挾说某焖俣ㄏ蜃则?qū)動(dòng), 最大速度可達(dá)73.8 m/s, 并給出了液滴在楔形基底上的最大位移與基底潤(rùn)濕性、基底角度和液滴表面張力之間的理論關(guān)系; 能量變化規(guī)律表明, 液滴的勢(shì)能與流-固界面相互作用能存在競(jìng)爭(zhēng)關(guān)系, 即液滴與基底之間的相互作用能有部分轉(zhuǎn)化為液滴的勢(shì)能,并非完全作用于液滴驅(qū)動(dòng); 在微納尺度下, 楔形基底寬度為非連續(xù)線性增加, 寬度增加的軸向間隔距離隨基底角度的減小而增大, 導(dǎo)致小角度基底的潤(rùn)濕梯度連續(xù)性較差, 小角度基底上的液滴速度波動(dòng)相對(duì)較大, 且液滴在小角度基底上的位移小于理論值; 在大角度基底上, 小液滴更易獲得較大的驅(qū)動(dòng)速度, 同時(shí)由于液滴與基底之間接觸面積較小, 小液滴在無潤(rùn)濕梯度階段會(huì)產(chǎn)生較長(zhǎng)距離的減速運(yùn)動(dòng), 導(dǎo)致其實(shí)際位移要大于理論值; 當(dāng)楔形銅結(jié)構(gòu)兩側(cè)凹槽深度較大, 無法將單層石墨烯吸附進(jìn)凹槽時(shí), 液滴擁有最大的自驅(qū)動(dòng)速度和最遠(yuǎn)的輸運(yùn)距離.