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

    玻璃微珠液滴碰撞分離過(guò)程數(shù)值研究

    2024-01-05 00:15:20李慧玲胡曉磊余子寒謝能剛
    關(guān)鍵詞:液橋矢量圖微珠

    李慧玲, 胡曉磊, 余子寒, 謝能剛

    (安徽工業(yè)大學(xué) 機(jī)械工程學(xué)院, 安徽 馬鞍山 243032)

    0 引 言

    玻璃微珠是近年來(lái)發(fā)展起來(lái)的一種用途廣泛、性能特殊的新型無(wú)機(jī)非金屬材料.玻璃微珠具有強(qiáng)度高、流動(dòng)性好、光學(xué)性能優(yōu)異等顯著特點(diǎn),因此被廣泛應(yīng)用于軍事、航空、建材及其他高科技領(lǐng)域作光反射材料、凈化拋光材料、填充材料、載體材料等.隨著科學(xué)技術(shù)的不斷發(fā)展和用戶對(duì)產(chǎn)品性能要求的不斷提高,玻璃微珠的直徑越來(lái)越小,范圍可達(dá)到數(shù)微米到數(shù)百微米[1].玻璃微珠的生產(chǎn)方法分為物理法和化學(xué)法,采用物理法制備玻璃微珠的技術(shù)發(fā)展較為成熟.生產(chǎn)中常用的固相玻璃粉末法就是一種物理法,該方法常將玻璃粉料和燃料噴入球化爐內(nèi),使玻璃粉末顆粒充分受熱熔化,獲得熔融的玻璃微珠液滴,在此過(guò)程中玻璃微珠液滴之間大多發(fā)生碰撞分離現(xiàn)象,但由于爐內(nèi)溫度過(guò)高,該微觀運(yùn)動(dòng)過(guò)程難以捕捉[2].

    為準(zhǔn)確預(yù)測(cè)玻璃微珠液滴碰撞過(guò)程,需要對(duì)玻璃微珠液滴碰撞過(guò)程進(jìn)行數(shù)值模擬.液滴碰撞過(guò)程是一個(gè)涉及多尺度、非線性問(wèn)題的復(fù)雜流動(dòng)過(guò)程,對(duì)于一些時(shí)間和空間的微觀信息,實(shí)驗(yàn)難以捕捉.與實(shí)驗(yàn)研究相比,數(shù)值模擬在流場(chǎng)及參數(shù)設(shè)置方面具有更大柔性.在此方面,夏盛勇等[3]采用流體體積(VOF)方法對(duì)直徑為10~200 μm的三氧化二鋁液滴碰撞過(guò)程進(jìn)行了數(shù)值模擬,獲得了反彈、大變形后聚合和自反分離等結(jié)果類型,研究了不同Weber數(shù)下的結(jié)果,并得到了它們之間的臨界Weber數(shù).Fostiropoulos等[4]研究了水乳化重油燃料爆炸液滴的碰撞行為,用3個(gè)VOF傳輸方程分離共存水液體和水蒸氣流體界面,分析模擬結(jié)果,得到液滴破碎只發(fā)生在燃料-空氣界面的一部分,類似于噴氣的特征.尹強(qiáng)等[5]將VOF方法用于研究循環(huán)冷卻系統(tǒng)中海水液滴碰撞過(guò)程,得到了聚合和自反分離兩種結(jié)果類型及二者的臨界Weber數(shù),分析了流場(chǎng)結(jié)構(gòu)、液滴直徑和海水濃度對(duì)碰撞結(jié)果的影響.

    Amani等[6]使用守恒的水平集(level set)方法模擬了所有體系中正面和偏離中心的二元液滴碰撞,介紹了一種新的薄層穩(wěn)定方法來(lái)數(shù)值解析薄膜,并對(duì)每個(gè)案例進(jìn)行了深入的能量分析,涵蓋了廣泛的碰撞范圍,為碰撞過(guò)程提供了更深入的見(jiàn)解.Ashna等[7]采用LBM方法研究了液滴的蒸發(fā)過(guò)程,引入了有效表面張力測(cè)量蒸發(fā)速率對(duì)碰撞結(jié)果的影響,得到了其隨Stefan數(shù)或Weber數(shù)增加而減少,隨時(shí)間的增加而增加的規(guī)律.Liu等[8]通過(guò)將實(shí)驗(yàn)可視化與數(shù)值模擬相結(jié)合,研究了雙乳化液滴的水動(dòng)力剪驅(qū)動(dòng)二元碰撞,發(fā)現(xiàn)了過(guò)渡和反轉(zhuǎn)這兩種典型碰撞運(yùn)動(dòng),并得到了它們之間的決定影響因素.Zhang等[9]用分子動(dòng)力學(xué)的方法研究了尺寸對(duì)二元碰撞納米液滴動(dòng)力學(xué)的影響,對(duì)不同直徑液滴的拉伸因子、能量耗散和碰撞結(jié)果進(jìn)行分析,得到了納米尺度上原子相互作用引起的動(dòng)能耗散對(duì)液滴分離的發(fā)生或其他方面有顯著影響.

    Li等[10]采用CLSVOF方法的數(shù)值模型,研究了空心液滴對(duì)干燥平面沖擊的動(dòng)力學(xué)和傳熱問(wèn)題,比較了高空化氣泡壓力的空心液滴與連續(xù)致密液滴的動(dòng)力學(xué)和傳熱特性,得到了液滴在這兩種條件下的碰撞規(guī)律.Qian等[11]利用CLSVOF方法和自適應(yīng)網(wǎng)格細(xì)化,對(duì)Newton液體和聚合物液體的碰撞行為進(jìn)行了數(shù)值研究.對(duì)液滴的變形過(guò)程、內(nèi)部流場(chǎng)和能量演化進(jìn)行了詳細(xì)分析,得到了二者在液滴碰撞過(guò)程中的不同之處,主要分析了聚合物液滴碰撞的一些變化特征.由上述文獻(xiàn)可知,國(guó)內(nèi)外對(duì)液滴碰撞機(jī)理的研究已然十分成熟.然而,截至目前,國(guó)內(nèi)外尚未見(jiàn)公開(kāi)發(fā)表針對(duì)玻璃微珠液滴碰撞物理規(guī)律及模型的研究.

    鑒于對(duì)熔融玻璃微珠液滴碰撞過(guò)程實(shí)驗(yàn)觀測(cè)的困難性,本文采用CLSVOF方法[12]開(kāi)展兩個(gè)相同尺寸玻璃微珠液滴對(duì)心碰撞的數(shù)值模擬.首先闡述了液滴碰撞的基本模型以及相關(guān)參數(shù)計(jì)算公式,隨后建立了玻璃微珠碰撞過(guò)程數(shù)值模型和計(jì)算方法,采用建立的數(shù)值方法對(duì)Qian等[13]的實(shí)驗(yàn)開(kāi)展數(shù)值驗(yàn)證,最后研究了玻璃微珠液滴對(duì)心碰撞分離的物理規(guī)律.

    1 液滴碰撞基本模型

    圖1給出了兩個(gè)相同尺寸液滴對(duì)心碰撞模型示意圖,其中D0是液滴直徑,v是液滴速度,L為兩液滴的初始間距.計(jì)算區(qū)域?yàn)?D0×2.5D0.由于模型具有軸對(duì)稱特性,本文采用二維軸對(duì)稱計(jì)算條件,計(jì)算域底部為軸對(duì)稱邊界條件,其他三邊均為壓力出口,計(jì)算區(qū)域如圖1所示.

    圖1 物理模型Fig. 1 The physical model

    2 數(shù)值計(jì)算方法

    液滴碰撞是兩相流問(wèn)題,目前,兩相界面的捕獲方法主要有:VOF法[14]、水平集[15]以及 CLSVOF方法.VOF法難于計(jì)算相界面曲率,會(huì)導(dǎo)致解的振蕩以及界面的陡峭變化被抹平.水平集法在計(jì)算時(shí)不可避免地引入數(shù)值耗散和舍入誤差,造成封閉界面內(nèi)流體質(zhì)量損失,從而導(dǎo)致質(zhì)量不守恒.VOF法和水平集法具有各自的優(yōu)勢(shì)和劣勢(shì),但是這兩種方法的優(yōu)、劣勢(shì)正好可以互補(bǔ).因此,Sussman等[16]提出使用VOF法和水平集法相耦合的方法,即CLSVOF方法進(jìn)行相界面追蹤,該方法能準(zhǔn)確捕獲相界面的同時(shí),保證質(zhì)量守恒.

    對(duì)于流體力學(xué)問(wèn)題,涉及到的基本控制方程如下所示:

    ?·v=0,

    (1)

    (2)

    式中v為速度矢量,m/s2;p為流體壓力,Pa;F為表面張力源項(xiàng);g為重力加速度,m/s2;μ為流體動(dòng)力黏度.CLSVOF方法在捕獲兩相界面時(shí)涉及到的水平集函數(shù)為

    (3)

    式中Ω1,Ω2分別是液相和氣相,Γ為兩相界面,L為區(qū)域內(nèi)點(diǎn)到界面的距離.

    通過(guò)引入Heaviside函數(shù)來(lái)使界面的密度和黏度光滑過(guò)渡:

    (4)

    式中a為1.5倍的最小網(wǎng)格尺寸.

    式(2)的表面張力源項(xiàng)F為

    F=σκ(φ)?H(φ)n,

    (5)

    式中,σ是液相的表面張力系數(shù),n是兩相界面的法線,κ(φ)是兩相界面的曲率.法線n和曲率κ(φ)的表達(dá)式如下:

    (6)

    (7)

    光滑后的密度和黏度為

    ρ(φ)=ρg+(ρl-ρg)H(φ),

    (8)

    μ(φ)=μg+(μl-μg)H(φ),

    (9)

    式中ρ(φ)為各相的平均密度,kg/m3;ρg為氣相的密度,kg/m3;ρl為液相的密度,kg/m3;μ(φ)為各相的平均動(dòng)力黏度,N·s/m2;μg為氣相的動(dòng)力黏度,N·s/m2;μl為液相的動(dòng)力黏度,N·s/m2.

    對(duì)于碰撞結(jié)果描述涉及到的Weber數(shù)和Reynolds數(shù)公式定義如下:

    (10)

    (11)

    式中,ρ為液滴密度,kg/m3;d為液滴直徑,m;u為液滴速度,m/s;σ為表面張力系數(shù),N/m;μ為液滴動(dòng)力黏度,N·s/m2.

    采用PISO方法耦合求解速度場(chǎng)和壓力場(chǎng),時(shí)間項(xiàng)采用Euler方法進(jìn)行離散,壓力項(xiàng)采用PRESTO方法進(jìn)行求解,動(dòng)量方程采用二階迎風(fēng)格式進(jìn)行求解,水平集函數(shù)采用一階迎風(fēng)格式離散求解.

    3 數(shù)值算法的驗(yàn)證

    為驗(yàn)證數(shù)值方法是否適合計(jì)算液滴碰撞及其可靠性,本文針對(duì)文獻(xiàn)[13]正十四烷液滴在1 atm(1 atm=101 325 Pa)的常溫氮?dú)猸h(huán)境下的實(shí)驗(yàn)進(jìn)行數(shù)值模擬驗(yàn)證計(jì)算,計(jì)算流體物性參數(shù)見(jiàn)表1.

    針對(duì)文獻(xiàn)[13]中圖4(h)近對(duì)心碰撞工況開(kāi)展了數(shù)值計(jì)算,計(jì)算結(jié)果如圖2所示,其中實(shí)驗(yàn)結(jié)果為三維拍攝圖.對(duì)比計(jì)算及實(shí)驗(yàn)結(jié)果可以看出,計(jì)算最終結(jié)果是自反分離,取得了與實(shí)驗(yàn)基本一致的液滴變形過(guò)程,只是在時(shí)刻上略微有所不同.

    表1 正十四烷液滴的物理參數(shù)

    (a) 實(shí)驗(yàn)結(jié)果[13] (b) 計(jì)算結(jié)果 (a) Experimental results[13] (b) Calculation results

    4 計(jì)算結(jié)果及分析

    4.1 玻璃微珠液滴及氣相物性參數(shù)

    玻璃微珠液滴的尺寸一般介于幾微米到幾百微米之間,本文選擇兩個(gè)尺寸相同的玻璃微珠液滴計(jì)算,直徑均為20 μm.模擬中所采用的物性參數(shù)是參考彭壽等[2]研究玻璃微珠球化時(shí)所用的參數(shù),如表2所示.

    表2 玻璃微珠液滴的物理參數(shù)

    4.2 玻璃微珠液滴對(duì)心碰撞數(shù)值研究

    本文僅研究?jī)蓚€(gè)相同尺寸玻璃微珠液滴對(duì)心碰撞分離的規(guī)律特性.在Weber數(shù)為149的工況下,液滴碰撞分離形成了多個(gè)衛(wèi)星液滴,但最終只有一個(gè)較大的衛(wèi)星液滴保持穩(wěn)定形態(tài),因此選取它來(lái)進(jìn)行玻璃微珠液滴碰撞分離的規(guī)律特性研究.通過(guò)對(duì)兩玻璃微珠液滴碰撞過(guò)程中的形態(tài)變化、 能量變化、 中心處壓力變化和關(guān)鍵時(shí)刻的速度矢量圖以及壓力云圖的分析, 研究玻璃微珠液滴碰撞后發(fā)生分離的原因.模擬結(jié)果如圖3所示.

    圖3為兩相同尺寸玻璃微珠液滴碰撞后自反分離的計(jì)算結(jié)果.兩液滴突破氣膜后發(fā)生聚合,0.5 μs聚合后在表面張力的作用下發(fā)生徑向伸展變形,形成片狀射流;在2.0 μs時(shí)液滴徑向拉伸到極限,之后液滴徑向收縮,并同時(shí)向軸向伸展;在6.0 μs時(shí),液滴軸向拉伸到極限,不再延展,但徑向仍在不斷收縮,形成液橋;在11.5 μs時(shí)液滴發(fā)生分離,液橋分裂形成3個(gè)衛(wèi)星液滴,在表面張力的作用下,大衛(wèi)星液滴兩側(cè)的小衛(wèi)星液滴迅速與左右兩個(gè)大液滴融合,最后形成兩個(gè)完整的液滴.

    圖3 計(jì)算結(jié)果(We=149, Re=135.52, Ur=44 m/s)Fig. 3 Calculation results (We=149, Re=135.52, Ur=44 m/s)

    4.2.1 玻璃微珠液滴對(duì)心碰撞過(guò)程中的能量變化

    許多研究表明在液滴碰撞過(guò)程中,存在3種能量即液滴動(dòng)能Ek、表面能Es和黏性耗散能Ev在發(fā)生變化[17-18].為研究玻璃微珠液滴對(duì)心碰撞過(guò)程中的能量變化情況,本文對(duì)上述3種能量進(jìn)行計(jì)算分析.液滴碰撞過(guò)程能量守恒可表示為

    Ek0+Es0=Ekn+Esn+Ev,

    (12)

    式中Ek0,Es0分別是液滴初始動(dòng)能和初始表面能;Ekn,Esn和Ev分別是液滴碰撞過(guò)程中的動(dòng)能、表面能和黏性耗散能.初始狀態(tài)下的動(dòng)能和表面能計(jì)算公式如下:

    (13)

    Es0=σS=4πσR2,

    (14)

    式中ρ為液滴密度;R為液滴半徑;v0=1/(2Ur),Ur為液滴碰撞的相對(duì)速度;σ為表面張力系數(shù).

    碰撞過(guò)程中的動(dòng)能和表面能需要計(jì)算液滴的體積和表面積,由于模擬環(huán)境為二維平面,無(wú)法直接得到液滴的體積和表面積.為此本文使用當(dāng)量半徑Re,它是由現(xiàn)有的二維平面的液滴面積演化而來(lái)的.我們假設(shè)二維平面的液滴面積S為三維液滴的中心截面面積,并且把它等量成球的中心截面Se,則該球的半徑即為當(dāng)量半徑Re,計(jì)算公式為

    (15)

    液滴碰撞過(guò)程中的動(dòng)能計(jì)算公式為

    (16)

    式中v為液滴速度.

    液滴碰撞過(guò)程中的表面能計(jì)算公式為

    (17)

    液滴碰撞過(guò)程中的黏性耗散能計(jì)算公式為

    Ev=Ek0+Es0-Ekn-Esn.

    (18)

    由上述公式計(jì)算出各個(gè)時(shí)刻相應(yīng)的能量,并將它們與初始總能量即Ek0+Es0相除進(jìn)行無(wú)量綱化處理,得到We=149時(shí)玻璃微珠液滴碰撞分離過(guò)程的能量變化曲線,如圖4所示.圖中Es,Ek和Ev曲線分別代表表面能、動(dòng)能和黏性耗散能的變化情況.從圖中可以看出,在玻璃微珠液滴相互靠近,突破氣膜阻力,碰撞后軸向收縮,徑向拉伸的過(guò)程中,動(dòng)能先迅速下降,在液滴即將拉伸到極限時(shí),動(dòng)能緩慢減小,在2.0 μs動(dòng)能減小至零;由于液滴融合,表面積減小,表面能緩慢下降,在2.0 μs降到最小;由于動(dòng)能和表面能減小,黏性耗散能持續(xù)增大,在2.0 μs增至最大.2.0~2.5 μs動(dòng)能和表面能基本保持不變,這是因?yàn)橐旱屋S向繼續(xù)收縮,徑向變化不大,在2.5 μs軸向收縮到極限.2.5 μs后,液滴徑向收縮,軸向拉伸,在4.0 μs徑向收縮和軸向拉伸達(dá)到平衡,11.5 μs液滴發(fā)生分離,產(chǎn)生3個(gè)衛(wèi)星液滴.在這個(gè)過(guò)程中,動(dòng)能先增大,當(dāng)徑向收縮和軸向拉伸達(dá)到平衡時(shí),動(dòng)能增至最大,隨后在表面張力的作用下,液滴徑向收縮,軸向拉伸形成液橋,在二者的持續(xù)作用下,液橋逐漸變細(xì),在12.0 μs被夾斷,在此期間,動(dòng)能緩慢減少至零.液滴在2.5 μs后,由于徑向收縮、軸向拉伸,表面積增大,表面能逐漸增大,在6.0 μs液滴軸向拉伸到極限時(shí)增加到最大值.隨后由于液滴徑向不斷收縮,液橋變細(xì),導(dǎo)致表面能小幅減?。谝簶虮粖A斷,形成3個(gè)大小不一的衛(wèi)星液滴后,由于小衛(wèi)星液滴和大液滴的融合,表面能發(fā)生小幅振蕩,最后和初始表面能相當(dāng),符合VOF方法的質(zhì)量守恒特征.黏性耗散能由于動(dòng)能和表面能的增加,急劇下降,在液滴收縮拉伸達(dá)到平衡后緩慢增長(zhǎng).

    通過(guò)上述分析可以得出,液滴動(dòng)能是前期液滴發(fā)生碰撞融合的主要?jiǎng)恿?在整個(gè)碰撞過(guò)程中對(duì)于液滴的收縮拉伸變化有著重要作用,在收縮拉伸過(guò)程中,動(dòng)能會(huì)出現(xiàn)峰值.液滴在碰撞后徑向拉伸過(guò)程表面能減小;在軸向拉伸至分離過(guò)程中,表面能先增大后減小,恢復(fù)至初始表面能水平.液滴碰撞過(guò)程中黏性耗散能變化受動(dòng)能和表面能影響,但從最終結(jié)果來(lái)看,液滴初始動(dòng)能絕大部分轉(zhuǎn)化為黏性耗散能.由此可以看出,液滴碰撞過(guò)程中的分離結(jié)果受動(dòng)能的起伏變化和表面能的變化的影響較大,在軸向拉伸前期,動(dòng)能增大,表面能增大,為后面液滴拉伸分離提供了能量支持.

    4.2.2 玻璃微珠液滴總壓變化

    玻璃微珠液滴在發(fā)生碰撞的過(guò)程中,其靜壓和動(dòng)壓都在發(fā)生變化,只分析其中一個(gè)都比較片面.而總壓為兩者之和,可以較為全面地反映液滴的壓力變化規(guī)律.圖5為兩液滴碰撞過(guò)程液滴平均總壓變化折線圖.從圖中可以看出, 兩液滴從零時(shí)刻開(kāi)始相互靠近,壓縮氣膜,液滴總壓升高.0.5~2.0 μs為兩液滴碰撞融合,并徑向拉伸到極限的過(guò)程,液滴總壓持續(xù)下降,在2.0 μs降至最低.2.0~4.0 μs為液滴徑向收縮,軸向拉伸達(dá)到平衡的過(guò)程,液滴總壓持續(xù)增大,在4.0 μs達(dá)到平衡時(shí)增至最大.4.0~6.0 μs為液滴徑向收縮,軸向拉伸到極限的過(guò)程,液滴總壓減小.6.0 μs之后,液滴軸向不再拉伸,徑向仍在不斷收縮,中間液橋部分形態(tài)變化均勻,因此液滴總壓保持不變.在11.5 μs時(shí),可以明顯看出液橋形態(tài)發(fā)生變化,兩端變細(xì),液滴總壓也忽然增大,這與末端夾斷機(jī)制一致[19-20].液橋斷裂之后,小衛(wèi)星液滴和大液滴發(fā)生融合,形成穩(wěn)定形態(tài)的過(guò)程中,液滴總壓降至穩(wěn)定狀態(tài).

    4.2.3 玻璃微珠壓力云圖和速度矢量圖分析

    從上述玻璃微珠液滴能量變化和總壓變化分析可知,液滴碰撞分離過(guò)程有4個(gè)重要狀態(tài),即:2.0 μs徑向拉伸到極限,4.0 μs徑向收縮和軸向拉伸達(dá)到平衡狀態(tài),6.0 μs軸向拉伸到極限和11.5 μs液滴液橋夾斷分離.下面對(duì)這4種狀態(tài)的壓力云圖和速度矢量圖進(jìn)行分析.圖6給出的是兩相同尺寸玻璃微珠液滴發(fā)生碰撞2.0 μs時(shí),即液滴徑向拉伸到極限的速度矢量圖(圖6(a))(紅色箭頭代表液相速度矢量,藍(lán)色箭頭代表氣相速度矢量,后同)和液滴內(nèi)部壓力云圖(圖6(b)).從液滴的速度矢量圖可以看出液滴的速度分布上下對(duì)稱,上方的液滴拉伸部分速度向上,上端部分速度向下擴(kuò)散;外部的速度在液滴周圍呈大小不一的漩渦擴(kuò)散,總體趨勢(shì)向上,但上端部分漩渦有向下壓縮液滴的趨勢(shì).下方液滴速度分布和上方相同,但方向相反,具有很好的對(duì)稱性.液滴拉伸部分向上下兩側(cè)擴(kuò)散的速度矢量箭頭密集,但兩端的速度矢量箭頭稀疏,這是由于啞鈴形玻璃微珠液滴中間部分較細(xì),其界面曲率κ較大,在此處的表面張力較大,即式(2)中的動(dòng)力源項(xiàng)F較大,使此處產(chǎn)生指向液滴內(nèi)部的較大速度.這也表示液滴徑向拉伸已達(dá)到極限,沒(méi)有持續(xù)拉伸的趨勢(shì)了,與此時(shí)刻之后液滴的形態(tài)變化相符.此時(shí)的液滴形態(tài)呈啞鈴形,液滴壓力云圖與之相符,且壓力是從兩端向中心逐漸減小,這和圖5中此時(shí)刻的平均總壓變化相符.

    圖4 玻璃微珠液滴碰撞分離過(guò)程能量變化曲線

    (a) 速度矢量圖 (b) 壓力云圖 (a) Velocity vectors (b) Pressure contours

    圖7為玻璃微珠液滴徑向收縮和軸向拉伸達(dá)到平衡時(shí)刻4.0 μs的速度矢量圖和內(nèi)部壓力云圖.從圖7(a)可以看出,此時(shí)液滴的速度大多集中在左右兩端,并分別向兩端擴(kuò)散.液滴兩端的速度矢量箭頭明顯比較密集,這是由于液滴形態(tài)呈鈴鐺形狀,中間左右兩端部分突出,其界面曲率κ就較大,則在此處的表面張力,即式(2)中的動(dòng)力源項(xiàng)F也就越大,使此處液體的速度越大.圖7(b)為玻璃微珠液滴壓力云圖,圖中顯示液滴壓力從上下兩端向中心逐漸增大,從中心向左右逐漸增大.液滴壓力和速度的分布的共同作用促成了此刻的液滴形態(tài).

    圖8為玻璃微珠液滴軸向拉伸到極限時(shí)刻6.0 μs的速度矢量圖和內(nèi)部壓力云圖.從圖8(a)可以看出,此時(shí)液滴的速度分布均勻,并且從中間分別向左右兩端擴(kuò)散,具有良好的對(duì)稱性.此時(shí)的液滴形態(tài)呈柱狀,形態(tài)勻稱,界面曲率κ沒(méi)有明顯變化,所以液滴速度分布均勻.在端部分布有速度漩渦,與液滴內(nèi)部速度方向相反,阻礙液滴的軸向拉伸.圖8(b)為玻璃微珠液滴壓力云圖,兩端的壓力比中部壓力大,但兩者相差不大.速度矢量圖和壓力云圖都與液滴的形態(tài)變化特性相符.

    (a) 速度矢量圖 (b) 壓力云圖 (a) Velocity vectors (b) Pressure contours

    (a) 速度矢量圖 (b) 液滴內(nèi)部壓力云圖 (a) Velocity vectors (b) Pressure contours

    (a) 速度矢量圖(a) Velocity vectors

    (b) 液橋局部放大速度矢量圖 (c) 壓力云圖 (b) Partially amplified velocity vectors of the liquid bridge (c) Pressure contours

    圖9為兩相同尺寸玻璃微珠液滴碰撞分離時(shí)刻11.5 μs的速度矢量圖和壓力云圖.從圖9(a)可以看出,此時(shí)液滴的速度大多集中在液橋部分,液滴形態(tài)呈啞鈴狀.圖9(b)為圖9(a)液橋部分速度矢量圖的放大圖,從中可以清晰地看出液橋部分的速度從兩端向中部擴(kuò)散,液橋外部末端靠近大液滴處有速度漩渦形成,速度從中心向兩端擴(kuò)散.圖9(c)為玻璃微珠液滴壓力云圖,圖中顯示液滴壓力都集中在液橋的部位,并從液橋中心向左右兩邊增大,導(dǎo)致末端壓力最大.綜合速度和壓力分析可知兩者的分布與液滴的末端夾斷機(jī)制一致,這也是導(dǎo)致液滴碰撞分離的重要原因.

    5 結(jié) 論

    本文對(duì)玻璃微珠液滴碰撞分離過(guò)程進(jìn)行了數(shù)值模擬研究,闡述了計(jì)算方法,并驗(yàn)證了計(jì)算方法可靠性,主要結(jié)論如下:

    1) 在碰撞過(guò)程中,液滴的初始表面能和最終表面能相當(dāng),動(dòng)能則絕大部分轉(zhuǎn)化為黏性耗散能.

    2) 在液滴軸向拉伸分離前期,動(dòng)能增大、表面能增大為后期液滴分離供能,對(duì)玻璃微珠液滴碰撞后分離有著重要作用.

    3) 通過(guò)對(duì)玻璃微珠液滴能量和平均總壓變化的研究,得到了液滴碰撞分離過(guò)程有4種狀態(tài),即徑向拉伸到極限、徑向收縮和軸向拉伸達(dá)到平衡、軸向拉伸到極限和液滴液橋夾斷分離.

    4) 液滴的速度分布與液滴的形態(tài)密切相關(guān),不規(guī)則的液滴形態(tài)導(dǎo)致界面曲率κ較大,以此影響界面的表面張力,使得不規(guī)則的液體部分產(chǎn)生較大的不同指向的速度分布.

    5) 末端夾斷機(jī)制是玻璃微珠液滴碰撞分離的主要原因.

    猜你喜歡
    液橋矢量圖微珠
    搭起一座液橋
    空心微珠負(fù)載鈰或氮摻雜氧化亞銅光催化劑的制備方法
    基于圖像處理的平面-球面液橋力檢測(cè)
    硅酸鋁微珠在人造花崗石中的應(yīng)用
    石材(2022年1期)2022-05-23 12:48:34
    Analysis of the line current differential protection considering inverter-interfaced generation station and countermeasures
    兩狹長(zhǎng)平行板間液橋兩端體積的計(jì)算
    不等徑顆粒間液橋力學(xué)參數(shù)及形態(tài)的試驗(yàn)研究1)
    利用矢量圖對(duì)小物體從光滑斜面下滑運(yùn)動(dòng)探討
    空心玻璃微珠對(duì)PMMA/SAN共混體系相分離的影響
    空心玻璃微珠/PNHMPA/PEG復(fù)配保溫蓄熱乳膠漆的制備與表征
    久久天堂一区二区三区四区| 亚洲第一欧美日韩一区二区三区| 日日夜夜操网爽| 日韩精品中文字幕看吧| 亚洲自拍偷在线| 国产三级黄色录像| 亚洲色图 男人天堂 中文字幕| 欧美日韩亚洲国产一区二区在线观看| 色婷婷久久久亚洲欧美| 99re在线观看精品视频| 一夜夜www| 老司机午夜十八禁免费视频| 91九色精品人成在线观看| 12—13女人毛片做爰片一| 久久亚洲精品不卡| 亚洲av成人av| 久久久久久国产a免费观看| 满18在线观看网站| 99精品久久久久人妻精品| 国产97色在线日韩免费| 国产高清videossex| 亚洲视频免费观看视频| 国产亚洲欧美在线一区二区| 男女下面插进去视频免费观看| 午夜免费激情av| 欧美成狂野欧美在线观看| x7x7x7水蜜桃| 国产成人欧美在线观看| 国产精品亚洲美女久久久| 亚洲少妇的诱惑av| 国产av精品麻豆| 一级作爱视频免费观看| 成人精品一区二区免费| 午夜视频精品福利| 亚洲av日韩精品久久久久久密| 亚洲色图 男人天堂 中文字幕| 黄片播放在线免费| 日韩免费av在线播放| 91精品三级在线观看| 午夜免费激情av| 午夜久久久久精精品| 国产精品99久久99久久久不卡| 美女免费视频网站| 母亲3免费完整高清在线观看| 91麻豆av在线| 又紧又爽又黄一区二区| 色综合亚洲欧美另类图片| av网站免费在线观看视频| 国产单亲对白刺激| 香蕉久久夜色| 波多野结衣av一区二区av| 免费高清在线观看日韩| 女人精品久久久久毛片| 欧美性长视频在线观看| 中文字幕人成人乱码亚洲影| 老鸭窝网址在线观看| 人妻久久中文字幕网| 搞女人的毛片| avwww免费| 99国产精品一区二区蜜桃av| 黄色丝袜av网址大全| 777久久人妻少妇嫩草av网站| 国产真人三级小视频在线观看| 日本a在线网址| 狂野欧美激情性xxxx| 国产伦一二天堂av在线观看| 亚洲少妇的诱惑av| av网站免费在线观看视频| 国产精品久久电影中文字幕| 日本黄色视频三级网站网址| 亚洲熟妇中文字幕五十中出| 非洲黑人性xxxx精品又粗又长| 一区二区三区国产精品乱码| 久久国产乱子伦精品免费另类| 国内毛片毛片毛片毛片毛片| 亚洲精品国产一区二区精华液| 女性被躁到高潮视频| 亚洲国产看品久久| 高清毛片免费观看视频网站| 窝窝影院91人妻| 亚洲专区字幕在线| 村上凉子中文字幕在线| 国产亚洲欧美在线一区二区| 亚洲五月婷婷丁香| 国产精品野战在线观看| 好男人电影高清在线观看| 精品久久久久久久毛片微露脸| 美女 人体艺术 gogo| 国产成人精品久久二区二区免费| 亚洲自拍偷在线| 日日干狠狠操夜夜爽| 99久久综合精品五月天人人| 97人妻天天添夜夜摸| 两个人免费观看高清视频| 欧美成人一区二区免费高清观看 | 亚洲伊人色综图| 久久久久久大精品| 久久草成人影院| 少妇熟女aⅴ在线视频| 美女午夜性视频免费| cao死你这个sao货| 91老司机精品| 国产精品久久久久久亚洲av鲁大| 国产成人精品久久二区二区91| 人妻久久中文字幕网| 国产精品一区二区精品视频观看| 十八禁网站免费在线| 一级作爱视频免费观看| 在线十欧美十亚洲十日本专区| 国产真人三级小视频在线观看| 一级,二级,三级黄色视频| 正在播放国产对白刺激| 亚洲全国av大片| 亚洲性夜色夜夜综合| 在线av久久热| 亚洲成人免费电影在线观看| 亚洲性夜色夜夜综合| 国产av在哪里看| 黑丝袜美女国产一区| 国产97色在线日韩免费| 亚洲三区欧美一区| 国产精品自产拍在线观看55亚洲| 亚洲一区高清亚洲精品| 精品福利观看| 男女午夜视频在线观看| 人人妻人人爽人人添夜夜欢视频| 久久久久久久精品吃奶| 精品国产一区二区久久| 两个人看的免费小视频| av免费在线观看网站| 亚洲美女黄片视频| 他把我摸到了高潮在线观看| 国产极品粉嫩免费观看在线| 欧美激情久久久久久爽电影 | 一级毛片女人18水好多| 久久人人精品亚洲av| 亚洲av熟女| 国产片内射在线| 国产一区二区三区视频了| 国产亚洲欧美在线一区二区| 欧美成人免费av一区二区三区| 麻豆成人av在线观看| 亚洲熟妇熟女久久| 国产亚洲欧美精品永久| 亚洲国产高清在线一区二区三 | 免费看美女性在线毛片视频| 久久国产精品人妻蜜桃| 国产三级在线视频| 99国产精品一区二区蜜桃av| 欧美最黄视频在线播放免费| 欧美av亚洲av综合av国产av| 99在线视频只有这里精品首页| 欧美乱码精品一区二区三区| 1024香蕉在线观看| 操出白浆在线播放| 日韩国内少妇激情av| 国产av一区二区精品久久| 久久九九热精品免费| 一级毛片高清免费大全| 日韩欧美一区视频在线观看| 三级毛片av免费| 老司机深夜福利视频在线观看| 丰满的人妻完整版| 一进一出抽搐动态| 国产视频一区二区在线看| 国产精品日韩av在线免费观看 | 久久中文看片网| 欧美日韩福利视频一区二区| 国产精品98久久久久久宅男小说| 亚洲久久久国产精品| 三级毛片av免费| 亚洲最大成人中文| 亚洲熟妇熟女久久| 国产精品永久免费网站| 此物有八面人人有两片| 最近最新免费中文字幕在线| 国产亚洲精品综合一区在线观看 | 婷婷精品国产亚洲av在线| 满18在线观看网站| 久久国产亚洲av麻豆专区| 美女大奶头视频| 亚洲五月婷婷丁香| 亚洲va日本ⅴa欧美va伊人久久| 免费av毛片视频| 国产麻豆成人av免费视频| 国产伦一二天堂av在线观看| 久久精品国产综合久久久| 18禁观看日本| 最近最新免费中文字幕在线| 久久精品亚洲精品国产色婷小说| 日本精品一区二区三区蜜桃| 亚洲第一欧美日韩一区二区三区| 18禁观看日本| 桃红色精品国产亚洲av| 免费在线观看视频国产中文字幕亚洲| 久久香蕉国产精品| 国产av一区在线观看免费| 一区二区三区高清视频在线| 一边摸一边抽搐一进一小说| 亚洲一区二区三区不卡视频| www.999成人在线观看| 国产97色在线日韩免费| 在线观看www视频免费| 亚洲人成77777在线视频| 免费久久久久久久精品成人欧美视频| 乱人伦中国视频| 啪啪无遮挡十八禁网站| 国产精品亚洲美女久久久| 黄片小视频在线播放| 极品人妻少妇av视频| 亚洲欧美日韩另类电影网站| 亚洲七黄色美女视频| 中文字幕久久专区| 两个人免费观看高清视频| 午夜老司机福利片| 成人特级黄色片久久久久久久| videosex国产| 9色porny在线观看| 美女免费视频网站| 99久久精品国产亚洲精品| 黑人巨大精品欧美一区二区mp4| or卡值多少钱| 我的亚洲天堂| 两个人免费观看高清视频| 精品国产一区二区久久| 色哟哟哟哟哟哟| 男女之事视频高清在线观看| 亚洲国产欧美网| 亚洲激情在线av| 成人国语在线视频| 日本黄色视频三级网站网址| 国产精品 国内视频| 一区二区三区国产精品乱码| 亚洲熟妇中文字幕五十中出| 天堂影院成人在线观看| 午夜免费成人在线视频| 香蕉国产在线看| 欧美 亚洲 国产 日韩一| 亚洲久久久国产精品| 成熟少妇高潮喷水视频| www.自偷自拍.com| 香蕉国产在线看| 美女国产高潮福利片在线看| 久久久精品欧美日韩精品| 亚洲国产中文字幕在线视频| 亚洲国产毛片av蜜桃av| 欧美日韩福利视频一区二区| 999久久久精品免费观看国产| 日本免费a在线| 夜夜躁狠狠躁天天躁| 男女下面进入的视频免费午夜 | 亚洲国产中文字幕在线视频| 午夜影院日韩av| 久9热在线精品视频| 色婷婷久久久亚洲欧美| 制服丝袜大香蕉在线| 欧美一级a爱片免费观看看 | 丁香欧美五月| 久久久精品国产亚洲av高清涩受| 亚洲精品中文字幕在线视频| 久久精品亚洲熟妇少妇任你| 高清黄色对白视频在线免费看| 老司机午夜十八禁免费视频| 国产成人免费无遮挡视频| 激情视频va一区二区三区| 黄色a级毛片大全视频| 变态另类丝袜制服| 国产国语露脸激情在线看| 亚洲欧美日韩高清在线视频| 国产精品亚洲一级av第二区| av有码第一页| 欧美国产日韩亚洲一区| 女人被躁到高潮嗷嗷叫费观| 亚洲国产精品成人综合色| 国产1区2区3区精品| 男女做爰动态图高潮gif福利片 | √禁漫天堂资源中文www| 一进一出抽搐gif免费好疼| 国产一区二区三区在线臀色熟女| 精品午夜福利视频在线观看一区| 国产精品1区2区在线观看.| 久久精品国产亚洲av香蕉五月| 精品国产乱码久久久久久男人| 久久影院123| 日本免费一区二区三区高清不卡 | 国产精品久久电影中文字幕| 国产精品亚洲一级av第二区| 变态另类丝袜制服| 国产精品美女特级片免费视频播放器 | 亚洲成人免费电影在线观看| 国产一区二区激情短视频| aaaaa片日本免费| 亚洲av五月六月丁香网| 变态另类丝袜制服| 亚洲成国产人片在线观看| 欧美 亚洲 国产 日韩一| 亚洲精品久久成人aⅴ小说| 久久久久久久久久久久大奶| 色av中文字幕| 国产亚洲精品一区二区www| 国产精品综合久久久久久久免费 | 免费不卡黄色视频| 在线十欧美十亚洲十日本专区| 在线观看一区二区三区| 色精品久久人妻99蜜桃| 久久人人爽av亚洲精品天堂| 一区二区三区激情视频| 深夜精品福利| 在线观看免费午夜福利视频| 国产亚洲精品av在线| www.999成人在线观看| 国产人伦9x9x在线观看| 久久精品国产清高在天天线| 欧美大码av| 美女午夜性视频免费| 亚洲一区高清亚洲精品| 日韩成人在线观看一区二区三区| 国产精品亚洲av一区麻豆| 亚洲av成人不卡在线观看播放网| 国产精品美女特级片免费视频播放器 | 黄色视频,在线免费观看| 国产精品野战在线观看| 成人免费观看视频高清| 在线国产一区二区在线| 国产成人精品在线电影| 在线国产一区二区在线| 国产精品自产拍在线观看55亚洲| 国产亚洲精品av在线| videosex国产| 精品卡一卡二卡四卡免费| 女生性感内裤真人,穿戴方法视频| 50天的宝宝边吃奶边哭怎么回事| 在线永久观看黄色视频| 亚洲精品国产一区二区精华液| 国产精品久久久av美女十八| 久久草成人影院| 91成人精品电影| 精品国产一区二区久久| 国产av在哪里看| 久久久精品国产亚洲av高清涩受| 18禁观看日本| 很黄的视频免费| 久久人人97超碰香蕉20202| 在线观看日韩欧美| 老熟妇仑乱视频hdxx| 久9热在线精品视频| 国产一区二区三区视频了| 禁无遮挡网站| 麻豆成人av在线观看| 日本 欧美在线| 精品福利观看| 最近最新免费中文字幕在线| 精品不卡国产一区二区三区| www.www免费av| 久久国产精品影院| 高清黄色对白视频在线免费看| 一区二区三区国产精品乱码| 日韩视频一区二区在线观看| 成人永久免费在线观看视频| 国产亚洲精品久久久久5区| 一级毛片精品| 成人免费观看视频高清| 91九色精品人成在线观看| 亚洲人成77777在线视频| 麻豆一二三区av精品| 老汉色av国产亚洲站长工具| 麻豆一二三区av精品| 亚洲av日韩精品久久久久久密| 久久国产亚洲av麻豆专区| 国产精品二区激情视频| 免费观看精品视频网站| 亚洲av日韩精品久久久久久密| 精品一区二区三区av网在线观看| 中文字幕人妻丝袜一区二区| 久久热在线av| 99在线人妻在线中文字幕| 99精品在免费线老司机午夜| 国产精品日韩av在线免费观看 | 91av网站免费观看| 欧美国产精品va在线观看不卡| 两个人视频免费观看高清| 国产成+人综合+亚洲专区| 久久中文字幕人妻熟女| 欧美日韩瑟瑟在线播放| 最近最新免费中文字幕在线| 精品国产亚洲在线| 久久久国产成人免费| 国产亚洲精品久久久久久毛片| 国产精品 国内视频| 变态另类成人亚洲欧美熟女 | 身体一侧抽搐| 中文字幕人妻熟女乱码| av天堂在线播放| 国产97色在线日韩免费| 免费少妇av软件| 亚洲avbb在线观看| 国产精品,欧美在线| 国内久久婷婷六月综合欲色啪| 丰满的人妻完整版| xxx96com| 久久久水蜜桃国产精品网| 伦理电影免费视频| 亚洲视频免费观看视频| 人人澡人人妻人| 欧美日本中文国产一区发布| 精品久久久久久久人妻蜜臀av | 国产日韩一区二区三区精品不卡| 变态另类丝袜制服| 精品卡一卡二卡四卡免费| 亚洲九九香蕉| 激情在线观看视频在线高清| 黄片大片在线免费观看| 午夜精品国产一区二区电影| 十分钟在线观看高清视频www| 日本欧美视频一区| 久久伊人香网站| 亚洲第一青青草原| 久久国产亚洲av麻豆专区| 天天添夜夜摸| 亚洲午夜精品一区,二区,三区| 免费搜索国产男女视频| 亚洲三区欧美一区| 国产国语露脸激情在线看| 久久久久亚洲av毛片大全| 亚洲成国产人片在线观看| 久久狼人影院| 丰满人妻熟妇乱又伦精品不卡| 亚洲av美国av| 丝袜美足系列| 国产成人精品久久二区二区91| 国产又爽黄色视频| 精品久久久久久成人av| 岛国视频午夜一区免费看| 欧美一级a爱片免费观看看 | 亚洲男人天堂网一区| 亚洲精品av麻豆狂野| 大陆偷拍与自拍| 欧美精品啪啪一区二区三区| 一本久久中文字幕| 嫩草影视91久久| 国产精品亚洲av一区麻豆| 一区福利在线观看| 精品第一国产精品| 一个人免费在线观看的高清视频| 岛国在线观看网站| 中国美女看黄片| 亚洲专区国产一区二区| 一级作爱视频免费观看| 无遮挡黄片免费观看| 欧美中文综合在线视频| 十八禁网站免费在线| 看黄色毛片网站| 亚洲七黄色美女视频| 精品不卡国产一区二区三区| 操出白浆在线播放| 丁香欧美五月| 国产日韩一区二区三区精品不卡| 90打野战视频偷拍视频| 曰老女人黄片| 国产色视频综合| 日韩精品青青久久久久久| 久久草成人影院| 久久天堂一区二区三区四区| 日日爽夜夜爽网站| 欧美日韩亚洲综合一区二区三区_| 国产精品国产高清国产av| 两个人视频免费观看高清| 午夜两性在线视频| 久久这里只有精品19| 亚洲第一av免费看| 精品第一国产精品| 成在线人永久免费视频| 黄色成人免费大全| 国产极品粉嫩免费观看在线| 国产真人三级小视频在线观看| av天堂久久9| 在线天堂中文资源库| 欧美乱色亚洲激情| 久久精品91无色码中文字幕| 悠悠久久av| 欧美丝袜亚洲另类 | 亚洲在线自拍视频| 久久中文字幕人妻熟女| 精品不卡国产一区二区三区| 人人妻人人澡人人看| 亚洲欧美日韩无卡精品| 国产aⅴ精品一区二区三区波| 国产精品日韩av在线免费观看 | 亚洲成av人片免费观看| 日韩欧美国产一区二区入口| 法律面前人人平等表现在哪些方面| 99热只有精品国产| 亚洲午夜理论影院| 母亲3免费完整高清在线观看| 黄色 视频免费看| 欧美黄色淫秽网站| 色综合亚洲欧美另类图片| 亚洲精品av麻豆狂野| 亚洲欧美一区二区三区黑人| 不卡一级毛片| 国产区一区二久久| 在线永久观看黄色视频| 国产成人精品久久二区二区免费| 国产精品 国内视频| 国产高清视频在线播放一区| 成年女人毛片免费观看观看9| 精品国产美女av久久久久小说| 日日摸夜夜添夜夜添小说| 满18在线观看网站| 每晚都被弄得嗷嗷叫到高潮| 婷婷丁香在线五月| 亚洲精品美女久久av网站| 欧美绝顶高潮抽搐喷水| 午夜精品在线福利| 村上凉子中文字幕在线| 黑人欧美特级aaaaaa片| 国产精品精品国产色婷婷| 麻豆成人av在线观看| 精品无人区乱码1区二区| 男女做爰动态图高潮gif福利片 | 国产成年人精品一区二区| 女人爽到高潮嗷嗷叫在线视频| 日本 欧美在线| 亚洲 国产 在线| 精品久久久久久成人av| 国产一区二区激情短视频| 亚洲午夜精品一区,二区,三区| 欧美人与性动交α欧美精品济南到| 丁香欧美五月| 啦啦啦观看免费观看视频高清 | 亚洲中文字幕日韩| e午夜精品久久久久久久| av视频免费观看在线观看| 久久久久久久久中文| 一级,二级,三级黄色视频| 亚洲av成人不卡在线观看播放网| 国产av一区在线观看免费| or卡值多少钱| 午夜两性在线视频| 成人国语在线视频| 91精品三级在线观看| 亚洲最大成人中文| 在线观看日韩欧美| 国产高清激情床上av| 成人欧美大片| 亚洲五月天丁香| 国产精品影院久久| 成人18禁高潮啪啪吃奶动态图| 国内毛片毛片毛片毛片毛片| 999久久久精品免费观看国产| 狂野欧美激情性xxxx| 十八禁网站免费在线| 欧美激情久久久久久爽电影 | 99re在线观看精品视频| 中文字幕av电影在线播放| 母亲3免费完整高清在线观看| 自线自在国产av| 国产高清videossex| 俄罗斯特黄特色一大片| 精品一区二区三区视频在线观看免费| 大码成人一级视频| 亚洲精品国产精品久久久不卡| 午夜激情av网站| 黄色片一级片一级黄色片| 免费在线观看黄色视频的| 国产精品久久久久久亚洲av鲁大| 午夜久久久久精精品| 91av网站免费观看| 免费在线观看影片大全网站| 亚洲欧美精品综合久久99| 久久精品影院6| 国产乱人伦免费视频| 很黄的视频免费| 日韩欧美三级三区| 老司机深夜福利视频在线观看| 69精品国产乱码久久久| 日韩欧美一区视频在线观看| 黑人操中国人逼视频| 国产成人影院久久av| 欧美av亚洲av综合av国产av| 日日夜夜操网爽| 最近最新免费中文字幕在线| 最好的美女福利视频网| 50天的宝宝边吃奶边哭怎么回事| cao死你这个sao货| 人人妻,人人澡人人爽秒播| 午夜久久久久精精品| 性色av乱码一区二区三区2| 一进一出抽搐gif免费好疼| 亚洲第一电影网av| 99国产综合亚洲精品| 身体一侧抽搐| 国产精品 欧美亚洲| 叶爱在线成人免费视频播放| 亚洲精品在线美女| 亚洲男人的天堂狠狠| or卡值多少钱| 国产一级毛片七仙女欲春2 | 国产精品1区2区在线观看.| 两人在一起打扑克的视频| 一二三四社区在线视频社区8| 日韩精品中文字幕看吧| 亚洲人成伊人成综合网2020| 国产精品免费一区二区三区在线| 在线十欧美十亚洲十日本专区| 人人妻人人澡人人看| 日韩欧美三级三区| 亚洲午夜精品一区,二区,三区| 免费在线观看亚洲国产| 国产精品一区二区三区四区久久 | 麻豆久久精品国产亚洲av| 日本在线视频免费播放| 午夜亚洲福利在线播放| 自线自在国产av| 91老司机精品|