• <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ī)院感染嗜水氣單胞菌的臨床治療分析
    午夜福利视频1000在线观看| 免费在线观看亚洲国产| 丁香六月欧美| 99热只有精品国产| 日日摸夜夜添夜夜添小说| 一个人免费在线观看电影 | 国产亚洲欧美在线一区二区| 欧美在线一区亚洲| 午夜福利免费观看在线| 国产探花在线观看一区二区| 国产一区二区在线av高清观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲一区二区三区色噜噜| 亚洲人成网站在线播放欧美日韩| 国模一区二区三区四区视频 | 亚洲片人在线观看| 网址你懂的国产日韩在线| 国产亚洲精品久久久久久毛片| 欧美日韩综合久久久久久 | 久久久久久九九精品二区国产| 亚洲国产精品久久男人天堂| 国产欧美日韩一区二区精品| 99久久综合精品五月天人人| 国产乱人伦免费视频| 亚洲成人中文字幕在线播放| 色综合亚洲欧美另类图片| АⅤ资源中文在线天堂| 国产一区二区在线观看日韩 | 少妇丰满av| 麻豆成人av在线观看| 757午夜福利合集在线观看| 久久欧美精品欧美久久欧美| 欧美av亚洲av综合av国产av| 欧美日韩综合久久久久久 | 精品久久久久久,| 999久久久国产精品视频| 亚洲av成人av| 久久久久性生活片| 中文字幕av在线有码专区| 国产精品电影一区二区三区| 麻豆av在线久日| 午夜视频精品福利| 国产成人精品久久二区二区免费| 国产亚洲精品综合一区在线观看| 婷婷亚洲欧美| 国产成+人综合+亚洲专区| 久久久色成人| 亚洲国产欧洲综合997久久,| 亚洲成人精品中文字幕电影| 男人舔奶头视频| 97超视频在线观看视频| 18美女黄网站色大片免费观看| 精品久久蜜臀av无| 九色国产91popny在线| 99久国产av精品| 嫩草影院精品99| 99久久精品热视频| 亚洲真实伦在线观看| 成年女人看的毛片在线观看| 丰满的人妻完整版| 久久精品91无色码中文字幕| 99热这里只有是精品50| 少妇熟女aⅴ在线视频| 黑人欧美特级aaaaaa片| 日本一二三区视频观看| 亚洲欧美一区二区三区黑人| 久久精品综合一区二区三区| 亚洲欧美日韩卡通动漫| 91久久精品国产一区二区成人 | 一卡2卡三卡四卡精品乱码亚洲| 亚洲aⅴ乱码一区二区在线播放| 免费大片18禁| 亚洲一区高清亚洲精品| 亚洲午夜精品一区,二区,三区| 亚洲国产精品999在线| 综合色av麻豆| 国产成人福利小说| 亚洲精品一卡2卡三卡4卡5卡| 色综合婷婷激情| 久久精品国产综合久久久| 不卡av一区二区三区| 欧美精品啪啪一区二区三区| 国产精品一区二区三区四区免费观看 | 日本在线视频免费播放| 一本久久中文字幕| 亚洲黑人精品在线| 国产精品自产拍在线观看55亚洲| 国产精品,欧美在线| 国产精品99久久久久久久久| 欧美不卡视频在线免费观看| 国产极品精品免费视频能看的| 午夜福利欧美成人| 舔av片在线| 国产美女午夜福利| 99热精品在线国产| 久久国产精品人妻蜜桃| 免费观看精品视频网站| 国产亚洲精品av在线| 成年版毛片免费区| 999精品在线视频| 女人高潮潮喷娇喘18禁视频| 国产精品一区二区三区四区久久| 欧美日韩中文字幕国产精品一区二区三区| 免费大片18禁| 精品熟女少妇八av免费久了| 久99久视频精品免费| 精品日产1卡2卡| 日韩高清综合在线| 一个人免费在线观看的高清视频| 国产精品女同一区二区软件 | 欧美成狂野欧美在线观看| 久久久久国产一级毛片高清牌| xxxwww97欧美| 麻豆成人午夜福利视频| 亚洲色图 男人天堂 中文字幕| 一区二区三区激情视频| 黑人巨大精品欧美一区二区mp4| 9191精品国产免费久久| 午夜日韩欧美国产| 校园春色视频在线观看| 午夜福利视频1000在线观看| 国产淫片久久久久久久久 | 亚洲专区中文字幕在线| 亚洲激情在线av| 亚洲国产看品久久| 搡老妇女老女人老熟妇| 又爽又黄无遮挡网站| 国产伦在线观看视频一区| 在线播放国产精品三级| 亚洲一区二区三区不卡视频| 国产亚洲精品一区二区www| 国产熟女xx| 黄色 视频免费看| 波多野结衣高清无吗| 热99在线观看视频| 少妇裸体淫交视频免费看高清| 久久九九热精品免费| 亚洲人成电影免费在线| 国产成人欧美在线观看| 欧美中文综合在线视频| 一级毛片女人18水好多| 亚洲av成人一区二区三| 国产极品精品免费视频能看的| 久久久水蜜桃国产精品网| 久久香蕉国产精品| 亚洲精华国产精华精| 国产成人影院久久av| 久久草成人影院| 久久久久久久久免费视频了| 亚洲va日本ⅴa欧美va伊人久久| 两个人的视频大全免费| 色在线成人网| 99国产综合亚洲精品| 亚洲五月婷婷丁香| 美女高潮的动态| www.精华液| 美女扒开内裤让男人捅视频| 真人做人爱边吃奶动态| 欧美日韩综合久久久久久 | 欧美日本视频| 国内少妇人妻偷人精品xxx网站 | 国产久久久一区二区三区| 国产亚洲精品av在线| 大型黄色视频在线免费观看| 久久久成人免费电影| 亚洲欧美日韩卡通动漫| 国内精品久久久久久久电影| 人人妻人人看人人澡| 成年女人看的毛片在线观看| 国产精品一及| 国产精品自产拍在线观看55亚洲| 亚洲七黄色美女视频| e午夜精品久久久久久久| 高潮久久久久久久久久久不卡| 亚洲va日本ⅴa欧美va伊人久久| 色在线成人网| 免费高清视频大片| 综合色av麻豆| 99精品久久久久人妻精品| 热99re8久久精品国产| 亚洲天堂国产精品一区在线| 两性午夜刺激爽爽歪歪视频在线观看| 国产乱人视频| 欧美+亚洲+日韩+国产| 欧美午夜高清在线| 欧美性猛交╳xxx乱大交人| 国产成人aa在线观看| 日日摸夜夜添夜夜添小说| 国产精品99久久99久久久不卡| 99re在线观看精品视频| 综合色av麻豆| 成年人黄色毛片网站| 99精品久久久久人妻精品| 最新美女视频免费是黄的| 999久久久精品免费观看国产| 精品国产美女av久久久久小说| 九九久久精品国产亚洲av麻豆 | 久9热在线精品视频| 天堂网av新在线| 一进一出抽搐gif免费好疼| 亚洲av成人av| 99精品久久久久人妻精品| 国产 一区 欧美 日韩| 亚洲欧美日韩高清在线视频| 啦啦啦韩国在线观看视频| 国产高清videossex| 男女那种视频在线观看| 亚洲欧美日韩高清在线视频| 亚洲av中文字字幕乱码综合| 91九色精品人成在线观看| 不卡av一区二区三区| 亚洲欧美精品综合一区二区三区| 亚洲人成网站在线播放欧美日韩| 日本黄大片高清| 嫩草影视91久久| 国产野战对白在线观看| 日韩 欧美 亚洲 中文字幕| 他把我摸到了高潮在线观看| 亚洲熟妇中文字幕五十中出| 久久久久久久午夜电影| 国产成人精品久久二区二区91| 久久精品91无色码中文字幕| 久久久久九九精品影院| 欧美色视频一区免费| 亚洲人成伊人成综合网2020| 国产精品亚洲av一区麻豆| 国产aⅴ精品一区二区三区波| 99精品在免费线老司机午夜| 日日摸夜夜添夜夜添小说| 成人永久免费在线观看视频| 曰老女人黄片| 国产午夜福利久久久久久| 日韩精品青青久久久久久| 亚洲乱码一区二区免费版| 一区福利在线观看| 黄片小视频在线播放| 欧美日韩中文字幕国产精品一区二区三区| 91av网站免费观看| 99久久国产精品久久久| 99久久成人亚洲精品观看| 观看免费一级毛片| 国产伦精品一区二区三区四那| 欧美高清成人免费视频www| 又黄又爽又免费观看的视频| 亚洲人成网站在线播放欧美日韩| 午夜视频精品福利| 国产成人系列免费观看| 国产91精品成人一区二区三区| 亚洲精品国产精品久久久不卡| av在线天堂中文字幕| av视频在线观看入口| 久久精品91无色码中文字幕| 久9热在线精品视频| 免费看日本二区| 久久国产精品人妻蜜桃| 人人妻人人看人人澡| 国语自产精品视频在线第100页| 国产高清视频在线播放一区| 在线观看一区二区三区| 五月玫瑰六月丁香| 免费观看人在逋| 美女cb高潮喷水在线观看 | 变态另类成人亚洲欧美熟女| 免费在线观看影片大全网站| 亚洲真实伦在线观看| 香蕉av资源在线| or卡值多少钱| 亚洲黑人精品在线| 亚洲成人久久爱视频| 99热这里只有精品一区 | 欧美日韩综合久久久久久 | 99久久精品一区二区三区| 又黄又爽又免费观看的视频| 人妻夜夜爽99麻豆av| 精品久久久久久,| 99精品欧美一区二区三区四区| 免费av毛片视频| 国产精品一及| 日韩欧美三级三区| 99精品欧美一区二区三区四区| 欧美日韩乱码在线| 免费观看的影片在线观看| 欧美中文日本在线观看视频| 日本一本二区三区精品| 亚洲自偷自拍图片 自拍| 亚洲人成网站在线播放欧美日韩| 人人妻人人澡欧美一区二区| 午夜免费观看网址| 香蕉丝袜av| 真人一进一出gif抽搐免费| 婷婷丁香在线五月| 免费观看精品视频网站| 国产亚洲av高清不卡| 国产精品一及| 日本免费a在线| 精品一区二区三区视频在线 | 夜夜夜夜夜久久久久| 久久人妻av系列| 国产精品乱码一区二三区的特点| 老汉色∧v一级毛片| 少妇丰满av| 最近视频中文字幕2019在线8| 老鸭窝网址在线观看| 两人在一起打扑克的视频| 久久久精品大字幕| 又黄又爽又免费观看的视频| 日韩精品青青久久久久久| 日韩欧美三级三区| 亚洲av熟女| 亚洲欧美激情综合另类| 亚洲av成人精品一区久久| 少妇熟女aⅴ在线视频| 午夜福利欧美成人| 90打野战视频偷拍视频| 午夜免费激情av| 欧美成人免费av一区二区三区| 一级毛片精品| 日本成人三级电影网站| 在线视频色国产色| 久久久久久大精品| 给我免费播放毛片高清在线观看| 成熟少妇高潮喷水视频| 看黄色毛片网站| 999精品在线视频| 他把我摸到了高潮在线观看| 黄频高清免费视频| x7x7x7水蜜桃| 看免费av毛片| 亚洲人成网站在线播放欧美日韩| 亚洲成人中文字幕在线播放| 九九热线精品视视频播放| 18禁裸乳无遮挡免费网站照片| 欧美黄色淫秽网站| 国产精品98久久久久久宅男小说| 神马国产精品三级电影在线观看| www.999成人在线观看| a级毛片在线看网站| 亚洲第一欧美日韩一区二区三区| 亚洲av成人一区二区三| 99久久久亚洲精品蜜臀av| 美女黄网站色视频| 亚洲五月天丁香| 国产单亲对白刺激| 亚洲成av人片在线播放无| 亚洲精品在线美女| 男女做爰动态图高潮gif福利片| 性色av乱码一区二区三区2| 国产成人av教育| 国产激情欧美一区二区| 亚洲成人免费电影在线观看| 久久精品国产亚洲av香蕉五月| 久久中文字幕一级| 无人区码免费观看不卡| 国产精品一区二区免费欧美| 全区人妻精品视频| 国产精品香港三级国产av潘金莲| 亚洲九九香蕉| 国产日本99.免费观看| 国产激情偷乱视频一区二区| 亚洲一区二区三区色噜噜| 美女黄网站色视频| 亚洲精品久久国产高清桃花| 国产成人啪精品午夜网站| 国产伦一二天堂av在线观看| 此物有八面人人有两片| 国产午夜精品久久久久久| 午夜精品在线福利| 亚洲天堂国产精品一区在线| 欧美日韩精品网址| 国产真人三级小视频在线观看| 免费大片18禁| 国产69精品久久久久777片 | 天堂动漫精品| 国产亚洲av嫩草精品影院| 久久国产精品人妻蜜桃| 午夜福利在线观看免费完整高清在 | 中文字幕av在线有码专区| 亚洲一区二区三区色噜噜| 国产成年人精品一区二区| 香蕉国产在线看| 中文字幕久久专区| 亚洲成人久久性| 99在线人妻在线中文字幕| 亚洲国产色片| 日韩欧美国产在线观看| 欧美日韩精品网址| av国产免费在线观看| 免费大片18禁| 免费在线观看成人毛片| 看免费av毛片| 在线观看66精品国产| 午夜福利18| 亚洲精品粉嫩美女一区| 国产精品一及| av视频在线观看入口| 国产亚洲欧美98| 国产成人精品久久二区二区免费| 精品国产超薄肉色丝袜足j| 婷婷六月久久综合丁香| 亚洲精品中文字幕一二三四区| 非洲黑人性xxxx精品又粗又长| 日日夜夜操网爽| 精品国产超薄肉色丝袜足j| 99热只有精品国产| 久久久国产成人免费| 窝窝影院91人妻| 国产三级中文精品| 久久中文字幕一级| 婷婷丁香在线五月| 国产精品1区2区在线观看.| 久久天堂一区二区三区四区| 俺也久久电影网| 亚洲欧美日韩高清在线视频| 日本免费一区二区三区高清不卡| 舔av片在线| 成年女人毛片免费观看观看9| 免费在线观看视频国产中文字幕亚洲| 巨乳人妻的诱惑在线观看| 亚洲黑人精品在线| 欧美大码av| 国产一区二区三区视频了| 亚洲中文av在线| 亚洲一区高清亚洲精品| 曰老女人黄片| 97超视频在线观看视频| 别揉我奶头~嗯~啊~动态视频| 亚洲中文日韩欧美视频| 亚洲欧美日韩无卡精品| 亚洲黑人精品在线| 真实男女啪啪啪动态图| 精品电影一区二区在线| 91av网站免费观看| 午夜成年电影在线免费观看| 全区人妻精品视频| 久久久久久国产a免费观看| 夜夜夜夜夜久久久久| 国产精品一区二区三区四区免费观看 | 欧美日韩乱码在线| 午夜成年电影在线免费观看| 亚洲精品在线美女| 久久久国产精品麻豆| 99久久综合精品五月天人人| 国产精品亚洲一级av第二区| 变态另类丝袜制服| 美女高潮的动态| 国产精品99久久99久久久不卡| 五月玫瑰六月丁香| 真人一进一出gif抽搐免费| 色综合婷婷激情| 最新中文字幕久久久久 | 国产激情久久老熟女| 精品一区二区三区视频在线观看免费| 国产亚洲欧美在线一区二区| 亚洲一区高清亚洲精品| 亚洲专区国产一区二区| 国产精品久久久人人做人人爽| 久久久国产成人免费| 欧美日韩国产亚洲二区| 99久久精品热视频| 国产美女午夜福利| 日韩精品中文字幕看吧| 日本精品一区二区三区蜜桃| 亚洲欧美激情综合另类| 国产精品一区二区精品视频观看| 久久久精品欧美日韩精品| 欧美中文综合在线视频| 每晚都被弄得嗷嗷叫到高潮| 免费观看人在逋| 九色国产91popny在线| 亚洲精品在线观看二区| 啪啪无遮挡十八禁网站| 怎么达到女性高潮| 在线视频色国产色| 久久性视频一级片| 欧美日本视频| 一进一出好大好爽视频| 中文字幕精品亚洲无线码一区| 激情在线观看视频在线高清| 老汉色av国产亚洲站长工具| www国产在线视频色| 夜夜看夜夜爽夜夜摸| 国产一区二区三区在线臀色熟女| 久久伊人香网站| 男人舔女人的私密视频| 免费搜索国产男女视频| 麻豆一二三区av精品| 亚洲国产看品久久| 在线观看免费午夜福利视频| 白带黄色成豆腐渣| 视频区欧美日本亚洲| 日本在线视频免费播放| 男人舔奶头视频| 欧美xxxx黑人xx丫x性爽| 一区二区三区国产精品乱码| 亚洲片人在线观看| 国产麻豆成人av免费视频| 黄色成人免费大全| 最近视频中文字幕2019在线8| 国产人伦9x9x在线观看| 免费看a级黄色片| 欧美乱妇无乱码| 又黄又爽又免费观看的视频| 国产精品九九99| 国产久久久一区二区三区| 久久久久亚洲av毛片大全| 国模一区二区三区四区视频 | 九色成人免费人妻av| 亚洲国产中文字幕在线视频| 亚洲专区中文字幕在线| 久久亚洲真实| 亚洲国产欧美网| 黄片小视频在线播放| 免费看日本二区| 国产爱豆传媒在线观看| 欧美xxxx黑人xx丫x性爽| 亚洲国产日韩欧美精品在线观看 | 亚洲国产精品999在线| 在线观看日韩欧美| av天堂中文字幕网| 偷拍熟女少妇极品色| 精品国内亚洲2022精品成人| 久久精品夜夜夜夜夜久久蜜豆| 精品无人区乱码1区二区| 国产亚洲精品久久久久久毛片| 搡老妇女老女人老熟妇| 日日摸夜夜添夜夜添小说| 香蕉av资源在线| 特大巨黑吊av在线直播| 日韩精品青青久久久久久| 给我免费播放毛片高清在线观看| 国产毛片a区久久久久| 国产精品自产拍在线观看55亚洲| 高潮久久久久久久久久久不卡| 久久亚洲精品不卡| 欧美丝袜亚洲另类 | 亚洲国产高清在线一区二区三| 亚洲精品中文字幕一二三四区| 中文字幕高清在线视频| 国产高清视频在线观看网站| 精品99又大又爽又粗少妇毛片 | 国产av在哪里看| 麻豆成人午夜福利视频| 国产视频一区二区在线看| 18禁美女被吸乳视频| 亚洲国产欧美网| 国产av不卡久久| 观看免费一级毛片| 亚洲片人在线观看| 精品久久蜜臀av无| 国产av在哪里看| 亚洲人与动物交配视频| 999久久久精品免费观看国产| 亚洲精华国产精华精| www日本在线高清视频| 99国产综合亚洲精品| 精品久久久久久久人妻蜜臀av| 中文字幕人妻丝袜一区二区| 在线播放国产精品三级| 在线观看一区二区三区| 国产成人欧美在线观看| 99国产极品粉嫩在线观看| 18禁黄网站禁片午夜丰满| 老汉色∧v一级毛片| 美女cb高潮喷水在线观看 | 国产精品女同一区二区软件 | 两个人的视频大全免费| 国产亚洲精品久久久com| 一级黄色大片毛片| 色播亚洲综合网| 国内精品一区二区在线观看| 观看美女的网站| 国产在线精品亚洲第一网站| 91麻豆精品激情在线观看国产| www.精华液| 午夜a级毛片| 亚洲avbb在线观看| 欧美在线黄色| 特大巨黑吊av在线直播| 亚洲成人中文字幕在线播放| 女人被狂操c到高潮| 午夜影院日韩av| 成年免费大片在线观看| 国产男靠女视频免费网站| 三级毛片av免费| 国产三级中文精品| 黄色片一级片一级黄色片| 欧美一区二区精品小视频在线| 久久久国产成人免费| 日韩欧美 国产精品| 99精品久久久久人妻精品| 国产高清视频在线播放一区| 午夜精品久久久久久毛片777| 999精品在线视频| 欧美色视频一区免费| 12—13女人毛片做爰片一| 欧美+亚洲+日韩+国产| 精品久久蜜臀av无| 男女做爰动态图高潮gif福利片| 精品久久久久久成人av| 欧美黑人欧美精品刺激| 亚洲欧美日韩高清在线视频| 国产精品自产拍在线观看55亚洲| 亚洲最大成人中文| 一本久久中文字幕| 老汉色∧v一级毛片| 国产精品美女特级片免费视频播放器 | 亚洲精品色激情综合| 啦啦啦韩国在线观看视频| 视频区欧美日本亚洲| 国产伦一二天堂av在线观看| 91在线精品国自产拍蜜月 | 亚洲欧美精品综合一区二区三区| 日韩高清综合在线| 国产高潮美女av| 国产视频内射|