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

    橢圓管外液膜流動的數(shù)值模擬及傳熱實驗研究Ⅰ:數(shù)值模擬研究

    2012-04-09 09:15:16齊春華馮厚軍呂慶春
    化學工業(yè)與工程 2012年4期
    關鍵詞:管外管型圓管

    齊春華,徐 克,馮厚軍,呂慶春

    (1.國家海洋局天津海水淡化與綜合利用研究所,天津 300192;2.國家海水利用工程技術(shù)研究中心,天津 300192)

    低溫多效蒸餾海水淡化主要采用水平管降膜蒸發(fā)技術(shù),水平管降膜蒸發(fā)器傳熱傳質(zhì)性能直接影響著海水淡化成本。長期以來,作為熱交換的核心元件,傳熱管均是以圓形截面為主,隨著計算流體力學和微觀傳熱技術(shù)的發(fā)展,人們越來越認識到橢圓管對液膜流動和強化傳熱的影響,并成功應用于空調(diào)與制冷、化工等行業(yè),表現(xiàn)出較好的節(jié)能效果。國內(nèi)外許多專家學者對橢圓管換熱開展了多方面的研究[1-4],2003年,Hasan[5]對橢圓管和圓管進行了試驗研究,表明橢圓管的傳熱與流阻性能是圓管的1.93~1.96倍。在換熱設計中,殼程阻力降是很重要的,它直接影響工藝流程上下游的穩(wěn)定性、泵功耗和操作費用[6]。在低溫多效海水淡化過程中,管外液膜分布狀態(tài)將直接影響傳熱效率以及噴淋供液泵的耗能。采用橢圓管換熱元件,對改變管外薄膜流動狀態(tài)、提升降膜蒸發(fā)效率有積極作用,但橢圓管降膜過程中的管外液膜流動相關研究較少,其在低溫多效海水淡化中的應用研究也未見報道。

    本研究采用開發(fā)后的FLUENT對橢圓管外液膜流動進行了三維數(shù)值模擬研究,得出管型對液膜分布、液膜厚度的影響,為開發(fā)基于橢圓管的降膜蒸發(fā)海水淡化技術(shù)做探索性研究。

    1 機理分析

    水平管降膜蒸發(fā)是20世紀70年代發(fā)展起來的一項蒸餾淡化技術(shù),液體從分布器噴淋到管子頂部,即成薄膜狀沿管子兩側(cè)流下進行降膜蒸發(fā),直到管底匯聚成液滴滴到下一排管子上,如圖1所示。因為傳熱管兩側(cè)均有相變發(fā)生,換熱系數(shù)較高,而管外壁降膜流動狀態(tài)是蒸發(fā)側(cè)傳熱膜系數(shù)的決定因素。

    圖1 圓管噴淋降膜物理模型Fig.1 Laminar flow physical-model of spraying liquid films on horizontal tube

    圖1為水平管外噴淋降膜的物理模型,在忽略二次蒸汽對液膜的剪切力作用時,可以認為液體流動的推動力是重力沿管壁的切向分力Gg,Gg隨著周向角θ增大有明顯的變化。在0<θ<π/2區(qū)域內(nèi),Gg隨著θ增大逐漸增大,液膜運動速度逐漸加快,造成上游液體供應不上,液膜減薄;在θ=π/2處,Gg等于重力G,液體受到推動力最大,流速達到最大值,液膜厚度出現(xiàn)最小值;在π/2<θ<π區(qū)域內(nèi),Gg隨著θ增大逐漸減小,液膜的運動速度也不斷減小,由于下游液體流速比上游慢,故在下游出現(xiàn)液體堆積,液膜厚度不斷增加,直至θ=π處,Gg=0,液體脫離管壁向下滴??梢姡讦?π/2處的附近區(qū)域,液膜流速大,膜厚小,擾動加劇,故熱阻較小,傳熱系數(shù)較高。若使這個較高傳熱系數(shù)的區(qū)域所占比例變大,傳熱效果會得到強化,所以截面形狀為橢圓型的傳熱管就是根據(jù)上述基礎理論開發(fā)的。在形狀因子E(長短軸之比)大于1的橢圓管周長中的較大區(qū)域(簡稱為B區(qū))上,作用在液膜上的推動力較大,而且這個區(qū)域隨形狀因子E的增大而增大。圖2是Moalen等[7]通過數(shù)學推導得出的結(jié)果。在相同液膜雷諾數(shù)和傳熱溫差條件下,E=2.0的橢圓管B區(qū)占半周長的75%,而E=0.5的橢圓管和E=1.0的圓管分別只占14%和54%;當E>2.0時,B區(qū)范圍很大,噴淋蒸發(fā)傳熱膜系數(shù)he很大,而E<2.0(E=0.5和1.0)的傳熱管噴淋蒸發(fā)傳熱性能要差得多。

    圖2 傳熱膜系數(shù)的周向分布[7]Fig.2 Variation of the heat transfer coefficient in circumferential direction[7]

    2 數(shù)值模擬方法

    采用計算流體力學方法對橢圓管型建立液膜流動模型進行數(shù)值模擬,確定管外液膜流動狀態(tài)較好地管型結(jié)構(gòu)尺寸。

    2.1 物理模型建立

    水平管降膜蒸發(fā)器橢圓管單管外流動換熱模型如圖3所示,傳熱管上方有小孔布液器。傳熱管頂部正對布液器開孔中心,過布液器開孔中心的縱剖面作為研究截面。利用scheme語言技術(shù)進行參數(shù)化設計;在Gambit中針對橢圓管型分別進行參數(shù)化建模、網(wǎng)格劃分和邊界條件設置;在Fluent中進行了參數(shù)化初始條件設置、邊界條件設置、預處理、求解和后處理。

    圖3 橢圓管單管外流動換熱模型Fig.3 Flow and heat transfer model outside the single elliptical tube

    2.2 基本假設及流體熱物性

    假設所選取的計算區(qū)域內(nèi)流體的流動為湍流,計算流體為自來水,處于蒸發(fā)器內(nèi)壓力為0.027 MPa的工作環(huán)境,假設蒸發(fā)器內(nèi)部除了計算流體外充滿了空氣,流體是不存在換熱情況下的純流動,其物性參數(shù)為常數(shù),如表1所示。

    表1 流體物性參數(shù)Table 1 Thermophysical parameters of liquid

    2.3 控制方程和VOF方法

    傳熱管外汽液兩相流動遵循的基本方程組,包括動量、質(zhì)量和能量守恒方程和雷諾方程,以及為封閉雷諾方程提出的各種湍流模型等。通過對一定區(qū)域內(nèi)基本方程組進行分析,就可以得到在特定情況下流體流動的規(guī)律。其通用形式:

    (1)

    Φ=1,ΓΦ=0,SΦ=0

    (2)

    得出連續(xù)性方程。

    (3)

    得出動量方程。

    (4)

    得出雷諾方程。

    采用模擬兩相流通用的VOF方法來追蹤氣液相界面。在VOF表達式中,動量方程與前面的表達式相同,而且求解方法也沒有變化。但連續(xù)方程需要稍作修正。對于某一相有如下表達式:

    ?α2/?t+v·α2=0

    (5)

    式(5)中,α2為第2相得體積分數(shù)。第1相的體積分數(shù)為:α1=1-α2。流體性質(zhì)的估算就可以通過體積分數(shù)來表達,例如,在任何網(wǎng)格單元的物質(zhì)密度可用式(6)來表達:

    ρ=α2ρ2+(1-α2)ρ1

    (6)

    式(6)中,ρ是在某一計算單元格內(nèi)的估算密度,ρ1和ρ2分別為第1相和第2相的密度。

    2.4 網(wǎng)格劃分和邊界條件

    為了提高計算精度,對計算區(qū)域網(wǎng)格的處理采用分區(qū)劃分、局部加密的方式。整個計算域長(x)×高(y)×寬(z)=18 mm × 50 mm × 40 mm,分成壁面區(qū)、臨界壁面區(qū)、入流區(qū)、軸向區(qū)、出流區(qū)和非主要區(qū),共6個區(qū),全部采用六面體網(wǎng)格,網(wǎng)格總數(shù)為1 420 000個,計算區(qū)域網(wǎng)格如圖4所示。

    圖4 網(wǎng)格劃分圖Fig.4 Meshed configuration of the bath

    建立相應的邊界條件如下:1)計算區(qū)域邊線包括2種邊界:速度入口(Velocity inlet)邊界及Wall邊界,速度入口邊界位于計算區(qū)域上邊線,寬度為1.5 mm; 2)Wall邊界為傳熱管外壁,其余部分為對稱邊界。模擬過程中計算區(qū)域初始給定為全氣相,布液孔入流液體為全液相。

    3 模擬結(jié)果與分析

    3.1 管型對液膜分布的影響

    在入口流速分別為0.2、0.5和0.8 m/s[液體負荷分別為0.037、0.093和0.149 kg/(m·s)]的工況下,對相同截面周長的圓管型EP01和橢圓管型EP03進行管外液體流動的數(shù)值模擬計算。

    表2 EP01和EP03管型參數(shù)Table 2 The structural parameters of the EP01 and EP03 tubes

    圖5為圓管型不同入口流速下的液膜分布情況,其中圖5a)為EP01圓管型不同入口流速下的液膜分布情況,圖5b)為EP03橢圓管型不同入口流速下的液膜分布情況。

    圖5 不同入口流速下的液膜分布情況Fig.5 Falling film distribution chart with the velocity

    從圖5a)可見,液滴在每個管型頂部完成滴落、鋪展、下落和脫離再到恢復初始狀態(tài)的全過程,其脫落時間隨著入口流速的增大而減少。從圖5b)可見,橢圓管上的液膜在流速為0.2 m/s時,液體呈不連續(xù)的滴狀模式流動,傳熱管表面存在大量干壁區(qū)域;流速為0.5 m/s時,管外液體流型呈柱狀模式,但液體不能完全潤濕傳熱管外壁,仍存在干壁區(qū)域;流速為0.8 m/s時,管外液體流型為滴-柱狀模式,一層薄液膜包裹在傳熱管外壁,而當流速增加至1.0 m/s時,管外液體流型呈現(xiàn)出柱狀模式,管外包裹了一層較厚的液膜,并且2個入流口中間區(qū)域的液膜厚度明顯增加,形成明顯的隆起。

    由圖5a)和圖5b)對比可見,液膜在圓管底端脫落的分布較為零亂,液膜在管底部聚集,脫落較慢,使管底部液膜厚度大大增加,從圖5b)入口流速0.8 m/s時可見,液膜在橢圓管底端脫落呈間隔柱狀,均勻程度比圓管有較大程度的改善。橢圓管外液膜局部流速增大,薄膜效果增強,且很快從管底部脫落,基本沒有液膜聚集現(xiàn)象。

    由圖5可見,橢圓管外液膜的形成、鋪展、下落和脫離所用時間都比圓管短,入口流速0.8 m/s時,橢圓管用時335 ms,而圓管用時375 ms,表明橢圓管外液體流速快,更新快,擾動增強,利于換熱。

    3.2 管型對液膜厚度的影響

    液體在管表面的分布決定著管外液膜厚度,而傳熱管管型對管外液體的分布有重要影響。在液體負荷為0.149 kg/(m·s),即入口流速0.8 m/s時,根據(jù)模擬結(jié)果對管外壁液膜厚度沿管周和管軸方向的變化情況進行了測量,測量結(jié)果見圖6,其中圖6a)為管外液膜厚度分布,圖6b)為管長方向液膜厚度分布(θ=90°),θ為測點和管重心的連線與過管重心的垂線之間的夾角。

    由圖6a)可見,液膜在傳熱管表面的開始分布及流動階段,在大部分區(qū)域,圓管的膜厚比橢圓管厚10%~20%,而只有在管周液體流速最大處(θ=90°),圓管外液膜厚度才與橢圓管外液膜厚度基本相當;液膜快脫落時,在管底端的液膜厚度,圓管達3.6~3.8 mm,而橢圓管也達2.7~2.9 mm,分別是θ=90 °時的11和8倍;橢圓管僅在最低點厚度較厚,在其他區(qū)域則沒有任何液膜聚集現(xiàn)象,橢圓管液膜厚度在整個管周區(qū)域均低于圓管。

    由圖6b)為在過管軸的平面與管周相交處(θ=90 °)液膜厚度沿管長方向的變化,可見橢圓管外的液膜較圓管外薄,在L=9 mm處膜厚最厚,這是因為2個入流口(0 mm處和18 mm處)中間區(qū)域會因液膜的相向鋪展而造成液體堆積,但由圖6可見橢圓管外的液體局部堆積現(xiàn)象明顯減輕。

    圖6 液膜厚度分布Fig.6 Film thickness vs.variation

    可見,橢圓管相對于圓管,其管外壁液膜脫落速度加快,分布均勻,易形成穩(wěn)定柱狀流,管表面的液膜薄,更利于傳熱。

    3.2.1同一流速下不同橢圓管管外液膜厚度對比

    在入口流速為0.8 m/s的工況下,對物理模型就截面周長相等而截面形狀因子不同的幾種橢圓管進行管外液體流動的數(shù)值模擬計算,結(jié)構(gòu)尺寸如下表:

    表3 5種不同橢圓管管型參數(shù)Table 3 The structural parameters of five kind of elliptical tubes

    由圖7可見,截面形狀因子E越大,橢圓管外膜厚越小,但隨著E的不斷增大,液膜的不均勻性也越來越嚴重(見圖8)。如當E=2.1時,θ=170 °和θ=20 °處,橢圓管外液膜厚度分別為0.55 mm、0.16 mm,相差2倍多。E=1.5要比E=1.2的液膜薄且更加均勻。

    圖7 不同橢圓管外的液膜厚度變化Fig.7 The falling film thickness around the different elliptical tube

    圖8 EP05橢圓管外液膜厚度分布圖Fig.8 Vaper-liquid phase distribution chart of the EP05 elliptical tube

    通過模擬研究,E=1.5的橢圓管具有液膜薄,均勻性好的特點,管外液體流動形態(tài)更易形成柱狀流,有利于傳熱。

    3.2.2同一橢圓管不同入口流速橢圓管外液膜厚度對比

    對E=1.5的橢圓管,在入口流速為0.2~1.4 m/s(θ=90 °)時進行管外液體流動的數(shù)值模擬計算。

    圖9 入口流速的變化對橢圓管外側(cè)液膜水平處厚度影響Fig.9 Variation of falling film thickness around elliptical tube with the velocity

    由圖9可見,隨著入口流速的增大,橢圓管外的液體負荷也增加,管外液膜的厚度相應增大。入口流速為0.2 m/s時,模擬結(jié)果顯示管外存在較大面積的干壁,降低蒸發(fā)效率;當入口流速大于0.8 m/s時,液膜厚度增加1.5倍左右,影響傳熱,且提高供應液體泵的功率消耗,不利于節(jié)能。

    模擬結(jié)果表明,橢圓管形狀因子E=1.5,入口流速V=0.8 m/s時,管外液膜狀態(tài)最好,液膜薄且均勻,利于換熱。

    4 結(jié)論

    1)首先建立橢圓管外降膜流動的物理模型,通過模擬研究管外液體流動狀態(tài),發(fā)現(xiàn)橢圓管型相對于圓管型,液體流動耗費時間縮短40 ms,說明橢圓管型上的表面液體更新速度比圓管加快,利于換熱。

    2)在傳熱管表面大部分區(qū)域的液膜,圓管的膜厚比橢圓管厚10%~20%;在管長方向上,橢圓管外的液膜厚度較圓管外薄,且局部聚集現(xiàn)象明顯減輕;截面形狀因子E越大,橢圓管外膜厚越小,但液膜的不均勻性也越嚴重。

    3)模擬結(jié)果表明,E=1.5的橢圓管具有液膜薄,均勻性好的特點,管表面更易形成柱狀流,當入口流速V=0.8 m/s時,此橢圓管外液膜狀態(tài)最好。

    參考文獻:

    [1]劉希女,楊永安,鄒同華,等.叉排橢圓管強迫對流換熱的數(shù)值模擬[J].制冷與空調(diào),2004(1):16-19

    [2]LI Q,MA L,QI S,etal. Experimental investigation on heat transfer and f low resistance of an 1.6/1.0 elliptic cylinder in cross-flow[J].Journal of Qingdao University of Science and Technology,2004,25(1):61-54

    [3]LI Q,QI S,ZHANG X,etal.Heat transfer characteristics and flow behaviors of elliptic cylinders in cross-flow[J].Journal of Qingdao University of Science and Technology,2004,25(5):434-440

    [4]SABOYA S,SABOYA F E M.Experiments on elliptic sections in one-and tow-row arrangements of plate fin and tube heat exchangers [J].Experimental Thermal and Fluid Science,2001,24:67-75

    [5]HASAN A,SIREN K.Performance investigation of plain and finned tube evaporatively cooled heat exchangers [J].Applied Thermal Engineering,2003,23 (3):325-340

    [6]HORVAT A,LESKOVAR M,MAVKO B.Comparison of heat transfer conditions in tube bundle cross-flow for different tube shapes[J].International Journal of Heat and Mass Transfer,2006,49(6):1 027-1 038

    [7]MOALEM D,SIDEMAN S.Theoretical analysis of a horizontal condenser-evaporator tube[J].International Journal of Heat and Mass Transfer,1976,19(3):259-270

    [8]齊春華,李炎,俞永江,等.水平管降膜蒸發(fā)海水淡化的傳熱試驗研究[J].化學工業(yè)與工程,2011,28(5):1-5

    猜你喜歡
    管外管型圓管
    溴化鋰制冷蒸發(fā)器中鈦橢圓管外降膜流動及傳熱特性
    一種方便連接的涂塑鋼管
    鋼管(2021年2期)2021-11-30 02:11:01
    頸動脈管外口及周邊骨性結(jié)構(gòu)的解剖特點*
    解剖學雜志(2021年4期)2021-09-17 05:51:12
    絕緣管型母線改半絕緣管型母線的工程應用
    電線電纜(2018年3期)2018-06-29 07:41:04
    一種圓管內(nèi)孔自動打磨機的設計
    尿沉渣中蠟樣管型與不同類型腎小球疾病的相關性研究
    柔性圓管在渦激振動下的模態(tài)響應分析
    圓管帶式輸送機最佳懸垂度研究
    管型結(jié)構(gòu)物維修摩擦疊焊設備研制
    焊接(2015年11期)2015-07-18 11:12:46
    低射頻下管型線圈感應耦合探究
    我的亚洲天堂| 成人黄色视频免费在线看| 伊人亚洲综合成人网| 大陆偷拍与自拍| 久久精品人人爽人人爽视色| 日韩,欧美,国产一区二区三区| 亚洲精品一二三| 在线观看免费视频网站a站| 精品亚洲乱码少妇综合久久| 999久久久精品免费观看国产| 一级黄色大片毛片| 我要看黄色一级片免费的| xxxhd国产人妻xxx| 制服人妻中文乱码| 婷婷丁香在线五月| 人人妻人人澡人人看| 波多野结衣一区麻豆| 色精品久久人妻99蜜桃| 国产免费福利视频在线观看| 叶爱在线成人免费视频播放| 美女脱内裤让男人舔精品视频| 中文字幕人妻丝袜一区二区| 久久狼人影院| 亚洲国产看品久久| 久久狼人影院| 国产精品香港三级国产av潘金莲| 妹子高潮喷水视频| 99热全是精品| www日本在线高清视频| 欧美性长视频在线观看| 午夜福利影视在线免费观看| 老司机深夜福利视频在线观看 | 在线精品无人区一区二区三| 美女高潮喷水抽搐中文字幕| 欧美黄色片欧美黄色片| 精品人妻熟女毛片av久久网站| 日韩三级视频一区二区三区| 国产精品九九99| 黑丝袜美女国产一区| 国产av精品麻豆| 久久中文字幕一级| 精品免费久久久久久久清纯 | 丝瓜视频免费看黄片| 女人高潮潮喷娇喘18禁视频| 精品亚洲成国产av| 啦啦啦免费观看视频1| 欧美黑人精品巨大| 热99久久久久精品小说推荐| 黄片小视频在线播放| 日本黄色日本黄色录像| av有码第一页| 大陆偷拍与自拍| 永久免费av网站大全| 视频在线观看一区二区三区| 亚洲av欧美aⅴ国产| 18在线观看网站| 中文欧美无线码| 女警被强在线播放| 国产精品秋霞免费鲁丝片| 成人国产av品久久久| 又紧又爽又黄一区二区| 秋霞在线观看毛片| 99热全是精品| 纯流量卡能插随身wifi吗| 亚洲精品一区蜜桃| 涩涩av久久男人的天堂| tube8黄色片| avwww免费| 国产欧美日韩一区二区三 | 青草久久国产| 在线观看免费午夜福利视频| 国产麻豆69| 中文字幕色久视频| 亚洲精品国产一区二区精华液| 欧美另类一区| 天堂8中文在线网| 久久国产亚洲av麻豆专区| 一区二区三区乱码不卡18| 热99国产精品久久久久久7| 老熟妇乱子伦视频在线观看 | 黄片小视频在线播放| 人人妻人人爽人人添夜夜欢视频| 欧美激情极品国产一区二区三区| 亚洲av电影在线进入| 黄频高清免费视频| 国产成人精品无人区| 超色免费av| 国产欧美日韩综合在线一区二区| 菩萨蛮人人尽说江南好唐韦庄| 爱豆传媒免费全集在线观看| 大片电影免费在线观看免费| 国产高清国产精品国产三级| av在线播放精品| 精品国产国语对白av| 久久久久网色| 国产视频一区二区在线看| 老司机靠b影院| 国产日韩欧美亚洲二区| 男女无遮挡免费网站观看| 国产不卡av网站在线观看| 国产在线免费精品| 国内毛片毛片毛片毛片毛片| 天堂中文最新版在线下载| 久久久久久久久久久久大奶| 亚洲五月婷婷丁香| bbb黄色大片| 亚洲第一欧美日韩一区二区三区 | 两性夫妻黄色片| 国产成人一区二区三区免费视频网站| 爱豆传媒免费全集在线观看| av一本久久久久| 亚洲成人手机| 下体分泌物呈黄色| 另类精品久久| 国产野战对白在线观看| 久久久精品免费免费高清| 国产欧美日韩精品亚洲av| 免费高清在线观看日韩| 在线看a的网站| 久久九九热精品免费| 极品人妻少妇av视频| 亚洲avbb在线观看| 天堂8中文在线网| 天天躁日日躁夜夜躁夜夜| 美女国产高潮福利片在线看| 一区二区三区激情视频| 黑人巨大精品欧美一区二区蜜桃| 精品亚洲成国产av| 日韩欧美一区视频在线观看| 亚洲第一青青草原| 久久女婷五月综合色啪小说| 性色av乱码一区二区三区2| 1024视频免费在线观看| 一区在线观看完整版| 老司机午夜福利在线观看视频 | 99精品久久久久人妻精品| 91麻豆av在线| 欧美久久黑人一区二区| 精品卡一卡二卡四卡免费| 欧美日韩av久久| 91字幕亚洲| 精品人妻在线不人妻| 国产一区二区三区综合在线观看| 国内毛片毛片毛片毛片毛片| 国产色视频综合| 少妇的丰满在线观看| 看免费av毛片| 水蜜桃什么品种好| 丝瓜视频免费看黄片| 国产av一区二区精品久久| 十八禁网站网址无遮挡| 永久免费av网站大全| 在线观看免费午夜福利视频| 国产男女内射视频| 日日摸夜夜添夜夜添小说| 一级黄色大片毛片| 国产日韩欧美亚洲二区| 国产精品二区激情视频| 狠狠狠狠99中文字幕| 免费在线观看黄色视频的| 丰满人妻熟妇乱又伦精品不卡| 99久久99久久久精品蜜桃| 一区二区三区激情视频| 国产真人三级小视频在线观看| www.av在线官网国产| 午夜免费鲁丝| 欧美日韩国产mv在线观看视频| 日韩一区二区三区影片| 国产成人系列免费观看| 19禁男女啪啪无遮挡网站| 亚洲一区二区三区欧美精品| 免费不卡黄色视频| 成人影院久久| 99国产精品一区二区三区| 三上悠亚av全集在线观看| 国产欧美日韩一区二区三 | 少妇的丰满在线观看| 国产av国产精品国产| 国产精品1区2区在线观看. | 最黄视频免费看| 伦理电影免费视频| 亚洲精品国产av成人精品| 欧美日韩黄片免| 国产淫语在线视频| 久久国产精品大桥未久av| 久久性视频一级片| 日本猛色少妇xxxxx猛交久久| 大片免费播放器 马上看| 亚洲成人免费电影在线观看| 午夜影院在线不卡| 国产精品一区二区在线观看99| 1024视频免费在线观看| 久久久久视频综合| 亚洲成人手机| 12—13女人毛片做爰片一| 欧美变态另类bdsm刘玥| 叶爱在线成人免费视频播放| 亚洲人成电影免费在线| 成人影院久久| 成年人午夜在线观看视频| 亚洲精品一区蜜桃| 国产欧美亚洲国产| 18在线观看网站| 中文精品一卡2卡3卡4更新| 亚洲av片天天在线观看| 老司机靠b影院| 久久亚洲国产成人精品v| 国产精品久久久久久人妻精品电影 | 欧美97在线视频| 国产精品av久久久久免费| tocl精华| 欧美日韩视频精品一区| 欧美中文综合在线视频| 丁香六月欧美| 日韩 亚洲 欧美在线| 九色亚洲精品在线播放| 亚洲精品一区蜜桃| 极品少妇高潮喷水抽搐| 精品人妻熟女毛片av久久网站| 18在线观看网站| 国产成人啪精品午夜网站| 亚洲自偷自拍图片 自拍| 国产在线一区二区三区精| 黄色视频不卡| 国产av又大| 久久影院123| 亚洲精品一卡2卡三卡4卡5卡 | 纵有疾风起免费观看全集完整版| 99久久国产精品久久久| 国产成人av激情在线播放| 一边摸一边做爽爽视频免费| 天天影视国产精品| 免费高清在线观看日韩| 久久久国产成人免费| 久久热在线av| 黄色视频在线播放观看不卡| 韩国精品一区二区三区| kizo精华| 高清黄色对白视频在线免费看| 成人国语在线视频| 亚洲欧美一区二区三区久久| 中文字幕制服av| 人妻一区二区av| 亚洲精品久久久久久婷婷小说| 欧美中文综合在线视频| 岛国在线观看网站| 波多野结衣av一区二区av| 国产精品自产拍在线观看55亚洲 | 美女午夜性视频免费| 欧美xxⅹ黑人| 亚洲成国产人片在线观看| 99九九在线精品视频| 韩国精品一区二区三区| 国产一区二区激情短视频 | 少妇粗大呻吟视频| 亚洲久久久国产精品| 中文欧美无线码| 老汉色av国产亚洲站长工具| 51午夜福利影视在线观看| 久久久久久亚洲精品国产蜜桃av| 国产成人精品久久二区二区免费| 999精品在线视频| 亚洲精品一二三| 超色免费av| 高潮久久久久久久久久久不卡| 亚洲精品国产区一区二| 国产一卡二卡三卡精品| 这个男人来自地球电影免费观看| 91大片在线观看| 国产精品1区2区在线观看. | 91国产中文字幕| 亚洲av电影在线观看一区二区三区| 精品少妇一区二区三区视频日本电影| 午夜福利一区二区在线看| 电影成人av| 欧美少妇被猛烈插入视频| 欧美日韩亚洲综合一区二区三区_| bbb黄色大片| 亚洲成人手机| 99国产极品粉嫩在线观看| 老司机影院毛片| 黄色片一级片一级黄色片| 久久青草综合色| 99热全是精品| 亚洲视频免费观看视频| 国产成人精品久久二区二区91| 国产男人的电影天堂91| 无限看片的www在线观看| 国产亚洲精品第一综合不卡| 成人亚洲精品一区在线观看| 美女主播在线视频| 91老司机精品| 一区二区三区激情视频| 亚洲,欧美精品.| 一区二区三区精品91| 欧美久久黑人一区二区| 亚洲精品在线美女| 精品少妇黑人巨大在线播放| 免费人妻精品一区二区三区视频| 亚洲av电影在线观看一区二区三区| 国产日韩欧美视频二区| 亚洲,欧美精品.| 制服诱惑二区| 免费在线观看日本一区| 亚洲国产欧美日韩在线播放| 免费在线观看完整版高清| 久久精品亚洲av国产电影网| 精品第一国产精品| 国产无遮挡羞羞视频在线观看| 人妻一区二区av| 如日韩欧美国产精品一区二区三区| 淫妇啪啪啪对白视频 | 亚洲av美国av| 久久性视频一级片| 黄片播放在线免费| 热re99久久精品国产66热6| 精品一品国产午夜福利视频| 中文字幕人妻丝袜制服| 操出白浆在线播放| 99久久99久久久精品蜜桃| av在线app专区| 国内毛片毛片毛片毛片毛片| 少妇的丰满在线观看| 亚洲一区中文字幕在线| 欧美97在线视频| 中文字幕人妻熟女乱码| 亚洲国产中文字幕在线视频| 精品福利观看| 一区福利在线观看| 天天操日日干夜夜撸| 国产人伦9x9x在线观看| 精品国产国语对白av| 999久久久国产精品视频| 18禁国产床啪视频网站| 中文字幕另类日韩欧美亚洲嫩草| 人妻久久中文字幕网| 国产精品一区二区精品视频观看| a级片在线免费高清观看视频| 亚洲色图 男人天堂 中文字幕| 日日摸夜夜添夜夜添小说| 涩涩av久久男人的天堂| 99国产精品一区二区蜜桃av | 国产在线免费精品| 免费不卡黄色视频| 曰老女人黄片| 精品乱码久久久久久99久播| 国产国语露脸激情在线看| 精品久久久精品久久久| avwww免费| 如日韩欧美国产精品一区二区三区| 国产成人免费无遮挡视频| 国产日韩欧美亚洲二区| 男女免费视频国产| 男女高潮啪啪啪动态图| 国产视频一区二区在线看| 中文字幕人妻丝袜一区二区| av福利片在线| 国产亚洲精品一区二区www | 一级毛片电影观看| 韩国高清视频一区二区三区| 性高湖久久久久久久久免费观看| 国产97色在线日韩免费| 一二三四在线观看免费中文在| 捣出白浆h1v1| 精品国产一区二区久久| 在线 av 中文字幕| 亚洲成人国产一区在线观看| 热99re8久久精品国产| 国产精品熟女久久久久浪| 亚洲,欧美精品.| 久久久精品国产亚洲av高清涩受| 五月天丁香电影| 中文字幕人妻丝袜一区二区| 操美女的视频在线观看| 少妇 在线观看| 亚洲中文av在线| 亚洲精品乱久久久久久| 侵犯人妻中文字幕一二三四区| 国产三级黄色录像| 少妇的丰满在线观看| 日韩一区二区三区影片| 91成年电影在线观看| 午夜成年电影在线免费观看| 欧美日韩亚洲高清精品| 久久99一区二区三区| 久久影院123| 精品亚洲乱码少妇综合久久| cao死你这个sao货| 天堂8中文在线网| 在线观看人妻少妇| 久久久久国内视频| 亚洲av片天天在线观看| 久久精品国产亚洲av香蕉五月 | 一区二区日韩欧美中文字幕| 国产精品成人在线| 9191精品国产免费久久| avwww免费| 国产成人av教育| 下体分泌物呈黄色| 国产在线免费精品| 美女扒开内裤让男人捅视频| 飞空精品影院首页| 大型av网站在线播放| 精品国产乱子伦一区二区三区 | 亚洲第一欧美日韩一区二区三区 | www日本在线高清视频| 欧美午夜高清在线| a在线观看视频网站| 91精品国产国语对白视频| 国产日韩欧美在线精品| 亚洲五月婷婷丁香| 久久精品国产综合久久久| 80岁老熟妇乱子伦牲交| 国产极品粉嫩免费观看在线| 亚洲一区中文字幕在线| 老汉色av国产亚洲站长工具| 久久久久国产一级毛片高清牌| 另类精品久久| kizo精华| 国产精品偷伦视频观看了| 亚洲欧美一区二区三区久久| 1024香蕉在线观看| 菩萨蛮人人尽说江南好唐韦庄| 两性午夜刺激爽爽歪歪视频在线观看 | 超碰成人久久| 日本精品一区二区三区蜜桃| 精品欧美一区二区三区在线| 少妇粗大呻吟视频| a级毛片在线看网站| 黄色怎么调成土黄色| tube8黄色片| 国产欧美日韩一区二区精品| 色老头精品视频在线观看| 亚洲中文av在线| 国产片内射在线| 少妇精品久久久久久久| 99精品久久久久人妻精品| 免费在线观看视频国产中文字幕亚洲 | 免费看十八禁软件| 国产在视频线精品| 十八禁网站免费在线| 欧美少妇被猛烈插入视频| 波多野结衣av一区二区av| 国产成人av教育| 激情视频va一区二区三区| 老鸭窝网址在线观看| 99九九在线精品视频| 18禁国产床啪视频网站| 大片电影免费在线观看免费| 99国产极品粉嫩在线观看| 日本撒尿小便嘘嘘汇集6| 精品少妇一区二区三区视频日本电影| 丝袜美腿诱惑在线| 免费在线观看日本一区| 日本欧美视频一区| 别揉我奶头~嗯~啊~动态视频 | 人成视频在线观看免费观看| 99久久99久久久精品蜜桃| 一二三四社区在线视频社区8| 国产精品 欧美亚洲| 丝袜喷水一区| 日本一区二区免费在线视频| 超碰97精品在线观看| 亚洲激情五月婷婷啪啪| av不卡在线播放| 天天添夜夜摸| 母亲3免费完整高清在线观看| svipshipincom国产片| 国产野战对白在线观看| 在线观看免费午夜福利视频| 亚洲精品成人av观看孕妇| 亚洲九九香蕉| 真人做人爱边吃奶动态| 美女高潮喷水抽搐中文字幕| 久久久国产一区二区| 亚洲伊人色综图| 国产欧美日韩综合在线一区二区| 亚洲va日本ⅴa欧美va伊人久久 | 9色porny在线观看| 国产精品国产三级国产专区5o| 成在线人永久免费视频| 777米奇影视久久| 十八禁高潮呻吟视频| 亚洲精品国产精品久久久不卡| 亚洲精品乱久久久久久| 在线观看免费高清a一片| 亚洲成av片中文字幕在线观看| 男人操女人黄网站| 午夜成年电影在线免费观看| tube8黄色片| 桃红色精品国产亚洲av| 一边摸一边抽搐一进一出视频| 欧美变态另类bdsm刘玥| 午夜久久久在线观看| 在线观看人妻少妇| 国产黄频视频在线观看| 午夜激情av网站| 日韩熟女老妇一区二区性免费视频| 国产成人av激情在线播放| 丝袜喷水一区| 性少妇av在线| 秋霞在线观看毛片| 老鸭窝网址在线观看| 中文字幕高清在线视频| 国产欧美亚洲国产| 熟女少妇亚洲综合色aaa.| 国产一区有黄有色的免费视频| 久久精品成人免费网站| 成人18禁高潮啪啪吃奶动态图| 十八禁高潮呻吟视频| 18在线观看网站| 久久精品久久久久久噜噜老黄| 美国免费a级毛片| 欧美黄色片欧美黄色片| 两人在一起打扑克的视频| 欧美成人午夜精品| 18在线观看网站| 国产视频一区二区在线看| 色精品久久人妻99蜜桃| 啦啦啦视频在线资源免费观看| 国产精品二区激情视频| 日本撒尿小便嘘嘘汇集6| 一级a爱视频在线免费观看| 免费女性裸体啪啪无遮挡网站| 国产精品一区二区在线不卡| 丝瓜视频免费看黄片| 51午夜福利影视在线观看| 又黄又粗又硬又大视频| 国产不卡av网站在线观看| 成人国产av品久久久| 男女免费视频国产| 精品熟女少妇八av免费久了| 欧美国产精品va在线观看不卡| 久久精品国产亚洲av香蕉五月 | 亚洲伊人色综图| 国产成人影院久久av| 一本色道久久久久久精品综合| 精品亚洲成a人片在线观看| 亚洲欧美日韩另类电影网站| 老司机靠b影院| 后天国语完整版免费观看| 亚洲欧美一区二区三区久久| av视频免费观看在线观看| 2018国产大陆天天弄谢| 女人被躁到高潮嗷嗷叫费观| 一级,二级,三级黄色视频| 少妇 在线观看| 精品国产一区二区三区久久久樱花| 老司机午夜十八禁免费视频| 动漫黄色视频在线观看| 天天添夜夜摸| 91大片在线观看| 欧美黑人欧美精品刺激| 亚洲国产av新网站| 亚洲欧美精品综合一区二区三区| 肉色欧美久久久久久久蜜桃| 男女国产视频网站| 嫩草影视91久久| 91av网站免费观看| 国产深夜福利视频在线观看| 亚洲国产精品一区二区三区在线| av国产精品久久久久影院| 成人亚洲精品一区在线观看| 99久久99久久久精品蜜桃| kizo精华| 精品第一国产精品| 久热爱精品视频在线9| 国产黄频视频在线观看| 妹子高潮喷水视频| 亚洲av日韩在线播放| 色老头精品视频在线观看| 色94色欧美一区二区| 久久综合国产亚洲精品| 天天躁狠狠躁夜夜躁狠狠躁| 久久这里只有精品19| 亚洲欧美激情在线| 欧美精品一区二区免费开放| a级毛片黄视频| 中亚洲国语对白在线视频| 亚洲成人免费电影在线观看| 交换朋友夫妻互换小说| 侵犯人妻中文字幕一二三四区| 亚洲av日韩在线播放| 欧美精品亚洲一区二区| 天堂中文最新版在线下载| 丝袜脚勾引网站| 五月开心婷婷网| 国产成人av激情在线播放| 国产主播在线观看一区二区| 成人影院久久| 久热这里只有精品99| 国产一区二区 视频在线| avwww免费| 中国国产av一级| 少妇被粗大的猛进出69影院| 天堂8中文在线网| 久久国产精品男人的天堂亚洲| 欧美国产精品一级二级三级| 美女午夜性视频免费| 亚洲欧洲精品一区二区精品久久久| 青春草视频在线免费观看| 亚洲欧洲精品一区二区精品久久久| 国产又爽黄色视频| 国产精品 欧美亚洲| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品秋霞免费鲁丝片| 欧美精品高潮呻吟av久久| 日本wwww免费看| 免费少妇av软件| 亚洲精品久久久久久婷婷小说| 婷婷色av中文字幕| 香蕉丝袜av| 亚洲五月色婷婷综合| 美女主播在线视频| 欧美亚洲日本最大视频资源|