久久高清超碰?V热热久久|久久久综合精品|国产精品日韩专区欧美专区|亚洲精品另类欧美|欧美精品AⅤ一区二区三区|绯色av一区二区三区亚洲人妻|久久久久人妻精品一区|亚洲中文久久精品无码软件

芬蘭Kibron專注表面張力儀測(cè)量技術(shù),快速精準(zhǔn)測(cè)量動(dòng)靜態(tài)表面張力

熱線:021-66110810,66110819,66110690,13564362870 Email: info@vizai.cn

合作客戶/

拜耳公司.jpg

拜耳公司

同濟(jì)大學(xué)

同濟(jì)大學(xué)

聯(lián)合大學(xué).jpg

聯(lián)合大學(xué)

寶潔公司

美國(guó)保潔

強(qiáng)生=

美國(guó)強(qiáng)生

瑞士羅氏

瑞士羅氏

當(dāng)前位置:首頁(yè) > 新聞中心

表面張力變化對(duì)含氣泡液體射流破裂的影響

來(lái)源:化工學(xué)報(bào) 瀏覽 2478 次 發(fā)布時(shí)間:2022-10-11

表面張力變化對(duì)含氣泡液體射流破裂的影響


吳兆偉 施浙杭 趙輝 周騖 蔡小舒 劉海峰


摘要

使用高速相機(jī)研究了表面張力變化對(duì)含氣泡液體射流破裂過(guò)程的影響。通過(guò)改變表面活性劑濃度獲得了不同表面張力的液體射流。實(shí)驗(yàn)發(fā)現(xiàn)當(dāng)液體射流速度保持不變時(shí),減小液體表面張力會(huì)增加射流破裂長(zhǎng)度。表面活性劑一方面降低了液體動(dòng)態(tài)表面張力,減小了射流表面不穩(wěn)定波的增長(zhǎng)率,增大了射流破裂長(zhǎng)度;另一方面表面活性劑在射流表面的非均勻分布會(huì)產(chǎn)生Marangoni應(yīng)力,促使液體向射流變形區(qū)運(yùn)動(dòng),從而推遲了射流破裂的發(fā)生,增大了射流破裂長(zhǎng)度。通過(guò)理論分析得到了液體射流破裂長(zhǎng)度表達(dá)式。發(fā)現(xiàn)射流內(nèi)部氣泡會(huì)顯著縮短含表面活性劑射流的破裂長(zhǎng)度。通過(guò)氣泡擾動(dòng)射流速度和吸附表面活性劑的分析,揭示了內(nèi)部氣泡對(duì)含表面活性劑射流破裂的影響規(guī)律。


引言

液體射流常見(jiàn)于日常生活以及工業(yè)生產(chǎn)領(lǐng)域,如化工生產(chǎn)、航空航天及交通運(yùn)輸?shù)萚1-3]。在煤化工中,多噴嘴對(duì)置式水煤漿氣化技術(shù)由于其易于大型化以及先進(jìn)的工藝指標(biāo)等特點(diǎn)得到了廣泛的應(yīng)用[4]。在多噴嘴對(duì)置式水煤漿氣化技術(shù)中,混合含表面活性劑廢水的水煤漿噴射進(jìn)入氣化爐,在環(huán)隙氣流的作用下發(fā)生破裂并霧化形成大量小液滴。液體射流破裂過(guò)程是液體霧化的基礎(chǔ),研究含表面活性劑射流的破裂過(guò)程具有重要意義[5-7]。

自Rayleigh[8]和Weber[9]使用線性不穩(wěn)定分析理論研究液體射流破裂過(guò)程以來(lái),關(guān)于液體射流的研究層出不窮[10-13]。液體射流破裂過(guò)程不僅受射流速度和液體黏度的影響,還受射流液體表面張力的影響[14-17]。Reitz[18]發(fā)現(xiàn),在射流破裂Rayleigh模式下,牛頓流體射流不穩(wěn)定波的最大增長(zhǎng)率可以表示為表面張力的函數(shù),提高液體表面張力可以促進(jìn)射流破裂的發(fā)生。Chang等[19]研究了冪律流體射流不穩(wěn)定性問(wèn)題,發(fā)現(xiàn)在Rayleigh模式下,表面張力會(huì)促進(jìn)液體射流的破裂,而液體黏度則會(huì)阻礙射流的破裂。Martínez-Calvo等[20]使用數(shù)值分析方法研究了表面吸附有不可溶表面活性劑的液絲的破裂過(guò)程,得到了不同工況下衛(wèi)星液滴體積的表達(dá)式。Hameed等[21]發(fā)現(xiàn)不可溶表面活性劑的存在能夠延緩液絲的破裂。Timmermans等[22]研究了表面吸附有不可溶表面活性劑的液絲的穩(wěn)定性并指出表面活性劑主要通過(guò)改變液絲毛細(xì)壓力以及液絲表面剪切應(yīng)力影響液絲的穩(wěn)定性。如何獲得適合的射流破裂長(zhǎng)度一直是射流破裂研究的熱點(diǎn)。Lad等[23]發(fā)現(xiàn)可以通過(guò)對(duì)液體射流施加額外的擾動(dòng)實(shí)現(xiàn)更小的射流破裂長(zhǎng)度。Wu等[24]發(fā)現(xiàn)在液體射流內(nèi)部注入氣泡能夠顯著縮短射流破裂長(zhǎng)度。關(guān)于射流內(nèi)部氣泡對(duì)含表面活性劑射流破裂過(guò)程影響的研究則少見(jiàn)報(bào)道。

前人在研究液體性質(zhì)對(duì)射流破裂過(guò)程的影響時(shí),多采用Ohnesorge數(shù)Oh和Weber數(shù)We等以描述黏度相對(duì)表面張力的大小或氣動(dòng)力相對(duì)表面張力的大小。較少的實(shí)驗(yàn)研究單獨(dú)關(guān)注液體表面張力對(duì)射流破裂過(guò)程的影響。這是由于當(dāng)采用不同表面張力的流體作為實(shí)驗(yàn)流體時(shí),往往流體的其他性質(zhì),例如黏度或密度也會(huì)不同。表面活性劑是一類只要少量加入即能顯著改變液體表面張力的物質(zhì)。這提供了一個(gè)極好的獨(dú)立考察表面張力影響液體射流破裂過(guò)程的視角。前人關(guān)于表面活性劑對(duì)液體射流或液絲破裂的研究多數(shù)基于數(shù)值模擬或理論推導(dǎo),關(guān)于含表面活性劑射流或液絲破裂過(guò)程的實(shí)驗(yàn)研究較少[20,25-26]。本文采用十二烷基苯磺酸鈉(SDBS)作為表面活性劑,研究了表面張力變化對(duì)液體射流破裂過(guò)程及射流內(nèi)部氣泡對(duì)射流破裂長(zhǎng)度的影響。對(duì)不同表面張力液體射流的破裂過(guò)程展開(kāi)研究有利于加深對(duì)射流破裂過(guò)程的理解,為含表面活性劑射流以及含氣泡射流的工業(yè)應(yīng)用提供理論指導(dǎo)。


1、實(shí)驗(yàn)流程與介質(zhì)

實(shí)驗(yàn)流程如圖1所示。來(lái)自注射泵的氣體通過(guò)充氣通道進(jìn)入噴嘴并在噴嘴內(nèi)形成氣泡,并與來(lái)自液體流量計(jì)的液體相混合形成氣液混合物。該氣液混合物離開(kāi)噴嘴后,在不穩(wěn)定波的作用下發(fā)生破裂,使用高速相機(jī)記錄這一過(guò)程,并使用開(kāi)源軟件ImageJ分析處理得到的圖片。實(shí)驗(yàn)通過(guò)控制注入噴嘴的氣體流量從而得到不同液體射流速度下直徑近似相等的氣泡。實(shí)驗(yàn)使用的噴嘴如圖2所示。其中,噴嘴內(nèi)徑d 0=2.96 mm,充氣通道外徑d 1=0.50 mm。

圖1實(shí)驗(yàn)流程Fig.1Sketch of experimental setup

圖2實(shí)驗(yàn)噴嘴Fig.2Sketch of nozzle

實(shí)驗(yàn)采用的表面活性劑為十二烷基苯磺酸鈉(SDBS),是一種陰離子型表面活性劑。在水中加入SDBS會(huì)顯著減小水的表面張力。當(dāng)SDBS濃度較小時(shí),水中的十二烷基苯磺酸鈉呈分子狀態(tài)分散在溶液中。提高SDBS濃度至一定值時(shí),溶液內(nèi)部的SDBS分子會(huì)結(jié)合成較大的基團(tuán),形成膠束。表面活性劑在水中形成膠束所需的最低濃度稱為臨界膠束濃度,即CMC。Zhao等[27]實(shí)驗(yàn)表明SDBS所對(duì)應(yīng)的CMC,使用質(zhì)量分?jǐn)?shù)c表示時(shí),c≈1.00×10-3。

由于表面活性劑溶液的性質(zhì)在CMC附近會(huì)發(fā)生突變,因此實(shí)驗(yàn)分別配制了六種不同濃度的表面活性劑溶液。當(dāng)含表面活性劑射流離開(kāi)射流噴嘴后,射流表面與環(huán)境氣體接觸,形成氣液界面。此時(shí)界面處不存在表面活性劑,溶液表面張力表現(xiàn)為純水的表面張力。溶液內(nèi)部與氣液界面處存在表面活性劑濃度梯度。在這一濃度梯度的作用下,溶解在液體中的表面活性劑分子自發(fā)地從溶液內(nèi)部向射流表面擴(kuò)散并分布在氣液界面上,從而導(dǎo)致溶液表面張力的持續(xù)減小。當(dāng)氣液界面對(duì)表面活性劑分子的吸附達(dá)到動(dòng)態(tài)平衡或氣液界面處完全布滿表面活性劑分子時(shí),氣液界面處的表面張力達(dá)到穩(wěn)定值。因此,獲得SDBS溶液在不同時(shí)刻下的表面張力,即溶液的動(dòng)態(tài)表面張力,是十分必要的。

實(shí)驗(yàn)使用SITA t100型動(dòng)態(tài)表面張力儀(SITA Messtechnik GmbH,Germany)測(cè)量了不同濃度下SDBS溶液的動(dòng)態(tài)表面張力,測(cè)量結(jié)果如圖3所示。該表面張力儀使用氣泡壓力法測(cè)量溶液的動(dòng)態(tài)表面張力。從圖3中可以看出,隨著氣泡壽命的增加,動(dòng)態(tài)表面張力逐漸減小并最終趨于穩(wěn)定。這表明在表面活性劑濃度梯度的作用下,SDBS分子逐漸從溶液內(nèi)部向氣液界面擴(kuò)散,液體動(dòng)態(tài)表面張力逐漸減小。當(dāng)氣泡壽命達(dá)到某一值時(shí),氣液界面對(duì)表面活性劑分子的吸附達(dá)到動(dòng)態(tài)平衡,氣液界面處的表面張力取得穩(wěn)定值,繼續(xù)增大氣泡壽命并不會(huì)導(dǎo)致動(dòng)態(tài)表面張力的進(jìn)一步減小。

圖3不同濃度表面活性劑溶液的動(dòng)態(tài)表面張力Fig.3Dynamic surface tension at various surfactant concentrations

Rosen等[28]根據(jù)表面活性劑溶液動(dòng)態(tài)表面張力隨時(shí)間變化的不同趨勢(shì)將動(dòng)態(tài)表面張力的時(shí)間演變劃分為四個(gè)區(qū)域,分別為誘導(dǎo)區(qū)、快速下降區(qū)、介平衡區(qū)和平衡區(qū)。他們發(fā)現(xiàn)前三個(gè)區(qū)域的動(dòng)態(tài)表面張力隨時(shí)間的演變規(guī)律可以采用式(1)表示:

(1)

式中,σ0為表面活性劑濃度為零時(shí)表面張力,即水的表面張力,mN/m;σt為t時(shí)刻溶液的表面張力,mN/m;σm為介平衡區(qū)溶液的表面張力,一般取動(dòng)態(tài)表面張力變化較小時(shí)的溶液表面張力[29],mN/m;n,t*為修正系數(shù)。對(duì)式(1)兩側(cè)取對(duì)數(shù),有:

(2)

(3)

此時(shí),lg S為lg t的一次函數(shù),通過(guò)擬合實(shí)驗(yàn)數(shù)據(jù)可以得到直線的斜率n與截距-n lg t*。

圖4展示了不同SDBS濃度下lg S隨lg t的變化規(guī)律,圖中實(shí)線對(duì)應(yīng)各濃度下動(dòng)態(tài)表面張力的擬合曲線,此時(shí)圖中時(shí)間t的單位為ms。從圖中可以看出,擬合結(jié)果與實(shí)驗(yàn)數(shù)據(jù)吻合較好,擬合曲線能夠較好地表征動(dòng)態(tài)表面張力的時(shí)間演變趨勢(shì)。擬合得到的系數(shù)見(jiàn)表1。

圖4不同濃度表面活性劑溶液的lgS-lgt曲線圖Fig.4lgS-lgtcurves of various aqueous

(surfactant solutions)

表1動(dòng)態(tài)表面張力擬合參數(shù)Table 1 Fitting parameters of dynamic surface tension

2、實(shí)驗(yàn)結(jié)果與討論

2.1、表面張力變化對(duì)射流破裂特征的影響

本文首先研究了不含氣泡的表面活性劑射流破裂過(guò)程以探究表面張力變化對(duì)射流破裂特征的影響。實(shí)驗(yàn)使用高速相機(jī)記錄了六種不同濃度表面活性劑溶液的射流破裂過(guò)程。實(shí)驗(yàn)獲得的射流破裂過(guò)程如圖5所示,其中Δt是相鄰兩張圖片間的時(shí)間間隔。從圖5中可以看出,雖然表面活性劑濃度并不相同,但射流的破裂過(guò)程在形態(tài)上是相似的。當(dāng)液體射流離開(kāi)噴嘴后,射流表面擾動(dòng)迅速增長(zhǎng),表面波動(dòng)在射流表面呈對(duì)稱分布,這是射流破裂Rayleigh區(qū)的典型特征。隨著表面波動(dòng)的進(jìn)一步發(fā)展,在射流方向上,射流表面形成瘤狀結(jié)構(gòu),射流破裂發(fā)生在兩個(gè)瘤狀結(jié)構(gòu)之間。

圖5含表面活性劑射流的破裂過(guò)程Fig.5Breakup processes of various surfactant-laden jets

為了定量表征表面張力變化對(duì)射流破裂過(guò)程的影響,實(shí)驗(yàn)測(cè)量了含表面活性劑射流的射流破裂長(zhǎng)度。射流破裂長(zhǎng)度定義為噴嘴出口至射流首次發(fā)生破裂處的距離。圖6是含表面活性劑射流破裂長(zhǎng)度隨射流速度u l的變化趨勢(shì),其中L s(mm)表示含表面活性劑射流的射流破裂長(zhǎng)度,L l(mm)為水的射流破裂長(zhǎng)度。從圖6(a)中可以看出,當(dāng)表面活性劑濃度保持不變時(shí),提高液體射流速度會(huì)增大射流破裂長(zhǎng)度。從圖6(b)中可以看出,當(dāng)液體射流速度保持不變時(shí),隨著表面活性劑濃度的增大,射流破裂長(zhǎng)度總體上呈現(xiàn)逐漸增大的趨勢(shì);射流速度越大,射流破裂長(zhǎng)度增大的趨勢(shì)就越明顯。使用相同工況下水的射流破裂長(zhǎng)度無(wú)量綱化含表面活性劑射流的破裂長(zhǎng)度后,可以得到射流破裂長(zhǎng)度比隨射流速度以及表面活性劑濃度的變化趨勢(shì),如圖6(c)、(d)所示。從圖中可以看出,表面活性劑的加入顯著增大了射流破裂長(zhǎng)度比。在表面活性劑濃度不變時(shí),提高液體射流速度對(duì)射流破裂長(zhǎng)度比的影響并不顯著。在射流速度保持不變時(shí),提高表面活性劑濃度會(huì)導(dǎo)致射流破裂長(zhǎng)度比的增大。

圖6含表面活性劑射流破裂長(zhǎng)度Fig.6Breakup lengths of various surfactant-laden jets

在氣液界面處,表面活性劑分子的親水基團(tuán)指向水溶液,親油基團(tuán)指向周圍氣體。在氣液界面處表面活性劑的表面吸附量越大,溶液的表面張力越小。當(dāng)液體射流離開(kāi)噴嘴后,氣液界面開(kāi)始形成。溶解在液體射流內(nèi)部的表面活性劑分子在濃度梯度的作用下向射流表面擴(kuò)散并吸附在射流表面。隨著射流的進(jìn)行,射流氣液界面表面活性劑的表面吸附量逐漸增大,射流表面張力逐漸減小。因此,與表面張力為定值的液體射流不同,處理含表面活性劑射流的破裂過(guò)程時(shí),液體動(dòng)態(tài)表面張力必須被考慮。

液體射流離開(kāi)噴嘴到發(fā)生破裂所需時(shí)間可以定義為射流破裂時(shí)間t b,忽略重力作用時(shí),其可以表示為:

(4)

不同濃度表面活性劑溶液的射流破裂時(shí)間如圖7所示。

圖7含表面活性劑射流的射流破裂時(shí)間Fig.7Breakup time of various surfactant-laden jets

從圖7(a)中可以看出,當(dāng)表面活性劑濃度保持不變時(shí),射流破裂時(shí)間幾乎不隨射流速度的改變而改變。在射流速度保持不變時(shí),隨著表面活性劑濃度的增大,射流破裂時(shí)間呈現(xiàn)逐漸增長(zhǎng)的趨勢(shì)。在實(shí)驗(yàn)參數(shù)范圍內(nèi),射流破裂時(shí)間均小于0.4 s??紤]重力對(duì)液體射流的加速作用,液體射流實(shí)際破裂時(shí)間要小于0.4 s。

對(duì)比圖3中表面活性劑溶液氣泡壽命與圖7中射流破裂時(shí)間可以發(fā)現(xiàn),在時(shí)間小于0.4 s時(shí),動(dòng)態(tài)表面張力尚未發(fā)展至平衡區(qū),動(dòng)態(tài)表面張力隨著時(shí)間的增長(zhǎng)逐漸減小。當(dāng)表面活性劑濃度c<0.1%時(shí),在射流破裂時(shí)間內(nèi),溶液動(dòng)態(tài)表面張力緩慢減小,且與純水的表面張力幾乎相等。當(dāng)表面活性劑濃度c≥0.1%時(shí),在射流破裂時(shí)間內(nèi),動(dòng)態(tài)表面張力呈現(xiàn)先迅速減小而后緩慢降低的趨勢(shì)。

當(dāng)含表面活性劑射流離開(kāi)噴嘴后,其表面會(huì)形成不穩(wěn)定波。射流表面不穩(wěn)定波隨著射流的前進(jìn)不斷發(fā)展并最終導(dǎo)致了液體射流的破裂。從圖5中可以看出,實(shí)驗(yàn)工況下,射流表面波動(dòng)為正對(duì)稱波;從圖6中可以看出,射流破裂長(zhǎng)度隨著液體射流速度的增大而增大。因此,可以認(rèn)為在實(shí)驗(yàn)工況下,含表面活性劑射流處于Rayleigh模式。在此射流破裂模式下,射流表面最不穩(wěn)定波的增長(zhǎng)率可以表示為[18]:

(5)

射流表面波振幅的時(shí)間演變規(guī)律可以表示為[9]:

(6)

式中,η和η0分別為射流表面波動(dòng)振幅和初始振幅,m。假設(shè)表面波振幅發(fā)展至射流半徑時(shí),射流發(fā)生破裂。因此,射流破裂長(zhǎng)度可以表示為[9]:

(7)

在含表面活性劑射流中,表面活性劑濃度較低,對(duì)除液體表面張力以外的其他液體性質(zhì)影響較小。因此,當(dāng)射流剛離開(kāi)射流噴嘴,新的氣液界面開(kāi)始形成時(shí),射流表面初始擾動(dòng)振幅與水的射流初始振幅相同,即≈10.14[24]。

當(dāng)液體射流其他性質(zhì)保持不變時(shí),減小液體表面張力會(huì)導(dǎo)致射流表面最不穩(wěn)定波增長(zhǎng)率的減小。當(dāng)液體射流內(nèi)部含有表面活性劑時(shí),液體射流表面張力隨著時(shí)間的增長(zhǎng)逐漸減小,射流表面最不穩(wěn)定波的增長(zhǎng)率逐漸減小[30-32],表面波振幅發(fā)展至射流半徑所耗費(fèi)時(shí)間大于同工況下水的射流所耗費(fèi)時(shí)間,射流破裂長(zhǎng)度大于同工況下水的射流破裂長(zhǎng)度。當(dāng)表面活性劑濃度c<0.1%時(shí),在射流破裂時(shí)間內(nèi),溶液動(dòng)態(tài)表面張力與水的表面張力幾乎相等,因此含表面活性劑射流的破裂長(zhǎng)度與水的射流破裂長(zhǎng)度幾乎相等。當(dāng)表面活性劑濃度c≥0.1%時(shí),在射流破裂時(shí)間內(nèi),動(dòng)態(tài)表面張力顯著小于水的表面張力。在相同時(shí)刻,動(dòng)態(tài)表面張力隨著表面活性劑濃度的增大而減小。因此,含表面活性劑射流破裂長(zhǎng)度大于水的射流破裂長(zhǎng)度,且隨著表面活性劑濃度的增大,射流破裂長(zhǎng)度逐漸增大。

在實(shí)驗(yàn)工況范圍內(nèi),含表面活性劑射流的動(dòng)態(tài)表面張力隨著時(shí)間的推移不斷減小。為了定量表征表面活性劑對(duì)射流破裂的影響,使用t=0與t=t b時(shí)刻動(dòng)態(tài)表面張力的平均值作為射流平均表面張力:

(8)

式中,、分別為t=0與t=t b時(shí)刻的動(dòng)態(tài)表面張力,mN/m。在t=0時(shí)刻,射流開(kāi)始離開(kāi)噴嘴,氣液界面開(kāi)始形成,氣液界面尚未吸附表面活性劑分子。因此,采用表面活性劑濃度c=0時(shí)的液體表面張力作為此時(shí)的動(dòng)態(tài)表面張力。

因此,射流表面最不穩(wěn)定波的平均增長(zhǎng)率為:

(9)

從式(2)~式(9)可以看出,隨著表面活性劑濃度的增大以及時(shí)間的推移,動(dòng)態(tài)表面張力逐漸減小,射流表面平均表面張力逐漸減小,因而導(dǎo)致射流表面最不穩(wěn)定波平均增長(zhǎng)率降低,射流破裂長(zhǎng)度逐漸增大。這與實(shí)驗(yàn)觀察到的射流破裂長(zhǎng)度變化趨勢(shì)一致。

射流表面的表面活性劑分子不僅會(huì)減小射流表面最不穩(wěn)定波的增長(zhǎng)率,其在射流表面的不均勻分布還會(huì)引發(fā)Marangoni效應(yīng),從而影響射流破裂過(guò)程[22,32-33]。

如圖5所示,隨著射流的發(fā)展,處于兩個(gè)瘤狀結(jié)構(gòu)中間的液體被不斷拉伸,其表面積逐漸增大,從而導(dǎo)致表面活性劑表面吸附量的減小。同時(shí),由于射流內(nèi)部與射流表面存在表面活性劑濃度梯度,在這一濃度梯度的作用下,表面活性劑分子從溶液內(nèi)部向射流表面擴(kuò)散,從而導(dǎo)致表面活性劑表面吸附量增大。在射流破裂時(shí)間內(nèi),不同濃度表面活性劑溶液的動(dòng)態(tài)表面張力均未達(dá)到平衡區(qū),表面活性劑尚在從射流內(nèi)部向射流表面擴(kuò)散。

Peclet數(shù)Pe常被用于表征對(duì)流速率與擴(kuò)散速率的相對(duì)大小,本文使用Peclet數(shù)表征液體拉伸導(dǎo)致的表面活性劑表面吸附量減小趨勢(shì)與擴(kuò)散導(dǎo)致的表面活性劑表面吸附量增大趨勢(shì)的相對(duì)強(qiáng)弱:

(10)

式中,L c為特征長(zhǎng)度,m;u l為液體射流速度,m/s,u l~1;D s為質(zhì)量擴(kuò)散系數(shù),m 2/s,溶液的擴(kuò)散系數(shù)D s~10-9。本文使用最不穩(wěn)定波的半波長(zhǎng)λ/2作為特征長(zhǎng)度L c。Rayleigh[8]指出,促使液體射流發(fā)生破裂的最不穩(wěn)定波的波長(zhǎng)滿足:

(11)

聯(lián)合式(10)和式(11)可以得到Pe>>1。因此,液體拉伸所導(dǎo)致的表面活性劑表面吸附量減小作用占據(jù)主導(dǎo)地位。隨著液柱的拉伸,表面活性劑的表面吸附量逐漸減小。

在毛細(xì)壓力的作用下,射流液柱直徑逐漸減小并最終導(dǎo)致了射流的破裂。在接近射流破裂處,毛細(xì)壓力驅(qū)動(dòng)流體向破裂點(diǎn)兩側(cè)運(yùn)動(dòng),破裂點(diǎn)附近表面活性劑的表面吸附量接近零,表面張力接近水的表面張力[22]。在瘤狀結(jié)構(gòu)處,表面活性劑的表面吸附量最大,表面張力最小。因此,射流表面存在自表面波波谷指向波峰的表面張力梯度,從而產(chǎn)生Marangoni效應(yīng)。在Marangoni應(yīng)力的作用下,液體自低表面張力的波峰向高表面張力的波谷運(yùn)動(dòng),從而導(dǎo)致了液體射流破裂的延遲和射流破裂長(zhǎng)度的增大。

Marangoni應(yīng)力τM可以表示為:

(12)

在射流破裂過(guò)程中促使射流直徑收縮的毛細(xì)壓力pσ可以表示為:

(13)

因此,定義Marangoni應(yīng)力與毛細(xì)壓力的比值Z以表征Marangoni效應(yīng)對(duì)射流破裂的阻礙作用:

(14)

由于射流破裂發(fā)生在Rayleigh區(qū),聯(lián)合式(11)并簡(jiǎn)化可以得到:

(15)

因此,綜合考慮表面張力變化對(duì)射流表面最不穩(wěn)定波增長(zhǎng)率的影響以及Marangoni效應(yīng),液體射流破裂長(zhǎng)度可以表示為:

(16)

(17)

式中,a和b為關(guān)于Marangoni效應(yīng)的修正系數(shù)。

通過(guò)擬合實(shí)驗(yàn)數(shù)據(jù)得到a=1.91,b=1.58。不同濃度表面活性劑射流的射流破裂長(zhǎng)度擬合曲線如圖8所示。從圖8中可以看出,擬合結(jié)果和實(shí)驗(yàn)數(shù)據(jù)吻合較好。擬合曲線能夠較好地描述不同表面張力下射流破裂長(zhǎng)度隨射流速度的變化趨勢(shì)。圖9是擬合公式計(jì)算值與實(shí)驗(yàn)測(cè)量值之間的比較。從圖中可以看出,理論預(yù)測(cè)值與實(shí)驗(yàn)測(cè)量值吻合較好。

圖8含表面活性劑射流的破裂長(zhǎng)度(圖中實(shí)線為對(duì)應(yīng)的不同濃度表面活性劑射流的擬合曲線)Fig.8Breakup lengths of various surfactant-laden jets(solid lines correspond to fitting curves)

圖9射流破裂長(zhǎng)度實(shí)驗(yàn)測(cè)量值與擬合計(jì)算值的比較Fig.9Comparison between experimental results and calculated results

2.2、表面張力變化對(duì)含氣泡射流破裂特征的影響

圖10展示了三種不同濃度表面活性劑射流在內(nèi)部氣泡作用下發(fā)生破裂的過(guò)程。從圖10中可以看出,不同濃度表面活性劑溶液射流受內(nèi)部氣泡影響發(fā)生破裂的過(guò)程是相似的。隨著射流的發(fā)展,表面波振幅逐漸增大并最終導(dǎo)致了液體射流的破裂。射流破裂發(fā)生在氣泡附近。雖然不同濃度表面活性劑射流的破裂過(guò)程在形態(tài)上是相似的,但射流破裂長(zhǎng)度并不相同。內(nèi)部氣泡對(duì)射流破裂長(zhǎng)度存在影響。實(shí)驗(yàn)測(cè)量了含氣泡的表面活性劑射流破裂長(zhǎng)度L sb以定量表征射流內(nèi)部氣泡對(duì)含表面活性劑射流破裂過(guò)程的影響。射流內(nèi)氣泡尺寸的測(cè)量見(jiàn)文獻(xiàn)[24]。

圖10含氣泡表面活性劑溶液射流破裂過(guò)程Fig.10Breakup process of surfactant-laden jet with inner bubbles

圖11是不同濃度表面活性劑射流破裂長(zhǎng)度L sb隨氣泡無(wú)量綱直徑D=d b/d 0(d b為氣泡直徑)的變化。從圖中可以看出,不同表面活性劑濃度下,射流破裂長(zhǎng)度隨氣泡無(wú)量綱直徑D的變化趨勢(shì)是相似的。射流破裂長(zhǎng)度隨著氣泡直徑的增大逐漸減??;當(dāng)氣泡直徑相同時(shí),提高液體射流速度會(huì)導(dǎo)致更大的射流破裂長(zhǎng)度。

圖11內(nèi)部氣泡直徑對(duì)含表面活性劑射流破裂長(zhǎng)度的影響Fig.11Surfactant-laden jet breakup length versusbubble diameter

實(shí)驗(yàn)還研究了含氣泡射流破裂長(zhǎng)度比隨氣泡直徑的變化趨勢(shì),如圖12所示。射流破裂長(zhǎng)度比定義為含氣泡射流破裂長(zhǎng)度對(duì)同工況下不含氣泡射流破裂長(zhǎng)度的比值。從圖12可以看出,隨著氣泡直徑的增大,破裂長(zhǎng)度比逐漸減小,這一趨勢(shì)在不同表面活性劑濃度和液體射流速度下均可被觀察到。這是因?yàn)闅馀莸囊雽?dǎo)致了液體射流速度的波動(dòng)從而縮短了射流破裂長(zhǎng)度。氣泡直徑越大,引起的速度波動(dòng)越強(qiáng)烈,射流破裂長(zhǎng)度比越小[24]。

圖12含表面活性劑射流的破裂長(zhǎng)度比隨氣泡直徑的變化Fig.12Surfactant-laden jet breakup length ratiosversus bubble diameter

當(dāng)氣泡被注入射流后,氣泡不僅會(huì)引起射流速度擾動(dòng),而且會(huì)產(chǎn)生新的氣液界面并吸附表面活性劑分子。氣泡表面對(duì)表面活性劑分子的吸附會(huì)降低表面活性劑濃度,減小射流內(nèi)部與射流表面間的表面活性劑濃度梯度。因此,相對(duì)于不含氣泡射流,相同時(shí)刻含氣泡射流表面的表面活性劑表面吸附量較小,動(dòng)態(tài)表面張力較大,射流破裂長(zhǎng)度較小。動(dòng)態(tài)表面張力的增大還會(huì)減小波谷和波峰間的表面張力梯度,從而削弱Marangoni效應(yīng),導(dǎo)致更小的射流破裂長(zhǎng)度。因此,射流內(nèi)部氣泡的存在會(huì)導(dǎo)致射流破裂長(zhǎng)度的減小。

隨著氣泡直徑的增大,氣泡表面能夠吸附的表面活性劑分子增多,射流內(nèi)部與射流表面的表面活性劑濃度梯度減小,相同時(shí)刻射流表面的表面活性劑表面吸附量減小,動(dòng)態(tài)表面張力增大,射流破裂長(zhǎng)度減小;同時(shí),Marangoni效應(yīng)也被削弱,射流破裂長(zhǎng)度減小。因此,增大氣泡直徑會(huì)導(dǎo)致更小的射流破裂長(zhǎng)度比。

如圖12所示,當(dāng)氣泡直徑較大時(shí),破裂長(zhǎng)度比基本不隨射流速度的改變而改變。這是由于此時(shí)氣泡引起的射流速度擾動(dòng)對(duì)射流破裂的影響占主導(dǎo)地位。在氣泡直徑較小時(shí),當(dāng)表面活性劑濃度c<1.00×10-2,不同射流速度下射流破裂長(zhǎng)度比基本相同。當(dāng)表面活性劑濃度c=1.00×10-2,破裂長(zhǎng)度比隨著射流速度的增大而逐漸減小。這可能與被削弱的Marangoni效應(yīng)有關(guān)。由式(10)可知,液體射流速度越大,液體對(duì)流造成的表面活性劑表面吸附量減少的趨勢(shì)越強(qiáng),波峰波谷間的表面張力差值越小,Marangoni效應(yīng)越弱,射流破裂長(zhǎng)度增加的趨勢(shì)越弱。因此,在氣泡較小時(shí),隨著射流速度的增大,含氣泡射流的破裂長(zhǎng)度比逐漸減小。在表面活性劑濃度較小時(shí),氣泡的存在進(jìn)一步減小了射流內(nèi)部表面活性劑的濃度,從而使得不同射流速度下Marangoni效應(yīng)的差別不明顯,不同射流速度下射流破裂長(zhǎng)度比基本相同。

實(shí)驗(yàn)還得到了不同射流速度下,射流破裂長(zhǎng)度比隨氣泡直徑的變化曲線,如圖13所示。從圖中可以看出,當(dāng)射流速度相同時(shí),不同濃度表面活性劑溶液的射流破裂長(zhǎng)度比是接近的,且?guī)缀醵夹∮谕葪l件下水的射流破裂長(zhǎng)度比。這表明射流內(nèi)部氣泡對(duì)含表面活性劑射流破裂過(guò)程的影響更為顯著。

圖13不同射流速度下的射流破裂長(zhǎng)度比(圖中實(shí)線為水的含氣泡射流破裂曲線)Fig.13Surfactant-laden jet breakup length ratios at various jet velocities(solid line corresponds to curve of water jet breakup)


3、結(jié)論

本文研究了表面張力變化下,處于Rayleigh模式的液體射流破裂過(guò)程,揭示了射流內(nèi)部氣泡對(duì)含表面活性劑射流破裂過(guò)程的影響,得到的主要結(jié)論如下。

(1)在液體射流速度保持不變時(shí),隨著表面活性劑濃度的提高,液體表面張力逐漸減小,液體射流破裂長(zhǎng)度逐漸增大。表面活性劑對(duì)射流破裂的影響主要通過(guò)動(dòng)態(tài)表面張力及表面活性劑在射流表面的非均勻分布實(shí)現(xiàn)。表面活性劑的加入使得溶液動(dòng)態(tài)表面張力隨著時(shí)間的推移逐漸減小,射流表面不穩(wěn)定波增長(zhǎng)率逐漸降低,表面波振幅發(fā)展到射流半徑所需時(shí)間延長(zhǎng),射流破裂長(zhǎng)度增大。表面活性劑在射流表面的不均勻分布會(huì)引發(fā)Marangoni效應(yīng)。液體在Marangoni應(yīng)力的作用下向射流拉伸區(qū)運(yùn)動(dòng),從而延緩了射流破裂的發(fā)生,增大了射流破裂長(zhǎng)度。得出了液體射流破裂長(zhǎng)度的表達(dá)式[式(16)],擬合結(jié)果與實(shí)驗(yàn)測(cè)量結(jié)果吻合較好。

(2)實(shí)驗(yàn)發(fā)現(xiàn)射流內(nèi)部氣泡會(huì)顯著縮短含表面活性劑射流的破裂長(zhǎng)度。射流內(nèi)部氣泡的引入會(huì)導(dǎo)致射流速度的波動(dòng),從而縮短射流破裂長(zhǎng)度。射流內(nèi)部氣泡對(duì)表面活性劑分子的吸附會(huì)降低表面活性劑濃度,減小表面活性劑的表面吸附量,提高射流表面的動(dòng)態(tài)表面張力,從而縮短了射流破裂長(zhǎng)度。此外,射流內(nèi)部氣泡導(dǎo)致的較大的動(dòng)態(tài)表面張力還會(huì)減小Marangoni應(yīng)力,削弱Marangoni效應(yīng),從而縮短了射流破裂長(zhǎng)度。射流內(nèi)部氣泡越大,射流破裂長(zhǎng)度越小。與水的含氣泡射流相比,射流內(nèi)部氣泡對(duì)含表面活性劑射流破裂過(guò)程的影響更為顯著。

国产无码.con| 国产精品国产三级国产普通话一| 自拍偷拍第十页| 久久精品无码一区二区三区| 91福利导| 亚洲无码二区| 无码精品人妻| 97国精产品无人区一码二码 | 久久国产综合| 国产精品久久国产精品99无码| 久久亚洲免费视频| 91aaa| 麻豆性爱视频| 国产精品国精产品一二三| 中文字幕日韩精品无码内射| 久久性爱视频| 国产高清无码免费| 亚洲人妻视频| 国产手机视频在线| 黄色片毛片| 久久无码电影| 免费一级大黄片| 中文字幕免费在线| 嘿嘿嘿视频免费网站| 亚洲免费在线| 99国产精品99久久久久久 | 黄色免费无码视频网站| 日本色色网| 91精品国产人妻女教师| 国产精品无码在线观看| 乱伦老女人一区二区| 亚洲AV丰满熟妇在线播放| 超碰97在线免费观看| 狠狠操影院| 天天日天天操天天干| 午夜福利国产| 无码在线观看一区| 欧美碰碰| 国模一区二区| 国产精品1区| 国产免费一级| 国产一级性爱| 亚洲国产二区| 中文字幕人妻无码| 午夜欧美精品久久久久久久| 国产视频不卡| 成人做爰A片一区二区| 亚洲第一网站| 亚洲狠狠干| 黄片免费视频| 97视频在线| 欧美日一区二区三区| 色爱综合网| 曰韩无码视频| 成人精品一区二区| 免费观看黄色网址| 国产农村高清无套内谢视频| 国产毛片在线| 中文字幕在线视频观看| 欧美污视频| 天天综合久久综合| 老女人chinese肥臀老女人| 91麻豆精品国产91久久久久久| 一本一道久久综合狠狠躁牛牛影视 | 国产00粉嫩馒头一线天91| 亚洲啪啪视频| 九色国产| 精品人妻少妇一区二区三区在线| 97超蹦在线人艹人| 国产成人精品久久| 国内精品久久久久久久影视4| 亚洲男人天堂| 亚洲三级片在线播放| www操笔网站| 国产精品天堂一区二区在线观看| 色欲精品人妻AV一区| 4438xx亚洲五月最大丁香| 亚洲国产熟妇伦| 不卡欧美| 中文字幕日产A片在线看| 国产激情一区二区三区| av无码aV天天aV天天爽| 欧美高清一区二区| 一级特黄大片69| 热久久这里只有精品| 乱伦精品| 天天日天天干天天操| 影音先锋男人资源网| 91性视频| 久久久999| 免费看的av| 成人免费观看网站| 色综合88| 日日操日日爽| 色妺妺视频网| 国产一级A片夜天码免费看| 91麻豆精品久久久久蜜臀| 欧美视频一区二区| 日韩无码专区| 波多野结衣中文字幕久久| 我跟闺蜜公交车被弄到高潮| 秋霞影院午夜丰满少妇在线视频| 草一次黄色av| 色妞综合网| 丁香五月社区| 91久久电影| A毛片网站| 成人精品视频在线| 中文字幕一区二区三区日韩精品| 涩涩视频网站| 91精品国产高清91久久久久久| 国产精品一区二区黑人巨大| 色无码视频| 一本一道久久a久久精品综合蜜臀| 熟女VS乱伦| 超碰在线伊人| 一级a一级a爰片免费免水l软件| 91这里只有精品| 国产91在线拍揄自揄拍无码九色| 国产一区二区三区免费视频| 一级黄色电影在线观看| 高清无码在线观看av| 人人操这里只有精品| 最新国产无码| 久久精品91| 亚洲av网站| 婷婷五月天视频| 97资源超碰| 午夜久久久| 污污网站在线观看| 日韩一级黄色电影| 欧洲-级毛片内射| 偷拍自拍网| 精品人妻久久| 香蕉色a片| 一级黄片免费视频| 日韩性爱无码| 久久激情综合| 又做又爱视频免费| 亚洲无码性爱| 国产1区2区3区中文字幕| 西西图吧| 国产精品免费观看| 逼操逼操逼操逼操| 午夜福利视频| 成人区精品一区二区婷婷| 无码深夜AAA片在线观看| aaa无码| 日本久草| 日韩欧美一级大片| 亚洲操逼网| 久久久无码精品人妻二区| 亚洲精品国产suv一区| 91网站在线播放| 日本天堂在线| 91性爱网站| MM1313又粗又大受不了| 国产一二精品| 欧美色图一区二区三区| 色综合久久久| 日日夜夜网站| 亚洲图片欧美日韩| 亚洲AV午夜精品无码专区在线 | 久久天天躁狠狠躁夜夜躁2014| 五月丁香五月婷婷| 国产一区二区三区在线视频| 日木精品人妻| 后入内射欧美99二区视频| 中文字幕人妻一区二区…| 久久久18禁一区二区三区精品| 啪啪导航| 无码人妻aⅴ一区二区三区69堂| 91丝袜视频| av在线一区二区| 天天操夜夜操免费视频| 国产精品极品白嫩在线| 成人午夜福利| 九九超碰| 亚洲精品日韩激情在线电影| 中文字幕黄色| 在线观看一区| 巨大巨粗巨长 黑人长吊| 亚洲一区免费| 国产日本精品| 亚洲精品自拍| 欧美亚洲天堂| √8天堂资源地址中文在线| 久久午夜视频| 日韩精品影院| 无码少妇一二三区免费| 国产精品亚洲无码| 极品视频在线| 激淫少妇被插视频在线观看| 国产精品日韩无码| 黄色成人在线| 国产一级A片在线观看免费视频| 国产精品久久久久久人妻黑料| 色窝窝无码一区二区三区成人网站 | 日韩精品免费一区二区夜夜嗨| 一级a一级a爰片免免免下载| 天天日天天干天天操| 蜜桃av在线播放| 秋霞无码视频| 亚洲AV无码国产精品久久不卡嫖娼| JlZZJlZZ亚洲日本少妇| 狼友91精品一区二区三区| 免费看黄网址| 91人妻无码精品一区二区毛片| 国产操逼视频免费观看| 四虎色播| 又粗又爽又猛高潮的在线视频| 极品91尤物被啪到呻吟喷水| 免费观看黄片| 亚洲欧美一级特黄大片| 久久久久久三级片| 一级毛片网址| 日韩黄片一区| 欧美视频三区| 在线中文字幕| 欧美怡春院| 国产美女裸体永久免费观看网站 | 亚洲AV成人www新版精品久久| 国产一区二区电影| 欧美性爱日韩高清| 好看的操逼视频| 国产免费无码av| 中文字幕国产| 国产激情在线观看| 奇米狠狠| 中文字幕无码av| 国产在线视频网站| 欧美一级日韩一级| 国产毛片在线| 欧美午夜激情| 欧美日韩日逼| 蜜桃狠狠干网| 亚洲中文国产精品| 欧美1区2区| 一区二区www| 欧美午夜在线视频| 中文字幕有码视频| 亚洲午夜无码AV毛片久久| 成人在线毛片| 精品福利| 水果派解说一区二区三区在线观看| 亚洲精品无码视频| 青娱乐极品视频| 亚洲无码国产精品| 欧美日韩黄色| 一级黄色网址| 亚洲无毛| 国产一码二码三码四码无码| 欧美一级黄色大片| 久久精品精品无码一区三区| 午夜无码片在线观看影院| 日韩精品久久久| 91无码高清视频| 国产一级a毛一级a看免费人娇| 午夜视频一区| 欧美激情一区| 伊人热久久| 国产a级免费| 永久黄网站色视频免费直播| 日韩美女在线| 欧美精品一区二区在线| 91丨九色丨熟女高潮| 天天干,夜夜操| av中文字幕一区| 亚欧艹逼| 高清无码黄| 亚洲少妇一区二区| 深夜成人视频在线| 天天干狠狠干| 粉嫩av久久一区二区三区小说| 久久精品不卡| 一级毛片久久久久久久18| 中文字幕乱偷无码av一区二区| 国产性爱一级片| 国产一级自拍| 国产日韩欧美在线| 在线免费看黄| 波多野结衣精品视频| 久久这里有精品| 免费国产网站| 精品日韩在线| 国产精品色哟哟| 在线观看无码视频| 日韩精品在线播放| 日韩免费一区二区三区| 国产成人精品亚洲日本在线观看 | 国产白浆视频| 疯狂的交换1—6真实交换3和2| 午夜精品久久久久久久| 美女污网站| 欧美成人性爱视频在线观看| 欧美多毛熟妇| 黄色无码大片| 亚洲国产精品自拍| 中文无码熟妇人妻AV在线| 国产a区| 日本护士高潮大叫| 成人大香蕉| 岛国高清无码| 久久成人麻豆午夜电影| 亚洲久草| 中韩XXX抄逼| 成人无码视频在线观看| 99大香蕉| 精灵梦叶罗丽第八季| 精品自拍视频| 狠狠操观看视频| 国产三级片在线看| 一级黄色片在线观察| 日本精品视频| 91人妻在线| 亚洲自拍一区| 国产 丝袜 另类 精品 综合| 国产色拍| 影音先锋国产资源| 欧美老熟妇一区二区三区| 亚洲黄色电影在线观看| 极品美女一区二区三区| 欧美成人精品一区二区男人看 | 亚洲国产片| 超碰亚洲| 蜜桃久久av无码牛牛影视| 在线观看91| 一本久道久久综合狠狠爱| 激情五月天在线| 精品一级黄片| 日韩高清无码一区| 人妻少妇精品视频免费看蜜桃| 亚洲女同视频| 无码白丝强行免费| 国产成人无码专区| 亚洲成人无码在线| 成人乱人伦一区二区三区| 日韩欧美三级在线| 日韩91| 无码在线免费| 日本视频一区二区三区| 乱伦熟妇| 伊人影院在线观看| 淫荡网站| 精品人妻一区| 高清无码一区二区三区| 日韩精品免费视频| 国产精品无码在线观看| 国内自拍偷拍视频| 国产成人在线视频观看| 亚洲高清在线无码| 人人摸人人看| 午夜天堂一区二区三区| 亚洲高清一区二区三区| 中文字幕不卡在线观看| 色哟呦AV永久免费| 天堂AV国产一区二区熟女人妻| 欧美日韩成人影院| 国产成人三区| 琪琪午夜福利| 日韩无码一区二区三区| 自拍偷拍欧美亚洲| 亚洲AV午夜精品一区二区三区| 久久99精品久久久久婷婷| AV无码波多野结衣| 污视频网站在线观看| 免费观看黄色网| 日韩乱伦一区| 久久99精品久久久水蜜桃| 国产在线观看一区二区| 久久久艹| 九九九国产视频| 午夜视频免费在线观看| 欧美人妻曰韩精品| 国产性爱一级| 日韩精品久久久久久免费| 亚洲熟妇无码AV| 日韩免费在线观看视频| 黄色免费AV| 天天草av| 免费A片三p视频| 亚洲少妇无套内射激情视频| 亚洲无码免费在线视频| 久久综合久| 久久亚洲电影| 好吊视频一区二区三区| 一区二区三区无码视频| 久久久久久影院| 中文字幕亚洲天堂| 国产成人在线视频播放| 青青在线| 嗯啊不要在线观看| 国产精品免费在线| 亚洲激情AV| 亚洲AA| 欧美性爱第1页| 亚洲一级AV无码毛片| 欧美午夜影院| 天堂а在线中文在线新版| 欧美天堂在线观看| h片在线免费观看| 国产精品天天狠天天看| 免费无码国产www| 日本三级网站| 一本久久精品久久综合桃色| 免费高潮视频| 久久91视频| 人人愛人人操| 免费无码视频| 91久久久精品| 成人一区视频| 国产精品乱伦视频| 乱熟女高潮一区二区在线| 国产老熟女伦老熟妇露脸| 色无码在线| а√天堂资源国产精品| 久久九九99| 无码不卡在线| 精品爆乳一区二区三区无码AV| 高清性色生活片| 国产亚洲AV永久无码国产天堂| 色欲久久久| 一区二区激情| 日本a视频| 欧美精品久久久久| 99毛片| 国产做受69高潮精品王| 少妇AV一区二区三区无码按摩| 超碰毛片| 操逼无码免费视频| 国产视频二区| 捷克视频一区二区三区无码| 91性爱视频| 亚洲欧美中文字幕| 不卡无码AV| 蜜桃伊人| 亚洲精品自拍| 无码一级毛片| 污视频在线| 伊人婷婷| 免费操逼| 成年人在线视频| 妖精视频黄色| 无码精品一区二区免费JIZZ| 中文字幕日韩人妻在线视频| 国产黄片免费在线观看| 久久精品影视| 91久久精品无码一区二区三区| 亚洲精品久久无码77777| av无码aV天天aV天天爽| 日韩在线观看AV| 91蜜桃在线免费观看| www欧美| 成人性生交大片免费看5| 欧美性爰一二三区| 国产又粗又猛又黄| 欧美日韩在线观看视频| 99亚洲精品| 日韩黄视频| 国产永久精品| 亚洲精选在线| 国产精品久久成人网站水多多| www毛片| 日日干天天干| 天天操天天透| 日本三级黄色麻豆| 日本巜侵犯人妻人伦| 无码喷水| 国产无码精品电影| 碰碰人人| 三级片中文字幕在线观看| 久久精品1| 无码人妻一区二区三区在线 | 91免费看视频| 视频一区二区在线观看| 久久99视频精品| 久久精品嫩草影院| 无码在线免费看| 一男一女一级一片| 一α一α在线看| 欧美一区二区无码三区有限公司 | 熟女肥臀白浆大屁股一区二区| 日本中文字幕在线播放| 天天插天天日| 欧美日韩一区二区三区四区五区| 少妇人妻一区二区三区| 亚洲国产精品一区二区三区| 欧美人与性动交α欧美精品 | 8050午夜| 国产中出| 久久成人影视| 成人精品一区二区三区| 欧美在线视频免费观看| 亚洲精品www| 亚洲成人AV在线| 亚洲国产精品无码观看久久| 在线无码播放| 久久91视频| 久久无码影视| 国产精品爽爽久久久久久豆腐| 亚洲精品一区23p| 久久欧美国产伦子伦精品按摩| 黄色网在线播放| 国产性色视频| 一级a一级a爰片免费免水l软件| 精品伊人| 作爱网站| 美日韩一区二区| 久久午夜影院| 无码人妻久久一区二区三区免费人妻| 女人高潮抽搐喷液30分钟视频| 欧美一级免费| 成人午夜福利| 亚洲无码aaa| 亚洲国产精品自拍| 亚洲无码免费网站| 欧美精品一卡二卡| 日本久久久久久| 亚洲国产AV自拍| 在线观看a视频| 视频一区二区在线观看| 国产在线成人| 99精品无码| 黄色高清无码性爱| 婷婷超碰| 中文字幕第一页在线| 99久久精品免费看国产免费软件| 97久久超碰| 精品国产一区二区三区性色AV| 久久精品久久精品| av亚洲欧洲日产国码无码苍井空| 国产精品久久久久久久AV超碰| 欧美天天干| 北条麻妃视频在线观看| 国产成人午夜视频| 波多野结衣网址| 欧美精品少妇| 九九香蕉视频| 久久久精品一区二区三区| 久久亚洲一区二区| 成人网站在线看| 一级黄片在线播放| 亚洲无码一区在线观看| 日本三级不卡| 国产高清视频在线| 无码视频在线| 成人欧美一区二区三区黑人免费 | 日本人妻在线播放| 久久国产一区二区深田咏美| 国产在线精品免费aaa片| 国产人妻精品一区二区三水牛| 少妇人妻偷人精品无码视频新浪| 国产黄色自拍视频| 国产精品婷婷| 嘿嘿射在线| 久久国产精品无码一级毛片| 91人妻无码一区二区久久| 超碰99在线| 国产乱码精品1区2区3区| 色婷婷影视| 亚洲国产成人精品久久久国产成人一区 | 奇米网| 久久久久久99| 久热中文字幕| 国产黄色在线播放| 拳交网| 国产69Av| 国产精品一区二区三| 老司机午夜影院| 婷婷色一二三区波多野结衣| 国产av无码片毛片一级流奶水| japanese老熟妇乱子伦视频| 黄色片免费网址| 污污污免费网站| 精品69| 国产精品av久久久| 三人成全免费观看电视剧高清| 91无码人妻精品一区二区三区四| 天天草天天爽| 人妻日韩中文字幕| 免费无码国产免费| 成人写真福利网| 中文字幕人妻系列| 亚洲蜜桃妇女| 综合久久久久| 国产第二页| 性色一区| 国产日韩亚洲欧美| 三级三级久久三级久久18| 青青草91| 狼友导航| 丁香激情五月| 国产成人网站在线观看| 亚欧免费视频| 久久午夜影院| 99久久精品国产一区二区三区| 一级特黄60分钟免费| 亚洲熟女乱综合一区二区三区| 欧美日韩久久久久| 伊人成人电影| 一起草视频免费观看无码| 99人妻碰碰碰久久久久禁片| 亚洲乱妇| 91精品久久久久久粉嫩| 粉嫩aⅴ一区二区三区四区五区 | 在线观看一区| 亚洲激情视频在线| 无套内谢波多野结衣| 特一级黄色片| 91中文在线| 亚洲日逼视频| 欧美日韩毛| 久操网站| 色鬼网站| 狠狠干av| 久久久久中文字幕| 免费日逼视频| 国产精品99| 人妻毛片| 日屁视频| 香蕉视频国产| 91视频黄| 超碰免费在线| 99人妻碰碰碰久久久久禁片| _中国一级特黄大片在线看| 久久日本无码中文字幕三级伦| 人人操人人插人人性| 人人操人人早| 中文字幕乱码亚洲精品一区| 99久久99久久久精品棕色圆| 91精品国产高清一区二区三区蜜臀| 久久精品国产亚洲AV高清色欲| 久久99精品久久久久婷婷| 欧洲av在线| 91AV视频在线观看| 精品久久99| 无码人妻一区| 婷婷五月综合在线| 在线观看a v| 欧美精品高清| 亚洲无码一区二区在线| 国产小视频在线| 最新国产在线| 无码成人黄网站在线观看| 国洲 一区二区| 少妇高潮毛片免费看欧美| 亚洲性爱无码| 三级国产| 色无码在线| 图片区偷拍区小说区| 欧美精品一区二区三区久久久竹菊 | 亚洲熟人妇一区二区三区| 九九香蕉视频| 一区二区三区四区无码| 亚洲在线视频| 99国产一区| 亚洲高清成人| 影音先锋男人资源网| 不卡av在线| www夜夜操| 色色欧美| 日本一区二区三区视频在线| 黄污视频| 亚洲AV永久无码精品国产精 | 久久免费小视频| 国产又粗又长又硬| 国产一区二区三区免费观看| 免费精品视频一区二区三区| 国产AV一区二区三区| 69久久| 欧美人人操人人舔| 久久性视频| 欧美色逼| 欧美乱码精品一区二区三| 一区二区高清| 青青免费在线视频| 波多野结衣一区二区| 国产精品久久久久久一级毛片探花| 狠狠干成人| 国产中文字幕视频| 白嫩少妇激情无码| 亚洲精品中文字幕| 一级黄色小视频| 亚洲午夜视频| 日韩中文在线| 久久无码一区二区三区| 欧美日韩一区二区三区不卡视频| 乱伦精品| 日韩少妇人妻| 无码H乳在线看| 三级黄视频| 国产电影一区二区三曲| 一本色道久久综合亚洲精品小说 | 三级片免费网址| 一级a一级a爰片免费啪啪女女| 亚洲免费精品| 免费A级视频| 欧洲熟妇的性久久久久久| 人人草人人摸| 久久久久久亚洲综合影院红桃| 亚洲无码aaa| 天天躁日日躁AAAAXXXX| 国产91色在线观看| 大香蕉99| 一区二区三区四区免费视频| 国产乱伦黄片| 狼友91精品一区二区三区| 香蕉三级片| 97成人站| 天天夜夜爽| 无码人妻一区二区三区免费九色| 在线观看国产黄片| 国产精品偷伦免费观看视频| 无码喷水| 蜜桃久久久| 久久精品视频6| 国产精品久久久久久亚洲调教| 国产suv精品一区二区| 中文字幕手机在线视频| 欧美在线一二三区| 国产精品三级| 丰满人妻一区二区三区免费视频棣 | 亚洲精品无码av牛牛影视| 国产视频一区二区三区四区| 亚洲三级网| 国产精品久久久久久模特 | 国产精品一二| 国产精品国产三级国产普通话一| 国产三级日本三级在线播放| 五月丁香综合| 欧美日本在线观看| 国产操b| 亚洲黄色片免费看| 国产精品久久久久久久久久| 亚洲精品91| 亚洲AV性爱电影| 日韩黄色视屏| 久久久五月天| 天天综合久久综合| av无码在线不卡| 国产1区2区3区| 91久久久久国产一区二区| 无码天堂| 国产精品毛片一区二区在线看| 91福利网| 特级毛片绝黄A片免费播冫| 久久久国产精品视频| 欧美一级黄片免费观看| 国产精品伦一区二区三级视频| 91网站入口| 成人精品视频在线| 屁屁影院第一页| 人人妻人人澡人人爽精品日本| 无码国产一区二区三区| 日韩成人片在线观看| 午夜成人在线| 人妖一区二区| 欧美视频第二页| 宅男午夜影院| 99国产一区| 无码在线电影| 亚洲图片欧美另类| 久久久久久伊人| 国产特级片| 日本无码高清| 久久精品一区二区三区四区| 亚洲综合无码| 强奸乱伦亚洲综合| blacked精品一区国产99| 日韩性爱免费网| 亚洲熟女乱伦| 免费无码又爽又黄又刺激网站| 99久久99| 国产一区不卡在线| 婷婷精品视频| 成人片在线观看| 精品国产乱码久久久久电车痴汉久| 欧美91| 五月丁香中文字幕| 自拍偷拍av| 综合色av| 无码国产精品一区| 亚洲丰满少妇在线播放| 免费黄色大片| A级黄片免费视频| 91精品人妻| 日韩欧美色图| 国产婷婷色| 在线看91| 三级久久| 久久亚洲精品成人AV| 国产精品毛片一区二区在线看| 日韩无码免费视频| 国产精品爱久久久久久久威尼斯| 91最新在线视频| 成人欧美一区二区三区黑人免费| 97中文字幕在线观看| 日韩欧美国产精品| 国产欧美一区二区精品97| 成人A片无码水蜜桃免费网站软件| 亚洲AV导航| 国产精品一区二区视频| 国产精品固产视频| 99久久99久久精品国产片果冻 | 色网在线| 女人18片毛片90分钟| 日韩激情网| 亚洲无码在线一区| 精品少妇人妻| 无码专区AV| 亚洲熟女一区| 五月婷婷综合| 成人黄色一级视频| 欧美91视频| 国产成人在线视频播放| 91在线精品一区二区三区| 亚洲AV永久纯肉无码精品动漫| 色天堂在线| 91精品在线视频观看| 丁香五月v国产| 久久精品国产亚洲AV无码娇色 | 国产做a爰片毛片A片美国| 国产精品久久欧美久久一区| 亚洲欧美乱伦| 日韩欧美国产综合| 91久久国产综合| 天天射天天干天天日| 国产露脸91国语对白| 人妻视频在线| 五月天狠狠爱| 国产精品成人一区二区网站软件 | 91免费在线播放| 亚洲熟女乱综合一区二区三区| 深山熟女Av| 欧美日本一区二区三区| 亚洲熟女久久| www..com操老师| 九九热视频在线| 男女高潮又爽又黄又无遮挡| 一级香蕉,黄色片| 91人妻中文字幕在线精品| 国产小视频91| 草视频黄在线| 亚洲熟妇无码AV无码| 99久久久国产精品免费蜜臀| 视频国产精品| 色色色婷婷| 欧美精品一区二区三区| 久久精彩视频| 一级亚洲| 欧美日韩偷拍视频| 91无码偷拍精品一区二区三区| 一本大道无码| 国产美女视频| 免费一级全黄少妇性色生活片| 亚洲精品免费在线观看| 欧美一级性爱视频| 国产手机视频在线| 久久人人爽人人| 日本免费高清视频| 第一福利视频导航| 欧美精品第一页| 亚洲AV导航| 丝袜灬啊灬快灬高潮了AV| 产国传媒91一区久久无码| 亚洲视频在线一区二区| 成人高清无码| 色悠悠在线| 亚洲1区2区| 国产在线观看一区二区| 91人人操| 无码人妻精品一区二区中文| 牲欲强的熟妇农村老妇女视频| 西欧毛片| 永久免费黄片| 老司机午夜福利视频| 女人高潮抽搐喷液30分钟视频| 少妇又紧又色又爽又刺激视频 | 亚洲国产日韩三级av探花| 久久精品噜噜噜成人| 久久e热| 人人弄人人摸| 翔田千里在线播放AV101| 久久九九精品99国产精品| 无码国产精品一区二区免费网站| 日韩无码性爱视频| 家庭乱伦网站国产| 在线观看av的网站| 三级片免费网址| 亚洲视频入口| 国产成人网站在线观看| 友田真希一区| 久久精品久久久久久久| 久久一区二区视频| 日韩美女一区二区三区| blacked精品一区国产99| 91久久久精品国产一区二区爱豆 | 黄色爱爱视频| 视频一区在线播放| 国产毛片久久久久| 精东粉嫩av免费一区二区三区| 国产成人精品一区二区三区 | 日韩av电影在线观看| 91sese| 亚洲无码在线一区| 成人av免费在线观看| 最新av导航| 伊人久久一区| 激情婷婷五月天| 这里只有精品视频| 97人妻蜜臀中文字幕| 九九精品久久| 后入内射欧美99二区视频| 久久久久一区| 久久国产精品久久久| av水蜜桃| 啄木乌欧美一区二区三区| 韩国无码在线观看| 精品人妻一区二区| 国产精品福利在线| 成人午夜在线| 国产三级网站| 国产成人三区| 亚洲AV成人无码久久精品| 一级片免费视频| xxxx18一20岁hd| 青青操免费在线视频| 国产91视频|