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

    流動(dòng)聚焦中液體錐形形態(tài)和流動(dòng)結(jié)構(gòu)實(shí)驗(yàn)研究

    2022-05-21 02:19:52康鵬郭鑒鋒穆愷司廷
    實(shí)驗(yàn)流體力學(xué) 2022年2期

    康鵬,郭鑒鋒,穆愷,司廷

    中國(guó)科學(xué)技術(shù)大學(xué) 近代力學(xué)系,合肥 230027

    0 引 言

    微納米尺度的液滴、顆粒和膠囊在科學(xué)研究和工程實(shí)踐中具有重要價(jià)值,廣泛應(yīng)用于醫(yī)學(xué)、生物學(xué)、材料科學(xué)、農(nóng)業(yè)、食品、化工等諸多領(lǐng)域。在醫(yī)學(xué)領(lǐng)域,微納膠囊可以實(shí)現(xiàn)藥物的靶向運(yùn)輸、固化和儲(chǔ)存、緩釋或可控釋放等;在農(nóng)業(yè)領(lǐng)域,利用經(jīng)過(guò)霧化的單分散農(nóng)藥液滴可以實(shí)現(xiàn)藥物的均勻噴灑并節(jié)省用量;在食品工業(yè)中,微納膠囊被普遍應(yīng)用,以保護(hù)生物活性、遮掩不良?xì)馕丁⒐潭ㄊ称芳庸ぶ械募?xì)胞或酶等。此外,微納米尺度的液滴、顆粒和膠囊在增材制造、組織工程、再生醫(yī)學(xué)等重大需求領(lǐng)域及其涉及的新興交叉學(xué)科領(lǐng)域也發(fā)揮著至關(guān)重要的作用。

    隨著科學(xué)技術(shù)的進(jìn)步,制備微納米尺度的液滴、顆粒和膠囊的方法也有了顯著的發(fā)展。在界面聚合法、乳化法等傳統(tǒng)制備技術(shù)基礎(chǔ)上,近年來(lái)逐漸出現(xiàn)了以電霧化技術(shù)、流動(dòng)聚焦技術(shù)、微流控芯片技術(shù)、玻璃微毛細(xì)管技術(shù)等為代表的新型制備技術(shù)。流動(dòng)聚焦(flow focusing)是一種典型的毛細(xì)流動(dòng)現(xiàn)象:流體從毛細(xì)管流出后,被另一種高速運(yùn)動(dòng)的流體驅(qū)動(dòng),在毛細(xì)管下方正對(duì)的小孔處聚焦,形成穩(wěn)定的錐-射流結(jié)構(gòu);在多種物理過(guò)程的共同作用下,界面擾動(dòng)逐漸發(fā)展,射流最終破碎形成微納米尺度液滴。流動(dòng)聚焦技術(shù)穩(wěn)定性好、操作簡(jiǎn)單、可控性強(qiáng)、實(shí)驗(yàn)條件寬松、材料適用性強(qiáng),且液滴制備效率高,收集的液滴經(jīng)過(guò)后處理可應(yīng)用于不同領(lǐng)域。根據(jù)核心裝置結(jié)構(gòu)的不同,流動(dòng)聚焦可以分為單軸流動(dòng)聚焦和復(fù)合流動(dòng)聚焦;根據(jù)提供驅(qū)動(dòng)力的方式不同,流動(dòng)聚焦又可分為氣驅(qū)流動(dòng)聚焦和液驅(qū)流動(dòng)聚焦。此外,在流動(dòng)聚焦技術(shù)的基礎(chǔ)上,還衍生出了將流動(dòng)聚焦與電霧化相結(jié)合的電流動(dòng)聚焦技術(shù);將流動(dòng)聚焦技術(shù)應(yīng)用于微流控芯片中,則形成了微流控流動(dòng)聚焦技術(shù)。

    一般而言,流動(dòng)聚焦中被聚焦液體在演化過(guò)程中會(huì)經(jīng)歷3種形態(tài):錐形、射流和液滴。自1998年流動(dòng)聚焦被提出以來(lái),國(guó)內(nèi)外諸多學(xué)者針對(duì)其演化過(guò)程中涉及的關(guān)鍵力學(xué)問(wèn)題在過(guò)程控制、流動(dòng)模式、尺度律和不穩(wěn)定性等方面進(jìn)行了廣泛研究。例如,在氣驅(qū)單軸流動(dòng)聚焦中,針對(duì)液體流量、速度和氣體壓差,Si等得到了流動(dòng)聚焦在這3個(gè)外部控制參數(shù)平面內(nèi)的6種流動(dòng)模式及工作區(qū)域。在液驅(qū)同軸流動(dòng)聚焦中,Mu等探討了復(fù)合錐-射流的流動(dòng)模態(tài)及模態(tài)轉(zhuǎn)換,并針對(duì)射流模態(tài)下兩層界面破碎的耦合規(guī)律及物理機(jī)理開(kāi)展了研究。在單界面的單軸流動(dòng)聚焦和雙界面的同軸流動(dòng)聚焦的基礎(chǔ)上,Si等發(fā)展出了具有三界面、四界面和五界面的復(fù)合流動(dòng)聚焦裝置。針對(duì)近年來(lái)興起的電流動(dòng)聚焦技術(shù),考慮到真實(shí)的物理環(huán)境,本團(tuán)隊(duì)在流動(dòng)聚焦技術(shù)中引入了電場(chǎng)和熱場(chǎng)對(duì)射流不穩(wěn)定性進(jìn)行了系統(tǒng)分析。為使液滴制備更加可控,Yang等設(shè)計(jì)了通過(guò)對(duì)流動(dòng)系統(tǒng)施加外部主動(dòng)激勵(lì)來(lái)控制液滴生成的流動(dòng)聚焦裝置,Mu等探討了擾動(dòng)相、激勵(lì)頻率、振幅等多種因素對(duì)液滴生成的影響。

    從流動(dòng)聚焦中被聚焦液體的流動(dòng)過(guò)程可以看出,生成穩(wěn)定的錐形是產(chǎn)生射流并收集液滴的前提條件。由于錐形的失穩(wěn)涉及全局不穩(wěn)定性問(wèn)題,對(duì)錐形的穩(wěn)定形成條件及流場(chǎng)結(jié)構(gòu)特征開(kāi)展研究具有重要的意義。值得注意的是,錐形和射流有著尺度上的巨大差距。以單軸流動(dòng)聚焦為例,錐形尺寸從管口附近的宏觀尺度(毫米量級(jí))逐漸減小到小孔附近的微觀尺度(微米量級(jí))。此外,從針對(duì)單軸流動(dòng)聚焦開(kāi)展的大量研究中發(fā)現(xiàn),當(dāng)流動(dòng)聚焦的幾何參數(shù)產(chǎn)生變化時(shí),錐形的形貌和不穩(wěn)定性都會(huì)有較大變化。當(dāng)錐形在某些特定參數(shù)條件下處于全局不穩(wěn)定性模式時(shí),可以觀察到錐形內(nèi)部會(huì)出現(xiàn)回流區(qū),因此有必要對(duì)回流區(qū)結(jié)構(gòu)進(jìn)行流場(chǎng)顯示和細(xì)致研究。

    現(xiàn)有針對(duì)流動(dòng)聚焦的研究多重點(diǎn)關(guān)注于射流結(jié)構(gòu)的變化情況及射流不穩(wěn)定性方面,而關(guān)于流動(dòng)聚焦中錐形結(jié)構(gòu)變化以及錐形不穩(wěn)定性的研究目前開(kāi)展較少。在錐形內(nèi)部精細(xì)流場(chǎng)結(jié)構(gòu)的觀測(cè)方面開(kāi)展實(shí)驗(yàn)研究,能夠細(xì)致觀察流動(dòng)聚焦的整個(gè)過(guò)程,深刻理解錐形結(jié)構(gòu)變化以及錐形不穩(wěn)定性的內(nèi)在機(jī)理。本文針對(duì)氣驅(qū)流動(dòng)聚焦開(kāi)展被聚焦液體錐形的實(shí)驗(yàn)研究,考察流動(dòng)控制參數(shù)和結(jié)構(gòu)參數(shù)對(duì)錐形形態(tài)及穩(wěn)定性的影響;研究液體錐形內(nèi)回流區(qū)的產(chǎn)生條件并探究主要過(guò)程參數(shù)對(duì)回流區(qū)的影響,為流動(dòng)聚焦技術(shù)的進(jìn)一步應(yīng)用提供實(shí)驗(yàn)依據(jù)與理論支撐。

    1 實(shí)驗(yàn)裝置與方法

    為便于對(duì)流動(dòng)聚焦中的液體錐形進(jìn)行實(shí)時(shí)觀測(cè),本文采用吸氣式流動(dòng)聚焦裝置將錐形暴露在氣壓室腔體外面。整套實(shí)驗(yàn)裝置系統(tǒng)如圖1所示,主要由3部分組成:流動(dòng)聚焦核心裝置系統(tǒng)、參數(shù)控制系統(tǒng)和圖像采集系統(tǒng)。

    圖1 實(shí)驗(yàn)裝置示意圖Fig.1 Schematic diagram of experimental device

    1)流動(dòng)聚焦核心裝置系統(tǒng)。裝置采用吸氣式設(shè)計(jì),包括位于上方的毛細(xì)管及其固定裝置和位于下方的氣壓室腔體,氣壓室正上方薄片中心位置雕刻有聚焦小孔,正對(duì)毛細(xì)管的出口。毛細(xì)管內(nèi)流體在高速氣體的聚焦作用下穿過(guò)小孔,進(jìn)入氣壓室形成射流并破碎。吸氣式流動(dòng)聚焦裝置的優(yōu)勢(shì)在于,實(shí)驗(yàn)中毛細(xì)管的位置相對(duì)固定,氣壓室固定于一個(gè)能夠三維移動(dòng)的升降平臺(tái)上,將二者分離便于改變毛細(xì)管與小孔的距離和同軸度。吸氣式設(shè)計(jì)使得錐形暴露在氣壓室之外,比吹氣式設(shè)計(jì)更利于觀測(cè)錐形結(jié)構(gòu)。氣壓室側(cè)面設(shè)計(jì)有3個(gè)小孔,分別連接至真空抽氣泵、壓力表(由于采用吸氣方式,壓力表為負(fù)壓力表)和抽取廢液的注射器。

    2)參數(shù)控制系統(tǒng)。主要包括氣閥、真空抽氣泵、負(fù)壓力表、微量注射泵、注射器、升降臺(tái)等。氣壓室固定于升降平臺(tái)上,抽氣泵與氣閥相連,再連接至氣壓室,注射器固定于微型注射泵上,并連接至毛細(xì)管。在流動(dòng)聚焦中,影響參數(shù)可分為流動(dòng)控制參數(shù)和結(jié)構(gòu)參數(shù)2部分:流動(dòng)控制參數(shù)主要包括被聚焦液體流量(毛細(xì)管流量)Q和驅(qū)動(dòng)氣體壓差Δp。結(jié)構(gòu)參數(shù)主要包括毛細(xì)管內(nèi)徑d、小孔直徑D、管孔距離H等。本實(shí)驗(yàn)主要關(guān)注被聚焦液體流量Q、氣體壓差Δp和管孔距離H的變化對(duì)錐形形態(tài)和回流區(qū)結(jié)構(gòu)的影響。實(shí)驗(yàn)中,通過(guò)微量注射泵控制毛細(xì)管中液體流量Q,流量大小可以從注射泵上直接讀取;通過(guò)調(diào)節(jié)進(jìn)氣閥門(mén)改變氣體壓差Δp,具體數(shù)值可從壓力表上讀出;管孔距離H的控制通過(guò)調(diào)節(jié)升降平臺(tái)改變氣壓室高度來(lái)實(shí)現(xiàn)。

    3)圖像采集系統(tǒng)。在流動(dòng)聚焦中,錐形尺寸從毛細(xì)管口附近的毫米量級(jí)逐漸減小到小孔附近的微米量級(jí),同時(shí)錐形內(nèi)部流場(chǎng)流速較高,因此,拍攝錐形形態(tài)和內(nèi)部回流區(qū)結(jié)構(gòu)的圖像采集系統(tǒng)需要具有較高的空間和時(shí)間分辨率。本實(shí)驗(yàn)圖像采集系統(tǒng)包含顯微鏡(Leica Z16 apo)和高速攝影相機(jī)(Phantom VEO 710S)。高速攝影相機(jī)的拍攝幀率較高(14 000 幀/s),因此采用高功率的連續(xù)LED光源(榮峰光電RF-100W LED冷光源)。不同于液驅(qū)流動(dòng)聚焦,氣驅(qū)流動(dòng)聚焦存在氣液彎曲界面的折射,導(dǎo)致液體錐形內(nèi)部流場(chǎng)的拍攝難度較大,因此使用菲涅爾透鏡將光源發(fā)出的光聚焦至流動(dòng)聚焦的錐形所在區(qū)域。菲涅爾透鏡是由聚乙烯塑料制成的薄鏡片,一面為光滑面,另一面為粗糙面,刻有數(shù)十圈同心圓,具有聚焦與濾光作用,可以使錐形內(nèi)部圖像更加明亮清晰,以便獲得較好的圖像顯示效果。顯微鏡可將錐形區(qū)域放大,調(diào)節(jié)顯微鏡與光源、菲涅爾透鏡和流動(dòng)聚焦錐形區(qū)域,使三者中心處于同一直線上。高速攝影相機(jī)與顯微鏡相連,將拍攝的圖像傳輸至計(jì)算機(jī)進(jìn)行圖像的后處理。

    在探究參數(shù)變化對(duì)錐形形態(tài)的影響時(shí),采用控制變量法,保持其他參數(shù)不變并調(diào)節(jié)一個(gè)參數(shù)來(lái)觀察錐形形態(tài)的變化。由于錐形在穩(wěn)定形態(tài)與不穩(wěn)定形態(tài)之間轉(zhuǎn)換存在遲滯效應(yīng),錐形是否穩(wěn)定與其初始狀態(tài)參數(shù)密切相關(guān),因此在進(jìn)行實(shí)驗(yàn)時(shí),參數(shù)的變化方向選擇為錐形從穩(wěn)定形態(tài)過(guò)渡到不穩(wěn)定形態(tài)的方向,最大程度避免遲滯效應(yīng)帶來(lái)的實(shí)驗(yàn)誤差。

    實(shí)驗(yàn)時(shí),在被聚焦液體中加入示蹤粒子使流場(chǎng)可視化,在錐形穩(wěn)定的條件下,考察過(guò)程參數(shù)變化對(duì)錐形內(nèi)部回流區(qū)結(jié)構(gòu)的影響。采用控制變量法,探究某一參數(shù)影響時(shí)保持其他參數(shù)不變,使用高速攝影相機(jī)拍攝各個(gè)參數(shù)條件下的回流區(qū)結(jié)構(gòu)。必須指出,由于錐形內(nèi)的回流區(qū)是三維結(jié)構(gòu),使用高速攝影相機(jī)拍攝的回流區(qū)結(jié)構(gòu)圖像會(huì)因光學(xué)折射而產(chǎn)生變形,實(shí)驗(yàn)中應(yīng)盡量縮小拍攝景深進(jìn)行彌補(bǔ),景深ΔL的計(jì)算公式為:

    式中:f為鏡頭焦距;δ為容許彌散圓直徑,取決于相機(jī)的感光底片;L為拍攝時(shí)的對(duì)焦距離;F為光圈值。在大光圈及高放大倍率(即小對(duì)焦距離)的情況下,景深ΔL較小。在本實(shí)驗(yàn)中,拍攝景深約為40 μm,比毛細(xì)管內(nèi)徑410 μm小了一個(gè)數(shù)量級(jí)。因此,較小的景深起到了近似于在焦平面上打片光的效果。同時(shí),在實(shí)驗(yàn)拍攝中調(diào)節(jié)相機(jī)及顯微鏡位置,使得錐形邊界附近的粒子清晰可見(jiàn),確保拍攝時(shí)相機(jī)焦平面位于錐形中心軸位置。拍攝后,對(duì)高速攝影相機(jī)獲取的連續(xù)變化的實(shí)驗(yàn)圖像進(jìn)行預(yù)處理,使用基于Matlab的PIVlab軟件對(duì)圖像進(jìn)行分析,在較短時(shí)間間隔內(nèi)計(jì)算相鄰2幅圖像的相關(guān)性,得到速度場(chǎng)和流線圖等,進(jìn)而刻畫(huà)出液體錐形內(nèi)部的回流區(qū)結(jié)構(gòu)。

    實(shí)驗(yàn)的驅(qū)動(dòng)氣體為空氣,被聚焦液體為去離子水,兩者的物性參數(shù)如表1所示。

    表1 實(shí)驗(yàn)材料的物理屬性(20 ℃)Table 1 Physical properties of the experimental materials(20 ℃)

    2 實(shí)驗(yàn)結(jié)果與分析

    2.1 錐形形態(tài)及穩(wěn)定性

    在流動(dòng)聚焦中,毛細(xì)管口與小孔之間的液體錐形保持穩(wěn)定是小孔下游穩(wěn)定射流形成的前提條件。錐形的形態(tài)主要有振動(dòng)形態(tài)和穩(wěn)定形態(tài)2種。本節(jié)將探究被聚焦液體流量Q、氣體壓差Δp和管孔距離H對(duì)錐形形狀、尺寸和穩(wěn)定性的影響。

    2.1.1 量綱分析

    采用量綱分析的方法探究參數(shù)變化對(duì)錐形穩(wěn)定性的影響。在單軸流動(dòng)聚焦中,錐形的形成是小孔處驅(qū)動(dòng)氣體帶來(lái)的黏性剪切力和被聚焦液體的表面張力之間相互競(jìng)爭(zhēng)的結(jié)果。其中,黏性剪切力τ促進(jìn)界面變形,可表示為:

    式中, ρ為 氣體密度。

    與剪切力相反,被聚焦液體的表面張力τ會(huì)阻止界面變形,其在小孔位置處可表示為:

    式中,σ為表面張力系數(shù), d為射流直徑:

    式中, ρ為液體密度。

    如果錐形要保持穩(wěn)定形態(tài),剪切力必須大于表面張力。將黏性剪切力與表面張力之比表示為:

    可以看出,隨著內(nèi)外氣體壓差 Δp和被聚焦液體流量Q的增大以及管孔距離H的減小,驅(qū)動(dòng)氣體的黏性剪切力作用將變強(qiáng),最終克服表面張力并促進(jìn)錐形穩(wěn)定。

    2.1.2 實(shí)驗(yàn)結(jié)果

    為探究液體流量對(duì)錐形形態(tài)及穩(wěn)定性的影響,選取去離子水在內(nèi)外壓差Δp=10.0 kPa、管孔距離H=0.9 mm的條件下進(jìn)行實(shí)驗(yàn),其他結(jié)構(gòu)參數(shù)保持不變(小孔直徑D=0.4 mm;針管型號(hào)為22g,外徑0.71 mm,內(nèi)徑0.41 mm)。實(shí)驗(yàn)中,選取流量Q=30 mL/h為初始狀態(tài),此時(shí)錐形形態(tài)保持穩(wěn)定;然后逐漸減小流量,流量減少至Q=27 mL/h時(shí),錐形失穩(wěn),進(jìn)入振動(dòng)形態(tài),如圖2(a)所示(相鄰2張圖片時(shí)間間隔為2.2 ms);之后逐漸增大流量,圖2 (b)展示了不同流量(28、30、60、100、140 mL/h)下的錐形穩(wěn)定形態(tài)。本文圖中的標(biāo)度尺都為100 μm。

    圖2 不同液體流量下的錐形形態(tài)圖Fig.2 Cone morphology diagram at different liquid flow rates

    可見(jiàn),在其他參數(shù)條件一定時(shí),存在一個(gè)臨界流量Q:當(dāng)Q大于Q時(shí),隨著Q的增大,錐形始終保持穩(wěn)定,錐形下游在聚焦小孔處的尺寸(即射流直徑)逐漸增大,但錐形上游的形態(tài)沒(méi)有明顯變化,而隨著Q的減小,表面張力逐漸增大;當(dāng)Q小于Q后,驅(qū)動(dòng)氣體的黏性剪切力無(wú)法克服聚焦液體的表面張力,此時(shí)錐形失穩(wěn),進(jìn)入振動(dòng)形態(tài)。值得注意的是,當(dāng)流動(dòng)聚焦裝置的其他參數(shù)改變時(shí),臨界流量Q也會(huì)隨之改變。

    為探究?jī)?nèi)外氣體壓差對(duì)錐形形態(tài)及穩(wěn)定性的影響,在Q=30 mL/h、H=0.9 mm的條件下進(jìn)行實(shí)驗(yàn)。與前文所述方法類(lèi)似,從Δp=2.0 kPa時(shí)錐穩(wěn)定形態(tài)開(kāi)始減小Δp,Δp=1.4 kPa時(shí),錐形開(kāi)始振動(dòng)失穩(wěn),如圖3 (a)所示(相鄰2張圖片時(shí)間間隔為2.1 ms)??梢?jiàn),在一定的被聚焦流體流量(大于Q)下,存在一個(gè)臨界內(nèi)外氣體壓差Δp,當(dāng)內(nèi)外壓差小于臨界值時(shí),錐形失穩(wěn)。之后逐漸增大Δp,得到如圖3 (b)所示的錐形穩(wěn)定形態(tài)(Δp分別為1.5、5.0、10.0、20.0和30.0 kPa)??梢园l(fā)現(xiàn),隨著Δp逐漸增大,射流直徑逐漸減小。進(jìn)一步的實(shí)驗(yàn)發(fā)現(xiàn),當(dāng)Q很小時(shí)(小于Q),無(wú)論Δp如何增大,錐形都始終處于振動(dòng)形態(tài),由此可以看出,足夠大的被聚焦液體流量是維持錐形結(jié)構(gòu)穩(wěn)定的前提條件。

    圖3 不同氣體壓力差下的錐形形態(tài)圖Fig.3 Cone morphology diagram at different gas pressure differences

    為探究管孔距離H對(duì)錐形形態(tài)及穩(wěn)定性的影響,在Q=30 mL/h、Δp=10.0 kPa的條件下進(jìn)行實(shí)驗(yàn)。管孔距離H從0.5 mm開(kāi)始,依次增加至 0.7、0.9和1.0 mm,如圖4 (a)所示??梢钥吹?,在上述結(jié)構(gòu)參數(shù)和控制參數(shù)下,錐形都為穩(wěn)定狀態(tài)。當(dāng)H較小時(shí),從毛細(xì)管口流出的液體呈匯聚狀,這是因?yàn)槊?xì)管口與聚焦小孔距離較近,驅(qū)動(dòng)氣流在錐形上下游表面均有較高速度,剪切拉伸作用較強(qiáng),使得錐形呈匯聚狀;隨著H的增大,毛細(xì)管口附近的氣體流速逐漸減小,剪切作用減小,錐形上游部分開(kāi)始向外擴(kuò)張并逐漸隆起,直至直徑與毛細(xì)管內(nèi)徑相當(dāng),而錐形在小孔附近的下游部分仍因驅(qū)動(dòng)氣流速度較快而被剪切拉伸呈匯聚狀,且射流直徑不變。當(dāng)管孔距離增大至1.1 mm時(shí),錐形已進(jìn)入振動(dòng)形態(tài),如圖4 (b)所示(相鄰2張圖片時(shí)間間隔為2.2 ms),此時(shí),毛細(xì)管口距小孔過(guò)遠(yuǎn),導(dǎo)致管口附近的氣流速度非常小,無(wú)法對(duì)被聚焦液體提供足夠的剪切力,使得錐形失穩(wěn),進(jìn)入振動(dòng)形態(tài)。

    圖4 不同管孔距離下的錐形形態(tài)圖Fig.4 Cone morphology diagram at different heights between the tube and the hole

    2.2 回流區(qū)結(jié)構(gòu)

    在氣驅(qū)流動(dòng)聚焦中,當(dāng)液體錐形處在穩(wěn)定形態(tài)時(shí),氣液界面內(nèi)部流場(chǎng)結(jié)構(gòu)會(huì)因界面上氣流剪切力分布不均勻而隨參數(shù)發(fā)生變化,錐形內(nèi)部的軸向速度分布直接受剪切力影響。在一定條件下,當(dāng)速度梯度變化不劇烈時(shí),能夠保證被聚焦液體都向下游流動(dòng),流線呈均勻分布。但在某些條件下,錐形內(nèi)部會(huì)產(chǎn)生回流區(qū),有效流動(dòng)截面積(沿水平方向垂直于錐形對(duì)稱(chēng)軸軸線的液體截面)減小,界面附近的液體速度梯度增大,從而維持界面上的剪切力平衡。為使錐形內(nèi)部流場(chǎng)可視化,在被聚焦液體中加入直徑約20 μm的玻璃微球作為示蹤粒子。用高速攝影相機(jī)采集圖像后,使用基于Matlab的PIVlab軟件對(duì)圖像進(jìn)行分析,得到速度場(chǎng)、流線圖和渦量場(chǎng)等。如圖5所示,在Δp=10.0 kPa、H=0.9 mm、Q=30 mL/h時(shí),可以在錐形中清晰地觀察到回流區(qū)的存在。

    圖5 回流區(qū)結(jié)構(gòu)圖Fig.5 Recirculation cell structure

    由于光在氣液界面折射引發(fā)的透鏡效應(yīng),實(shí)驗(yàn)拍攝的回流區(qū)存在一定的變形,計(jì)算的渦量場(chǎng)并不準(zhǔn)確。因此,本文主要定性考察不同參數(shù)下錐形內(nèi)部回流區(qū)的結(jié)構(gòu)和存在情況。定量研究需要關(guān)注多尺度下錐形內(nèi)部流場(chǎng)顯示的光路問(wèn)題,這是實(shí)驗(yàn)研究的難點(diǎn),留待進(jìn)一步解決。

    在Δp=10.0 kPa、H=0.9 mm條件下,實(shí)驗(yàn)探究了液體流量因素對(duì)回流區(qū)結(jié)構(gòu)的影響。在此條件下,錐形穩(wěn)定的臨界流量Q=28 mL/h。將毛細(xì)管中液體流量從28 mL/h開(kāi)始逐漸增大,觀察錐形內(nèi)部有無(wú)回流區(qū)出現(xiàn),用高速攝影相機(jī)拍攝流量變化時(shí)回流區(qū)結(jié)構(gòu)的變化。結(jié)果發(fā)現(xiàn)Q=Q、錐形穩(wěn)定時(shí)內(nèi)部即存在回流區(qū),如圖6 (a)所示?;亓鲄^(qū)的不對(duì)稱(chēng)也證明了由于光的折射而產(chǎn)生的圖像變形。進(jìn)一步用高速攝影相機(jī)拍攝了Q=30 mL/h與Q=40 mL/h時(shí)的回流區(qū)圖像,其流線圖如圖6 (b)、(c)所示。從圖中可以看出,在該組參數(shù)條件下,回流區(qū)結(jié)構(gòu)在錐形穩(wěn)定的臨界條件Q=Q時(shí)就已存在,且回流區(qū)大小隨Q的增大而減小。當(dāng)Q大于特定值Q=42 mL/h時(shí),回流區(qū)消失,如圖6(d)所示。

    圖6 不同毛細(xì)管流量的回流區(qū)結(jié)構(gòu)Fig.6 Recirculation cell structure at different liquid flow rates

    回流區(qū)的產(chǎn)生是錐形界面兩側(cè)剪切力平衡的結(jié)果,毛細(xì)管流量減小時(shí),被聚焦液體的平均速度隨之降低。因此,回流區(qū)尺寸需要增大以減小錐形內(nèi)被聚焦液體的有效流動(dòng)截面積,從而使界面內(nèi)側(cè)液體的速度梯度增大,以達(dá)到界面兩側(cè)剪切力平衡。而毛細(xì)管流量較大時(shí),錐形內(nèi)部液體的平均流速增大,回流區(qū)會(huì)逐漸減小甚至消失。

    圖7 不同氣體壓差下的回流區(qū)結(jié)構(gòu)Fig.7 Recirculation cell structure at different gas pressure differences

    為探究管孔距離對(duì)回流區(qū)結(jié)構(gòu)的影響,在Q=30 mL/h、Δp=10.0 kPa的條件下進(jìn)行實(shí)驗(yàn)。已知在此條件下,H=1.1 mm時(shí)錐形進(jìn)入振動(dòng)形態(tài),因此本組實(shí)驗(yàn)中管孔距離分別設(shè)定為0.7 、0.8、0.9和1.0 mm,處理得到的錐形內(nèi)部流線圖如圖8所示??梢园l(fā)現(xiàn),H=0.7~0.9 mm時(shí),回流區(qū)隨著H的增大而增大,而H=1.0 mm時(shí)的回流區(qū)卻略小于H=0.9 mm時(shí)的回流區(qū)。其原因是:管孔距離增大時(shí),錐形尺寸逐漸變大,其有效流動(dòng)截面積也變大,因此回流區(qū)增大,使錐形內(nèi)界面附近液體達(dá)到界面剪切力平衡所需要的速度;管孔距離較大時(shí),繼續(xù)增大管孔距離,有效流動(dòng)截面積的增大會(huì)減慢,且錐形上游的氣體流速也會(huì)因距小孔太遠(yuǎn)而減小,導(dǎo)致剪切作用減弱,進(jìn)而導(dǎo)致回流區(qū)變小且只存在于錐形下游。

    圖8 不同管孔距離下的回流區(qū)結(jié)構(gòu)Fig.8 Recirculation cell structure at different heights between the tube and the hole

    3 結(jié) 論

    本文采用控制變量法在吸氣式流動(dòng)聚焦平臺(tái)上定性研究了不同過(guò)程參數(shù)對(duì)錐形動(dòng)力學(xué)行為以及錐形內(nèi)部流場(chǎng)結(jié)構(gòu)的影響規(guī)律。結(jié)果表明,增大被聚焦液體流量、驅(qū)動(dòng)氣體壓差或減小管孔距離,都能提高錐形的穩(wěn)定性,顯著改變錐形的形態(tài)。此外,在特定的參數(shù)條件下,錐形內(nèi)部會(huì)產(chǎn)生回流區(qū),其尺寸隨著聚焦液體流量的減小或驅(qū)動(dòng)氣體壓差的增大而增大。管孔距離的變化對(duì)回流區(qū)大小的影響則是非單調(diào)的。雖然粒子示蹤方法為研究錐形的內(nèi)部流場(chǎng)結(jié)構(gòu)打開(kāi)了一扇窗,但光學(xué)折射導(dǎo)致的多尺度錐形內(nèi)部流場(chǎng)變形影響了實(shí)驗(yàn)結(jié)果的精度,因此實(shí)驗(yàn)方法仍需改進(jìn),這也為下一步研究工作指明了方向。

    嫩草影视91久久| 国产一级毛片七仙女欲春2| 欧美日韩瑟瑟在线播放| 成人永久免费在线观看视频| 国产毛片a区久久久久| 欧美三级亚洲精品| 亚洲精品亚洲一区二区| 丰满的人妻完整版| 嫩草影院精品99| or卡值多少钱| 久久人人精品亚洲av| 色综合亚洲欧美另类图片| 国产不卡一卡二| 亚洲内射少妇av| 午夜福利在线观看吧| 国产黄色小视频在线观看| 国产一区二区在线av高清观看| 精品免费久久久久久久清纯| 国产一区二区在线av高清观看| 亚洲av熟女| 欧美成狂野欧美在线观看| 午夜激情欧美在线| 国产色爽女视频免费观看| 欧美最黄视频在线播放免费| 亚洲七黄色美女视频| 色综合亚洲欧美另类图片| 搡老岳熟女国产| 美女免费视频网站| 精华霜和精华液先用哪个| 免费看a级黄色片| 欧美又色又爽又黄视频| 亚洲国产精品sss在线观看| 中文字幕人成人乱码亚洲影| 亚洲美女黄片视频| 999久久久精品免费观看国产| 免费av不卡在线播放| 成人欧美大片| 少妇的逼水好多| 午夜日韩欧美国产| 国产在线精品亚洲第一网站| 直男gayav资源| 日本一本二区三区精品| 人人妻人人看人人澡| 国产成+人综合+亚洲专区| 好看av亚洲va欧美ⅴa在| 麻豆国产av国片精品| 两人在一起打扑克的视频| 欧洲精品卡2卡3卡4卡5卡区| 欧美乱妇无乱码| 国产日本99.免费观看| 成年免费大片在线观看| 国产av麻豆久久久久久久| 日日摸夜夜添夜夜添小说| 日韩欧美在线乱码| 成人一区二区视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产69精品久久久久777片| 欧美黄色淫秽网站| 能在线免费观看的黄片| 免费在线观看成人毛片| 国内毛片毛片毛片毛片毛片| 日韩欧美国产一区二区入口| 亚洲精品亚洲一区二区| 亚洲自偷自拍三级| 51午夜福利影视在线观看| 99riav亚洲国产免费| 三级毛片av免费| 国产美女午夜福利| 五月伊人婷婷丁香| 国产黄a三级三级三级人| 免费搜索国产男女视频| .国产精品久久| 好看av亚洲va欧美ⅴa在| 女同久久另类99精品国产91| 亚洲经典国产精华液单 | netflix在线观看网站| 欧美最新免费一区二区三区 | 男女做爰动态图高潮gif福利片| www日本黄色视频网| 自拍偷自拍亚洲精品老妇| 热99在线观看视频| 国产精品自产拍在线观看55亚洲| 精品午夜福利在线看| 亚洲av日韩精品久久久久久密| 69av精品久久久久久| 久久性视频一级片| 中文字幕精品亚洲无线码一区| 精品人妻一区二区三区麻豆 | 好看av亚洲va欧美ⅴa在| 少妇的逼水好多| 桃色一区二区三区在线观看| 精品一区二区三区av网在线观看| 亚洲av五月六月丁香网| 成人三级黄色视频| 有码 亚洲区| 成人无遮挡网站| 中文字幕免费在线视频6| 国产色婷婷99| 国产伦一二天堂av在线观看| 成年女人看的毛片在线观看| 一区福利在线观看| 国产精品女同一区二区软件 | 亚洲综合色惰| 香蕉av资源在线| 天天躁日日操中文字幕| 男人狂女人下面高潮的视频| netflix在线观看网站| 1000部很黄的大片| 波多野结衣高清作品| 免费搜索国产男女视频| 非洲黑人性xxxx精品又粗又长| 午夜久久久久精精品| 在线a可以看的网站| 亚洲国产精品999在线| 国产成人福利小说| 最近中文字幕高清免费大全6 | 日韩欧美一区二区三区在线观看| 日韩欧美三级三区| 日日干狠狠操夜夜爽| 99精品久久久久人妻精品| 久久久国产成人精品二区| 国产色婷婷99| 欧美日韩综合久久久久久 | 免费观看人在逋| 久久久久九九精品影院| 亚洲久久久久久中文字幕| 国产精品不卡视频一区二区 | 亚洲精品影视一区二区三区av| 精品无人区乱码1区二区| 日本一本二区三区精品| 国产乱人伦免费视频| 日韩精品青青久久久久久| 免费av毛片视频| 十八禁人妻一区二区| 老司机午夜十八禁免费视频| 看黄色毛片网站| 日韩精品青青久久久久久| 我要看日韩黄色一级片| 能在线免费观看的黄片| 国产精品电影一区二区三区| 国产精品一区二区性色av| 美女黄网站色视频| 亚洲国产精品成人综合色| 国产单亲对白刺激| 99热精品在线国产| 国产乱人伦免费视频| av欧美777| 国产三级黄色录像| 两个人视频免费观看高清| 欧美丝袜亚洲另类 | 97碰自拍视频| 久久人人爽人人爽人人片va | 一级av片app| 欧洲精品卡2卡3卡4卡5卡区| netflix在线观看网站| 99视频精品全部免费 在线| 99热精品在线国产| 亚洲久久久久久中文字幕| 看片在线看免费视频| 亚洲,欧美,日韩| 色噜噜av男人的天堂激情| 精品乱码久久久久久99久播| 一个人免费在线观看电影| 午夜福利视频1000在线观看| 人妻久久中文字幕网| 小蜜桃在线观看免费完整版高清| 午夜福利在线观看免费完整高清在 | 禁无遮挡网站| 观看免费一级毛片| 精品福利观看| 欧美又色又爽又黄视频| 国产伦精品一区二区三区视频9| 亚洲av美国av| 毛片一级片免费看久久久久 | 久久久久国内视频| 免费黄网站久久成人精品 | 长腿黑丝高跟| 俺也久久电影网| 色尼玛亚洲综合影院| 婷婷精品国产亚洲av| 欧美在线黄色| 观看美女的网站| 首页视频小说图片口味搜索| 国产精品人妻久久久久久| 一区二区三区高清视频在线| АⅤ资源中文在线天堂| 给我免费播放毛片高清在线观看| 国产综合懂色| 成年女人看的毛片在线观看| 国产私拍福利视频在线观看| 欧美xxxx黑人xx丫x性爽| 日本三级黄在线观看| 欧美最新免费一区二区三区 | 十八禁人妻一区二区| 天堂影院成人在线观看| 久久国产精品影院| 亚洲激情在线av| 中亚洲国语对白在线视频| 国产成人a区在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲最大成人手机在线| 亚洲中文日韩欧美视频| 亚洲,欧美精品.| 在线看三级毛片| 又紧又爽又黄一区二区| 日本a在线网址| 国产在线精品亚洲第一网站| 深夜a级毛片| 青草久久国产| 精品福利观看| 国产伦人伦偷精品视频| av专区在线播放| 色噜噜av男人的天堂激情| 国产真实乱freesex| 在线十欧美十亚洲十日本专区| 国产精品久久久久久精品电影| 国产中年淑女户外野战色| av福利片在线观看| 日本免费a在线| 又黄又爽又刺激的免费视频.| 最好的美女福利视频网| 99在线视频只有这里精品首页| 国产精品国产高清国产av| 亚洲内射少妇av| 精品日产1卡2卡| 免费在线观看影片大全网站| 国产精品久久视频播放| 99国产精品一区二区三区| 哪里可以看免费的av片| 中文在线观看免费www的网站| 男人舔奶头视频| 俺也久久电影网| 欧美绝顶高潮抽搐喷水| 综合色av麻豆| 亚洲经典国产精华液单 | 色5月婷婷丁香| 国产精品,欧美在线| 久久99热6这里只有精品| 国产亚洲精品久久久久久毛片| 俄罗斯特黄特色一大片| 精品久久久久久久人妻蜜臀av| 他把我摸到了高潮在线观看| 精品99又大又爽又粗少妇毛片 | 舔av片在线| 中文字幕人妻熟人妻熟丝袜美| 日本五十路高清| 老鸭窝网址在线观看| 99热这里只有是精品50| 一区二区三区高清视频在线| 91麻豆精品激情在线观看国产| 91在线观看av| 成人毛片a级毛片在线播放| 色综合站精品国产| 在线看三级毛片| 欧美+日韩+精品| 国产精品影院久久| 偷拍熟女少妇极品色| 国产精品久久久久久久久免 | 欧美又色又爽又黄视频| 深爱激情五月婷婷| 人妻久久中文字幕网| 99视频精品全部免费 在线| 精品午夜福利在线看| 老司机午夜十八禁免费视频| 婷婷亚洲欧美| 五月玫瑰六月丁香| 美女高潮的动态| av天堂中文字幕网| 亚洲激情在线av| 久久99热6这里只有精品| 亚洲 欧美 日韩 在线 免费| 有码 亚洲区| 夜夜夜夜夜久久久久| 人妻制服诱惑在线中文字幕| 一二三四社区在线视频社区8| xxxwww97欧美| 看十八女毛片水多多多| 国产精品三级大全| 久久精品国产99精品国产亚洲性色| 黄色日韩在线| 久久久久久大精品| 久久久久精品国产欧美久久久| 最好的美女福利视频网| 中亚洲国语对白在线视频| 日韩欧美在线乱码| 香蕉av资源在线| 欧美丝袜亚洲另类 | 每晚都被弄得嗷嗷叫到高潮| 亚洲自偷自拍三级| 亚洲成a人片在线一区二区| 女同久久另类99精品国产91| 久久久久久久久久黄片| 老司机午夜十八禁免费视频| 国产乱人伦免费视频| 久久精品综合一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 身体一侧抽搐| 天美传媒精品一区二区| 别揉我奶头 嗯啊视频| 亚洲精品亚洲一区二区| 97热精品久久久久久| 少妇的逼水好多| 村上凉子中文字幕在线| 亚洲精品影视一区二区三区av| 久久久久久久久中文| avwww免费| 好男人在线观看高清免费视频| 真人做人爱边吃奶动态| 国内揄拍国产精品人妻在线| 久9热在线精品视频| 伊人久久精品亚洲午夜| 亚洲18禁久久av| 国产视频内射| 99视频精品全部免费 在线| 深夜精品福利| 国产人妻一区二区三区在| 国语自产精品视频在线第100页| 亚洲国产欧洲综合997久久,| 女人十人毛片免费观看3o分钟| 一卡2卡三卡四卡精品乱码亚洲| 在线观看免费视频日本深夜| 我的女老师完整版在线观看| 久久久久九九精品影院| 日日夜夜操网爽| 色在线成人网| 观看美女的网站| 女人十人毛片免费观看3o分钟| av黄色大香蕉| 国产真实伦视频高清在线观看 | 男女床上黄色一级片免费看| 成年女人毛片免费观看观看9| 亚洲最大成人中文| a级毛片免费高清观看在线播放| 中出人妻视频一区二区| 国产免费一级a男人的天堂| 国产精品久久久久久人妻精品电影| 99热精品在线国产| 亚洲国产精品合色在线| 精品久久久久久,| 国产欧美日韩一区二区精品| 免费无遮挡裸体视频| 国产蜜桃级精品一区二区三区| 99久久精品国产亚洲精品| 国产亚洲精品久久久久久毛片| 丰满乱子伦码专区| 老司机午夜福利在线观看视频| 国产精品98久久久久久宅男小说| 有码 亚洲区| 如何舔出高潮| 99久久久亚洲精品蜜臀av| 日日干狠狠操夜夜爽| 亚洲欧美日韩高清在线视频| 精品人妻视频免费看| 一本一本综合久久| 国产乱人伦免费视频| 欧美性猛交╳xxx乱大交人| 日韩欧美在线二视频| 男人舔奶头视频| 午夜福利视频1000在线观看| 欧美+日韩+精品| АⅤ资源中文在线天堂| 欧美另类亚洲清纯唯美| 精品人妻偷拍中文字幕| 色精品久久人妻99蜜桃| 午夜日韩欧美国产| 久久婷婷人人爽人人干人人爱| 国内精品久久久久久久电影| 日韩中文字幕欧美一区二区| 免费大片18禁| 久久人妻av系列| 精品久久久久久,| 一级作爱视频免费观看| 波野结衣二区三区在线| 久久精品影院6| 国产在线精品亚洲第一网站| 亚洲内射少妇av| 丁香六月欧美| 国产一区二区三区在线臀色熟女| 特级一级黄色大片| 亚洲欧美日韩卡通动漫| 亚洲一区二区三区不卡视频| 国产伦人伦偷精品视频| 看十八女毛片水多多多| 亚洲欧美清纯卡通| 国产高清视频在线播放一区| 热99re8久久精品国产| 国产黄片美女视频| 国产精品爽爽va在线观看网站| 天堂网av新在线| 色噜噜av男人的天堂激情| 国产人妻一区二区三区在| 亚洲,欧美,日韩| 国产成年人精品一区二区| 1024手机看黄色片| 在线观看66精品国产| 亚洲av熟女| 国产精品电影一区二区三区| 最新中文字幕久久久久| 无人区码免费观看不卡| 国产亚洲精品久久久久久毛片| 69人妻影院| 可以在线观看毛片的网站| 久久久久久久久大av| 亚洲一区高清亚洲精品| 麻豆久久精品国产亚洲av| 啪啪无遮挡十八禁网站| 国产亚洲精品av在线| 一二三四社区在线视频社区8| 亚洲精品乱码久久久v下载方式| 国产野战对白在线观看| 狂野欧美白嫩少妇大欣赏| 欧美三级亚洲精品| 久久午夜福利片| 亚洲国产精品合色在线| 久久久精品欧美日韩精品| 精品久久久久久久久久免费视频| 亚洲精品在线美女| 18禁在线播放成人免费| 亚洲三级黄色毛片| 欧美+亚洲+日韩+国产| 欧美在线一区亚洲| 9191精品国产免费久久| 99在线人妻在线中文字幕| 最近最新中文字幕大全电影3| 午夜福利欧美成人| 成人鲁丝片一二三区免费| 好男人电影高清在线观看| 久久九九热精品免费| 国产黄片美女视频| 国产黄色小视频在线观看| 在线国产一区二区在线| 天堂动漫精品| 九色国产91popny在线| www.色视频.com| 亚洲男人的天堂狠狠| 成人av一区二区三区在线看| 午夜亚洲福利在线播放| 国产精品电影一区二区三区| 一进一出好大好爽视频| 免费看日本二区| 嫩草影视91久久| 亚洲欧美日韩高清专用| 一级黄色大片毛片| 亚洲黑人精品在线| 伊人久久精品亚洲午夜| 黄色视频,在线免费观看| 丁香六月欧美| 黄色女人牲交| 色综合欧美亚洲国产小说| 色在线成人网| xxxwww97欧美| 午夜福利在线在线| 亚洲人成电影免费在线| 永久网站在线| 99riav亚洲国产免费| 国产白丝娇喘喷水9色精品| 国产精品女同一区二区软件 | 亚洲人成网站高清观看| 国产精品一区二区三区四区免费观看 | 69av精品久久久久久| 中文字幕av在线有码专区| 变态另类丝袜制服| 99精品在免费线老司机午夜| 99热只有精品国产| 级片在线观看| 精品一区二区免费观看| 亚洲国产色片| 窝窝影院91人妻| 久久精品夜夜夜夜夜久久蜜豆| 亚洲男人的天堂狠狠| 97人妻精品一区二区三区麻豆| 国产精品伦人一区二区| 精品久久久久久,| 夜夜夜夜夜久久久久| 久久国产乱子伦精品免费另类| 国产亚洲欧美98| 又粗又爽又猛毛片免费看| 欧美xxxx性猛交bbbb| 日韩精品青青久久久久久| 亚洲人成伊人成综合网2020| 国产午夜精品久久久久久一区二区三区 | 99热这里只有精品一区| 国产高清激情床上av| 亚洲第一电影网av| 宅男免费午夜| 国产成人a区在线观看| 亚洲熟妇熟女久久| 窝窝影院91人妻| 欧美一区二区国产精品久久精品| 99久久精品一区二区三区| 舔av片在线| 黄色配什么色好看| 亚洲aⅴ乱码一区二区在线播放| 天天一区二区日本电影三级| 午夜福利在线观看吧| 国产亚洲精品久久久久久毛片| 欧美成人一区二区免费高清观看| av天堂中文字幕网| 一区二区三区免费毛片| 亚洲电影在线观看av| 一级作爱视频免费观看| 99国产极品粉嫩在线观看| 免费在线观看成人毛片| 亚洲乱码一区二区免费版| 亚洲一区高清亚洲精品| av国产免费在线观看| 亚洲avbb在线观看| 久久99热这里只有精品18| 国产av一区在线观看免费| 香蕉av资源在线| 成人美女网站在线观看视频| 国产国拍精品亚洲av在线观看| 成人一区二区视频在线观看| 九九热线精品视视频播放| 日日摸夜夜添夜夜添av毛片 | www.www免费av| 永久网站在线| 可以在线观看的亚洲视频| 欧美+亚洲+日韩+国产| 国产视频内射| 97人妻精品一区二区三区麻豆| 一卡2卡三卡四卡精品乱码亚洲| 国产欧美日韩一区二区三| 直男gayav资源| 国产精品久久久久久精品电影| 国产大屁股一区二区在线视频| 天堂动漫精品| 欧美潮喷喷水| 可以在线观看的亚洲视频| 嫩草影视91久久| 在线免费观看不下载黄p国产 | 亚洲美女视频黄频| 日本三级黄在线观看| 国产免费av片在线观看野外av| 99久久无色码亚洲精品果冻| 国产单亲对白刺激| 搡老岳熟女国产| 色av中文字幕| 精品国产三级普通话版| 天堂影院成人在线观看| 久久中文看片网| 最新在线观看一区二区三区| av中文乱码字幕在线| 国内精品久久久久久久电影| 国产精品影院久久| 日本黄色片子视频| 国产真实伦视频高清在线观看 | 亚洲欧美激情综合另类| 日韩欧美精品v在线| 夜夜夜夜夜久久久久| 亚洲自拍偷在线| 极品教师在线视频| 久久精品国产自在天天线| 免费黄网站久久成人精品 | 午夜福利成人在线免费观看| 搡女人真爽免费视频火全软件 | 久久久成人免费电影| 久久久精品大字幕| 有码 亚洲区| 好看av亚洲va欧美ⅴa在| 91久久精品国产一区二区成人| 国产老妇女一区| 久99久视频精品免费| 国产野战对白在线观看| 蜜桃久久精品国产亚洲av| 美女黄网站色视频| 国产av麻豆久久久久久久| 听说在线观看完整版免费高清| 亚洲欧美日韩东京热| 久久久久国内视频| 国产高潮美女av| 欧美色欧美亚洲另类二区| 亚洲国产精品999在线| 色尼玛亚洲综合影院| 国产欧美日韩精品亚洲av| 高潮久久久久久久久久久不卡| 他把我摸到了高潮在线观看| 亚洲精品乱码久久久v下载方式| 美女高潮喷水抽搐中文字幕| 精品久久久久久久人妻蜜臀av| 久久久精品欧美日韩精品| 人妻制服诱惑在线中文字幕| 一级毛片久久久久久久久女| 波多野结衣巨乳人妻| or卡值多少钱| 又黄又爽又免费观看的视频| 我要搜黄色片| 国产亚洲精品av在线| 国产黄色小视频在线观看| 精品福利观看| bbb黄色大片| 国产精品,欧美在线| 少妇熟女aⅴ在线视频| 无人区码免费观看不卡| 在线免费观看的www视频| 国内精品一区二区在线观看| 无人区码免费观看不卡| 亚洲av美国av| 亚洲av五月六月丁香网| xxxwww97欧美| АⅤ资源中文在线天堂| 在线观看66精品国产| 一级av片app| 久久久久久久午夜电影| 欧美精品啪啪一区二区三区| 午夜免费激情av| 亚洲 国产 在线| 欧美成人免费av一区二区三区| 久久精品国产亚洲av香蕉五月| 成人美女网站在线观看视频| 综合色av麻豆| 欧美+日韩+精品| 久久久久国内视频| 久久欧美精品欧美久久欧美| 国产成人啪精品午夜网站|