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

    液滴碰撞親-疏水性組合壁面的數(shù)值分析

    2016-10-13 18:52:10陳圓圓徐進良李季巍
    化工學(xué)報 2016年12期
    關(guān)鍵詞:量綱液柱慣性力

    陳圓圓,徐進良,李季巍

    ?

    液滴碰撞親-疏水性組合壁面的數(shù)值分析

    陳圓圓,徐進良,李季巍

    (華北電力大學(xué)能源動力與機械工程學(xué)院,北京102206)

    采用VOF方法對液滴碰撞親-疏水性組合壁面過程進行數(shù)值模擬。通過量綱1化,在Reynolds數(shù)不變的情況下,研究Weber數(shù)和親水性區(qū)域直徑對液滴碰撞動態(tài)特性的影響,繪制了液滴碰撞結(jié)果分區(qū)圖。結(jié)果表明:隨著Weber數(shù)從1增加至100,液滴呈現(xiàn)出附著、短液滴一次破碎和長液滴多次破碎3種狀態(tài);液滴破碎過程受到表面張力、壁面黏附力及Rayleigh不穩(wěn)定性的影響,親水性區(qū)域直徑增加時,液滴破碎所對應(yīng)的臨界Weber數(shù)隨之增大;此外,Weber數(shù)較大時,液滴最大鋪展因子不受親水性區(qū)域直徑的影響,Weber數(shù)較小時,最大鋪展因子隨親水性區(qū)域直徑的增大而增大。

    液滴碰撞;親-疏水性;數(shù)值模擬;動態(tài)特性;鋪展因子

    引 言

    液滴碰撞現(xiàn)象廣泛存在于自然界、農(nóng)業(yè)及工業(yè)應(yīng)用中,如雨滴滴落、農(nóng)藥噴灑、噴墨打印等。液滴鋪展特性和碰撞結(jié)果與很多的實際應(yīng)用相關(guān),掌握其規(guī)律具有重要的意義,如在農(nóng)藥噴灑過程中,較大的鋪展面積可以提高農(nóng)藥利用率并節(jié)約成本;在法醫(yī)鑒定中,通過血液的鋪展或濺射情況可以對罪犯的位置進行推斷[1];在相變傳熱裝置中,液滴彈跳現(xiàn)象可用于先進換熱器的設(shè)計[2]。碰撞結(jié)果由液滴和壁面性質(zhì)共同決定,受慣性力、黏性力和表面張力控制,用量綱1參數(shù)Reynolds數(shù)和Weber數(shù)進行表征。?ikalo等[3]綜合研究了各因素對碰撞特性的影響。Marengo等[4]對液滴碰撞復(fù)雜壁面的特性進行了分析。李大樹等[5]、梁超等[6]研究了壁面浸潤性對碰撞結(jié)果及最大鋪展因子的影響。Palacios等[7]對液滴附著與飛濺的臨界條件進行了劃分。

    隨著材料制造技術(shù)的發(fā)展,人們已經(jīng)不再滿足于單一的均勻浸潤性壁面。當(dāng)壁面具有梯度浸潤性時,液滴會自發(fā)地由疏水性區(qū)域向親水性區(qū)域移動[8]。親-疏水性材料的組合產(chǎn)生了一些新的應(yīng)用,對于微流體和生物技術(shù)的研究具有重要意義,如控制壁面上濕潤和非濕潤區(qū)域,對液滴進行幾何塑形[9];平板印刷中,控制油性墨水和水性潤版液的潤濕和擴散,提高打印的準(zhǔn)確性[10];生物學(xué)上,利用生物活性分子在親水性與疏水性區(qū)域附著能力的差異,選擇性地控制蛋白質(zhì)、細胞及細菌的黏附[11]。彭本利等[12]利用親-疏水性間隔分布的設(shè)計強化了蒸汽冷凝傳熱。Mock等[13]使液滴偏心撞擊疏水性壁面上的圓形親水性區(qū)域,發(fā)現(xiàn)當(dāng)偏心距較小時,液滴會自動向親水性區(qū)域移動,最終完全位于親水性區(qū)域。Kim等[14]在超疏水性壁面上設(shè)計了一個環(huán)形超親水性區(qū)域,觀察到液滴碰撞壁面后中心區(qū)域飛濺,其余部分在親水性區(qū)域內(nèi)形成水環(huán)。Schutzius等[15]設(shè)計了雙圓環(huán)親水性結(jié)構(gòu),通過調(diào)整兩個圓環(huán)的寬度,使液滴中心形成一個軸對稱的“氣穴”,形成水環(huán)的同時又不會造成中心液滴的飛濺。當(dāng)親水性表面上具有一個條形的疏水性區(qū)域時,液滴可能被分割成兩部分。Son等[16]采用格子-玻耳茲曼的方法對親-疏水性區(qū)域接觸角的大小進行了研究,發(fā)現(xiàn)當(dāng)1?2≥80°時,液滴被分割,且差值越大,分割需要的時間越短。Song等[17]對碰撞速度和條形區(qū)域?qū)挾冗M行了研究,發(fā)現(xiàn)當(dāng)≥20時,液滴被分割,且分割時間約為12.2 ms,與碰撞速度無關(guān)。

    目前,大多數(shù)研究主要關(guān)注于親-疏水性組合壁面的制造技術(shù),而對其特性的研究依然有限。為了進一步掌握液滴碰撞組合壁面的規(guī)律,采用VOF方法,對其進行了模擬。壁面中心為圓形親水性區(qū)域,其余部分為超疏水性區(qū)域。將親水性區(qū)域直徑i與液滴初始直徑0之比定義為親水性區(qū)域相對直徑i,采用量綱1化計算,在壁面浸潤性組合和Reynolds數(shù)()不變的情況下,系統(tǒng)地研究了Weber數(shù)()和親水性區(qū)域相對直徑i對碰撞結(jié)果的影響。

    1 模型及方法

    1.1 控制方程及量綱1化

    模擬液滴碰撞親-疏水性組合壁面的過程,采用VOF方法捕捉氣液界面,兩相界面通過每個網(wǎng)格內(nèi)各相所占的體積分數(shù)來追蹤,當(dāng)體積分數(shù)=1時,為氣體;=0時,為液體;0<<1時,為兩相分界線?;究刂品匠虨?/p>

    式中,為速度矢量;為時間;為壓力;為密度;為動力黏度;為體積分數(shù);為重力加速度;為表面張力系數(shù);為曲率半徑。

    計算網(wǎng)格內(nèi)的平均密度和動力黏度分別為

    表面張力模型是由Brackbill等[18]提出的連續(xù)表面力(CSF)模型,將VOF計算中附加的表面張力處理為動量方程的源項。

    其中

    =?.(7)

    式中,下角標(biāo)和g分別代表了液相和氣相;為壁面接觸角。

    量綱1參數(shù)、均采用碰撞初始時刻液相參數(shù)進行定義,以液滴初始直徑0為特征長度,液滴與壁面接觸時的速度0為特征速度,則=ρu00/1,=ρu200/。普通重力條件下,采用空氣-水為工質(zhì),g=1.225 kg·m-3,ρ=103kg·m-3,g=1.79×10-5Pa·s,u=10-3Pa·s,=0.073 N·m-1,=9.8 m·s-2。

    對基本控制方程進行量綱1化,令

    將式(9)代入式(1)~式(3)得到量綱1控制方程

    其中,計算網(wǎng)格內(nèi)的平均密度和動力黏度分別為

    表面張力的計算采用連續(xù)表面力模型

    1.2 量綱1參數(shù)、邊界條件及初始條件

    與有量綱計算相比,量綱1計算只需將物性參數(shù)、初始條件和邊界條件轉(zhuǎn)化成對應(yīng)的量綱1值即可,其他求解方式與有量綱方法相同。根據(jù)控制方程的量綱1化推導(dǎo)可知,量綱1條件下密度、動力黏度、表面張力系數(shù)分別為ρ*=1,g*=g/ρμ*=1/,g*=g/(×μ),*=1/。在已知和的條件下,即可得出上述物性參數(shù),以及對應(yīng)的0和0,代入式(9)可得對應(yīng)的量綱1壓力和量綱1重力。

    同時,由式(9)可知,量綱1條件下=/0,=/0,因此初始條件下液滴直徑0=0/0=1,碰撞速度0=0/0=1??紤]到液滴下降過程會帶動周圍氣體的流動,初始時刻使液滴具有一定的高度,在重力加速度的條件下,保證液滴與壁面接觸時速度為1。采用二維軸對稱模型,選取5×5的計算域,如圖1所示:左邊為軸對稱邊界條件;上邊和右邊為壓力入口邊界條件,大小為量綱1大氣壓力atm;下邊為無滑移壁面邊界條件,壁面分為兩個部分,0≤≤i為親水性區(qū)域,>i為超疏水性區(qū)域,采用靜態(tài)接觸角模型,分別設(shè)為1=60°和2=150°。

    利用GAMBIT對網(wǎng)格進行劃分。網(wǎng)格尺寸采用量綱1值,分為兩種:在液滴的主要活動區(qū)域加密,加密區(qū)域為2×2的正方形,網(wǎng)格邊長為0.01;遠離液滴的區(qū)域網(wǎng)格邊長為0.02,網(wǎng)格總數(shù)為92500個。對網(wǎng)格獨立性進行檢驗,結(jié)果表明采用該網(wǎng)格既能保證計算精度,又可節(jié)約計算資源。連續(xù)性方程和動量方程采用二階隱式格式求解,對時間一階離散。壓力速度的耦合采用PISO算法,壓力求解采用Body Force Weighted,單元液相體積分數(shù)的計算采用Geo-Reconstruct方法。計算的時間步長為10-3,每個時間步長內(nèi)迭代次數(shù)為30次,收斂標(biāo)準(zhǔn)為連續(xù)性方程及動量方程殘差低于10-5。

    2 計算模型驗證

    為了驗證數(shù)值計算的合理性,分別針對文獻[19]中實驗4和實驗7進行了對比,分別對應(yīng)于液滴碰撞親水性壁面和疏水性壁面。以鋪展因子=/0隨量綱1時間的變化為驗證依據(jù),對比結(jié)果如圖2所示。雖然存在一定的誤差,但在一定程度上可以較好地模擬液滴碰撞過程。

    3 結(jié)果與分析

    3.1 不同浸潤性對液滴碰撞結(jié)果的影響

    液滴在=1000,=100的條件下,分別碰撞親水性(1=60°)、超疏水性(2=150°)和親水性-超疏水性組合壁面。其中,組合壁面上親水性區(qū)域直徑與液滴初始直徑相等,i=1。圖3、圖4分別是液滴碰撞3種壁面后鋪展因子及液滴形態(tài)隨時間的變化。根據(jù)圖3將碰撞過程分為3個階段:碰撞初期(<0.2),處于運動階段,由慣性力控制,不受壁面浸潤性影響[20],因此3條曲線完全重合,鋪展因子隨時間按指數(shù)規(guī)律變化,=3.08e0.61;之后浸潤性表現(xiàn)出影響,液滴在親水性壁面鋪展速度大于超疏水性壁面,在組合壁面上液滴已經(jīng)鋪展至超疏水性區(qū)域,因此與超疏水性壁面上鋪展因子繼續(xù)重合;當(dāng)液滴回縮至=1時,二者表現(xiàn)出差異,在超疏水性壁面上液滴繼續(xù)回縮直至脫離壁面(=0),在組合壁面上,液滴回縮至親-疏水性分界線,接觸線被親水性區(qū)域捕獲,不再回縮。

    由圖4可以看出,液滴形態(tài)變化與鋪展因子的變化相對應(yīng)。在親水性壁面上液滴快速鋪展形成薄液膜,達到最大鋪展后回縮,高度略微增加,但始終附著在壁面上。在超疏水性壁面上,液滴鋪展程度較小,因此與壁面摩擦所損耗的能量減少,具有足夠多的剩余動能可使液滴發(fā)生反彈,由于慣性力遠大于表面張力,液滴形成長液柱,完全脫離壁面,同時頂部形成二次液滴。液滴在組合壁面上形態(tài)類似于超疏水性壁面,同樣向上運動形成長液柱,并在頂部形成二次液滴。但不同的是,受底部親水性區(qū)域的影響,>5之后液滴接觸線不再回縮,上部繼續(xù)向上運動,在慣性力和壁面黏附力的反向作用下,液滴被進一步拉長,且壁面附近越來越細,最終斷裂,形成多個小液滴。

    3.2對液滴碰撞親-疏水性組合壁面的影響

    代表液滴慣性力與表面張力的比值,在不同下液滴會呈現(xiàn)出不同的現(xiàn)象。在低下,液滴慣性力≈表面張力,液滴鋪展程度較小,回縮過程獲得的動能相對較少,因此變形程度較低;高下,液滴慣性力?表面張力,液滴鋪展程度較大,有更多的表面能可轉(zhuǎn)化為動能,因此回縮過程可獲得足夠的動能向上運動。從Richard等[21]的研究可以看出,當(dāng)水滴碰撞超疏水性壁面時,根據(jù)不同,反彈過程會呈現(xiàn)兩種狀態(tài):≈1時,液滴近似為橢球形;≈18時,液滴被拉長形成液柱,并產(chǎn)出二次液滴,其形態(tài)與圖4(b)相似。Tsai等[22]分別針對低和高下水滴碰撞超疏水性壁面進行了研究,同樣可以看出,在低下(=1~10),液滴變形程度較小,近似呈橢球形。在高下(=120),液滴被拉長為液柱,并破碎產(chǎn)生多個小液滴,同時底部仍然殘留部分液體,與圖4(c)相似。

    由此可見,在均勻超疏水性壁面上對其形態(tài)具有重要影響。接下來將研究親-疏水性組合壁面上的影響。圖5是i=0.5時,分別對應(yīng)于=1、10、100液滴形態(tài)隨時間的變化。當(dāng)=1時,液滴變形程度較小,在經(jīng)歷了鋪展過程后略微回縮,并未出現(xiàn)反彈和破碎的情況。其原因有兩點:①=1,慣性力≈表面張力,液滴運動的慣性力被表面張力克服,在回縮過程中沒有足夠的動力向上運動。②Bond數(shù)()代表重力與表面張力之比,等于ρ*2/*,對應(yīng)于=1、10、100的3種工況,分別是6.3、6.3×10-2、6.3×10-4。對于后兩種工況<<1,重力可以忽略不計,而=1時,重力大于表面張力,大于慣性力,起主導(dǎo)作用,使液滴在垂直方向受力增加,抑制了液滴的向上運動。當(dāng)=10時,由文獻[22]可知,在均勻超疏水性壁面上,液滴呈近橢圓形反彈。觀察圖5(b)可以發(fā)現(xiàn),液滴變形程度介于圖5(a)、(c)之間,形成短液柱。同時由于底部親水性區(qū)域的存在,液滴未完全反彈,而是破碎成了兩部分,底部殘留在壁面上。當(dāng)=100時,慣性力?表面張力,液滴受慣性力主導(dǎo),使得表面張力幾乎可以忽略不計。由圖5(c)可以看出,在鋪展階段,液滴鋪展直徑遠大于圖5(a)、(b),使得液滴在回縮過程中,有更多的表面能轉(zhuǎn)化為動能。液滴回縮后,液滴在慣性力的作用下向上運動,同時底部受壁面黏附力的作用無法脫離,使得液滴逐漸被拉長成為液柱,并且破碎產(chǎn)生多個小液滴。

    Fig .5 Influence of Weber number on outcomes of drop impact dynamic on hybrid surface withi=0.5

    圖5(b)、(c)中兩種破碎雖然都是表面張力作用的結(jié)果,但具體原因卻不相同。短液滴的斷裂位置離壁面較近,受壁面影響較大。圖6是液滴一次破碎過程中壓強和速度的分布,在=2時刻,液滴頂端在表面張力下收縮,速度轉(zhuǎn)向,頂端受到擠壓壓力增大為液滴內(nèi)部最大值。隨著液滴不斷向上運動,球形頂部半徑增大,壓強減小,同時“頸部”壓強增大,在=4時達到最大值。在壓強差下,“頸部”附近液體加速向上運動。同時,壁面附近的液滴受到親水性區(qū)域的黏附力,在=3時刻速度方向發(fā)生改變,向下運動,導(dǎo)致液體快速運動至液滴頂部的同時,“頸部”以下得不到及時的補充,因此頂部區(qū)域半徑進一步增加,而“頸部”越來越細,最終在表面張力的作用下發(fā)生斷裂。

    對于長液柱來說,頂端收縮成球,產(chǎn)生二次液滴的機理與短液滴一致。但是由于液柱較長,產(chǎn)生二次液滴后會繼續(xù)收縮產(chǎn)生新的頂端,再次脫離。液柱向上運動的同時產(chǎn)生毛細波向下傳遞,并相互干涉,使液滴受到Rayleigh不穩(wěn)定性的影響[23]。液柱底部受到親水性壁面黏附力的作用,速度轉(zhuǎn)向,在壁面附近產(chǎn)生凹陷,加劇了Rayleigh不穩(wěn)定,最終導(dǎo)致液柱破碎,產(chǎn)生多個小液滴。

    3.3 親水性區(qū)域相對直徑i對液滴碰撞組合壁面的影響

    圖5介紹了i=0.5時,不同下液滴的形態(tài),圖7~圖9將對i增大時,對應(yīng)工況的變化進行研究。圖7~圖9分別對應(yīng)于=1、10、100,針對每個工況選取4個i,分別為0.5、1.0、1.5、2.0。首先觀察圖7,在i從0.5增大至2.0的過程中,液滴形態(tài)變化不大,在經(jīng)歷了鋪展過程后都略微回縮,并未出現(xiàn)反彈和破碎的情況。由3.2節(jié)的分析可知,其原因在于該工況下較大,重力的作用超過了表面張力和慣性力的作用,使得浸潤性的影響減小并抑制了液滴的向上運動。當(dāng)增大時,慣性力的作用增強,使得液滴具有更多的動能向上運動,可以看出圖8中液滴在豎直方向的運動明顯高于圖7,且隨著i的增加,液滴經(jīng)歷了一次破碎到不破碎的過渡過程:當(dāng)i=0.5時,液滴在向上運動過程中受表面張力的作用,頂端收縮形成“頸部”,并在慣性力和壁面黏附力的反向拉扯下,逐漸變細,最終發(fā)生斷裂。當(dāng)i增加時,親水性區(qū)域增加,壁面黏附力的作用增強,消耗了慣性力的作用。從圖8(b)可以看出,在=3時液滴同樣收縮形成了“頸部”,但是底部與壁面接觸的面積增大,且在=4時液柱高度不再增加,說明液滴開始向下運動,不再對“頸部”產(chǎn)生拉扯作用,使得“頸部”直徑再次增大,而沒能發(fā)生斷裂。當(dāng)i繼續(xù)增加時,液滴向上運動程度明顯減弱,甚至不再產(chǎn)生明顯的“頸部”,因此附著在壁面上不再斷裂。當(dāng)=100時,慣性力遠遠大于表面張力,處于主導(dǎo)地位,液滴鋪展程度較大,且回縮過程具有足夠的剩余動能使液滴發(fā)生反彈。由于底部受親水性區(qū)域的黏附力無法脫離壁面,隨著液滴不斷向上運動,液滴被拉長成為液柱,如圖9所示。當(dāng)i=0.5時,親水性區(qū)域較小,液柱較為均勻,在表面張力、慣性力和壁面黏附力的作用下斷裂形成小液珠。當(dāng)i=1.0時,形態(tài)與i=0.5相似,但是底部與壁面接觸的面積增加,斷裂后底部殘留的液滴體積增加。i=1.5時,液滴底部體積進一步增大,使得液柱上邊部分變細,產(chǎn)生的小液滴體積減小,數(shù)量增加。尤其是當(dāng)i=2.0時,液滴形態(tài)與i=0.5時已大不相同,大部分液滴都位于壁面附近,=8開始頂部液柱已經(jīng)破碎成多個小液珠,而底部向下收縮,逐漸穩(wěn)定在親水性區(qū)域。

    3.4 不同和i下的碰撞結(jié)果

    在固定的前提下,液滴是否破碎主要受和親水性區(qū)域相對直徑i影響。決定了液滴回縮后向上運動的程度,i決定了液滴向上運動過程中與壁面接觸面積的大小。當(dāng)較大而i較小時,液滴向上運動趨勢明顯,與壁面接觸面積較小,在慣性力和壁面黏附力的相反作用下,容易被拉長發(fā)生破碎;當(dāng)較小而i較大時,液滴向上運動趨勢減弱,與壁面接觸面積增大,慣性力和壁面黏附力之間的反向作用減弱,抑制了破碎的發(fā)生。為了進一步掌握液滴破碎所對應(yīng)的臨界條件,進行了多組計算,結(jié)果如圖10所示:對于不同的i,隨著的增大,液滴均經(jīng)歷了附著、一次破碎到多次破碎的過渡;當(dāng)i=0.5時,親水性區(qū)域直徑小于液滴初始直徑,極易發(fā)生破碎,i增大時,破碎對應(yīng)的臨界隨之增大。當(dāng)i=0.5、1.0、1.5時,一次破碎與多次破碎的分界線相同,臨界小于i=2.0的情況。

    ▲ for multi-time pinch-off; ■ for one-time pinch-off; ● for non-breakup

    3.5 液滴最大鋪展因子max

    液滴碰撞均勻浸潤性壁面時,最大鋪展因子max隨著和的增大而增大。在和較小時,浸潤性對max的影響較大,隨著和的增加,浸潤性的影響逐漸減弱[24]。液滴碰撞親-疏水性組合壁面時,除了和之外,還受親水性區(qū)域相對直徑i的影響。圖11是=1000,=100、10、1的液滴碰撞組合壁面時,最大鋪展因子max隨i的變化。

    i=0時,壁面為均勻超疏水性。可以看出,當(dāng)=100時,max不受i的影響,與均勻超疏水性壁面上結(jié)果一致;當(dāng)=10、1,、0.5、1.0時,max依然不受i的影響,i=1.5、2.0時,max隨i的增加而增加。從總體來看,max并沒有隨著的增加而增加,=1所對應(yīng)的max遠遠大于=10的結(jié)果。Wu等[25]曾通過數(shù)值模擬研究過均勻浸潤性壁面上對液滴鋪展的影響,發(fā)現(xiàn)隨著的增加,重力作用增強,液滴覆蓋面積會隨著增大。該研究為本文觀察到的反?,F(xiàn)象提供了依據(jù):由3.2節(jié)可知,對應(yīng)于=1、10、100的3種工況,6.3、6.3×10-2、6.3×10-4。當(dāng)=1時,較大,重力作用大于表面張力,使液滴在垂直方向受力增加,影響了液滴的鋪展規(guī)律,出現(xiàn)了=1所對應(yīng)的max遠遠大于=10所對應(yīng)的max的反?,F(xiàn)象。

    4 結(jié) 論

    采用VOF方法模擬了液滴碰撞親-疏水性組合壁面的動態(tài)特性。通過量綱1化,研究了=1000的條件下,和親水性區(qū)域相對直徑i對碰撞結(jié)果的影響。具體結(jié)論如下:

    (1)=100時,液滴分別碰撞親水性、超疏水性和親-疏水性組合壁面,依次表現(xiàn)為附著、反彈和液柱破碎。

    (2)液滴碰撞親-疏水性組合壁面時,隨著從1增加至100,液滴分別表現(xiàn)為附著、短液柱一次破碎和長液柱多次破碎3種結(jié)果。破碎的主要原因是表面張力,短液柱的破碎受壁面黏附力影響,而長液柱的破碎受Rayleigh不穩(wěn)定性影響。

    (3)親水性區(qū)域相對直徑i=0.5時,液滴易發(fā)生破碎,對應(yīng)的臨界值為5;i增大時,破碎對應(yīng)的臨界隨之增大。

    (4)在=100時,最大鋪展因子max不受i的影響。較小時,若i>1.0,則max隨i的增大而增大。且=1時,對應(yīng)的較大,重力的影響大于表面張力,液滴在垂直方向受到擠壓,使得max大于=10的情況。

    符 號 說 明

    Bo——Bond數(shù) D——量綱1直徑 d——直徑,m F——量綱1表面張力 f——表面張力,N G——量綱1重力 g——重力,N K——量綱1曲率半徑 k——曲率半徑,m?1 n——壁面法向量 P——量綱1壓力 p——壓力,Pa R——量綱1半徑 Re——Reynolds數(shù) r——半徑,m t——時間,s U——垂直方向量綱1速度 u——垂直方向速度,m·s?1 V——量綱1速度矢量 v——速度矢量,m·s?1 We——Weber數(shù) Z——量綱1高度 z——高度,m β——鋪展因子 θ——接觸角,(°) μ——動力黏度,Pa·s μ*——量綱1動力黏度 ρ——密度,kg·m?3 ρ*——量綱1密度 σ——表面張力系數(shù),N·m?1 σ*——量綱1表面張力系數(shù) τ——量綱1時間 φ——氣相體積分數(shù) 下角標(biāo) g——氣相 i——親-疏水性分界線 l——液相 max——最大值 0——初始值 1——親水性壁面 2——超疏水性壁面

    [1] LAAN N, DE BRUIN K G, BARTOLO D,Maximum diameter of impacting liquid droplets [J]. Physical Review Applied, 2014, 2(4): 1-7.

    [2] BOREYKO J B, CHEN C H. Vapor chambers with jumping-drop liquid return from superhydrophobic condensers [J]. International Journal of Heat and Mass Transfer, 2013, 61(1): 409-418.

    [3] ?IKALO ?, GANIC E N. Phenomena of droplet -surface interactions [J]. Experimental ThermalFluid Science, 2006, 31(2): 97-110.

    [4] MARENGO M, ANTONINI C, ROISMAN I V,. Drop collisions with simple and complex surfaces [J]. Current Opinion in Colloid & Interface Science, 2011, 16(4): 292-302.

    [5] 李大樹, 仇性啟, 鄭志偉, 等. 液滴沖擊不同浸潤性壁面的數(shù)值分析[J]. 農(nóng)業(yè)機械學(xué)報, 2015, 46(7): 294-302. LI D S, QIU X Q, ZHENG Z W,Numerical analysis of droplet impact on surfaces with different wettabilities [J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(7): 294-302.

    [6] 梁超, 王宏, 朱恂, 等. 液滴撞擊不同浸潤性壁面動態(tài)過程的數(shù)值模擬[J]. 化工學(xué)報, 2013, 64(8): 2745-2751. LIANG C H, WANG H, ZHU X,. Numerical simulation of droplet impact on surfaces with different wettabilities [J]. CIESC Journal, 2013, 64(8): 2745-2751.

    [7] PALACIOS J, HERNANDEZ J, GOMEZ P,Experimental study of splashing patterns and the splashing/deposition threshold in drop impacts onto dry smooth solid surfaces [J]. Experimental Thermal & Fluid Science, 2013, 44(44): 571-582.

    [8] 王宏, 廖強, 朱恂. 梯度表面能材料上液滴運動機理[J]. 化工學(xué)報, 2007, 58(9): 2313-2320. WANG H, LIAO Q, ZHU X. Mechanism of liquid droplet movement on surface with gradient surface energy [J]. Journal of Chemical Industry and Engineering (China), 2007, 58(9): 2313-2320.

    [9] HANCOCK M J, YANAGAWA F, JANGY H,Designer hydrophilic regions regulate droplet shape for controlled surface patterning and 3D microgel synthesis [J]. Small, 2012, 8(3): 393-403.

    [10] 卞燁, 王卉. 超親/超疏水微圖案的制作及應(yīng)用[J]. 現(xiàn)代化工, 2014, 34(9): 12-16. BIAN Y, WANG H. Fabrication and applications of superhydrophilic/superhydrophobic micropatterns [J]. Modern Chemical Industry, 2014, 34(9): 12-16.

    [11] GALOPIN E, PIRET G, SZUNERITS S,Selective adhesion ofspores on heterogeneously wetted silicon nanowires [J]. Langmuir, 2010, 26(5): 3479-3484.

    [12] 彭本利, 馬學(xué)虎, 蘭忠, 等. 組合表面調(diào)控液滴特性強化蒸汽冷凝傳熱[J]. 化工學(xué)報, 2015, 66(10): 3826-3833. PENG B L, MA X H, LAN H,. Steam condensation heat transfer enhancement through droplet properties manipulation with hybrid surfaces [J]. CIESC Journal, 2015, 66(10): 3826-3833.

    [13] MOCK U, MICHEL T, TROPEA C,. Drop impact on chemically structured arrays [J]. Journal of Physics Condensed Matter, 2005, 17(9): S595-S605.

    [14] KIM S, MOON M W, KIM H Y. Drop impact on super-wettability-contrast annular patterns [J]. Journal of Fluid Mechanics, 2013, 730(5): 328-342.

    [15] SCHUTZIUS T M, GRAEBER G, ELSHARKAWY M,. Morphing and vectoring impacting droplets by means of wettability-engineered surfaces [J]. Scientific Reports, 2014, 4: 7029.

    [16] SON S W, CHAN H J, HA M Y,. A numerical study on the dynamic behavior of the liquid droplet located on heterogeneous surface [J]. Computers & Fluids, 2014, 105: 294-306.

    [17] SONG D, SONG B, HU H,. Selectively splitting a droplet using superhydrophobic stripes on hydrophilic surfaces [J]. Physical Chemistry Chemical Physics, 2015, 17(21): 13800-13803.

    [18] BRACKBILL J U, KOTHE D B, ZEMACH C. A continuum method for modeling surface tension [J]. Journal of Computational Physics, 1992, 100(2): 335-354.

    [19] ?IKALO ?, WILHELM H D, ROISMAN I V,. Dynamic contact angle of spreading droplets: experiments and simulations [J]. Physics of Fluids, 2005, 17(6): 1-13.

    [20] RIOBOO R, MARENGO M, TROPEA C. Time evolution of liquid drop impact onto solid, dry surfaces [J]. Experiments in Fluids, 2002, 33(1): 112-124.

    [21] RICHARD D, QUERE D. Bouncing water drops [J]. Europhysics Letters, 2000, 50(6): 769-775.

    [22] TSAI P, PACHECO S, PIRAT C,Drop impact upon micro- and nanostructured superhydrophobic surfaces [J]. Langmuir, 2009, 25(20): 12293-12298.

    [23] SHINJO J, UMERNMURA A. Simulation of liquid jet primary breakup: dynamics of ligament and droplet formation [J]. International Journal of Multiphase Flow, 2010, 36(7): 513-532.

    [24] ?IKALO ?, MARENGO M, TROPEA C,. Analysis of impact of droplets on horizontal surfaces [J]. Experimental Thermal & Fluid Science, 2002, 25(7): 503-510.

    [25] WU J, LI Y D. Dynamic performance of a static or throwing droplet impact onto a solid substrate with different properties [J]. Physica A Statistical Mechanics & Its Applications, 2016, 446: 158-170.

    Numerical simulation of droplet impact on hybrid surfaces with different wettabilities

    CHEN Yuanyuan, XU Jinliang, LI Jiwei

    (College of Energy Power and Mechanical Engineering, North China Electric Power University, Beijing 102206, China)

    Surfaces with patterned domains of different wettabilities can be used for controlling the impact process. The dynamic behavior of droplet impact on hybrid surfaces with different wettabilities is simulated using the VOF method by dimensionless. Effects of Weber number and hydrophilic region diameter are studied with constant Reynold number. A drop impact regime map is generated, in which the impact dynamics is characterized as a function of Weber number and relative diameter of hydrophilic region. With the increase of Weber number, the droplet presents three different status: deposition, pinch-off as a short liquid columnand break up as a long liquid column. The threshold Weber number where pinch-off occurs will increase with the increase of hydrophilic region diameter. For high Weber number, the maximum spreading factors are the same for different hydrophilic region diameter. While for low Weber number, the maximum spreading factors will increase with it.

    droplet impact; wettabilities; simulation; dynamic behavior; spreading factors

    date: 2016-05-12.

    Prof. XU Jinliang, xjl@ncepu.edu.cn

    10.11949/j.issn.0438-1157.20160662

    TK 124

    A

    0438—1157(2016)12—5006—09

    國家自然科學(xué)基金項目(51436004);國家自然科學(xué)基金國際合作與交流項目(51210011)。

    supported by the National Natural Science Foundation of China (51436004) and the National Natural Science Foundation of International Cooperation and Exchange Programs (51210011).

    2016-05-12收到初稿,2016-09-22收到修改稿。

    聯(lián)系人:徐進良。第一作者:陳圓圓(1991—),女,碩士研究生。

    猜你喜歡
    量綱液柱慣性力
    巧用“形象思維”,速解液柱(活塞)移動問題
    量綱分析在熱力學(xué)統(tǒng)計物理中的教學(xué)應(yīng)用*
    物理通報(2022年6期)2022-05-27 09:43:14
    新型模切機滑塊運動特性及其慣性力特性研究
    淺談量綱法推導(dǎo)物理公式的優(yōu)勢
    ——以勻加速直線運動公式為例
    培養(yǎng)科學(xué)思維 落實核心素養(yǎng)
    從“慣性力”到“洛希極限”
    科技論文編輯加工中的量綱問題
    豎直液柱與水平液面作用激起毛細波探究
    液柱在激波沖擊下RM不穩(wěn)定性和破裂過程的數(shù)值計算
    量綱分析及其應(yīng)用
    物理與工程(2012年6期)2012-07-05 05:31:48
    av在线播放精品| 六月丁香七月| 99九九在线精品视频 | 精品人妻熟女毛片av久久网站| 久热这里只有精品99| 少妇的逼水好多| 国产白丝娇喘喷水9色精品| 日本黄大片高清| 日韩av不卡免费在线播放| 久久99蜜桃精品久久| 日产精品乱码卡一卡2卡三| 亚洲av.av天堂| 国产亚洲午夜精品一区二区久久| 久久久久网色| 亚洲图色成人| 久久这里有精品视频免费| 亚洲成色77777| 国产真实伦视频高清在线观看| 精品一区二区免费观看| 黄色怎么调成土黄色| 国产片特级美女逼逼视频| 免费黄网站久久成人精品| 日韩欧美一区视频在线观看 | 欧美日韩视频高清一区二区三区二| 一区二区三区免费毛片| 国产乱来视频区| 伊人久久国产一区二区| 伦精品一区二区三区| www.av在线官网国产| 午夜日本视频在线| 亚洲精品国产色婷婷电影| 精品久久久久久久久av| 看非洲黑人一级黄片| 高清视频免费观看一区二区| 国产精品久久久久久久久免| 97精品久久久久久久久久精品| 免费黄色在线免费观看| 欧美精品亚洲一区二区| 国产精品一区www在线观看| 丝袜脚勾引网站| 欧美精品亚洲一区二区| 啦啦啦啦在线视频资源| 黑人巨大精品欧美一区二区蜜桃 | 日本爱情动作片www.在线观看| 热re99久久国产66热| 美女cb高潮喷水在线观看| 国产片特级美女逼逼视频| 美女视频免费永久观看网站| 国产精品国产av在线观看| 高清午夜精品一区二区三区| 99久久精品热视频| 人人澡人人妻人| av女优亚洲男人天堂| 蜜桃久久精品国产亚洲av| 日日摸夜夜添夜夜添av毛片| 日韩一本色道免费dvd| 国产成人精品婷婷| 午夜福利影视在线免费观看| 亚洲精品456在线播放app| 午夜福利网站1000一区二区三区| 搡女人真爽免费视频火全软件| 嫩草影院新地址| 久久久久精品性色| 日韩电影二区| 亚洲av欧美aⅴ国产| 国产极品粉嫩免费观看在线 | 国产精品女同一区二区软件| 亚洲国产欧美日韩在线播放 | 一区二区三区免费毛片| 国产在线一区二区三区精| 97超碰精品成人国产| 亚洲国产精品一区二区三区在线| 国产男女内射视频| 久久精品国产亚洲网站| 一级毛片aaaaaa免费看小| 成人亚洲欧美一区二区av| 亚洲精品,欧美精品| 亚洲内射少妇av| 亚洲va在线va天堂va国产| 少妇人妻精品综合一区二区| 国产色爽女视频免费观看| av天堂中文字幕网| 91aial.com中文字幕在线观看| 亚洲精品国产色婷婷电影| 99热6这里只有精品| av在线app专区| 国产在线视频一区二区| 日本午夜av视频| 欧美少妇被猛烈插入视频| 国产亚洲精品久久久com| 亚洲国产精品国产精品| 午夜福利在线观看免费完整高清在| av福利片在线| 91精品伊人久久大香线蕉| 免费看日本二区| 国产欧美日韩精品一区二区| 日韩欧美 国产精品| 在线观看三级黄色| 国产成人aa在线观看| 日韩一区二区视频免费看| 国产免费福利视频在线观看| 爱豆传媒免费全集在线观看| 国产 精品1| 成人综合一区亚洲| 免费观看在线日韩| 亚洲精品色激情综合| 99九九在线精品视频 | 菩萨蛮人人尽说江南好唐韦庄| 精品少妇久久久久久888优播| 建设人人有责人人尽责人人享有的| 国产精品国产av在线观看| 欧美亚洲 丝袜 人妻 在线| 精品久久久久久久久av| 亚洲欧美一区二区三区国产| 国产欧美日韩精品一区二区| av在线观看视频网站免费| 日本av手机在线免费观看| 亚洲欧美中文字幕日韩二区| 欧美高清成人免费视频www| 日日啪夜夜撸| 国产成人免费无遮挡视频| 国产深夜福利视频在线观看| av.在线天堂| 国产精品秋霞免费鲁丝片| 国语对白做爰xxxⅹ性视频网站| 国产一区有黄有色的免费视频| 两个人的视频大全免费| 国产日韩一区二区三区精品不卡 | 又黄又爽又刺激的免费视频.| 午夜福利网站1000一区二区三区| 校园人妻丝袜中文字幕| 免费黄网站久久成人精品| 黄色日韩在线| 亚洲国产欧美在线一区| 色视频www国产| 岛国毛片在线播放| 亚洲精品色激情综合| 简卡轻食公司| 少妇高潮的动态图| 国产精品三级大全| 三上悠亚av全集在线观看 | 亚洲精品日本国产第一区| 美女中出高潮动态图| 日韩人妻高清精品专区| 天堂8中文在线网| 久久99热6这里只有精品| 国产精品人妻久久久久久| 欧美日韩视频精品一区| 热99国产精品久久久久久7| 赤兔流量卡办理| 亚洲美女黄色视频免费看| 国产av国产精品国产| 亚洲精品乱码久久久v下载方式| 国产伦精品一区二区三区四那| 成年人免费黄色播放视频 | 久久99一区二区三区| 国产永久视频网站| 国产一区二区在线观看av| 丰满乱子伦码专区| 最近手机中文字幕大全| 特大巨黑吊av在线直播| 欧美人与善性xxx| 亚洲欧美精品自产自拍| 亚洲av成人精品一二三区| 成人漫画全彩无遮挡| 国产精品熟女久久久久浪| 在线亚洲精品国产二区图片欧美 | 黄色配什么色好看| 日韩av不卡免费在线播放| av免费在线看不卡| 性色avwww在线观看| 国产一级毛片在线| 欧美日韩一区二区视频在线观看视频在线| 9色porny在线观看| 久久久a久久爽久久v久久| 天天操日日干夜夜撸| 青春草亚洲视频在线观看| 国产91av在线免费观看| 日韩av不卡免费在线播放| 久久精品久久精品一区二区三区| 99热这里只有是精品50| 亚洲天堂av无毛| 国产在视频线精品| a级片在线免费高清观看视频| 久久精品国产自在天天线| 高清欧美精品videossex| 亚洲无线观看免费| 又大又黄又爽视频免费| 日韩一区二区三区影片| 九九在线视频观看精品| 在线观看国产h片| 18禁动态无遮挡网站| 欧美日韩一区二区视频在线观看视频在线| 少妇裸体淫交视频免费看高清| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品乱码久久久久久按摩| 日日摸夜夜添夜夜爱| 亚洲国产毛片av蜜桃av| 99热国产这里只有精品6| 最近手机中文字幕大全| 久久精品国产亚洲av涩爱| 我的女老师完整版在线观看| 另类精品久久| 亚洲av日韩在线播放| 免费大片18禁| 乱码一卡2卡4卡精品| 亚洲精品亚洲一区二区| 亚洲精品国产av成人精品| 蜜桃在线观看..| 久久久久精品久久久久真实原创| 欧美精品亚洲一区二区| 精品国产一区二区三区久久久樱花| 男女边吃奶边做爰视频| 熟妇人妻不卡中文字幕| 亚洲真实伦在线观看| 久久午夜福利片| 午夜免费观看性视频| 国产极品粉嫩免费观看在线 | 麻豆精品久久久久久蜜桃| 国产白丝娇喘喷水9色精品| 国产国拍精品亚洲av在线观看| 亚洲精品中文字幕在线视频 | 毛片一级片免费看久久久久| 尾随美女入室| 在线观看三级黄色| 乱人伦中国视频| 国产精品国产三级专区第一集| 国产免费一区二区三区四区乱码| 国产成人freesex在线| 啦啦啦在线观看免费高清www| 婷婷色av中文字幕| 永久免费av网站大全| 国产成人免费无遮挡视频| 两个人免费观看高清视频 | 免费少妇av软件| 美女主播在线视频| 国产精品99久久99久久久不卡 | 亚洲国产av新网站| 久久人人爽人人片av| 你懂的网址亚洲精品在线观看| 偷拍熟女少妇极品色| 欧美一级a爱片免费观看看| 国产精品偷伦视频观看了| 一本大道久久a久久精品| 精品少妇黑人巨大在线播放| 又大又黄又爽视频免费| 国产精品伦人一区二区| 亚洲一区二区三区欧美精品| 男人爽女人下面视频在线观看| 女人久久www免费人成看片| 亚洲欧洲精品一区二区精品久久久 | 亚洲内射少妇av| 高清毛片免费看| 免费人妻精品一区二区三区视频| 久久久久网色| 欧美xxⅹ黑人| 夜夜骑夜夜射夜夜干| 国产在线视频一区二区| 美女xxoo啪啪120秒动态图| 在线观看一区二区三区激情| 桃花免费在线播放| 熟妇人妻不卡中文字幕| 69精品国产乱码久久久| 性色av一级| 交换朋友夫妻互换小说| 亚洲四区av| 亚洲中文av在线| 中国国产av一级| 男人添女人高潮全过程视频| 一级毛片 在线播放| 中文字幕亚洲精品专区| 亚洲精品日韩av片在线观看| 国产视频内射| 91午夜精品亚洲一区二区三区| 插阴视频在线观看视频| 久久久久久久久久人人人人人人| 日本爱情动作片www.在线观看| 80岁老熟妇乱子伦牲交| xxx大片免费视频| 成年av动漫网址| 国产精品偷伦视频观看了| 久久久久精品性色| 在线天堂最新版资源| 免费观看无遮挡的男女| 亚洲欧美日韩东京热| 噜噜噜噜噜久久久久久91| 国产免费又黄又爽又色| 色吧在线观看| 极品少妇高潮喷水抽搐| 99久久精品热视频| 黑人巨大精品欧美一区二区蜜桃 | 亚洲精品自拍成人| 黄色视频在线播放观看不卡| 大话2 男鬼变身卡| 美女cb高潮喷水在线观看| 成人免费观看视频高清| 久久99热这里只频精品6学生| 免费看不卡的av| 亚洲av福利一区| 97在线人人人人妻| 亚洲综合色惰| 蜜桃久久精品国产亚洲av| 国产午夜精品久久久久久一区二区三区| 国产老妇伦熟女老妇高清| 黑丝袜美女国产一区| 亚洲欧美日韩卡通动漫| 欧美日本中文国产一区发布| 国产高清国产精品国产三级| 欧美高清成人免费视频www| 国产国拍精品亚洲av在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产极品粉嫩免费观看在线 | 伊人久久国产一区二区| 九九久久精品国产亚洲av麻豆| 美女视频免费永久观看网站| 国产黄频视频在线观看| 久久久久国产精品人妻一区二区| 久久免费观看电影| 国产成人aa在线观看| 在线看a的网站| 亚洲成人手机| 黑丝袜美女国产一区| 高清不卡的av网站| 新久久久久国产一级毛片| 成人美女网站在线观看视频| 新久久久久国产一级毛片| 日日啪夜夜撸| 最近手机中文字幕大全| 亚洲人成网站在线观看播放| 日本午夜av视频| 亚洲精品日本国产第一区| 国产成人一区二区在线| freevideosex欧美| 91久久精品国产一区二区三区| 青春草国产在线视频| 十八禁网站网址无遮挡 | 日本欧美国产在线视频| 国产伦理片在线播放av一区| 精品一区二区免费观看| 中文字幕人妻熟人妻熟丝袜美| 欧美老熟妇乱子伦牲交| 菩萨蛮人人尽说江南好唐韦庄| 国产日韩欧美视频二区| 美女脱内裤让男人舔精品视频| 在线观看国产h片| 国产日韩欧美亚洲二区| 日韩人妻高清精品专区| 男男h啪啪无遮挡| 亚洲第一av免费看| 午夜日本视频在线| 大香蕉97超碰在线| 免费大片黄手机在线观看| 国产精品国产三级国产av玫瑰| 美女视频免费永久观看网站| 亚洲av欧美aⅴ国产| 啦啦啦中文免费视频观看日本| 大又大粗又爽又黄少妇毛片口| 日韩中字成人| 亚洲精品第二区| 国产精品一区二区性色av| 久久韩国三级中文字幕| 少妇被粗大的猛进出69影院 | 一级毛片我不卡| 久热这里只有精品99| 少妇人妻一区二区三区视频| 午夜91福利影院| 欧美一级a爱片免费观看看| 国产一区二区在线观看日韩| 我要看黄色一级片免费的| 亚洲av成人精品一区久久| 两个人的视频大全免费| 久久久精品94久久精品| 18+在线观看网站| 成人影院久久| 人人妻人人爽人人添夜夜欢视频 | 欧美一级a爱片免费观看看| 精品午夜福利在线看| 中文字幕久久专区| 激情五月婷婷亚洲| av视频免费观看在线观看| 精品一区在线观看国产| 国产伦理片在线播放av一区| 国产高清有码在线观看视频| 日韩熟女老妇一区二区性免费视频| 美女福利国产在线| 欧美区成人在线视频| 一本久久精品| 高清黄色对白视频在线免费看 | 亚洲天堂av无毛| 久久久久久久大尺度免费视频| av免费观看日本| 中文字幕免费在线视频6| 久久久久精品久久久久真实原创| 精品亚洲成a人片在线观看| 一本色道久久久久久精品综合| 观看免费一级毛片| 精品久久国产蜜桃| 国产精品人妻久久久久久| 久久ye,这里只有精品| 亚洲无线观看免费| 日日摸夜夜添夜夜爱| 国产男女超爽视频在线观看| 大香蕉97超碰在线| 男女边摸边吃奶| 99热6这里只有精品| 日本与韩国留学比较| 久热久热在线精品观看| 国产精品秋霞免费鲁丝片| 精品亚洲成a人片在线观看| 黑人高潮一二区| 在线观看一区二区三区激情| 激情五月婷婷亚洲| 精品人妻熟女毛片av久久网站| 嫩草影院新地址| 国产白丝娇喘喷水9色精品| 最后的刺客免费高清国语| 成年美女黄网站色视频大全免费 | 少妇被粗大的猛进出69影院 | 日韩中字成人| 国产日韩欧美亚洲二区| 免费久久久久久久精品成人欧美视频 | 多毛熟女@视频| 91精品一卡2卡3卡4卡| 一级二级三级毛片免费看| 91精品一卡2卡3卡4卡| 久久人人爽av亚洲精品天堂| 性高湖久久久久久久久免费观看| 在线观看免费高清a一片| 五月开心婷婷网| 亚洲国产最新在线播放| 国产成人精品久久久久久| 日韩亚洲欧美综合| 一本色道久久久久久精品综合| 亚洲精品456在线播放app| 视频区图区小说| 乱码一卡2卡4卡精品| 大片电影免费在线观看免费| a级片在线免费高清观看视频| 黄色怎么调成土黄色| h日本视频在线播放| 精品久久久精品久久久| 天天躁夜夜躁狠狠久久av| 人妻系列 视频| 日韩中字成人| 国产精品麻豆人妻色哟哟久久| 多毛熟女@视频| 在线精品无人区一区二区三| 午夜91福利影院| 国产免费一级a男人的天堂| 在线免费观看不下载黄p国产| 久久国产乱子免费精品| 又爽又黄a免费视频| 少妇的逼好多水| 亚洲国产精品999| 少妇猛男粗大的猛烈进出视频| 插逼视频在线观看| 51国产日韩欧美| 亚洲天堂av无毛| 极品少妇高潮喷水抽搐| 黑人猛操日本美女一级片| 十八禁网站网址无遮挡 | 中文天堂在线官网| 男人狂女人下面高潮的视频| 极品人妻少妇av视频| 美女中出高潮动态图| 18禁在线播放成人免费| 波野结衣二区三区在线| 亚洲国产精品成人久久小说| 日本av免费视频播放| 伊人久久国产一区二区| 欧美少妇被猛烈插入视频| 国产成人freesex在线| 一级毛片我不卡| a级毛片在线看网站| 免费看日本二区| 女人精品久久久久毛片| 久热这里只有精品99| 最近中文字幕高清免费大全6| 夫妻性生交免费视频一级片| 国产精品熟女久久久久浪| 只有这里有精品99| 97在线人人人人妻| 国产视频首页在线观看| 欧美日本中文国产一区发布| 夜夜骑夜夜射夜夜干| 午夜日本视频在线| 久久精品国产亚洲av天美| 一级片'在线观看视频| 亚洲,欧美,日韩| 国产免费一区二区三区四区乱码| 久久久亚洲精品成人影院| 美女视频免费永久观看网站| 日韩三级伦理在线观看| 国产毛片在线视频| 国产在视频线精品| 亚洲精品,欧美精品| 日日啪夜夜爽| 欧美国产精品一级二级三级 | 五月天丁香电影| 国产午夜精品久久久久久一区二区三区| 亚洲av在线观看美女高潮| 国产亚洲午夜精品一区二区久久| 精品酒店卫生间| 国产老妇伦熟女老妇高清| 午夜av观看不卡| 极品少妇高潮喷水抽搐| 久久久久人妻精品一区果冻| av福利片在线观看| 亚洲欧美日韩另类电影网站| 视频区图区小说| 久久婷婷青草| 91精品国产国语对白视频| 边亲边吃奶的免费视频| 女的被弄到高潮叫床怎么办| 高清av免费在线| 看免费成人av毛片| 亚洲欧美中文字幕日韩二区| 99热国产这里只有精品6| 日韩电影二区| 国产高清有码在线观看视频| 九草在线视频观看| 久久久久精品性色| 免费播放大片免费观看视频在线观看| av在线老鸭窝| 亚洲av电影在线观看一区二区三区| 国产精品一区二区在线观看99| 亚洲人与动物交配视频| 国产伦精品一区二区三区视频9| 久热这里只有精品99| 在线观看免费高清a一片| 22中文网久久字幕| 蜜桃久久精品国产亚洲av| 国国产精品蜜臀av免费| 欧美最新免费一区二区三区| 久久国内精品自在自线图片| 亚洲熟女精品中文字幕| tube8黄色片| 我要看日韩黄色一级片| .国产精品久久| 熟女人妻精品中文字幕| 精品久久久久久电影网| 97超碰精品成人国产| 性色avwww在线观看| 国产黄频视频在线观看| 免费少妇av软件| av在线播放精品| av网站免费在线观看视频| 精品99又大又爽又粗少妇毛片| 精品久久久久久久久av| 亚洲成色77777| 久久久久国产精品人妻一区二区| 亚洲美女视频黄频| 国产精品不卡视频一区二区| 我要看黄色一级片免费的| 如日韩欧美国产精品一区二区三区 | 亚洲av综合色区一区| 18禁裸乳无遮挡动漫免费视频| 欧美+日韩+精品| 国产精品熟女久久久久浪| 精品久久久久久久久亚洲| 秋霞伦理黄片| 99九九在线精品视频 | 乱码一卡2卡4卡精品| 一级毛片久久久久久久久女| 在线看a的网站| 国产91av在线免费观看| 亚洲精品国产色婷婷电影| 熟女电影av网| 色视频www国产| 国产欧美另类精品又又久久亚洲欧美| 在线精品无人区一区二区三| 亚洲国产精品成人久久小说| 精品国产乱码久久久久久小说| 精品一品国产午夜福利视频| 一区二区三区免费毛片| 亚洲精品久久久久久婷婷小说| 涩涩av久久男人的天堂| 校园人妻丝袜中文字幕| 美女大奶头黄色视频| 久久久久久久久久久丰满| 一级毛片久久久久久久久女| 永久免费av网站大全| av线在线观看网站| 在线观看国产h片| 国产成人精品久久久久久| 国产91av在线免费观看| av天堂久久9| 欧美xxⅹ黑人| 内地一区二区视频在线| 狂野欧美激情性xxxx在线观看| 全区人妻精品视频| 麻豆精品久久久久久蜜桃| 国产男人的电影天堂91| 国产精品一区二区在线观看99| 久久精品久久久久久久性| 又粗又硬又长又爽又黄的视频| 三级经典国产精品| 国产日韩欧美视频二区| 欧美精品国产亚洲| av福利片在线观看| 亚洲精品456在线播放app| 欧美国产精品一级二级三级 | 99热全是精品| 亚洲av福利一区| .国产精品久久| 男人狂女人下面高潮的视频| 男女边吃奶边做爰视频| 一级毛片黄色毛片免费观看视频| 晚上一个人看的免费电影| 噜噜噜噜噜久久久久久91| 在线看a的网站| 国产乱人偷精品视频| 久久青草综合色| 最新的欧美精品一区二区| 亚洲三级黄色毛片|