• <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
    低射頻下管型線圈感應耦合探究
    欧美丝袜亚洲另类| 丝袜喷水一区| 少妇人妻一区二区三区视频| 国产淫片久久久久久久久| 色94色欧美一区二区| 久久久久久久精品精品| 国产成人一区二区在线| 免费大片18禁| a级片在线免费高清观看视频| 激情五月婷婷亚洲| 看免费成人av毛片| 国产午夜精品久久久久久一区二区三区| 日本欧美国产在线视频| 寂寞人妻少妇视频99o| 在现免费观看毛片| 三级国产精品欧美在线观看| 一本一本综合久久| 国产成人精品无人区| 美女内射精品一级片tv| 色网站视频免费| av不卡在线播放| 一级黄片播放器| 在线免费观看不下载黄p国产| 亚洲av欧美aⅴ国产| 免费黄色在线免费观看| 汤姆久久久久久久影院中文字幕| 午夜福利影视在线免费观看| 内地一区二区视频在线| 亚洲精品日韩av片在线观看| 亚洲人成网站在线播| 精品一区二区免费观看| 日日啪夜夜爽| 777米奇影视久久| 午夜影院在线不卡| 校园人妻丝袜中文字幕| 欧美+日韩+精品| 久久 成人 亚洲| 国产av码专区亚洲av| 一本一本综合久久| 婷婷色麻豆天堂久久| 麻豆乱淫一区二区| 欧美日韩视频精品一区| 国产91av在线免费观看| www.色视频.com| 99re6热这里在线精品视频| 精品一区二区三卡| 成人18禁高潮啪啪吃奶动态图 | 亚洲激情五月婷婷啪啪| 国产亚洲一区二区精品| av在线播放精品| 亚洲欧美日韩卡通动漫| 日韩中文字幕视频在线看片| 黄片无遮挡物在线观看| 夫妻性生交免费视频一级片| 99久久精品一区二区三区| 美女福利国产在线| 日韩av不卡免费在线播放| 纯流量卡能插随身wifi吗| 精品午夜福利在线看| 中国美白少妇内射xxxbb| 91精品国产国语对白视频| 男女免费视频国产| 亚洲精品中文字幕在线视频 | 大片电影免费在线观看免费| 大片电影免费在线观看免费| 91久久精品国产一区二区三区| 少妇人妻 视频| 最近的中文字幕免费完整| 日本欧美国产在线视频| 亚洲精品成人av观看孕妇| 一本色道久久久久久精品综合| 国产成人午夜福利电影在线观看| 精品少妇内射三级| 久久免费观看电影| 亚洲图色成人| av一本久久久久| av天堂中文字幕网| 99热6这里只有精品| 五月伊人婷婷丁香| 日韩精品免费视频一区二区三区 | 日本av手机在线免费观看| 寂寞人妻少妇视频99o| 久热久热在线精品观看| 自拍欧美九色日韩亚洲蝌蚪91 | 各种免费的搞黄视频| 欧美激情国产日韩精品一区| 热99国产精品久久久久久7| av免费观看日本| 国产免费又黄又爽又色| 最近中文字幕高清免费大全6| 高清黄色对白视频在线免费看 | 亚洲电影在线观看av| 欧美丝袜亚洲另类| 五月玫瑰六月丁香| 精品少妇黑人巨大在线播放| 日韩亚洲欧美综合| 国产真实伦视频高清在线观看| videos熟女内射| 久久人人爽av亚洲精品天堂| 日本-黄色视频高清免费观看| 久久久久久久大尺度免费视频| 高清av免费在线| 欧美日本中文国产一区发布| 国产美女午夜福利| 日本黄色日本黄色录像| 国产片特级美女逼逼视频| 美女主播在线视频| 大又大粗又爽又黄少妇毛片口| 亚洲怡红院男人天堂| 女人精品久久久久毛片| 777米奇影视久久| a级毛色黄片| 狠狠精品人妻久久久久久综合| 日韩在线高清观看一区二区三区| .国产精品久久| 国产欧美日韩综合在线一区二区 | 亚洲色图综合在线观看| 久久久久网色| 国产69精品久久久久777片| 日本欧美视频一区| 丝袜喷水一区| 啦啦啦啦在线视频资源| 欧美成人午夜免费资源| 中文字幕亚洲精品专区| 欧美日韩在线观看h| 国产在线一区二区三区精| 欧美人与善性xxx| 久久久久国产精品人妻一区二区| 丰满人妻一区二区三区视频av| 内射极品少妇av片p| 最近的中文字幕免费完整| 自拍欧美九色日韩亚洲蝌蚪91 | 成人国产麻豆网| 久久鲁丝午夜福利片| 秋霞在线观看毛片| 高清av免费在线| 午夜福利在线观看免费完整高清在| 99久久人妻综合| 啦啦啦视频在线资源免费观看| 亚洲欧美精品专区久久| 国产av一区二区精品久久| 亚洲精品日本国产第一区| 国产精品国产三级国产av玫瑰| 日韩熟女老妇一区二区性免费视频| 久久人人爽人人片av| 男女国产视频网站| 永久免费av网站大全| 国产免费又黄又爽又色| 我要看黄色一级片免费的| 美女国产视频在线观看| 久久女婷五月综合色啪小说| 18禁动态无遮挡网站| 国产高清有码在线观看视频| 午夜av观看不卡| 亚洲性久久影院| 性色av一级| 国产综合精华液| 国产黄频视频在线观看| 亚洲精品一二三| 韩国av在线不卡| 日韩强制内射视频| 久久久午夜欧美精品| 看非洲黑人一级黄片| 人人澡人人妻人| 国产精品久久久久久av不卡| 免费av中文字幕在线| 2018国产大陆天天弄谢| 国产成人精品婷婷| 国产精品免费大片| 欧美精品一区二区免费开放| 欧美激情极品国产一区二区三区 | 人人妻人人澡人人爽人人夜夜| 久久精品国产亚洲av天美| 欧美日韩国产mv在线观看视频| 中国三级夫妇交换| 国产成人精品无人区| 少妇被粗大的猛进出69影院 | 日韩,欧美,国产一区二区三区| 大话2 男鬼变身卡| 国产免费一级a男人的天堂| 成人毛片60女人毛片免费| 高清黄色对白视频在线免费看 | 国产成人精品无人区| 亚洲久久久国产精品| 老熟女久久久| 日本欧美国产在线视频| 亚洲精品国产成人久久av| 亚洲天堂av无毛| 黑丝袜美女国产一区| 国产日韩欧美在线精品| 欧美日韩av久久| 成人漫画全彩无遮挡| 3wmmmm亚洲av在线观看| 亚洲经典国产精华液单| 国产一区有黄有色的免费视频| 黄色视频在线播放观看不卡| 久久影院123| 精品国产乱码久久久久久小说| .国产精品久久| 欧美97在线视频| 韩国av在线不卡| 国产免费视频播放在线视频| 日本vs欧美在线观看视频 | av在线播放精品| 色吧在线观看| 亚洲精品乱码久久久久久按摩| 精品人妻熟女毛片av久久网站| 久久精品国产a三级三级三级| 99久久精品国产国产毛片| 男女国产视频网站| 3wmmmm亚洲av在线观看| av.在线天堂| a级片在线免费高清观看视频| 女的被弄到高潮叫床怎么办| 99九九在线精品视频 | 亚洲国产成人一精品久久久| 夜夜骑夜夜射夜夜干| 亚洲欧洲精品一区二区精品久久久 | 色视频在线一区二区三区| 国产一区二区在线观看av| 永久免费av网站大全| 国产精品99久久久久久久久| 亚洲国产精品国产精品| 国产黄色视频一区二区在线观看| 嫩草影院新地址| 桃花免费在线播放| 这个男人来自地球电影免费观看 | 亚洲真实伦在线观看| 久久女婷五月综合色啪小说| 日韩一区二区视频免费看| 极品教师在线视频| 日韩亚洲欧美综合| 久久精品久久精品一区二区三区| 97超视频在线观看视频| 最新的欧美精品一区二区| 国产男女超爽视频在线观看| 在线观看免费视频网站a站| 狠狠精品人妻久久久久久综合| 女性被躁到高潮视频| 嫩草影院新地址| 又粗又硬又长又爽又黄的视频| 五月伊人婷婷丁香| 国产免费一级a男人的天堂| 国产成人freesex在线| 亚洲三级黄色毛片| 极品教师在线视频| 精品人妻熟女av久视频| 精品少妇内射三级| av天堂久久9| 18禁裸乳无遮挡动漫免费视频| 日产精品乱码卡一卡2卡三| 欧美区成人在线视频| 在线精品无人区一区二区三| 男女边摸边吃奶| 日韩电影二区| 99精国产麻豆久久婷婷| 日韩不卡一区二区三区视频在线| 亚洲高清免费不卡视频| 丝袜在线中文字幕| 国产一区二区三区综合在线观看 | 国产精品国产三级国产专区5o| av专区在线播放| 久久99精品国语久久久| 久久久久久久精品精品| 国产精品久久久久久久久免| 一个人看视频在线观看www免费| 精品国产一区二区久久| 91精品伊人久久大香线蕉| av在线老鸭窝| 中文天堂在线官网| 久久人人爽av亚洲精品天堂| 在线观看www视频免费| 两个人免费观看高清视频 | 国产成人a∨麻豆精品| 久久精品国产鲁丝片午夜精品| 伦理电影大哥的女人| 日韩制服骚丝袜av| 色婷婷久久久亚洲欧美| 在线播放无遮挡| 只有这里有精品99| 黄色怎么调成土黄色| 亚洲国产欧美日韩在线播放 | 国产精品一区二区在线不卡| 99热全是精品| 国产成人freesex在线| 久久久久久伊人网av| 日本欧美国产在线视频| 国产av一区二区精品久久| 国产极品天堂在线| 综合色丁香网| 亚洲人成网站在线播| 国产av码专区亚洲av| 国产亚洲欧美精品永久| 国产日韩欧美视频二区| 九草在线视频观看| av又黄又爽大尺度在线免费看| 少妇 在线观看| 欧美日本中文国产一区发布| av免费观看日本| 亚洲欧美日韩卡通动漫| 亚洲人成网站在线播| 黄色日韩在线| 久久久久视频综合| 97精品久久久久久久久久精品| 亚洲三级黄色毛片| 97在线人人人人妻| 成年美女黄网站色视频大全免费 | 免费高清在线观看视频在线观看| 国产综合精华液| 中文字幕亚洲精品专区| 国产精品一区二区在线观看99| 日本av免费视频播放| 女性被躁到高潮视频| h日本视频在线播放| 精品久久国产蜜桃| 高清在线视频一区二区三区| 22中文网久久字幕| 日本av手机在线免费观看| h视频一区二区三区| 一区在线观看完整版| 中文字幕精品免费在线观看视频 | 五月天丁香电影| 国产伦在线观看视频一区| 在线观看美女被高潮喷水网站| 美女视频免费永久观看网站| 中文字幕人妻丝袜制服| 精品一区在线观看国产| 蜜桃久久精品国产亚洲av| 寂寞人妻少妇视频99o| 9色porny在线观看| 97在线人人人人妻| 2022亚洲国产成人精品| 免费观看性生交大片5| 国产色爽女视频免费观看| 久久国产精品大桥未久av | av线在线观看网站| 看免费成人av毛片| 国产一区二区在线观看日韩| 精品一区二区三区视频在线| 亚洲久久久国产精品| 在线观看三级黄色| 在线亚洲精品国产二区图片欧美 | 汤姆久久久久久久影院中文字幕| 久久97久久精品| av黄色大香蕉| 成人毛片a级毛片在线播放| 综合色丁香网| 日本欧美视频一区| 日本午夜av视频| 美女cb高潮喷水在线观看| 九草在线视频观看| 啦啦啦在线观看免费高清www| 久久99热6这里只有精品| 亚洲成人av在线免费| 在线观看国产h片| 涩涩av久久男人的天堂| 久久99精品国语久久久| 亚洲综合精品二区| 中国美白少妇内射xxxbb| 中国三级夫妇交换| 妹子高潮喷水视频| 欧美精品一区二区大全| 亚洲婷婷狠狠爱综合网| 欧美老熟妇乱子伦牲交| 人人妻人人澡人人看| 久久久久国产精品人妻一区二区| 日韩av不卡免费在线播放| 熟女人妻精品中文字幕| 美女大奶头黄色视频| 国产成人aa在线观看| 久久免费观看电影| 国产精品国产三级国产专区5o| 亚洲精品久久久久久婷婷小说| 丝袜脚勾引网站| 在线观看av片永久免费下载| 亚洲精品第二区| 男女啪啪激烈高潮av片| 久久人人爽人人片av| 汤姆久久久久久久影院中文字幕| 免费看av在线观看网站| 精品国产乱码久久久久久小说| 亚洲人成网站在线播| 免费播放大片免费观看视频在线观看| 乱人伦中国视频| 国产伦精品一区二区三区视频9| 亚洲av在线观看美女高潮| 精品国产一区二区久久| 国产亚洲av片在线观看秒播厂| 亚洲一级一片aⅴ在线观看| 99热这里只有是精品50| 男女边吃奶边做爰视频| 在线播放无遮挡| 青青草视频在线视频观看| 久久精品国产亚洲av涩爱| 久久久久久久精品精品| 色94色欧美一区二区| 成人无遮挡网站| 欧美三级亚洲精品| 亚洲欧美中文字幕日韩二区| 国产极品天堂在线| 嘟嘟电影网在线观看| av福利片在线| av在线app专区| 亚州av有码| 美女中出高潮动态图| 欧美国产精品一级二级三级 | 亚洲成人手机| 欧美xxⅹ黑人| 在线免费观看不下载黄p国产| 一区二区三区精品91| 美女视频免费永久观看网站| 亚洲国产av新网站| 国产精品麻豆人妻色哟哟久久| 大码成人一级视频| 亚洲人成网站在线观看播放| 黄色怎么调成土黄色| av有码第一页| 在线观看美女被高潮喷水网站| 国产欧美日韩精品一区二区| 最近最新中文字幕免费大全7| 51国产日韩欧美| 一级毛片黄色毛片免费观看视频| 我的女老师完整版在线观看| 嘟嘟电影网在线观看| 日日啪夜夜撸| 老司机亚洲免费影院| 美女脱内裤让男人舔精品视频| 午夜福利视频精品| 亚洲久久久国产精品| 青春草亚洲视频在线观看| 日韩精品免费视频一区二区三区 | 高清欧美精品videossex| 日韩,欧美,国产一区二区三区| 亚洲精品国产成人久久av| 欧美97在线视频| 国产成人免费观看mmmm| 丝袜在线中文字幕| 2022亚洲国产成人精品| 国产精品久久久久久精品电影小说| 成人亚洲欧美一区二区av| 一区二区三区四区激情视频| 亚洲欧美成人精品一区二区| 欧美精品国产亚洲| av黄色大香蕉| 少妇丰满av| 香蕉精品网在线| 亚洲精品,欧美精品| 欧美xxxx性猛交bbbb| 在线看a的网站| 亚洲va在线va天堂va国产| h日本视频在线播放| 国产在线男女| 亚洲美女视频黄频| 中文在线观看免费www的网站| 国产亚洲欧美精品永久| 午夜福利,免费看| 少妇猛男粗大的猛烈进出视频| 中文字幕人妻丝袜制服| 少妇熟女欧美另类| 人人妻人人澡人人看| 一级爰片在线观看| 国产日韩一区二区三区精品不卡 | 男女免费视频国产| 久久精品国产亚洲网站| 一区二区av电影网| 国产色爽女视频免费观看| 黑人猛操日本美女一级片| 久久精品国产亚洲av涩爱| 久久ye,这里只有精品| 国产伦理片在线播放av一区| 嫩草影院入口| 国产成人精品无人区| 久久99热6这里只有精品| 成人国产麻豆网| 国产免费又黄又爽又色| 亚洲精品久久午夜乱码| 久久久国产精品麻豆| 在现免费观看毛片| 91aial.com中文字幕在线观看| 亚洲图色成人| 国产精品欧美亚洲77777| 在线精品无人区一区二区三| 最后的刺客免费高清国语| 边亲边吃奶的免费视频| 国产精品伦人一区二区| 日韩大片免费观看网站| 国产午夜精品一二区理论片| 美女福利国产在线| 男的添女的下面高潮视频| 男人狂女人下面高潮的视频| 一本—道久久a久久精品蜜桃钙片| 久久这里有精品视频免费| 国产亚洲最大av| 精品酒店卫生间| 热99国产精品久久久久久7| 日产精品乱码卡一卡2卡三| 高清毛片免费看| 午夜激情久久久久久久| 国产在线视频一区二区| 男人和女人高潮做爰伦理| 国产精品福利在线免费观看| 91成人精品电影| 亚洲精品国产成人久久av| 亚洲中文av在线| 国产老妇伦熟女老妇高清| 亚洲美女黄色视频免费看| 91成人精品电影| 亚洲欧洲日产国产| 伦理电影大哥的女人| 亚洲国产色片| 日韩中字成人| 日韩精品免费视频一区二区三区 | 美女主播在线视频| 日韩视频在线欧美| 91在线精品国自产拍蜜月| av福利片在线观看| 国产中年淑女户外野战色| 成人国产麻豆网| 日韩一区二区视频免费看| 亚洲,一卡二卡三卡| 国产成人精品久久久久久| 免费看av在线观看网站| 国产一区亚洲一区在线观看| 欧美+日韩+精品| 高清欧美精品videossex| 国产黄频视频在线观看| 久久久a久久爽久久v久久| 一级二级三级毛片免费看| 插逼视频在线观看| 国产日韩欧美视频二区| 女人精品久久久久毛片| 亚洲无线观看免费| 久久国内精品自在自线图片| 日本与韩国留学比较| 久久久国产欧美日韩av| 少妇裸体淫交视频免费看高清| 欧美三级亚洲精品| 丝袜在线中文字幕| 国产国拍精品亚洲av在线观看| 久久精品国产亚洲av涩爱| 国产精品熟女久久久久浪| 哪个播放器可以免费观看大片| 黄色欧美视频在线观看| 高清不卡的av网站| 日韩人妻高清精品专区| 老熟女久久久| 在线观看国产h片| 高清午夜精品一区二区三区| 超碰97精品在线观看| 日韩三级伦理在线观看| 免费看光身美女| 亚洲国产色片| 精品国产一区二区久久| 亚洲精品一二三| 日韩欧美精品免费久久| 欧美精品一区二区免费开放| 夜夜爽夜夜爽视频| 一二三四中文在线观看免费高清| 22中文网久久字幕| 伊人久久国产一区二区| 成人18禁高潮啪啪吃奶动态图 | 蜜臀久久99精品久久宅男| 老女人水多毛片| 久久人人爽av亚洲精品天堂| 三级国产精品欧美在线观看| 在线亚洲精品国产二区图片欧美 | 国产中年淑女户外野战色| 亚洲av.av天堂| 国产精品99久久久久久久久| 国产片特级美女逼逼视频| 国内揄拍国产精品人妻在线| 午夜免费男女啪啪视频观看| 一本久久精品| 69精品国产乱码久久久| av专区在线播放| av国产久精品久网站免费入址| 久久99蜜桃精品久久| 欧美97在线视频| 国产成人a∨麻豆精品| 日产精品乱码卡一卡2卡三| 国产精品麻豆人妻色哟哟久久| av国产精品久久久久影院| 国内少妇人妻偷人精品xxx网站| 伊人久久国产一区二区| 少妇裸体淫交视频免费看高清| 午夜福利网站1000一区二区三区| 国产深夜福利视频在线观看| 9色porny在线观看| 亚洲精品aⅴ在线观看| 高清毛片免费看| 欧美性感艳星| 午夜福利网站1000一区二区三区| 老熟女久久久| 欧美日韩视频高清一区二区三区二| 永久免费av网站大全| 乱人伦中国视频| 综合色丁香网| 精品人妻偷拍中文字幕| 91成人精品电影| 久久97久久精品| 久久久午夜欧美精品| 在线观看一区二区三区激情| 久久国产亚洲av麻豆专区| av.在线天堂| 久久精品国产亚洲网站| 人人妻人人澡人人看| 国产免费又黄又爽又色| 久久影院123| 免费av不卡在线播放| 国产精品熟女久久久久浪| 最新的欧美精品一区二区| 天堂俺去俺来也www色官网| 视频区图区小说| 久久精品熟女亚洲av麻豆精品| 校园人妻丝袜中文字幕|