• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    納米結(jié)構(gòu)及浸潤性對液滴潤濕行為的影響*

    2021-07-01 09:42:36李文馬驍婧徐進(jìn)良2王艷雷俊鵬
    物理學(xué)報(bào) 2021年12期
    關(guān)鍵詞:潤濕浸潤性液滴

    李文 馬驍婧? 徐進(jìn)良2) 王艷 雷俊鵬

    1) (華北電力大學(xué), 低品位能源多相流與傳熱北京市重點(diǎn)實(shí)驗(yàn)室, 北京 102206)

    2) (華北電力大學(xué), 電站能量傳遞轉(zhuǎn)化與系統(tǒng)教育部重點(diǎn)實(shí)驗(yàn)室, 北京 102206)

    1 引 言

    液滴在固體表面的潤濕行為是自然界中普遍存在的物理現(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軟件.

    2 模型建立及模擬細(xì)節(jié)

    2.1 模型構(gòu)建

    本文采用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.

    2.2 模擬體系控制方程

    不同種原子間的相互作用采用修正后的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.

    2.3 接觸角的標(biāo)定

    圖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.

    3 結(jié)果與討論

    3.1 表面浸潤性及納米結(jié)構(gòu)傾角對接觸角的影響

    納米液滴在潤濕過程中, 表面張力和分子間作用力是液滴變形與移動的驅(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 理論計(jì)算與MD模擬比較

    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ù)目隨φ的增加而增加.

    4 結(jié) 論

    采用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í),可以容納更多的流體原子.

    猜你喜歡
    潤濕浸潤性液滴
    基于低場核磁共振表征的礦物孔隙潤濕規(guī)律
    液滴間相互碰撞融合與破碎的實(shí)驗(yàn)研究
    噴淋液滴在空氣環(huán)境下的運(yùn)動特性
    浸潤性乳腺癌超聲及造影表現(xiàn)與P63及Calponin的相關(guān)性
    乙醇潤濕對2種全酸蝕粘接劑粘接性能的影響
    預(yù)潤濕對管道潤濕性的影響
    乳腺浸潤性微乳頭狀癌的研究進(jìn)展
    乳腺浸潤性導(dǎo)管癌組織β-catenin、cyclinD1、CDK4蛋白表達(dá)及臨床意義
    潤濕反轉(zhuǎn)劑的優(yōu)選及在氣體鉆井中的應(yīng)用
    氣井多液滴攜液理論模型研究
    河南科技(2013年10期)2013-08-12 05:55:10
    国产无遮挡羞羞视频在线观看| 成人特级黄色片久久久久久久| 999久久久国产精品视频| 亚洲片人在线观看| 叶爱在线成人免费视频播放| 国产精品成人在线| 午夜精品在线福利| 久久人妻熟女aⅴ| 国产精品自产拍在线观看55亚洲 | 多毛熟女@视频| 亚洲精品久久午夜乱码| 高潮久久久久久久久久久不卡| 一区二区三区国产精品乱码| 女性被躁到高潮视频| 动漫黄色视频在线观看| 久久久国产欧美日韩av| 亚洲成人国产一区在线观看| 成年人免费黄色播放视频| 久久久久久久午夜电影 | av天堂久久9| 超碰成人久久| av在线播放免费不卡| 性色av乱码一区二区三区2| x7x7x7水蜜桃| 久久精品亚洲精品国产色婷小说| 天天影视国产精品| 欧美在线一区亚洲| 美女扒开内裤让男人捅视频| 97人妻天天添夜夜摸| 精品国产一区二区三区四区第35| 欧美在线一区亚洲| 国产精品欧美亚洲77777| 97人妻天天添夜夜摸| 精品视频人人做人人爽| 亚洲av熟女| 大码成人一级视频| 国产精品99久久99久久久不卡| 日本欧美视频一区| 亚洲av欧美aⅴ国产| 亚洲情色 制服丝袜| 91在线观看av| 午夜影院日韩av| 大陆偷拍与自拍| 午夜老司机福利片| 亚洲av片天天在线观看| 久久久国产欧美日韩av| 精品福利永久在线观看| 亚洲欧洲精品一区二区精品久久久| 美国免费a级毛片| 国产精品久久久久久精品古装| 亚洲全国av大片| 欧美成人午夜精品| av中文乱码字幕在线| av视频免费观看在线观看| 亚洲av日韩精品久久久久久密| 亚洲aⅴ乱码一区二区在线播放 | 高清视频免费观看一区二区| 国产av精品麻豆| 视频在线观看一区二区三区| 国产极品粉嫩免费观看在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲中文日韩欧美视频| 嫩草影视91久久| 成人免费观看视频高清| 黄色片一级片一级黄色片| 少妇猛男粗大的猛烈进出视频| 欧美 日韩 精品 国产| 老汉色av国产亚洲站长工具| 性少妇av在线| 亚洲国产精品合色在线| 99re在线观看精品视频| 欧美精品人与动牲交sv欧美| 国内久久婷婷六月综合欲色啪| 国产精品98久久久久久宅男小说| 欧美黄色片欧美黄色片| 18禁黄网站禁片午夜丰满| 真人做人爱边吃奶动态| 久久人人97超碰香蕉20202| 女性生殖器流出的白浆| 国产精品亚洲av一区麻豆| 中文字幕精品免费在线观看视频| 久久精品aⅴ一区二区三区四区| 精品福利永久在线观看| 亚洲欧美日韩高清在线视频| 国产精品影院久久| 夜夜爽天天搞| 91大片在线观看| 国产亚洲精品久久久久久毛片 | 啦啦啦在线免费观看视频4| 三上悠亚av全集在线观看| 五月开心婷婷网| 三级毛片av免费| 香蕉久久夜色| 亚洲视频免费观看视频| 亚洲人成电影免费在线| 99国产精品一区二区三区| avwww免费| 一级,二级,三级黄色视频| 黄色丝袜av网址大全| 日本黄色视频三级网站网址 | 高清毛片免费观看视频网站 | 亚洲精品粉嫩美女一区| а√天堂www在线а√下载 | 精品高清国产在线一区| 一级黄色大片毛片| 啦啦啦免费观看视频1| 亚洲中文av在线| 母亲3免费完整高清在线观看| 18禁美女被吸乳视频| 国产精品国产av在线观看| 国产亚洲精品久久久久5区| 操出白浆在线播放| 一级毛片女人18水好多| 国产精品亚洲一级av第二区| 午夜免费成人在线视频| 久久香蕉国产精品| 最近最新免费中文字幕在线| 成熟少妇高潮喷水视频| 美女福利国产在线| 婷婷成人精品国产| 91精品三级在线观看| 黄片大片在线免费观看| svipshipincom国产片| 成人永久免费在线观看视频| 热re99久久精品国产66热6| 欧美午夜高清在线| 久久亚洲真实| 国产精品久久电影中文字幕 | 夫妻午夜视频| 久久青草综合色| 女性被躁到高潮视频| 欧美+亚洲+日韩+国产| 精品久久久久久,| 亚洲人成77777在线视频| 后天国语完整版免费观看| 日韩一卡2卡3卡4卡2021年| 高清在线国产一区| 高清欧美精品videossex| 欧美日韩亚洲高清精品| 欧美午夜高清在线| 亚洲av熟女| aaaaa片日本免费| 精品一区二区三区av网在线观看| 欧美精品啪啪一区二区三区| 亚洲伊人色综图| 99国产精品免费福利视频| 在线观看午夜福利视频| 亚洲第一欧美日韩一区二区三区| 午夜两性在线视频| 在线天堂中文资源库| 日本撒尿小便嘘嘘汇集6| 免费日韩欧美在线观看| 亚洲精品美女久久久久99蜜臀| 无人区码免费观看不卡| 变态另类成人亚洲欧美熟女 | 久久精品熟女亚洲av麻豆精品| 亚洲精品自拍成人| 91在线观看av| 国产成人精品久久二区二区免费| 午夜影院日韩av| 亚洲五月天丁香| 人人妻人人添人人爽欧美一区卜| 久久精品国产a三级三级三级| 中文字幕另类日韩欧美亚洲嫩草| 国产日韩一区二区三区精品不卡| videos熟女内射| 国产成人精品久久二区二区91| av天堂久久9| 韩国精品一区二区三区| 亚洲人成77777在线视频| 欧美人与性动交α欧美精品济南到| 麻豆国产av国片精品| 中文字幕最新亚洲高清| 亚洲欧美激情在线| 国产在线观看jvid| 亚洲中文av在线| 久久精品亚洲av国产电影网| 欧美国产精品一级二级三级| 欧美日韩精品网址| 国产又色又爽无遮挡免费看| 精品国产乱码久久久久久男人| 久99久视频精品免费| 免费观看精品视频网站| 热re99久久精品国产66热6| 国产欧美日韩一区二区三| 最近最新免费中文字幕在线| 国产激情久久老熟女| a级毛片黄视频| 久久亚洲精品不卡| 母亲3免费完整高清在线观看| av网站在线播放免费| 妹子高潮喷水视频| 国产伦人伦偷精品视频| 国产激情欧美一区二区| 久久精品熟女亚洲av麻豆精品| 如日韩欧美国产精品一区二区三区| 欧美日韩国产mv在线观看视频| 人妻 亚洲 视频| 亚洲熟女精品中文字幕| 激情在线观看视频在线高清 | 亚洲国产欧美网| 国产蜜桃级精品一区二区三区 | 久久国产精品影院| 91成人精品电影| 一二三四社区在线视频社区8| 精品免费久久久久久久清纯 | 日本黄色视频三级网站网址 | 国产精品永久免费网站| 国产精品免费视频内射| 国产高清激情床上av| 久久99一区二区三区| 欧美一级毛片孕妇| av不卡在线播放| 91成人精品电影| 久久人妻熟女aⅴ| 婷婷丁香在线五月| 欧美日韩av久久| 国产免费av片在线观看野外av| 亚洲一区二区三区欧美精品| 国产片内射在线| 亚洲精品在线美女| 丝袜美腿诱惑在线| 国产精品电影一区二区三区 | 久久热在线av| 国产高清激情床上av| 在线观看日韩欧美| 亚洲精品一卡2卡三卡4卡5卡| 五月开心婷婷网| 黑人猛操日本美女一级片| 天天躁夜夜躁狠狠躁躁| 成年人黄色毛片网站| 热99国产精品久久久久久7| 天堂√8在线中文| 成熟少妇高潮喷水视频| 亚洲成人免费av在线播放| 国产色视频综合| 一本大道久久a久久精品| 国产欧美日韩精品亚洲av| 成熟少妇高潮喷水视频| 亚洲成人免费av在线播放| 如日韩欧美国产精品一区二区三区| 又黄又粗又硬又大视频| 久久久久久久午夜电影 | 精品少妇一区二区三区视频日本电影| 99精国产麻豆久久婷婷| 久久国产亚洲av麻豆专区| 另类亚洲欧美激情| 99热网站在线观看| 伊人久久大香线蕉亚洲五| 最近最新免费中文字幕在线| 女性被躁到高潮视频| 久久久久久久午夜电影 | 亚洲欧美激情在线| 日韩人妻精品一区2区三区| 亚洲国产欧美网| 欧美黑人欧美精品刺激| 国产片内射在线| 天天躁狠狠躁夜夜躁狠狠躁| 三级毛片av免费| 亚洲中文字幕日韩| 久久天躁狠狠躁夜夜2o2o| 91在线观看av| 日韩熟女老妇一区二区性免费视频| 亚洲国产中文字幕在线视频| 婷婷精品国产亚洲av在线 | 久久精品国产a三级三级三级| 这个男人来自地球电影免费观看| 亚洲专区字幕在线| 老司机影院毛片| 精品第一国产精品| 国产欧美日韩一区二区三| 久久精品国产99精品国产亚洲性色 | 一区二区三区国产精品乱码| 国产欧美日韩精品亚洲av| 欧美国产精品一级二级三级| 久久国产乱子伦精品免费另类| 午夜精品在线福利| 建设人人有责人人尽责人人享有的| 中文字幕另类日韩欧美亚洲嫩草| 色播在线永久视频| 人人妻人人爽人人添夜夜欢视频| 麻豆国产av国片精品| 亚洲精品av麻豆狂野| 高清欧美精品videossex| 免费黄频网站在线观看国产| 亚洲人成77777在线视频| 黄色片一级片一级黄色片| 91在线观看av| 久久国产精品男人的天堂亚洲| 国产成人精品无人区| av一本久久久久| 飞空精品影院首页| 丝袜美腿诱惑在线| 法律面前人人平等表现在哪些方面| 麻豆乱淫一区二区| 99精品欧美一区二区三区四区| 精品熟女少妇八av免费久了| 美国免费a级毛片| 国产深夜福利视频在线观看| 久久中文字幕人妻熟女| 精品久久蜜臀av无| 欧美成人午夜精品| 黄色视频不卡| x7x7x7水蜜桃| 91九色精品人成在线观看| 一区福利在线观看| 久久精品亚洲av国产电影网| 久久国产亚洲av麻豆专区| 精品久久久久久,| 男人的好看免费观看在线视频 | 欧美激情极品国产一区二区三区| 这个男人来自地球电影免费观看| 丁香六月欧美| 韩国av一区二区三区四区| 亚洲美女黄片视频| 老鸭窝网址在线观看| 人妻丰满熟妇av一区二区三区 | 人人妻人人澡人人爽人人夜夜| 精品久久蜜臀av无| 黑人巨大精品欧美一区二区mp4| 99re6热这里在线精品视频| 夜夜夜夜夜久久久久| 侵犯人妻中文字幕一二三四区| 国产精品二区激情视频| e午夜精品久久久久久久| 亚洲色图av天堂| 亚洲熟妇熟女久久| 两人在一起打扑克的视频| 国产精品98久久久久久宅男小说| 国产av一区二区精品久久| 叶爱在线成人免费视频播放| 777米奇影视久久| 50天的宝宝边吃奶边哭怎么回事| 亚洲欧美激情在线| 成人亚洲精品一区在线观看| av天堂久久9| 亚洲美女黄片视频| 成人手机av| 脱女人内裤的视频| 国产亚洲精品第一综合不卡| e午夜精品久久久久久久| 黑人巨大精品欧美一区二区蜜桃| 久久这里只有精品19| 欧美日韩亚洲高清精品| 18禁黄网站禁片午夜丰满| 一进一出抽搐gif免费好疼 | xxx96com| 久久久久精品人妻al黑| 啦啦啦视频在线资源免费观看| 99精品欧美一区二区三区四区| 国产人伦9x9x在线观看| 午夜福利,免费看| 成人18禁高潮啪啪吃奶动态图| 在线观看免费视频网站a站| av网站在线播放免费| 啦啦啦免费观看视频1| 可以免费在线观看a视频的电影网站| 女人高潮潮喷娇喘18禁视频| 亚洲精品乱久久久久久| 国产又爽黄色视频| 高清av免费在线| 99re在线观看精品视频| 夫妻午夜视频| 久久中文看片网| 久久影院123| 在线永久观看黄色视频| 两个人看的免费小视频| 欧美日本中文国产一区发布| 国产精品电影一区二区三区 | 成人永久免费在线观看视频| 亚洲七黄色美女视频| 欧美日韩瑟瑟在线播放| 天天操日日干夜夜撸| 手机成人av网站| 高清黄色对白视频在线免费看| 两性午夜刺激爽爽歪歪视频在线观看 | 精品国产美女av久久久久小说| 婷婷丁香在线五月| 久久久精品免费免费高清| 色尼玛亚洲综合影院| 国产精品国产av在线观看| av一本久久久久| 国产一卡二卡三卡精品| 久久人人爽av亚洲精品天堂| 精品人妻熟女毛片av久久网站| 丰满迷人的少妇在线观看| 高清视频免费观看一区二区| 啦啦啦视频在线资源免费观看| 亚洲精品国产精品久久久不卡| 久久青草综合色| 大香蕉久久网| 免费一级毛片在线播放高清视频 | 一进一出好大好爽视频| 久久人妻av系列| 日本精品一区二区三区蜜桃| 丁香六月欧美| 99国产综合亚洲精品| 午夜福利影视在线免费观看| 满18在线观看网站| 亚洲专区国产一区二区| 黄片大片在线免费观看| 国产精品综合久久久久久久免费 | 岛国毛片在线播放| 在线观看舔阴道视频| 精品国产一区二区三区久久久樱花| 欧美乱码精品一区二区三区| 51午夜福利影视在线观看| 久久精品亚洲熟妇少妇任你| 人成视频在线观看免费观看| 亚洲熟女毛片儿| 在线观看午夜福利视频| 国产成人欧美| 老司机午夜十八禁免费视频| 欧美在线黄色| 久久午夜综合久久蜜桃| 日韩三级视频一区二区三区| 欧美乱码精品一区二区三区| 天天操日日干夜夜撸| 又紧又爽又黄一区二区| 麻豆国产av国片精品| 好男人电影高清在线观看| e午夜精品久久久久久久| 搡老岳熟女国产| 多毛熟女@视频| 欧美激情 高清一区二区三区| 国产成人av激情在线播放| 老汉色av国产亚洲站长工具| 精品国产国语对白av| 99国产综合亚洲精品| 欧美精品啪啪一区二区三区| 国产精品久久视频播放| 999久久久精品免费观看国产| 99精品久久久久人妻精品| 自线自在国产av| 亚洲一区高清亚洲精品| 久久精品国产亚洲av香蕉五月 | 成人亚洲精品一区在线观看| 不卡av一区二区三区| 亚洲五月天丁香| 热re99久久精品国产66热6| 欧美日韩福利视频一区二区| 色在线成人网| www.熟女人妻精品国产| 捣出白浆h1v1| 成人永久免费在线观看视频| 久久亚洲精品不卡| 久久ye,这里只有精品| 99精品欧美一区二区三区四区| 黄色毛片三级朝国网站| 久久人妻熟女aⅴ| 国产精品1区2区在线观看. | 变态另类成人亚洲欧美熟女 | 国产亚洲精品久久久久5区| 91九色精品人成在线观看| 免费高清在线观看日韩| 国产高清videossex| 韩国精品一区二区三区| 99热国产这里只有精品6| 岛国在线观看网站| 亚洲av熟女| 欧洲精品卡2卡3卡4卡5卡区| 一a级毛片在线观看| 首页视频小说图片口味搜索| 美女视频免费永久观看网站| 99精国产麻豆久久婷婷| 亚洲色图 男人天堂 中文字幕| 久久 成人 亚洲| 国产成人精品在线电影| 亚洲一区二区三区不卡视频| 中文字幕另类日韩欧美亚洲嫩草| 99国产精品99久久久久| 19禁男女啪啪无遮挡网站| 啦啦啦视频在线资源免费观看| 捣出白浆h1v1| 中文字幕av电影在线播放| 亚洲成av片中文字幕在线观看| 夫妻午夜视频| 黄色女人牲交| 久久国产精品人妻蜜桃| 老汉色av国产亚洲站长工具| 脱女人内裤的视频| 亚洲精品美女久久久久99蜜臀| 怎么达到女性高潮| 亚洲精品美女久久av网站| 一级毛片女人18水好多| 80岁老熟妇乱子伦牲交| 成人18禁高潮啪啪吃奶动态图| 午夜免费观看网址| www.熟女人妻精品国产| 少妇被粗大的猛进出69影院| av网站在线播放免费| 国产不卡一卡二| 久久久国产一区二区| 99re在线观看精品视频| 久久久国产一区二区| 亚洲国产看品久久| 欧美日韩亚洲国产一区二区在线观看 | 咕卡用的链子| 亚洲人成77777在线视频| 王馨瑶露胸无遮挡在线观看| 岛国毛片在线播放| 啦啦啦在线免费观看视频4| 在线天堂中文资源库| 国产99久久九九免费精品| 亚洲精品美女久久久久99蜜臀| 国产成人av激情在线播放| 天堂动漫精品| 脱女人内裤的视频| 岛国毛片在线播放| 美女福利国产在线| 免费高清在线观看日韩| 韩国av一区二区三区四区| 久久中文字幕人妻熟女| 欧美乱色亚洲激情| 性色av乱码一区二区三区2| 一本综合久久免费| 天堂中文最新版在线下载| 女性被躁到高潮视频| 交换朋友夫妻互换小说| 国产精品98久久久久久宅男小说| av超薄肉色丝袜交足视频| 午夜激情av网站| 麻豆av在线久日| 成人特级黄色片久久久久久久| 久久午夜亚洲精品久久| 国产高清视频在线播放一区| 久久久精品国产亚洲av高清涩受| 日本五十路高清| 国产区一区二久久| 天天影视国产精品| 黄色丝袜av网址大全| 国产不卡av网站在线观看| 黄片播放在线免费| 99热只有精品国产| 午夜免费鲁丝| 人人妻人人澡人人爽人人夜夜| 91在线观看av| 亚洲av日韩在线播放| 精品国产国语对白av| 精品国产乱码久久久久久男人| 欧美人与性动交α欧美软件| а√天堂www在线а√下载 | 日韩视频一区二区在线观看| 亚洲专区字幕在线| 满18在线观看网站| 国产精品98久久久久久宅男小说| 1024视频免费在线观看| 在线观看www视频免费| 99精品久久久久人妻精品| 嫩草影视91久久| 最新在线观看一区二区三区| 日本欧美视频一区| 国产成人啪精品午夜网站| 欧美午夜高清在线| 亚洲欧美日韩高清在线视频| 女人久久www免费人成看片| 成人精品一区二区免费| 香蕉丝袜av| 亚洲国产精品合色在线| 国产精品久久视频播放| 动漫黄色视频在线观看| 在线观看免费午夜福利视频| 乱人伦中国视频| 91麻豆av在线| 在线永久观看黄色视频| 淫妇啪啪啪对白视频| 成年人免费黄色播放视频| 在线观看舔阴道视频| 99国产精品一区二区蜜桃av | av网站免费在线观看视频| 久久香蕉精品热| 男人操女人黄网站| 久久婷婷成人综合色麻豆| 午夜老司机福利片| 欧美av亚洲av综合av国产av| 熟女少妇亚洲综合色aaa.| 久久天堂一区二区三区四区| 99精品在免费线老司机午夜| 国产精品自产拍在线观看55亚洲 | 99精国产麻豆久久婷婷| 又黄又爽又免费观看的视频| 中亚洲国语对白在线视频| 亚洲欧美日韩另类电影网站| 黄片大片在线免费观看| 999久久久国产精品视频| 亚洲成a人片在线一区二区| 99热只有精品国产| 又黄又爽又免费观看的视频| 国产成人一区二区三区免费视频网站| 两个人免费观看高清视频| 女人爽到高潮嗷嗷叫在线视频| 欧美成狂野欧美在线观看| avwww免费| 咕卡用的链子| 一区在线观看完整版| 国产av一区二区精品久久| 成人永久免费在线观看视频| 中文字幕人妻熟女乱码| 欧美午夜高清在线| 久久精品亚洲熟妇少妇任你| 国产成人av教育| 亚洲男人天堂网一区| 曰老女人黄片| 夜夜夜夜夜久久久久| 亚洲熟妇熟女久久| 麻豆成人av在线观看| 精品久久蜜臀av无| 午夜福利免费观看在线| 精品一区二区三区视频在线观看免费 | av网站免费在线观看视频| 国产精品影院久久| 黄色怎么调成土黄色| 国产精品久久视频播放| 欧美精品高潮呻吟av久久| 一区二区三区国产精品乱码|