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

    三維挑流水舌及其碰撞的LBM數(shù)值模擬研究

    2022-04-18 10:09:44彭燕祥吳國華
    關(guān)鍵詞:水舌外緣風場

    彭燕祥,張 華,吳國華

    (華北電力大學 水利與水電工程學院,北京 102206)

    1 研究背景

    隨著中國水電事業(yè)的不斷發(fā)展,水電站向高水頭、大流量的方向發(fā)展。泄洪霧化問題越發(fā)突出,尤其是采用挑流消能方式泄洪時,會在水電站下游局部區(qū)域形成大面積的雨霧。相關(guān)的研究表明:挑流泄洪霧化與水舌自身摻氣擴散、水舌之間相互碰撞和水舌入水噴濺密切相關(guān),要對挑流泄洪霧化進行更深入的研究,就必須對挑流水舌的運動特性,以及水舌與空氣的耦合作用引起的水舌風場進行研究。

    針對水舌的運動特性,許多學者建立不同的數(shù)學模型進行了研究,劉宣烈等[1-3]通過實驗對水舌運動過程的摻氣和擴散過程進行研究,得到了水舌運動的軌跡方程以及水舌橫、縱斷面的擴散規(guī)律。梁在潮等[4]對水舌運動軌跡、水舌摻氣量以及霧源量等參數(shù)進行了研究,并對這些參數(shù)進行擬合。張華等[5]通過建立水舌運動微分方程,并用數(shù)值方法進行求解,來研究空中水舌運動特性。劉士和等[6]運用相似理論在考慮空氣阻力及卷吸摻氣的作用下,建立了挑流水舌的數(shù)學模型。劉繼廣等[7]依據(jù)摻氣濃度相似律準則進行氣水兩相流模擬研究。張貝貝等[8]提出一種在排放口近區(qū)考慮射流垂向效應、在全計算域使用二維平面模型的計算方法。但是這些研究僅僅是對水舌進行了研究,由水舌引起的水舌風場并沒有涉及。

    為了研究水舌與空氣的耦合作用引起的水舌風場。李敏等[9]采用非靜力閉合邊界層模式,模擬了水電站泄洪條件下水舌風的發(fā)展過程。彭燕祥等[10]用CFD方法模擬了二灘水電站泄洪時水舌的運動及水舌風場的空間分布。彭莘仔等[11]采用N-S方程模擬射流的水動力過程。但是這些方法都是基于有限元網(wǎng)格的歐拉方法,這個方法對于流體運動的范圍有其局限性,網(wǎng)格的劃分密度以及范圍都對誤差的影響很大。

    格子玻爾茲曼方法(Lattice Boltzmann Method,LBM)是一種介觀流體系統(tǒng)模擬方法,兼具宏觀流體連續(xù)性和微觀分子動力學方法的優(yōu)點,而且計算效率較高,現(xiàn)已被廣泛應用于多組分多相流、化學反應流、氣固兩相流的研究[12-14]。在水汽兩相流兩方面,宋歌等[15]應用LBM方法對離心式霧化噴嘴出口處的液膜厚度、霧化錐角以及其他霧化特性進行了研究。葉永等[16]采用LBM對三維潰壩流進行了數(shù)值模擬與驗證。兩者都沒有對空氣流場進行研究。張華等[17]應用格子Boltzmann兩相流模型,對水舌在空氣中的運動進行了數(shù)值模擬,但是模擬的水舌是二維的,且研究的水舌尺度較小。

    為了分析三維挑流水舌的運動特性,以及兩股挑流水舌的碰撞特性,本文采用LBM方法對挑流水舌進行三維數(shù)值模擬研究。

    2 格子玻爾茲曼方法

    粒子運動的速度矢量見圖1。

    圖1 D3Q19網(wǎng)格模型與粒子速度向量

    3 挑流水舌空中運動數(shù)值模擬

    3.1 計算模型與參數(shù)設(shè)置 采用長江水科院[20]對泄洪霧化概化模型,進行數(shù)值模擬研究,如圖2所示,泄水建筑物長2.3 m,高1.7 m,寬0.28 m,出口挑坎挑角20°。數(shù)值模擬工況采用了典型的實驗工況,即水頭0.244 m,水頭落差2.05 m。

    圖2 數(shù)值模擬幾何模型

    3.2 模型參數(shù)設(shè)置 模型采用xflow兩相流計算模型,分別為液體相和空氣相,液體相密度為998.3 kg/m3,動力黏度為0.001 Pa·s??諝庀嗝芏葹?.225 kg/m3,動力黏度為1.7894 10-5Pa·s。

    由于時間步長和解析度直接關(guān)系到計算結(jié)果的精度和計算效率。所以在保證計算結(jié)果精度的同時還必須考慮計算時間。經(jīng)多方面考慮以及多次試算之后,選擇時間步長為0.001 s,解析度為0.03 m最為合適。

    環(huán)境參考壓力取一個標準大氣壓,邊界條件1、2截面采用壓力進口條件,3截面采用壓力出口,相對壓力都設(shè)置為0。其余面均設(shè)置為自由滑移壁面邊界條件。

    3.3 水舌風場計算結(jié)果 圖3為t=2.98 s時的水舌風場的空間分布計算圖,從圖中可以比較清楚的看出水舌與空氣運動交界面,挑流水舌隨著水舌的運動不斷的摻氣擴大。水舌的運動帶動周圍的空氣沿水舌方向運動,并且在水舌入水處,水舌表面空氣流速與水舌的流速相近風速達到6 m/s,風向沿水舌速度方向。

    圖3 水舌風場的空間分布

    由于水舌下緣在臨近下游水面區(qū)域與下游水面形成一個相對閉塞的空間,特別時水舌入水附近,如圖3(c)所示,水舌下緣的風場情況相較于水舌上緣相對更加的復雜,在水舌下沿出現(xiàn)較多的漩渦流場,如圖3(d)所示。從圖中可以看出水舌入水與下游水面碰撞,引起下游水面液滴飛濺,在下游水面形成一個與水舌入水面為邊界的水滴噴濺區(qū),并且沿水舌外緣入水區(qū)域的噴濺水滴,由于水舌入水夾角的關(guān)系,沿x方向的速度分量比其他區(qū)域的大。這也是泄洪霧化的主要霧源。

    3.4 水舌運動特性驗證分析 數(shù)值模擬采用了典型工況:水頭0.244 m,水頭落差2.05 m。對水舌的內(nèi)緣和外緣軌跡和挑距,以及水舌的厚度沿水舌軌跡的變化等進行分析。

    采用2次多項式對水舌的內(nèi)緣、外緣以及中心內(nèi)線進行擬合,分別得到:

    內(nèi)緣曲線表達式為:

    外緣曲線表達式為:

    中心線曲線表達式:

    挑流水舌的外緣和內(nèi)緣,分別與文獻[5]做了對比,如圖4所示。

    圖4 挑流水舌外緣和內(nèi)緣的空間變化

    從圖4中可以看出,兩者的水舌軌跡基本一致。在水舌的中間段,本文的外緣線比文獻[5]的結(jié)果偏高,中間段水舌的擴散速度比文獻[5]更快。

    圖5為水舌流速與下游水位截面位置的體積分布計算結(jié)果,水舌入水的速度在6 m/s左右。水舌內(nèi)緣挑距約為2.8 m,外緣挑距約為3.5 m,橫向?qū)挾燃s為0.75 m。與實測數(shù)據(jù)對比誤差分別為9.6%、7.9%、7.1%,如表1所示,可以看出挑距的模擬值相對實測數(shù)據(jù)偏低。水舌入水的形狀如圖5(b)所示,大致為一個圓角的矩形趨近于橢圓形,與文獻[21]中單孔泄洪實驗的水舌落水形狀基本一致,這可以為水舌入水噴濺霧源的進一步研究提供參考。

    表1 水舌挑距實測值[17]與模擬值的對比

    圖5 挑流水舌俯視圖與水舌入水截面圖

    圖6為水舌模擬與理論計算的厚度與長度關(guān)系變化圖,其中H為水舌厚度/m,H0為水舌初始斷面厚度/m,S為水舌長度/m。通過挑流水舌厚度與水舌長度的關(guān)系圖,可以將水舌的厚度變化可以分成三個階段,如圖的Ⅰ、Ⅱ、Ⅲ階段。

    圖6 挑流水舌厚度與水舌長度的關(guān)系

    Ⅰ階段水舌擴散緩慢,且水舌的厚度的增長與水舌長度呈線性關(guān)系,這是由于Ⅰ階段水舌尚有核心區(qū)存在,水舌摻氣量少,水體比較密實。

    Ⅱ階段水舌的厚度擴散速度先加快后慢慢變緩,這是因為由于水舌速度的增加,空氣對水舌的阻力增加,這樣水舌內(nèi)部的紊動也就隨之增加,引起水舌的破碎和大量摻氣,使水舌厚度快速的增加,隨著水舌摻氣的增加逐漸趨于飽和,使得摻氣速度逐漸趨于變緩,水舌的厚度增加速度變緩。

    Ⅲ階段為水舌的均勻擴散段,水舌的厚度的增長與水舌長度呈現(xiàn)近似線性變化。這與文獻[1]中通過實驗數(shù)據(jù)分析,將水舌厚度的變化分為三個不同變化區(qū)段,即:(1)出口緩擴段;(2)水舌急劇擴散段;(3)水舌均勻擴散段;結(jié)論基本一致。

    如圖7為水舌速度與水舌長度的關(guān)系對比圖,其中的文獻[5]-v、文獻[5]-vx、文獻[5]-vy為采用文獻[5]的數(shù)學模型,忽略空氣阻力情況下的計算結(jié)果。從圖中可以看出本文-v和文獻[5]-v隨著S增加先減小后增加、本文-v和文獻[5]-v的上升階段,本文-v相對與文獻[5]-v的值偏小,這是因為采用文獻[5]數(shù)學模型計算水舌時,沒有考慮空氣阻力以及水舌風場的影響。

    圖7 水舌速度與水舌長度的關(guān)系

    本文-vx值隨著水舌挑S的增加而減小,而文獻[5]-vx水舌由于沒有考慮空氣阻力的影響,水舌速度文獻[5]-vx的值隨著水舌挑距S的增加幾乎不變。

    本文-vy和文獻[5]-vy隨S增加先減小到0后增加,當本文-vy和文獻[5]-vy等于0時,水舌上升到最高點。本文-vy與文獻[5]-vy隨S變化趨勢基本是一致,這是由于在y方向水舌雖然受重力和空氣阻力以及水舌風場的共同作用,但是重力起主導作用。

    3.5 水舌渦量分析 渦量為流體速度的旋度:

    渦量是描寫旋渦運動的重要物理量之一,在xy平面取水舌的中心平面得水舌的中間截面的渦量云圖,即渦量在z軸的投影用Ωz表示,如圖8所示。從圖中可以看出水舌的上下表面在初始階段Ωz遠遠高于內(nèi)部水舌,這是因為水舌初始階段摻氣并不充分,僅僅外部水舌與空氣作用,使得水舌的上下表面的Ωz遠遠高于內(nèi)部水舌。這也是水舌摻氣的重要影響因素,隨著水舌的摻氣量不斷變大,Ωz向水舌的內(nèi)部輸運,在水舌的末端Ωz達到最大,Ωz達到139 s-1。

    圖8 水舌中心截面渦量云圖

    4 水舌碰撞的數(shù)值模擬

    格子Boltzmann方法可以更加容易的實現(xiàn)兩股水舌的碰撞,以及由于水舌碰撞而形成的相對復雜的水汽交界面,如圖9所示。

    圖9 兩股水舌碰撞水舌速度分布圖

    圖9和圖10分別為水舌碰撞的數(shù)值模擬圖與文獻[22]中水舌碰撞的實驗圖象。從圖中可以看出兩者的碰撞水舌形態(tài)基本一致,水舌經(jīng)過碰撞后,水花四濺,水舌在碰撞后形成一個錐形的射流形態(tài),并且由于兩股水舌的碰撞形成的壓力波影響,使得上股水舌的下表面破碎形成一股反向微射流,微射流以液滴為主且運動狀態(tài)較為復雜,但其大多數(shù)液滴的運動方向與水舌的運動方向相反。這為進一步研究多股水舌泄流,以及水舌入水噴濺等復雜情況提供參考。

    圖10 兩股水舌碰撞的試驗結(jié)果圖[22]

    5 結(jié)論

    本文采用格子Boltzmann方法針對三維挑流水舌以及兩股水舌的碰撞進行數(shù)值模擬研究,得到如下結(jié)論:

    (1)格子Boltzmann方法能很好的模擬泄洪挑流水舌和水舌風場,水舌的運動界面清晰可見。并且不用對水舌的周圍的大片區(qū)域進行網(wǎng)格劃分,計算更加便捷,效率更高。水舌內(nèi)緣挑距約為2.8 m,外緣挑距約為3.5 m,橫向?qū)挾燃s為0.75 m。與實測數(shù)據(jù)對比誤差分別為9.6%、7.9%、7.1%,挑距的模擬值相對實測數(shù)據(jù)偏低。為泄洪挑流水舌的流場研究提供新的方法。

    (2)挑流水舌在初始階段由于水舌與空氣的相互作用,使得水舌的上下表面的Ωz遠遠高于內(nèi)部水舌,這也是水舌摻氣的重要影響因素,隨著水舌的摻氣量不斷變大,Ωz向水舌的內(nèi)部輸運,在水舌的末端渦量達到最大,渦量達到139 s-1。

    (3)格子Boltzmann方法可以對兩股水舌的碰撞進行模擬,水舌碰撞的數(shù)值模擬與實驗條件下水舌運動形態(tài)基本一致,水舌經(jīng)過碰撞后,水花四濺,水舌在碰撞后形成一個錐形的射流形態(tài),并且由于兩股水舌的碰撞形成的壓力波影響,使得上股水舌的下表面破碎形成一股反向微射流。為多孔泄流的復雜情況,水舌的運動特性以及流場研究提供參考。

    猜你喜歡
    水舌外緣風場
    燕尾挑坎挑射水流數(shù)值模擬
    基于FLUENT的下?lián)舯┝魅S風場建模
    中國白肩天蛾屬記述與一新記錄種(鱗翅目:天蛾科)
    從外緣到內(nèi)蘊:《紅樓夢》異文的詮釋空間與小說特質(zhì)
    紅樓夢學刊(2019年2期)2019-04-12 05:49:56
    淺析藏醫(yī)學中的疾病產(chǎn)生機理
    “最美風場”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    渤、黃、東海冬季海表冷暖水舌的時空變化及機理分析
    海洋學研究(2017年1期)2017-04-21 03:33:06
    側(cè)向風場中無人機的飛行研究
    泄洪洞出口扭曲斜切挑坎挑流數(shù)值模擬
    從建筑文化傳承中解讀建筑史
    草草在线视频免费看| 免费黄色在线免费观看| 国产伦精品一区二区三区四那| 三级经典国产精品| 午夜福利在线在线| 舔av片在线| 下体分泌物呈黄色| 国产日韩欧美在线精品| 欧美xxⅹ黑人| 久久人人爽人人爽人人片va| 久久国内精品自在自线图片| 青春草视频在线免费观看| 在线天堂最新版资源| 日本色播在线视频| 91久久精品国产一区二区成人| 国产精品伦人一区二区| 欧美激情在线99| 51国产日韩欧美| 中文字幕久久专区| 久久久久久九九精品二区国产| 日韩三级伦理在线观看| 午夜爱爱视频在线播放| 美女被艹到高潮喷水动态| 国产成人精品久久久久久| av免费在线看不卡| 免费观看的影片在线观看| 久久久久久久久久久丰满| 成人免费观看视频高清| 99久久精品热视频| 婷婷色综合www| 伊人久久国产一区二区| 熟女电影av网| 丰满人妻一区二区三区视频av| 亚洲高清免费不卡视频| tube8黄色片| 国产高清不卡午夜福利| 精品午夜福利在线看| 一本久久精品| 在线观看人妻少妇| 97超视频在线观看视频| 少妇猛男粗大的猛烈进出视频 | 国产有黄有色有爽视频| av又黄又爽大尺度在线免费看| 国产在视频线精品| 免费人成在线观看视频色| 在线看a的网站| 五月玫瑰六月丁香| 最近2019中文字幕mv第一页| 国产成人一区二区在线| 少妇人妻 视频| 亚洲欧美日韩卡通动漫| 新久久久久国产一级毛片| av天堂中文字幕网| 国产成年人精品一区二区| 99久久精品国产国产毛片| 尤物成人国产欧美一区二区三区| 亚洲欧美一区二区三区黑人 | 男女下面进入的视频免费午夜| 国产av国产精品国产| 晚上一个人看的免费电影| 搞女人的毛片| 成人国产麻豆网| 亚洲精品影视一区二区三区av| 九九爱精品视频在线观看| 色综合色国产| 久久久久久久久大av| 日韩欧美精品免费久久| 亚洲三级黄色毛片| 汤姆久久久久久久影院中文字幕| 69人妻影院| 女的被弄到高潮叫床怎么办| 亚洲色图av天堂| 精品一区二区免费观看| 亚洲综合精品二区| 国产午夜精品一二区理论片| 日本色播在线视频| 性插视频无遮挡在线免费观看| 中文字幕免费在线视频6| 色综合色国产| 中国三级夫妇交换| 在线 av 中文字幕| 亚洲人成网站在线播| 日韩电影二区| 亚洲性久久影院| 国产精品福利在线免费观看| av卡一久久| 亚洲国产精品国产精品| 久久久久久九九精品二区国产| 欧美少妇被猛烈插入视频| 一级黄片播放器| 日韩在线高清观看一区二区三区| 街头女战士在线观看网站| 国产精品女同一区二区软件| 欧美 日韩 精品 国产| 亚洲国产日韩一区二区| 草草在线视频免费看| 亚洲综合精品二区| 亚洲精华国产精华液的使用体验| 少妇被粗大猛烈的视频| 国产精品嫩草影院av在线观看| 可以在线观看毛片的网站| 日韩av不卡免费在线播放| 激情 狠狠 欧美| 国产视频内射| 伦理电影大哥的女人| 日本wwww免费看| 午夜老司机福利剧场| 男女那种视频在线观看| 2022亚洲国产成人精品| 校园人妻丝袜中文字幕| 国产精品麻豆人妻色哟哟久久| 制服丝袜香蕉在线| 少妇裸体淫交视频免费看高清| 极品教师在线视频| 九九在线视频观看精品| 男女国产视频网站| 亚洲久久久久久中文字幕| 精品国产露脸久久av麻豆| 最近中文字幕高清免费大全6| 久久精品国产亚洲av天美| 你懂的网址亚洲精品在线观看| 五月玫瑰六月丁香| 狂野欧美激情性bbbbbb| 国产精品.久久久| 成人国产麻豆网| 高清欧美精品videossex| 精品久久久噜噜| 国内少妇人妻偷人精品xxx网站| 成人漫画全彩无遮挡| 香蕉精品网在线| 成人亚洲欧美一区二区av| 免费不卡的大黄色大毛片视频在线观看| 好男人在线观看高清免费视频| 国产午夜精品一二区理论片| 久久99热这里只有精品18| 高清欧美精品videossex| 美女被艹到高潮喷水动态| 国产精品成人在线| 天天躁日日操中文字幕| 国产在线男女| 亚洲欧洲国产日韩| 免费黄色在线免费观看| 人体艺术视频欧美日本| 欧美性猛交╳xxx乱大交人| 极品教师在线视频| 国产高潮美女av| 久久影院123| 熟妇人妻不卡中文字幕| 亚洲一区二区三区欧美精品 | 免费av不卡在线播放| 午夜激情久久久久久久| 欧美zozozo另类| 亚洲av在线观看美女高潮| 制服丝袜香蕉在线| 18禁在线播放成人免费| 久久久久久久亚洲中文字幕| 哪个播放器可以免费观看大片| 精品久久久久久久末码| 国产又色又爽无遮挡免| 色播亚洲综合网| 欧美日韩亚洲高清精品| 国产免费一级a男人的天堂| 久久久精品欧美日韩精品| 久久久久久久大尺度免费视频| av福利片在线观看| 亚洲欧洲日产国产| 日本三级黄在线观看| 午夜福利网站1000一区二区三区| 午夜福利视频精品| 欧美老熟妇乱子伦牲交| 日本wwww免费看| 国产人妻一区二区三区在| 十八禁网站网址无遮挡 | 亚洲欧美精品专区久久| 18禁裸乳无遮挡免费网站照片| 人妻 亚洲 视频| 免费观看a级毛片全部| 最后的刺客免费高清国语| 美女被艹到高潮喷水动态| 97人妻精品一区二区三区麻豆| 王馨瑶露胸无遮挡在线观看| 一级毛片aaaaaa免费看小| 国产一区二区三区av在线| 一区二区三区免费毛片| 搞女人的毛片| 国产亚洲午夜精品一区二区久久 | 97超视频在线观看视频| 成人二区视频| 日韩av不卡免费在线播放| 国产成人免费无遮挡视频| 亚洲三级黄色毛片| 亚洲伊人久久精品综合| 美女高潮的动态| 精品国产一区二区三区久久久樱花 | 哪个播放器可以免费观看大片| 国语对白做爰xxxⅹ性视频网站| 免费人成在线观看视频色| av国产免费在线观看| 久久久色成人| 亚洲最大成人中文| 亚洲自偷自拍三级| av线在线观看网站| freevideosex欧美| 亚洲精品日本国产第一区| 99re6热这里在线精品视频| 高清日韩中文字幕在线| 成人二区视频| 身体一侧抽搐| 国产免费一区二区三区四区乱码| 亚洲无线观看免费| 亚洲欧洲日产国产| 男女那种视频在线观看| 亚洲av成人精品一二三区| 亚洲精品日本国产第一区| 天天一区二区日本电影三级| 欧美国产精品一级二级三级 | 成年免费大片在线观看| 身体一侧抽搐| 欧美日韩国产mv在线观看视频 | 欧美亚洲 丝袜 人妻 在线| videossex国产| 国产精品国产av在线观看| 亚洲欧美成人综合另类久久久| 亚洲人与动物交配视频| 亚洲欧美精品自产自拍| av在线播放精品| 成人欧美大片| 国产欧美另类精品又又久久亚洲欧美| 蜜桃久久精品国产亚洲av| 日韩av在线免费看完整版不卡| 精品久久久久久久末码| 麻豆成人av视频| 大码成人一级视频| 国产欧美日韩一区二区三区在线 | 嘟嘟电影网在线观看| 我的老师免费观看完整版| 日日摸夜夜添夜夜爱| 久久久精品免费免费高清| 精品久久久久久电影网| 男人狂女人下面高潮的视频| 国产黄色免费在线视频| 国产精品久久久久久精品电影小说 | 在线天堂最新版资源| 超碰av人人做人人爽久久| 亚洲在线观看片| 美女主播在线视频| 天堂俺去俺来也www色官网| 网址你懂的国产日韩在线| 亚洲精品国产色婷婷电影| 老师上课跳d突然被开到最大视频| 欧美日韩在线观看h| 热99国产精品久久久久久7| 国产探花极品一区二区| 综合色av麻豆| 可以在线观看毛片的网站| 菩萨蛮人人尽说江南好唐韦庄| 一级黄片播放器| 亚洲电影在线观看av| 亚洲色图综合在线观看| 人人妻人人爽人人添夜夜欢视频 | 久久久亚洲精品成人影院| 国产免费一区二区三区四区乱码| 大香蕉久久网| 欧美日韩视频高清一区二区三区二| 中文字幕久久专区| 色吧在线观看| 丝瓜视频免费看黄片| 亚洲精品自拍成人| 一本一本综合久久| 夜夜看夜夜爽夜夜摸| 日本av手机在线免费观看| 久久精品熟女亚洲av麻豆精品| 菩萨蛮人人尽说江南好唐韦庄| 欧美丝袜亚洲另类| av在线蜜桃| 国产成人a∨麻豆精品| 亚洲国产精品999| 午夜激情久久久久久久| 亚洲在线观看片| 99久久精品热视频| 九九爱精品视频在线观看| av播播在线观看一区| 国产 精品1| 国产熟女欧美一区二区| 亚洲精品影视一区二区三区av| 欧美+日韩+精品| 久久久久久久午夜电影| 欧美97在线视频| 久久精品久久久久久噜噜老黄| 99热这里只有是精品在线观看| 纵有疾风起免费观看全集完整版| 99久久精品国产国产毛片| 成人亚洲精品一区在线观看 | 插阴视频在线观看视频| 黄色日韩在线| 最近中文字幕高清免费大全6| 少妇的逼水好多| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久久久久大尺度免费视频| 神马国产精品三级电影在线观看| 国产 精品1| 精品熟女少妇av免费看| 2018国产大陆天天弄谢| 精品一区二区三卡| 日韩三级伦理在线观看| 久久久a久久爽久久v久久| 精品国产三级普通话版| 午夜精品国产一区二区电影 | 狂野欧美激情性bbbbbb| 日韩人妻高清精品专区| 一区二区三区乱码不卡18| 国产成人精品福利久久| 久久精品久久精品一区二区三区| 在线精品无人区一区二区三 | 国产精品国产三级国产专区5o| 成人毛片a级毛片在线播放| 精品久久久精品久久久| 精品国产露脸久久av麻豆| av卡一久久| 欧美激情久久久久久爽电影| 久久久久久久精品精品| 特级一级黄色大片| 成人漫画全彩无遮挡| 国产午夜精品一二区理论片| 男女边吃奶边做爰视频| 精品久久久噜噜| 日韩大片免费观看网站| 老师上课跳d突然被开到最大视频| 亚洲av成人精品一二三区| 国产精品蜜桃在线观看| 精品国产乱码久久久久久小说| 国产69精品久久久久777片| 亚洲一区二区三区欧美精品 | 久久精品国产自在天天线| 欧美丝袜亚洲另类| 国产午夜精品久久久久久一区二区三区| 国产美女午夜福利| 午夜日本视频在线| 狂野欧美白嫩少妇大欣赏| 欧美精品一区二区大全| 熟女人妻精品中文字幕| 精品人妻偷拍中文字幕| 99热6这里只有精品| 在线观看av片永久免费下载| 精品人妻熟女av久视频| 成人毛片60女人毛片免费| 免费观看的影片在线观看| 国产av码专区亚洲av| 欧美日韩国产mv在线观看视频 | av免费在线看不卡| 在线观看人妻少妇| 精品午夜福利在线看| 国产精品麻豆人妻色哟哟久久| 免费电影在线观看免费观看| 亚洲天堂国产精品一区在线| 日韩三级伦理在线观看| 青春草视频在线免费观看| 中文精品一卡2卡3卡4更新| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧洲国产日韩| 中文字幕人妻熟人妻熟丝袜美| 亚洲伊人久久精品综合| 久久久a久久爽久久v久久| 色吧在线观看| 国产高清三级在线| 日韩av免费高清视频| 男女无遮挡免费网站观看| 欧美精品国产亚洲| 欧美潮喷喷水| 97在线视频观看| 亚洲av中文字字幕乱码综合| 在线观看国产h片| 看黄色毛片网站| 久久久久久国产a免费观看| 久久这里有精品视频免费| 久久久久久国产a免费观看| 久久久精品免费免费高清| 国国产精品蜜臀av免费| 日韩成人伦理影院| av国产精品久久久久影院| 亚洲欧美日韩东京热| 亚洲,一卡二卡三卡| 久久鲁丝午夜福利片| 国精品久久久久久国模美| 色综合色国产| 王馨瑶露胸无遮挡在线观看| 亚洲电影在线观看av| 亚洲av免费高清在线观看| 亚洲第一区二区三区不卡| 久久久色成人| 男人舔奶头视频| 韩国av在线不卡| 免费高清在线观看视频在线观看| 三级国产精品片| av播播在线观看一区| 亚洲精品中文字幕在线视频 | 五月开心婷婷网| 性插视频无遮挡在线免费观看| 国产高清不卡午夜福利| 男人和女人高潮做爰伦理| 日本黄大片高清| 久久久久精品性色| 蜜臀久久99精品久久宅男| 成年免费大片在线观看| 26uuu在线亚洲综合色| 亚洲美女搞黄在线观看| 亚洲av.av天堂| 亚洲人成网站高清观看| 亚洲欧洲日产国产| 麻豆成人av视频| 国产av国产精品国产| 国产黄片美女视频| 亚洲最大成人手机在线| 亚洲精品亚洲一区二区| 欧美丝袜亚洲另类| 最新中文字幕久久久久| 国产色爽女视频免费观看| 国产成年人精品一区二区| 国产淫语在线视频| 激情 狠狠 欧美| 在线精品无人区一区二区三 | 国产久久久一区二区三区| 久久久久久久久大av| 免费少妇av软件| 97超视频在线观看视频| 午夜视频国产福利| 亚洲精品久久午夜乱码| 欧美日韩综合久久久久久| 男人和女人高潮做爰伦理| 天天躁夜夜躁狠狠久久av| 免费在线观看成人毛片| 亚洲精品日韩av片在线观看| 97在线人人人人妻| 国产成人福利小说| 国产精品99久久久久久久久| 久久韩国三级中文字幕| 亚洲欧美一区二区三区黑人 | 一本久久精品| 天天躁夜夜躁狠狠久久av| 韩国高清视频一区二区三区| 黄片wwwwww| 久久精品国产a三级三级三级| 国产精品国产三级国产av玫瑰| 亚洲丝袜综合中文字幕| 观看美女的网站| 少妇高潮的动态图| 高清欧美精品videossex| 男的添女的下面高潮视频| 青春草国产在线视频| av免费观看日本| 日日啪夜夜爽| 国产精品三级大全| 不卡视频在线观看欧美| 97热精品久久久久久| av在线蜜桃| 国产乱来视频区| eeuss影院久久| 亚洲一区二区三区欧美精品 | 搡女人真爽免费视频火全软件| 97在线人人人人妻| 国国产精品蜜臀av免费| 久久鲁丝午夜福利片| 精品久久久精品久久久| 亚洲激情五月婷婷啪啪| 国产欧美日韩一区二区三区在线 | 久久亚洲国产成人精品v| 女人十人毛片免费观看3o分钟| 人体艺术视频欧美日本| 97精品久久久久久久久久精品| 亚洲欧美精品专区久久| 日日啪夜夜爽| 欧美少妇被猛烈插入视频| 国产探花极品一区二区| 婷婷色综合www| 日韩av在线免费看完整版不卡| 久热久热在线精品观看| 爱豆传媒免费全集在线观看| 欧美日本视频| 精品久久久久久久末码| 中文字幕av成人在线电影| 一级av片app| 久久久久久久久大av| 午夜免费观看性视频| 午夜免费鲁丝| 国产有黄有色有爽视频| 少妇熟女欧美另类| 久久久久精品久久久久真实原创| 国产免费又黄又爽又色| 超碰av人人做人人爽久久| 夫妻性生交免费视频一级片| 超碰av人人做人人爽久久| 亚洲经典国产精华液单| 黄色视频在线播放观看不卡| 黄片wwwwww| 国产成人a区在线观看| 久久影院123| 国产精品一区二区在线观看99| 国产淫片久久久久久久久| 黄色怎么调成土黄色| 国产老妇女一区| 亚洲人与动物交配视频| 视频区图区小说| 日韩电影二区| 亚洲欧美日韩卡通动漫| 两个人的视频大全免费| 18禁动态无遮挡网站| 五月玫瑰六月丁香| 亚洲在久久综合| 一级毛片我不卡| 亚洲精品自拍成人| 久久精品国产亚洲av天美| 有码 亚洲区| 十八禁网站网址无遮挡 | 少妇猛男粗大的猛烈进出视频 | 国产精品久久久久久久久免| av又黄又爽大尺度在线免费看| 91狼人影院| 亚洲va在线va天堂va国产| 汤姆久久久久久久影院中文字幕| 高清在线视频一区二区三区| 久久ye,这里只有精品| 免费av不卡在线播放| 久久精品夜色国产| 精品一区二区免费观看| 欧美高清成人免费视频www| 91精品伊人久久大香线蕉| 国产精品一区二区性色av| 美女视频免费永久观看网站| 欧美精品国产亚洲| 九草在线视频观看| 熟女人妻精品中文字幕| 国产精品av视频在线免费观看| 欧美潮喷喷水| 男男h啪啪无遮挡| 中文天堂在线官网| 欧美高清性xxxxhd video| 国产欧美另类精品又又久久亚洲欧美| 特级一级黄色大片| 五月天丁香电影| 亚洲精品日韩av片在线观看| 欧美97在线视频| 亚洲欧美日韩卡通动漫| 国产永久视频网站| 国产久久久一区二区三区| 一级毛片电影观看| 亚洲成人av在线免费| 在线免费观看不下载黄p国产| www.色视频.com| 日韩在线高清观看一区二区三区| 国语对白做爰xxxⅹ性视频网站| 自拍偷自拍亚洲精品老妇| 国产色爽女视频免费观看| 午夜免费鲁丝| 日韩欧美 国产精品| 亚洲久久久久久中文字幕| 欧美少妇被猛烈插入视频| 欧美最新免费一区二区三区| 99热全是精品| 亚洲精品乱久久久久久| 日本熟妇午夜| 午夜爱爱视频在线播放| 99久久精品热视频| 亚洲无线观看免费| 在现免费观看毛片| 日韩精品有码人妻一区| 国产精品国产av在线观看| 最后的刺客免费高清国语| 中文欧美无线码| 欧美97在线视频| 少妇猛男粗大的猛烈进出视频 | 九九爱精品视频在线观看| 综合色av麻豆| 肉色欧美久久久久久久蜜桃 | 大码成人一级视频| 岛国毛片在线播放| 水蜜桃什么品种好| 免费av毛片视频| av在线亚洲专区| 在线 av 中文字幕| 一区二区av电影网| 91狼人影院| 麻豆成人午夜福利视频| 亚洲国产av新网站| 欧美老熟妇乱子伦牲交| 一本一本综合久久| 精品国产三级普通话版| 日韩成人伦理影院| 久久久久网色| 精品久久久久久电影网| 国产精品99久久久久久久久| 国产人妻一区二区三区在| 精品一区二区免费观看| 在线精品无人区一区二区三 | 久久久欧美国产精品| 最近手机中文字幕大全| 亚洲成人中文字幕在线播放| 精品久久久久久久人妻蜜臀av| 18禁裸乳无遮挡免费网站照片| av免费观看日本| 午夜免费男女啪啪视频观看| 亚洲欧美清纯卡通| 黄片wwwwww| 91精品一卡2卡3卡4卡| 在现免费观看毛片| 99久久精品热视频| 久久精品国产亚洲av涩爱| 国产成人aa在线观看| 丰满人妻一区二区三区视频av| 久久久久国产网址| 99热6这里只有精品| 精品午夜福利在线看| 久久女婷五月综合色啪小说 | 久久久亚洲精品成人影院| 免费黄频网站在线观看国产| 日本爱情动作片www.在线观看|