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

    封凍河道冰下流速分布規(guī)律研究

    2022-09-08 13:05:00藺建銘張寶森謝志剛田治宗于國卿
    人民黃河 2022年9期
    關(guān)鍵詞:糙率冰蓋床面

    蘇 磊,藺建銘,張寶森,謝志剛,田治宗,鄧 宇,于國卿

    (1.黃河水利委員會 黃河水利科學(xué)研究院,河南 鄭州 450003;2.鄂爾多斯市水利事業(yè)發(fā)展中心,內(nèi)蒙古 鄂爾多斯 017000)

    1 引言

    高寒地區(qū)冬季氣溫低,河渠易結(jié)冰而形成冰蓋。冰蓋顯著改變河渠斷面的水流結(jié)構(gòu),使封凍河道各水力要素變得十分復(fù)雜[1],河道邊界條件、水動(dòng)力條件變化導(dǎo)致冰下水位、流量、流速變化,將影響冬季河渠輸水,甚至誘發(fā)凌汛災(zāi)害,因此研究封凍河道水流特性具有重要的理論及實(shí)際意義[2]。

    國內(nèi)外學(xué)者對冰下流速分布問題進(jìn)行了大量的研究,冰凌數(shù)學(xué)模型方面,國內(nèi)外開發(fā)了一系列一維動(dòng)/靜態(tài)模型(RICE、ICEJAM、RIVER1D、HEC-RAS等)、二維動(dòng)/靜態(tài)模型(CRISSP2-D、DYNARICEDE[3]等)。王軍等[4]、楊開林等[5]、茅澤育等[6]結(jié)合水力學(xué)、熱力學(xué)原理改進(jìn)了河道冰凌數(shù)學(xué)模型,但針對冰蓋下水流垂向流速分布特性的研究相對較少。

    本文對N-S方程進(jìn)行簡化、離散,采用k-ε紊流模型建立冰蓋下水流垂向二維數(shù)學(xué)模型,模擬封凍河道冰蓋下流速分布情況,并與實(shí)測資料進(jìn)行比較,分析、探討其變化規(guī)律。

    2 控制方程與離散

    2.1 控制方程

    如圖1所示,x表示流動(dòng)方向坐標(biāo),z表示垂直于河床方向坐標(biāo)。對于冰蓋下的明渠紊流,考慮靜壓假定的不可壓縮流體的控制方程無量綱化后,可以表示如下。

    圖1 坐標(biāo)系示意

    連續(xù)方程:

    x方向動(dòng)量方程:

    在z方向上引入σ坐標(biāo),某一坐標(biāo)處的垂向質(zhì)點(diǎn)流速w可用垂向場流速ω表示為

    式中:u、w分別為x、z方向的時(shí)均流速;t為時(shí)間;ζ為自由水面高程(ζ=zb+H),zb為床面高程;R0為羅斯比數(shù);為科氏系數(shù),g為重力加速度;ρ為密度;p為壓強(qiáng);AH為水平向紊動(dòng)渦黏系數(shù);AV為垂向紊動(dòng)渦黏系數(shù);H為水深。

    此外,式中定義如下無量綱組合:①Rossby數(shù),R0=Ur/fXr;Froude數(shù),Fr=Ur/(gZr)1/2。②密度Froude數(shù),。③水平向Ekman數(shù),EH=AHr/;垂向Ekman數(shù),EV=AVr/。其中下標(biāo)r表示對應(yīng)變量的參考值。

    垂向紊動(dòng)渦黏系數(shù)AV采用垂向k-ε模型進(jìn)行計(jì)算,該模型可以較好地反映風(fēng)剪切應(yīng)力、床面剪切應(yīng)力、流速梯度產(chǎn)生項(xiàng)、耗散、擴(kuò)散和分層等的影響。k-ε模型的基本思想是將垂向紊動(dòng)渦黏系數(shù)AV與紊流動(dòng)能k和紊流耗散率ε聯(lián)系起來,即

    紊流動(dòng)能k和紊流耗散率ε的輸運(yùn)方程表示如下:

    式中:cν、c1、c2、σk、σε為常數(shù),取值cν=0.09、c1=1.44、c2=1.92、σk=1.0、σε=1.3;源項(xiàng)P為流速垂向梯度引起的紊流產(chǎn)生項(xiàng);G為浮力產(chǎn)生或耗散項(xiàng)。

    河床阻力的確定,以x方向?yàn)槔?河床切應(yīng)力在x方向表示為,其中:U、V分別為x、y方向的平均流速;cf為(無量綱)床面摩阻系數(shù),cf=gn2/H1/3=g/C2,n為Manning糙率系數(shù),C為Chézy系數(shù);mb反映了河床底坡的影響。

    冰蓋阻力的確定,采用H.A.Einstein提出的阻力劃分方法可將流體分為兩層等效明流層(冰面層和床面層),假定上下兩層互不影響,分別只受冰蓋底部和床面粗糙度的影響,與表示河床阻力相同的方式,冰蓋底部的切應(yīng)力在x方向表示為。

    2.2 離散求解策略

    將x方向流速u表示為如下形式:

    式中:u′為x方向流速相較于垂向平均流速的差值,即攝動(dòng)流速。

    u′所滿足的方程可通過將三維控制方程減去二維控制方程得到,這樣的處理移除了原始方程中的自由水面重力波項(xiàng),因此求解u′后再通過式(7)間接計(jì)算u,能使模型獲得更高的穩(wěn)定性,保證計(jì)算質(zhì)量。

    對于攝動(dòng)動(dòng)量方程,除時(shí)間項(xiàng)外,只有垂直擴(kuò)散項(xiàng)被視為隱式的,其余各項(xiàng)(包括科氏力項(xiàng)、對流項(xiàng)、水平擴(kuò)散項(xiàng)和斜壓項(xiàng))均被視為顯式的,也就是說,攝動(dòng)動(dòng)量方程可以簡寫為如下形式:

    式中:Rx為x方向的水平源項(xiàng),包括科氏力項(xiàng)、對流項(xiàng)、水平擴(kuò)散項(xiàng)和斜壓項(xiàng)4部分。

    將式(8)采用有限體積法進(jìn)行離散,將其在一個(gè)時(shí)間步長Δt和一個(gè)控制體積Δσ上進(jìn)行積分可得

    如圖2所示,假設(shè)變量在中心點(diǎn)的值代表該變量在整個(gè)控制體積Δσ上的積分平均值,并認(rèn)為隱式項(xiàng)在n、n+1兩個(gè)時(shí)刻的加權(quán)平均值代表其在整個(gè)時(shí)間步長Δt上的積分平均值,顯式項(xiàng)在n時(shí)刻的值代表其在整個(gè)時(shí)間步長Δt上的積分平均值,則上式可表示為

    圖2 垂向網(wǎng)格示意

    式中:θ為權(quán)重系數(shù);k為垂向分層數(shù)。

    代入公式(9),得到

    上式即式(8)的最終離散形式,并對k=2,3,…,m-1成立,其中m為垂向分層總數(shù),并可寫成通用形式:

    其中

    2.3 定解條件及方程組的求解

    2.3.1 邊界條件

    (1)床面邊界條件(k=1),如圖3所示。

    圖3 床面邊界條件示意

    式(8)在床面附近所需滿足的邊界條件為

    代入式(9),得到

    式(13)寫成通用形式為

    其中

    (2)冰蓋邊界條件(k=n),如圖4所示。

    圖4 冰蓋邊界條件示意

    式(8)在冰蓋附近所需滿足的邊界條件為

    代入式(9),得到

    式(9)寫成通用形式為

    其中

    2.3.2 初始條件

    2.3.3 方程組求解

    離散方程[式(12)]結(jié)合床面邊界條件[式(14)]、自由水面邊界條件[式(16)],最終形成如下三對角方程組,該方程組可用“追趕法”求解。

    求解上述三對角方程組獲得攝動(dòng)流速u′后,可由式(7)計(jì)算得到每層流速u。整個(gè)計(jì)算的迭代收斂條件為,其中ε取10-4。

    3 計(jì)算結(jié)果分析

    應(yīng)用實(shí)測資料對上述模型計(jì)算結(jié)果進(jìn)行了驗(yàn)證,如圖5、圖6所示,圖中h為水深、y為相對水深、為斷面平均流速。其中,圖5為計(jì)算值與文獻(xiàn)[7]中現(xiàn)場實(shí)測資料(River Moskva上長為400 m的封凍順直河段第4號量測斷面垂線流速分布實(shí)測成果)的比較情況,圖6(a)為計(jì)算值與文獻(xiàn)[8]中室內(nèi)水槽實(shí)測資料(槽寬80 cm,水深8.92 cm)的比較情況,圖6(b)為計(jì)算值與文獻(xiàn)[8]中第12號水流條件(水深200 cm,冰蓋糙率與床面糙率相同)的比較情況。結(jié)果顯示,模型的計(jì)算值與實(shí)測值均吻合較好。同時(shí),流速分布對于冰蓋糙率ni和床面糙率nb的影響較敏感,最大流速Umax的位置偏向糙率較小的一側(cè),且在相同的水流條件下,ni/nb的值越小,最大流速點(diǎn)的位置越靠近糙率較小的一側(cè);當(dāng)冰蓋糙率和床面糙率相等時(shí),流速分布以最大流速Umax為軸呈對稱分布,且在最大流速Umax處流速梯度幾乎為零[9]。

    圖5 計(jì)算值與文獻(xiàn)[7]流速分布比較

    圖6 計(jì)算值與文獻(xiàn)[8]流速分布比較

    將冰蓋下水流以最大流速為界分為床面層和冰面層兩層,分層示意如圖7所示,hb為床面層水深,yb為床面層相對水深,hi為冰面層水深,yi為冰面層相對水深。以文獻(xiàn)[8]中兩組工況為例子,工況1為床面糙率和冰蓋糙率不相等時(shí)(ni/nb=0.76),將其計(jì)算結(jié)果繪于圖8所示的對數(shù)坐標(biāo)系下;工況2為床面糙率和冰蓋糙率相等時(shí)(ni/nb=1),將其計(jì)算結(jié)果繪于圖9所示的對數(shù)坐標(biāo)系下。由圖8可知,工況1當(dāng)ln(yb/hb)<-0.4時(shí)(即yb/hb<0.67時(shí)),床面層流速分布符合對數(shù)分布,當(dāng)ln(yb/hb)>-0.4時(shí)(即yb/hb>0.67時(shí)),床面層流速分布將不再符合對數(shù)分布,且流速值略小于符合對數(shù)分布的情況;在冰面層有同樣的規(guī)律但影響區(qū)域略有不同,當(dāng)ln(yi/hi)<-0.32時(shí)(即yi/hi<0.72時(shí)),冰面層流速分布符合對數(shù)分布;ln(yi/hi)>-0.32時(shí)(即yi/hi>0.72時(shí)),冰面層流速分布將不再符合對數(shù)分布,且流速值略小于符合對數(shù)分布的情況。同時(shí),由圖9可知,工況2當(dāng)ln(yb/hb)<-0.41時(shí)(即yb/hb<0.66時(shí)),床面層和冰面層流速分布符合對數(shù)分布;ln(yb/hb)>-0.41時(shí)(即yb/hb>0.66時(shí)),床面層和冰面層流速分布將不再符合對數(shù)分布。

    圖8 對數(shù)坐標(biāo)系下流速分布結(jié)果(文獻(xiàn)[8]中,ni/nb=0.76)

    圖9 對數(shù)坐標(biāo)系下流速分布結(jié)果(文獻(xiàn)[8]中,ni/nb=1)

    計(jì)算結(jié)果表明,封凍河道的斷面流速分布主要與床面糙率、冰蓋糙率及河道的水流條件等影響因素相關(guān)。床面糙率與冰蓋糙率相等和床面糙率與冰蓋糙率不相等兩種工況,冰面層和床面層均有70%左右的流速符合對數(shù)分布,其余靠近最大流速點(diǎn)約30%范圍內(nèi)的流速不再符合對數(shù)分布且流速值略小于對數(shù)分布的情況。工況1當(dāng)床面糙率與冰蓋糙率相等時(shí),冰面阻力和床面阻力干擾影響的區(qū)域相等,最大流速Umax在水深中點(diǎn)的位置。工況2當(dāng)床面糙率與冰蓋糙率不相等時(shí),冰面阻力和床面阻力干擾影響的區(qū)域不相等,且糙率越大阻力越大,阻力干擾影響的區(qū)域越大,從而導(dǎo)致最大流速Umax位置偏向相對光滑的一側(cè);Umax的位置與ni/nb比值相關(guān),且ni/nb的值越小Umax位置越靠近相對光滑的一側(cè)。床面糙率與冰蓋糙率相等和不相等兩種情況,計(jì)算結(jié)果均表現(xiàn)出相同的規(guī)律,這與前人的試驗(yàn)成果吻合。

    4 結(jié)論

    本文基于淺水方程,采用k-ε紊流模型建立冰蓋下水流二維數(shù)值模型,并將其應(yīng)用于封凍河道冰下流速分布的計(jì)算。結(jié)果表明,冰蓋下水流流速分布對于冰蓋糙率ni和床面糙率nb的影響較敏感,當(dāng)床面糙率與冰蓋糙率不相等時(shí),冰面阻力和床面阻力干擾影響的區(qū)域不相等,從而導(dǎo)致最大流速Umax位置偏向相對光滑的一側(cè),Umax位置的偏離程度與ni/nb值直接相關(guān);當(dāng)冰蓋糙率和床面糙率相等時(shí),流速分布以最大流速Umax為軸呈對稱分布。床面糙率與冰蓋糙率相等和不相等兩種情況,計(jì)算結(jié)果均顯示在最大流速點(diǎn)附近30%的區(qū)域內(nèi)流速分布并不遵循對數(shù)分布規(guī)律,且流速值略小于對數(shù)分布的情況。

    猜你喜歡
    糙率冰蓋床面
    軍事文摘(2024年6期)2024-04-30 03:13:59
    梯形河道行洪能力與護(hù)岸糙率的關(guān)系研究
    陜西水利(2023年7期)2023-07-28 09:30:26
    魚鱗狀床面粗糙特性
    格陵蘭島的冰蓋悄悄融化
    參花(下)(2022年1期)2022-01-15 00:45:01
    對瓦里安碳纖維治療床面模型的評估
    淹沒植物明渠床面沖淤及其對水流運(yùn)動(dòng)的影響
    改進(jìn)的投影覆蓋方法對遼河河道粗糙床面分維量化研究
    新疆阿勒泰哈巴河縣養(yǎng)殖渠人工渠道糙率的試驗(yàn)分析
    陜西水利(2019年5期)2019-06-26 06:27:34
    復(fù)式河道整治設(shè)計(jì)中綜合糙率研究
    大口徑玻璃鋼管道糙率及過流能力分析
    波多野结衣巨乳人妻| 亚洲欧美成人精品一区二区| 日韩一区二区视频免费看| 国产亚洲精品av在线| 少妇高潮的动态图| 日韩欧美 国产精品| 2021天堂中文幕一二区在线观| 91麻豆精品激情在线观看国产| 国产精品三级大全| 中国美女看黄片| 少妇熟女欧美另类| 十八禁国产超污无遮挡网站| 色5月婷婷丁香| 亚洲中文字幕日韩| 美女大奶头视频| 天堂√8在线中文| 99国产精品一区二区蜜桃av| 成人二区视频| 欧美日本视频| 日韩人妻高清精品专区| 亚洲av.av天堂| 国产精品爽爽va在线观看网站| 免费搜索国产男女视频| 色视频www国产| 国产精品麻豆人妻色哟哟久久 | 亚洲欧洲日产国产| 搡老妇女老女人老熟妇| 久久九九热精品免费| 如何舔出高潮| 日韩国内少妇激情av| 99久久久亚洲精品蜜臀av| 成人国产麻豆网| 久久久久久久午夜电影| 好男人视频免费观看在线| 日日摸夜夜添夜夜爱| 1024手机看黄色片| 亚洲人成网站在线播| 尤物成人国产欧美一区二区三区| 最近最新中文字幕大全电影3| 麻豆一二三区av精品| 可以在线观看毛片的网站| 婷婷精品国产亚洲av| 日日啪夜夜撸| 有码 亚洲区| 午夜免费激情av| 欧美成人精品欧美一级黄| 啦啦啦韩国在线观看视频| 免费av不卡在线播放| 日韩欧美精品v在线| 亚洲欧洲日产国产| 好男人在线观看高清免费视频| 别揉我奶头 嗯啊视频| 欧美+日韩+精品| 精品人妻偷拍中文字幕| 草草在线视频免费看| 超碰av人人做人人爽久久| 久99久视频精品免费| 观看免费一级毛片| 免费av不卡在线播放| 亚洲成人精品中文字幕电影| 国产精品日韩av在线免费观看| 亚洲欧美清纯卡通| 国产乱人偷精品视频| 国内精品一区二区在线观看| avwww免费| 综合色丁香网| 干丝袜人妻中文字幕| 国产精品野战在线观看| av在线蜜桃| 久久精品国产清高在天天线| 在线观看午夜福利视频| 桃色一区二区三区在线观看| 晚上一个人看的免费电影| 国产亚洲精品久久久com| 丰满乱子伦码专区| 欧美精品国产亚洲| 久久久久久久久久黄片| 国产男人的电影天堂91| 亚洲国产欧洲综合997久久,| 国产精品一区二区三区四区久久| 免费不卡的大黄色大毛片视频在线观看 | 国产免费一级a男人的天堂| 在线观看免费视频日本深夜| 国产极品精品免费视频能看的| 亚洲第一区二区三区不卡| 91久久精品电影网| 99久久精品热视频| 岛国在线免费视频观看| 国产精品伦人一区二区| 国产成人一区二区在线| 亚洲精品亚洲一区二区| 久久精品综合一区二区三区| 久久人人爽人人片av| 久久久久久久久中文| 国产女主播在线喷水免费视频网站 | 精品久久久久久久人妻蜜臀av| 亚洲国产高清在线一区二区三| 看非洲黑人一级黄片| 嫩草影院新地址| 免费看av在线观看网站| 久久久久免费精品人妻一区二区| 婷婷色综合大香蕉| 综合色丁香网| 亚洲人成网站在线播| 成年版毛片免费区| 久久久成人免费电影| 性插视频无遮挡在线免费观看| 别揉我奶头 嗯啊视频| av福利片在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 一边摸一边抽搐一进一小说| 十八禁国产超污无遮挡网站| 欧美一区二区精品小视频在线| 国产黄a三级三级三级人| 午夜精品在线福利| 99精品在免费线老司机午夜| 亚洲精品亚洲一区二区| 国产av不卡久久| 国产精品嫩草影院av在线观看| 亚洲天堂国产精品一区在线| 一进一出抽搐gif免费好疼| 长腿黑丝高跟| 99热全是精品| 亚洲在久久综合| 欧美日韩精品成人综合77777| 欧美性猛交╳xxx乱大交人| 亚洲成a人片在线一区二区| 中文在线观看免费www的网站| 亚洲人成网站在线播放欧美日韩| 天堂√8在线中文| 一本久久精品| 乱系列少妇在线播放| 日韩亚洲欧美综合| 亚洲在线观看片| 国产一级毛片七仙女欲春2| 国产一区亚洲一区在线观看| 成人av在线播放网站| 欧美潮喷喷水| 久久人人精品亚洲av| 黄色视频,在线免费观看| 国产视频内射| 老熟妇乱子伦视频在线观看| 亚洲国产欧洲综合997久久,| 在线免费十八禁| 亚洲在线自拍视频| 日韩一本色道免费dvd| 亚洲熟妇中文字幕五十中出| 日本一二三区视频观看| 老女人水多毛片| 欧美三级亚洲精品| 长腿黑丝高跟| 一级黄片播放器| 亚洲无线观看免费| 亚洲成人中文字幕在线播放| 男女视频在线观看网站免费| 亚洲av成人精品一区久久| 日本黄大片高清| 少妇的逼好多水| 熟女人妻精品中文字幕| 中文资源天堂在线| 少妇人妻一区二区三区视频| 好男人视频免费观看在线| 亚洲天堂国产精品一区在线| 精品一区二区三区人妻视频| 欧美性感艳星| 精品99又大又爽又粗少妇毛片| 网址你懂的国产日韩在线| 久久99热这里只有精品18| 免费观看的影片在线观看| 国产在视频线在精品| 精品熟女少妇av免费看| 床上黄色一级片| 麻豆乱淫一区二区| 免费av观看视频| 天天一区二区日本电影三级| 午夜a级毛片| 国产精品久久视频播放| 国产色婷婷99| 又粗又爽又猛毛片免费看| 嫩草影院精品99| 日本撒尿小便嘘嘘汇集6| 国产免费男女视频| 国产三级在线视频| 国产黄色小视频在线观看| 啦啦啦啦在线视频资源| 国产精品一区二区三区四区免费观看| 亚洲精品色激情综合| 午夜精品在线福利| 久久久久免费精品人妻一区二区| 久久婷婷人人爽人人干人人爱| 免费观看的影片在线观看| 国产探花在线观看一区二区| 99国产极品粉嫩在线观看| 综合色丁香网| 久久人妻av系列| 亚洲中文字幕一区二区三区有码在线看| 99久久精品热视频| or卡值多少钱| 18禁在线播放成人免费| 国产亚洲av嫩草精品影院| 免费看日本二区| 成人一区二区视频在线观看| 欧美不卡视频在线免费观看| 精品人妻偷拍中文字幕| 国产成人91sexporn| 欧美日韩乱码在线| 人妻少妇偷人精品九色| 国产在视频线在精品| 伊人久久精品亚洲午夜| 日韩人妻高清精品专区| 免费av毛片视频| 综合色丁香网| 乱系列少妇在线播放| 国产精品野战在线观看| 黄片无遮挡物在线观看| 能在线免费看毛片的网站| 亚洲丝袜综合中文字幕| 欧美日韩在线观看h| 少妇猛男粗大的猛烈进出视频 | 精华霜和精华液先用哪个| 最近中文字幕高清免费大全6| 高清毛片免费看| 国内精品美女久久久久久| 欧美一区二区国产精品久久精品| 久久久国产成人精品二区| 亚洲av二区三区四区| 亚洲18禁久久av| 在线免费观看不下载黄p国产| 久久久欧美国产精品| 婷婷色综合大香蕉| 精品一区二区三区视频在线| 欧美+日韩+精品| 丰满的人妻完整版| 国产成人aa在线观看| 内射极品少妇av片p| 日韩制服骚丝袜av| 精品一区二区三区视频在线| 国产精品一区二区三区四区免费观看| av免费观看日本| 九九热线精品视视频播放| 久久精品国产99精品国产亚洲性色| 亚洲成人av在线免费| 国产男人的电影天堂91| 蜜桃亚洲精品一区二区三区| 国内精品美女久久久久久| 国产一区二区三区在线臀色熟女| 在线播放国产精品三级| 波多野结衣高清无吗| 一级毛片久久久久久久久女| 久久久久久久亚洲中文字幕| 18禁黄网站禁片免费观看直播| 插阴视频在线观看视频| 看黄色毛片网站| 26uuu在线亚洲综合色| 久久九九热精品免费| 在线播放无遮挡| www.色视频.com| 国产淫片久久久久久久久| 精品少妇黑人巨大在线播放 | 日本av手机在线免费观看| 亚洲国产精品久久男人天堂| 男女啪啪激烈高潮av片| 六月丁香七月| 欧美日韩国产亚洲二区| 久久久久久久久久黄片| 免费看av在线观看网站| 午夜福利在线观看免费完整高清在 | 国产单亲对白刺激| 国产熟女欧美一区二区| 一级二级三级毛片免费看| 九九爱精品视频在线观看| 欧美一区二区精品小视频在线| 狠狠狠狠99中文字幕| 国产在线精品亚洲第一网站| 一区福利在线观看| 欧美xxxx性猛交bbbb| .国产精品久久| a级毛片a级免费在线| 内地一区二区视频在线| 日日啪夜夜撸| 久久国内精品自在自线图片| 国产午夜精品论理片| 久久久久久久久久黄片| 成熟少妇高潮喷水视频| 日韩av不卡免费在线播放| 少妇猛男粗大的猛烈进出视频 | 乱人视频在线观看| 久久人妻av系列| 国产精品无大码| 禁无遮挡网站| 白带黄色成豆腐渣| 国产一区二区在线av高清观看| 欧洲精品卡2卡3卡4卡5卡区| 欧美日韩精品成人综合77777| 美女 人体艺术 gogo| 高清午夜精品一区二区三区 | 国内精品美女久久久久久| 性插视频无遮挡在线免费观看| 亚洲,欧美,日韩| 色哟哟·www| 真实男女啪啪啪动态图| 男女做爰动态图高潮gif福利片| 久久精品久久久久久久性| 成人毛片60女人毛片免费| 18禁在线无遮挡免费观看视频| 精品少妇黑人巨大在线播放 | 久久精品国产亚洲网站| 国产女主播在线喷水免费视频网站 | 最新中文字幕久久久久| 女人被狂操c到高潮| 欧美一区二区精品小视频在线| 亚洲最大成人av| 欧美另类亚洲清纯唯美| 欧美日韩在线观看h| 在线免费十八禁| 亚洲中文字幕日韩| 国产成人a∨麻豆精品| 99久久人妻综合| 日韩av不卡免费在线播放| 国内精品宾馆在线| 国产伦在线观看视频一区| 欧美日韩精品成人综合77777| 51国产日韩欧美| 日韩成人av中文字幕在线观看| 一级毛片aaaaaa免费看小| 国产精品国产三级国产av玫瑰| 国产精品蜜桃在线观看 | 我要看日韩黄色一级片| 国产欧美日韩精品一区二区| 99热只有精品国产| 亚洲国产精品久久男人天堂| 亚洲欧美日韩无卡精品| 美女 人体艺术 gogo| 成人毛片60女人毛片免费| 亚洲欧美清纯卡通| 在线天堂最新版资源| av视频在线观看入口| 国产成人一区二区在线| 国产成人福利小说| 在线免费观看不下载黄p国产| 国产69精品久久久久777片| 中文字幕人妻熟人妻熟丝袜美| 又黄又爽又刺激的免费视频.| 日韩视频在线欧美| 黄色配什么色好看| 中文字幕av成人在线电影| 国产精品久久久久久精品电影| 中文字幕熟女人妻在线| 五月伊人婷婷丁香| 国产成人精品婷婷| 性插视频无遮挡在线免费观看| 国产69精品久久久久777片| 成人欧美大片| 欧美+日韩+精品| 麻豆久久精品国产亚洲av| 美女国产视频在线观看| 成人午夜高清在线视频| 国产精品一二三区在线看| 国产黄片美女视频| 国产黄片视频在线免费观看| 亚洲欧美日韩高清在线视频| 自拍偷自拍亚洲精品老妇| 亚洲人与动物交配视频| 在线观看免费视频日本深夜| 久久午夜福利片| 欧美最黄视频在线播放免费| 亚洲精品乱码久久久久久按摩| 亚洲欧美日韩高清专用| 国产免费男女视频| 亚洲自偷自拍三级| 免费观看精品视频网站| 国产高清激情床上av| 内地一区二区视频在线| 91午夜精品亚洲一区二区三区| 亚洲av一区综合| 亚洲欧美成人精品一区二区| 美女黄网站色视频| 好男人在线观看高清免费视频| 亚洲中文字幕一区二区三区有码在线看| 我要看日韩黄色一级片| 久久久久久久久久成人| 美女国产视频在线观看| 成人鲁丝片一二三区免费| 久久99蜜桃精品久久| 国产精品久久视频播放| av天堂中文字幕网| 中国国产av一级| 99热精品在线国产| 日韩制服骚丝袜av| 久久久久网色| 久久久久久九九精品二区国产| 我要看日韩黄色一级片| 国产av一区在线观看免费| 国产成人91sexporn| 国产中年淑女户外野战色| 欧美变态另类bdsm刘玥| 成人特级av手机在线观看| 精品99又大又爽又粗少妇毛片| 免费无遮挡裸体视频| 欧美在线一区亚洲| 午夜精品一区二区三区免费看| 亚洲精品久久久久久婷婷小说 | 欧美成人一区二区免费高清观看| 天堂av国产一区二区熟女人妻| 看非洲黑人一级黄片| 亚洲人成网站在线播| 国产女主播在线喷水免费视频网站 | 中文字幕av在线有码专区| 精品国产三级普通话版| 久久精品夜色国产| 黄色日韩在线| 久久九九热精品免费| 婷婷亚洲欧美| 色播亚洲综合网| 波多野结衣巨乳人妻| 亚洲美女视频黄频| 国产久久久一区二区三区| 日韩人妻高清精品专区| 亚洲欧美日韩卡通动漫| 波野结衣二区三区在线| 欧美性感艳星| 亚洲最大成人中文| 麻豆国产97在线/欧美| 高清毛片免费看| 久久久久久伊人网av| 久久草成人影院| 男女下面进入的视频免费午夜| 久久人人爽人人片av| 天堂中文最新版在线下载 | 日韩av在线大香蕉| h日本视频在线播放| 国国产精品蜜臀av免费| 精品久久久久久久久久久久久| 亚洲在线观看片| 青春草亚洲视频在线观看| 少妇熟女欧美另类| 日韩视频在线欧美| 国产黄片美女视频| 欧美日韩精品成人综合77777| 成人无遮挡网站| 精品99又大又爽又粗少妇毛片| 亚洲国产日韩欧美精品在线观看| 国产精品人妻久久久久久| 国产午夜福利久久久久久| 亚洲自偷自拍三级| 亚洲高清免费不卡视频| 欧美最黄视频在线播放免费| 国产精品一区二区性色av| 国产成人福利小说| 日本欧美国产在线视频| 国产在线精品亚洲第一网站| 日韩一本色道免费dvd| 日日撸夜夜添| 如何舔出高潮| 夜夜看夜夜爽夜夜摸| 91久久精品国产一区二区成人| 日本在线视频免费播放| 欧美日本视频| 小说图片视频综合网站| 黑人高潮一二区| 国产大屁股一区二区在线视频| 美女被艹到高潮喷水动态| 人人妻人人澡人人爽人人夜夜 | 又粗又硬又长又爽又黄的视频 | av又黄又爽大尺度在线免费看 | av免费在线看不卡| 欧美日韩乱码在线| 国产单亲对白刺激| 欧美潮喷喷水| 亚洲乱码一区二区免费版| 国产精品爽爽va在线观看网站| 成人鲁丝片一二三区免费| 欧美人与善性xxx| 26uuu在线亚洲综合色| 性色avwww在线观看| 欧美一级a爱片免费观看看| 国产一级毛片七仙女欲春2| 乱人视频在线观看| 亚洲av成人精品一区久久| 九九在线视频观看精品| 波多野结衣高清作品| 国产精品久久久久久精品电影小说 | 国产亚洲av片在线观看秒播厂 | 成人性生交大片免费视频hd| av.在线天堂| 精品欧美国产一区二区三| 熟女电影av网| 最近中文字幕高清免费大全6| 亚洲av二区三区四区| 一级毛片久久久久久久久女| 国产欧美日韩精品一区二区| 在线观看美女被高潮喷水网站| 天堂av国产一区二区熟女人妻| 久久这里只有精品中国| 成人综合一区亚洲| 日韩欧美三级三区| 国产精品人妻久久久久久| АⅤ资源中文在线天堂| 色尼玛亚洲综合影院| 身体一侧抽搐| 国产精品人妻久久久影院| 高清在线视频一区二区三区 | 欧美一级a爱片免费观看看| 哪个播放器可以免费观看大片| 在线免费观看不下载黄p国产| 国产一区二区在线观看日韩| 麻豆av噜噜一区二区三区| 国产亚洲5aaaaa淫片| 精品99又大又爽又粗少妇毛片| 97人妻精品一区二区三区麻豆| 99热6这里只有精品| 成人国产麻豆网| 精品99又大又爽又粗少妇毛片| 亚洲av成人av| 乱码一卡2卡4卡精品| 国产精品爽爽va在线观看网站| av专区在线播放| 亚洲综合色惰| 赤兔流量卡办理| 欧美精品国产亚洲| 欧美人与善性xxx| 免费观看人在逋| 中文字幕久久专区| 一区二区三区四区激情视频 | 18禁在线无遮挡免费观看视频| 九草在线视频观看| 黄色视频,在线免费观看| 在线观看66精品国产| 一区二区三区免费毛片| 成人毛片a级毛片在线播放| 精品欧美国产一区二区三| 国产精品一区二区在线观看99 | 国产成人午夜福利电影在线观看| 亚洲国产欧洲综合997久久,| 久久精品国产亚洲av涩爱 | 国产精品伦人一区二区| 亚洲一级一片aⅴ在线观看| 99久久成人亚洲精品观看| 亚洲国产欧美人成| 久久这里有精品视频免费| 国产成人精品婷婷| 久久鲁丝午夜福利片| 91在线精品国自产拍蜜月| 我的老师免费观看完整版| av女优亚洲男人天堂| 久久精品国产自在天天线| 久久精品久久久久久久性| 高清午夜精品一区二区三区 | av福利片在线观看| 好男人视频免费观看在线| 波多野结衣高清无吗| 国产精品一区www在线观看| 久久精品91蜜桃| 久久午夜福利片| 老熟妇乱子伦视频在线观看| 一进一出抽搐动态| 91午夜精品亚洲一区二区三区| 九九爱精品视频在线观看| 国产女主播在线喷水免费视频网站 | 黄片无遮挡物在线观看| 国产精品麻豆人妻色哟哟久久 | 国产色婷婷99| 91在线精品国自产拍蜜月| 久久精品国产亚洲av涩爱 | 色综合色国产| 91久久精品国产一区二区成人| 中文字幕av成人在线电影| 亚洲高清免费不卡视频| 国产伦理片在线播放av一区 | 丝袜喷水一区| 国产男人的电影天堂91| 麻豆一二三区av精品| 97热精品久久久久久| 女同久久另类99精品国产91| 夜夜爽天天搞| av在线观看视频网站免费| 国产美女午夜福利| 婷婷六月久久综合丁香| 在线播放无遮挡| 日本av手机在线免费观看| 永久网站在线| 97热精品久久久久久| 黄片无遮挡物在线观看| 青春草视频在线免费观看| 亚洲内射少妇av| 欧美最新免费一区二区三区| 深爱激情五月婷婷| 国产午夜福利久久久久久| 人妻久久中文字幕网| 色综合亚洲欧美另类图片| 在线播放无遮挡| 精品久久久久久久久av| 国产激情偷乱视频一区二区| 国产伦精品一区二区三区视频9| 亚洲欧美日韩高清在线视频| 成年女人永久免费观看视频| 国产老妇女一区| 男插女下体视频免费在线播放| 欧美最新免费一区二区三区| 麻豆精品久久久久久蜜桃| 免费不卡的大黄色大毛片视频在线观看 | 国产午夜精品论理片| 成人午夜高清在线视频| 一区二区三区高清视频在线| 青春草国产在线视频 | 国产精品精品国产色婷婷| 啦啦啦韩国在线观看视频| 国产成人a∨麻豆精品| 欧美人与善性xxx| 激情 狠狠 欧美| 熟妇人妻久久中文字幕3abv| 免费观看人在逋| 亚洲无线观看免费| 99九九线精品视频在线观看视频| 久久精品国产亚洲网站|