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

    非牛頓流體中在線氣泡聚并的混沌特征

    2010-05-10 06:26:42姜韶堃馬友光范文元朱春英
    關(guān)鍵詞:牛頓流體氣速相空間

    姜韶堃 ,馬友光,范文元,楊 珂,朱春英

    (天津大學(xué)化工學(xué)院化學(xué)工程聯(lián)合國(guó)家重點(diǎn)實(shí)驗(yàn)室,天津 300072)

    非牛頓流體中在線氣泡聚并的混沌特征

    姜韶堃 ,馬友光,范文元,楊 珂,朱春英

    (天津大學(xué)化工學(xué)院化學(xué)工程聯(lián)合國(guó)家重點(diǎn)實(shí)驗(yàn)室,天津 300072)

    利用光電信號(hào)采集系統(tǒng)測(cè)定了非牛頓流體中伴有聚并現(xiàn)象的在線連續(xù)上升氣泡的通過(guò)頻率信號(hào),并采用基于相空間重構(gòu)技術(shù)的混沌時(shí)間序列法對(duì)其進(jìn)行了分析.相圖和 Cao算法計(jì)算結(jié)果定性表明了非牛頓流體中伴有聚并現(xiàn)象的在線連續(xù)氣泡上升行為具有混沌吸引子特征,服從確定性的混沌運(yùn)動(dòng)機(jī)理.最大 Lyapunov指數(shù)量化分析表明在線連續(xù)氣泡上升行為的混沌特征與聚并程度密切相關(guān):氣泡行為的混沌程度隨聚并程度增大而增大,增大氣速易于導(dǎo)致氣泡行為混沌程度增大.

    氣泡;聚并;非牛頓流體;相空間重構(gòu);最大Lyapunov指數(shù)

    非牛頓流體中的氣液兩相流在氣液接觸、氣液分離、氣體吸收、沸騰、鼓泡塔、發(fā)酵和氣蝕等眾多領(lǐng)域廣泛存在.由于非牛頓流體的流變特性,其中的氣泡行為比在牛頓流體中更加復(fù)雜,例如氣泡間的聚并過(guò)程.關(guān)于非牛頓流體中氣泡聚并現(xiàn)象的研究報(bào)道較少.Acharya等[1]和 De Kee等[2-3]采用高速攝像機(jī)對(duì)黏彈性流體中氣泡形狀及氣泡聚并過(guò)程進(jìn)行了圖像采集和實(shí)驗(yàn)研究.其后,Li等[4-5]和 Frank等[6]采用混沌理論分析了黏彈性流體中伴有聚并現(xiàn)象的氣泡上升過(guò)程,得到其混沌運(yùn)動(dòng)特征.此外,Lin等[7-9]通過(guò)粒子圖像分析儀測(cè)量了黏彈性流體中氣泡聚并過(guò)程時(shí)的流線和剪切應(yīng)力分布,分析了剪切變稀效應(yīng)和黏彈性在聚并過(guò)程的作用.迄今,人們對(duì)非牛頓流體中的氣泡聚并機(jī)理的認(rèn)識(shí)還不是很清楚,更難以進(jìn)行精確的理論描述.

    對(duì)于非牛頓流體中在線連續(xù)氣泡上升行為,在某一固定高度,由于氣泡的不斷聚并,連續(xù)氣泡上升行為表現(xiàn)出明顯的非周期特征,因此很難建立有效的理論預(yù)測(cè)模型.采用混沌時(shí)間序列分析研究非牛頓流體中在線連續(xù)氣泡上升行為,是一種獲得系統(tǒng)非線性動(dòng)力學(xué)特性的有效途徑.理論上,混沌是指確定的非周期行為,對(duì)初始條件具有很強(qiáng)的敏感依賴性.一般利用相空間重構(gòu)方法,通過(guò)相圖、龐加萊截面、關(guān)聯(lián)維數(shù)、Kolmogrov熵和最大 Lyapunov指數(shù)等重要指標(biāo)來(lái)定性和定量判斷系統(tǒng)混沌行為.筆者主要分析了非牛頓流體中伴有聚并現(xiàn)象的在線連續(xù)氣泡上升行為,鑒于過(guò)程的復(fù)雜性,采用不干擾流體流動(dòng)的無(wú)接觸式光電信號(hào)系統(tǒng)采集氣泡通過(guò)頻率信號(hào),以相空間重構(gòu)技術(shù)為基礎(chǔ),利用相圖、最大Lyapunov指數(shù)等混沌時(shí)間序列分析指標(biāo)進(jìn)行了研究,討論了非牛頓流體中伴有聚并現(xiàn)象的在線連續(xù)氣泡上升行為的非線性混沌特征.

    1 實(shí) 驗(yàn)

    實(shí)驗(yàn)采用底部中央帶有一個(gè)噴嘴(直徑:1×10-3m)的方形有機(jī)玻璃槽(長(zhǎng)×寬×高:0.15,m×0.15,m×1.50,m)作為氣泡發(fā)生裝置,槽內(nèi)液面距噴嘴高度恒定為 1.10,m,如圖 1所示,裝置下方連接一個(gè)較大氣室來(lái)抑制氣泡生成過(guò)程引起的壓力波動(dòng).調(diào)節(jié)轉(zhuǎn)子流量計(jì)至適當(dāng)氣速,氮?dú)馀萃ㄟ^(guò)噴嘴在槽內(nèi)非牛頓流體中連續(xù)生成并不斷上升.采用槽外光電信號(hào)采集系統(tǒng)對(duì)在線連續(xù)上升氣泡進(jìn)行信號(hào)捕捉,該系統(tǒng)主要由激光器、光電二極管、運(yùn)算放大器以及數(shù)據(jù)采集卡組成,采樣頻率為 200,Hz.激光光束垂直于玻璃槽一側(cè)壁面入射,穿過(guò)氣泡上升軌跡,當(dāng)氣泡通過(guò)激光光束時(shí),由于氣泡表面的反射,只有部分光線通過(guò)氣泡,位于玻璃槽另一側(cè)的光電二極管接收該光強(qiáng)變化,并以電壓形式表現(xiàn)出來(lái),從而形成相應(yīng)的脈沖信號(hào).所有實(shí)驗(yàn)均在恒定室溫297,K下完成.

    圖1 氣泡聚并實(shí)驗(yàn)裝置示意Fig.1 Schematic representation of experimental apparatus Fig.1 for bubble coalescence

    實(shí)驗(yàn)中采用的非牛頓流體的質(zhì)量分?jǐn)?shù)分別為1.0%、1.3%和1.6%的羧甲基纖維素鈉(carboxymethylcellulose sodium,CMC)水溶液,其流變性質(zhì)由 DV-III ULTRA PROGRAMMABLE RHEOMETER (Brookfield Engineering Laboratories,Inc.,U.S.A.)測(cè)定.從測(cè)定結(jié)果(見(jiàn)圖2)可以看出,3種CMC水溶液的黏度均隨著剪切速率的升高而降低,表現(xiàn)出明顯的剪切變稀效應(yīng).

    圖2 黏度隨剪切速率變化Fig.2 Variation of viscosity vs shear rate

    2 混沌時(shí)間序列分析

    在 3種 CMC水溶液中,氣泡運(yùn)動(dòng)行為十分相似,因此僅以 1.0% CMC(質(zhì)量分?jǐn)?shù))的實(shí)驗(yàn)結(jié)果為例進(jìn)行分析.

    首先采用雙光束法和照相法分別估測(cè)了氣泡上升速度和氣泡體積.結(jié)果顯示,氣泡上升速度范圍為0.20~ 0.40,m/s,氣 泡 體 積 范 圍 為 0.52×10-6~27.16×10-6,m3.按照文獻(xiàn)[10]對(duì)氣泡剪切速率的描述

    式中:v表示氣泡上升速度;deq表示氣泡當(dāng)量直徑;VB為氣泡體積.計(jì)算得到氣泡上升剪切速率范圍為10.54~25,s-1,恰好位于剪切變稀速率范圍內(nèi).

    圖 3表示了不同位置的氣泡通過(guò)頻率信號(hào)片段.各個(gè)片段中都存在不同的信號(hào)脈沖,每一個(gè)脈沖都對(duì)應(yīng)于一個(gè)單獨(dú)通過(guò)的氣泡.在噴嘴上方較低的位置,信號(hào)脈沖均勻密集,表示氣泡的周期性上升過(guò)程;隨著位置的升高,由于氣泡聚并的發(fā)生,逐漸產(chǎn)生無(wú)規(guī)則的脈沖分布;當(dāng)位置接近液面時(shí),聚并頻率大大降低,單位時(shí)間內(nèi)氣泡通過(guò)數(shù)量也明顯減少,信號(hào)脈沖變得稀疏無(wú)序.

    采用混沌時(shí)間序列分析方法來(lái)研究氣泡通過(guò)頻率信號(hào)的非線性動(dòng)力學(xué)特性,最關(guān)鍵的步驟就是相空間重構(gòu)參數(shù)的選擇,包括時(shí)間延遲τ 和嵌入維數(shù)m.以圖3中h=0.3,m時(shí)的中信號(hào)為例,將實(shí)際測(cè)得的氣泡通過(guò)頻率信號(hào)歸一化為均值為 0振幅為 1的時(shí)間序列xi,i=1,2,…,N.圖4顯示出采用不同的τ對(duì)時(shí)間序列xi重構(gòu)得到的二維相軌跡演變.當(dāng)τ選擇過(guò)小時(shí),如圖 4(a)所示,相軌跡都集中于對(duì)角線附近,數(shù)據(jù)之間關(guān)聯(lián)程度太強(qiáng),不易區(qū)分;當(dāng) τ選擇過(guò)大時(shí),如圖 4(b)所示,相軌跡呈現(xiàn)出明顯的折疊現(xiàn)象,數(shù)據(jù)之間毫不相關(guān).因此,必須選擇適當(dāng)?shù)?τ,使得重構(gòu)的相空間既不線性相關(guān),又非完全獨(dú)立.去偏復(fù)自相關(guān)法[11]在實(shí)際信號(hào)處理上具有較好的魯棒性,筆者采用該法計(jì)算τ.時(shí)間序列xi的m維去偏復(fù)自相關(guān)法存在基本關(guān)系為

    圖3 氣泡通過(guò)頻率信號(hào)片段Fig.3 Part of frequency signal of bubble passage

    圖4 不同時(shí)間延遲的相圖Fig.4 Phase portraits with different time delays

    圖5 去偏復(fù)自相關(guān)法計(jì)算結(jié)果Fig.5 Results of nonbias multiple autocorrelation method

    對(duì)于嵌入維數(shù) m的選擇,已有大量文獻(xiàn)報(bào)道采用 Cao算法[12-15]能夠準(zhǔn)確獲得許多實(shí)際信號(hào)的m.其基本計(jì)算方法為

    圖6 不同位置的相圖Fig.6 Phase portraits of different positions

    式中:Xi(m)表示實(shí)測(cè)氣泡通過(guò)頻率信號(hào)經(jīng)相空間重構(gòu)得到的向量(xi,xi+τ,…,xi+(m-1)τ),i=1,2,…,Nm;下標(biāo) n(i,m)(1≤n(i,m)≤N-mτ)使得 Xn(i,m)(m)是在按上述距離定義的 m維重構(gòu)相空間中 Xi(m)的最近鄰點(diǎn).隨著m的不斷增大,當(dāng)m增大到某一確定值時(shí),E1(m)基本保持恒定,這時(shí)的 m 即為嵌入維數(shù) m;若m 由小到大變化,E2(m)始終保持在 1左右,那么信號(hào)就是隨機(jī)的;反之,就是確定性的.圖7為采用Cao算法對(duì)氣泡通過(guò)頻率信號(hào)的計(jì)算結(jié)果.從圖中可以看出,E2隨著 m的增大并不是始終為 1,由此表明擬塑性流體中的氣泡連續(xù)上升過(guò)程是確定性的,而非隨機(jī)運(yùn)動(dòng);此外,當(dāng)E1趨于收斂時(shí)所對(duì)應(yīng)的m值均比較大,而 m標(biāo)志著能夠足夠描述相軌跡狀態(tài)所需的最小相空間維數(shù),即高氣速條件下 CMC水溶液中氣泡運(yùn)動(dòng)行為的系統(tǒng)自由度也比較大,與文獻(xiàn)[16]中黏彈性流體低氣速條件下的結(jié)果相比,氣泡行為更加復(fù)雜.

    圖7 Cao算法計(jì)算結(jié)果Fig.7 Results of Cao′s method

    當(dāng)τ和m確定后,即可通過(guò)Lyapunov指數(shù)量化氣泡通過(guò)頻率信號(hào)的混沌特性.在許多實(shí)際信號(hào)處理過(guò)程中,采用 Wolf法[17]能夠有效地獲得最大Lyapunov指數(shù).首先在重構(gòu)的相空間中選中初始基點(diǎn) X(t0)及其最近鄰點(diǎn) X(t0′),隨著時(shí)間演化,當(dāng)?shù)竭_(dá)t1時(shí)刻,兩點(diǎn)間距離由 L0變化為 L0′,并且L0′>ε(ε>0),保留基點(diǎn) X(t1),重新選取與之夾角盡可能小的最近鄰點(diǎn) X(t1′),使得兩點(diǎn)之間距離 L1<ε,繼續(xù)上述演化跟蹤過(guò)程,直到相軌跡終點(diǎn),此時(shí)演化跟蹤過(guò)程迭代次數(shù)為 M,則最大 Lyapunov指數(shù)λ1可表示為一般只要λ1大于 0,就可判斷系統(tǒng)處于混沌狀態(tài);而對(duì)于周期運(yùn)動(dòng),λ1小于或等于0.

    表 1列出了不同氣速條件下λ1隨高度的分布.對(duì)于如Qg=5×10-6m3/s、h=0.1 m時(shí)的氣泡運(yùn)動(dòng)行為,位置較低,氣速較小,氣泡生成間隔較大,氣泡間存在線性相互作用,形成較均勻的氣泡上升過(guò)程,因此λ1較小,接近于 0;隨著氣速增大,氣泡生成間隔變小,氣泡間相互作用逐漸由線性轉(zhuǎn)變?yōu)榉蔷€性,破壞了氣泡上升的周期性規(guī)律,形成非周期性上升過(guò)程,λ1明顯增大.隨著位置升高,不同氣速條件下的氣泡均加速上升,CMC水溶液的剪切變稀效應(yīng)開(kāi)始逐漸影響氣泡行為:氣泡通過(guò)后,局部黏度降低,使得下方氣泡在局部阻力減小的情況下加速上升.受剪切變稀效應(yīng)和先行氣泡尾流等因素共同影響,氣泡行為仍表現(xiàn)為非周期上升特征,甚至當(dāng)位置升高到某一高度后,發(fā)生氣泡聚并現(xiàn)象,氣泡行為非周期特征加劇,λ1則表現(xiàn)出隨著位置升高而增大,當(dāng)位置繼續(xù)升高,氣泡頻繁聚并,λ1達(dá)到最大值;隨后,由于先前頻繁聚并生成體積較大的新氣泡,嚴(yán)重破壞了氣泡原有的受力分布,使得氣泡間非線性相互作用減弱,聚并頻率大大減小,λ1明顯減小.此外,隨著氣速增大,氣泡生成間隔變小,上升速度增大,剪切速率增大,流體的剪切變稀效應(yīng)促使氣泡更易加速上升,在隨后的上升過(guò)程中,更易強(qiáng)化非線性相互作用,導(dǎo)致聚并程度增大,從表 1中可以看出,氣速增大更易得到較大的λ1.通常,由λ1的數(shù)值可以定量判斷系統(tǒng)運(yùn)動(dòng)的混沌程度,表中 1所有條件下的λ1均大于 0,結(jié)合相圖和Cao算法等定性分析,CMC水溶液中在線氣泡上升過(guò)程表現(xiàn)為確定性的混沌運(yùn)動(dòng),頻繁聚并現(xiàn)象使得氣泡行為混沌程度增大,而增大氣速也易造成氣泡運(yùn)動(dòng)混沌程度增大.

    表1 不同氣速條件下λ1隨高度分布Tab.1 Variation of λ1 with height at different gas flowrates

    3 結(jié) 論

    利用光電信號(hào)采集系統(tǒng)測(cè)量了非牛頓流體 CMC水溶液中在線連續(xù)上升氣泡的通過(guò)頻率信號(hào),采用基于相空間重構(gòu)技術(shù)的混沌時(shí)間序列分析方法對(duì)該信號(hào)進(jìn)行了分析,得到如下結(jié)論:

    (1) 相空間重構(gòu)二維相圖清晰地顯示,伴有聚并現(xiàn)象的在線氣泡連續(xù)上升行為表現(xiàn)出明顯的混沌吸引子特征.

    (2) Cao算法計(jì)算結(jié)果表明,CMC水溶液中的在線氣泡上升過(guò)程具有較高的系統(tǒng)自由度,是一個(gè)復(fù)雜的確定性非隨機(jī)過(guò)程.

    (3) 最大 Lyapunov指數(shù)均大于 0,表明了 CMC水溶液中在線氣泡連續(xù)上升行為的混沌特征.

    (4) 由最大Lyapunov指數(shù)定量描述的氣泡行為混沌程度與氣泡聚并程度密切相關(guān),前者隨后者增大而增大,相同高度下,增大氣速導(dǎo)致氣泡行為混沌程度增大.

    [1]Acharya A,Ulbrecht J. Note on the influence of viscoelasticity on the coalescence rate of bubbles and drops[J].AIChE Journal,1978,24(2):348-351.

    [2]De Kee D,Chhabra R P,Dajan A. Motion and coalescence of gas bubbles in non-Newtonian polymer solutions[J].Journal of Non-Newtonian Fluid Mechanics,1990,37(1):1-18.

    [3]Dekee D,Chan M,F(xiàn)ong C F,et al. Bubble shape in non-Newtonian fluids[J].Journal of Applied Mechanics,2002,69(5) :703-708.

    [4]Li H Z,F(xiàn)rank X,F(xiàn)unfschilling D,et al. Bubbles′ rising dynamics in polymeric solutions[J].Physics Letter A,2004,325(1):43-50.

    [5]Li H Z. Bubbles in non-Newtonian fluids:Formation,interactions and coalescence[J].Chemical Engineering Science,1999,54(13/14):2247-2254.

    [6]Frank X,Li H Z. Route to chaos in the rising dynamics of a bubble chain in a polymeric fluid[J].Physics Letter A,2008,372(40):6155-6160.

    [7]Lin T J,Lin G M. The mechanisms of bubble coalescence in a non-newtonian fluid[J].The Canadian Journal of Chemical Engineering,2003,81(3/4):476-482.

    [8]Lin T J,Lin G M. An experimental study on flow structures of a single bubble rising in a shear-thinning viscoe-lastic fluid with a new measurement technique[J].International Journal of Multiphase Flow,2005,31(2):239-252.

    [9]Lin T J,Lin G M. Mechanisms of in-line coalescence of two-unequal bubbles in a non-Newtonian fluid[J].Chemical Engineering Journal,2009,15(3):750-756.

    [10]林嘉宇,王躍科,黃芝平,等. 語(yǔ)音信號(hào)相空間重構(gòu)中時(shí)間延遲的選擇:復(fù)自相關(guān)法[J]. 信號(hào)處理,1999,15(3):220-225.

    Lin Jiayu,Wang Yueke,Huang Zhiping,et al. Selection of proper time-delay in phase space reconstruction of speech signals[J].Signal Processing,1999,15(3):220-225(in Chinese).

    [11]Ma H G,Han C Z. Selection of embedding dimension and delay time in phase space reconstruction[J].Frontiers of Electrical and Electronic Engineering in China,2006,1(1):111-114.

    [12]Cao Liangyue. Practical method for determining the minimum embedding dimension of a scalar time series[J].Physica D:Nonlinear Phenomena,1997,110(1/2):43-50.

    [13]Kumar U,Prakasha A,Jain V K. Characterization of chaos in air pollutants:A Volterra-Wiener-Korenberg series and numerical titration approach[J].Atmospheric Environment,2008,42(7):1537-1551.

    [14]Wang J,Sun L,F(xiàn)ei X Y,et al. Chaos analysis of the electrical signal time series evoked by acupuncture[J].Chaos,Solitons and Fractals,2007,33(3):901-907.

    [15]Li X L,Ouyang G. Nonlinear similarity analysis for epileptic seizures prediction[J].Nonlinear Analysis:Theory,Methods and Applications,2006,64(8):1666-1678.

    [16]Li H Z,Mouline Y,Choplin L,et al. Chaotic bubble coalescence in non-Newtonian fluids[J].International Journal of Multiphase Flow,1997,23(4):713-723.

    [17]Wolf A,Swift J B,Swinney H L,et al. Determining Lyapunov exponents from a time series[J].Physica D:Nonlinear Phenomena,1985,16(3):285-317.

    Chaotic Characteristics of in Line Bubble Coalescence in Non-Newtonian Fluid

    JIANG Shao-kun,MA You-guang,F(xiàn)AN Wen-yuan,YANG Ke,ZHU Chun-ying
    (State Key Laboratory of Chemical Engineering,School of Chemical Engineering and Technology,Tianjin University,Tianjin 300072,China)

    The frequency signal of continuous rising bubbles with coalescences in non-Newtonian fluid was measured by means of a photoelectrical signal collection system and analyzed using chaotic time series analysis based on the phase space reconstruction. The results of phase portrait and Cao's method show qualitatively that the in-line continuous bubble rising behavior with coalescence in non-Newtonianfluid follows a deterministic and chaotic motion mechanism with the characteristics of chaotic attractor. The largest Lyapunov exponents quantify that the chaotic characteristics of the in-line continuous bubble rising behavior are closely related with the coalescence degree. The chaotic degree of bubble behavior increases with the coalescence degree and increasing the gas flowrate tends to increase the chaotic degree of bubble behavior.

    bubble;coalescence;non-Nnewtonian fluid;phase space reconstruction;largest Lyapunov exponent

    TQ021

    A

    0493-2137(2010)11-1025-06

    2009-08-24;

    2010-06-23.

    國(guó)家自然科學(xué)基金資助項(xiàng)目(20476073);化學(xué)工程聯(lián)合國(guó)家重點(diǎn)實(shí)驗(yàn)室課題(SKL-ChE-08B03).

    堃姜韶 (1980— ),男,博士,j_sk45@yahoo.com.cn.

    馬友光,ygma@tju.edu.cn.

    猜你喜歡
    牛頓流體氣速相空間
    束團(tuán)相空間分布重建技術(shù)在西安200 MeV質(zhì)子應(yīng)用裝置的應(yīng)用
    非牛頓流體
    什么是非牛頓流體
    少兒科技(2019年3期)2019-09-10 07:22:44
    區(qū)別牛頓流體和非牛頓流體
    濕法煙氣脫硫吸收塔阻力特性實(shí)驗(yàn)研究
    新型折板除霧器的流場(chǎng)和壓降數(shù)值模擬
    首款XGEL非牛頓流體“高樂(lè)高”系列水溶肥問(wèn)世
    非對(duì)易空間中的三維諧振子Wigner函數(shù)
    基于相空間重構(gòu)的電磁繼電器電性能參數(shù)預(yù)測(cè)研究
    D120 mm流化床冷模實(shí)驗(yàn)研究
    化工科技(2014年5期)2014-06-09 05:17:22
    a级毛片免费高清观看在线播放| 新久久久久国产一级毛片| 欧美一级a爱片免费观看看| 岛国毛片在线播放| 国产 一区精品| 有码 亚洲区| 热99国产精品久久久久久7| 国产精品蜜桃在线观看| 亚洲性久久影院| 两个人的视频大全免费| 中文精品一卡2卡3卡4更新| av国产久精品久网站免费入址| 如何舔出高潮| 伦精品一区二区三区| 免费看光身美女| 国产极品天堂在线| 日韩强制内射视频| 少妇人妻精品综合一区二区| 久久久欧美国产精品| 亚洲第一区二区三区不卡| 国产一区二区在线观看av| av国产精品久久久久影院| 麻豆精品久久久久久蜜桃| 欧美变态另类bdsm刘玥| 男的添女的下面高潮视频| 亚洲精品色激情综合| 亚洲色图综合在线观看| 99久国产av精品国产电影| 我要看黄色一级片免费的| 国语对白做爰xxxⅹ性视频网站| 美女内射精品一级片tv| 久久久久网色| 成年人午夜在线观看视频| 国产中年淑女户外野战色| 99热6这里只有精品| 亚洲精品中文字幕在线视频 | 国产免费又黄又爽又色| 国产伦精品一区二区三区四那| 久久国产精品男人的天堂亚洲 | 久久 成人 亚洲| 99热6这里只有精品| 成人午夜精彩视频在线观看| 国产精品国产三级专区第一集| 欧美少妇被猛烈插入视频| 国产亚洲精品久久久com| 伊人久久国产一区二区| 大陆偷拍与自拍| 搡女人真爽免费视频火全软件| 久久精品久久精品一区二区三区| 内地一区二区视频在线| 国产熟女欧美一区二区| 少妇的逼好多水| 久久久久网色| 国产精品不卡视频一区二区| 日本色播在线视频| 热re99久久国产66热| 免费黄网站久久成人精品| 伦理电影大哥的女人| av福利片在线| 六月丁香七月| 欧美一级a爱片免费观看看| 大话2 男鬼变身卡| 亚洲欧洲日产国产| 国产精品一区www在线观看| 欧美97在线视频| 两个人的视频大全免费| 丰满迷人的少妇在线观看| 日韩一区二区三区影片| 黄色配什么色好看| 国产精品国产三级国产av玫瑰| 人妻人人澡人人爽人人| 男女无遮挡免费网站观看| 永久免费av网站大全| 又粗又硬又长又爽又黄的视频| 老司机亚洲免费影院| 夜夜爽夜夜爽视频| 蜜臀久久99精品久久宅男| av天堂久久9| 久久精品夜色国产| 国产伦在线观看视频一区| 久久99蜜桃精品久久| 亚洲情色 制服丝袜| 亚洲av综合色区一区| 内射极品少妇av片p| av免费观看日本| 蜜臀久久99精品久久宅男| 久久人人爽人人爽人人片va| 乱系列少妇在线播放| 久久久精品免费免费高清| a级毛片在线看网站| 国产视频首页在线观看| 国产无遮挡羞羞视频在线观看| 一个人看视频在线观看www免费| 午夜av观看不卡| 免费观看无遮挡的男女| 久久久久久久亚洲中文字幕| 国产女主播在线喷水免费视频网站| 国产亚洲午夜精品一区二区久久| 国产中年淑女户外野战色| 国产毛片在线视频| 九九久久精品国产亚洲av麻豆| av国产久精品久网站免费入址| 免费黄频网站在线观看国产| 91午夜精品亚洲一区二区三区| 我要看日韩黄色一级片| 蜜桃久久精品国产亚洲av| 亚洲成人av在线免费| 成人国产av品久久久| 一级a做视频免费观看| 大片电影免费在线观看免费| 性色avwww在线观看| 欧美 日韩 精品 国产| 国产成人午夜福利电影在线观看| 国产精品一区二区性色av| 久久久久久久久大av| 亚洲高清免费不卡视频| 九九爱精品视频在线观看| www.av在线官网国产| 街头女战士在线观看网站| 哪个播放器可以免费观看大片| xxx大片免费视频| 久久6这里有精品| 国产精品99久久99久久久不卡 | 热99国产精品久久久久久7| 熟女人妻精品中文字幕| 黄色欧美视频在线观看| 日本免费在线观看一区| av福利片在线| 午夜免费鲁丝| 狂野欧美激情性bbbbbb| 久久毛片免费看一区二区三区| 91久久精品电影网| 黑人猛操日本美女一级片| 曰老女人黄片| 黄色欧美视频在线观看| 天天操日日干夜夜撸| 亚洲av.av天堂| 亚洲,一卡二卡三卡| 亚洲国产欧美日韩在线播放 | 久久久久人妻精品一区果冻| 亚洲四区av| 观看免费一级毛片| 女人久久www免费人成看片| 亚洲一区二区三区欧美精品| 亚洲国产av新网站| 中国三级夫妇交换| 免费大片18禁| 一级二级三级毛片免费看| 日本爱情动作片www.在线观看| 日韩,欧美,国产一区二区三区| 午夜激情福利司机影院| 六月丁香七月| 久久人人爽人人片av| 99热这里只有是精品50| 精品少妇黑人巨大在线播放| av福利片在线观看| 天天躁夜夜躁狠狠久久av| 亚洲国产av新网站| 国产高清国产精品国产三级| 91在线精品国自产拍蜜月| 在线免费观看不下载黄p国产| 亚洲国产欧美日韩在线播放 | 免费观看a级毛片全部| 九色成人免费人妻av| 亚洲国产最新在线播放| 午夜免费观看性视频| 一区在线观看完整版| 国产在视频线精品| 亚洲久久久国产精品| 亚洲国产av新网站| 妹子高潮喷水视频| 边亲边吃奶的免费视频| 少妇猛男粗大的猛烈进出视频| 亚洲三级黄色毛片| 99热国产这里只有精品6| 亚洲精品国产av成人精品| av视频免费观看在线观看| 午夜av观看不卡| 春色校园在线视频观看| 麻豆成人av视频| 亚洲国产精品专区欧美| 97精品久久久久久久久久精品| 五月开心婷婷网| 国产精品人妻久久久影院| 日本黄色片子视频| 成人毛片a级毛片在线播放| 91在线精品国自产拍蜜月| 亚洲av成人精品一区久久| 日产精品乱码卡一卡2卡三| 美女福利国产在线| 在线播放无遮挡| 日本vs欧美在线观看视频 | 日本午夜av视频| 一本久久精品| 亚洲,欧美,日韩| 亚洲av男天堂| 欧美成人精品欧美一级黄| 亚洲电影在线观看av| 国产精品一区二区在线不卡| 久久国产精品男人的天堂亚洲 | av网站免费在线观看视频| 啦啦啦中文免费视频观看日本| 国产高清三级在线| 久久久久精品久久久久真实原创| 在线免费观看不下载黄p国产| 亚洲成色77777| √禁漫天堂资源中文www| 亚洲国产成人一精品久久久| 国产亚洲午夜精品一区二区久久| 亚洲精品456在线播放app| 天堂中文最新版在线下载| 国产日韩一区二区三区精品不卡 | 丝袜喷水一区| 午夜av观看不卡| 老司机亚洲免费影院| 国产在线视频一区二区| av天堂久久9| 乱码一卡2卡4卡精品| 2022亚洲国产成人精品| 日本欧美国产在线视频| 伊人久久精品亚洲午夜| 久久午夜综合久久蜜桃| 国产欧美日韩精品一区二区| 国产成人免费观看mmmm| 纯流量卡能插随身wifi吗| 美女大奶头黄色视频| 啦啦啦视频在线资源免费观看| 国产伦精品一区二区三区四那| 国产伦精品一区二区三区视频9| 亚洲欧美日韩另类电影网站| 天堂中文最新版在线下载| 久久青草综合色| 男人添女人高潮全过程视频| 久久久久精品久久久久真实原创| 亚洲av综合色区一区| 成人毛片60女人毛片免费| 国产av国产精品国产| 日韩电影二区| 18禁裸乳无遮挡动漫免费视频| 性色av一级| 亚洲电影在线观看av| 最近中文字幕高清免费大全6| 国产免费又黄又爽又色| 九九久久精品国产亚洲av麻豆| 国产视频内射| 久久久亚洲精品成人影院| av国产久精品久网站免费入址| 亚洲性久久影院| 久久人人爽人人爽人人片va| 内地一区二区视频在线| 观看av在线不卡| 欧美xxxx性猛交bbbb| 久久国产精品大桥未久av | 夫妻性生交免费视频一级片| 精品少妇久久久久久888优播| 成年女人在线观看亚洲视频| 丰满人妻一区二区三区视频av| 亚洲精品456在线播放app| 欧美精品人与动牲交sv欧美| 国语对白做爰xxxⅹ性视频网站| 亚洲在久久综合| 免费av中文字幕在线| 久久精品国产a三级三级三级| 亚洲真实伦在线观看| 夜夜骑夜夜射夜夜干| 国产在线视频一区二区| 天堂8中文在线网| 国产伦精品一区二区三区四那| 99热网站在线观看| 久久精品国产亚洲av天美| 亚洲精品久久久久久婷婷小说| 一本大道久久a久久精品| av有码第一页| 黄色怎么调成土黄色| 亚洲丝袜综合中文字幕| 99久久中文字幕三级久久日本| 热re99久久精品国产66热6| 在线看a的网站| 精品少妇内射三级| 免费大片黄手机在线观看| 日韩强制内射视频| 久久毛片免费看一区二区三区| 亚洲精品色激情综合| 亚洲精品乱码久久久v下载方式| 欧美精品一区二区免费开放| 狂野欧美激情性bbbbbb| 久久热精品热| 亚洲欧美中文字幕日韩二区| 大话2 男鬼变身卡| 热99国产精品久久久久久7| 成人亚洲精品一区在线观看| 久久青草综合色| 91精品一卡2卡3卡4卡| 自拍欧美九色日韩亚洲蝌蚪91 | 日韩av在线免费看完整版不卡| 欧美另类一区| 99九九在线精品视频 | 男人舔奶头视频| 国产欧美另类精品又又久久亚洲欧美| 亚洲婷婷狠狠爱综合网| 极品教师在线视频| 国产成人精品福利久久| 亚洲精品日韩av片在线观看| 久久韩国三级中文字幕| 国产高清三级在线| 黑丝袜美女国产一区| 久久av网站| 日本vs欧美在线观看视频 | av视频免费观看在线观看| 美女主播在线视频| 国产亚洲精品久久久com| 亚洲久久久国产精品| 视频区图区小说| 国产在视频线精品| 国产成人freesex在线| av福利片在线| 女性被躁到高潮视频| 日本黄色日本黄色录像| 欧美bdsm另类| 人妻人人澡人人爽人人| 日韩成人伦理影院| 国内揄拍国产精品人妻在线| 麻豆成人午夜福利视频| 国产黄片美女视频| 久久影院123| 丁香六月天网| 亚洲美女黄色视频免费看| 亚洲综合精品二区| 国产精品国产av在线观看| 日韩精品有码人妻一区| 水蜜桃什么品种好| 大香蕉久久网| 大片电影免费在线观看免费| 91精品国产九色| 日本黄大片高清| 久久久久久久精品精品| 在线观看人妻少妇| 国产精品无大码| 能在线免费看毛片的网站| 欧美高清成人免费视频www| 视频区图区小说| 欧美日韩av久久| 国产在线视频一区二区| 一级片'在线观看视频| 亚洲国产成人一精品久久久| 国产免费一级a男人的天堂| 最近中文字幕高清免费大全6| 久久婷婷青草| 久久精品国产自在天天线| 一本久久精品| 国产又色又爽无遮挡免| 26uuu在线亚洲综合色| 久久久久网色| 岛国毛片在线播放| 国产亚洲欧美精品永久| 另类精品久久| 欧美日韩综合久久久久久| 国产精品欧美亚洲77777| 极品教师在线视频| 欧美xxⅹ黑人| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品国产av成人精品| 国产一区二区在线观看av| 亚洲,一卡二卡三卡| 成人美女网站在线观看视频| 久久久午夜欧美精品| 人妻人人澡人人爽人人| 人人妻人人爽人人添夜夜欢视频 | 最近手机中文字幕大全| 熟女av电影| 国产成人免费观看mmmm| 狂野欧美白嫩少妇大欣赏| 国产探花极品一区二区| 男女边吃奶边做爰视频| 极品少妇高潮喷水抽搐| 在线观看免费日韩欧美大片 | 国产精品一区www在线观看| 99国产精品免费福利视频| 国产一区亚洲一区在线观看| 丝袜在线中文字幕| 三级国产精品欧美在线观看| 国产黄片视频在线免费观看| 久久毛片免费看一区二区三区| 亚洲一区二区三区欧美精品| 最近手机中文字幕大全| 国产伦在线观看视频一区| 日本猛色少妇xxxxx猛交久久| 亚洲图色成人| 亚洲欧美精品自产自拍| 久久影院123| 日本欧美视频一区| 国产亚洲5aaaaa淫片| 免费观看性生交大片5| 中文字幕亚洲精品专区| 国产高清不卡午夜福利| 欧美国产精品一级二级三级 | 免费观看a级毛片全部| 卡戴珊不雅视频在线播放| 亚洲av成人精品一区久久| 午夜日本视频在线| 国产av国产精品国产| 国产午夜精品一二区理论片| 中文欧美无线码| 国模一区二区三区四区视频| 青春草视频在线免费观看| 亚洲美女黄色视频免费看| 亚洲欧美清纯卡通| 中国美白少妇内射xxxbb| 搡女人真爽免费视频火全软件| 亚洲av在线观看美女高潮| 色吧在线观看| a级一级毛片免费在线观看| 国产亚洲午夜精品一区二区久久| 午夜福利,免费看| 91在线精品国自产拍蜜月| 久久久国产精品麻豆| 午夜老司机福利剧场| 欧美三级亚洲精品| 最近的中文字幕免费完整| 国产黄色免费在线视频| 婷婷色av中文字幕| 国产一区二区在线观看av| av在线老鸭窝| 国国产精品蜜臀av免费| 亚洲在久久综合| 毛片一级片免费看久久久久| 精品熟女少妇av免费看| 国产美女午夜福利| 91aial.com中文字幕在线观看| 中文字幕亚洲精品专区| 97超视频在线观看视频| 另类亚洲欧美激情| 最近中文字幕高清免费大全6| 久久久久久久久久久免费av| 夜夜骑夜夜射夜夜干| 成人毛片a级毛片在线播放| 日韩一区二区三区影片| 日韩精品有码人妻一区| 日本欧美视频一区| 亚洲av电影在线观看一区二区三区| 亚洲性久久影院| 在线观看免费高清a一片| 久久国产精品大桥未久av | 欧美精品一区二区免费开放| 久久精品国产亚洲av天美| 国产欧美另类精品又又久久亚洲欧美| 亚洲人成网站在线观看播放| 亚洲色图综合在线观看| 欧美日韩一区二区视频在线观看视频在线| 婷婷色综合www| 欧美日韩综合久久久久久| 毛片一级片免费看久久久久| 99热这里只有是精品50| 日韩av在线免费看完整版不卡| 一本大道久久a久久精品| 黄色一级大片看看| 久久ye,这里只有精品| 精品亚洲乱码少妇综合久久| 亚洲av综合色区一区| 欧美精品人与动牲交sv欧美| 国产精品久久久久成人av| 久久久久久久久久久免费av| 国产黄片视频在线免费观看| 欧美+日韩+精品| 一级二级三级毛片免费看| 人人妻人人看人人澡| 一级毛片电影观看| 少妇被粗大的猛进出69影院 | 国产真实伦视频高清在线观看| 老熟女久久久| 男女边吃奶边做爰视频| 亚洲欧美日韩东京热| 一区二区三区乱码不卡18| 国内少妇人妻偷人精品xxx网站| 超碰97精品在线观看| 亚洲av成人精品一二三区| 国产免费又黄又爽又色| 青青草视频在线视频观看| 黄色日韩在线| av女优亚洲男人天堂| 亚洲美女黄色视频免费看| 亚洲国产精品国产精品| 欧美少妇被猛烈插入视频| 久久鲁丝午夜福利片| 国语对白做爰xxxⅹ性视频网站| 在线观看av片永久免费下载| 在线精品无人区一区二区三| 精品国产国语对白av| 成人毛片60女人毛片免费| 热re99久久精品国产66热6| 菩萨蛮人人尽说江南好唐韦庄| 成人美女网站在线观看视频| 在线观看美女被高潮喷水网站| 久久国产精品大桥未久av | av福利片在线| 精品一区二区三卡| 最近中文字幕高清免费大全6| 夜夜看夜夜爽夜夜摸| 免费在线观看成人毛片| 精品一区二区免费观看| 成人午夜精彩视频在线观看| 在线观看免费高清a一片| 永久免费av网站大全| 嘟嘟电影网在线观看| 色吧在线观看| 久久av网站| 性色avwww在线观看| 亚洲av国产av综合av卡| 乱人伦中国视频| 国产精品欧美亚洲77777| 免费大片黄手机在线观看| 老司机亚洲免费影院| 国产真实伦视频高清在线观看| 在线观看人妻少妇| 9色porny在线观看| 人体艺术视频欧美日本| 99热这里只有是精品在线观看| 亚洲伊人久久精品综合| 18+在线观看网站| 精品卡一卡二卡四卡免费| 男女啪啪激烈高潮av片| 国产一级毛片在线| 国产亚洲91精品色在线| 亚洲国产日韩一区二区| 精品一区二区三区视频在线| 69精品国产乱码久久久| 国产淫语在线视频| 国产成人免费无遮挡视频| 欧美精品一区二区大全| 欧美 日韩 精品 国产| 国产男人的电影天堂91| 人妻制服诱惑在线中文字幕| 国产淫片久久久久久久久| 美女主播在线视频| h视频一区二区三区| 91久久精品国产一区二区成人| 亚洲国产精品一区二区三区在线| 成人二区视频| 久久久国产一区二区| 伊人久久精品亚洲午夜| 精品一区二区三区视频在线| 日韩伦理黄色片| 嫩草影院入口| 亚洲av国产av综合av卡| 中文字幕人妻熟人妻熟丝袜美| 青春草亚洲视频在线观看| 在线天堂最新版资源| 啦啦啦啦在线视频资源| 国产精品人妻久久久影院| 一个人看视频在线观看www免费| 在现免费观看毛片| 日韩伦理黄色片| 两个人免费观看高清视频 | 久久久亚洲精品成人影院| 黑人巨大精品欧美一区二区蜜桃 | 日韩中文字幕视频在线看片| 不卡视频在线观看欧美| 简卡轻食公司| 草草在线视频免费看| 亚洲精华国产精华液的使用体验| 精品久久久精品久久久| 只有这里有精品99| 久久久久久久亚洲中文字幕| 久久青草综合色| av网站免费在线观看视频| 亚洲图色成人| 欧美性感艳星| 天堂中文最新版在线下载| 欧美日韩亚洲高清精品| 夫妻午夜视频| 最近中文字幕2019免费版| 国产中年淑女户外野战色| 国产成人91sexporn| 亚洲无线观看免费| 老司机影院成人| 国产极品粉嫩免费观看在线 | 成人免费观看视频高清| 少妇熟女欧美另类| 日日撸夜夜添| 日本黄色日本黄色录像| 亚洲欧美一区二区三区国产| 女人精品久久久久毛片| 久久av网站| 男人爽女人下面视频在线观看| 亚洲欧美成人综合另类久久久| 97超碰精品成人国产| 我要看黄色一级片免费的| 欧美日韩视频精品一区| 久久婷婷青草| 一级毛片我不卡| 欧美日韩在线观看h| 中文字幕人妻熟人妻熟丝袜美| 色视频在线一区二区三区| 伊人亚洲综合成人网| av在线app专区| 日本爱情动作片www.在线观看| 在线观看av片永久免费下载| 涩涩av久久男人的天堂| 啦啦啦中文免费视频观看日本| 男女边摸边吃奶| 观看免费一级毛片| 高清不卡的av网站| 亚洲精品亚洲一区二区| 久久久久久久亚洲中文字幕| 91精品国产国语对白视频| 亚洲精品,欧美精品| 免费人成在线观看视频色| 嫩草影院入口| 国产真实伦视频高清在线观看| 国产伦精品一区二区三区四那| 午夜福利影视在线免费观看| 看非洲黑人一级黄片| 91精品一卡2卡3卡4卡| 国产淫片久久久久久久久| 日日撸夜夜添|