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

    植被對波浪作用下床面切應(yīng)力影響的數(shù)值模擬分析

    2022-11-14 06:33:58葛昭佩唐軍趙楚嫣
    海洋學(xué)報 2022年11期
    關(guān)鍵詞:波流床面邊界層

    葛昭佩,唐軍*,趙楚嫣

    ( 1. 大連理工大學(xué) 海岸和近海工程國家重點實驗室,遼寧 大連 116023)

    1 引言

    波浪作用下的床面切應(yīng)力對海床泥沙起動具有關(guān)鍵性作用,是模擬分析近岸泥沙輸移、岸灘演變的重要參數(shù)。近岸植被具有防浪固灘護岸的作用,研究近岸植被對波浪作用下床面切應(yīng)力的影響有助于理解其防護海岸的作用機制。

    目前關(guān)于波浪作用下床面切應(yīng)力的研究取得了一定成果。Jonsson[1]從理論公式出發(fā),推導(dǎo)了層流邊界層的床面切應(yīng)力及邊界層厚度公式,并引入波浪底摩阻系數(shù)及邊界層流態(tài)判別方法,給出了紊流邊界層的波浪摩阻系數(shù)經(jīng)驗公式;孔令雙等[2]基于邊界層理論,導(dǎo)出了波流環(huán)境中床面切應(yīng)力,并建立了泥沙起動及水體挾沙力公式;蔡翠蘇[3]基于水槽試驗,探討了規(guī)則波和不規(guī)則波作用下底摩阻系數(shù)計算方法;齊富康等[4]基于渤海海峽觀測數(shù)據(jù),分別應(yīng)用Soulsby和Van Rijn 的模型計算波流環(huán)境中床面切應(yīng)力;Lin和Zhang[5]建立三維數(shù)值模型,模擬計算線性波、Stokes 波、孤立波等條件下層流邊界層床面切應(yīng)力,并分別與解析解比較分析;滕涌等[6]基于ECOMSED波流耦合底邊界層模型,分析了波高、近底流速及水深對波流條件下床面切應(yīng)力的影響;Larsen 和Fuhrman[7–8]基于waves2Foam 模擬海嘯波的傳播及爬升過程,著重討論了海嘯波作用下的邊界層及床面切應(yīng)力并提出了一種預(yù)測海嘯波作用下瞬時邊界層厚度和床面切應(yīng)力的方法。相對裸床,植被影響下波浪水動力特性復(fù)雜,目前關(guān)于植被水域波浪作用下床面切應(yīng)力的計算以經(jīng)驗公式為主。Wang 等[9]基于水槽試驗分析了水流和波流條件下植被水域平均流速及湍流動能,并由湍流動能推算植被水域床面切應(yīng)力,發(fā)現(xiàn)純水流時植被水域床面切應(yīng)力沿程降低,而波流條件下植被水域床面切應(yīng)力增大并在離開植被區(qū)后顯著大幅降低;Reidenbach 和Thomas[10]觀測了弗吉尼亞海岸保護區(qū)水動力,并應(yīng)用經(jīng)驗公式分別計算波、流的床面切應(yīng)力,發(fā)現(xiàn)與純水流時相比波流條件下床面切應(yīng)力大幅增大,位于大葉藻水域測點處的床面切應(yīng)力始終小于裸床處且小于泥沙起動臨界切應(yīng)力,有效保護了海床免受侵蝕;陳家貴和沈小雄[11]基于Jonsson[1]提出的層流邊界層公式,結(jié)合水槽試驗數(shù)據(jù),分析了入射波高及柔性植被對最大床面切應(yīng)力的影響;李勰等[12]基于Jonsson[1]提出的邊界層最大切應(yīng)力公式及Luhar 等[13]提出的植被帶中邊界層流速公式,分析了剛性植被根莖對床面切應(yīng)力沿程變化的影響。

    總體來看,目前關(guān)于植被對波浪或波流作用下床面切應(yīng)力影響的研究主要基于水槽試驗或現(xiàn)場觀測數(shù)據(jù)并結(jié)合經(jīng)驗公式分析計算,但由于植被影響下的流速剖面及雷諾應(yīng)力分布復(fù)雜,經(jīng)驗公式未能很好地描述波浪條件下植被水域床面切應(yīng)力[9,14–15]。本文以O(shè)penFOAM 中的waves2Foam 求解器[16]為基礎(chǔ),在其中引入植被源項,建立植被水域波流數(shù)值水槽,通過計算全水深速度場及邊界層內(nèi)速度梯度而直接給出床面切應(yīng)力分布。在應(yīng)用已有實驗數(shù)據(jù)及理論公式對數(shù)值模型驗證的基礎(chǔ)上,著重分析了入射波高、植被密度、植被淹沒高度及水流對植被水域波浪作用下床面切應(yīng)力的影響特性。

    2 數(shù)值模型

    waves2Foam 求解器是由Jacobsen 等[16]基于Open-FOAM 開發(fā)的一個用于模擬波浪的第三方模塊庫,提供了水流、線性波、Stokes 波、孤立波等造波文件,亦可自定義造波文件實現(xiàn)其他類型波浪模擬。

    2.1 控制方程

    waves2Foam 求解器控制方程基于OpenFOAM 標準,采用氣液兩相流方程。為考慮波浪與植被相互作用,在動量方程和湍流模型中引入植被源項,即拖曳阻力項和慣性力項[17](稱之為宏觀數(shù)值模型,以下未特別說明的均指該數(shù)值模型)??刂品匠虨?/p>

    傳統(tǒng)湍流模型應(yīng)用于氣液兩相流時,由于氣液界面處較大的密度梯度[19]及湍流動能的過度評估[20]會造成大波陡波傳播過程中的波高非物理衰減。為抑制此種現(xiàn)象,本文采用Larsen 和Fuhrman[20]提出的修正k-ε湍流模型。

    表1 經(jīng)驗系數(shù)取值Table 1 Default values for the closure coefficient

    waves2Foam 為兩相流求解器,其定義體積分數(shù)α ( 0 ≤α ≤1)來描述每個網(wǎng)格中各相的體積比重。各相的運動滿足的控制方程為

    式(8)第三項是用于保持氣液界面清晰的人工壓縮項;ur,i為 壓縮速度;壓縮系數(shù)Cα∈[0,1],其默認值為1。

    2.2 植被拖曳力系數(shù)

    植被拖曳力系數(shù)CD是對植被水域水動力特性進行宏觀數(shù)值模擬時的關(guān)鍵參數(shù),與波浪、水流、植被密切相關(guān),直接影響數(shù)值模擬的準確性,但其取值尚無普適方法。本文將上述控制方程中的植被源項去除(令CD、CM、Ckp、Cεp4 項經(jīng)驗系數(shù)均取為0),通過直接考慮剛性植被外形對波流的擾動,詳細模擬純波及波流條件下植被水域的水動力特征(稱之為精細數(shù)值模型,在本文中僅用于計算各工況代表CD值),通過提取植被上的波流力及其所在截面的沿植被高度平均流速并使用Morison 公式(式(10))直接計算植被拖曳力系數(shù)。

    式中,F(xiàn)D為 植株上的拖曳力;FI為植株上的慣性力;F為植株上的波流力,Dalrymple 等[22]指出阻力與流速之間無相位差,阻力的做功可由阻力和流速的乘積在周期內(nèi)積分獲得,然而慣性力與流速之間存在 π/2的相位差,使得慣性力在周期內(nèi)的做功為0,故上式計算過程中可忽略慣性力[18],F(xiàn)根據(jù)植被上壓力及附近流場由式(11)計算;u為植株所在截面的平均流速;ρ為水的密度;hv為 植被淹沒高度;bv為植被直徑。

    2.3 邊界條件

    數(shù)值水槽入口與出口采用松弛區(qū)造波與消波,入口處速度場、壓力場和自由波面由波浪理論公式給出,出口采用自由出流條件;床面采用無滑移邊界條件;兩側(cè)壁采用周期性邊界條件,該邊界條件能夠消除兩側(cè)壁面的影響,可將水槽簡化為較窄的數(shù)值水槽以提高計算效率,且已被成功應(yīng)用于波浪的數(shù)值模擬[23–24];頂端采用壓力出口邊界。時間步長由庫朗數(shù)自動調(diào)節(jié),為保證計算穩(wěn)定,最大庫朗數(shù)設(shè)置為0.25。

    3 模型驗證

    由于缺乏波浪作用下植被水域床面切應(yīng)力實測數(shù)據(jù),本文將分別驗證植被水域波面演化和無植被時波浪作用下床面切應(yīng)力。采用海岸和近海工程國家重點實驗室開展的植被影響下波浪傳播物理水槽試驗結(jié)果[25]驗證本文模型模擬植被水域波浪傳播的有效性;采用理論公式驗證層流邊界層時的床面切應(yīng)力,采用徐華等[26]的試驗驗證紊流邊界層時的床面切應(yīng)力。

    3.1 植被水域波浪傳播驗證

    選取表2 中的兩種試驗工況驗證模型模擬植被水域波浪傳播的有效性,其中CD值由精細數(shù)值模型計算,CD值計算過程以表2 中工況1 為例,數(shù)值模擬中設(shè)置數(shù)值水槽長為18 m,高為0.6 m,寬為0.12 m,取植被水域長為1.08 m,其中以圓柱代表剛性植被,整體網(wǎng)格尺寸 Δx=Δy=0 .03 m, Δz=0.01 m,為更好地貼合植被外形,分別加密近植被域及單株植被附近網(wǎng)格。邊界條件同2.3 節(jié),此外,植被表面設(shè)為無滑移邊界條件。

    表2 植被及波浪參數(shù)Table 2 Parameters of vegetation and waves

    圖1 為精細模擬中工況1 的波面演化和x=0.6 m(x=0 為植被域起點)處植被上的波浪力及截面平均速度,使用式(10)計算該處周期平均CD值,并利用同樣的方法計算每一排植株處的周期平均CD值后,做平均得出空間–周期平均CD值作為該工況的代表CD值。

    圖1 考慮植被外形擾動的模擬結(jié)果Fig. 1 Simulation results considering the disturbance of vegetation shape

    將利用精細模型計算的各工況代表CD值代入宏觀模型中模擬計算。波浪經(jīng)過植被水域的波面演化及流速衰減的模擬值與試驗值對比如圖2 和圖3 所示,由圖可知,模型計算結(jié)果和試驗數(shù)據(jù)[25]符合性較好,表明利用精細數(shù)值模型計算CD值準確并且使用宏觀數(shù)值模型能夠準確模擬植被水域波浪及波流水動力變化。

    圖2 植被水域波面演化Fig. 2 Free surface evolution along the vegetation zones

    圖3 植被水域流速衰減驗證Fig. 3 Verification of velocity attenuation in vegetation zones

    3.2 波浪作用下床面切應(yīng)力驗證

    本文數(shù)值模型通過求解雷諾平均Navier-Stokes 方程得出波浪作用下全水深速度場及邊界層內(nèi)速度梯度,并由公式(12)計算床面切應(yīng)力。

    為驗證模型計算床面切應(yīng)力的準確性,此處設(shè)置3 組驗證工況,其中工況1、工況2 為本試驗工況,工況3 為徐華等[26]的試驗工況,波浪參數(shù)如表3 所示。工況1 和工況2 的邊界層雷諾數(shù)均小于1.26×104,依據(jù)Jonsson[1]提出的判定標準,兩者均屬于層流邊界層,工況3 則屬于光滑紊流邊界層。

    表3 驗證工況參數(shù)Table 3 Parameters of verification conditions

    邊界層流速剖面及床面切應(yīng)力驗證結(jié)果如圖4至圖6 所示。其中,工況1 采用線性波理論解驗證,工況2 采用五階Stokes 波理論解驗證,工況3 采用試驗測量值[26]驗證。從圖中可以看出,待波浪場穩(wěn)定后,工況1 模擬值與理論值符合性較好;工況2 床面切應(yīng)力在波谷時模擬值較理論值偏大,但差值小于10%;工況3 模擬值與試驗測量值在一個完整周期內(nèi)總體符合良好,表明該數(shù)值模型能夠準確計算波浪作用下床面切應(yīng)力。

    圖4 理論值與模擬值對比(工況1)Fig. 4 Comparison between theoretical and simulated values (case 1)

    圖5 理論值與模擬值對比(工況2)Fig. 5 Comparison between theoretical and simulated values (case 2)

    圖6 床面切應(yīng)力對比(工況3)Fig. 6 Comparison of bed shear stress (case 3)

    4 模擬結(jié)果分析

    為分析不同入射波高、植被密度、植被淹沒高度及水流流速條件下,植被水域床面切應(yīng)力的分布特征,本文共模擬計算22 種工況,見表4。其中工況1、2、3、12、14、15、18、19、20 為文獻[25]的試驗工況,這9 種工況的波面演化及流速衰減模擬值與試驗值均符合良好。

    表4 模擬工況參數(shù)Table 4 Parameters of numerical simulation conditions

    純波時,波浪經(jīng)過植被水域后波高衰減的同時床面切應(yīng)力也隨之減小,見圖7(x=0 m 處為植被起點),這是由于植被的阻水作用,近底流速沿植被水域衰減使得邊界層流速梯度減小所致(圖8)。圖7 中植被前后的波面與床面切應(yīng)力均存在 π/4的相位差,說明植被對波高和床面切應(yīng)力的衰減作用是同步的。此外,該相位差與式(13)、式(14)相一致,進一步驗證了該數(shù)值模型的準確性。

    圖7 波面及床面切應(yīng)力(工況1)Fig. 7 Free surface and bed shear stress (case 1)

    圖8 近底流速剖面(工況1)Fig. 8 The near-bottom velocity profile (case 1)

    當植被處于淹沒狀態(tài)時,植被的阻水效應(yīng)僅發(fā)生在冠層以下,使這一區(qū)域流速減小且阻水作用隨著淹沒高度的增大而增大,而冠層以上由于過水斷面的減小使得流速增大且超過無植被時的流速,故在植被頂端形成較大的速度梯度(圖9);此外,由于邊界層的作用在近底處存在一速度峰值。上述水平速度沿水深的分布規(guī)律與Liu 等[29]的試驗結(jié)果一致。

    圖9 淹沒植被沿水深速度剖面Fig. 9 Longitudinal velocity profile of submerged vegetation

    同向水流的存在會增加波浪水質(zhì)點速度的非線性,同時增大正向床面切應(yīng)力幅值,減小負向床面切應(yīng)力幅值,使得一個波周期內(nèi)正向切應(yīng)力時間延長,負向切應(yīng)力時間縮短,如圖10 所示。

    圖10 不同流速下床面切應(yīng)力變化Fig. 10 Variation of bed shear stress under different current velocity

    5 結(jié)論

    本文OpenFOAM 中的waves2Foam 求解器建立了含植被水域的波流數(shù)值水槽,模擬分析了波浪、水流、植被要素對植被水域床面切應(yīng)力衰減的影響。模擬結(jié)果表明:當波浪傳播到植被前端時,床面切應(yīng)力會出現(xiàn)增大現(xiàn)象,并隨著入射波高增大、水流流速增大、植被密度增大此種現(xiàn)象愈加明顯;純波時,波高和植被密度的增大均會導(dǎo)致植被水域床面切應(yīng)力衰減幅度的增大,并且床面切應(yīng)力在植被水域前段的衰減幅度較大,后段衰減幅度減??;對于完全淹沒植向床面切應(yīng)力幅值減小,使得一個波周期內(nèi)正向切應(yīng)力時間延長負向切應(yīng)力時間縮短;當水流流速較小時,床面切應(yīng)力的衰減率與無水流時相近,之后隨著流速增大衰減率隨之增大。

    圖11 不同工況下植被水域最大床面切應(yīng)力分布Fig. 11 Distribution of maximum bed shear stress in vegetation zones under different conditions

    圖12 不同流速下最大床面切應(yīng)力衰減率Fig. 12 Decay rate of maximum bed shear stress at different current velocities

    猜你喜歡
    波流床面邊界層
    魚鱗狀床面粗糙特性
    波流耦合下樁周珊瑚砂沖刷機理研究
    海洋通報(2022年5期)2022-11-30 12:04:02
    對瓦里安碳纖維治療床面模型的評估
    波流聯(lián)合作用下海上輸油漂浮軟管動力響應(yīng)分析
    淹沒植物明渠床面沖淤及其對水流運動的影響
    基于HIFiRE-2超燃發(fā)動機內(nèi)流道的激波邊界層干擾分析
    改進的投影覆蓋方法對遼河河道粗糙床面分維量化研究
    偶感
    槽道內(nèi)渦波流場展向渦的分布特征
    一類具有邊界層性質(zhì)的二次奇攝動邊值問題
    99国产综合亚洲精品| 欧美日韩中文字幕国产精品一区二区三区 | 人人妻人人添人人爽欧美一区卜| 成人手机av| 久久人妻熟女aⅴ| 亚洲精品粉嫩美女一区| 日韩欧美一区视频在线观看| 日韩欧美国产一区二区入口| 热re99久久精品国产66热6| 国产老妇伦熟女老妇高清| 午夜福利影视在线免费观看| 国产精品一区二区在线不卡| 一级片'在线观看视频| 中文精品一卡2卡3卡4更新| 国产麻豆69| 国产黄色免费在线视频| cao死你这个sao货| 久久人妻福利社区极品人妻图片| 色94色欧美一区二区| 亚洲精品美女久久av网站| 午夜福利乱码中文字幕| 日韩精品免费视频一区二区三区| av天堂在线播放| 超碰成人久久| 两个人看的免费小视频| 黑人猛操日本美女一级片| 亚洲欧美成人综合另类久久久| 亚洲人成电影免费在线| 香蕉国产在线看| 男女之事视频高清在线观看| 免费在线观看完整版高清| 美女中出高潮动态图| 搡老岳熟女国产| 亚洲伊人色综图| 美女视频免费永久观看网站| 久久国产精品影院| 亚洲精品国产av成人精品| 久久青草综合色| 看免费av毛片| 人人妻人人添人人爽欧美一区卜| 久久久精品94久久精品| 亚洲免费av在线视频| 成人国语在线视频| 亚洲av男天堂| 极品人妻少妇av视频| 国产成人av教育| 亚洲人成电影免费在线| av福利片在线| 波多野结衣一区麻豆| 免费人妻精品一区二区三区视频| 在线观看免费日韩欧美大片| 免费看十八禁软件| 亚洲av美国av| 日韩视频在线欧美| 91九色精品人成在线观看| 高潮久久久久久久久久久不卡| 人人澡人人妻人| 国产成人av激情在线播放| 亚洲国产av新网站| 亚洲av男天堂| 爱豆传媒免费全集在线观看| 婷婷丁香在线五月| 午夜激情久久久久久久| 99久久精品国产亚洲精品| 菩萨蛮人人尽说江南好唐韦庄| 国产一区二区在线观看av| 国产成人免费观看mmmm| 国产97色在线日韩免费| 色婷婷av一区二区三区视频| 国产精品亚洲av一区麻豆| 国产极品粉嫩免费观看在线| 欧美亚洲 丝袜 人妻 在线| 性色av乱码一区二区三区2| 19禁男女啪啪无遮挡网站| 性高湖久久久久久久久免费观看| 日韩,欧美,国产一区二区三区| 午夜福利,免费看| 一区二区日韩欧美中文字幕| 巨乳人妻的诱惑在线观看| 国产亚洲精品第一综合不卡| 精品亚洲乱码少妇综合久久| 自线自在国产av| 亚洲 国产 在线| 成年av动漫网址| 精品国产乱码久久久久久小说| 99国产精品免费福利视频| 19禁男女啪啪无遮挡网站| 捣出白浆h1v1| 天天添夜夜摸| 亚洲一码二码三码区别大吗| 黄片大片在线免费观看| 波多野结衣一区麻豆| 成人手机av| 麻豆av在线久日| 欧美成狂野欧美在线观看| 青草久久国产| 18禁黄网站禁片午夜丰满| 亚洲人成电影免费在线| 精品国产一区二区久久| 一区在线观看完整版| 在线观看免费午夜福利视频| 日韩大片免费观看网站| e午夜精品久久久久久久| 狂野欧美激情性xxxx| 一本—道久久a久久精品蜜桃钙片| 国产精品亚洲av一区麻豆| av网站在线播放免费| 法律面前人人平等表现在哪些方面 | 亚洲天堂av无毛| 亚洲欧美成人综合另类久久久| 操美女的视频在线观看| 免费久久久久久久精品成人欧美视频| 久久国产精品人妻蜜桃| 国产成人a∨麻豆精品| 欧美乱码精品一区二区三区| 久久精品亚洲熟妇少妇任你| 亚洲欧美精品综合一区二区三区| 少妇 在线观看| 每晚都被弄得嗷嗷叫到高潮| 无遮挡黄片免费观看| 亚洲 国产 在线| 飞空精品影院首页| 国产亚洲av高清不卡| 一边摸一边做爽爽视频免费| 欧美成人午夜精品| 国产亚洲一区二区精品| 欧美日韩亚洲国产一区二区在线观看 | 日本av免费视频播放| 色精品久久人妻99蜜桃| 90打野战视频偷拍视频| 亚洲精品一二三| 日韩电影二区| 一区二区日韩欧美中文字幕| 国产日韩欧美在线精品| www.自偷自拍.com| av在线老鸭窝| 日韩电影二区| tube8黄色片| 韩国精品一区二区三区| 午夜视频精品福利| 婷婷色av中文字幕| 嫩草影视91久久| 狠狠精品人妻久久久久久综合| 每晚都被弄得嗷嗷叫到高潮| 国产精品 欧美亚洲| 亚洲五月婷婷丁香| 视频区欧美日本亚洲| 亚洲全国av大片| 不卡一级毛片| 啦啦啦啦在线视频资源| 蜜桃国产av成人99| 久久人人97超碰香蕉20202| 午夜免费成人在线视频| 亚洲少妇的诱惑av| 午夜福利视频在线观看免费| 精品国产一区二区三区四区第35| 精品亚洲成a人片在线观看| videos熟女内射| 国产精品麻豆人妻色哟哟久久| 国产精品香港三级国产av潘金莲| 男人添女人高潮全过程视频| 男女床上黄色一级片免费看| 啦啦啦 在线观看视频| 国产免费现黄频在线看| 黄色视频,在线免费观看| 欧美日韩亚洲综合一区二区三区_| 亚洲国产日韩一区二区| 国产精品.久久久| 亚洲欧美一区二区三区黑人| 久久精品成人免费网站| 欧美成狂野欧美在线观看| 午夜精品久久久久久毛片777| av福利片在线| 别揉我奶头~嗯~啊~动态视频 | 亚洲欧美清纯卡通| 亚洲自偷自拍图片 自拍| 王馨瑶露胸无遮挡在线观看| 精品第一国产精品| 国产精品 欧美亚洲| 久久久久久久精品精品| 下体分泌物呈黄色| 欧美日韩黄片免| 国产一区有黄有色的免费视频| 亚洲国产日韩一区二区| 国产av精品麻豆| 久久香蕉激情| 热99re8久久精品国产| 亚洲成国产人片在线观看| 亚洲第一欧美日韩一区二区三区 | 久久人人爽av亚洲精品天堂| 国产一区二区 视频在线| 精品少妇内射三级| 亚洲国产精品成人久久小说| 韩国高清视频一区二区三区| 999精品在线视频| 国产欧美日韩一区二区三区在线| 日韩制服丝袜自拍偷拍| 新久久久久国产一级毛片| 成年动漫av网址| 亚洲第一青青草原| 蜜桃在线观看..| 国产有黄有色有爽视频| 视频在线观看一区二区三区| 久久香蕉激情| 老司机午夜福利在线观看视频 | 亚洲第一av免费看| 国产野战对白在线观看| 极品少妇高潮喷水抽搐| 国产精品麻豆人妻色哟哟久久| 日韩欧美一区二区三区在线观看 | 欧美 日韩 精品 国产| av天堂在线播放| 精品一品国产午夜福利视频| 99国产精品免费福利视频| 亚洲av日韩在线播放| cao死你这个sao货| 桃花免费在线播放| 欧美精品av麻豆av| 午夜福利视频精品| 不卡av一区二区三区| 亚洲三区欧美一区| 在线观看人妻少妇| 在线永久观看黄色视频| 国产av一区二区精品久久| 欧美激情极品国产一区二区三区| 欧美激情高清一区二区三区| 国产欧美日韩一区二区精品| 精品久久蜜臀av无| 亚洲精品久久午夜乱码| 亚洲色图综合在线观看| 高潮久久久久久久久久久不卡| 9191精品国产免费久久| 久久久精品国产亚洲av高清涩受| 中文字幕色久视频| 日日夜夜操网爽| 一级a爱视频在线免费观看| 啦啦啦免费观看视频1| 一本一本久久a久久精品综合妖精| 欧美日韩亚洲综合一区二区三区_| 一区二区三区精品91| 欧美日韩一级在线毛片| 久久久久国产一级毛片高清牌| 亚洲av片天天在线观看| 97人妻天天添夜夜摸| 午夜视频精品福利| 久久毛片免费看一区二区三区| 久久久国产精品麻豆| 男女之事视频高清在线观看| 日本a在线网址| 日本av免费视频播放| 日本撒尿小便嘘嘘汇集6| 搡老乐熟女国产| 咕卡用的链子| 欧美国产精品va在线观看不卡| 超碰97精品在线观看| 日韩欧美一区视频在线观看| 一本一本久久a久久精品综合妖精| 日韩熟女老妇一区二区性免费视频| 国产亚洲精品久久久久5区| 免费高清在线观看视频在线观看| tocl精华| 日韩一卡2卡3卡4卡2021年| 正在播放国产对白刺激| 欧美人与性动交α欧美软件| 国产欧美日韩一区二区三区在线| 男女无遮挡免费网站观看| 91麻豆精品激情在线观看国产 | 美女国产高潮福利片在线看| 久久人人97超碰香蕉20202| 精品人妻在线不人妻| 叶爱在线成人免费视频播放| 一区二区三区精品91| 人妻久久中文字幕网| 久久国产亚洲av麻豆专区| 巨乳人妻的诱惑在线观看| 亚洲全国av大片| 欧美变态另类bdsm刘玥| 欧美日本中文国产一区发布| 大陆偷拍与自拍| 亚洲美女黄色视频免费看| 欧美另类一区| 亚洲色图综合在线观看| 97人妻天天添夜夜摸| 久久久国产欧美日韩av| 中文字幕制服av| 成年动漫av网址| 黄色视频不卡| 人人妻人人爽人人添夜夜欢视频| 免费高清在线观看日韩| 国产精品偷伦视频观看了| 少妇人妻久久综合中文| 久久久精品国产亚洲av高清涩受| www.精华液| 精品卡一卡二卡四卡免费| 成人av一区二区三区在线看 | 啪啪无遮挡十八禁网站| 亚洲天堂av无毛| 一本色道久久久久久精品综合| 久久人人爽av亚洲精品天堂| 一边摸一边抽搐一进一出视频| 亚洲精品一区蜜桃| 国产av一区二区精品久久| 欧美xxⅹ黑人| 三级毛片av免费| 亚洲欧美一区二区三区久久| 亚洲国产毛片av蜜桃av| 蜜桃国产av成人99| 丰满饥渴人妻一区二区三| 日韩精品免费视频一区二区三区| 黑丝袜美女国产一区| av福利片在线| 亚洲国产中文字幕在线视频| 国精品久久久久久国模美| 亚洲熟女毛片儿| 巨乳人妻的诱惑在线观看| 99re6热这里在线精品视频| 麻豆av在线久日| 午夜激情久久久久久久| 99国产精品99久久久久| 老司机靠b影院| 国产亚洲精品久久久久5区| 久热爱精品视频在线9| 国产xxxxx性猛交| www.av在线官网国产| 大香蕉久久成人网| 欧美中文综合在线视频| 欧美 亚洲 国产 日韩一| 亚洲伊人色综图| 国产一区二区三区综合在线观看| 香蕉丝袜av| 亚洲欧美成人综合另类久久久| 91精品国产国语对白视频| 婷婷色av中文字幕| 国产精品久久久av美女十八| 最近中文字幕2019免费版| 国产av又大| 日日摸夜夜添夜夜添小说| 超碰97精品在线观看| 18禁国产床啪视频网站| 欧美精品一区二区大全| 操美女的视频在线观看| 中文字幕av电影在线播放| 亚洲熟女毛片儿| 黄片小视频在线播放| 母亲3免费完整高清在线观看| 99热国产这里只有精品6| 亚洲精品久久午夜乱码| 国产无遮挡羞羞视频在线观看| 日本黄色日本黄色录像| 91国产中文字幕| 啦啦啦视频在线资源免费观看| 97精品久久久久久久久久精品| 久久久欧美国产精品| 母亲3免费完整高清在线观看| 国产亚洲午夜精品一区二区久久| 岛国在线观看网站| 亚洲精品国产av蜜桃| kizo精华| 欧美日韩亚洲综合一区二区三区_| 日本vs欧美在线观看视频| 免费一级毛片在线播放高清视频 | 久久精品aⅴ一区二区三区四区| 国产伦人伦偷精品视频| 亚洲免费av在线视频| 国产精品一区二区在线不卡| 亚洲自偷自拍图片 自拍| 精品亚洲成a人片在线观看| 黄色视频,在线免费观看| 99热网站在线观看| 50天的宝宝边吃奶边哭怎么回事| 国产老妇伦熟女老妇高清| 日本欧美视频一区| 母亲3免费完整高清在线观看| 一级片'在线观看视频| 丰满人妻熟妇乱又伦精品不卡| 欧美成人午夜精品| 午夜视频精品福利| 亚洲九九香蕉| 桃花免费在线播放| 久久 成人 亚洲| 国产精品影院久久| 亚洲中文av在线| 少妇被粗大的猛进出69影院| 午夜免费观看性视频| av网站免费在线观看视频| 精品欧美一区二区三区在线| 在线十欧美十亚洲十日本专区| 一个人免费看片子| 女人高潮潮喷娇喘18禁视频| 成人亚洲精品一区在线观看| 国产麻豆69| 王馨瑶露胸无遮挡在线观看| 老熟女久久久| 久久久国产精品麻豆| 亚洲免费av在线视频| 亚洲欧洲日产国产| 9191精品国产免费久久| 国产日韩欧美亚洲二区| 丝袜美腿诱惑在线| 亚洲欧美成人综合另类久久久| 久久国产亚洲av麻豆专区| 美女视频免费永久观看网站| 免费在线观看完整版高清| 午夜福利在线观看吧| 狂野欧美激情性xxxx| 热99国产精品久久久久久7| 涩涩av久久男人的天堂| 日本a在线网址| 亚洲avbb在线观看| 少妇 在线观看| 18禁国产床啪视频网站| 在线观看免费日韩欧美大片| 国产男女内射视频| 色婷婷久久久亚洲欧美| 天天影视国产精品| 欧美在线一区亚洲| svipshipincom国产片| a级毛片黄视频| 国产精品国产av在线观看| 亚洲性夜色夜夜综合| 欧美一级毛片孕妇| 欧美乱码精品一区二区三区| 丁香六月天网| 国产福利在线免费观看视频| svipshipincom国产片| 麻豆乱淫一区二区| 国产精品九九99| 亚洲国产欧美在线一区| 成年人午夜在线观看视频| 免费黄频网站在线观看国产| 永久免费av网站大全| 欧美精品av麻豆av| 国产成人av激情在线播放| 丰满饥渴人妻一区二区三| 国产精品二区激情视频| av福利片在线| 国产亚洲精品一区二区www | 欧美 亚洲 国产 日韩一| 91成人精品电影| 男女免费视频国产| 大码成人一级视频| 久热爱精品视频在线9| 老熟妇仑乱视频hdxx| 午夜福利在线免费观看网站| 人人妻人人澡人人看| 精品久久久精品久久久| 一区福利在线观看| 男人操女人黄网站| 亚洲精品国产av成人精品| 国产av一区二区精品久久| 亚洲天堂av无毛| 欧美日韩国产mv在线观看视频| 日韩电影二区| 久久精品成人免费网站| 精品一区在线观看国产| 韩国高清视频一区二区三区| 国产一区二区激情短视频 | 99国产精品一区二区三区| 黄片小视频在线播放| 一区二区av电影网| www.熟女人妻精品国产| 少妇粗大呻吟视频| 国内毛片毛片毛片毛片毛片| 亚洲精品乱久久久久久| 中文字幕高清在线视频| 蜜桃国产av成人99| av天堂在线播放| 欧美在线一区亚洲| 少妇的丰满在线观看| 91九色精品人成在线观看| 搡老岳熟女国产| 亚洲av美国av| 欧美日韩福利视频一区二区| 欧美 亚洲 国产 日韩一| a级片在线免费高清观看视频| 成人黄色视频免费在线看| 日日摸夜夜添夜夜添小说| 久久精品国产a三级三级三级| 欧美激情极品国产一区二区三区| 18禁裸乳无遮挡动漫免费视频| 建设人人有责人人尽责人人享有的| 久久精品国产亚洲av香蕉五月 | 美女扒开内裤让男人捅视频| 亚洲精品在线美女| 中文字幕另类日韩欧美亚洲嫩草| 欧美少妇被猛烈插入视频| 另类精品久久| 最近中文字幕2019免费版| www.熟女人妻精品国产| 久久久久久人人人人人| 亚洲欧美精品自产自拍| 天堂8中文在线网| 波多野结衣av一区二区av| 香蕉国产在线看| 亚洲色图综合在线观看| 国产精品麻豆人妻色哟哟久久| 最近最新免费中文字幕在线| 窝窝影院91人妻| 日韩视频在线欧美| 欧美日韩成人在线一区二区| 久久人人爽人人片av| 99国产综合亚洲精品| 国产一区二区 视频在线| 国产一区有黄有色的免费视频| 色精品久久人妻99蜜桃| 色播在线永久视频| 久久亚洲精品不卡| 两性夫妻黄色片| a 毛片基地| 国产免费现黄频在线看| 久久精品久久久久久噜噜老黄| 王馨瑶露胸无遮挡在线观看| 考比视频在线观看| 久久香蕉激情| 777久久人妻少妇嫩草av网站| 老司机影院毛片| 天堂8中文在线网| 亚洲国产看品久久| 久久精品久久久久久噜噜老黄| 老司机亚洲免费影院| 99re6热这里在线精品视频| 国产男女内射视频| 国产精品影院久久| 捣出白浆h1v1| 777米奇影视久久| 性色av一级| 久久女婷五月综合色啪小说| 久久青草综合色| 精品一品国产午夜福利视频| 亚洲色图 男人天堂 中文字幕| 欧美激情久久久久久爽电影 | 精品免费久久久久久久清纯 | www.自偷自拍.com| 久久天躁狠狠躁夜夜2o2o| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产成人一精品久久久| 免费女性裸体啪啪无遮挡网站| av天堂久久9| 啦啦啦免费观看视频1| 久久久国产成人免费| 美女国产高潮福利片在线看| 国产精品国产av在线观看| 日韩中文字幕视频在线看片| 18禁观看日本| 嫁个100分男人电影在线观看| 一本大道久久a久久精品| 国产精品一区二区在线观看99| 人妻久久中文字幕网| 成年av动漫网址| 亚洲精品自拍成人| 91精品伊人久久大香线蕉| 亚洲成av片中文字幕在线观看| 国产男女内射视频| 老汉色∧v一级毛片| 母亲3免费完整高清在线观看| 曰老女人黄片| 久久免费观看电影| 欧美精品啪啪一区二区三区 | 国产97色在线日韩免费| 午夜影院在线不卡| 免费久久久久久久精品成人欧美视频| 久久精品aⅴ一区二区三区四区| 男女国产视频网站| 国产不卡av网站在线观看| 欧美久久黑人一区二区| 一二三四社区在线视频社区8| 精品一区二区三区av网在线观看 | 精品亚洲乱码少妇综合久久| 最新在线观看一区二区三区| 欧美av亚洲av综合av国产av| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久久久精品古装| 国产福利在线免费观看视频| 欧美国产精品一级二级三级| 亚洲全国av大片| tocl精华| 久久国产精品影院| 国产精品久久久久久人妻精品电影 | av在线app专区| 后天国语完整版免费观看| 99国产精品一区二区三区| 亚洲精品一二三| 2018国产大陆天天弄谢| 午夜福利视频在线观看免费| 精品国产乱码久久久久久小说| 日韩一区二区三区影片| 亚洲精品中文字幕一二三四区 | 男女边摸边吃奶| 国产精品熟女久久久久浪| 成年美女黄网站色视频大全免费| av在线播放精品| 亚洲,欧美精品.| 久久这里只有精品19| 国产欧美日韩一区二区三区在线| 青草久久国产| 超碰成人久久| 亚洲性夜色夜夜综合| 精品第一国产精品| 国产亚洲精品第一综合不卡| 久久久国产成人免费| 水蜜桃什么品种好| 啦啦啦 在线观看视频| 亚洲人成电影免费在线| 一级黄色大片毛片| 久久99一区二区三区| 亚洲欧美清纯卡通| 国产在线视频一区二区| 大片电影免费在线观看免费| 国产精品久久久久久精品电影小说| 欧美日本中文国产一区发布| 色综合欧美亚洲国产小说| 久久久久国产精品人妻一区二区| 免费少妇av软件| 亚洲精品久久久久久婷婷小说| av国产精品久久久久影院|