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

    用于建筑繞流預測的非線性渦粘模型改良

    2014-06-15 17:15:12邵建濤趙加寧
    哈爾濱工業(yè)大學學報 2014年4期
    關(guān)鍵詞:預測出尾跡粘性

    邵建濤,劉 京,趙加寧,李 彪

    用于建筑繞流預測的非線性渦粘模型改良

    邵建濤1,2,劉 京1,3,趙加寧1,李 彪1

    (1.哈爾濱工業(yè)大學市政環(huán)境工程學院,150090哈爾濱;2.華東建筑設(shè)計研究院有限公司,200002上海;3.哈爾濱工業(yè)大學城市水資源與水環(huán)境國家重點實驗室,150090哈爾濱)

    為了改良非線性渦粘性模型模擬建筑繞流問題的表現(xiàn),首先介紹了采用渦粘性模型預測建筑繞流問題的現(xiàn)狀,分析了模擬中存在的問題,然后基于Craft非線性渦粘性模型提出了一種用于預測建筑繞流的改良的非線性渦粘性模型,并利用日本建筑學會提供的風洞實驗數(shù)據(jù)對改良的非線性模型進行分析驗證.結(jié)果表明:改良的非線性渦粘性模型一方面改善了標準k-ε模型建筑前端湍動動能預測過大的問題,預測出了建筑頂部的分離和再附著;另一方面通過增大尾跡區(qū)的湍動動能,改善了渦粘性模型在預測建筑尾跡區(qū)流動中的表現(xiàn).通過改良,非線性渦粘性模型可以較好地預測建筑風環(huán)境.

    風環(huán)境;建筑繞流;非線性渦粘性模型;計算流體力學;改良

    據(jù)聯(lián)合國人口調(diào)查局統(tǒng)計,2000年時,全球47%的人口居住在城市;預計到2015年,全球城市人口比例將會達52%,且大部分城市化發(fā)生在發(fā)展中國家[1].伴隨我國經(jīng)濟的發(fā)展,城市化進程的加快,城市人口占總?cè)丝诘谋戎卦絹碓酱螅鞘忻娣e越來越大.城市的熱環(huán)境和風環(huán)境直接影響著城市居民的健康.為此,人們投入了大量的精力研究城市氣候.隨著計算機計算能力的提高和CFD模擬軟件的發(fā)展,在設(shè)計階段,CFD技術(shù)越來越多地用來預測建筑風環(huán)境.建筑風環(huán)境模擬的核心是建筑繞流問題的模擬.多年來,許多學者對建筑繞流問題作了大量的研究[2-6];然而,建筑繞流問題的模擬長期以來一直是計算風工程領(lǐng)域的難點和熱點.對模擬建筑繞流的難點進行了分析,認為流動的高雷諾數(shù)、建筑前端的撞擊、鈍體尖銳的邊緣和出入口邊界條件是建筑繞流計算中的難點.由于對建筑前端湍流動能預測過大,標準kε模型被認為不適合用來模擬建筑繞流.文獻[7]通過研究發(fā)現(xiàn)修正的k-ε模型(如LK、MMK模型)能夠很好地修正標準k-ε模型在建筑前端高估湍動動能的問題.然而,LK模型和MMK模型存在對建筑再附著長度預測過長的缺點.

    由于非線性模擬能夠克服常規(guī)湍流模型的湍流各向同性的缺點,近年來,有研究者嘗試采用非線性湍流模式來模擬建筑繞流問題.Wright等嘗試用非線性湍流模型模擬建筑繞流問題,只有一個非線性的雷諾應(yīng)力模擬得到收斂解,其他模型很難得到收斂解[8].邵建濤等嘗試利用非線性kε模型對建筑繞流進行了模擬[9],發(fā)現(xiàn)非線性kε模型對建筑繞流問題不能得到穩(wěn)態(tài)收斂解;而利用非穩(wěn)態(tài)求解能夠得到建筑尾跡區(qū)橫向速度的周期性波動.盡管,利用Craft等提出的非線性kε模型采用URANS方法,可以得到較標準k-ε模型、LK模型、MMK模型更為準確的風環(huán)境預測.但是,URANS方法計算時間較長,很難在存在較大計算區(qū)域的現(xiàn)實工程中廣泛應(yīng)用.本文將嘗試通過風洞實驗數(shù)據(jù)對現(xiàn)有的非線性k-ε模型進行校準和改良,以使其能夠更為快速和準確的預測建筑風環(huán)境.

    1 CFD控制方程和湍流模型

    1.1 控制方程

    不可壓縮流體流動的控制方程為

    其中:ui為xi方向的瞬時速度,m/s;xi為空間坐標,m;i,j=1,2,3表示空間坐標系的3個方向;p為瞬時壓力,Pa;ν為粘性系數(shù),m2/s;ρ為流體密度,kg/m3.

    在雷諾平均模型中,瞬時速度ui可以雷諾分解成平均速度和脈動速度兩部分,即

    其中:Ui為平均速度,m/s;ui′為脈動速度,m/s.

    把式(3)分別代入式(1)、(2)后,可以得到

    1.2 非線性渦粘性模型

    為了更好地在渦粘性模型內(nèi)描述雷諾應(yīng)力,非線性渦粘性模型通常將雷諾應(yīng)力表示為平均變形率張量和平均渦量的高階形式.下式是一個典型的三階非線性渦粘性模型的雷諾應(yīng)力表達式.等號的右端前兩項為線性模型的表達式,等號右端前5項為二階非線性模型的表達式.

    式中:νt采用標準k-ε模型同樣的計算方法,即

    其中:k為湍動動能,m/s;ε為湍動動能耗散率,m2/s3.

    非線性k-ε模型中湍動動能和湍動動能耗散率的控制方程及其經(jīng)驗系數(shù)的取值和標準k-ε模型相應(yīng)的控制方程及其經(jīng)驗系數(shù)的取值相同.Sij和Ωij分別為平均變形率張量和平均渦量,1/s,可由下式計算得到.

    非線性渦粘性湍流模型大都十分類似,如Suga模型和Craft模型[10],Shih模型[11].各種非線性渦粘性模型的主要區(qū)別在于式(6)中C1到C7的非線性項系數(shù)選擇和式(7)中的cμ不同.

    1.3 非線性渦粘性模型的改進

    建筑繞流計算的難點在于正確預測建筑周圍的湍動動能分布.由Craft和Shih非線性渦粘性湍流模型的表達式可知,這兩個模型都對標準k-ε模型中的cμ系數(shù)進行了修正,修正了撞擊區(qū)湍動動能預測過大的問題.但根據(jù)前期研究結(jié)果可知,這兩種非線性渦粘性湍流模型對cμ的修正的同時會導致計算得到的尾跡區(qū)湍流粘性過低,穩(wěn)態(tài)計算很難達到收斂[9].本研究的主要工作是對cμ進行進一步適當?shù)母倪M,以達到改進撞擊區(qū)湍動動能預測過大的問題,同時增大尾跡區(qū)的湍動動能,一方面補償渦脫落造成的湍動增強,另一方面加強計算的穩(wěn)定性.通過對實驗數(shù)據(jù)的分析和大量的數(shù)值試驗,在Craft模型的基礎(chǔ)上,本文提出cμ表達式為

    本文改良的非線性渦粘性模型計算雷諾應(yīng)力的經(jīng)驗系數(shù)和Craft非線性模型相同.

    本模型中C1~C7分別為-0.1、-0.1、

    2 算例分析

    2.1 算例概況

    本文將采用文獻中最為常見的日本建筑學會標準算例中的1∶1∶2型建筑為例對本文改良的非線性渦粘性模型進行驗證,該建筑如圖1所示[12];建筑物的長度和寬度b均為0.08 m,高度H為0.16 m.以建筑高度H為特征長度和建筑高度處風速UH為特征速度的雷諾數(shù)約為24 000.

    為了避免計算區(qū)域邊界對模擬結(jié)果產(chǎn)生影響,本算例所采用的計算區(qū)域大小為21b×11b× 11.25b,如圖1所示.由于一階迎風差分的數(shù)值粘性較大,容易造成較大的計算誤差,計算中所有方程中的對流項均采用二階迎風差分.壓力速度耦合采用SIMPLE算法.計算采用立方體網(wǎng)格,如圖2所示.網(wǎng)格的最小間距為0.05H,出現(xiàn)在建筑壁面處,總網(wǎng)格數(shù)約為30萬.計算采用實驗數(shù)據(jù)和平衡湍流假設(shè)設(shè)置入口邊界;其中入口邊界的速度和湍動動能分布,按照風洞實驗數(shù)據(jù)擬合曲線給定;湍動動能耗散率邊界采用假設(shè)平衡態(tài)湍流.出口邊界采用壓力出口,即各物理量沿x方向的梯度為0.建筑表面和地面采用壁面函數(shù),其中建筑表面按光滑表面處理,地面的粗糙長度z0約為0.005.上表面和側(cè)面采用對稱邊界條件,各物理量在邊界的法向方向均為0,在切線方向梯度為0.

    圖1 計算區(qū)域

    圖2 計算網(wǎng)格

    為了使實驗結(jié)果更加可靠,本研究對網(wǎng)格無關(guān)性進行了檢驗.

    2.2 計算結(jié)果

    為了便于討論,本文對標準k-ε模型、Craft非線性渦粘模型、本文改良模型3種模型的計算結(jié)果和實驗結(jié)果進行對比討論.

    本文首先對預測結(jié)果中的速度場進行討論.圖3為建筑繞流再附著長度示意圖.表1為不同湍流模型的再附著長度模擬結(jié)果.圖4為y/b=0截面的流線圖.從圖中可以看出,除標準k-ε模型外,其他模型由于改善了撞擊區(qū)湍動動能預測過大的問題,都可不同程度地預測出來建筑頂部的渦.和其他研究一樣,標準k-ε模型沒能預測出建筑頂部氣流剝離和再貼附形成的渦旋.

    圖3 再附著長度示意圖

    表1 不同湍流模型的再附著長度模擬結(jié)果

    圖4 y/b=0截面的流線圖

    對建筑尾跡區(qū)而言,標準k-ε模型和Craft模型預測的建筑尾跡區(qū)再附著長度過長.本文改良模型預測的尾跡區(qū)在附著長度和風洞實驗值較為接近.

    圖5為y/b=0截面上多條豎線上的速度分布.在建筑前,x/b=-0.75處,標準k-ε模型、Craft模型和本文改良模型都能較好地預測建筑前的速度分布,本文改進模型和風洞實驗數(shù)據(jù)更為接近.

    圖5 y/b=0截面上多條豎線上的速度分布

    而在建筑頂部x/b=0,由于標準k-ε模型沒有預測出建筑頂部的渦,所以標準k-ε模型沒有預測到建筑頂部的回流.而Craft模型和本文改進模型都成功預測出建筑頂部的建筑回流.同時,標準k-ε模型和Craft模型都低估了z=0.2 m處附近的風速,而本文改良模型成功捕捉到了該處的最大風速.在建筑尾跡區(qū),x/b=1.25處,由于Craft模型預測建筑回流速度過大,標準k-ε模型和本文改良模型預測該處的風速和風洞實驗值較為接近.在再附著點附近,x/b=2處,標準k-ε模型和Craft模型都還能預測出地面附近的回流,這和風洞實驗的結(jié)果不符;而本文改良模型預測結(jié)果和風洞實驗結(jié)果非常接近,沒有出現(xiàn)地面附近的回流,成功預測除了流動進入到邊界層的再發(fā)展區(qū).

    圖6 y/b=0截面的湍動動能分布云圖(m2/s2)

    圖6為y/b=0截面的湍動動能分布云圖.從圖中可以看出,如文獻中指出,標準k-ε模型在建筑前端預測湍動動能過大,在建筑尾跡區(qū)預測湍動過小.通過對計算湍動應(yīng)力計算方法的改進,考慮到湍動的各項異性,Craft模型在一定程度上修正了在標準k-ε模型在建筑前端湍動動能預測過大的影響;然而這種修正也同時造成了尾跡區(qū)湍動動能預測更小.這也導致Craft模型成功預測了建筑頂部的渦,而惡化了建筑尾跡區(qū)渦大小的預測.本文改良模型由于在Craft模型的基礎(chǔ)上改進了湍動應(yīng)力計算方法,同時成功地修正了建筑前端湍動動能預測過大的問題和尾跡區(qū)湍動動能預測過小的問題,提高了湍動動能的預測精度.

    圖7為y/b=0截面上多條豎線上的湍動動能分布.從圖中可以看出,在建筑前端,x/b= -0.75處,標準k-ε模型預測建筑前端湍動動能最大;Craft模型改善了標準k-ε模型的預測結(jié)果,但預測的湍動動能較實驗值仍然較大;本文改良模型預測的建筑前端湍動動能值和風洞實驗值最為接近,有較好的預測精度.在建筑頂部,x/b=-0.75處,3種湍流模型都不能很好的預測湍動動能的分布.從圖中可以看出,在建筑屋面附近,標準k-ε模型能夠預測湍動動能的最大值,而Craft模型和本文改良模型都不能預測出該處的湍動動能最大值;在z=0.2 m處附近,本文改良模型能夠很好的預測湍動動能大小,而標準k -ε模型和Craft模型都低估了該處湍動動能的大小.在建筑尾跡區(qū),x/b=1.25處,Craft模型嚴重低估了該處的湍動動能大小,而標準k-ε模型在一定程度上低估了該處的湍動動能大小,本文改良模型稍微高估了該處的湍動動能大小.在x/b=2處,標準k-ε模型和Craft模型都低估了該處的湍動動能大小,而本文改良模型由于采用了新的湍動應(yīng)力的計算方法,在該處成功準確的預測了該處的湍動動能大小.

    圖7 y/b=0截面上多條豎線上的湍動動能分布

    2.3 結(jié)果分析

    從以上的結(jié)果可以看出,由于本文改良的非線性模型考慮了湍流的各項異性,同時采用了新的湍流應(yīng)力計算方法,本文改良模型一方面修正了標準k-ε模型建筑前端湍動動能預測過大的問題,同時又增大了建筑尾跡區(qū)的湍動動能.因此,本文改良模型一方面預測出了建筑頂部漩渦,另一方面較為準確地預測了建筑尾跡區(qū)的再附著長度.

    從建筑繞流的流動特性而言,標準k-ε模型對建筑前端的湍動動能預測過大是由于標準k-ε模型對湍流應(yīng)力的預測具有各項同性特點.很多對標準k-ε模型的修正都可以解決這個問題.但是,大部分修正的k-ε模型都會惡化對建筑尾跡區(qū)的預測,使預測再附著長度更長.事實上,建筑尾跡區(qū)存在著劇烈的非穩(wěn)態(tài)流動,大量的渦脫落使得該部分混合較為劇烈.大部分渦粘性不能預測出該部分的渦脫落,因此對再附著長度的預測都顯得過長.本文改良模型通過增大尾跡區(qū)的湍動動能,加強湍動擴散以求達到和渦脫落相似的混合效果,取得較為準確的尾跡區(qū)的流場.

    3 結(jié) 論

    1)建筑繞流問題是建筑環(huán)境模擬中的重要問題.基于Craft模型提出了一種用于建筑繞流預測的改良的非線性渦粘性模型,并利用日本建筑學會提供的實驗數(shù)據(jù)對其進行了驗證.

    2)改良了模型,一方面修正了標準k-ε模型建筑前端湍動動能預測過大的問題,預測出了建筑頂部的分離和再附著;另一方面通過增大尾跡區(qū)的湍動動能,改善了渦粘性模型在預測建筑尾跡區(qū)流動中的表現(xiàn).改良模型預測的建筑尾跡區(qū)再附著長度和風洞實驗數(shù)據(jù)非常接近.

    [1]United Nations,Department of Economic and Social Affairs.World urbanization prospects:the 2009 revision[EB/OL].[2009-03-01].http://esa.un.org/unpd/wup/CD-ROM-2009/WUP2009-F03-Urban-Population. xls.

    [2]STATHOPOULOST,BASKARANB.Computer simulation of wind environmental conditions around buildings[J].Engineering Structures,1996,18(11):876-885.

    [3]TSUCHIYA M,MURAKAMI S,MOCHIDA A,et al. Development of a new k-ε model for flow and pressure fields aroundbluffbody[J].JournalofWind Engineering and Industrial Aerodynamics,1997,67-68:169-182.

    [4]MURAKAMI S,OOKA R,MOCHIDA A,et al.CFD analysis of wind climate from human scale to urban scale[J].Journal of Wind Engineering and Industrial Aerodynamics,1999,81(1/2/3):57-81.

    [5]WANG X,MCNAMARA K F.Evaluation of CFD simulation using RANS turbulence models for building effects on pollution dispersion[J].Environmental Fluid Mechanics,2006,6(2):181-202.

    [6]村上周三.CFD與建筑環(huán)境[M].朱清宇,譯.北京:中國建筑工業(yè)出版社,2007.

    [7]MURAKAMI S.Current status and future trends in computational wind engineering[J].Journal of Wind Engineering and Industrial Aerodynamics,1997,67-68:3-34.

    [8]WRIGHT N G,EASOM G J.Non-linear turbulence model results for flow over a building at full-scale[J]. Applied Mathematical Modelling,2003,27(12): 1013-1033.

    [9]SHAO Jiantao,LIU Jing,ZHAO Jianing.Evaluation of various non-linear k-ε models for predicting wind flow around an isolated high-rise building within the surface boundary layer[J].Building and Environment,2012,57(11):145-155.

    [10]CRAFT T J,LAUNDER B E,SUGA K.Development and application of a cubic eddy-viscosity model of turbulence[J].International Journal of Heat and Fluid Flow,1996,17(2):108-115.

    [11]SHIH T H,ZHU J,LUMLEY J L.A realizable reynolds stress algebraicequationmodel[R].Cleveland: NASA,1993.

    [12]TOMINAGA Y,MOCHIDA A,MURAKAMI S,et al. Comparison of various revised k-ε models and LES applied to flow around a high-rise building model with 1∶1∶2 shape placed within the surface boundary layer[J].Journal of Wind Engineering and Industrial Aerodynamics,2008,96(4):389-411.

    (編輯 魏希柱)

    Improvement of the non-linear eddy viscosity model applied to predicting wind flow around building

    SHAO Jiantao1,2,LIU Jing1,3,ZHAO Jianing1,LI Biao1
    (1.School of Municipal and Environmental Engineering,Harbin Institute of Technology,150090 Harbin,China;2.East China Architectural Design&Research Institute Co.Ltd.,200002 Shanghai,China;3.State Key Laboratory of Urban Water Resource and Environment,Harbin Institute of Technology,150090 Harbin,China)

    The aim of this paper is to improve the performance of the non-linear eddy viscosity model for simulating the wind flow around the building.Firstly,the state of art of the predicting wind flow around building using RANS model was introduced,and the problems in the simulation were analyzed.Then an improved non-linear eddy viscosity for predicting the wind flow around buildings was proposed based on Craft model.The improved non-linear eddy viscosity was validated and analyzed through the wind tunnel data provided by AIJ.The results showed that the proposed non-linear eddy viscosity improved the overestimation of turbulent kinetic energy in impingement region by the standard k-ε model,and predicted better results in the wake region behind buildings simultaneously through strengthening the eddy viscosity in the wake region.After the improvement,the non-linear eddy viscosity model can predict the wind environment around buildings better.

    wind environment;flow around building;non-linear eddy viscosity model;computational fluid dynamics;improvement

    TU111.19+3

    A

    0367-6234(2014)04-0050-07

    2013-03-30.

    國家自然科學基金資助項目(40505025);

    城市氣象科學研究基金資助項目(UMRF201004).

    邵建濤(1983—),男,博士,高級工程師;

    劉 京(1972—),男,教授,博士生導師;

    趙加寧(1956—),女,教授,博士生導師.

    邵建濤,shaojiantao@gmail.com.

    猜你喜歡
    預測出尾跡粘性
    一類具有粘性項的擬線性拋物型方程組
    一種基于Radon 變換和尾跡模型的尾跡檢測算法
    帶粘性的波動方程組解的逐點估計
    基于EEMD-Hilbert譜的渦街流量計尾跡振蕩特性
    吸引力無法則
    粘性非等熵流體方程平衡解的穩(wěn)定性
    量體薦衣
    家庭醫(yī)生增強基層首診粘性
    基于FABEMD和Goldstein濾波器的SAR艦船尾跡圖像增強方法
    雷達學報(2012年4期)2012-10-03 12:25:02
    閱讀理解沙場點兵
    国产精品一区二区三区四区免费观看| 91久久精品电影网| 色噜噜av男人的天堂激情| 久久久精品大字幕| 精品人妻偷拍中文字幕| 在线天堂最新版资源| 国产极品精品免费视频能看的| 国产精品久久久久久久电影| 中文字幕人妻熟人妻熟丝袜美| 国产亚洲av嫩草精品影院| 国产一区二区三区在线臀色熟女| 国产av不卡久久| 亚洲美女视频黄频| av国产免费在线观看| 国产av在哪里看| 成人毛片60女人毛片免费| 在线观看免费视频日本深夜| 女人被狂操c到高潮| 久久久色成人| 在线观看午夜福利视频| 乱人视频在线观看| 成年女人看的毛片在线观看| 热99在线观看视频| 久久精品久久久久久噜噜老黄 | 免费电影在线观看免费观看| 女人被狂操c到高潮| 有码 亚洲区| 久久久久久久久大av| 99久久精品热视频| 免费av毛片视频| 久久精品国产清高在天天线| 高清日韩中文字幕在线| 久久精品人妻少妇| 不卡一级毛片| 三级男女做爰猛烈吃奶摸视频| 免费av不卡在线播放| 搞女人的毛片| 熟妇人妻久久中文字幕3abv| 国产成人一区二区在线| 日韩在线高清观看一区二区三区| 在线国产一区二区在线| 欧美+日韩+精品| 久久久久久久久中文| 国产在线男女| 老女人水多毛片| 直男gayav资源| 国产精品永久免费网站| 国产亚洲av片在线观看秒播厂 | 中文资源天堂在线| 18+在线观看网站| 日韩欧美在线乱码| 91狼人影院| 欧美最黄视频在线播放免费| 黄色欧美视频在线观看| 欧美潮喷喷水| 欧美变态另类bdsm刘玥| 联通29元200g的流量卡| 天堂影院成人在线观看| 久久久a久久爽久久v久久| 亚洲av成人精品一区久久| 中文字幕人妻熟人妻熟丝袜美| 又黄又爽又刺激的免费视频.| 级片在线观看| 亚洲精品乱码久久久v下载方式| 国产久久久一区二区三区| 亚洲欧美精品综合久久99| 老司机福利观看| 又爽又黄a免费视频| 九九爱精品视频在线观看| 日本在线视频免费播放| 日本在线视频免费播放| 高清在线视频一区二区三区 | 亚洲欧洲日产国产| 日韩强制内射视频| 国产一区二区三区在线臀色熟女| 噜噜噜噜噜久久久久久91| 淫秽高清视频在线观看| 亚洲精品国产成人久久av| 国产精品三级大全| 淫秽高清视频在线观看| 日韩强制内射视频| 乱人视频在线观看| 国产精品人妻久久久影院| 日韩欧美精品v在线| 自拍偷自拍亚洲精品老妇| 国产精品一区二区性色av| 精品久久久久久久人妻蜜臀av| 欧美成人精品欧美一级黄| 国产视频内射| 啦啦啦啦在线视频资源| 夜夜夜夜夜久久久久| 久久草成人影院| 一级av片app| 亚洲av不卡在线观看| 午夜精品国产一区二区电影 | 能在线免费看毛片的网站| 男女视频在线观看网站免费| 99国产精品一区二区蜜桃av| 欧美xxxx黑人xx丫x性爽| 最后的刺客免费高清国语| 亚洲四区av| 国模一区二区三区四区视频| 18禁裸乳无遮挡免费网站照片| 97热精品久久久久久| 久久久久九九精品影院| 国产淫片久久久久久久久| 青春草国产在线视频 | 国产在线精品亚洲第一网站| 久久这里有精品视频免费| 少妇的逼水好多| 欧美激情国产日韩精品一区| 白带黄色成豆腐渣| 亚洲中文字幕一区二区三区有码在线看| 寂寞人妻少妇视频99o| 99在线视频只有这里精品首页| 一本久久中文字幕| 男女视频在线观看网站免费| 18禁裸乳无遮挡免费网站照片| 午夜免费激情av| 人妻少妇偷人精品九色| 国产一级毛片在线| 直男gayav资源| 十八禁国产超污无遮挡网站| 人人妻人人看人人澡| 成年版毛片免费区| 国产极品精品免费视频能看的| av天堂在线播放| 国产一区二区在线观看日韩| 国产精品免费一区二区三区在线| 亚洲国产欧美在线一区| 国产精品三级大全| 日本av手机在线免费观看| 亚洲高清免费不卡视频| 最近手机中文字幕大全| 美女 人体艺术 gogo| 桃色一区二区三区在线观看| 中文资源天堂在线| 国产日韩欧美在线精品| av专区在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 国产高清不卡午夜福利| 日韩av不卡免费在线播放| 久久热精品热| 哪个播放器可以免费观看大片| 九九久久精品国产亚洲av麻豆| 亚洲国产欧美在线一区| 日日撸夜夜添| 免费观看a级毛片全部| 国产一区二区在线av高清观看| 一夜夜www| 男人的好看免费观看在线视频| 麻豆乱淫一区二区| 成人亚洲精品av一区二区| 欧美最黄视频在线播放免费| 日韩在线高清观看一区二区三区| 网址你懂的国产日韩在线| 中文亚洲av片在线观看爽| 久久精品91蜜桃| 国内少妇人妻偷人精品xxx网站| 深夜a级毛片| 九九久久精品国产亚洲av麻豆| 免费不卡的大黄色大毛片视频在线观看 | 亚洲欧美成人精品一区二区| 欧美日韩国产亚洲二区| 久久久久久久久中文| 天天一区二区日本电影三级| 国产亚洲精品久久久久久毛片| av在线播放精品| av视频在线观看入口| 国产淫片久久久久久久久| 啦啦啦啦在线视频资源| 一夜夜www| 久久久色成人| 免费看av在线观看网站| 国产亚洲av片在线观看秒播厂 | 人妻夜夜爽99麻豆av| 亚洲av中文av极速乱| 草草在线视频免费看| 99热只有精品国产| 成人漫画全彩无遮挡| 嫩草影院入口| 国产美女午夜福利| 中出人妻视频一区二区| 国产成人精品一,二区 | 国产男人的电影天堂91| 激情 狠狠 欧美| 成人无遮挡网站| 99在线视频只有这里精品首页| 免费在线观看成人毛片| 91精品国产九色| 成人一区二区视频在线观看| 久久亚洲国产成人精品v| 99视频精品全部免费 在线| 日韩av不卡免费在线播放| 免费观看a级毛片全部| 国产精品国产三级国产av玫瑰| 久久久久久久午夜电影| 久久午夜福利片| 国产中年淑女户外野战色| 国产精品麻豆人妻色哟哟久久 | 乱人视频在线观看| 久久久久九九精品影院| 久久99热6这里只有精品| 小说图片视频综合网站| 人妻夜夜爽99麻豆av| av免费观看日本| 亚洲av成人av| 天天躁日日操中文字幕| 久久精品久久久久久久性| 成年免费大片在线观看| 日本一本二区三区精品| 欧美日韩国产亚洲二区| 国产亚洲精品久久久久久毛片| 日本色播在线视频| 97在线视频观看| 欧美xxxx黑人xx丫x性爽| 99在线视频只有这里精品首页| 国产视频内射| 精品久久久久久久久久免费视频| 国产伦一二天堂av在线观看| 久久热精品热| 国产精品一区二区三区四区久久| 亚洲高清免费不卡视频| 51国产日韩欧美| 日本一本二区三区精品| 啦啦啦观看免费观看视频高清| 国产免费男女视频| 日产精品乱码卡一卡2卡三| 国产av不卡久久| 国产精品一区www在线观看| 女人被狂操c到高潮| 床上黄色一级片| 99国产精品一区二区蜜桃av| 性欧美人与动物交配| 精品久久久久久成人av| 午夜视频国产福利| 伦理电影大哥的女人| 麻豆乱淫一区二区| 国产又黄又爽又无遮挡在线| av福利片在线观看| 久久精品久久久久久噜噜老黄 | 18禁黄网站禁片免费观看直播| 国产精品久久电影中文字幕| 色综合色国产| 成人国产麻豆网| 美女内射精品一级片tv| 国产乱人偷精品视频| 国产精品一区二区性色av| 欧美一区二区国产精品久久精品| a级毛色黄片| 国产午夜精品论理片| 亚洲欧美精品专区久久| 国产成人福利小说| 亚洲va在线va天堂va国产| 成人亚洲精品av一区二区| 美女被艹到高潮喷水动态| 久久精品91蜜桃| 亚洲欧美成人综合另类久久久 | 成人鲁丝片一二三区免费| 成年版毛片免费区| 国产白丝娇喘喷水9色精品| 欧美丝袜亚洲另类| 韩国av在线不卡| 国产精华一区二区三区| 久久精品久久久久久久性| 老司机福利观看| 国产精品久久久久久久久免| 亚洲欧美精品综合久久99| 久久久a久久爽久久v久久| 国产探花在线观看一区二区| 欧美潮喷喷水| 日韩,欧美,国产一区二区三区 | 亚洲自偷自拍三级| 少妇裸体淫交视频免费看高清| 在线观看美女被高潮喷水网站| 综合色丁香网| 日本欧美国产在线视频| 麻豆精品久久久久久蜜桃| 日本黄色片子视频| 亚洲国产精品久久男人天堂| 久久久精品大字幕| 青春草视频在线免费观看| 18+在线观看网站| 精品不卡国产一区二区三区| 国产爱豆传媒在线观看| 边亲边吃奶的免费视频| 一区二区三区高清视频在线| 久久精品国产自在天天线| 中文精品一卡2卡3卡4更新| 九色成人免费人妻av| 国产精品爽爽va在线观看网站| av卡一久久| 欧美成人精品欧美一级黄| 少妇裸体淫交视频免费看高清| 一本久久中文字幕| 桃色一区二区三区在线观看| 真实男女啪啪啪动态图| a级毛片a级免费在线| 国产成人a∨麻豆精品| 日本三级黄在线观看| 国产亚洲精品久久久久久毛片| 丰满的人妻完整版| 变态另类丝袜制服| 人妻夜夜爽99麻豆av| 在线免费观看的www视频| 丝袜美腿在线中文| 麻豆国产av国片精品| 欧美极品一区二区三区四区| 亚洲av成人精品一区久久| 亚洲在久久综合| 亚洲国产高清在线一区二区三| 91av网一区二区| 日韩强制内射视频| av视频在线观看入口| 少妇熟女欧美另类| 亚洲国产精品久久男人天堂| 日韩大尺度精品在线看网址| 少妇的逼水好多| 亚洲电影在线观看av| 免费看a级黄色片| 久久精品国产亚洲av香蕉五月| 噜噜噜噜噜久久久久久91| 美女内射精品一级片tv| 好男人在线观看高清免费视频| 免费黄网站久久成人精品| 国产亚洲精品av在线| 高清午夜精品一区二区三区 | 99久久无色码亚洲精品果冻| 亚洲精品影视一区二区三区av| 国产中年淑女户外野战色| 成年免费大片在线观看| 国产伦精品一区二区三区四那| 国产精品久久久久久久电影| 亚洲国产欧洲综合997久久,| av卡一久久| 日韩精品青青久久久久久| 亚洲精品日韩av片在线观看| 久久久成人免费电影| 日本爱情动作片www.在线观看| 精品一区二区三区视频在线| 又爽又黄无遮挡网站| 99精品在免费线老司机午夜| 亚洲最大成人中文| 变态另类成人亚洲欧美熟女| 亚洲精品久久久久久婷婷小说 | 国产伦精品一区二区三区四那| av福利片在线观看| 最好的美女福利视频网| 97在线视频观看| 欧美成人a在线观看| 秋霞在线观看毛片| 国产激情偷乱视频一区二区| 小说图片视频综合网站| 日日啪夜夜撸| 麻豆国产97在线/欧美| 能在线免费看毛片的网站| 亚洲人与动物交配视频| 丝袜喷水一区| 看十八女毛片水多多多| 久久精品人妻少妇| 男女下面进入的视频免费午夜| 两个人视频免费观看高清| 国产片特级美女逼逼视频| 国产91av在线免费观看| 国产精品乱码一区二三区的特点| 男的添女的下面高潮视频| 国内揄拍国产精品人妻在线| 免费看美女性在线毛片视频| 精品一区二区三区视频在线| 欧美高清性xxxxhd video| 网址你懂的国产日韩在线| 91aial.com中文字幕在线观看| 精品国产三级普通话版| av视频在线观看入口| 欧美丝袜亚洲另类| 十八禁国产超污无遮挡网站| 国产成人精品久久久久久| 99热网站在线观看| 搡老妇女老女人老熟妇| 又爽又黄无遮挡网站| 日韩在线高清观看一区二区三区| 精品熟女少妇av免费看| 国产精品免费一区二区三区在线| 久久6这里有精品| 日韩一区二区三区影片| 亚洲成av人片在线播放无| av在线蜜桃| 日韩欧美一区二区三区在线观看| 嫩草影院入口| 女的被弄到高潮叫床怎么办| 久久这里有精品视频免费| 国产爱豆传媒在线观看| 国产成人精品婷婷| 亚洲婷婷狠狠爱综合网| 欧美高清性xxxxhd video| 国内精品一区二区在线观看| 2022亚洲国产成人精品| 欧美一区二区亚洲| 国产视频内射| 国产激情偷乱视频一区二区| 最好的美女福利视频网| 麻豆国产97在线/欧美| 在线免费十八禁| 国产精品电影一区二区三区| 亚洲成人av在线免费| 国产精品av视频在线免费观看| 亚洲综合色惰| 99久久无色码亚洲精品果冻| 亚洲av成人精品一区久久| av免费观看日本| 亚洲五月天丁香| 国产又黄又爽又无遮挡在线| 精品不卡国产一区二区三区| 国产精品av视频在线免费观看| 毛片女人毛片| 中文在线观看免费www的网站| 美女脱内裤让男人舔精品视频 | a级毛片免费高清观看在线播放| 欧美色欧美亚洲另类二区| 国产精品人妻久久久久久| 久久精品人妻少妇| 精品免费久久久久久久清纯| 欧美色欧美亚洲另类二区| 亚洲人成网站在线观看播放| 欧美+亚洲+日韩+国产| 亚洲国产欧美人成| 国产美女午夜福利| 国产亚洲精品av在线| 91午夜精品亚洲一区二区三区| 舔av片在线| 成人鲁丝片一二三区免费| 尾随美女入室| 在线观看一区二区三区| 最近中文字幕高清免费大全6| 又爽又黄无遮挡网站| 国产熟女欧美一区二区| 亚洲av免费在线观看| 国内精品宾馆在线| 精品久久久久久久人妻蜜臀av| 悠悠久久av| 国产黄色小视频在线观看| 国产v大片淫在线免费观看| 搞女人的毛片| 99热只有精品国产| 欧美成人一区二区免费高清观看| 欧美潮喷喷水| 亚洲av中文av极速乱| 丰满乱子伦码专区| 亚洲人成网站在线播| 欧美激情久久久久久爽电影| 欧美bdsm另类| 午夜福利视频1000在线观看| 国产大屁股一区二区在线视频| 中文字幕av成人在线电影| 高清毛片免费看| 亚洲最大成人中文| 亚洲av成人av| 久久午夜亚洲精品久久| av专区在线播放| 久久人人精品亚洲av| 国产一区二区三区在线臀色熟女| 美女高潮的动态| 少妇高潮的动态图| 婷婷六月久久综合丁香| av女优亚洲男人天堂| 欧美高清成人免费视频www| 欧美+亚洲+日韩+国产| 卡戴珊不雅视频在线播放| 国产久久久一区二区三区| 老熟妇乱子伦视频在线观看| 国产精品一区二区三区四区久久| 日韩一区二区三区影片| 婷婷精品国产亚洲av| 国产成人精品婷婷| 国产黄色视频一区二区在线观看 | 亚洲精品乱码久久久v下载方式| 99视频精品全部免费 在线| 欧美精品国产亚洲| 成熟少妇高潮喷水视频| .国产精品久久| 哪里可以看免费的av片| 啦啦啦啦在线视频资源| 国产精品乱码一区二三区的特点| 精品午夜福利在线看| 国产日本99.免费观看| 成人高潮视频无遮挡免费网站| 波多野结衣高清作品| 乱系列少妇在线播放| www.av在线官网国产| 丰满人妻一区二区三区视频av| 久久久久网色| 亚洲在久久综合| 色吧在线观看| 久久人人爽人人片av| 欧美区成人在线视频| 国产黄a三级三级三级人| 亚洲精品乱码久久久v下载方式| 九九在线视频观看精品| 青青草视频在线视频观看| 国产av一区在线观看免费| av在线蜜桃| 性插视频无遮挡在线免费观看| 99久久成人亚洲精品观看| 亚洲av中文av极速乱| 国模一区二区三区四区视频| 国产伦精品一区二区三区四那| 成人国产麻豆网| 午夜激情福利司机影院| 成年版毛片免费区| 久久综合国产亚洲精品| 久久精品国产鲁丝片午夜精品| 日本与韩国留学比较| 久久九九热精品免费| 亚洲欧美成人精品一区二区| 99热只有精品国产| 夜夜夜夜夜久久久久| 人人妻人人澡人人爽人人夜夜 | 男女边吃奶边做爰视频| 夜夜夜夜夜久久久久| 悠悠久久av| 狠狠狠狠99中文字幕| 日韩大尺度精品在线看网址| 午夜a级毛片| 日韩亚洲欧美综合| h日本视频在线播放| 少妇的逼好多水| 国产av一区在线观看免费| 波多野结衣高清作品| 精品日产1卡2卡| 国产爱豆传媒在线观看| 国产成人影院久久av| 看黄色毛片网站| av又黄又爽大尺度在线免费看 | 丰满的人妻完整版| 51国产日韩欧美| 变态另类成人亚洲欧美熟女| 91久久精品国产一区二区三区| 91精品国产九色| 蜜臀久久99精品久久宅男| 毛片一级片免费看久久久久| 亚洲人成网站在线观看播放| 免费不卡的大黄色大毛片视频在线观看 | or卡值多少钱| av在线老鸭窝| 18禁黄网站禁片免费观看直播| 亚洲最大成人手机在线| 国产精品国产高清国产av| 亚洲成人中文字幕在线播放| 欧美xxxx黑人xx丫x性爽| 久久久久国产网址| 97在线视频观看| 啦啦啦韩国在线观看视频| 在线观看一区二区三区| 嘟嘟电影网在线观看| 亚洲一区二区三区色噜噜| 日韩高清综合在线| 少妇熟女aⅴ在线视频| 中文字幕熟女人妻在线| 日本一本二区三区精品| 美女国产视频在线观看| 一进一出抽搐动态| 身体一侧抽搐| 午夜福利在线观看吧| 欧美激情国产日韩精品一区| 亚洲精品国产成人久久av| 国产在视频线在精品| eeuss影院久久| 午夜激情欧美在线| 日韩欧美 国产精品| 高清毛片免费观看视频网站| 两性午夜刺激爽爽歪歪视频在线观看| 最近视频中文字幕2019在线8| 国产亚洲av片在线观看秒播厂 | 精品午夜福利在线看| 国产精品久久电影中文字幕| 少妇被粗大猛烈的视频| 狂野欧美白嫩少妇大欣赏| 99久久久亚洲精品蜜臀av| 能在线免费看毛片的网站| 国产精品一区二区三区四区免费观看| 深爱激情五月婷婷| 免费观看精品视频网站| 五月玫瑰六月丁香| 亚洲国产精品久久男人天堂| 久久精品久久久久久噜噜老黄 | 好男人视频免费观看在线| 精品少妇黑人巨大在线播放 | 丰满人妻一区二区三区视频av| 晚上一个人看的免费电影| 最好的美女福利视频网| 精品无人区乱码1区二区| 我的老师免费观看完整版| 欧美色视频一区免费| 美女脱内裤让男人舔精品视频 | 国产视频内射| 欧美一区二区精品小视频在线| 国产人妻一区二区三区在| 成人特级黄色片久久久久久久| 久久国内精品自在自线图片| 边亲边吃奶的免费视频| 午夜激情福利司机影院| 久久精品人妻少妇| 99热全是精品| 久久久午夜欧美精品| 久久久久性生活片| 三级国产精品欧美在线观看| 欧美色视频一区免费| 免费观看在线日韩| 99在线视频只有这里精品首页| 夜夜看夜夜爽夜夜摸| 在现免费观看毛片| 精品不卡国产一区二区三区| 久久久久九九精品影院| 一个人看的www免费观看视频| 一级毛片我不卡|