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

    風(fēng)速引起的湍流對揮發(fā)性污染物傳質(zhì)的影響

    2012-10-12 09:35:56陳麗萍程璟濤蔣軍成鄧廣發(fā)
    關(guān)鍵詞:交界面水氣傳質(zhì)

    陳麗萍,程璟濤,蔣軍成,鄧廣發(fā)

    (1.南京工業(yè)大學(xué)城市建設(shè)與安全工程學(xué)院,江蘇南京 210009;2.江蘇方天電力技術(shù)有限公司,江蘇南京 211102)

    準確預(yù)測揮發(fā)性污染物泄漏后在水體和大氣中的時空分布,能為水氣生態(tài)環(huán)境健康診斷、控制管理提供重要依據(jù)。揮發(fā)性污染物泄漏后在水氣中的時空分布受到諸多因素影響,如污染物種類,水、氣及其交界面的流態(tài)等。由于水氣界面上揮發(fā)性污染物揮發(fā)傳質(zhì)的復(fù)雜性,迄今為止,人們對其機理尚不清楚。本文探討了由風(fēng)速引起的水、氣及其交界面的湍流對揮發(fā)性污染物揮發(fā)傳質(zhì)的作用。水氣交界面附近的流動結(jié)構(gòu)性質(zhì)很獨特:如分界面通常不平整;在水流速度非常高時會形成自然摻氣,水流速度較高時水面上有波出現(xiàn),水流速度較低時,分界面本身變形不大,氣側(cè)的流動性質(zhì)與氣流繞過固體邊界的流動性質(zhì)相似,而水側(cè)卻并非如此[1-2]。這主要是因為水氣交界面水側(cè)的時均流速梯度很小,因而在水氣交界面附近耗散的湍動能只能來源于自由液面以下水側(cè)產(chǎn)生的湍動能及水氣交界面以上通過壓強場傳遞的湍動能[3]。從這一點來看,水氣交界面附近的流動結(jié)構(gòu)極大地依賴于水側(cè)的渦與水氣交界面的碰撞。但當(dāng)氣流速度較高時,水氣交界面附近的流動結(jié)構(gòu)也與氣側(cè)的湍流運動有關(guān)。

    目前還很難通過實驗得到水氣交界面附近流場的準確數(shù)據(jù),數(shù)值模擬成為研究水氣交界面湍流結(jié)構(gòu)的有效方法。根據(jù)數(shù)值模擬的精細程度,湍流數(shù)值模擬可分為直接模擬(DNS)、大渦模擬(LES)和雷諾平均數(shù)值模擬(RANS)[4]。Kawamura等[5-6]將水氣交界面上的速度和剪切應(yīng)力耦合,用DNS方法模擬了水氣交界面的湍流結(jié)構(gòu),但由于DNS方法對計算機內(nèi)存和計算速度的要求非常高,目前該方法無法用于工程實際。Magnaudet等[7]用LES方法研究了水氣交界面上湍流結(jié)構(gòu)及傳質(zhì)過程。Hasegawa等[8-9]提出了DNS和LES相混合的方法:在水氣交界面上的模擬用DNS方法;而水、氣交界面附近的水、氣模擬用 LES方法。Hasegawa等[8-9]雖然考慮了水、氣流動的動量耦合,但沒有考慮污染物濃度耦合。

    本文用氣液兩相流VOF(volume of fluid)方法捕捉水氣交界面[10],用可實現(xiàn)k-ε模型計算水氣交界面湍流結(jié)構(gòu),在此基礎(chǔ)上利用揮發(fā)性污染物水氣耦合擴散模型預(yù)測揮發(fā)性污染物點源泄漏后的時空分布,分析風(fēng)速引起的水氣湍流對水氣交界面上揮發(fā)性污染物揮發(fā)傳質(zhì)的作用機理。

    1 數(shù)學(xué)模型

    1.1 VOF動力學(xué)模型

    氣液兩相流VOF方法是常見的捕捉水氣交界面的方法,該方法中質(zhì)量守恒、動量方程及捕捉水氣交界面體積百分比函數(shù)α方程如下[11]:

    式中:ui——網(wǎng)格單元總速度;ρ——網(wǎng)格單元總密度,由液體、氣體的體積加權(quán)平均得出;uLIQi——液體速度;uGASi——氣體速度;ρGAS——氣體密度;p——壓強;ηt——湍流運動黏性系數(shù);fj——質(zhì)量力,本文中質(zhì)量力僅有重力。水氣交界面上,0<α<1。

    1.2 湍流計算模型

    其中

    式中:S——應(yīng)變率;k——湍動能;ε——湍動能耗散率;σm,σn,C2,C1——湍流系數(shù),σm=1.0,σn=1.2,C2=1.92;η——分子黏性系數(shù)。

    ηt由下式確定:

    其中

    式中:A0,As,U*——系數(shù);Ωi,ji,j——渦通量;ωk—— 角速度;δi,j,k——里奇(Ricci)符號,j——從角速度為ωk的參考系中觀察到的渦通量。

    1.3 揮發(fā)性污染物水氣耦合擴散模型

    由亨利定律表達式知,污染物在氣相中的質(zhì)量濃度CGAS與液相中的質(zhì)量濃度CLIQ之比為量綱為1的亨利常數(shù)Haw,即Haw=CGAS/CLIQ。忽略化學(xué)反應(yīng),揮發(fā)性污染物水氣耦合擴散模型表達式為

    式中:C——網(wǎng)格單元總質(zhì)量濃度;Ei,D——污染物在液相、氣相中湍流質(zhì)擴散系數(shù),為ηt與Schmidt數(shù)之比。

    2 模擬分析

    2.1 模擬條件

    計算域長直水槽14.0 m×0.4 m×0.4 m,空氣部14.0 m×6.4 m×2.6 m。坐標(biāo)原點取在水槽中部靠岸,水槽底部。水流動方向為x正方向,水槽寬度方向為y,垂直向上為z方向。進口邊界劃分為水流速度進口和空氣速度進口,進口的k,ε值由經(jīng)驗公式給出。假設(shè)出口斷面為充分發(fā)展的湍流,所有物理量一階法向梯度為零,50 mL質(zhì)量濃度為2 g/L的三氯乙烯溶液在x=1.25 m,y=0.20 m處瞬時泄漏。三氯乙烯的Haw為0.49,在水和空氣中的Schmidt數(shù)分別取754和1.03[13]。在交錯網(wǎng)格上,數(shù)值離散用有限體積法,速度壓力耦合用PISO算法,時間差分采用Crank-Nicholson格式。風(fēng)速太高,水氣交界面將變?yōu)榉沁B續(xù)界面,而本文研究連續(xù)水氣交界面上風(fēng)速引起的湍流對揮發(fā)性污染物揮發(fā)傳質(zhì)的影響,故模擬風(fēng)速較低,分別取0.4 m/s,0.8 m/s和1.2 m/s。模擬平臺為開源計算流體力學(xué)軟件OpenFOAM。

    2.2 k和ε與風(fēng)速的關(guān)系

    圖1是風(fēng)速為0.8 m/s,1.2 m/s時,x=3.45 m截面上水氣交界面附近水側(cè)k,ε相對值隨時間變化的曲線。如果考察的截面太近,風(fēng)速對揮發(fā)傳質(zhì)的作用時間太短,分析無意義。如果考察的截面太遠,揮發(fā)傳質(zhì)太少,研究意義不大。

    圖1 x=3.45 m截面水中 k/kv air=0.4 m/s,ε/εv air=0.4 m/s隨時間的變化Fig.1 k/kv air=0.4m/s and ε/εv air=0.4m/s as functions of time at cross-section x=3.45 m

    圖1中vair=0.4 m/s表示風(fēng)速為0.4 m/s的工況。圖中k/kvair=0.4m/s和ε/εvair=0.4m/s數(shù)值基本大于1,說明風(fēng)速為0.8 m/s,1.2 m/s時,水側(cè)k,ε大于風(fēng)速0.4 m/s時水側(cè)k,ε。這是因為風(fēng)速增加,水氣交界面附近時均速度梯度增加,進而導(dǎo)致湍動能平均值增加,同時湍動能波動幅度提高,氣流誘導(dǎo)的切應(yīng)力及表面變形隨著風(fēng)速增加而促進湍能耗散率增加。圖1中(a),(b)對應(yīng)的曲線變化趨勢一致,大多數(shù)時刻,ε/εvair=0.4m/s大于k/kvair=0.4m/s,說明湍能耗散率隨風(fēng)速增加而增加的程度高于湍動能隨風(fēng)速增加而增加的程度。

    2.3 揮發(fā)性污染物濃度分布

    圖2是不同風(fēng)速下,污染物泄漏8 s后,縱截面y=0.225m上污染物質(zhì)量濃度的空間分布。由圖2可見,污染物質(zhì)量濃度分布區(qū)域在縱向上得到了延伸,并且濃度中心由水中經(jīng)過水氣交界面向空氣中轉(zhuǎn)移,進而使空氣中的污染物越來越多。水流速度決定污染團前鋒的分布狀況,由于3種情況的水流速度相同,故他們的濃度前鋒相差無幾。3種風(fēng)速下污染團尾部區(qū)域相差甚遠。圖2(a)情況下,污染團幾乎沒有濃度尾巴;圖2(b)的污染團后部拖著長長的濃度尾巴,質(zhì)量濃度值較小;圖2(c)的污染團縱向延伸程度最大,濃度尾巴長且數(shù)值高。說明風(fēng)速的增加對污染團尾部的揮發(fā)傳質(zhì)過程影響巨大。由此說明,風(fēng)速引起的k,ε的增加不僅有利于水氣交界面的傳質(zhì)而且有助于污染物在空氣中的擴散。這種水氣中污染物濃度相互耦合的特性使得污染團尾部的濃度分布占據(jù)的空間隨風(fēng)速增加而增大。傳統(tǒng)非耦合模型認為空氣中污染物濃度為零,水中污染物揮發(fā)不受空氣中污染物濃度影響。事實上,空氣中污染物濃度抑制了水中污染物的揮發(fā),使得整個揮發(fā)過程得到延遲。

    圖3是x=3.45 m截面水氣交界面附近水側(cè)和氣側(cè)污染物質(zhì)量濃度隨時間變化的曲線。

    圖2 縱截面y=0.225 m污染物質(zhì)量濃度空間分布Fig.2 Pollutant concentration distribution at longitudinal section y=0.225 m

    圖3 x=3.45 m截面處不同風(fēng)速下水側(cè)、氣側(cè)污染物質(zhì)量濃度隨時間的變化Fig.3 Pollutant concentration distribution at cross-section x=3.45 m in water and air at different wind velocities

    圖3(a)中,0~4 s期間3條曲線基本沒變化,說明泄漏初期,風(fēng)速對水中污染物的揮發(fā)傳質(zhì)作用很小。但是從第4秒開始,曲線開始發(fā)生分離,3種風(fēng)速對應(yīng)的污染物質(zhì)量濃度開始增長,vair=0.8 m/s和vair=1.2 m/s對應(yīng)的質(zhì)量濃度增長較快。第10秒時,3種風(fēng)速下污染物質(zhì)量濃度均達到峰值,vair=1.2 m/s的峰值明顯高于其他2種。在質(zhì)量濃度達到峰值前,風(fēng)速越高,水中污染物質(zhì)量濃度越高。質(zhì)量濃度達到峰值后,同一時刻,vair=0.4 m/s與vair=0.8 m/s對應(yīng)質(zhì)量濃度幾乎相等??梢妚air=0.8 m/s時質(zhì)量濃度下降速率大于vair=0.4 m/s時質(zhì)量濃度的下降速率。vair=1.2 m/s對應(yīng)的質(zhì)量濃度稍高些,但其下降速率最大。說明達到峰值濃度后,風(fēng)速越大,ε越大,水中污染物質(zhì)量濃度下降越快。第6秒時,圖3(b)中3條曲線開始分離,氣側(cè)污染物質(zhì)量濃度達到峰值的時間要早于水側(cè)污染物,說明風(fēng)速對氣側(cè)污染物擴散的影響要先于對水側(cè)污染物擴散的影響。風(fēng)速越大,氣側(cè)污染物質(zhì)量濃度越高,質(zhì)量濃度達到峰值以后的衰減速率也越大,說明風(fēng)速引起的k,ε的增加有助于水中污染物向空氣中揮發(fā)傳質(zhì)及空氣中污染物的遷移擴散。

    3 結(jié) 語

    將揮發(fā)性污染物在水氣交界面上揮發(fā)傳質(zhì)的耦合特性與可實現(xiàn)k-ε模型相結(jié)合,揭示湍流特性對揮發(fā)傳質(zhì)的作用機理。

    風(fēng)速增加,水氣交界面附近時均速度梯度增加,進而導(dǎo)致湍動能平均值增加,同時湍動能波動幅度提高,氣流誘導(dǎo)的切應(yīng)力及表面變形隨著風(fēng)速增加而促進湍能耗散率增加。湍能耗散率隨風(fēng)速增加而增加的程度高于湍動能隨風(fēng)速增加而增加的程度。

    風(fēng)速引起的湍動能和湍能耗散率的提高增強了污染團尾部的揮發(fā)傳質(zhì),風(fēng)速越大,污染團尾部的濃度分布占據(jù)的空間越大。風(fēng)速對氣側(cè)污染物擴散的影響要先于對水側(cè)污染物擴散的影響。

    [1]LIU S,KERMANI A,SHEN L,et al.Investigation of coupled air-water turbulent boundary layers using direct numerical simulations[J].Phys Fluids,2009,21:062108.

    [2]YAMASHITA S,CHEN C G,TAKAHASHI K,et al.Large scale numerical simulations for multi-phase fluid dynamics with moving interfaces[J].International Journal of Computation Fluicl Dynamics,2008,22:405-410.

    [3]李佳佳,陳春剛,肖鋒.風(fēng)波水氣界面湍流渦旋結(jié)構(gòu)研究[J].中國科學(xué):物理學(xué) 力學(xué) 天文學(xué),2011,41:980-994.(LI Jiajia,CHEN Chungang,XIAO Feng.Study on vortex structure near wind-driven air-water interfaces[J].Sci Sin Phys Mech Astron,2011,41:980-994.(in Chinese))

    [4]張兆順,崔桂香,許春曉.湍流理論與模擬[M].北京:清華大學(xué)出版社,2005:163-165.

    [5]KAWAMURA H,OHSAKE K,ABE H,et al.DNS of turbulent heat transfer in channel flow with low to medium-high Prandtl number[J].International Journal Heat Fluid Flow,1998,19(5):482-491.

    [6]LAKEHAL D,F(xiàn)ULGOSI M,YADIGAROGLU G,et al.Direct numerical simulation of turbulent heat transfer across a mobile,sheared gas-liquid interface[J].ASME Journal Heat Transfer,2003,125(6):1129-1139.

    [7]MAGNAUDET J,CALMET I.Turbulent mass transfer through a flat shear-free surface[J].Fluid Mech Journal,2006,553:155-185.

    [8]HASEGAWA Y,KASAGI N.Effects of interfacial velocity boundary condition on turbulent mass transfer at high Schmidt numbers[J].International Journal Heat Fluid Flow,2007,28(6):1192-1203.

    [9]HASEGAWA Y,KASAGI N.Hybrid DNS/LES of high Schmidt number mass transfer across turbulent air-water interface[J].International Journal of Heat and Mass Transfer,2009,52:1012-1022.

    [10]DOLBOW JT,MOSSO S,ROBBINSJ,et al.Coupling volume-of-fluid based interface reconstructions with the extended finite element method [J].Computer Methods in Applied Mechanics and Engineering,2008,197:439-447.

    [11]CHEN Lipin,JIANG Juncheng.Experiments and numerical simulations on transport of dissolved pollutants around spur dike[J].Water Science and Engineering,2010,3(3):343-355.

    [12]陶文銓.數(shù)值傳熱學(xué)[M].西安:西安交通大學(xué)出版社,2001:278-280.

    [13]施瓦茨巴赫 R P,施格文 P M,英博登 D M.環(huán)境有機化學(xué)[M].北京:化學(xué)工業(yè)出版社,2004:425-432.

    猜你喜歡
    交界面水氣傳質(zhì)
    鋼-混凝土交界面法向粘結(jié)性能研究
    遼中區(qū)患病草魚體內(nèi)嗜水氣單胞菌分離、鑒定與致病力測定
    海上邊水氣藏利用試井資料確定水侵狀況研究
    海洋石油(2021年3期)2021-11-05 07:42:26
    高速公路機電工程相關(guān)交界面管理組織建設(shè)探討
    雙塊式無砟軌道軌枕與道床交界面損傷特性分析
    中國鐵路(2019年1期)2019-03-23 01:11:58
    氨基酸鹽吸收二氧化碳過程的傳質(zhì)特性
    PTFE膜吸收CO2的工藝參數(shù)對傳質(zhì)性能的影響
    清潔轉(zhuǎn)向酸H+表面?zhèn)髻|(zhì)行為實驗研究
    煤顆粒熱解的傳熱傳質(zhì)分析
    計算物理(2014年1期)2014-03-11 17:00:36
    醫(yī)院感染嗜水氣單胞菌的臨床治療分析
    欧美黄色片欧美黄色片| 两个人免费观看高清视频| 1024视频免费在线观看| 男女午夜视频在线观看| 午夜91福利影院| 久久久久视频综合| 久久国产精品大桥未久av| 日本欧美视频一区| 国产成+人综合+亚洲专区| bbb黄色大片| 男女之事视频高清在线观看| 日韩欧美国产一区二区入口| 9191精品国产免费久久| 老司机午夜十八禁免费视频| 在线永久观看黄色视频| 脱女人内裤的视频| 国产亚洲午夜精品一区二区久久| 黄片大片在线免费观看| 男女之事视频高清在线观看| 精品人妻熟女毛片av久久网站| 国产真人三级小视频在线观看| 一区二区三区激情视频| 曰老女人黄片| 黄色视频在线播放观看不卡| 不卡一级毛片| 国产精品熟女久久久久浪| 老司机深夜福利视频在线观看 | 美女扒开内裤让男人捅视频| 男人操女人黄网站| 国产熟女午夜一区二区三区| 国产精品成人在线| 亚洲精品中文字幕一二三四区 | 搡老岳熟女国产| 乱人伦中国视频| 丝瓜视频免费看黄片| 真人做人爱边吃奶动态| 黄色怎么调成土黄色| 欧美一级毛片孕妇| 99久久国产精品久久久| 久久久久久久久久久久大奶| 国产欧美日韩精品亚洲av| 久久中文看片网| 考比视频在线观看| 国产在线一区二区三区精| 建设人人有责人人尽责人人享有的| 国产在视频线精品| 少妇人妻久久综合中文| 精品福利永久在线观看| 最近中文字幕2019免费版| 国产成人精品在线电影| 在线永久观看黄色视频| 亚洲精品一二三| 亚洲国产看品久久| 久久精品aⅴ一区二区三区四区| 久久精品熟女亚洲av麻豆精品| 国产又爽黄色视频| 久久久精品免费免费高清| 久久毛片免费看一区二区三区| 一区二区日韩欧美中文字幕| 久久亚洲精品不卡| 美女国产高潮福利片在线看| av电影中文网址| 高清欧美精品videossex| 欧美精品亚洲一区二区| 精品久久久精品久久久| 91九色精品人成在线观看| 丝袜美足系列| 成年人午夜在线观看视频| 欧美精品啪啪一区二区三区 | 精品少妇一区二区三区视频日本电影| 五月天丁香电影| 久久精品亚洲熟妇少妇任你| 国产成+人综合+亚洲专区| 啦啦啦在线免费观看视频4| 搡老乐熟女国产| 欧美另类亚洲清纯唯美| 亚洲精品乱久久久久久| 免费在线观看黄色视频的| tocl精华| 国产精品成人在线| 欧美成狂野欧美在线观看| 这个男人来自地球电影免费观看| 天天操日日干夜夜撸| 国产伦人伦偷精品视频| 国产av又大| 高清av免费在线| 91精品国产国语对白视频| 久久久久国产精品人妻一区二区| 午夜影院在线不卡| 久久久水蜜桃国产精品网| 国产精品av久久久久免费| 国产av精品麻豆| 十八禁网站网址无遮挡| 成人亚洲精品一区在线观看| 国产精品二区激情视频| 91av网站免费观看| 丝袜人妻中文字幕| 国产精品成人在线| 欧美日韩成人在线一区二区| 精品熟女少妇八av免费久了| 日本av免费视频播放| tocl精华| 中文字幕人妻熟女乱码| 女警被强在线播放| 国产一区二区在线观看av| 日本wwww免费看| 99热国产这里只有精品6| 成人黄色视频免费在线看| 两性夫妻黄色片| 国产色视频综合| 两性夫妻黄色片| 成人手机av| 成人手机av| 欧美日韩一级在线毛片| 国产精品偷伦视频观看了| av超薄肉色丝袜交足视频| 天堂中文最新版在线下载| 大香蕉久久成人网| 极品人妻少妇av视频| 777米奇影视久久| 一级毛片女人18水好多| 亚洲午夜精品一区,二区,三区| 999久久久国产精品视频| 在线永久观看黄色视频| 母亲3免费完整高清在线观看| 久久久久久亚洲精品国产蜜桃av| 欧美大码av| 色老头精品视频在线观看| 动漫黄色视频在线观看| 日本91视频免费播放| 国产精品香港三级国产av潘金莲| 精品乱码久久久久久99久播| 欧美日韩亚洲综合一区二区三区_| 久久久久久久久久久久大奶| 国产欧美日韩综合在线一区二区| 精品一区在线观看国产| 亚洲 欧美一区二区三区| 精品卡一卡二卡四卡免费| 亚洲专区字幕在线| 亚洲专区国产一区二区| 三上悠亚av全集在线观看| 巨乳人妻的诱惑在线观看| 国产精品国产av在线观看| 精品人妻熟女毛片av久久网站| 99精国产麻豆久久婷婷| 伊人亚洲综合成人网| 久久国产精品大桥未久av| 国产欧美日韩精品亚洲av| 黄色视频不卡| 欧美日韩亚洲综合一区二区三区_| 69av精品久久久久久 | 捣出白浆h1v1| 国产精品自产拍在线观看55亚洲 | 每晚都被弄得嗷嗷叫到高潮| 在线永久观看黄色视频| 在线观看免费视频网站a站| 亚洲视频免费观看视频| 人妻一区二区av| 午夜福利影视在线免费观看| 亚洲av男天堂| 男女高潮啪啪啪动态图| 黄色怎么调成土黄色| 精品卡一卡二卡四卡免费| 欧美精品啪啪一区二区三区 | 国产成人欧美| 啦啦啦免费观看视频1| 在线观看一区二区三区激情| 男女无遮挡免费网站观看| 一区二区日韩欧美中文字幕| 1024香蕉在线观看| 国产伦理片在线播放av一区| 美女高潮到喷水免费观看| 亚洲av欧美aⅴ国产| 69精品国产乱码久久久| 亚洲伊人久久精品综合| 最新在线观看一区二区三区| 一级片'在线观看视频| 国产片内射在线| 日日摸夜夜添夜夜添小说| 美国免费a级毛片| 五月天丁香电影| 又紧又爽又黄一区二区| 在线亚洲精品国产二区图片欧美| 国产日韩欧美亚洲二区| 女人高潮潮喷娇喘18禁视频| 美女中出高潮动态图| 国产免费av片在线观看野外av| 人人妻人人爽人人添夜夜欢视频| 欧美激情极品国产一区二区三区| 国产精品一区二区在线观看99| 欧美精品一区二区大全| 在线观看免费日韩欧美大片| 精品国产乱码久久久久久小说| 老司机午夜福利在线观看视频 | 俄罗斯特黄特色一大片| 永久免费av网站大全| 久久香蕉激情| 欧美精品av麻豆av| 午夜福利视频精品| 在线十欧美十亚洲十日本专区| 国产亚洲av高清不卡| 伊人久久大香线蕉亚洲五| 国产一区有黄有色的免费视频| 久久久欧美国产精品| 亚洲精品国产精品久久久不卡| 国产伦人伦偷精品视频| 免费看十八禁软件| 国产精品秋霞免费鲁丝片| 人人妻人人澡人人看| 男女边摸边吃奶| 亚洲人成电影免费在线| 在线av久久热| 亚洲精品国产色婷婷电影| 精品国产国语对白av| 美女国产高潮福利片在线看| 999精品在线视频| 丁香六月天网| 国产精品久久久人人做人人爽| 成年人免费黄色播放视频| 永久免费av网站大全| 别揉我奶头~嗯~啊~动态视频 | 另类精品久久| 天天躁狠狠躁夜夜躁狠狠躁| 宅男免费午夜| 人人妻人人澡人人爽人人夜夜| 十八禁人妻一区二区| 亚洲av男天堂| 青青草视频在线视频观看| 国产成+人综合+亚洲专区| 亚洲成人免费电影在线观看| 日韩制服丝袜自拍偷拍| 国产精品久久久人人做人人爽| 成人亚洲精品一区在线观看| 电影成人av| 亚洲国产av影院在线观看| www.熟女人妻精品国产| 国产成人a∨麻豆精品| 欧美日韩av久久| 亚洲av成人不卡在线观看播放网 | av片东京热男人的天堂| 91大片在线观看| 美国免费a级毛片| 久久久久久免费高清国产稀缺| 国产日韩欧美亚洲二区| 久久精品aⅴ一区二区三区四区| 国产亚洲欧美在线一区二区| 国产精品 国内视频| 亚洲欧美成人综合另类久久久| 欧美另类一区| 夜夜骑夜夜射夜夜干| 亚洲欧美精品综合一区二区三区| 国产在线一区二区三区精| 亚洲精品一二三| 日韩欧美免费精品| 国产亚洲欧美精品永久| 亚洲精品久久午夜乱码| 男人爽女人下面视频在线观看| 不卡一级毛片| 亚洲第一青青草原| 久久久国产精品麻豆| 2018国产大陆天天弄谢| 曰老女人黄片| 一区二区日韩欧美中文字幕| 亚洲欧美一区二区三区黑人| 女警被强在线播放| av一本久久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 男女高潮啪啪啪动态图| netflix在线观看网站| 激情视频va一区二区三区| 男女无遮挡免费网站观看| 久久久国产精品麻豆| 免费人妻精品一区二区三区视频| 熟女少妇亚洲综合色aaa.| 精品国产国语对白av| 蜜桃国产av成人99| 久久精品国产亚洲av香蕉五月 | 亚洲欧美日韩高清在线视频 | 女性被躁到高潮视频| 91字幕亚洲| 欧美老熟妇乱子伦牲交| 性色av乱码一区二区三区2| 99热国产这里只有精品6| 久久久久精品国产欧美久久久 | 热re99久久精品国产66热6| www日本在线高清视频| 动漫黄色视频在线观看| 一级a爱视频在线免费观看| 纯流量卡能插随身wifi吗| av欧美777| 国产不卡av网站在线观看| 男女之事视频高清在线观看| 热99re8久久精品国产| 欧美大码av| 亚洲精品一卡2卡三卡4卡5卡 | 成在线人永久免费视频| 国产欧美日韩一区二区精品| 人人澡人人妻人| 精品亚洲成国产av| 超色免费av| 亚洲av成人一区二区三| 青春草视频在线免费观看| 波多野结衣一区麻豆| 少妇猛男粗大的猛烈进出视频| 亚洲精品国产区一区二| 考比视频在线观看| 永久免费av网站大全| 国产97色在线日韩免费| 亚洲成人手机| 好男人电影高清在线观看| 十八禁网站网址无遮挡| 免费人妻精品一区二区三区视频| 亚洲精品国产一区二区精华液| 亚洲av成人一区二区三| 欧美性长视频在线观看| 久久人人97超碰香蕉20202| 深夜精品福利| 日韩大码丰满熟妇| 男女国产视频网站| 亚洲午夜精品一区,二区,三区| 欧美av亚洲av综合av国产av| 久久中文字幕一级| 老司机午夜十八禁免费视频| 在线观看舔阴道视频| videos熟女内射| 中文字幕人妻丝袜制服| 女人精品久久久久毛片| 一区二区三区乱码不卡18| 91成人精品电影| 亚洲欧美一区二区三区黑人| 亚洲欧洲日产国产| 久久久久国产精品人妻一区二区| 国产一卡二卡三卡精品| 国产黄频视频在线观看| 91老司机精品| 亚洲精品中文字幕在线视频| 欧美日本中文国产一区发布| 一区在线观看完整版| 9色porny在线观看| 亚洲精品国产区一区二| 一二三四社区在线视频社区8| 日韩大码丰满熟妇| 中文字幕最新亚洲高清| 亚洲欧美日韩另类电影网站| 亚洲视频免费观看视频| 午夜精品久久久久久毛片777| 亚洲自偷自拍图片 自拍| 精品人妻在线不人妻| 一级片免费观看大全| 免费久久久久久久精品成人欧美视频| 美女国产高潮福利片在线看| 免费少妇av软件| 精品免费久久久久久久清纯 | 国产野战对白在线观看| 国产伦人伦偷精品视频| 自拍欧美九色日韩亚洲蝌蚪91| 老司机亚洲免费影院| 少妇被粗大的猛进出69影院| 午夜福利在线观看吧| 成人影院久久| 亚洲三区欧美一区| 黄片小视频在线播放| 亚洲五月婷婷丁香| 人人澡人人妻人| 亚洲精品中文字幕在线视频| 最近中文字幕2019免费版| 99久久精品国产亚洲精品| 天天躁日日躁夜夜躁夜夜| 国产日韩欧美亚洲二区| 国产av精品麻豆| 久久香蕉激情| 久久 成人 亚洲| 午夜免费成人在线视频| 99国产综合亚洲精品| 狂野欧美激情性xxxx| 一本大道久久a久久精品| 国产精品久久久av美女十八| 亚洲欧美精品自产自拍| 高清av免费在线| 999久久久国产精品视频| 久久久久久久久久久久大奶| 国产一区二区 视频在线| 亚洲欧美激情在线| 99精品欧美一区二区三区四区| 免费观看a级毛片全部| 91精品伊人久久大香线蕉| 久久人人97超碰香蕉20202| 国产男女超爽视频在线观看| 久久久久久亚洲精品国产蜜桃av| 亚洲专区字幕在线| 国产又爽黄色视频| 麻豆av在线久日| 91老司机精品| 中文字幕制服av| 午夜福利免费观看在线| 国产麻豆69| 无限看片的www在线观看| 国产在线观看jvid| 人成视频在线观看免费观看| av一本久久久久| 激情视频va一区二区三区| 久久久久精品人妻al黑| 女警被强在线播放| 在线观看舔阴道视频| 天天影视国产精品| 色精品久久人妻99蜜桃| 国产一区二区 视频在线| 悠悠久久av| 天天添夜夜摸| 亚洲综合色网址| 国产视频一区二区在线看| 99re6热这里在线精品视频| 日韩有码中文字幕| 一本大道久久a久久精品| 在线观看免费视频网站a站| 汤姆久久久久久久影院中文字幕| 青青草视频在线视频观看| 成人国产一区最新在线观看| 中文字幕人妻丝袜制服| 国产亚洲午夜精品一区二区久久| 国产精品一二三区在线看| 岛国毛片在线播放| 日韩大片免费观看网站| 精品福利观看| 一本一本久久a久久精品综合妖精| 日韩有码中文字幕| 我要看黄色一级片免费的| 一区二区日韩欧美中文字幕| 国产色视频综合| 精品少妇内射三级| 国产日韩欧美在线精品| 1024视频免费在线观看| xxxhd国产人妻xxx| 亚洲伊人色综图| 纵有疾风起免费观看全集完整版| 成人18禁高潮啪啪吃奶动态图| 999久久久国产精品视频| 免费女性裸体啪啪无遮挡网站| 搡老乐熟女国产| 午夜福利视频在线观看免费| 成人手机av| 精品福利观看| 精品久久久久久久毛片微露脸 | 色播在线永久视频| 人妻 亚洲 视频| 精品熟女少妇八av免费久了| 精品久久久精品久久久| 国产精品国产av在线观看| 老司机午夜十八禁免费视频| 少妇人妻久久综合中文| 中国国产av一级| h视频一区二区三区| 美女主播在线视频| 免费日韩欧美在线观看| 捣出白浆h1v1| 大片电影免费在线观看免费| 悠悠久久av| cao死你这个sao货| 人人妻人人添人人爽欧美一区卜| 欧美av亚洲av综合av国产av| 男人操女人黄网站| av不卡在线播放| 亚洲av电影在线进入| 国产99久久九九免费精品| 国产精品国产av在线观看| 欧美黄色淫秽网站| 欧美亚洲日本最大视频资源| 久久久久久久久久久久大奶| 亚洲精品国产av成人精品| 日韩三级视频一区二区三区| 制服人妻中文乱码| h视频一区二区三区| 色婷婷av一区二区三区视频| 国产无遮挡羞羞视频在线观看| 电影成人av| 我要看黄色一级片免费的| 亚洲av日韩精品久久久久久密| 成人影院久久| 久久亚洲精品不卡| 18禁裸乳无遮挡动漫免费视频| 一级a爱视频在线免费观看| 国产一区二区三区在线臀色熟女 | av在线播放精品| 欧美激情 高清一区二区三区| 日韩一卡2卡3卡4卡2021年| 亚洲午夜精品一区,二区,三区| 日韩欧美免费精品| 日日摸夜夜添夜夜添小说| 国产欧美日韩精品亚洲av| 精品久久蜜臀av无| 亚洲精品国产区一区二| 男女下面插进去视频免费观看| 午夜福利视频精品| 欧美黑人精品巨大| 精品久久久久久电影网| 一级a爱视频在线免费观看| 国产精品自产拍在线观看55亚洲 | 亚洲激情五月婷婷啪啪| 女人爽到高潮嗷嗷叫在线视频| 亚洲av美国av| 欧美中文综合在线视频| 日本91视频免费播放| 亚洲一卡2卡3卡4卡5卡精品中文| 九色亚洲精品在线播放| 人成视频在线观看免费观看| 老司机亚洲免费影院| 日本a在线网址| 亚洲自偷自拍图片 自拍| 美女脱内裤让男人舔精品视频| 精品一区二区三区av网在线观看 | 大片免费播放器 马上看| 十八禁网站免费在线| 高潮久久久久久久久久久不卡| 久久中文看片网| 亚洲欧美日韩高清在线视频 | 19禁男女啪啪无遮挡网站| 日本av手机在线免费观看| 国产一卡二卡三卡精品| 一级毛片电影观看| 搡老岳熟女国产| 菩萨蛮人人尽说江南好唐韦庄| 在线观看一区二区三区激情| 美女扒开内裤让男人捅视频| 日韩视频在线欧美| 老司机午夜福利在线观看视频 | 黄网站色视频无遮挡免费观看| 中国美女看黄片| 侵犯人妻中文字幕一二三四区| 午夜福利在线观看吧| 又紧又爽又黄一区二区| 麻豆av在线久日| 亚洲精品成人av观看孕妇| 久久精品成人免费网站| 日韩 亚洲 欧美在线| 久久免费观看电影| 亚洲欧美日韩另类电影网站| 麻豆av在线久日| 亚洲精品第二区| 51午夜福利影视在线观看| 亚洲成人手机| 国产三级黄色录像| 黄色a级毛片大全视频| 国产一卡二卡三卡精品| 欧美少妇被猛烈插入视频| 午夜影院在线不卡| 成人国产av品久久久| 少妇猛男粗大的猛烈进出视频| 国产一卡二卡三卡精品| 国产精品免费大片| 亚洲国产看品久久| 亚洲天堂av无毛| av电影中文网址| 中文字幕另类日韩欧美亚洲嫩草| 国产高清国产精品国产三级| 午夜免费成人在线视频| 日本一区二区免费在线视频| 秋霞在线观看毛片| 伊人亚洲综合成人网| a级毛片黄视频| 纵有疾风起免费观看全集完整版| 久久精品亚洲熟妇少妇任你| 国产亚洲精品久久久久5区| 国产免费一区二区三区四区乱码| 免费高清在线观看日韩| 亚洲国产欧美网| 久久久国产精品麻豆| 窝窝影院91人妻| 在线天堂中文资源库| 国产精品一区二区在线不卡| 国产91精品成人一区二区三区 | 日本猛色少妇xxxxx猛交久久| 日本欧美视频一区| 老司机福利观看| 久久亚洲精品不卡| 亚洲国产中文字幕在线视频| 久久亚洲精品不卡| 久久久久国产一级毛片高清牌| av在线app专区| 99热网站在线观看| av线在线观看网站| 菩萨蛮人人尽说江南好唐韦庄| 男人舔女人的私密视频| 别揉我奶头~嗯~啊~动态视频 | 满18在线观看网站| 人妻人人澡人人爽人人| 国产成人精品在线电影| 午夜日韩欧美国产| 黄色视频,在线免费观看| 精品少妇内射三级| 精品久久久久久久毛片微露脸 | 天堂俺去俺来也www色官网| 日韩大码丰满熟妇| 建设人人有责人人尽责人人享有的| 69av精品久久久久久 | 91精品国产国语对白视频| 岛国毛片在线播放| 免费一级毛片在线播放高清视频 | 亚洲 欧美一区二区三区| 国产一区二区三区综合在线观看| 黄色a级毛片大全视频| 中文字幕人妻丝袜一区二区| 黑人欧美特级aaaaaa片| 亚洲成av片中文字幕在线观看| 黄片播放在线免费| 亚洲男人天堂网一区| 成人18禁高潮啪啪吃奶动态图| 日本精品一区二区三区蜜桃| 两个人看的免费小视频| 精品少妇黑人巨大在线播放| 啦啦啦中文免费视频观看日本| 美女高潮到喷水免费观看| 午夜日韩欧美国产| 中文字幕另类日韩欧美亚洲嫩草| 两性午夜刺激爽爽歪歪视频在线观看 | 国产麻豆69|