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

    基于振幅方程組的行波對流的數(shù)值模擬

    2014-03-26 08:20:12王卓運寧利中王娜李開繼
    西安理工大學(xué)學(xué)報 2014年2期
    關(guān)鍵詞:群速度瑞利臨界點

    王卓運,寧利中,王娜,李開繼

    (西安理工大學(xué) 陜西省西北旱區(qū)生態(tài)水利工程國家重點實驗室,陜西 西安 710048)

    對流現(xiàn)象不僅在日常生活中很常見,同時也是自然界中很普遍的現(xiàn)象,比如:引起天氣、氣候變化的大氣環(huán)流、海洋環(huán)流,引起大陸板塊漂移的地幔對流,決定地磁場產(chǎn)生和變化的外地核對流等[1]。這些對流現(xiàn)象都涉及到了Rayleigh-Benard對流系統(tǒng)。Rayleigh-Benard對流則是研究熱對流現(xiàn)象的一個經(jīng)典模型,是指在一個封閉的腔體內(nèi),其上表面溫度恒定,下表面加熱,從而形成溫度差導(dǎo)致腔體內(nèi)流體運動的流動現(xiàn)象。由于該系統(tǒng)試驗簡單且易于控制,以及描述該系統(tǒng)的精確方程已知,便于理論分析,故該模型已被普遍接受,是用來研究對流穩(wěn)定性、時空結(jié)構(gòu)和非線性特性的模型之一。

    目前,對于Rayleigh-Benard對流問題研究手段主要有實驗研究、模型方程探討、流體力學(xué)方程數(shù)值模擬等方法[2-9]。王濤、葛永斌等人對腔體內(nèi)的對流進行了高精度數(shù)值模擬,對雙局部進行波的研究進一步進行了擴充,并研究了具有較弱Soret效應(yīng)下,附加一個微小的溫度擾動作為擾動源的中等長高比腔體內(nèi)混合流體對流系統(tǒng)中時空結(jié)構(gòu)的發(fā)展[10-12]。郝建武在長高比Γ=40、分離比ψ=-0.6的條件下,得到了一種新的現(xiàn)象,即單側(cè)缺陷擺動對傳行波,并對其形成機理與特性作了進一步探討[13]。在弱非線性的假定下,在分歧點附近通過級數(shù)展開等方法,人們已建立了各種模型方程或振幅方程。這些方程包括:Ginzbarg-Landau(GL)方程式、復(fù)數(shù)GL方程式、耦合的GL方程式、Kuramoto-Sivashinsky方程式、Swift-Hohenberg方程式等[14-21]。

    目前,基于對流振幅來表征行波對流斑圖的研究報道還很少,本文利用Ginzbarg-Landau(GL)方程組,通過研究對流的振幅變化來探討行波對流的斑圖特性及非平衡體對流的非線性動力特性。

    1 基本方程組

    1.1 Ginzbarg-Landau(GL)方程組

    本文研究行波對流的斑圖時選取如下GL方程組[15-16]:

    (1)

    (2)

    式中A、B為行波對流的振幅,τ0、ξ0、s分別反映對流運動的特性時間、長度、行進波的群速度;g1、g2、k反映非線性系數(shù)。

    定義參數(shù):

    整理式(1)、(2),然后去掉方程(1)、(2)參數(shù)的“′”,則GL方程組可簡化為:

    (3)

    (4)

    方程(3)、(4)中,v反映行波群速度,ε反映相對瑞利數(shù),g代表非線性系數(shù)。

    兩端采用反射邊界條件,方程形式為[16]:

    (5)

    (6)

    其中 ,α、β為反射系數(shù),計算中取α=-0.06,β=-0.06。

    本文在計算時給出左行波振幅A和右行波振幅B的初值分別為:

    式中L代表腔體長度。

    1.2 振幅方程組的差分式

    GL方程組簡化式經(jīng)過一階差分和二階差分,式子為:

    (7)

    (8)

    左邊界方程的差分式為:

    (9)

    (10)

    右邊界方程的差分式為:

    (11)

    (12)

    2 結(jié)果與分析

    通過對有限腔體長度L=20、L=40、L=60內(nèi)的對流進行模擬,分析了在不同腔體長度下,行波振幅隨相對瑞利數(shù)和群速度的變化,發(fā)現(xiàn)了行波振幅存在三種狀態(tài):TI(Time independent),即行波的振幅不隨時間變化;TD(Time dependent),即行波的振幅對時間有依賴性,隨時間發(fā)生變化;Null,即行波的振幅為零。

    將振幅狀態(tài)與實際流體斑圖狀態(tài)一一對應(yīng),發(fā)現(xiàn)存在行波(Traveling wave(TW))、局部行波(Localized traveling wave (LTW))、對傳波(Counterpropagating waves (CPW)) 、Blinking traveling wave(BTW)斑圖。并在ε-v平面中根據(jù)行波對流振幅的變化情況進行了分區(qū)。

    2.1 行波對流的斑圖

    2.1.1 振幅狀態(tài)為TI的對傳波CPW

    對于在腔體總長度為60個單位長度,群速度v=1.0,相對瑞利數(shù)ε=0.29的情況,得到了振幅狀態(tài)為TI的對傳波CPW,如圖1所示。通過對圖1的分析,得到以下特性。

    經(jīng)過250個單位時間左右行波對流狀態(tài)都達到穩(wěn)定。在長度為60個單位長度的區(qū)間內(nèi),左行波從腔體的最左端即L=0處到L=20處附近振幅為0.5,然后在L=20處振幅開始衰減,到L=35附近處開始衰減為零,一直延續(xù)到L=50附近振幅仍為零,而后一直到腔體的右邊壁L=60處行波振幅有小幅度的回升,其數(shù)量級是10-2級,大小可以忽略,在整個長度區(qū)間左行波保持這種狀態(tài)不隨時間變化而持續(xù)前進;右行波從腔體的最右端L=60處開始一直到L=45處附近行波振幅為0.5,隨后振幅開始出現(xiàn)連續(xù)性衰減,直到L=30個單位長度附近振幅為零,這種狀態(tài)一直延續(xù)發(fā)展到L=10處附近振幅仍為零,而后行波振幅出現(xiàn)“抬頭”并保持增長直到腔體的最左端L=0處,其數(shù)量級也是10-2級,大小可以忽略,右行波在整個腔體長度區(qū)間內(nèi)保持這種狀態(tài)且不隨時間發(fā)生變化而發(fā)展下去。

    圖1 L=60, v=1.0, ε=0.29左右行波振幅

    左右行波的振幅狀態(tài)達到穩(wěn)定后都不隨時間發(fā)生變化即振幅狀態(tài)為TI狀態(tài),此種波在實際中對應(yīng)對傳波。

    2.1.2 振幅狀態(tài)為TD的BTW

    腔體長度L=60、群速度v=1.5、相對瑞利數(shù)ε=0.2的振幅情況如圖2所示。通過對圖2的分析,得到以下特性。

    左右行波的對流狀態(tài)在達到穩(wěn)定后,其振幅均表現(xiàn)出了“雜亂”性,且振幅大小隨時間發(fā)生變化,即振幅狀態(tài)都為TD。

    左行波振幅在腔體的左端開始發(fā)展,在L=10處附近振幅達到最大值0.4,而后出現(xiàn)連續(xù)性衰減直到腔體的最右端為零,隨著時間的發(fā)展,出現(xiàn)最大振幅區(qū)域的長度與整個腔體長度的比例時而變小,如250個單位時間到300個單位時間段;時而在整個腔體中振幅為零,如350個單位時間到390個單位時間段;時而變大,如425個單位時間到450個單位時間段。右行波也展現(xiàn)出了類似的特性。

    此種波在實際中對應(yīng)Blinking traveling wave(BTW)。

    圖2 L=60, v=1.5, ε=0.2左右行波振幅

    2.1.3 振幅狀態(tài)為Null的行波

    腔體長度L=60、群速度v=1.0、相對瑞利數(shù)ε=0.0的振幅情況如圖3所示。通過對圖3的分析,得到以下特性。

    圖3 L=60, v=1.0, ε=0.0左右行波振幅

    經(jīng)過300個單位時間對流達到穩(wěn)定。

    達到穩(wěn)定狀態(tài)后左行波與右行波振幅都為零,且都不隨時間發(fā)生變化,即振幅狀態(tài)為Null。

    左右行波振幅都為零,且不隨時間發(fā)生變化即振幅狀態(tài)為TI。

    2.1.4 振幅狀態(tài)為TI的行波TW

    腔體長度L=40、群速度v=0.2、相對瑞利數(shù)ε=0.6的振幅情況如圖4所示。通過對圖4的分析,得到以下特性。

    圖4 L=40, v=0.2, ε=0.6左右行波振幅

    經(jīng)過250個單位時間,左右行波對流狀態(tài)都達到穩(wěn)定。

    在長度為40個單位長度的區(qū)間內(nèi),左行波從腔體內(nèi)部L=5處到L=35處附近振幅為零,但在腔體的兩端L=0到L=5附近及L=35到L=40附近,振幅出現(xiàn)“抬頭”并保持增長,其數(shù)量級是10-2級,大小可以忽略。

    右行波在整個腔體內(nèi)振幅大小恒為0.6,保持這種狀態(tài)且不隨時間發(fā)生變化地發(fā)展下去。

    左右行波的振幅狀態(tài)達到穩(wěn)定后都不隨時間發(fā)生變化即振幅狀態(tài)為TI狀態(tài),此種波在實際中對應(yīng)行波。

    2.1.5 振幅狀態(tài)為TD且具有一定周期性的行波

    腔體長度L=40、群速度v=1.5、相對瑞利數(shù)ε=0.3的振幅情況如圖5所示。通過對圖5的分析,得到以下特性。

    圖5 L=40, v=1.5, ε=0.3左右行波振幅

    左右行波振幅狀態(tài)表現(xiàn)為TD,但區(qū)別于BTW,即出現(xiàn)了一定的周期性和對稱性。左右行波的變化周期均為100個單位時間左右。

    2.1.6 振幅狀態(tài)為TI的局部行波LTW

    腔體長度L=20、群速度v=0.5、相對瑞利數(shù)ε=0.2的振幅情況如圖6所示。通過對圖6的分析,得到以下特性。

    圖6 L=20, v=0.5, ε=0.2左右行波振幅

    右行波在腔體兩端同時連續(xù)增長,最終在腔體中間L=4到L=12區(qū)域振幅保持最大值0.4;左行波振幅在整個腔體內(nèi)一直為零,但在腔體的左端有輕微“抬頭”,數(shù)值可以忽略。腔體中僅僅留下右行波。

    左右行波振幅狀態(tài)表現(xiàn)為TI狀態(tài),對應(yīng)局部行波。

    2.2 振幅狀態(tài)分區(qū)圖

    2.2.1L=60的行波振幅狀態(tài)分區(qū)

    在長度L=60的腔體內(nèi),筆者研究了行波振幅隨群速度v在0.0~2.0、相對瑞利數(shù)ε在-0.20~0.60的范圍內(nèi)的分區(qū)情況,如圖7所示。

    圖7 L=60的行波振幅狀態(tài)分區(qū)

    群速度v=0時: 臨界狀態(tài)出現(xiàn)在相對瑞利數(shù)ε=-0.01處,此時當(dāng)相對瑞利數(shù)ε再增大0.01即ε=0.00時行波振幅隨時間的變化情況就發(fā)生明顯變化:-0.20≤ε≤-0.01時,左右行波穩(wěn)定后,行波振幅都變?yōu)榱?,即處于Null區(qū)域;0.00≤ε≤0.60時,左右行波穩(wěn)定后, 左右行波振幅不隨時間變化,即處于TI區(qū)域。同理,在v=0.25時臨界點ε=0.00;在v=0.5時臨界點ε=0.02。

    群速度v=0.75:-0.20≤ε≤0.03時,行波振幅處于Null區(qū)域;0.14≤ε≤0.15時,處于TD區(qū)域;0.16≤ε≤0.60時,行波振幅處于TI區(qū)域。同理,在v=1.0時有兩個臨界點ε=0.04和ε=0.28;在v=1.25時有臨界點ε=0.05和ε=0.42。

    群速度v=1.5:-0.20≤ε≤0.05時,行波振幅處于Null區(qū)域;0.06≤ε≤0.60時,行波振幅處于TD區(qū)域。同理,在v=1.75時臨界點是ε=0.07;在v=2.0時臨界點是ε=0.07。

    在群速度v=0.0到v=0.55,行波振幅狀態(tài)劃分為Null與TI兩個區(qū)域; 在群速度v=0.55到v=2.0,行波振幅狀態(tài)劃分為Null、TD及TI三個區(qū)域。

    2.2.2L=40的行波振幅狀態(tài)分區(qū)

    L=40的行波振幅分區(qū)圖如圖8所示。v=0、 -0.20≤ε≤ -0.02時,行波振幅處于Null區(qū)域;v=0,-0.01≤ε≤0.60時,行波振幅處于TI區(qū)域。同理,在v=0.25時臨界點ε=0.01,在v=0.5時臨界點ε=0.03。

    圖8 L=40的行波振幅狀態(tài)分區(qū)

    群速度v=0.75:-0.20≤ε≤ 0.06時,行波振幅處于Null區(qū)域;0.07≤ε≤ 0.13時,即處于TD區(qū)域;0.14≤ε≤ 0.60時,行波振幅處于TI區(qū)域。同理,在v=1.0時有兩個臨界點ε=0.06和ε=0.27;在v=1.25時有臨界點ε=0.06和ε=0.43。

    群速度v=1.5:-0.20≤ε≤ 0.08時,行波振幅處于Null區(qū)域;0.09≤ε≤ 0.60時,行波振幅處于TD區(qū)域。同理,在v=1.75時臨界點ε=0.09;在v=2.0時臨界點ε=0.10。

    在群速度v=0.0到v=0.75,行波振幅狀態(tài)劃分為Null與TI兩個區(qū)域,在群速度v=0.75到v=2.0,行波振幅狀態(tài)劃分為Null、TD及TI三個區(qū)域。

    2.2.3L=20的行波振幅狀態(tài)分區(qū)

    在長度L=20的腔體內(nèi),左右行波振幅的劃分區(qū)域如圖9所示。

    圖9 L=20的行波振幅狀態(tài)分區(qū)

    群速度v=0:相對瑞利數(shù)-0.20≤ε≤ 0.00時,處于Null區(qū)域;0.01≤ε≤ 0.60時,處于TI區(qū)域。同理,在v=0.25時有臨界點ε=0.02;在v=0.5時有臨界點ε=0.06;v=0.75時有臨界點ε=0.10。

    群速度v=1.0:-0.20≤ε≤0.13時,行波振幅處于Null區(qū)域;0.14≤ε≤0.32時,處于TD區(qū)域;0.33≤ε≤0.60時,行波振幅處于TI區(qū)域。同理,在v=1.25時有臨界點ε=0.15和ε=0.50。

    群速度v=1.5:相對瑞利數(shù)-0.20≤ε≤0.17時,處于Null區(qū)域;0.18≤ε≤0.60時,行波振幅處于TD區(qū)域。同理,在v=1.75時臨界點ε=0.20;在v=2.0時臨界點ε=0.22。

    在群速度v=0.0到v=0.8,行波振幅狀態(tài)劃分為Null與TI兩個區(qū)域; 在群速度v=0.8到v=2.0,行波振幅狀態(tài)劃分為Null、TD及TI三個區(qū)域。

    3 結(jié) 語

    在腔體長度一定的情況下影響振幅狀態(tài)的參數(shù)主要有兩個:群速度與相對瑞利數(shù)。

    隨著群速度v與相對瑞利數(shù)ε變化,腔體長度L=20、L=40、L=60的振幅都被劃分為“TI、TD、Null”三個區(qū)域,即時間獨立區(qū)、時間依賴區(qū)、零振幅區(qū)。

    在TI區(qū),行波振幅不依賴于時間,即行波振幅達到穩(wěn)定狀態(tài)后不隨時間而改變。它包括兩種情況:①均勻行波,即左行波振幅穩(wěn)定且不為零,右行波振幅穩(wěn)定且為零,或者右行波振幅穩(wěn)定且不為零,左行波振幅穩(wěn)定且為零;②對傳波,即左右行波振幅同時穩(wěn)定存在且都不為零。

    在TD區(qū),行波振幅依賴于時間,即行波振幅達到穩(wěn)定狀態(tài)后隨時間而變化,它包括兩種情況:①BTW,即左右行波振幅都隨時間在發(fā)生變化且雜亂無章;②空間周期性,即左右行波振幅隨時間在發(fā)生變化但出現(xiàn)一定的對稱性,具有周期性。

    腔體長度L=20、L=40、L=60,在TI區(qū)域振幅都是經(jīng)過了250個單位時間達到穩(wěn)定;在Null區(qū)域振幅都是經(jīng)過了300個時間單位達到穩(wěn)定。

    參考文獻:

    [1] 余荔,寧利中,魏炳乾,等. Rayleigh-Benard對流及其在工程中的應(yīng)用[J]. 水資源與水工程學(xué)報,2008,19(3):52-54.

    Yu Li, Ning Lizhong, Wei Binqian, et al. Rayleigh-Benard convection and application in engineering[J]. Journal o f Water Resources & Water Engineering, 2008,19(3):52-54.

    [2] Ning Lizhong,Harada Y, Yahata H. Localized traveling waves in binary fluid convection[J]. Progress of Theoretical Physics,1996,96(4):669-682.

    [3] Ning Lizhong, Harada Y, Yahata H. Dynamics of localized traveling wave in binary fluid mixtures [J]. Journal of Hydrodynamics, 1998,10(2):29-39.

    [4] Ning Lizhong, Harada Y, Yahata H. Formation process of the traveling wave state with a defect in binary fluid convection [J]. Progress of Theoretical Physics,1997,98(3):551-566.

    [5] Ning Lizhong, Harada Y, Yahata H. Modulated traveling waves in binary fluid convection in an intermediate-aspect-ratio rectangular [J]. Progress of Theoretical Physics,1997,97(6):831-848.

    [6] 寧利中,余荔, 袁喆, 等. 沿混合流體對流分叉曲線上部分支行波斑圖的演化[J].中國科學(xué)G,2009,39(5):746-751.

    Ning Lizhong, Yu Li, Yuan Zhe, et al. Evolution of traveling wave patterns along upper branch of bifurcation diagram in binary fluid convection [J]. Science in China G, 2009,39(5): 746-751.

    [7] 寧利中,齊昕, 周洋,等.混合流體行波對流中的缺陷結(jié)構(gòu)[J].物理學(xué)報,2009,58(4):2528-2534.

    Ning Lizhong, Qi Xin, Zhou Yang, et al. Defect structures of Rayleigh-Benard traveling wave convection in binary fluid mixtures [J]. Acta Phys Sin-cn ED, 2009,58(4):2528-2534.

    [8] 趙秉新, 田振夫.底部加熱平面Poiseuille 流中的局部行波結(jié)構(gòu)[J].水動力學(xué)研究與進展:A 輯,2012,27(6):649-658.

    Zhao Bingxin, Tian Zhenfu.Localized traveling wave convection in plan Poiseuille flow heated from below[J].Chinese Journal of Hydrodynamics, 2012,27(6):649-658.

    [9] 齊昕,寧利中,劉嘉夫,等. 中等長高比腔體內(nèi)的混合流體Undulation行進波[J].西安理工大學(xué)學(xué)報,2010,26(3):271-276.

    Qi Xin,Ning Lizhong,Liu Jiafu,et al.Undulation travelling wave convection in binary fluid mixtures in an intermediate-aspect-ratio cell[J].Journal of Xi’an University of Technology,2010,26(3):271-276.

    [10] 王濤,田振夫,葛永斌.長腔體內(nèi)混合流體行進波對流的高精度數(shù)值模擬[J].水動力學(xué)研究與進展:A 輯,2011,26(1):41-47.

    Wang Tao,Tian Zhenfu,Ge Yongbin.High accuracy numerical simulation of travelling wave convection in binary fluid mixture in large aspect ratio rectangular cell[J].Chinese Journal of Hydrodynamics,2011,26(1):41-47.

    [11] 王濤,葛永斌.雙局部行進波對流的高精度數(shù)值模擬[J].合肥工業(yè)大學(xué)學(xué)報:自然科學(xué)版,2012,35(6):842-847.

    Wang Tao,Ge Yongbin.High accuracy numerical simulation of doubly localized travelling wave convection [J].Journal of Hefei University of Technology(Natural Science),2012,35(6):842-847.

    [12] 王濤,葛永斌.微小擾動下中等長高比腔體內(nèi)行進波對流的高精度數(shù)值模擬[J].合肥工業(yè)大學(xué)學(xué)報:自然科學(xué)版,2012,35(11):1569-1573.

    Wang Tao,Ge Yongbin.High accuracy numerical simulation of travelling wave convection in an intermediate-aspect-ratio rectangular cell with small perturbation[J].Journal of Hefei University of Technology(Natural Science),2012,35(11):1569-1573.

    [13] 郝建武,寧利中,王卓運,等.長進行截面腔體內(nèi)具有缺陷的對傳行波斑圖[J].力學(xué)季刊,2013,34(1):139-146.

    Hao Jianwu,Ning Lizhong,Wang Zhuoyun,et al.Patterns of counter propagating wave with defects in a long rectangular cell[J].Chinese Quarterly of Mechanics,2013,34(1):139-146.

    [14] Aranson I S,Kramer L. The world of the complex Ginzburg-Landau equation[J].Review of Modern Physics, 2002,74(1):99-143.

    [15] Sullivan T S,Deissler R J. Elimination of hysteresis in a system of coupled Ginzburg-Landau equations [J]. Physical Review A,1989,40(11):6748-6751.

    [16] Cross M C. Structure of nonlinear traveling-wave states in finite geometries[J].Physical Review A,1988,38(7):3593-3600.

    [17] Cross M,Hohenberg P C.Pattern formation outside of equilibrium[J]. Reviews of Modern Physics,1993, 65 (3): 855-1112.

    [18] RIECKE. Self-trapping of traveling-wave pulses in binary mixture convection[J]. Physical Review Letter,1992,(68):301-304.

    [19] Cross M, Greenside H.Pattern formation and dynamics in nonequilibrium systems[M]. UK:Cambridge University Press,2009.

    [20] 蔵本由記,川崎恭治,山田道夫,等.パターン形成[M].東京:朝倉書店,1991:1-146.

    [21] 森肇,蔵本由記.散逸構(gòu)造とカオス[M].東京:巖波書店,2000:1-46.

    猜你喜歡
    群速度瑞利臨界點
    亞瑞利散斑場的二階累積量鬼成像
    基于臨界點的杭州灣水體富營養(yǎng)化多年變化研究
    海洋通報(2022年4期)2022-10-10 07:41:48
    VTI介質(zhì)地震波群速度三維擴展各向異性線性近似表征
    激光調(diào)制摻鉺光纖中光速及其在高靈敏光纖傳感領(lǐng)域的潛在應(yīng)用研究
    物質(zhì)波的波速與頻率公式中的能量
    小議超光速
    超越生命的臨界點
    超越生命的臨界點
    馬瑞利推出多項汽車零部件技術(shù)
    汽車零部件(2015年4期)2015-12-22 05:32:56
    瑞利波頻散成像方法的實現(xiàn)及成像效果對比研究
    我的女老师完整版在线观看| 天堂动漫精品| 成人特级黄色片久久久久久久| 全区人妻精品视频| 国产综合懂色| 精品99又大又爽又粗少妇毛片| 成人鲁丝片一二三区免费| 99热只有精品国产| 国产精品,欧美在线| 看黄色毛片网站| 校园春色视频在线观看| 欧美成人a在线观看| 狠狠狠狠99中文字幕| 日韩精品有码人妻一区| 国产三级在线视频| 国产一区二区三区av在线 | 六月丁香七月| 国产精品一区二区三区四区久久| 黄片wwwwww| 国产不卡一卡二| 美女大奶头视频| 中文字幕av成人在线电影| 日本五十路高清| 两个人的视频大全免费| 国产精品乱码一区二三区的特点| 熟妇人妻久久中文字幕3abv| 一个人看视频在线观看www免费| 国产精品电影一区二区三区| 国产欧美日韩精品一区二区| 大又大粗又爽又黄少妇毛片口| 欧美日韩国产亚洲二区| 日日摸夜夜添夜夜添av毛片| 91狼人影院| 黄色配什么色好看| 久久鲁丝午夜福利片| 国产精品亚洲美女久久久| 大香蕉久久网| 一个人看的www免费观看视频| 亚洲人成网站在线播放欧美日韩| 在线a可以看的网站| 国产熟女欧美一区二区| 99视频精品全部免费 在线| 婷婷亚洲欧美| av.在线天堂| 亚洲精品日韩av片在线观看| 嫩草影视91久久| 男人舔女人下体高潮全视频| 成人国产麻豆网| 久久久久久久久大av| 我的女老师完整版在线观看| 国产高清视频在线观看网站| 亚洲av一区综合| 美女cb高潮喷水在线观看| 99热全是精品| 秋霞在线观看毛片| 免费在线观看成人毛片| 别揉我奶头~嗯~啊~动态视频| 在现免费观看毛片| 国产中年淑女户外野战色| 国产精品美女特级片免费视频播放器| 欧美成人免费av一区二区三区| 日本一本二区三区精品| 成人欧美大片| 黄色一级大片看看| 熟女电影av网| 欧美日本亚洲视频在线播放| 午夜a级毛片| 国产高清视频在线观看网站| 一区二区三区免费毛片| 熟女人妻精品中文字幕| 麻豆av噜噜一区二区三区| 欧美另类亚洲清纯唯美| 精品久久久噜噜| 成人欧美大片| 变态另类成人亚洲欧美熟女| 成年女人看的毛片在线观看| 免费在线观看影片大全网站| 国产男靠女视频免费网站| 日韩精品青青久久久久久| 国产高清激情床上av| 麻豆国产97在线/欧美| 露出奶头的视频| 校园人妻丝袜中文字幕| 午夜爱爱视频在线播放| 可以在线观看的亚洲视频| 成年av动漫网址| 国内精品久久久久精免费| 成人鲁丝片一二三区免费| 精品久久久久久久久久久久久| 精品久久久久久久人妻蜜臀av| 少妇的逼水好多| 国产欧美日韩精品一区二区| 亚洲经典国产精华液单| 人妻夜夜爽99麻豆av| 日韩中字成人| 日韩强制内射视频| 禁无遮挡网站| 如何舔出高潮| 亚洲成av人片在线播放无| 波多野结衣高清无吗| 五月玫瑰六月丁香| 我要搜黄色片| 欧美xxxx黑人xx丫x性爽| 乱系列少妇在线播放| 国产精品美女特级片免费视频播放器| 最近2019中文字幕mv第一页| 最近最新中文字幕大全电影3| 国产 一区 欧美 日韩| 久久久久精品国产欧美久久久| 亚洲丝袜综合中文字幕| 日韩欧美精品v在线| 日本五十路高清| 亚洲熟妇中文字幕五十中出| av卡一久久| 精品人妻熟女av久视频| 在线播放国产精品三级| 国产成人freesex在线 | 久久综合国产亚洲精品| 性色avwww在线观看| 国产69精品久久久久777片| 美女高潮的动态| 男女那种视频在线观看| 亚洲国产日韩欧美精品在线观看| 啦啦啦观看免费观看视频高清| 欧美不卡视频在线免费观看| 一区二区三区高清视频在线| 国产精品久久久久久精品电影| 成年女人永久免费观看视频| 国产片特级美女逼逼视频| 真实男女啪啪啪动态图| 18+在线观看网站| 小说图片视频综合网站| 亚洲美女黄片视频| 亚洲精品色激情综合| 超碰av人人做人人爽久久| 可以在线观看的亚洲视频| 99热这里只有是精品在线观看| 免费黄网站久久成人精品| 国产 一区 欧美 日韩| 国产女主播在线喷水免费视频网站 | 波多野结衣高清作品| 欧美xxxx性猛交bbbb| 熟妇人妻久久中文字幕3abv| 午夜a级毛片| 国产亚洲精品久久久com| 国产在线男女| 九九久久精品国产亚洲av麻豆| 国产亚洲精品久久久com| 波多野结衣巨乳人妻| 97人妻精品一区二区三区麻豆| 夜夜爽天天搞| 赤兔流量卡办理| 97人妻精品一区二区三区麻豆| 一区二区三区免费毛片| 在线国产一区二区在线| 少妇丰满av| 俺也久久电影网| 天堂av国产一区二区熟女人妻| 2021天堂中文幕一二区在线观| av卡一久久| 青春草视频在线免费观看| 亚洲第一区二区三区不卡| 婷婷六月久久综合丁香| 高清日韩中文字幕在线| 亚洲国产欧美人成| 嫩草影院新地址| 搡老妇女老女人老熟妇| 日本一二三区视频观看| 国产v大片淫在线免费观看| 欧美最新免费一区二区三区| АⅤ资源中文在线天堂| 老司机福利观看| 国产一级毛片七仙女欲春2| 床上黄色一级片| 久久精品国产99精品国产亚洲性色| 日本成人三级电影网站| 亚洲最大成人av| 精品99又大又爽又粗少妇毛片| 三级国产精品欧美在线观看| 日日摸夜夜添夜夜爱| 成人美女网站在线观看视频| 午夜a级毛片| av女优亚洲男人天堂| 国产成人一区二区在线| 99热只有精品国产| 99久国产av精品| 午夜福利在线观看免费完整高清在 | 日韩三级伦理在线观看| 亚洲精品456在线播放app| 精品不卡国产一区二区三区| 乱人视频在线观看| 乱码一卡2卡4卡精品| 国产三级中文精品| 欧美xxxx黑人xx丫x性爽| 国产在线男女| 亚洲精品色激情综合| 欧美丝袜亚洲另类| 国产精品1区2区在线观看.| 女生性感内裤真人,穿戴方法视频| 亚洲欧美精品综合久久99| 亚洲第一区二区三区不卡| 俄罗斯特黄特色一大片| 亚洲在线自拍视频| 国产 一区精品| 男人的好看免费观看在线视频| 91在线观看av| 国产精品一区二区三区四区久久| 九九爱精品视频在线观看| 成人av一区二区三区在线看| 欧美xxxx黑人xx丫x性爽| 少妇的逼水好多| 九九久久精品国产亚洲av麻豆| 最近2019中文字幕mv第一页| 国产精品亚洲美女久久久| 美女大奶头视频| a级毛片免费高清观看在线播放| 欧美zozozo另类| 非洲黑人性xxxx精品又粗又长| 亚洲av电影不卡..在线观看| 中文字幕av在线有码专区| 大型黄色视频在线免费观看| 日韩欧美精品免费久久| 国产精品亚洲美女久久久| 日韩三级伦理在线观看| 校园春色视频在线观看| 人人妻人人澡欧美一区二区| 久久九九热精品免费| 国产探花极品一区二区| 欧美一区二区国产精品久久精品| 久久人妻av系列| av天堂在线播放| 一区二区三区四区激情视频 | 俺也久久电影网| 精品99又大又爽又粗少妇毛片| 日韩成人伦理影院| 男女做爰动态图高潮gif福利片| 亚洲成人久久爱视频| 婷婷精品国产亚洲av| 久久久久久大精品| 午夜精品一区二区三区免费看| 日韩欧美精品v在线| 最近在线观看免费完整版| 人妻制服诱惑在线中文字幕| 国产一区二区在线观看日韩| 美女 人体艺术 gogo| .国产精品久久| 亚洲精品国产成人久久av| 亚洲天堂国产精品一区在线| 精品国产三级普通话版| 成人二区视频| 亚洲精品亚洲一区二区| 久久久色成人| 国产探花极品一区二区| 搡女人真爽免费视频火全软件 | 国产色爽女视频免费观看| 尾随美女入室| 婷婷精品国产亚洲av在线| 亚洲av第一区精品v没综合| 国产91av在线免费观看| 亚洲成人av在线免费| 十八禁网站免费在线| 久99久视频精品免费| av黄色大香蕉| 国产一级毛片七仙女欲春2| 亚洲无线在线观看| 成人av在线播放网站| 欧美最黄视频在线播放免费| 国内少妇人妻偷人精品xxx网站| 久久综合国产亚洲精品| 国产伦精品一区二区三区四那| 亚洲美女视频黄频| 蜜桃亚洲精品一区二区三区| 亚洲精品一区av在线观看| 亚洲四区av| 少妇高潮的动态图| 亚洲av第一区精品v没综合| 免费一级毛片在线播放高清视频| 精华霜和精华液先用哪个| 一个人观看的视频www高清免费观看| 亚洲一区高清亚洲精品| 精品人妻一区二区三区麻豆 | 亚洲成人精品中文字幕电影| 国产真实乱freesex| 亚洲av二区三区四区| 露出奶头的视频| 亚洲精品乱码久久久v下载方式| 日本五十路高清| 亚洲国产精品成人久久小说 | 深夜a级毛片| 精华霜和精华液先用哪个| 亚洲国产欧美人成| 亚洲美女搞黄在线观看 | 99热这里只有精品一区| 村上凉子中文字幕在线| 在线免费观看不下载黄p国产| 成年女人看的毛片在线观看| 亚洲婷婷狠狠爱综合网| 亚洲色图av天堂| 综合色丁香网| 久久久久免费精品人妻一区二区| 国产国拍精品亚洲av在线观看| 久久久久精品国产欧美久久久| 精品久久久久久久久亚洲| 国产日本99.免费观看| 国产高清视频在线播放一区| 一本一本综合久久| 国产一级毛片七仙女欲春2| 久久午夜福利片| 久久精品国产亚洲av涩爱 | 桃色一区二区三区在线观看| 97碰自拍视频| 亚洲人与动物交配视频| 日韩av不卡免费在线播放| 亚洲欧美日韩高清专用| 成人特级黄色片久久久久久久| 国产成人a∨麻豆精品| 熟妇人妻久久中文字幕3abv| 在线观看美女被高潮喷水网站| 人人妻人人看人人澡| 十八禁网站免费在线| 最近最新中文字幕大全电影3| 久久亚洲精品不卡| 亚洲国产精品久久男人天堂| 国产蜜桃级精品一区二区三区| 午夜福利成人在线免费观看| 看非洲黑人一级黄片| 色视频www国产| 日韩人妻高清精品专区| 三级毛片av免费| 身体一侧抽搐| 精品一区二区三区视频在线| 国产精品一区二区性色av| 丝袜美腿在线中文| 国产av麻豆久久久久久久| 波多野结衣高清无吗| 免费看美女性在线毛片视频| 桃色一区二区三区在线观看| 大香蕉久久网| 草草在线视频免费看| 亚洲av熟女| 国产在视频线在精品| 欧美又色又爽又黄视频| 国产精品一区www在线观看| av.在线天堂| 毛片一级片免费看久久久久| 国产高清不卡午夜福利| 美女内射精品一级片tv| 赤兔流量卡办理| 国产伦精品一区二区三区视频9| 如何舔出高潮| 成人毛片a级毛片在线播放| 少妇裸体淫交视频免费看高清| 日韩成人伦理影院| 天美传媒精品一区二区| 日韩精品中文字幕看吧| 国产精品久久久久久av不卡| 人人妻人人澡人人爽人人夜夜 | 中国美白少妇内射xxxbb| 亚洲最大成人av| 国产亚洲精品综合一区在线观看| 亚洲最大成人av| 免费av观看视频| or卡值多少钱| 欧美激情国产日韩精品一区| 岛国在线免费视频观看| 国产黄片美女视频| 久久天躁狠狠躁夜夜2o2o| 免费av不卡在线播放| 性插视频无遮挡在线免费观看| 国产精品国产高清国产av| 69av精品久久久久久| 国产精品国产高清国产av| 99国产精品一区二区蜜桃av| 亚洲专区国产一区二区| 禁无遮挡网站| 久久久久久九九精品二区国产| 人妻少妇偷人精品九色| 亚洲激情五月婷婷啪啪| 美女高潮的动态| 亚洲欧美日韩无卡精品| 校园春色视频在线观看| 一进一出好大好爽视频| 久久久精品大字幕| 亚洲熟妇熟女久久| 成人精品一区二区免费| 俄罗斯特黄特色一大片| 黄色配什么色好看| 欧美成人一区二区免费高清观看| 国产午夜福利久久久久久| 1024手机看黄色片| 亚洲av中文av极速乱| 三级毛片av免费| 日产精品乱码卡一卡2卡三| 亚洲一级一片aⅴ在线观看| 亚洲久久久久久中文字幕| 在线观看一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 国产精品,欧美在线| 国产乱人偷精品视频| 精品人妻视频免费看| 国产精品伦人一区二区| 男人舔女人下体高潮全视频| 国产在视频线在精品| 97超视频在线观看视频| 国产在线精品亚洲第一网站| 精品久久国产蜜桃| 国产午夜精品论理片| 日本爱情动作片www.在线观看 | 禁无遮挡网站| 国产乱人视频| 成人二区视频| 美女内射精品一级片tv| 麻豆乱淫一区二区| 亚洲欧美日韩无卡精品| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久久久久黄片| 99热精品在线国产| 黄色视频,在线免费观看| 欧美绝顶高潮抽搐喷水| 精品久久久久久久末码| 色哟哟哟哟哟哟| 51国产日韩欧美| 婷婷精品国产亚洲av在线| 欧美中文日本在线观看视频| 久久精品夜夜夜夜夜久久蜜豆| 国产精品永久免费网站| 国产精品人妻久久久影院| 亚洲无线在线观看| 免费看日本二区| 无遮挡黄片免费观看| 久久久久性生活片| 免费观看的影片在线观看| 国产aⅴ精品一区二区三区波| 国产亚洲精品av在线| 国产一区二区亚洲精品在线观看| 免费不卡的大黄色大毛片视频在线观看 | 国产高清不卡午夜福利| 日本免费a在线| 午夜日韩欧美国产| 精品乱码久久久久久99久播| 国产aⅴ精品一区二区三区波| 高清日韩中文字幕在线| 国产在视频线在精品| 国产精品一及| 国产亚洲欧美98| 岛国在线免费视频观看| 内地一区二区视频在线| 自拍偷自拍亚洲精品老妇| 午夜激情欧美在线| 一进一出抽搐gif免费好疼| 国产精华一区二区三区| 欧美最新免费一区二区三区| 深爱激情五月婷婷| 午夜久久久久精精品| 成人三级黄色视频| 亚洲成a人片在线一区二区| 此物有八面人人有两片| 欧美zozozo另类| 国产精品99久久久久久久久| 中国国产av一级| 久久久欧美国产精品| 久久精品影院6| 国产大屁股一区二区在线视频| 亚洲精品成人久久久久久| 国产日本99.免费观看| 日韩av不卡免费在线播放| 51国产日韩欧美| 美女免费视频网站| 十八禁国产超污无遮挡网站| 舔av片在线| 国产白丝娇喘喷水9色精品| 男人的好看免费观看在线视频| 天天躁夜夜躁狠狠久久av| 老司机午夜福利在线观看视频| 亚洲综合色惰| 亚洲中文日韩欧美视频| 亚洲国产精品成人综合色| 午夜福利在线观看免费完整高清在 | 亚洲va在线va天堂va国产| 成人毛片a级毛片在线播放| 女人被狂操c到高潮| 色综合亚洲欧美另类图片| 可以在线观看毛片的网站| 观看美女的网站| 国产视频一区二区在线看| 久久久色成人| 亚洲国产精品sss在线观看| av黄色大香蕉| 激情 狠狠 欧美| 黄色一级大片看看| 国产欧美日韩一区二区精品| 欧美日本亚洲视频在线播放| 日韩欧美国产在线观看| 亚洲最大成人av| 亚洲av二区三区四区| 成人av在线播放网站| 亚洲精品日韩av片在线观看| 亚洲国产精品久久男人天堂| 香蕉av资源在线| 黄色欧美视频在线观看| 深爱激情五月婷婷| 亚洲人成网站在线播放欧美日韩| 麻豆乱淫一区二区| 免费av毛片视频| 亚洲av成人精品一区久久| 狂野欧美激情性xxxx在线观看| 免费无遮挡裸体视频| 久久午夜亚洲精品久久| 久久国产乱子免费精品| 久久久精品大字幕| 嫩草影院新地址| 中文亚洲av片在线观看爽| 亚洲欧美中文字幕日韩二区| 97人妻精品一区二区三区麻豆| 日本一本二区三区精品| 欧美日韩乱码在线| 亚洲中文字幕一区二区三区有码在线看| 少妇熟女欧美另类| 亚洲综合色惰| 亚洲国产日韩欧美精品在线观看| 一边摸一边抽搐一进一小说| 日产精品乱码卡一卡2卡三| 国产av一区在线观看免费| 九九热线精品视视频播放| h日本视频在线播放| 国产成人aa在线观看| 日韩精品有码人妻一区| 看片在线看免费视频| 99久久无色码亚洲精品果冻| 国产老妇女一区| 性色avwww在线观看| 特级一级黄色大片| 国产白丝娇喘喷水9色精品| 美女免费视频网站| 能在线免费观看的黄片| 日韩欧美 国产精品| 亚洲成人久久爱视频| av在线观看视频网站免费| 麻豆国产av国片精品| 一区二区三区高清视频在线| 老司机影院成人| 在线天堂最新版资源| 男女视频在线观看网站免费| 国产精品1区2区在线观看.| 欧美成人一区二区免费高清观看| 亚洲成人av在线免费| 插阴视频在线观看视频| 国产午夜福利久久久久久| 在线免费十八禁| 男人舔奶头视频| 久久精品夜夜夜夜夜久久蜜豆| 国产精品永久免费网站| 亚洲自偷自拍三级| 麻豆一二三区av精品| 中国国产av一级| 久久精品国产亚洲网站| 精华霜和精华液先用哪个| 亚洲欧美中文字幕日韩二区| av在线老鸭窝| 亚洲18禁久久av| 日韩一区二区视频免费看| 最近视频中文字幕2019在线8| 伊人久久精品亚洲午夜| 欧美成人一区二区免费高清观看| 成年女人看的毛片在线观看| 国产精品一区www在线观看| 日韩av不卡免费在线播放| 成人二区视频| 久久精品国产清高在天天线| 亚洲中文日韩欧美视频| 两个人视频免费观看高清| 日本黄色视频三级网站网址| 亚洲欧美日韩高清专用| 国产男靠女视频免费网站| 久久精品国产鲁丝片午夜精品| 国产精品三级大全| 免费看光身美女| 国产成人a∨麻豆精品| 国产精品一区二区性色av| 热99re8久久精品国产| 日韩中字成人| 国产 一区 欧美 日韩| 综合色av麻豆| 又黄又爽又免费观看的视频| 中国美女看黄片| 夜夜看夜夜爽夜夜摸| 精品欧美国产一区二区三| 男人舔奶头视频| 色av中文字幕| 久久久久久久久久黄片| 国产在线男女| 深爱激情五月婷婷| 99久久九九国产精品国产免费| 免费观看的影片在线观看| 午夜福利在线在线| 欧美成人免费av一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 国产av一区在线观看免费| 亚洲熟妇中文字幕五十中出| 亚洲婷婷狠狠爱综合网| 日日摸夜夜添夜夜添av毛片| 日韩av在线大香蕉| 插逼视频在线观看| 国产淫片久久久久久久久| 亚洲乱码一区二区免费版| 亚洲精品乱码久久久v下载方式| 亚洲av美国av| 看黄色毛片网站| 欧美激情在线99| 成年女人看的毛片在线观看| 国产麻豆成人av免费视频| 无遮挡黄片免费观看| 在线播放国产精品三级| 少妇的逼好多水| 国产成人一区二区在线|