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

    規(guī)整填料塔內氣-液兩相并流流動的CFD研究

    2017-04-09 06:55:27高騰飛羅文媛韓少強
    化學工業(yè)與工程 2017年1期
    關鍵詞:規(guī)整液量氣液

    何 杰,高騰飛,羅文媛,韓少強,王 昕

    (天津大學化工學院,天津 300072)

    規(guī)整填料塔目前已經廣泛應用與石油化工領域。因為它具有壓降低、通量大、效率高,而且持液量小的優(yōu)勢。同時,規(guī)整填料塔具有結構規(guī)整的特點,與板式塔相比,通過合理巧妙的設計,能使因工業(yè)化而產生的放大效應降到幾乎為0[1],但是由于技術水平的限制,學術界對規(guī)整填料塔內流體的流動機理尚無定論,由于機理的不明確,限制了其進一步的發(fā)展和更新。對規(guī)整填料機理的研究具有現(xiàn)實的難度,因為這種填料在結構上有規(guī)整性與復雜性并存的特點,即存在各項異性,這是結構上區(qū)別于散裝填料的顯著特點。隨著現(xiàn)代計算機理論和算法的進步,電算法解決工程問題已經越來越普遍。這些進步為CFD研究填料塔內的流體力學性能提供了技術支持。谷芳[2]對規(guī)整填料板表面上的液體降膜流動過程建立了三維模擬,對這一實際過程進行了研究,此法是通過與二維模型的對比,并在此基礎上構建波紋板和平板的三維降膜兩相流模型的基礎之上進行的。陳江波[3-4]提出的三維兩相流動模型中,對氣液相間作用力源項(FLG)和液相表面張力源項(Fvol)對結果的影響給予重視。在他的工作中,明確闡釋了規(guī)整填料特征結構單元內的潤濕面積、持液量和該單元結構之分離效率于不同操作條件下的變化,并得到了氣液相的速度和各個組分的濃度分布圖,以直觀的方式表現(xiàn)出氣液相在模型表面的運動狀態(tài)。值得注意的是,以上研究均是以規(guī)整填料內氣液兩相逆流流動現(xiàn)象做為研究背景。研究流體并流流動在當今具有非?,F(xiàn)實的意義,例如:汽液噴射反應器及環(huán)流反應器中氣液兩相垂直向下或垂直向上的并流流動;真空蒸發(fā)和真空干燥系統(tǒng)用于排除不凝性氣體并制造真空的水力噴射器的大氣腿中氣液兩相的垂直向下并流流動;管內氣液兩相混合的傾斜并流流動等。隨著經濟發(fā)展的需要以及綠色化工概念的指導,如何能節(jié)省設備投資,降低設備占地面積已變得非常重要。并流流動能有效解決上述問題,但遺憾的是,與并流流動相關的理論研究目前非常之少,特別是并流過程中的傳質過程和流體力學行為。

    為了填補氣液兩相在規(guī)整填料中并流流動時機理研究領域的空白。本論文通過構建三維模型,應用流體力學的計算方法對簡化的Mellapak350Y規(guī)整填料內的并流現(xiàn)象進行了研究,并通過相關實驗,以驗證研究成果。

    1 物理模型的建立

    受到當今技術條件的限制,并不能實現(xiàn)直接利用CFD方法模擬規(guī)整填料塔內氣液兩相流動的過程。主要是因為規(guī)整填料內部復雜的幾何結構,使填料內部出現(xiàn)溝流或者溪流等比較復雜的流動形式,而不僅僅是均勻降膜流動。其次,當今計算機的算法和運算速率還不能滿足這一復雜計算的需求。為解決這些不足,常用的方法是將規(guī)整填料床層的內部結構根據(jù)流體在其表面的傳質特點分成4類特征單元,并對不同類型的單元分別進行CFD計算。通過比較,發(fā)現(xiàn)交叉特征單元,即第二特征單元在填料層中占比高達70%之上[5],與其他3種相比,比例最高。為更好的描述這一單元,本論文在其基礎上構建了并流流動的CFD三維模型,如圖1。其中,把方向相反的2條三角形流道也包含在模擬區(qū)域內(以Mellapak350Y型號的板波紋規(guī)整填料的波紋為基礎模型)。

    為加快計算機對模型的處理速度,在可控的模擬范圍內對模型的物理結構稍作調整,假設其表面沒有開孔,而且不加紋。稱之為簡化的Mellapak填料。

    1.1 邊界條件

    模擬規(guī)整填料內氣液兩相并流時的流動過程,對氣液兩相,共有2個進口和2個出口(見圖1)。

    圖1 規(guī)整填料三維近似物理模型的建立Fig.1 Construction of three-dimensional physical model of structured packing

    1.2 網(wǎng)格劃分方法

    在計算之前,對模型進行網(wǎng)格劃分。由于本論文中計算模型的幾何形狀非常復雜,因此采用了非結構化四面體網(wǎng)格進行分網(wǎng)。所劃分的網(wǎng)格總數(shù)為29 474個,如圖2。

    圖2 計算模型的網(wǎng)格Fig.2 Grid of computational model

    有文獻報道[33-34]將規(guī)整填料的結構切割為4類特征單元(REU,representative elementary unit),如圖3所示,對各單元分別進行CFD模擬,而在這幾類特征單元中,第2類單元,即交叉特征單元在填料中所占的份額最多。因此,本論文將采用此結構單元對氣液兩相的并流流動進行近似描述。

    表1 邊界條件

    圖3 板波紋規(guī)整填料的四種典型的結構單元Fig.3 Typical representative elementary units approximating corrugated-sheet structured packing

    2 數(shù)學模型的建立

    2.1 數(shù)學模型

    本論文采用VOF模型,即一種追蹤氣液相界面模型,來模擬本研究體系中氣液兩相的非穩(wěn)態(tài)分層流動現(xiàn)象。經試驗證明,若利用常用的拉格朗日離散相模型或歐拉連續(xù)介質模型描述分散相都難以很好地描述整體的流動過程,究其本質,主要因為氣液兩相之間其實存在非常復雜的相互作用。

    2.2 控制方程

    假設氣液兩相均是牛頓流體,液相不可壓縮,氣相是理想氣體。當氣液兩相處于1個單元結構時,可認為兩相恒溫,此時,混合體系的物理性質僅為組成的函數(shù),控制方程分列于下。

    質量守恒方程:

    (1)

    動量守恒方程:

    通常,對于流體的湍流行為常用標準的k-ε模型,本體系中氣液均處于湍流流動狀態(tài)。通過Raynal等[6]對不同湍流模型的模擬效果的比較,發(fā)現(xiàn)描述規(guī)整填料內流體的流動最合理的模型是RNG(Renormalization group)k-ε模型,所以本論文利用RNGk-ε模型來描述該體系的湍流行為。

    湍動能k的輸運方程為:

    (3)

    耗散率ε的輸運方程為:

    方程中所涉及的流體物性參數(shù)(黏度μ和密度ρ)可以與各相的體積分率建立函數(shù)關系??赏ㄟ^公式(5)和(6)求得黏度和密度:

    ρ=αLρL+(1-αL)ρG(5)

    μ=αLμL+(1-αL)μG(6)

    αL為計算單元內液相的體積分率。

    在VOF法中,通過求解包含相體積分率的連續(xù)方程來完成氣液兩相的相界面的計算。本論文中,液相(分散相)體積分率的連續(xù)性方程為:

    ▽αL=0(7)

    氣相(主相)的體積分率由式(8)得到。

    αG+αL=1(8)

    2.3 動量源項的確定

    氣液兩相在填料塔中的流動有其特殊性,根據(jù)其特征,本論文主要關注了以下2個動力源項的作用,分別是氣液相間作用力動量源項FLG和表面張力動量源項Fvol,其中Fi為各方向的動量源項。方程如下:

    Fi=Fvol,i+FLG,i(9)

    表面張力在本體系中的影響非常重要,這主要是根據(jù)液相在填料塔內的特有流動形式決定的。液膜和液滴是液相在填料塔內的主要存在形式,其自上而下在塔內運動。根據(jù)表面張力的知識,液滴大小和液膜于填料表面的鋪展程度都受到該力的影響,所以,模型構建的過程中,必須考慮表面張力的影響。

    根據(jù)流程模擬系統(tǒng)Fluent的設置,可利用CSF模型對表面張力源項進行模擬(系統(tǒng)默認),其推導過程參見文獻[7]。

    最終,表面張力源項的表達形式如式(10)。

    (10)

    2.4 氣液相間作用力動量源項

    對氣液兩項之間作用的影響也構建了模型。根據(jù)谷芳[2]的研究,填料塔內液相表面很容易因為即便是微小的氣速而影響其表面的波動程度,這種波動對氣液相間的接觸面積影響較大,進而影響規(guī)整填料內的總傳質擴散系數(shù)(KOGae或KOLa),最終影響傳質速率。

    本論文采用Hewitt[8]提出的氣液相間作用力模型,其具體表達式為:

    (11)

    式(11)中,fLG為氣液相界面處的摩擦系數(shù),本論文的摩擦系數(shù)采用Riazi提出的模型[9]。

    (12)

    方程(12)描述的氣液相間作用力源項只是存在于氣液相界面周圍的計算單元。這不同于氣液單項流動區(qū)域計算單元的計算方式,因為在該計算區(qū)域并沒有相間作用力源項。氣液相間作用力模型在軟件中一般首先需要編寫用戶自定義函數(shù)UDF,然后從Fluent的主控制面板中導入,這和表面張力源項的導入方式有所不同。

    值得一提的是,在模型中加入相間作用力之前必須判斷該計算單元內的相含率,這主要是為了確定研究的計算單元是否處于相界面周圍。因為相間相互作用力的本質(兩相界面的摩擦力)告訴我們,該力一般出現(xiàn)在兩相界面附近,可是VOF模型中變量及屬性均為體積平均值,且只有1個動量方程。若計算單元在液相中,F(xiàn)D為正,否則為負。

    3 計算結果與討論

    對填料塔內并流流動持液量和壓降的變化情況作了研究。在該三維模型中,氣液相并流流動,從填料上部的氣相口進入。

    3.1 填料表面液相的流動受氣相Fv因子的影響

    不同F(xiàn)v因子[噴淋密度設定為12.286 m3·(m2·h)-1]時,填料表面的液體含有率分布如圖4所示(此時規(guī)整填料單元的流動基本達到穩(wěn)定)。

    如圖4所示,由于重力和表面張力對液相的影響不同,相比而言,前者作用力更大。這使得當液相在填料上表面形成連續(xù)液膜時,即氣相Fv因子較小時,液膜下表面要寬于上表面。隨著氣相Fv因子的逐漸增加,液相的存在形式也隨著逐漸發(fā)生變化。此時,在填料的表面,液相的存在形式是液膜和液滴。在Fv因子較小時,此時氣相對液相的推動力比較小,所以流動相能形成連續(xù)的液膜;隨著氣速增加,氣相的Fv因子增加到5及以上時,液相的存在形式有液膜轉變?yōu)橐旱?,流動相由液相為主轉變?yōu)闅庀嗪鸵合嗟幕旌衔铩kS著氣相的Fv因子持續(xù)提高,氣體對液膜的推動力也將持續(xù)增加,最終,液相在填料的表面將很難形成連續(xù)膜。同時可以發(fā)現(xiàn),流動形式主要以溪流或者溝流為主,且其在下表面的寬度均寬于上表面。

    圖4 填料表面液體分布受Fv因子的影響Fig.4 Liquid distribution on the surface of the structured packing at the different Fv factors

    3.2 不同氣相Fv因子時,壓降隨液相噴淋密度的變化關系

    研究了并流流動時填料層內的壓降隨氣液相流速的變化情況,并將模擬值與實驗值進行了對比。實驗是以空氣和水作為測試物系,在常壓和常溫下進行的并流操作,實驗中的測試填料為簡化的Mellapak350Y型金屬板波紋填料,實驗流程如圖5所示。模型選取了5個液相噴淋密度在4個氣相Fv因子下的壓降變化情況進行了比較驗證。

    CFD中的設置列于表2。

    表2 CFD中邊界條件設定

    當氣液兩相并流流動,且氣相Fv因子不同時,模型內壓降和液相噴淋密度的關系見于圖6。

    從圖6中可以看出,液相中噴淋密度與壓降成正相關關系,且在不同氣相Fv因子時皆成立,從實驗結果上看,模擬值與實驗室在20%的誤差范圍之內,吻合效果較好。出現(xiàn)此現(xiàn)象主要是因為填料內的積液對氣相流動的阻力比噴淋密度增大對氣相推動的影響要大的多,使得積液的影響占據(jù)主要地位,最終導致了壓降的增大;值得一提的是,在小氣相Fv因子下,壓降的變化隨噴淋密度增大比較反常。其機理尚不明確,原因可能與大氣相Fv因子時的流動機理正好相反,即液相流動對氣相的推動力要比填料內液相的積液阻力大,導致壓降下降,此時液相對氣相的流動推動力不能忽視。由圖6可以看出,同一噴淋密度條件下,壓降與氣相Fv因子成正相關關系。與逆流流動的壓降相比,順流流動具有較大的優(yōu)勢,其壓降值可以降低70%,具體數(shù)據(jù)可以通過與文獻值[10]的比較得到。

    3.3 不同氣相Fv因子下,持液量隨液相噴淋密度的變化

    研究了并流流動時填料層內的持液量隨氣液相流速的變化情況,模型選取了5個液相噴淋密度在4個氣相Fv因子下的壓降變化情況,并將模擬值與實驗值進行了對比。邊界條件設定如表2。

    圖5 試驗流程圖Fig.5 Experimental procedure chart

    圖6 壓降與液相噴淋密度受Fv因子的影響Fig.6 Relationship between pressure drop and liquid spray density at the different Fv factors

    當氣液兩相并流流動,且氣相Fv因子不同時,模型內持液量和液相噴淋密度的關系見于圖7。

    從圖7中可以看出,體系的持液量與液相的噴淋密度呈現(xiàn)正相關關系,且在不同的氣相Fv因子時都呈現(xiàn)相同的規(guī)律。出現(xiàn)這種情況的內在原因是噴淋密度的增大使得填料層的積液增加,但前提是氣相的速度一定。由于為該物理過程構建的計算機模型前期經過簡化,且該試驗的宏觀因素可能造成一定的不確定性,使得模擬結果與實驗值存在一定的允許誤差,但是在實驗值與計算值在變化趨勢上式一致的。如圖7所示,由于氣相流動對液相有推動作用,當氣相Fv因子的增大時,即氣相流量增大時,填料層的液相持液量明顯降低(在同一噴淋密度條件下)。這一現(xiàn)象明顯區(qū)別于逆流流動過程,在逆流流動過程中,如果體系的操作條件在載點以下,填料層的持液量隨氣相Fv因子增大變化較小[11],并流時,不會出現(xiàn)這一現(xiàn)象。且與逆流流動相比[12],并流過程的填料層持液量下降達50%之多。

    圖7 持液量與液相噴淋密度受Fv因子的影響Fig.7 Relationship between liquid holdup and liquid spray density at the different Fv factors

    4 結論

    研究了板波紋規(guī)整填料內部氣液兩相的并流流動現(xiàn)象。通過對實體結構內部的簡化,用流體力學軟件構建了簡化的Mellapak350Y規(guī)整填料的三維立體模型。通過該模型,對氣液并流流動時的持液量和壓降進行了研究,通過相關實驗驗證了模型的準確性,且對氣液并流流動在不同氣相Fv因子條件下填料表面的液體分布做了研究,給出其分布圖;并提出了并流流動的壓降和持液量的模擬計算方法,通過實驗驗證發(fā)現(xiàn)計算結果與實驗值吻合良好,與逆流流動相比,并流流動的壓降下降近70%,持液量下降50%左右;小氣相Fv因子的并流流動中,壓降和液相噴淋密度成反比例關系;持液量隨著氣相Fv因子的增大而減小。

    加強了對規(guī)整填料內流體流動機理的認識,從而為規(guī)整填料的進一步更新和發(fā)展提供了較好的分析方法。

    參考文獻:

    [1]張鵬, 劉春江, 成弘, 等. 規(guī)整填料塔內流體流動研究進展及展望[J]. 化學工程, 2001, 29(3): 66-70, 73

    Zhang Peng, Liu Chunjiang, Cheng Hong,etal. Progress in the studies on fluid flow in structured packing tower and its prospect[J]. Chemical Industry and Engineering, 2001, 29(3): 66-70, 73 (in Chinese)

    [2]Gu F, Liu C, Yuan X,etal. CFD simulation of liquid film flow on inclined plates[J]. Chem Eng Technol, 2004, 27(10): 1 099-1 104

    [3]陳江波. 高壓下規(guī)整填料塔的計算傳遞及傳質性能 [D].天津: 天津大學, 2006

    Chen Jiangbo. Numerical and experimental study of transport phenomena in a structured packed column at high pressure [D]. Tianjin: Tianjin University, 2006(in Chinese)

    [4]陳江波, 劉春江. 計算流體力學方法模擬規(guī)整填料塔內流體流動行為的研究進展[J]. 石油化工, 2007, 36(9): 962-965

    Chen Jiangbo, Liu Chunjiang. Progress of simulation of fluid flow behavior in structured packed column by computational fluid dynamics[J]. Petrochemical Technology, 2007, 36(9): 962-965 (in Chinese)

    [5]Petre C F, Larachi F, Iliuta I,etal. Pressure drop through structured packings: Breakdown into the contributing mechanisms by CFD modeling[J]. Chemical Engineering Science, 2003, 58: 163-177

    [6]Raynal L, Boyer C, Ballaguet J P. Liquid holdup and pressure drop determination in structured packing with CFD simulations[J]. Canadian Journal of Chemical Engineering, 2008, 82(5): 871-879

    [7]Brackbill J U, Kothe D B, Zemach C. A continuum method for modeling surface tension[J]. J Comput Phys, 1992, 100: 335-354

    [8]Bird B, Stewart W E, Lightfoot E N. Transport phenomena[M]. New York: John Wiley & Sons, Inc., 2002

    [9]Riazi M R, Faghri A. Effect of interfacial drag on gas absorption with chemical reaction in a vertical tube[J]. AIChE J, 1986, 32(4): 696-699

    [10]王樹楹. 現(xiàn)代填料塔技術指南[M]. 北京: 中國石化出版社, 1998

    [11]董誼仁, 裘俊紅, 侯章德. 現(xiàn)代填料塔技術(四) 填料塔的流體力學和傳質(上)[J]. 化工生產與技術, 1997, (1): 18-28

    [12]谷芳, 劉春江, 袁希鋼, 等. 板波紋規(guī)整填料液相分布模型[J]. 天津大學學報:自然科學與工程技術版, 2005, 38(7): 586-591

    Gu Fang, Liu Chunjiang, Yuan Xigang,etal. Liquid distribution model for sheet corrugated structured packing[J]. Journal of Tianjin University, 2005, 38(7): 586-591 (in Chinese)

    猜你喜歡
    規(guī)整液量氣液
    微重力下兩相控溫型儲液器內氣液界面仿真分析
    TAP閥壓開地層解決方案及其在蘇南的應用
    氣液分離罐液位計接管泄漏分析
    300kt/a硫酸系統(tǒng)規(guī)整填料使用情況簡介
    工程因素對頁巖氣井產水的影響分析
    CO2 驅低液量高氣液比井下氣錨模擬與優(yōu)化
    提高日用玻璃陶瓷規(guī)整度和表面光滑度的處理方法
    佛山陶瓷(2016年11期)2016-12-23 08:50:27
    電梯的建筑化藝術探索
    大觀(2016年9期)2016-11-16 10:31:30
    基于發(fā)音機制的貪婪自適應語音時長規(guī)整算法
    計算機工程(2015年8期)2015-07-03 12:20:34
    臨床問題 如何記錄出入液量更準確?
    or卡值多少钱| 国产av麻豆久久久久久久| 亚洲图色成人| 一级毛片久久久久久久久女| 精品乱码久久久久久99久播| 成人一区二区视频在线观看| 精品久久久久久久末码| 亚洲一级一片aⅴ在线观看| 国产一区二区三区av在线 | av在线播放精品| 亚洲精品日韩av片在线观看| 欧美不卡视频在线免费观看| 午夜精品国产一区二区电影 | 亚洲av熟女| 精华霜和精华液先用哪个| 欧美绝顶高潮抽搐喷水| 国产69精品久久久久777片| 欧美一区二区亚洲| 成人三级黄色视频| 伦精品一区二区三区| 一本一本综合久久| 亚洲av免费在线观看| 91在线观看av| 熟女人妻精品中文字幕| 最近视频中文字幕2019在线8| 国产午夜福利久久久久久| 国产乱人偷精品视频| 男插女下体视频免费在线播放| 菩萨蛮人人尽说江南好唐韦庄 | 97人妻精品一区二区三区麻豆| 久久久久免费精品人妻一区二区| 久久久久性生活片| 卡戴珊不雅视频在线播放| 国产精品三级大全| 久久欧美精品欧美久久欧美| 啦啦啦观看免费观看视频高清| 午夜福利视频1000在线观看| 18禁黄网站禁片免费观看直播| 亚洲精华国产精华液的使用体验 | 高清毛片免费观看视频网站| 国产av不卡久久| 成年女人看的毛片在线观看| 亚洲人成网站在线观看播放| 成人二区视频| 久久久久九九精品影院| 亚洲精品成人久久久久久| 国产亚洲精品久久久久久毛片| 亚洲人成网站高清观看| 欧美xxxx性猛交bbbb| 少妇裸体淫交视频免费看高清| 久久久午夜欧美精品| 国产黄片美女视频| 美女被艹到高潮喷水动态| 又爽又黄a免费视频| 亚洲人与动物交配视频| 日本a在线网址| 九色成人免费人妻av| 99热只有精品国产| 97碰自拍视频| 午夜福利在线观看吧| 熟妇人妻久久中文字幕3abv| 亚洲精品国产av成人精品 | 男人舔女人下体高潮全视频| 成年女人毛片免费观看观看9| 菩萨蛮人人尽说江南好唐韦庄 | 插阴视频在线观看视频| 亚洲精品日韩av片在线观看| 内地一区二区视频在线| 久久人人爽人人爽人人片va| 亚洲色图av天堂| 午夜福利在线观看吧| 日韩欧美 国产精品| 亚洲精华国产精华液的使用体验 | 亚洲精品456在线播放app| 六月丁香七月| 欧美性猛交黑人性爽| 日本 av在线| 精品熟女少妇av免费看| 一级黄色大片毛片| 91在线观看av| 久久久久久久亚洲中文字幕| 老司机福利观看| 麻豆国产av国片精品| 91在线精品国自产拍蜜月| 免费高清视频大片| 精华霜和精华液先用哪个| 波多野结衣高清无吗| 国产伦一二天堂av在线观看| 国产aⅴ精品一区二区三区波| 夜夜夜夜夜久久久久| 毛片一级片免费看久久久久| 丰满乱子伦码专区| 日韩欧美一区二区三区在线观看| ponron亚洲| 国内精品一区二区在线观看| 国产av麻豆久久久久久久| 亚洲,欧美,日韩| 精品国内亚洲2022精品成人| 成人亚洲精品av一区二区| 1024手机看黄色片| 性色avwww在线观看| 亚洲av一区综合| 亚洲性夜色夜夜综合| 成人三级黄色视频| 大型黄色视频在线免费观看| 看十八女毛片水多多多| 少妇的逼水好多| 中文字幕久久专区| 国产久久久一区二区三区| 中文字幕精品亚洲无线码一区| 午夜老司机福利剧场| 一级毛片电影观看 | 女生性感内裤真人,穿戴方法视频| 美女大奶头视频| 国产色婷婷99| 国产精品伦人一区二区| 精品日产1卡2卡| 在线免费观看不下载黄p国产| 免费人成视频x8x8入口观看| 亚洲自偷自拍三级| 丰满的人妻完整版| 成年女人看的毛片在线观看| 人妻制服诱惑在线中文字幕| 亚洲欧美日韩高清在线视频| 日本a在线网址| 欧美日本亚洲视频在线播放| 中文在线观看免费www的网站| 观看美女的网站| 欧美一区二区精品小视频在线| 国产v大片淫在线免费观看| 精品一区二区三区视频在线| 国产片特级美女逼逼视频| 亚洲四区av| 欧美在线一区亚洲| 亚洲最大成人av| 免费不卡的大黄色大毛片视频在线观看 | 我要看日韩黄色一级片| 亚洲精品久久国产高清桃花| 欧美xxxx性猛交bbbb| 国产综合懂色| 久久久欧美国产精品| 女同久久另类99精品国产91| 身体一侧抽搐| 国产精品福利在线免费观看| 亚洲综合色惰| 少妇裸体淫交视频免费看高清| 观看美女的网站| 女人十人毛片免费观看3o分钟| 久久久久九九精品影院| 人妻少妇偷人精品九色| 国产精品伦人一区二区| 欧美潮喷喷水| 亚洲av五月六月丁香网| 一个人看的www免费观看视频| 精品一区二区三区视频在线观看免费| 我的老师免费观看完整版| 亚洲国产精品sss在线观看| 搡女人真爽免费视频火全软件 | 听说在线观看完整版免费高清| 男女之事视频高清在线观看| videossex国产| 久久99热6这里只有精品| 熟女人妻精品中文字幕| 联通29元200g的流量卡| 不卡视频在线观看欧美| 亚洲精品在线观看二区| 国产不卡一卡二| 欧美精品国产亚洲| 精品免费久久久久久久清纯| 特大巨黑吊av在线直播| 非洲黑人性xxxx精品又粗又长| 观看免费一级毛片| 你懂的网址亚洲精品在线观看 | 亚洲av免费高清在线观看| 尤物成人国产欧美一区二区三区| 日日摸夜夜添夜夜添小说| 午夜福利在线观看免费完整高清在 | 熟女人妻精品中文字幕| 综合色av麻豆| 精品人妻熟女av久视频| 中文字幕av在线有码专区| 精品熟女少妇av免费看| 国产三级在线视频| 中文字幕人妻熟人妻熟丝袜美| 啦啦啦韩国在线观看视频| 男女啪啪激烈高潮av片| a级毛片a级免费在线| 一进一出抽搐动态| 亚洲国产精品国产精品| 女同久久另类99精品国产91| 久久国产乱子免费精品| 麻豆国产97在线/欧美| 亚洲无线在线观看| 亚洲欧美清纯卡通| 国产综合懂色| 久久精品91蜜桃| 日日啪夜夜撸| 一进一出好大好爽视频| 日韩成人av中文字幕在线观看 | 中文字幕精品亚洲无线码一区| 又黄又爽又刺激的免费视频.| 人妻丰满熟妇av一区二区三区| 亚洲性久久影院| av福利片在线观看| 黄片wwwwww| 亚洲内射少妇av| 日日干狠狠操夜夜爽| 嫩草影院新地址| 亚洲国产精品成人久久小说 | 亚洲成人中文字幕在线播放| 色吧在线观看| 狠狠狠狠99中文字幕| 亚洲高清免费不卡视频| 啦啦啦啦在线视频资源| 国产又黄又爽又无遮挡在线| 国产探花在线观看一区二区| 亚洲欧美日韩东京热| 国产老妇女一区| 亚洲激情五月婷婷啪啪| 男人舔女人下体高潮全视频| 在线免费观看的www视频| 一级毛片电影观看 | 欧美高清性xxxxhd video| 美女被艹到高潮喷水动态| 三级国产精品欧美在线观看| 美女免费视频网站| 舔av片在线| 可以在线观看的亚洲视频| 国产色婷婷99| 亚洲av二区三区四区| 亚洲精品日韩av片在线观看| 欧美性感艳星| 免费在线观看成人毛片| 国产精品野战在线观看| 大型黄色视频在线免费观看| 一级a爱片免费观看的视频| 少妇人妻一区二区三区视频| 久久精品国产清高在天天线| 男女边吃奶边做爰视频| 最近手机中文字幕大全| 如何舔出高潮| 亚洲最大成人中文| 日韩av不卡免费在线播放| 亚洲av免费在线观看| 日本免费a在线| 亚洲av免费高清在线观看| 中文在线观看免费www的网站| 真人做人爱边吃奶动态| 在线播放无遮挡| 国产精品一区二区三区四区免费观看 | 午夜影院日韩av| 国产精品日韩av在线免费观看| 国国产精品蜜臀av免费| 久久久久国内视频| videossex国产| 成人鲁丝片一二三区免费| 欧美国产日韩亚洲一区| 国产精品电影一区二区三区| 观看免费一级毛片| 精品少妇黑人巨大在线播放 | 欧美高清成人免费视频www| 我的老师免费观看完整版| 99久久中文字幕三级久久日本| АⅤ资源中文在线天堂| 91久久精品国产一区二区成人| 久久九九热精品免费| 老司机午夜福利在线观看视频| 亚洲国产精品国产精品| a级一级毛片免费在线观看| 亚洲经典国产精华液单| 亚洲精品一区av在线观看| av中文乱码字幕在线| 精品一区二区免费观看| 国产精品一区二区三区四区久久| 男女边吃奶边做爰视频| 一边摸一边抽搐一进一小说| 国产中年淑女户外野战色| 乱码一卡2卡4卡精品| 亚洲av二区三区四区| 日韩av不卡免费在线播放| 99热网站在线观看| 日韩欧美免费精品| 亚洲国产色片| 欧美最黄视频在线播放免费| 欧美日韩一区二区视频在线观看视频在线 | 人人妻人人澡人人爽人人夜夜 | 一级黄色大片毛片| 国产成人91sexporn| 久久天躁狠狠躁夜夜2o2o| 一进一出抽搐动态| 欧美一级a爱片免费观看看| 亚洲专区国产一区二区| 国产三级在线视频| 日本精品一区二区三区蜜桃| 99九九线精品视频在线观看视频| 97热精品久久久久久| 九九爱精品视频在线观看| 欧美日韩国产亚洲二区| 欧美日韩在线观看h| 亚洲婷婷狠狠爱综合网| 精品久久久噜噜| 亚洲电影在线观看av| 国产男人的电影天堂91| 有码 亚洲区| 男人和女人高潮做爰伦理| 国产aⅴ精品一区二区三区波| 亚洲性久久影院| 中文字幕av成人在线电影| 麻豆国产97在线/欧美| 成人二区视频| 国产大屁股一区二区在线视频| 欧美日本视频| 人人妻,人人澡人人爽秒播| 日本与韩国留学比较| 高清日韩中文字幕在线| 国产成人a∨麻豆精品| 老师上课跳d突然被开到最大视频| 免费人成在线观看视频色| 国产视频一区二区在线看| 男女做爰动态图高潮gif福利片| 欧美性感艳星| 欧美色欧美亚洲另类二区| 麻豆一二三区av精品| 国产熟女欧美一区二区| 久久99热这里只有精品18| 亚洲五月天丁香| 蜜桃久久精品国产亚洲av| 老熟妇乱子伦视频在线观看| 日本黄大片高清| 亚洲欧美精品自产自拍| 91麻豆精品激情在线观看国产| а√天堂www在线а√下载| 可以在线观看的亚洲视频| 国产在视频线在精品| 永久网站在线| 97人妻精品一区二区三区麻豆| 国产精品女同一区二区软件| 人妻久久中文字幕网| 中文字幕av在线有码专区| 国产精品无大码| 激情 狠狠 欧美| 国产黄色视频一区二区在线观看 | 亚洲国产精品合色在线| 欧美bdsm另类| 精品免费久久久久久久清纯| 一级毛片我不卡| 亚洲精品乱码久久久v下载方式| 免费av观看视频| 99国产极品粉嫩在线观看| 精品一区二区三区视频在线观看免费| 亚洲美女黄片视频| 99热只有精品国产| 国产一区二区三区在线臀色熟女| 一夜夜www| 可以在线观看的亚洲视频| 人妻丰满熟妇av一区二区三区| 黄色一级大片看看| 久久婷婷人人爽人人干人人爱| 岛国在线免费视频观看| 极品教师在线视频| 两性午夜刺激爽爽歪歪视频在线观看| 神马国产精品三级电影在线观看| 哪里可以看免费的av片| 久99久视频精品免费| 精品一区二区三区av网在线观看| h日本视频在线播放| 亚洲国产精品sss在线观看| 国产一区二区激情短视频| 国产亚洲精品久久久com| 国产乱人偷精品视频| 欧美激情久久久久久爽电影| 最新中文字幕久久久久| 一本久久中文字幕| 黄色欧美视频在线观看| 天堂网av新在线| 又黄又爽又免费观看的视频| 成人毛片a级毛片在线播放| 一夜夜www| 国产精品伦人一区二区| 黄色日韩在线| 九色成人免费人妻av| 国产精品伦人一区二区| 全区人妻精品视频| 一级毛片电影观看 | 少妇的逼水好多| 日韩欧美精品v在线| 欧美成人a在线观看| 老师上课跳d突然被开到最大视频| 欧美最新免费一区二区三区| 国产高清视频在线观看网站| 免费大片18禁| 成人二区视频| 熟女电影av网| 国产精品精品国产色婷婷| 国产v大片淫在线免费观看| 亚洲最大成人中文| 国产亚洲91精品色在线| 51国产日韩欧美| 婷婷亚洲欧美| 国产成年人精品一区二区| 亚洲成人久久爱视频| 日本熟妇午夜| 国产精品久久久久久精品电影| 日日摸夜夜添夜夜爱| 中文字幕久久专区| 欧美日韩在线观看h| videossex国产| 国产av麻豆久久久久久久| 少妇人妻一区二区三区视频| 久久韩国三级中文字幕| 免费人成在线观看视频色| 国产精品无大码| 亚洲精品影视一区二区三区av| 人妻丰满熟妇av一区二区三区| 中文亚洲av片在线观看爽| 亚洲熟妇中文字幕五十中出| 51国产日韩欧美| 日韩国内少妇激情av| 搡老岳熟女国产| 最后的刺客免费高清国语| 麻豆一二三区av精品| 成人三级黄色视频| 性插视频无遮挡在线免费观看| 亚洲图色成人| 日韩av不卡免费在线播放| 日日摸夜夜添夜夜爱| 97超碰精品成人国产| 久久久成人免费电影| 99久久精品国产国产毛片| 18+在线观看网站| 日韩人妻高清精品专区| 综合色丁香网| 亚洲人成网站高清观看| 欧美又色又爽又黄视频| 日本a在线网址| 欧美又色又爽又黄视频| 菩萨蛮人人尽说江南好唐韦庄 | 免费大片18禁| 99久国产av精品| 看十八女毛片水多多多| 男人舔女人下体高潮全视频| 午夜亚洲福利在线播放| 精品日产1卡2卡| 真人做人爱边吃奶动态| 国产精品福利在线免费观看| 亚洲国产精品成人综合色| 一卡2卡三卡四卡精品乱码亚洲| 亚洲欧美日韩东京热| 99久久精品一区二区三区| 高清日韩中文字幕在线| 久久久久性生活片| 亚洲欧美精品自产自拍| 欧美在线一区亚洲| 日韩欧美在线乱码| 搡女人真爽免费视频火全软件 | 性色avwww在线观看| 国产免费男女视频| 99久久中文字幕三级久久日本| 99久久精品国产国产毛片| 赤兔流量卡办理| 国产伦精品一区二区三区视频9| 亚洲婷婷狠狠爱综合网| 亚洲欧美日韩高清在线视频| 日日啪夜夜撸| 一边摸一边抽搐一进一小说| 中文字幕久久专区| www日本黄色视频网| av女优亚洲男人天堂| 亚洲欧美日韩东京热| 十八禁网站免费在线| 99热6这里只有精品| 亚洲av二区三区四区| 狂野欧美激情性xxxx在线观看| 亚洲不卡免费看| 日韩av不卡免费在线播放| 欧美成人a在线观看| 少妇的逼水好多| 国产蜜桃级精品一区二区三区| 一本一本综合久久| 少妇熟女aⅴ在线视频| 有码 亚洲区| 天堂影院成人在线观看| 自拍偷自拍亚洲精品老妇| 久久久国产成人免费| 亚洲最大成人中文| 成年女人永久免费观看视频| 亚洲图色成人| 欧美bdsm另类| 国产v大片淫在线免费观看| 亚洲美女黄片视频| 一个人观看的视频www高清免费观看| 啦啦啦韩国在线观看视频| 午夜爱爱视频在线播放| 黄色欧美视频在线观看| 色噜噜av男人的天堂激情| 午夜精品在线福利| 观看免费一级毛片| 久久人人精品亚洲av| 精品免费久久久久久久清纯| 亚洲第一电影网av| 天天躁夜夜躁狠狠久久av| 最近手机中文字幕大全| 免费看光身美女| 有码 亚洲区| 日韩欧美国产在线观看| 国产午夜精品久久久久久一区二区三区 | 国产黄片美女视频| 丰满乱子伦码专区| 日本三级黄在线观看| 99久久久亚洲精品蜜臀av| 一个人免费在线观看电影| 97超碰精品成人国产| 麻豆国产97在线/欧美| 大又大粗又爽又黄少妇毛片口| 国产 一区精品| 久99久视频精品免费| 亚洲经典国产精华液单| 国产精品精品国产色婷婷| 人妻丰满熟妇av一区二区三区| 久久亚洲精品不卡| 欧洲精品卡2卡3卡4卡5卡区| 久久国产乱子免费精品| 熟女人妻精品中文字幕| 国产亚洲精品久久久久久毛片| 亚洲精品国产成人久久av| av在线老鸭窝| 精品国产三级普通话版| 亚洲成a人片在线一区二区| 国产亚洲精品久久久久久毛片| 日韩强制内射视频| 亚洲在线观看片| 亚洲欧美成人综合另类久久久 | 国产毛片a区久久久久| 亚洲最大成人av| 免费人成在线观看视频色| 香蕉av资源在线| 人妻久久中文字幕网| 久久久久久九九精品二区国产| 久久韩国三级中文字幕| 国产成人freesex在线 | 成年av动漫网址| 一a级毛片在线观看| 最近2019中文字幕mv第一页| 悠悠久久av| 亚洲av免费高清在线观看| 亚洲专区国产一区二区| 日日摸夜夜添夜夜爱| 欧美日韩国产亚洲二区| 嫩草影院精品99| 久久久欧美国产精品| 女生性感内裤真人,穿戴方法视频| 欧美高清性xxxxhd video| 精品久久久久久久久亚洲| 在线看三级毛片| 亚洲精品亚洲一区二区| 国产淫片久久久久久久久| 久久精品91蜜桃| 夜夜爽天天搞| 一个人看视频在线观看www免费| 国产欧美日韩精品亚洲av| 十八禁国产超污无遮挡网站| 成人av在线播放网站| 在线天堂最新版资源| 亚洲va在线va天堂va国产| 黄色日韩在线| 99久久无色码亚洲精品果冻| 亚洲欧美日韩东京热| 日本免费a在线| 波多野结衣巨乳人妻| 亚洲精品粉嫩美女一区| 99热这里只有是精品50| 免费高清视频大片| 久久久久久久久久久丰满| 久久这里只有精品中国| 人妻久久中文字幕网| av中文乱码字幕在线| 非洲黑人性xxxx精品又粗又长| 两个人的视频大全免费| 日韩亚洲欧美综合| 欧美一区二区精品小视频在线| 简卡轻食公司| ponron亚洲| 欧美色视频一区免费| 成年版毛片免费区| 69人妻影院| 国产精品国产高清国产av| 日韩欧美免费精品| 日韩成人伦理影院| 嫩草影院精品99| 成人精品一区二区免费| 免费高清视频大片| 国产精品一区二区免费欧美| 免费av毛片视频| 校园春色视频在线观看| 精品久久久久久久久亚洲| 亚洲第一电影网av| 精品午夜福利在线看| 青春草视频在线免费观看| 国产黄a三级三级三级人| 国内精品久久久久精免费| АⅤ资源中文在线天堂| 日韩一区二区视频免费看| 国产精品日韩av在线免费观看| 欧美一级a爱片免费观看看| 九色成人免费人妻av| 国产激情偷乱视频一区二区| 欧美性猛交黑人性爽| 亚洲无线观看免费| 观看美女的网站| 国产精品一区二区三区四区免费观看 | 岛国在线免费视频观看| 一个人看的www免费观看视频| 日本免费a在线| 国产成人一区二区在线| 一区二区三区四区激情视频 |