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

    較大雷諾數(shù)下方柱繞流的數(shù)值模擬

    2017-09-03 09:26:48宋文武
    動(dòng)力工程學(xué)報(bào) 2017年8期
    關(guān)鍵詞:方柱尾緣雷諾數(shù)

    肖 姚, 秦 浩,2, 宋文武,2, 胡 帥

    (1.西華大學(xué) 能源與動(dòng)力工程學(xué)院,成都 610039;2.西華大學(xué) 流體及動(dòng)力機(jī)械教育部重點(diǎn)實(shí)驗(yàn)室,成都 610039)

    較大雷諾數(shù)下方柱繞流的數(shù)值模擬

    肖 姚1, 秦 浩1,2, 宋文武1,2, 胡 帥1

    (1.西華大學(xué) 能源與動(dòng)力工程學(xué)院,成都 610039;2.西華大學(xué) 流體及動(dòng)力機(jī)械教育部重點(diǎn)實(shí)驗(yàn)室,成都 610039)

    基于k-ωSST模型和γ-Reθ轉(zhuǎn)捩模型,分析了雷諾數(shù)為1.4×104和2.2×104時(shí)二維方柱的流場(chǎng)特征,計(jì)算方柱的阻力系數(shù)和斯特勞哈爾數(shù)Sr,并對(duì)方柱平均流場(chǎng)的速度剖面進(jìn)行分析,綜合考慮不同湍流模型的計(jì)算成本及計(jì)算精度.結(jié)果表明:γ-Reθ轉(zhuǎn)捩模型比k-ωSST模型計(jì)算結(jié)果更加精確且所需計(jì)算成本較?。沪?Reθ轉(zhuǎn)捩模型能較好地預(yù)測(cè)方柱尾緣近壁回流區(qū)的位置和大小,且能較好地觀測(cè)流場(chǎng)在方柱尾緣的再附著現(xiàn)象;在求解實(shí)際大雷諾數(shù)工程問題中γ-Reθ轉(zhuǎn)捩模型為優(yōu)秀的湍流模型.

    方柱繞流; 數(shù)值模擬; 阻力系數(shù); 湍流模型; 時(shí)間步長(zhǎng)

    鈍體繞流廣泛存在于流體機(jī)械、動(dòng)力工程和石油工程等領(lǐng)域中,如汽輪機(jī)葉片、換熱器管群及海洋鉆井平臺(tái)等均涉及鈍體繞流問題.鈍體繞流所引起的振動(dòng)會(huì)造成設(shè)備工作不穩(wěn)定,當(dāng)振動(dòng)頻率與結(jié)構(gòu)的固有頻率相接近時(shí)會(huì)產(chǎn)生共振,影響結(jié)構(gòu)的使用性能,振動(dòng)嚴(yán)重時(shí)甚至?xí)?duì)結(jié)構(gòu)造成破壞性損壞.方柱繞流存在復(fù)雜的渦結(jié)構(gòu),而渦的脫落對(duì)于方柱所受的升力和阻力均有較大影響.因此,深入研究方柱繞流對(duì)解決工程實(shí)際中相關(guān)問題有著重要意義.

    為了準(zhǔn)確高效地計(jì)算出方柱流場(chǎng)的分布以及了解渦的形成與脫落規(guī)律,很多學(xué)者對(duì)方柱繞流開展了研究工作[1-10].如Durao等[1]在水洞中利用多普勒激光測(cè)速儀(LDV)研究了方柱在雷諾數(shù)為1.4×104時(shí)的繞流問題,結(jié)果表明隨著湍流的波動(dòng),在顫振速度最大區(qū)域的能量可達(dá)總能量的40%.Lyn等[2]基于LDV實(shí)驗(yàn)測(cè)試了雷諾數(shù)為2.2×104時(shí)方柱的流場(chǎng),得出方柱尾緣的速度分布規(guī)律以及渦的分布與大小.張偉等[3]采用粒子圖像測(cè)速儀(PIV)和k-ωSST模型研究了雷諾數(shù)為2.2×104時(shí)的方柱繞流流場(chǎng),得出k-ωSST模型計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)吻合較好.Bouris等[4]利用大渦模型數(shù)值模擬了小雷諾數(shù)(Re=250)時(shí)的二維方柱繞流,通過優(yōu)化網(wǎng)格質(zhì)量得出二維方柱的大渦模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)吻合較好.梁思超等[6]基于分離渦(DES)模型研究了雷諾數(shù)為4×104時(shí)方柱繞流中馬蹄渦系統(tǒng)的發(fā)展過程.Franke等[7]采用k-ε渦黏模型和雷諾應(yīng)力模型研究了雷諾數(shù)為2.2×104時(shí)的二維方柱流場(chǎng),結(jié)果表明2種模型均不能較好的與實(shí)驗(yàn)結(jié)果吻合.謝志剛等[8]基于大渦模型,利用精細(xì)的邊界研究雷諾數(shù)為2.2×104時(shí)的方柱繞流,其結(jié)果可以媲美密網(wǎng)格的結(jié)果.

    研究者們從實(shí)驗(yàn)、數(shù)值模擬和理論分析3個(gè)方面同時(shí)研究了方柱的流場(chǎng),然而實(shí)驗(yàn)結(jié)果與數(shù)值模擬還存在一定的差異.文獻(xiàn)[7]的研究結(jié)果與實(shí)驗(yàn)結(jié)果相差較大,k-ε渦黏模型計(jì)算出的Sr和阻力因數(shù)CD均偏??;文獻(xiàn)[3]的結(jié)果雖然與文獻(xiàn)[2]的實(shí)驗(yàn)結(jié)果吻合較好,但其計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果在近方柱的尾流區(qū)存在一個(gè)明顯的跳躍,其對(duì)尾流區(qū)未能準(zhǔn)確預(yù)測(cè);文獻(xiàn)[4]的計(jì)算結(jié)果較好,但計(jì)算的雷諾數(shù)較小且由于該湍流模型計(jì)算成本較高,將其應(yīng)用在實(shí)際工程中不現(xiàn)實(shí).實(shí)際工程多涉及大雷諾數(shù)下渦街的形成與脫落,為此需在考慮成本的情況下得出適合于工程大雷諾數(shù)的湍流模型.γ-Reθ轉(zhuǎn)捩模型[11]基于k-ωSST輸運(yùn)方程、間歇系數(shù)γ方程和動(dòng)量厚度雷諾數(shù)Reθt方程,相比k-ωSST模型,該模型可以較好地捕捉流場(chǎng)分離過程和湍流強(qiáng)度的變化.筆者基于Fluent軟件分析方柱在不同Re下的流場(chǎng),通過對(duì)比分析方柱的Sr、阻力系數(shù)CD、升力系數(shù)CL和流場(chǎng)的速度剖面等參數(shù)來(lái)對(duì)比不同湍流模型的評(píng)估結(jié)果.

    1 數(shù)值方法

    1.1 控制方程

    二維黏性不可壓牛頓流體運(yùn)動(dòng)的基本控制方程為連續(xù)性方程和Navier-Stokes方程:

    (1)

    (2)

    1.2k-ωSST模型

    k-ωSST模型基于Menter[12]研究成果,該湍流模型在多種流態(tài)下較standardk-ω和BSLk-ω湍流模型更加精確和可靠,其表達(dá)式如下:

    Gk-Yk+Sk

    (3)

    Gw-Yw+Dw+Sw

    (4)

    式中:Gk、Gw分別代表湍動(dòng)能和w的產(chǎn)生;Γk、Γw分別為k和w的有效擴(kuò)散率;Yk、Yw分別為由湍動(dòng)產(chǎn)生的耗散率;Sk、Sw為源項(xiàng).

    1.3γ-Reθ轉(zhuǎn)捩模型

    Menter等[11]提出了γ-Reθ轉(zhuǎn)捩模型,該湍流模型由k-ωSST輸運(yùn)方程、間歇系數(shù)γ方程和動(dòng)量厚度雷諾數(shù)Reθt方程組成,其中間歇系數(shù)γ用來(lái)描述邊界層的狀態(tài),并控制過渡點(diǎn)下游湍動(dòng)能的生成項(xiàng),式(5)為間歇系數(shù)無(wú)量綱輸運(yùn)方程的守恒形式.動(dòng)量厚度雷諾數(shù)Reθt通過捕獲非局部影響的湍流強(qiáng)度來(lái)控制由邊界層內(nèi)間歇系數(shù)的增長(zhǎng)和當(dāng)?shù)販u量雷諾數(shù)構(gòu)成的轉(zhuǎn)捩觸發(fā)準(zhǔn)則,其表達(dá)式如下:

    (5)

    (6)

    式中:μt、σθ、σγ為模型系數(shù);Pγ1、Eγ1代表轉(zhuǎn)捩源項(xiàng);Pγ2、Eγ2分別為破裂系數(shù)和再附著系數(shù);Pθ為源項(xiàng).

    1.4 方柱幾何模型與網(wǎng)格劃分

    圖1 方柱流場(chǎng)計(jì)算域

    圖2 網(wǎng)格劃分

    采用Fluent軟件進(jìn)行數(shù)值模擬,設(shè)定初始條件:u1=U0,u2=0;方柱入口定義為速度邊界條件,來(lái)流為均勻速度U0;出口定義為自由流邊界條件,方柱壁面與計(jì)算域邊界定義為無(wú)滑移壁面.

    2 結(jié)果與討論

    2.1 湍流模型對(duì)比

    首先以文獻(xiàn)[1]和文獻(xiàn)[2]的實(shí)驗(yàn)結(jié)果作為數(shù)值模擬的參照對(duì)象,計(jì)算不同Re下各湍流模型最佳的模擬時(shí)間步長(zhǎng),得出方柱阻力系數(shù)和升力系數(shù)的時(shí)程曲線,將升力系數(shù)時(shí)程曲線進(jìn)行快速傅氏變換(FFT)得到渦街脫落頻率f.其中阻力系數(shù)CD、升力系數(shù)CL及Sr的定義分別為:

    (a) 方柱上壁面Y+

    (b) 方柱尾緣壁面Y+

    (7)

    式中:FD和FL分別為方柱所受的阻力和升力.

    表1 不同Re下方柱繞流的Sr和

    (a) k-ω SST模型

    (b) γ-Reθ轉(zhuǎn)捩模型

    2.2 流場(chǎng)分析

    由第2.1節(jié)可知γ-Reθ轉(zhuǎn)捩模型可以較好地評(píng)估實(shí)際工程大Re下的流場(chǎng)變化,為了與實(shí)驗(yàn)進(jìn)行對(duì)比,下文只分析Re為2.2×104時(shí)方柱尾緣渦的脫落過程.圖5為方柱平均流場(chǎng)速度分布曲線與水平湍動(dòng)能變化曲線.從圖5(a)可以看出,2種湍流模型計(jì)算結(jié)果在方柱前緣的速度場(chǎng)與文獻(xiàn)[1]的實(shí)驗(yàn)結(jié)果吻合.在1.8

    (a) 水平速度剖面

    (b) 豎向速度剖面

    (c) 水平剖面湍動(dòng)能

    綜上,γ-Reθ轉(zhuǎn)捩模型與三維大渦模擬結(jié)果以及實(shí)驗(yàn)結(jié)果較為接近,但由于尾流近壁回流區(qū)是三維流動(dòng),本文結(jié)果還存在一定的偏差.

    2.3 尾渦脫落

    由表1可知,雷諾數(shù)為2.2×104時(shí)采用γ-Reθ轉(zhuǎn)捩模型計(jì)算得出的方柱渦脫落頻率為17 Hz,周期為0.059 s.圖6給出了Re為2.2×104時(shí)方柱尾緣渦脫落的半個(gè)周期過程,其中τ為時(shí)間.從圖6可以看出,在方柱前緣-0.05 m左右流場(chǎng)開始分離,并在方柱兩側(cè)產(chǎn)生回流區(qū),在方柱尾緣也有短暫再附著現(xiàn)象.在這半個(gè)周期內(nèi)尾渦從方柱尾緣上方逐漸形成,并在尾緣擴(kuò)大,在一定程度時(shí)渦開始脫落,其形成尾渦的大小與方柱的特征長(zhǎng)度相近.

    (a) τ=0 s

    (b) τ=0.000 2 s

    (c) τ=0.000 4 s

    (d) τ=0.000 6 s

    (e) τ=0.000 8 s

    (f) τ=0.001 s

    (g) τ=0.001 2 s

    (h) τ=0.001 4 s

    (i) τ=0.001 6 s

    (j) τ=0.001 8 s

    圖6 方柱尾渦脫落過程

    Fig.6 Shedding process of vortex from the square cylinder

    3 結(jié) 論

    (1) 與其他湍流模型相比,γ-Reθ轉(zhuǎn)捩模型可以精確地評(píng)估大Re下方柱的Sr、CD和CL.

    (2) 與相關(guān)實(shí)驗(yàn)結(jié)果對(duì)比,γ-Reθ轉(zhuǎn)捩模型與實(shí)驗(yàn)結(jié)果的速度剖面吻合較好,可以較好地預(yù)測(cè)方柱尾緣回流區(qū)的大小和位置,計(jì)算結(jié)果與三維大渦模擬結(jié)果較為接近.

    (3)γ-Reθ轉(zhuǎn)捩模型可以準(zhǔn)確地分析方柱壁面渦的形成與脫落過程以及流場(chǎng)再附著現(xiàn)象,很好地揭示了繞流的流動(dòng)機(jī)理.

    [1] DURAO D, HEITOR M V, PEREIRA J C F. Measurements of turbulent and periodic flows around a square cross-section cylinder[J]. Experiments in Fluids, 1988, 6(5): 298-304.

    [2] LYN D A, EINAV S, RODI W, et al. A laser-Doppler velocimetry study of ensemble-averaged characteristics of the turbulent near wake of a square cylinder[J]. Journal of Fluid Mechanics, 1995, 304: 285-319.

    [3] 張偉, 葛耀君. 方柱繞流粒子圖像測(cè)速試驗(yàn)與數(shù)值模擬[J]. 同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2009, 37(7): 857-861, 892.

    ZHANG Wei, GE Yaojun. Particle image velocimetry study and numerical simulation of turbulent near wake of square cylinder[J]. Journal of Tongji University(Natural Science), 2009, 37(7): 857-861, 892.

    [4] BOURIS D, BERGELES G. 2D LES of vortex shedding from a square cylinder[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1999, 80(1/2): 31-46.

    [5] 李慧君, 王炯. 基于水平圓管外加裝排液板降膜流動(dòng)的數(shù)值模擬[J]. 動(dòng)力工程學(xué)報(bào), 2016, 36(4): 265-270.

    LI Huijun, WANG Jiong. Numerical simulation of falling film flow on horizontal circular tubes with drainage strips[J]. Journal of Chinese Society of Power Engineering, 2016, 36(4): 265-270.

    [6] 梁思超, 張曉東, 康順, 等. 方體繞流的DES模擬及其旋渦結(jié)構(gòu)分析[J]. 工程熱物理學(xué)報(bào), 2015, 36(3): 535-540.

    LIANG Sichao, ZHANG Xiaodong, KANG Shun, et al. DES simulation of the flow around a cube and analysis of the vortex structure[J]. Journal of Engineering Thermophysics, 2015, 36(3): 535-540.

    [7] FRANKE R, RODI W. Calculation of vortex shedding past a square cylinder with various turbulence models[M]//DURST F, FRIEDRICH R, LAUNDER BE, et al. Turbulent Shear Flows 8.Berlin Heidelberg:Springer, 1993.

    [8] 謝志剛, 許春曉, 崔桂香, 等. 方柱繞流大渦模擬[J]. 計(jì)算物理, 2007, 24(2): 171-180.

    XIE Zhigang, XU Chunxiao, CUI Guixiang, et al. Large eddy simulation of flows around a square cylinder[J]. Chinese Journal of Computational Physics, 2007, 24(2): 171-180.

    [9] OKAJIMA A. Strouhal numbers of rectangular cylinders[J]. Journal of Fluid Mechanics, 1982, 123: 379-398.

    [10] SAMANI M, BERGSTROM D J. Effect of a wall on the wake dynamics of an infinite square cylinder[J]. International Journal of Heat and Fluid Flow, 2015, 55: 158-166.

    [11] MENTER F R, LANGTRY R B, LIKKI S R, et al. A correlation-based transition model using local variables—Part I: model formulation[J]. Journal of Turbomachinery, 2006, 128(3): 413-422.

    [12] MENTER F R. Two-equation eddy-viscosity turbulence models for engineering applications[J]. AIAA Journal, 1994, 32(8): 1598-1605.

    [13] BREUER M, POURQUIé M. First experiences with LES of flows past bluff bodies[C]//Proceedings of the 3rd International Symposium of Engineering Turbulence Modelling and Measurements. Amsterdam,Holland:Elsevier,1996: 177-186.

    [14] SAHA A K, MURALIDHAR K, BISWAS G. Experimental study of flow past a square cylinder at high Reynolds numbers[J]. Experiments in Fluids, 2000, 29(6): 553-563.

    Numerical Simulation of the Flow Around a Square Cylinder at Higher Reynolds Numbers

    XIAOYao1,QINHao1,2,SONGWenwu1,2,HUShuai1

    (1.School of Energy and Power Engineering, Xihua University, Chengdu 610039, China; 2. Key Laboratory of Fluid and Power Machinery, Ministry of Education, Xihua University, Chengdu 610039, China)

    Based on thek-ωSST model andγ-Reθtransition model, a study was conducted on the flow characteristics around a square cylinder respectively at the Reynolds number of 1.4×104and 2.2×104, so as to calculate the drag coefficient and Strouhal number of the square cylinder, analyze its average velocity profiles, and finally to find out the calculation cost and precision of different turbulent models. Results show that compared withk-ωSST model, theγ-Reθtransition model is more superior in calculation precision and cost, which is more effective in predicting the position and magnitude of the recirculation region behind the square cylinder, and in analyzing the reattachment of separated flows on side walls near the rear corner, and therefore is a practical model for engineering applications at high Reynolds numbers.

    flow around square cylinder; numerical simulation; drag coefficient; turbulent model; time step

    1674-7607(2017)08-0679-06

    O357

    A

    470.10

    2016-08-04

    2016-10-24

    流體及動(dòng)力機(jī)械教育部重點(diǎn)實(shí)驗(yàn)室開放基金資助項(xiàng)目(szjj2016-003);西華大學(xué)自然科學(xué)重點(diǎn)基金資助項(xiàng)目(Z1510418);國(guó)家自然科學(xué)基金資助項(xiàng)目(51408505)

    肖 姚(1990-),男,四川資陽(yáng)人,碩士研究生,研究方向?yàn)榱黧w機(jī)械. 秦 浩(通信作者),男,講師,博士,電話(Tel.):13547904041;E-mail:qh.email.@foxmail.com.

    猜你喜歡
    方柱尾緣雷諾數(shù)
    上游切角倒角小間距比串列方柱大渦模擬研究
    串列多方柱氣動(dòng)特性的試驗(yàn)研究
    基于強(qiáng)化換熱的偏斜尾緣設(shè)計(jì)
    能源工程(2021年1期)2021-04-13 02:05:46
    基于Transition SST模型的高雷諾數(shù)圓柱繞流數(shù)值研究
    翼型湍流尾緣噪聲半經(jīng)驗(yàn)預(yù)測(cè)公式改進(jìn)
    具有尾緣襟翼的風(fēng)力機(jī)動(dòng)力學(xué)建模與恒功率控制
    2017年中考數(shù)學(xué)模擬試題(十)
    失穩(wěn)初期的低雷諾數(shù)圓柱繞流POD-Galerkin 建模方法研究
    基于轉(zhuǎn)捩模型的低雷諾數(shù)翼型優(yōu)化設(shè)計(jì)研究
    民機(jī)高速風(fēng)洞試驗(yàn)的阻力雷諾數(shù)效應(yīng)修正
    三级毛片av免费| 亚洲午夜理论影院| 在线观看免费视频日本深夜| 老熟妇仑乱视频hdxx| 国产私拍福利视频在线观看| 久久天躁狠狠躁夜夜2o2o| 久久久国产成人精品二区| 日韩大尺度精品在线看网址 | 国产精品 国内视频| 亚洲免费av在线视频| 男男h啪啪无遮挡| 老汉色∧v一级毛片| 美女午夜性视频免费| av在线天堂中文字幕| 可以在线观看的亚洲视频| 欧美久久黑人一区二区| 欧美日韩黄片免| 他把我摸到了高潮在线观看| 亚洲一码二码三码区别大吗| 欧美一级毛片孕妇| 国产精品国产高清国产av| 村上凉子中文字幕在线| 琪琪午夜伦伦电影理论片6080| 久久性视频一级片| 夜夜躁狠狠躁天天躁| 欧美日韩瑟瑟在线播放| 99re在线观看精品视频| 亚洲黑人精品在线| 日韩欧美国产在线观看| 动漫黄色视频在线观看| 日韩三级视频一区二区三区| 久久午夜亚洲精品久久| 中文字幕av电影在线播放| 国产午夜福利久久久久久| 少妇被粗大的猛进出69影院| 日本 欧美在线| 午夜福利成人在线免费观看| 成人永久免费在线观看视频| 美女免费视频网站| 一夜夜www| 老司机午夜福利在线观看视频| 欧美日本视频| 日本五十路高清| 亚洲无线在线观看| 国产av又大| 欧美一级a爱片免费观看看 | 黄网站色视频无遮挡免费观看| 两个人免费观看高清视频| 69av精品久久久久久| 久久久久国产一级毛片高清牌| 香蕉久久夜色| 村上凉子中文字幕在线| 国产成人精品久久二区二区91| 国产精品亚洲av一区麻豆| 久久精品国产综合久久久| 亚洲精品在线观看二区| 淫秽高清视频在线观看| 热99re8久久精品国产| 一级毛片精品| 国产三级在线视频| 日韩欧美免费精品| 婷婷丁香在线五月| 欧美精品啪啪一区二区三区| 国产av一区二区精品久久| 免费搜索国产男女视频| 亚洲天堂国产精品一区在线| 国产一区二区三区视频了| 久久天堂一区二区三区四区| 欧美色视频一区免费| 亚洲欧美日韩高清在线视频| 欧美一区二区精品小视频在线| 国产一区二区在线av高清观看| 国产欧美日韩一区二区精品| 一级毛片高清免费大全| 搡老岳熟女国产| 国产高清有码在线观看视频 | 国产高清有码在线观看视频 | 最近最新免费中文字幕在线| www.精华液| 侵犯人妻中文字幕一二三四区| 他把我摸到了高潮在线观看| 十八禁人妻一区二区| 久久影院123| 精品国产亚洲在线| 一个人免费在线观看的高清视频| 极品人妻少妇av视频| 黑人操中国人逼视频| 最好的美女福利视频网| 免费搜索国产男女视频| 九色国产91popny在线| 身体一侧抽搐| 久久中文字幕一级| 丰满人妻熟妇乱又伦精品不卡| 变态另类成人亚洲欧美熟女 | 波多野结衣巨乳人妻| 日韩精品青青久久久久久| 黑丝袜美女国产一区| 国产单亲对白刺激| 黑人巨大精品欧美一区二区mp4| 久久中文看片网| 午夜福利视频1000在线观看 | 久久人妻福利社区极品人妻图片| 在线播放国产精品三级| 亚洲 欧美 日韩 在线 免费| 久久久久久大精品| 日韩欧美在线二视频| 国产精品亚洲一级av第二区| 色综合站精品国产| 少妇的丰满在线观看| 成年版毛片免费区| 久久天堂一区二区三区四区| 精品一区二区三区四区五区乱码| 欧美 亚洲 国产 日韩一| 国产不卡一卡二| 最近最新免费中文字幕在线| 精品一区二区三区av网在线观看| 国产欧美日韩一区二区三区在线| 亚洲av电影在线进入| av片东京热男人的天堂| 国产精品免费视频内射| 亚洲中文av在线| 国产亚洲欧美98| 成年版毛片免费区| 午夜福利,免费看| 久久香蕉激情| 欧美性长视频在线观看| 亚洲视频免费观看视频| 操出白浆在线播放| 欧美精品啪啪一区二区三区| 天堂√8在线中文| 久久这里只有精品19| 中文字幕精品免费在线观看视频| 欧美黄色片欧美黄色片| 免费在线观看亚洲国产| 一二三四在线观看免费中文在| 亚洲,欧美精品.| 久9热在线精品视频| 男男h啪啪无遮挡| 精品久久久久久,| 亚洲国产精品成人综合色| 国产精品一区二区在线不卡| 一进一出抽搐gif免费好疼| 日韩视频一区二区在线观看| 色av中文字幕| 久久久久久人人人人人| 欧美日韩乱码在线| 久久久久久久久免费视频了| 亚洲国产精品sss在线观看| 老鸭窝网址在线观看| 精品第一国产精品| 日韩精品青青久久久久久| 国产高清有码在线观看视频 | 给我免费播放毛片高清在线观看| 欧美乱码精品一区二区三区| 岛国视频午夜一区免费看| x7x7x7水蜜桃| 亚洲av五月六月丁香网| 视频区欧美日本亚洲| 亚洲成国产人片在线观看| 精品福利观看| 女性生殖器流出的白浆| 99久久综合精品五月天人人| 国产成人欧美| 亚洲一码二码三码区别大吗| 看免费av毛片| 久久狼人影院| av福利片在线| 侵犯人妻中文字幕一二三四区| 欧美激情极品国产一区二区三区| www国产在线视频色| 国产精品一区二区精品视频观看| 成人免费观看视频高清| 成人18禁在线播放| 男人舔女人的私密视频| 久久天躁狠狠躁夜夜2o2o| 一级毛片精品| 一进一出好大好爽视频| 两人在一起打扑克的视频| 亚洲色图av天堂| 天堂√8在线中文| 亚洲成国产人片在线观看| 精品少妇一区二区三区视频日本电影| 久久精品人人爽人人爽视色| 欧美日韩亚洲国产一区二区在线观看| 麻豆久久精品国产亚洲av| 黑丝袜美女国产一区| 成人18禁在线播放| 亚洲自拍偷在线| 18禁裸乳无遮挡免费网站照片 | 色老头精品视频在线观看| 天堂动漫精品| 日本三级黄在线观看| 12—13女人毛片做爰片一| 欧美日韩亚洲综合一区二区三区_| 成人18禁在线播放| 亚洲伊人色综图| cao死你这个sao货| 亚洲国产欧美网| 国产高清激情床上av| 可以在线观看的亚洲视频| 激情在线观看视频在线高清| 天天添夜夜摸| 免费观看精品视频网站| 日韩欧美一区视频在线观看| 国产熟女午夜一区二区三区| 国产精品98久久久久久宅男小说| 国产精品av久久久久免费| 国产主播在线观看一区二区| 国产精品秋霞免费鲁丝片| 亚洲精品一卡2卡三卡4卡5卡| 亚洲久久久国产精品| 久久草成人影院| 黄色 视频免费看| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品国产色婷婷电影| 大型黄色视频在线免费观看| 一个人观看的视频www高清免费观看 | 12—13女人毛片做爰片一| 久久精品国产亚洲av高清一级| 久久久久久久久免费视频了| 两个人免费观看高清视频| 国产成人av教育| 亚洲色图av天堂| 90打野战视频偷拍视频| 日韩大尺度精品在线看网址 | 50天的宝宝边吃奶边哭怎么回事| 一级a爱片免费观看的视频| 成年人黄色毛片网站| 99国产精品一区二区三区| 亚洲熟妇中文字幕五十中出| www.精华液| 神马国产精品三级电影在线观看 | 18禁美女被吸乳视频| 黄色a级毛片大全视频| 亚洲国产精品sss在线观看| 男人操女人黄网站| 夜夜看夜夜爽夜夜摸| 欧美激情高清一区二区三区| 一级毛片精品| 国产精品免费一区二区三区在线| 精品久久蜜臀av无| 午夜福利高清视频| 叶爱在线成人免费视频播放| √禁漫天堂资源中文www| 一区二区三区国产精品乱码| 欧美日韩精品网址| 乱人伦中国视频| 国产亚洲精品第一综合不卡| 午夜免费观看网址| 国产三级黄色录像| 国产精品久久视频播放| 法律面前人人平等表现在哪些方面| 又黄又爽又免费观看的视频| 国产三级在线视频| 国产精品爽爽va在线观看网站 | av网站免费在线观看视频| 亚洲男人的天堂狠狠| 免费无遮挡裸体视频| 宅男免费午夜| 久久久国产成人精品二区| 一区二区三区高清视频在线| 午夜精品在线福利| 欧美成人午夜精品| av超薄肉色丝袜交足视频| 女人爽到高潮嗷嗷叫在线视频| 99国产精品免费福利视频| 色av中文字幕| 免费少妇av软件| 久99久视频精品免费| 一级a爱视频在线免费观看| 国产精品精品国产色婷婷| 午夜免费鲁丝| 免费在线观看完整版高清| 熟女少妇亚洲综合色aaa.| 韩国精品一区二区三区| 操出白浆在线播放| 每晚都被弄得嗷嗷叫到高潮| 国产成人影院久久av| 日韩国内少妇激情av| 国产精品日韩av在线免费观看 | 亚洲精品美女久久av网站| 在线播放国产精品三级| 日韩精品免费视频一区二区三区| 久久精品国产99精品国产亚洲性色 | 国产私拍福利视频在线观看| 亚洲第一电影网av| 无遮挡黄片免费观看| 国产麻豆成人av免费视频| 亚洲人成网站在线播放欧美日韩| 久久久久久久精品吃奶| 欧美 亚洲 国产 日韩一| 成人亚洲精品一区在线观看| 国产激情欧美一区二区| 精品国产美女av久久久久小说| 日本vs欧美在线观看视频| 麻豆久久精品国产亚洲av| 久久精品成人免费网站| 国产乱人伦免费视频| 18禁美女被吸乳视频| 国产亚洲精品综合一区在线观看 | 婷婷六月久久综合丁香| 午夜久久久在线观看| 涩涩av久久男人的天堂| 国产成人av教育| 久久婷婷成人综合色麻豆| 一级,二级,三级黄色视频| 熟妇人妻久久中文字幕3abv| 俄罗斯特黄特色一大片| 国产亚洲精品一区二区www| 国产亚洲欧美在线一区二区| 国产精品野战在线观看| 国产精品日韩av在线免费观看 | 免费观看人在逋| 久久 成人 亚洲| 在线免费观看的www视频| 亚洲精品国产一区二区精华液| 无遮挡黄片免费观看| 精品一区二区三区四区五区乱码| 亚洲av成人不卡在线观看播放网| 69精品国产乱码久久久| 99久久久亚洲精品蜜臀av| 亚洲人成电影观看| 国产乱人伦免费视频| 满18在线观看网站| 两性午夜刺激爽爽歪歪视频在线观看 | 精品电影一区二区在线| 青草久久国产| 亚洲人成网站在线播放欧美日韩| 视频区欧美日本亚洲| 午夜福利成人在线免费观看| 麻豆成人av在线观看| 免费观看人在逋| 我的亚洲天堂| 国产又爽黄色视频| 免费看十八禁软件| 搞女人的毛片| 日韩精品免费视频一区二区三区| x7x7x7水蜜桃| 一区二区三区激情视频| 国产午夜精品久久久久久| 欧美国产日韩亚洲一区| 97碰自拍视频| 夜夜爽天天搞| 一区二区三区国产精品乱码| 美女免费视频网站| 亚洲人成电影免费在线| 国产成+人综合+亚洲专区| 美女大奶头视频| 日韩大码丰满熟妇| 韩国av一区二区三区四区| 两个人免费观看高清视频| 国产国语露脸激情在线看| 欧美国产精品va在线观看不卡| 亚洲成av片中文字幕在线观看| 夜夜夜夜夜久久久久| 亚洲成人免费电影在线观看| 成在线人永久免费视频| 日韩一卡2卡3卡4卡2021年| 亚洲一区二区三区不卡视频| 他把我摸到了高潮在线观看| 国内精品久久久久久久电影| 夜夜看夜夜爽夜夜摸| 免费观看人在逋| 高潮久久久久久久久久久不卡| 久久精品国产亚洲av香蕉五月| 熟女少妇亚洲综合色aaa.| 日韩精品免费视频一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 久久精品aⅴ一区二区三区四区| 自拍欧美九色日韩亚洲蝌蚪91| 国产在线观看jvid| 国产成人啪精品午夜网站| 国产真人三级小视频在线观看| 色综合欧美亚洲国产小说| 免费少妇av软件| 精品午夜福利视频在线观看一区| 国产精品香港三级国产av潘金莲| 午夜成年电影在线免费观看| 亚洲激情在线av| 91大片在线观看| 午夜福利,免费看| 亚洲在线自拍视频| 老汉色av国产亚洲站长工具| 99国产精品99久久久久| 丝袜美腿诱惑在线| 天天一区二区日本电影三级 | 一边摸一边做爽爽视频免费| 亚洲成av人片免费观看| 一进一出抽搐gif免费好疼| 女警被强在线播放| 91成人精品电影| 日本vs欧美在线观看视频| 国产三级在线视频| 亚洲avbb在线观看| 一区二区日韩欧美中文字幕| 级片在线观看| 亚洲男人天堂网一区| 日韩大尺度精品在线看网址 | 亚洲一区中文字幕在线| 免费看十八禁软件| 欧美亚洲日本最大视频资源| 免费在线观看黄色视频的| 日韩精品青青久久久久久| 国产精品秋霞免费鲁丝片| 精品人妻1区二区| 亚洲成国产人片在线观看| 精品国产一区二区三区四区第35| 黄色丝袜av网址大全| 免费在线观看完整版高清| 怎么达到女性高潮| 嫩草影院精品99| 一本综合久久免费| 亚洲av熟女| 一级黄色大片毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 波多野结衣高清无吗| 亚洲欧美激情综合另类| 中文字幕人成人乱码亚洲影| 制服人妻中文乱码| 女生性感内裤真人,穿戴方法视频| 国产成人免费无遮挡视频| 久久久国产精品麻豆| 99久久99久久久精品蜜桃| 国产伦人伦偷精品视频| 一级a爱视频在线免费观看| 桃红色精品国产亚洲av| 亚洲全国av大片| 亚洲人成伊人成综合网2020| 免费一级毛片在线播放高清视频 | www.www免费av| 成年版毛片免费区| 欧美成狂野欧美在线观看| 一边摸一边做爽爽视频免费| av福利片在线| 波多野结衣av一区二区av| 国产激情久久老熟女| 国产精品电影一区二区三区| 亚洲第一青青草原| 亚洲色图综合在线观看| 大陆偷拍与自拍| 国产精品乱码一区二三区的特点 | 悠悠久久av| 18禁黄网站禁片午夜丰满| 婷婷精品国产亚洲av在线| www国产在线视频色| 婷婷六月久久综合丁香| 免费观看精品视频网站| 亚洲中文日韩欧美视频| 久久久久久久午夜电影| 欧美国产精品va在线观看不卡| 一级,二级,三级黄色视频| 国产午夜福利久久久久久| 黄色片一级片一级黄色片| 涩涩av久久男人的天堂| 成熟少妇高潮喷水视频| 1024香蕉在线观看| 精品高清国产在线一区| 久久精品国产亚洲av高清一级| 少妇 在线观看| 禁无遮挡网站| 婷婷精品国产亚洲av在线| 成人三级做爰电影| 亚洲av美国av| 国产成人免费无遮挡视频| 国产黄a三级三级三级人| xxx96com| 欧美老熟妇乱子伦牲交| 一进一出抽搐动态| 性欧美人与动物交配| 在线观看66精品国产| 狠狠狠狠99中文字幕| 国产成人免费无遮挡视频| 人人妻人人爽人人添夜夜欢视频| 午夜久久久久精精品| 国产精品免费一区二区三区在线| 制服人妻中文乱码| 亚洲精品中文字幕在线视频| 国产精品久久视频播放| 一二三四在线观看免费中文在| 91av网站免费观看| 一本大道久久a久久精品| 久久精品国产亚洲av高清一级| 国产伦一二天堂av在线观看| 国产亚洲精品久久久久5区| 一本大道久久a久久精品| 欧美日韩乱码在线| 久久精品亚洲精品国产色婷小说| 欧美日韩乱码在线| 在线十欧美十亚洲十日本专区| 男女之事视频高清在线观看| 无限看片的www在线观看| 国产精品精品国产色婷婷| 亚洲一区二区三区不卡视频| 日韩欧美免费精品| 成人国产综合亚洲| 中国美女看黄片| 身体一侧抽搐| 午夜精品在线福利| 少妇粗大呻吟视频| 热99re8久久精品国产| 国产av在哪里看| 免费在线观看黄色视频的| 亚洲 国产 在线| 在线观看日韩欧美| 亚洲专区中文字幕在线| 日韩有码中文字幕| 亚洲一区二区三区色噜噜| 久久九九热精品免费| 日韩中文字幕欧美一区二区| 这个男人来自地球电影免费观看| 欧美精品亚洲一区二区| 一夜夜www| 9191精品国产免费久久| 国产成人精品在线电影| 一区二区三区国产精品乱码| 色综合站精品国产| 在线国产一区二区在线| 天天添夜夜摸| 国产成+人综合+亚洲专区| 午夜亚洲福利在线播放| 在线十欧美十亚洲十日本专区| 高清黄色对白视频在线免费看| 两性午夜刺激爽爽歪歪视频在线观看 | 美女 人体艺术 gogo| 乱人伦中国视频| 女人精品久久久久毛片| av中文乱码字幕在线| 黄色丝袜av网址大全| 午夜福利免费观看在线| 麻豆成人av在线观看| 亚洲欧美激情综合另类| 国产精品,欧美在线| 男男h啪啪无遮挡| 国产精品 国内视频| 嫩草影视91久久| www.精华液| 久久精品人人爽人人爽视色| 黄色成人免费大全| 日本vs欧美在线观看视频| 亚洲激情在线av| 精品久久久精品久久久| 亚洲成av人片免费观看| 啦啦啦免费观看视频1| 久久青草综合色| 69av精品久久久久久| 18禁观看日本| 伦理电影免费视频| 久久中文字幕一级| 国产一区二区三区视频了| 欧美一级毛片孕妇| 国产一区二区激情短视频| av免费在线观看网站| 国内精品久久久久久久电影| 欧美老熟妇乱子伦牲交| 搞女人的毛片| 别揉我奶头~嗯~啊~动态视频| avwww免费| 亚洲中文日韩欧美视频| 日韩大码丰满熟妇| 亚洲国产精品999在线| 夜夜躁狠狠躁天天躁| 国产成人精品久久二区二区91| 午夜精品国产一区二区电影| 亚洲国产精品久久男人天堂| 成在线人永久免费视频| 国产一级毛片七仙女欲春2 | 亚洲一区中文字幕在线| 国产国语露脸激情在线看| 他把我摸到了高潮在线观看| 久久精品亚洲精品国产色婷小说| tocl精华| 国产主播在线观看一区二区| 91在线观看av| 如日韩欧美国产精品一区二区三区| 99re在线观看精品视频| 男男h啪啪无遮挡| 黄片大片在线免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 免费观看精品视频网站| a级毛片在线看网站| 男女午夜视频在线观看| 啦啦啦观看免费观看视频高清 | 老熟妇仑乱视频hdxx| 欧美国产日韩亚洲一区| 欧美日韩亚洲综合一区二区三区_| 精品国产超薄肉色丝袜足j| 人妻久久中文字幕网| 女性生殖器流出的白浆| 亚洲人成网站在线播放欧美日韩| 亚洲第一欧美日韩一区二区三区| 欧美日本中文国产一区发布| 欧美精品亚洲一区二区| 最新在线观看一区二区三区| 99在线人妻在线中文字幕| 变态另类成人亚洲欧美熟女 | 在线免费观看的www视频| 色婷婷久久久亚洲欧美| 精品一区二区三区四区五区乱码| 午夜福利,免费看| 手机成人av网站| e午夜精品久久久久久久| 99在线人妻在线中文字幕| 变态另类成人亚洲欧美熟女 | 日本精品一区二区三区蜜桃| 女性被躁到高潮视频| 精品一区二区三区四区五区乱码| 深夜精品福利| 一边摸一边抽搐一进一出视频| 亚洲欧美激情综合另类| 老汉色av国产亚洲站长工具| 久久国产精品影院| 亚洲第一欧美日韩一区二区三区| 亚洲成人国产一区在线观看| 久久国产亚洲av麻豆专区| 午夜亚洲福利在线播放|