• <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ì)的二次奇攝動邊值問題
    亚洲精品中文字幕在线视频| 男人添女人高潮全过程视频| av免费在线观看网站| 老司机影院毛片| 日韩精品免费视频一区二区三区| www.自偷自拍.com| 手机成人av网站| 欧美人与性动交α欧美软件| 丝袜美足系列| 一级毛片女人18水好多| 多毛熟女@视频| 男女高潮啪啪啪动态图| 亚洲欧洲日产国产| 可以免费在线观看a视频的电影网站| 国产一区二区在线观看av| 久久精品久久久久久噜噜老黄| 亚洲人成电影免费在线| 国产国语露脸激情在线看| 亚洲一码二码三码区别大吗| 中文字幕另类日韩欧美亚洲嫩草| 国产av精品麻豆| 亚洲av日韩精品久久久久久密| 日本黄色日本黄色录像| 免费在线观看完整版高清| 免费在线观看完整版高清| 午夜福利,免费看| h视频一区二区三区| 国产一区有黄有色的免费视频| 婷婷成人精品国产| 亚洲熟女毛片儿| 亚洲精品国产区一区二| 黑人巨大精品欧美一区二区mp4| 日韩 亚洲 欧美在线| 老司机在亚洲福利影院| 国产成+人综合+亚洲专区| 中文字幕精品免费在线观看视频| 中文字幕av电影在线播放| 亚洲欧洲精品一区二区精品久久久| 国产一区二区 视频在线| xxxhd国产人妻xxx| 精品少妇内射三级| 丝袜美足系列| 国产在视频线精品| 免费在线观看视频国产中文字幕亚洲 | 99久久人妻综合| 好男人电影高清在线观看| 999精品在线视频| 飞空精品影院首页| 建设人人有责人人尽责人人享有的| 亚洲欧美日韩另类电影网站| 韩国高清视频一区二区三区| 久久久久久免费高清国产稀缺| 老汉色∧v一级毛片| 国产一级毛片在线| 欧美日韩黄片免| 狠狠精品人妻久久久久久综合| 99久久人妻综合| 亚洲少妇的诱惑av| 欧美av亚洲av综合av国产av| 在线观看一区二区三区激情| 男人添女人高潮全过程视频| 人妻 亚洲 视频| 他把我摸到了高潮在线观看 | 欧美变态另类bdsm刘玥| www.av在线官网国产| 肉色欧美久久久久久久蜜桃| 亚洲av电影在线进入| 午夜福利在线免费观看网站| 岛国毛片在线播放| 色94色欧美一区二区| 国产精品一区二区在线观看99| 欧美xxⅹ黑人| a在线观看视频网站| 波多野结衣一区麻豆| 亚洲色图 男人天堂 中文字幕| 69精品国产乱码久久久| 国产一区二区在线观看av| 国产成人欧美在线观看 | 精品亚洲成a人片在线观看| 精品视频人人做人人爽| 另类亚洲欧美激情| 亚洲男人天堂网一区| 91av网站免费观看| 午夜福利在线观看吧| 日本av免费视频播放| 久久这里只有精品19| 2018国产大陆天天弄谢| 91大片在线观看| 高潮久久久久久久久久久不卡| 亚洲成人手机| 热99国产精品久久久久久7| 国产日韩欧美视频二区| 久久久久久人人人人人| www.自偷自拍.com| 亚洲第一av免费看| 国产又色又爽无遮挡免| 精品人妻熟女毛片av久久网站| 丝袜美腿诱惑在线| 超色免费av| 日韩制服骚丝袜av| √禁漫天堂资源中文www| 国产激情久久老熟女| 欧美国产精品一级二级三级| 日本av手机在线免费观看| 亚洲成人手机| 一边摸一边抽搐一进一出视频| 丝袜脚勾引网站| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品一区蜜桃| 日韩 欧美 亚洲 中文字幕| 国产主播在线观看一区二区| 亚洲精品美女久久久久99蜜臀| 欧美日本中文国产一区发布| 99国产综合亚洲精品| 男女之事视频高清在线观看| cao死你这个sao货| 免费人妻精品一区二区三区视频| 男女下面插进去视频免费观看| 99国产精品一区二区蜜桃av | 建设人人有责人人尽责人人享有的| 纯流量卡能插随身wifi吗| 极品少妇高潮喷水抽搐| 男女无遮挡免费网站观看| 深夜精品福利| 成在线人永久免费视频| 男女下面插进去视频免费观看| 亚洲av国产av综合av卡| 国产一区二区三区av在线| 女人精品久久久久毛片| 久久午夜综合久久蜜桃| 色播在线永久视频| 亚洲国产精品成人久久小说| 纯流量卡能插随身wifi吗| h视频一区二区三区| 岛国毛片在线播放| 久久女婷五月综合色啪小说| 免费观看av网站的网址| 欧美亚洲 丝袜 人妻 在线| 国产成人免费观看mmmm| 国产日韩一区二区三区精品不卡| 亚洲精品久久午夜乱码| 岛国在线观看网站| 欧美 日韩 精品 国产| 国产精品麻豆人妻色哟哟久久| av视频免费观看在线观看| av有码第一页| 美女大奶头黄色视频| 女警被强在线播放| www.av在线官网国产| 久久久水蜜桃国产精品网| 男女免费视频国产| 在线观看免费日韩欧美大片| 午夜老司机福利片| 久久久久国产一级毛片高清牌| 中国国产av一级| 男男h啪啪无遮挡| 欧美激情久久久久久爽电影 | 一进一出抽搐动态| 午夜成年电影在线免费观看| 高潮久久久久久久久久久不卡| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲欧美精品综合一区二区三区| 久久久久久免费高清国产稀缺| 97在线人人人人妻| 亚洲第一青青草原| 久久99一区二区三区| 欧美日韩成人在线一区二区| 在线观看免费午夜福利视频| 成人三级做爰电影| 青春草亚洲视频在线观看| 免费一级毛片在线播放高清视频 | 18禁国产床啪视频网站| 国产av精品麻豆| 成年人黄色毛片网站| 一本—道久久a久久精品蜜桃钙片| 9色porny在线观看| av不卡在线播放| 国产一级毛片在线| 久久亚洲精品不卡| 国产精品.久久久| 性色av乱码一区二区三区2| 欧美成狂野欧美在线观看| 成人国语在线视频| 中文字幕最新亚洲高清| 国产成+人综合+亚洲专区| 丰满饥渴人妻一区二区三| 国产精品麻豆人妻色哟哟久久| 捣出白浆h1v1| 欧美精品一区二区大全| tocl精华| 亚洲欧洲精品一区二区精品久久久| 搡老熟女国产l中国老女人| 婷婷成人精品国产| 欧美变态另类bdsm刘玥| 成人亚洲精品一区在线观看| 日韩视频一区二区在线观看| 免费观看人在逋| tocl精华| 国产精品秋霞免费鲁丝片| 免费久久久久久久精品成人欧美视频| 久久国产亚洲av麻豆专区| 九色亚洲精品在线播放| 韩国精品一区二区三区| 亚洲第一青青草原| 亚洲av日韩在线播放| 亚洲色图 男人天堂 中文字幕| 亚洲欧洲日产国产| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品成人av观看孕妇| 日本91视频免费播放| 天天添夜夜摸| 免费高清在线观看视频在线观看| 老司机午夜福利在线观看视频 | 亚洲一卡2卡3卡4卡5卡精品中文| 国产男女内射视频| 人妻人人澡人人爽人人| 国产一级毛片在线| 日韩欧美一区二区三区在线观看 | 少妇的丰满在线观看| 丁香六月欧美| 日韩中文字幕欧美一区二区| 成年动漫av网址| 日韩欧美免费精品| 热99re8久久精品国产| 国产亚洲午夜精品一区二区久久| 精品卡一卡二卡四卡免费| 亚洲熟女毛片儿| 国产日韩欧美在线精品| 国产男女超爽视频在线观看| 国产一区二区三区av在线| 人妻久久中文字幕网| 黄色毛片三级朝国网站| 99热全是精品| 精品一区二区三卡| 满18在线观看网站| 日韩一卡2卡3卡4卡2021年| 免费在线观看日本一区| 秋霞在线观看毛片| av欧美777| 天堂中文最新版在线下载| 午夜精品久久久久久毛片777| 一级毛片电影观看| 亚洲自偷自拍图片 自拍| 亚洲成人免费av在线播放| 国产在线视频一区二区| 亚洲一区二区三区欧美精品| 日韩中文字幕视频在线看片| 欧美大码av| a 毛片基地| 精品亚洲成a人片在线观看| 国产三级黄色录像| 国产精品偷伦视频观看了| 欧美另类亚洲清纯唯美| 超碰成人久久| 丁香六月天网| 国产成人精品在线电影| 热99re8久久精品国产| 人妻久久中文字幕网| www.熟女人妻精品国产| 超色免费av| 日韩视频一区二区在线观看| 午夜免费成人在线视频| 99九九在线精品视频| 视频区图区小说| 日韩三级视频一区二区三区| 日韩视频一区二区在线观看| tocl精华| 黑人欧美特级aaaaaa片| 老司机在亚洲福利影院| 亚洲精品国产av蜜桃| 18禁国产床啪视频网站| 中亚洲国语对白在线视频| www.精华液| 97精品久久久久久久久久精品| 国产av一区二区精品久久| 亚洲国产精品一区二区三区在线| 久久久久久亚洲精品国产蜜桃av| 咕卡用的链子| 波多野结衣av一区二区av| 国产精品av久久久久免费| 日韩一卡2卡3卡4卡2021年| 9191精品国产免费久久| 国产亚洲精品第一综合不卡| 国产精品一区二区在线观看99| 菩萨蛮人人尽说江南好唐韦庄| 99九九在线精品视频| 午夜福利视频在线观看免费| 亚洲国产精品999| 欧美亚洲 丝袜 人妻 在线| 亚洲精品久久午夜乱码| 美女视频免费永久观看网站| 亚洲国产日韩一区二区| 免费在线观看黄色视频的| 丁香六月欧美| 自线自在国产av| 丝袜美腿诱惑在线| 欧美黄色淫秽网站| a级毛片黄视频| 日韩欧美一区视频在线观看| 日本a在线网址| 亚洲欧洲日产国产| 亚洲欧美激情在线| 亚洲国产日韩一区二区| 亚洲国产av新网站| 老司机深夜福利视频在线观看 | 国产欧美日韩一区二区精品| 一级毛片精品| 自线自在国产av| 欧美激情高清一区二区三区| 中文字幕色久视频| 欧美日韩精品网址| 两个人免费观看高清视频| 久久久精品免费免费高清| 亚洲av日韩在线播放| 久久久久视频综合| 欧美日韩精品网址| 欧美激情 高清一区二区三区| 后天国语完整版免费观看| 免费一级毛片在线播放高清视频 | 久久国产精品人妻蜜桃| 精品国产超薄肉色丝袜足j| 男女床上黄色一级片免费看| 69精品国产乱码久久久| 真人做人爱边吃奶动态| 无遮挡黄片免费观看| av网站免费在线观看视频| 欧美少妇被猛烈插入视频| 午夜两性在线视频| 老汉色av国产亚洲站长工具| 国产麻豆69| 亚洲欧美清纯卡通| 久久狼人影院| 国产成人影院久久av| 看免费av毛片| av线在线观看网站| 亚洲七黄色美女视频| 正在播放国产对白刺激| 久久久久久久大尺度免费视频| 久久久国产精品麻豆| 大片电影免费在线观看免费| 动漫黄色视频在线观看| 亚洲av电影在线进入| 各种免费的搞黄视频| 久久人人爽人人片av| 视频区欧美日本亚洲| 精品欧美一区二区三区在线| 欧美成人午夜精品| 亚洲精品国产精品久久久不卡| 视频区图区小说| 欧美人与性动交α欧美软件| 亚洲欧美色中文字幕在线| 国产精品久久久久久人妻精品电影 | 国产91精品成人一区二区三区 | av片东京热男人的天堂| 手机成人av网站| 精品乱码久久久久久99久播| 视频区图区小说| 男人爽女人下面视频在线观看| 99精品欧美一区二区三区四区| 亚洲人成电影免费在线| 国产高清videossex| 12—13女人毛片做爰片一| 欧美日韩亚洲综合一区二区三区_| 国产精品 欧美亚洲| 飞空精品影院首页| 欧美日韩亚洲综合一区二区三区_| 日本猛色少妇xxxxx猛交久久| 亚洲国产毛片av蜜桃av| 法律面前人人平等表现在哪些方面 | 国产区一区二久久| 亚洲专区国产一区二区| 国产在线视频一区二区| 亚洲国产av影院在线观看| a 毛片基地| 欧美激情极品国产一区二区三区| 亚洲欧美成人综合另类久久久| 免费日韩欧美在线观看| 国产不卡av网站在线观看| 一区二区av电影网| 午夜福利影视在线免费观看| 无遮挡黄片免费观看| 中国美女看黄片| 午夜日韩欧美国产| 在线观看免费视频网站a站| 久久久国产成人免费| 欧美国产精品va在线观看不卡| kizo精华| 老司机午夜十八禁免费视频| 精品少妇一区二区三区视频日本电影| 岛国毛片在线播放| www.精华液| 老熟女久久久| 肉色欧美久久久久久久蜜桃| 青青草视频在线视频观看| 欧美av亚洲av综合av国产av| 欧美日韩成人在线一区二区| 夜夜骑夜夜射夜夜干| 日本精品一区二区三区蜜桃| 夜夜夜夜夜久久久久| 国产日韩欧美视频二区| 中文字幕人妻丝袜制服| 国产一区二区 视频在线| 国产野战对白在线观看| 国产成人影院久久av| 午夜福利影视在线免费观看| 人妻一区二区av| 18禁裸乳无遮挡动漫免费视频| 精品国产超薄肉色丝袜足j| 欧美黑人欧美精品刺激| 欧美人与性动交α欧美精品济南到| 免费在线观看视频国产中文字幕亚洲 | 男人添女人高潮全过程视频| 亚洲男人天堂网一区| 女人精品久久久久毛片| 亚洲熟女毛片儿| 天堂中文最新版在线下载| 色视频在线一区二区三区| 欧美在线黄色| 中文字幕人妻熟女乱码| 欧美精品亚洲一区二区| 久久久水蜜桃国产精品网| 亚洲精品国产一区二区精华液| 亚洲专区中文字幕在线| 电影成人av| 亚洲欧美一区二区三区久久| 国产又色又爽无遮挡免| 菩萨蛮人人尽说江南好唐韦庄| 黑人巨大精品欧美一区二区mp4| 天天添夜夜摸| 中国国产av一级| 亚洲国产看品久久| 18禁黄网站禁片午夜丰满| 国产精品一区二区精品视频观看| 一个人免费在线观看的高清视频 | 欧美 亚洲 国产 日韩一| 十八禁网站网址无遮挡| 高清av免费在线| 久久久久久久大尺度免费视频| 中文字幕av电影在线播放| 99九九在线精品视频| 欧美黄色片欧美黄色片| 日本wwww免费看| 99热全是精品| 精品国产乱码久久久久久男人| 久久综合国产亚洲精品| 国产一级毛片在线| 国产精品 国内视频| 99久久综合免费| 亚洲久久久国产精品| 777米奇影视久久| 成年av动漫网址| 国产黄频视频在线观看| 高清欧美精品videossex| 国产伦理片在线播放av一区| 一本色道久久久久久精品综合| 国产精品久久久久久精品古装| 美女大奶头黄色视频| av电影中文网址| 国产又爽黄色视频| 一区二区三区精品91| 99国产精品一区二区三区| 亚洲精品美女久久av网站| 午夜成年电影在线免费观看| 亚洲欧美日韩另类电影网站| 欧美黑人欧美精品刺激| 亚洲国产欧美在线一区| 成年人午夜在线观看视频| 国产精品久久久人人做人人爽| 少妇的丰满在线观看| 狠狠精品人妻久久久久久综合| 精品福利永久在线观看| 99久久综合免费| 捣出白浆h1v1| 久久久精品94久久精品| 男女午夜视频在线观看| 一区在线观看完整版| 一边摸一边做爽爽视频免费| 国产成人精品久久二区二区91| 日日夜夜操网爽| 亚洲男人天堂网一区| 多毛熟女@视频| 久久狼人影院| 国产成人一区二区三区免费视频网站| 啦啦啦啦在线视频资源| 欧美精品一区二区大全| 精品久久久精品久久久| 十八禁网站网址无遮挡| 亚洲欧洲精品一区二区精品久久久| 午夜福利视频在线观看免费| 欧美精品一区二区免费开放| 久热爱精品视频在线9| 亚洲一码二码三码区别大吗| 国产亚洲精品一区二区www | 窝窝影院91人妻| 国产伦理片在线播放av一区| www.精华液| 国产xxxxx性猛交| 桃花免费在线播放| 天天躁夜夜躁狠狠躁躁| 波多野结衣一区麻豆| 亚洲熟女毛片儿| 久久免费观看电影| 欧美精品啪啪一区二区三区 | 久久这里只有精品19| 每晚都被弄得嗷嗷叫到高潮| 亚洲五月婷婷丁香| 十分钟在线观看高清视频www| 中文字幕制服av| 日韩人妻精品一区2区三区| 成年人午夜在线观看视频| 精品国产乱码久久久久久男人| 精品少妇内射三级| 91麻豆精品激情在线观看国产 | 国产欧美亚洲国产| 亚洲avbb在线观看| 精品亚洲成a人片在线观看| 成人手机av| 久热爱精品视频在线9| 久久人妻福利社区极品人妻图片| 啦啦啦视频在线资源免费观看| 亚洲精品美女久久久久99蜜臀| 国产高清国产精品国产三级| 精品视频人人做人人爽| 各种免费的搞黄视频| 国产亚洲精品第一综合不卡| 日韩一卡2卡3卡4卡2021年| 午夜精品国产一区二区电影| 黄片大片在线免费观看| 黑人操中国人逼视频| 亚洲久久久国产精品| 精品国产乱子伦一区二区三区 | 久久人人97超碰香蕉20202| 久久性视频一级片| 韩国高清视频一区二区三区| 久久久久国内视频| 亚洲一区中文字幕在线| 免费观看av网站的网址| 国产亚洲精品一区二区www | 欧美日韩精品网址| 男女边摸边吃奶| 97精品久久久久久久久久精品| 一级片免费观看大全| 色老头精品视频在线观看| 69av精品久久久久久 | 狂野欧美激情性bbbbbb| 亚洲av日韩精品久久久久久密| 一区二区三区四区激情视频| 亚洲国产日韩一区二区| 国精品久久久久久国模美| 日韩三级视频一区二区三区| 窝窝影院91人妻| 女性生殖器流出的白浆| 午夜福利视频精品| 国产精品一区二区在线观看99| 9191精品国产免费久久| 99精国产麻豆久久婷婷| 美国免费a级毛片| 性色av乱码一区二区三区2| 超碰成人久久| 啪啪无遮挡十八禁网站| 亚洲国产精品成人久久小说| 男女下面插进去视频免费观看| 亚洲国产精品成人久久小说| av网站免费在线观看视频| 九色亚洲精品在线播放| a级毛片在线看网站| 99热全是精品| 久久精品国产亚洲av高清一级| 色播在线永久视频| 最近中文字幕2019免费版| 高清欧美精品videossex| 男女免费视频国产| a级毛片黄视频| 久久久久国内视频| 不卡一级毛片| 精品国产乱码久久久久久小说| 自线自在国产av| 十分钟在线观看高清视频www| 午夜影院在线不卡| 2018国产大陆天天弄谢| 日韩欧美一区视频在线观看| 国产高清国产精品国产三级| 日本精品一区二区三区蜜桃| 真人做人爱边吃奶动态| 亚洲九九香蕉| 亚洲一区中文字幕在线| 久久精品成人免费网站| 久久久精品区二区三区| 美女主播在线视频| 国产在线一区二区三区精| 亚洲免费av在线视频| 国产日韩欧美在线精品| www.精华液| 久久天堂一区二区三区四区| 91老司机精品| 久久久国产精品麻豆| 高清欧美精品videossex| 啪啪无遮挡十八禁网站| 日韩中文字幕欧美一区二区| 亚洲天堂av无毛| www.自偷自拍.com| 桃花免费在线播放| 国产精品香港三级国产av潘金莲| 夫妻午夜视频| 国产精品二区激情视频| 国产99久久九九免费精品| 国产成人精品在线电影| 老司机亚洲免费影院| 亚洲精华国产精华精| 性色av一级| 一区二区av电影网| 精品亚洲成a人片在线观看| 亚洲avbb在线观看| 1024香蕉在线观看| 日日夜夜操网爽|