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

    泵站虹吸式出水流道駝峰排氣過程氣液兩相流研究

    2022-10-27 09:42:38劉霞黃可燦趙亮常正柏蔣勁
    中國農(nóng)村水利水電 2022年10期
    關(guān)鍵詞:虹吸式氣團駝峰

    陳 奇,劉霞,黃可燦,趙亮,常正柏,蔣勁

    (1.湖北省水利水電科學(xué)研究院,湖北武漢 430072;2.湖北省國際灌排研究培訓(xùn)中心,湖北武漢 430072;3.武漢大學(xué)動力與機械學(xué)院,湖北武漢 430072)

    0 引言

    泵站是指通過泵進行抽水、增壓,再通過相關(guān)水工筑物進行輸水、配水、排水的工程設(shè)施[1]。泵站在實際工程中應(yīng)用較為廣泛,在農(nóng)業(yè)灌溉、輸水調(diào)水、城市排澇等方面起著重要作用,為國家經(jīng)濟建設(shè)和人民的幸福生活做出了重要貢獻。出水流道作為泵站的關(guān)鍵過流部件之一,其對泵站的效率、出水流態(tài)以及泵站的穩(wěn)定運行影響較大。出水流道的形式主要有虹吸式、直管式、貓背式、屈膝式等[2]。對于虹吸式出水流道而言,其頂部裝有真空破壞閥,在水泵停機時可用來破壞虹吸切斷水流,運行可靠方便,且不會影響到防洪堤,工程造價較低,適用于出口水位變化不太大的泵站。

    目前已經(jīng)有諸多學(xué)者對虹吸式流道的相關(guān)問題進行了研究。朱紅耕[3,4]采用了三維數(shù)值模擬方法分析了虹吸管內(nèi)部流道的流態(tài),發(fā)現(xiàn)在駝峰段和下降段水流速度變化較大時,容易發(fā)生脫流現(xiàn)象。譚淋露[5]采用三維數(shù)值模擬方法分析了虹吸式出水管內(nèi)部水力特性,通過對比計算虹吸管不同的上升角、駝峰斷面高寬比、下降角和出口斷面高寬比對其水力特性的影響,確定了更優(yōu)參數(shù),降低了水力損失。模擬馮建剛[6]通過模型試驗,研究了虹吸形成的3 個主要階段,分別為水力驅(qū)氣階段、水力挾氣階段、穩(wěn)定虹吸階段。其中水力挾氣階段對形成穩(wěn)定虹吸所需時間影響較大。汪尚紅、李澤等[7]對虹吸式出水流道的軸流泵過渡過程進行了研究,發(fā)現(xiàn)開機過程出水閘門的開啟時間不宜過長,以防止機組過載;關(guān)機過程出水閘門先快關(guān)后慢關(guān),且對于不同泵站開閘和關(guān)閘相關(guān)時間都有最優(yōu)值。王曉升、馮建剛等[8,9]采用RNG k-? 湍流模型以及歐拉多相流模型對虹吸式出水流道虹吸形成過程進行了數(shù)值模擬計算,發(fā)現(xiàn)其他參數(shù)結(jié)構(gòu)固定的情況下,對水力挾氣階段影響最大的是流速,流速越大所需時間越短。M.Cihan 等[10]在數(shù)值模擬時使用VOF模型,模擬出了排水過程中空氣進入虹吸管的現(xiàn)象,在其物理模型試驗中也觀察到了與數(shù)值模擬類似的空氣夾帶現(xiàn)象,空氣進入虹吸管時會降低排水效率并且引起振動。Babaeyan-Koopaei等[11]通過模型試驗,分階段對特定的虹吸溢流道的進氣口設(shè)計和流道結(jié)構(gòu)設(shè)計進行改善,在對實驗?zāi)P椭贿M行簡單修改的情況下,發(fā)現(xiàn)在虹吸式流道特定高度的上方開一定大小和數(shù)量的水平方形槽孔,可以獲得最佳的空氣調(diào)節(jié)穩(wěn)定性,并且提高了溢洪道的容量。Rahim等[12]采用RNG k-?湍流模型和VOF模型計算出了虹吸式出水流道的流量系數(shù),并與模型試驗相對比,數(shù)據(jù)表明一致性較高,吻合較好,驗證了數(shù)值模擬的可行性。

    虹吸式出水流道形成穩(wěn)定虹吸需要經(jīng)歷4 個階段:啟動階段,水力驅(qū)氣階段、水力挾氣階段、穩(wěn)定虹吸階段[8],其中水力挾氣階段最為關(guān)鍵,水流需要將駝峰內(nèi)的大氣團帶走,才能進入到穩(wěn)定階段。若將氣團整體吹走,則該階段經(jīng)歷的時間短,但需要的流速較大,對于大型的虹吸式出水流道,整體將氣團帶走所需要的流速更大;若設(shè)計流速較小,則挾氣階段的時間將增加,流速過小甚至不能進入穩(wěn)定虹吸階段。處于該階段的時間過長,會導(dǎo)致駝峰處水力損失偏大,且水流狀態(tài)不穩(wěn)定,使泵的工況不穩(wěn)定,這將引起機組產(chǎn)生振動、噪音等。若此時水泵工作在其特性曲線的馬鞍區(qū)內(nèi),則將產(chǎn)生更強烈的振動和噪聲,處于該工況時間過長甚至有可能造成事故。因此,研究流速大小以及出水流道的自身結(jié)構(gòu)參數(shù)對虹吸形成階段的影響,以及形成穩(wěn)定虹吸所需時間,對于泵站安全高效穩(wěn)定運行具有重大意義。

    1 瞬態(tài)數(shù)值模擬方法

    1.1 基本控制方程

    流體流動遵循的三大基本定理為:能量守恒、質(zhì)量守恒、動量守恒,這三大定理結(jié)合數(shù)學(xué)方法即可得到流體運動的基本控制方程。本文所研究的虹吸式出水流道中的虹吸形成過程的可壓縮性不明顯,可忽略能量方程,僅考慮連續(xù)性方程及動量方程。

    (1)連續(xù)性方程如下:

    式中:ρ為流體密度;u為流體運動速度。

    (2)動量方程如下:

    式中:fi為單位質(zhì)量流體所受的質(zhì)量力;σij為流體的應(yīng)力張量。

    1.2 VOF多相流模型及RNG k-?湍流模型

    1.2.1 VOF多相流模型

    VOF 模型通過求解流域內(nèi)的各相的體積分數(shù)及各相的動量方程來模擬兩種及以上的流體,求解出的速度場等屬性由各相共享,所有相在每個網(wǎng)格單元里的體積分數(shù)之和均為1。VOF 模型可應(yīng)用于水中大氣泡運動,含有自由液面的流動等,故可用于模擬虹吸式出水流道內(nèi)的挾氣過程。

    在VOF 模型中通過求解各相的體積分數(shù)連續(xù)方程來追蹤各相之間的界面,其方程如下所示:

    式中:下標(biāo)i表示模型中流體的相;ρ為密度;v為速度。

    1.2.2 RNG k-?湍流模型

    線性渦黏湍流模型是目前在實際工程湍流問題中使用最為廣泛的湍流模型,如k-?、k-ω 模型等。而k-?模型經(jīng)過不斷完善后,成為目前使用最廣泛的模型,包括Standardk-?、Realizablek-?模型和RNGk-?模型。

    在之前對虹吸式出水流道的挾氣過程研究中,已經(jīng)比較了Standardk-?、Realizablek-?和RNGk-?模型模擬結(jié)果與實驗的相似程度,發(fā)現(xiàn)選用RNGk-?模型的仿真結(jié)果與實驗有較高的相似性[13],且RNGk-?模型考慮到了湍流漩渦,對低雷諾數(shù)也適用,可信度和精度較高,故選用RNGk-?模型。

    RNGk-?模型不可壓縮流體的湍動能k方程及耗散率ε方程如下所示:

    k方程:

    ε方程:

    式中:μeff是有效湍流黏度,μeff=μ+μt;常數(shù)項取值:αk=αε=0.139,=1.42,C2ε=1.68。

    1.3 幾何模型及邊界條件

    本文涉及到的模型有3個,改變的參數(shù)為駝峰截面高寬比,分別為0.4和0.35,0.3,模型的駝峰截面尺寸均相等,面積與3 m直徑的圓管相同。為了更快速地帶走氣泡,下降段需要保持較高的流速,故下降段沒有采用擴散,快到底部時采用擴散管型以減小流速。

    3 個模型總體較為相似,正視圖如圖1 所示,入口截面高度為H1,駝峰截面高度為H2,駝峰頂距離水面H0,出水管淹沒在水面下的高度為H4,出口截面高度為H3,模型長度為L0,模型寬度為W0。模型具體尺寸如表1所示。

    圖1 虹吸式出水流道正模型正視圖Fig.1 Front view of the front model of the siphon outlet channel

    表1 模型幾何尺寸數(shù)據(jù)表Tab.1 Model geometry data

    3 個模型均采用速度入口和壓力出口,入口流速按照工況不同則設(shè)置不同,由駝峰截面平均流速作為衡量標(biāo)準,出口壓力由淹沒深度及出口高度共同確定,高寬比為0.4、0.35、0.3 的模型出口壓力分別為19 300、18 400、17 420 Pa,重力方向設(shè)置為豎直向下。

    1.4 網(wǎng)格劃分

    將高寬比為0.3,下降角為30度的模型分別按照0.03、0.04、0.05 m 的網(wǎng)格單元尺寸劃分網(wǎng)格,并在相同邊界條件與求解器設(shè)置下,計算至各項殘差均低于10-4,設(shè)立如圖2 所示的ABCD四個監(jiān)測點檢測流速。計算結(jié)果如表2所示。

    圖2 速度監(jiān)測點分布圖Fig.2 Distribution of velocity monitoring points

    從表2中的數(shù)據(jù)可以看出,A、B、C三點流速幾乎沒有差別,D 點流速0.03 m 和0.04 m 單元尺寸的網(wǎng)格較為接近,0.05 m 單元尺寸的網(wǎng)格與前兩者差異相對較大,故選用0.04 m 作為網(wǎng)格劃分的主要尺寸。

    表2 速度監(jiān)測點結(jié)果表Tab.2 Velocity monitoring point results

    采用網(wǎng)格單元尺寸為0.04 m 生成的網(wǎng)格如圖3 圖4 所示。圖3 為高寬比0.4 模型的網(wǎng)格,圖4 是高寬比0.3 模型僅僅縮短入口管道后的網(wǎng)格,縮短入口的目的在于減小網(wǎng)格數(shù)量,進而縮短模擬破碎帶走的整個過程所需時間。高寬比為0.35 和0.3的模型網(wǎng)格與圖3所示網(wǎng)格類似。

    圖3 高寬比為0.4的模型網(wǎng)格圖Fig.3 Mesh diagram of the model with an aspect ratio of 0.4

    圖4 高寬比為0.3的模型縮短入口管段后的網(wǎng)格圖Fig.4 Mesh diagram of the model with an aspect ratio of 0.3 after shortening the inlet pipe

    1.5 求解設(shè)置

    VOF 模型比較適合用來計算氣團的運動,且在先前的工作中已經(jīng)確定RNG k-? 模型用來計算挾氣過程比較適合,故選用上述兩種模型來模擬氣團帶走的過程。使用瞬態(tài)求解,速度入口的速度大小根據(jù)不同工況設(shè)置,壓力出口根據(jù)出口高度以及淹沒深度共同確定,上文已經(jīng)提及。管壁為無滑移邊界,粗糙度為1 mm,粗糙系數(shù)為0.5,選用SIMPLE算法,計算時間步長為0.01 s,殘差控制全部設(shè)置為10-4,各項殘差達到10-4以下視為收斂。

    2 數(shù)值模擬結(jié)果及分析

    2.1 短時間內(nèi)氣團一次性排出計算

    該部分計算了3 個不同駝峰高寬比的虹吸式出水流道,當(dāng)駝峰頂部存在大氣團時,一次性帶走所需要的最小流速,其目的在于快速形成虹吸。氣團初始狀態(tài)如圖5 所示,氣團底部比駝峰底部略低。

    圖5 氣團初始狀態(tài)圖Fig.5 Diagram of the initial state of the air masses

    2.1.1 高寬比為0.4的模型一次性排出氣團所需最小流速

    高寬比為0.4 的虹吸式出水流道的計算結(jié)果如圖6 和圖7所示。設(shè)定駝峰截面流速為2.6 m/s,模擬至10 s 時的氣液分布及速度云圖如圖6 所示。從圖6(a)中可以看出,水流未能充滿駝峰頂部,氣泡將持續(xù)破碎并被帶向出口,該流速下不能一次性排出氣團。

    圖6 駝峰截面平均流速為2.6 m/s時的結(jié)果圖Fig.6 The result diagram when the average velocity of the hump section is 2.6 m/s

    逐漸增大駝峰截面流速,重復(fù)計算,發(fā)現(xiàn)當(dāng)流速增大至4.3 m/s,駝峰頂部可以直接充滿水。如圖7 所示,模擬至6 s時,水流可以充滿駝峰頂部,后續(xù)帶走氣團會比較迅速,但此時駝峰截面流速為4.3 m/s,流速過大,實際工程較少采用。

    圖7 駝峰截面平均流速為4.3 m/s時的結(jié)果圖Fig.7 The result diagram when the average velocity of the hump section is 4.3 m/s

    2.1.2 高寬比為0.35的模型一次性排出氣團所需最小流速

    為減小駝峰截面平均流速,將駝峰高寬比降低至0.35,計算結(jié)果如圖8 和圖9 所示。當(dāng)駝峰截面平均流速為3.5 m/s時,模擬到6 s時發(fā)現(xiàn)水流未能充滿駝峰頂部,如圖8所示。當(dāng)駝峰截面流速增大到4.1 m/s時,模擬至6 s時水流充滿整個駝峰,如圖9所示,后續(xù)氣團排出會較為迅速。

    圖8 駝峰截面平均流速為3.5 m/s時的結(jié)果圖Fig.8 The result graph when the average velocity of the hump section is 3.5 m/s

    圖9 駝峰截面平均流速為4.1 m/s時的結(jié)果圖Fig.9 The result diagram when the average velocity of the hump section is 4.1 m/s

    2.1.3 高寬比為0.3的模型一次性排出氣團所需最小流速

    為了進一步減小駝峰截面平均流速,減小駝峰高寬比至0.3,計算結(jié)果如圖10 和圖11 所示。圖10 是駝峰截面流速為3.5 m/s,模擬至6 s時的氣液分布及速度云圖,從圖10(a)中可以看出,水流未能充滿駝峰頂部,需增大流速。

    圖10 駝峰截面平均流速為3.5 m/s時的結(jié)果圖Fig.10 Result diagram when the average velocity of the hump section is 3.5 m/s

    當(dāng)截面流速增大至3.8 m/s時,模擬至6 s 時水流可以充滿駝峰頂部,如圖11(a)所示。

    綜上,高寬比為0.4、0.35、0.3 的模型一次性排出的流速分別為4.3、4.1、3.8 m/s,降低高寬比后,一次性排出氣團所需要的臨界流速有所減小,但此流速仍顯偏大,所以大口徑的虹吸式出水流道想要短時間一次性整體將氣團從駝峰頂部排出需要較大的流速。

    2.2 氣團破碎帶走計算

    從上述結(jié)果來看,駝峰截面高寬比越小,一次性排出氣團所需流速越小,故選用高寬比為0.3 的模型來進行氣泡破碎帶走全過程模擬,有利于進一步降低流速。本節(jié)模擬出了120 s和180 s 排氣時間內(nèi)將氣團全部排出所對應(yīng)的駝峰平均流速以及排氣全過程。

    2.2.1 排氣時間為120 s的氣團排出過程

    當(dāng)排氣時間為120 s時,氣團破碎排出過程及形成穩(wěn)定虹吸時的速度云圖分別如圖12 和圖13 所示。圖12 所示的(a)~(f)的6 個狀態(tài)分別為0.3、8、34、50、80、120 s 時的氣液分布狀態(tài),圖12(f)所示的狀態(tài)為最終的虹吸形成狀態(tài),氣泡完全排出,整個過程耗時120 s,此時駝峰截面平均流速為3.0 m/s。從圖13 的速度云圖可以看出,形成穩(wěn)定虹吸后,駝峰截面在豎直方向有較大的速度梯度。

    圖12 駝峰截面平均流速為3.0 m/s時的排氣過程Fig.12 The exhaust process when the average velocity of the hump section is 3.0 m/s

    圖13 駝峰截面平均流速為3.0 m/s時形成虹吸后的速度云圖Fig.13 Velocity diagram after siphon formation when the average flow velocity of the hump section is 3.0 m/s

    2.2.2 排氣時間為180 s的氣團排出過程

    當(dāng)排氣時間為180 s時,氣團破碎排氣過程及形成穩(wěn)定虹吸時的速度云圖分別如圖14 和圖15 所示。圖14 所示的(a)~(g)的6個狀態(tài)分別為0.3、8、34、50、80、120、180 s時的氣液分布狀態(tài),圖14(g)所示的狀態(tài)為最終的虹吸形成狀態(tài),氣泡幾乎完全排出,整個過程耗時180 s,此時駝峰截面平均流速為2.5 m/s。相較于120 s 的排氣時間,駝峰截面流速從3.0 m/s 降低到了2.5 m/s,降低了16.67%。從圖14(e)和圖14(f)可以看出,排出駝峰頂部的少量氣體所耗的時間較長,即駝峰頂部的氣體相較于下降段內(nèi)的氣體更難排出,結(jié)合速度云圖可以看出,由于存在垂直方向上的速度梯度,駝峰頂部流速較小,這是駝峰頂部氣體更難排出的主要原因。

    圖14 駝峰截面平均流速為2.5 m/s時的排氣過程Fig.14 The exhaust process when the average velocity of the hump section is 2.5 m/s

    圖15 駝峰截面平均流速為2.5 m/s時形成虹吸后的速度云圖Fig.15 Velocity diagram after siphon formation when the average flow velocity of the hump section is 2.5 m/s

    2.3 排氣過程不同時間點湍流動能分布

    對駝峰截面流速為2.5 m/s 的排氣過程中t=10 s 和t=80 s 時刻進行湍流動能進行分析,其分布分別如圖16、17所示。

    圖16 t=10 s時的氣液分布及湍流動能分布Fig.16 Gas-liquid distribution and turbulent kinetic energy distribution at t=10 s

    從上面的結(jié)果可以看出,t=10 s時,下降段中氣液交界面附近的湍流動能較大,氣團破碎迅速;而當(dāng)t=80 s時,駝峰段中氣液交界面附近的湍流動能大大減小,氣團破碎速度變慢。這說明湍流動能的大小會影響氣團的破碎速度。

    3 結(jié)論

    圖17 t=80 s時的氣液分布及湍流動能Fig.17 Gas-liquid distribution and turbulent kinetic energy distribution at t=80 s

    基于CFD 研究了以兩種方式帶走大口徑矩形斷面駝峰頂部大氣團,分別是一次性快速排出氣團和一定時間內(nèi)破碎排出氣團,其中前一種方式計算了0.4、0.35、0.3 三種不同高寬比的虹吸管模型,而后一種方式采用0.3 高寬比的模型模擬了兩種不同流速下氣團運動、破碎帶走、形成虹吸的整個過程,結(jié)論如下。

    (1)降低高寬比可以降低一次性排出氣團需要流速va。模型駝峰高寬比為0.4、0.35、0.3時,對應(yīng)的va分別為4.3、4.1、3.8 m/s,0.3 高寬比的模型與0.4 高寬比的模型相比,va降低了11.6%。雖然有所降低,但流速仍然偏大,所以較大口徑的虹吸式出水流道不建議采用一次性排出氣團的方式。

    (2)允許排氣時間越長,對應(yīng)的駝峰截面流速越小,駝峰頂部的氣體相比于下降段中的更難排出。當(dāng)排氣時間tb為120 s時,對應(yīng)的駝峰截面流速vb為3.0 m/s,當(dāng)tb為180 s時,對應(yīng)vb為2.5 m/s,流速降低了16.67%。由于駝峰頂部存在垂直方向上的速度梯度,頂部流速較小,且排氣過程中駝峰段的湍流動能較小,使得駝峰頂部氣體比下降段中的氣體更難排出。

    猜你喜歡
    虹吸式氣團駝峰
    虹吸式雨水排水技術(shù)在建筑給排水施工中的應(yīng)用
    泗陽二站虹吸式出水流道水力優(yōu)化設(shè)計
    輸水管線啟動填充過程含滯留氣團瞬變流數(shù)值模擬
    噴射虹吸式坐便器水包裝置試驗與結(jié)構(gòu)優(yōu)化分析
    基于后向軌跡模式的合肥市臭氧來源分析
    駝峰場編發(fā)線向到發(fā)場交授權(quán)電路探討
    駝峰第三制動位減速器夾停鉤車問題的改進
    氣團來源對沿海城市PM2.5中二次水溶性無機離子形成特征的影響研究
    HZ330CL橫折缺陷原因分析與改進
    淺析虹吸式雨水排水系統(tǒng)與重力式雨水排水系統(tǒng)
    一级毛片高清免费大全| 亚洲欧美日韩高清在线视频| 女警被强在线播放| 亚洲精品久久国产高清桃花| 亚洲性夜色夜夜综合| 久久欧美精品欧美久久欧美| 国产成+人综合+亚洲专区| 夜夜夜夜夜久久久久| 性色av乱码一区二区三区2| 亚洲精品久久成人aⅴ小说| 好看av亚洲va欧美ⅴa在| 特大巨黑吊av在线直播 | 黑人欧美特级aaaaaa片| 午夜精品久久久久久毛片777| 亚洲精华国产精华精| 亚洲精品粉嫩美女一区| 两个人免费观看高清视频| 神马国产精品三级电影在线观看 | 欧美性长视频在线观看| 人人妻人人看人人澡| 精品国产乱子伦一区二区三区| 国产亚洲欧美精品永久| 精品日产1卡2卡| 亚洲人成电影免费在线| 亚洲精品久久国产高清桃花| 久久狼人影院| 怎么达到女性高潮| 亚洲国产中文字幕在线视频| 精品久久久久久久久久久久久 | 一级a爱视频在线免费观看| 久久国产精品影院| 18禁观看日本| 嫩草影视91久久| 波多野结衣av一区二区av| 亚洲片人在线观看| 免费在线观看完整版高清| 90打野战视频偷拍视频| 99国产精品一区二区蜜桃av| 久久久久九九精品影院| 欧美日韩一级在线毛片| 精品无人区乱码1区二区| 欧美在线一区亚洲| 久久久精品欧美日韩精品| 午夜免费观看网址| 日韩欧美国产一区二区入口| 亚洲欧美精品综合一区二区三区| 国产精品电影一区二区三区| √禁漫天堂资源中文www| 人人妻人人澡人人看| 国产伦人伦偷精品视频| 波多野结衣高清作品| 宅男免费午夜| 日韩 欧美 亚洲 中文字幕| 精品国产超薄肉色丝袜足j| 国产免费男女视频| 脱女人内裤的视频| 法律面前人人平等表现在哪些方面| 91国产中文字幕| 两个人视频免费观看高清| 啦啦啦 在线观看视频| 一a级毛片在线观看| 97碰自拍视频| 变态另类成人亚洲欧美熟女| 欧美乱色亚洲激情| 中文字幕av电影在线播放| 欧美另类亚洲清纯唯美| 热99re8久久精品国产| 男人的好看免费观看在线视频 | 女生性感内裤真人,穿戴方法视频| 亚洲aⅴ乱码一区二区在线播放 | 少妇被粗大的猛进出69影院| 好男人电影高清在线观看| 日本免费一区二区三区高清不卡| 国产精品98久久久久久宅男小说| 久久久久免费精品人妻一区二区 | 国产伦一二天堂av在线观看| 我的亚洲天堂| 亚洲av电影不卡..在线观看| 中文字幕人妻熟女乱码| 欧美在线一区亚洲| 男女之事视频高清在线观看| 久久香蕉国产精品| 亚洲第一青青草原| 国产精品香港三级国产av潘金莲| 日本五十路高清| 亚洲一区二区三区色噜噜| 老司机在亚洲福利影院| 久久久久免费精品人妻一区二区 | 日韩精品中文字幕看吧| 精品高清国产在线一区| 久久久久久久久免费视频了| 婷婷丁香在线五月| 91av网站免费观看| 国产又色又爽无遮挡免费看| 妹子高潮喷水视频| 亚洲专区国产一区二区| 日韩国内少妇激情av| 免费一级毛片在线播放高清视频| 亚洲人成网站高清观看| 欧美国产日韩亚洲一区| 精品午夜福利视频在线观看一区| 最近在线观看免费完整版| 啪啪无遮挡十八禁网站| 性色av乱码一区二区三区2| 欧美性猛交╳xxx乱大交人| www.999成人在线观看| 色在线成人网| 久久久久久九九精品二区国产 | 一级毛片高清免费大全| 黑人欧美特级aaaaaa片| 午夜久久久久精精品| 国产成人啪精品午夜网站| 国产成人欧美在线观看| 三级毛片av免费| 成年人黄色毛片网站| 亚洲精华国产精华精| 一夜夜www| 久久天躁狠狠躁夜夜2o2o| 免费观看人在逋| 亚洲狠狠婷婷综合久久图片| 国产精品九九99| 久久天堂一区二区三区四区| 免费av毛片视频| 国产av一区二区精品久久| 亚洲精品久久成人aⅴ小说| 美女大奶头视频| 男女午夜视频在线观看| 人人妻,人人澡人人爽秒播| 少妇粗大呻吟视频| 亚洲一区二区三区色噜噜| 亚洲精品在线美女| 亚洲精品一卡2卡三卡4卡5卡| 一本大道久久a久久精品| 精品卡一卡二卡四卡免费| 日日夜夜操网爽| 日韩大尺度精品在线看网址| 久久久久久久久免费视频了| 最近最新中文字幕大全电影3 | 两性午夜刺激爽爽歪歪视频在线观看 | 黄色女人牲交| 97人妻精品一区二区三区麻豆 | 国产av不卡久久| 亚洲av成人av| 亚洲熟妇熟女久久| 亚洲av成人av| 亚洲精品美女久久av网站| 99在线人妻在线中文字幕| 深夜精品福利| 宅男免费午夜| 国产精品,欧美在线| 午夜免费观看网址| 天堂动漫精品| 啦啦啦免费观看视频1| 午夜福利18| 91成人精品电影| 婷婷精品国产亚洲av在线| 国产野战对白在线观看| 亚洲国产中文字幕在线视频| 亚洲国产精品久久男人天堂| 九色国产91popny在线| 欧美久久黑人一区二区| 国产精品一区二区免费欧美| 怎么达到女性高潮| 久久久久久久久中文| 亚洲av片天天在线观看| 韩国av一区二区三区四区| 欧美乱色亚洲激情| 欧美性长视频在线观看| 久久国产精品影院| 国语自产精品视频在线第100页| 国产精品亚洲一级av第二区| 在线观看免费视频日本深夜| 91字幕亚洲| 日韩精品免费视频一区二区三区| 亚洲,欧美精品.| 叶爱在线成人免费视频播放| 欧美黑人巨大hd| 久久香蕉激情| 欧美一级毛片孕妇| 欧美一级a爱片免费观看看 | 国产片内射在线| 日日夜夜操网爽| 搞女人的毛片| 欧美日韩黄片免| 午夜福利欧美成人| 欧美+亚洲+日韩+国产| 午夜福利高清视频| 女同久久另类99精品国产91| 韩国精品一区二区三区| 亚洲专区中文字幕在线| 亚洲精品在线观看二区| 一a级毛片在线观看| 亚洲中文av在线| 99久久无色码亚洲精品果冻| 18禁黄网站禁片午夜丰满| 精品久久蜜臀av无| 中文在线观看免费www的网站 | 一级毛片精品| 欧美性长视频在线观看| 日本成人三级电影网站| 波多野结衣高清无吗| 一区福利在线观看| a级毛片a级免费在线| 免费观看人在逋| 两性夫妻黄色片| 婷婷丁香在线五月| 亚洲中文日韩欧美视频| 亚洲av五月六月丁香网| 日韩免费av在线播放| 日韩欧美一区视频在线观看| 免费高清在线观看日韩| 日韩高清综合在线| 夜夜夜夜夜久久久久| 国产高清videossex| 中文资源天堂在线| 国产视频一区二区在线看| 久久久国产欧美日韩av| 国产成年人精品一区二区| 亚洲九九香蕉| 国产精品 国内视频| aaaaa片日本免费| 国产aⅴ精品一区二区三区波| 麻豆一二三区av精品| 麻豆久久精品国产亚洲av| 男人舔女人下体高潮全视频| av有码第一页| 一本久久中文字幕| 成年女人毛片免费观看观看9| 正在播放国产对白刺激| 色播在线永久视频| 成人特级黄色片久久久久久久| 日本五十路高清| 中文字幕人妻丝袜一区二区| 久久久久久久精品吃奶| 首页视频小说图片口味搜索| 级片在线观看| 欧美国产精品va在线观看不卡| 麻豆久久精品国产亚洲av| 丁香欧美五月| 国产一区二区激情短视频| 免费搜索国产男女视频| 日本一区二区免费在线视频| 观看免费一级毛片| 亚洲五月色婷婷综合| 国产av在哪里看| 哪里可以看免费的av片| 免费搜索国产男女视频| 国产精品亚洲av一区麻豆| 看黄色毛片网站| 欧美又色又爽又黄视频| 琪琪午夜伦伦电影理论片6080| 亚洲成av人片免费观看| 欧美一级a爱片免费观看看 | 免费一级毛片在线播放高清视频| 国产亚洲精品久久久久久毛片| www日本黄色视频网| 国产精品av久久久久免费| 非洲黑人性xxxx精品又粗又长| 精品国产乱子伦一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 午夜影院日韩av| 啦啦啦 在线观看视频| 亚洲五月婷婷丁香| 国产精品久久电影中文字幕| 欧美 亚洲 国产 日韩一| 日韩免费av在线播放| 日韩欧美一区视频在线观看| 日韩欧美免费精品| 国产精品乱码一区二三区的特点| 天天一区二区日本电影三级| 悠悠久久av| 亚洲成人国产一区在线观看| 亚洲国产精品成人综合色| 国产日本99.免费观看| 啦啦啦观看免费观看视频高清| 天天一区二区日本电影三级| 亚洲成人精品中文字幕电影| 日韩视频一区二区在线观看| 亚洲成a人片在线一区二区| 高清在线国产一区| 国产精品99久久99久久久不卡| 99re在线观看精品视频| a级毛片a级免费在线| 麻豆国产av国片精品| 国产高清videossex| 在线十欧美十亚洲十日本专区| 嫩草影视91久久| 亚洲第一青青草原| 18美女黄网站色大片免费观看| 在线av久久热| 国产伦在线观看视频一区| 免费看a级黄色片| 免费观看人在逋| 狂野欧美激情性xxxx| 精品卡一卡二卡四卡免费| 老熟妇仑乱视频hdxx| 日本精品一区二区三区蜜桃| 中亚洲国语对白在线视频| 成人国语在线视频| 亚洲国产欧美网| 久久久久久人人人人人| 久久精品影院6| a在线观看视频网站| 免费看日本二区| 天天一区二区日本电影三级| 久久国产精品男人的天堂亚洲| 亚洲九九香蕉| 免费在线观看亚洲国产| 91九色精品人成在线观看| 好男人电影高清在线观看| 91成年电影在线观看| 三级毛片av免费| 国产99久久九九免费精品| 老汉色∧v一级毛片| 长腿黑丝高跟| 别揉我奶头~嗯~啊~动态视频| 我的亚洲天堂| 久久精品91蜜桃| 国产91精品成人一区二区三区| 欧美成人一区二区免费高清观看 | 12—13女人毛片做爰片一| 亚洲狠狠婷婷综合久久图片| 亚洲成av人片免费观看| 男男h啪啪无遮挡| 好男人电影高清在线观看| 女性生殖器流出的白浆| av福利片在线| 日日摸夜夜添夜夜添小说| 99久久久亚洲精品蜜臀av| 中文字幕人妻丝袜一区二区| 男人操女人黄网站| 久久国产亚洲av麻豆专区| 伊人久久大香线蕉亚洲五| 欧美性猛交黑人性爽| 99精品久久久久人妻精品| 夜夜躁狠狠躁天天躁| 999久久久精品免费观看国产| 国产精品二区激情视频| 国产精华一区二区三区| 久久午夜综合久久蜜桃| 午夜两性在线视频| 亚洲男人天堂网一区| 两性夫妻黄色片| 国产成人av激情在线播放| 国产亚洲精品久久久久久毛片| 一进一出好大好爽视频| 夜夜爽天天搞| 欧美中文综合在线视频| 又黄又爽又免费观看的视频| 国产不卡一卡二| 69av精品久久久久久| 欧美另类亚洲清纯唯美| 亚洲电影在线观看av| 国产一区二区在线av高清观看| 99热6这里只有精品| 日本熟妇午夜| av有码第一页| 精品第一国产精品| 国产精品久久久久久亚洲av鲁大| 午夜亚洲福利在线播放| 怎么达到女性高潮| 国产真人三级小视频在线观看| 日本 欧美在线| 夜夜躁狠狠躁天天躁| 波多野结衣高清无吗| 精品国产乱码久久久久久男人| 可以免费在线观看a视频的电影网站| 色老头精品视频在线观看| 国产av一区二区精品久久| 别揉我奶头~嗯~啊~动态视频| 欧美 亚洲 国产 日韩一| avwww免费| 无人区码免费观看不卡| 亚洲国产日韩欧美精品在线观看 | 国产高清视频在线播放一区| АⅤ资源中文在线天堂| 免费在线观看黄色视频的| 黑人操中国人逼视频| 国产亚洲精品第一综合不卡| 日日干狠狠操夜夜爽| 欧美另类亚洲清纯唯美| 成人特级黄色片久久久久久久| www.999成人在线观看| 免费在线观看亚洲国产| 国产成+人综合+亚洲专区| 久久精品夜夜夜夜夜久久蜜豆 | 日本 av在线| a在线观看视频网站| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产精品sss在线观看| 少妇裸体淫交视频免费看高清 | 最近最新中文字幕大全免费视频| 国产成人啪精品午夜网站| 国产成人欧美| 人人妻人人澡欧美一区二区| 日韩 欧美 亚洲 中文字幕| 母亲3免费完整高清在线观看| 国产在线精品亚洲第一网站| 国产精品久久久久久亚洲av鲁大| 亚洲国产精品成人综合色| 黄色毛片三级朝国网站| 男女做爰动态图高潮gif福利片| 香蕉国产在线看| 国产乱人伦免费视频| 午夜久久久在线观看| 亚洲av中文字字幕乱码综合 | 日韩大码丰满熟妇| 人人澡人人妻人| 午夜福利18| www.www免费av| 亚洲av片天天在线观看| 国产精品香港三级国产av潘金莲| 91国产中文字幕| 精品久久久久久久人妻蜜臀av| 国产精品久久久av美女十八| 日本 欧美在线| 脱女人内裤的视频| 久久人人精品亚洲av| 久久亚洲真实| 亚洲第一电影网av| 久久精品91蜜桃| 香蕉丝袜av| 国产aⅴ精品一区二区三区波| 丁香欧美五月| 日日干狠狠操夜夜爽| 淫妇啪啪啪对白视频| 搡老岳熟女国产| 欧美成狂野欧美在线观看| 少妇熟女aⅴ在线视频| 亚洲av电影在线进入| 日韩欧美国产在线观看| 国产精品久久久av美女十八| 日韩有码中文字幕| 亚洲 欧美 日韩 在线 免费| 中亚洲国语对白在线视频| 99国产精品一区二区三区| 国产精品日韩av在线免费观看| 亚洲精品国产一区二区精华液| 午夜福利高清视频| 日韩av在线大香蕉| 亚洲专区字幕在线| АⅤ资源中文在线天堂| 亚洲天堂国产精品一区在线| 国产男靠女视频免费网站| 色精品久久人妻99蜜桃| 九色国产91popny在线| 一区二区三区激情视频| 亚洲av电影在线进入| 精品久久久久久久久久久久久 | 久久久久精品国产欧美久久久| 精华霜和精华液先用哪个| 国产日本99.免费观看| 国产精品98久久久久久宅男小说| 日本免费a在线| 国产精品亚洲美女久久久| 变态另类丝袜制服| www.熟女人妻精品国产| 亚洲国产精品sss在线观看| 欧美中文日本在线观看视频| 男女午夜视频在线观看| 热99re8久久精品国产| 亚洲中文字幕一区二区三区有码在线看 | 热99re8久久精品国产| 一区二区三区高清视频在线| 午夜福利在线观看吧| 不卡av一区二区三区| 国产精品二区激情视频| 久久午夜综合久久蜜桃| 久久婷婷人人爽人人干人人爱| 色综合欧美亚洲国产小说| 波多野结衣av一区二区av| 99热6这里只有精品| 国产精品综合久久久久久久免费| 十八禁人妻一区二区| 日本成人三级电影网站| 最近最新免费中文字幕在线| 亚洲av五月六月丁香网| 久久精品国产亚洲av高清一级| 中文资源天堂在线| 久久久久久久午夜电影| 日日爽夜夜爽网站| 久久久水蜜桃国产精品网| 午夜免费鲁丝| 亚洲成av片中文字幕在线观看| 欧美久久黑人一区二区| 国产亚洲欧美在线一区二区| 岛国在线观看网站| 久久午夜亚洲精品久久| 在线视频色国产色| 精品第一国产精品| 叶爱在线成人免费视频播放| 亚洲人成网站在线播放欧美日韩| 非洲黑人性xxxx精品又粗又长| av在线播放免费不卡| 欧美一级a爱片免费观看看 | 久久久久国产一级毛片高清牌| 宅男免费午夜| 韩国av一区二区三区四区| av超薄肉色丝袜交足视频| 黄色 视频免费看| 亚洲成国产人片在线观看| 18禁国产床啪视频网站| 欧美日本视频| 99riav亚洲国产免费| 神马国产精品三级电影在线观看 | 国产精品98久久久久久宅男小说| 操出白浆在线播放| 男女下面进入的视频免费午夜 | 一区二区日韩欧美中文字幕| 黄色毛片三级朝国网站| 人成视频在线观看免费观看| 久久久久久人人人人人| 午夜福利成人在线免费观看| 变态另类丝袜制服| 婷婷亚洲欧美| 99精品久久久久人妻精品| 国产三级在线视频| 成人18禁高潮啪啪吃奶动态图| 国产精品爽爽va在线观看网站 | 免费看美女性在线毛片视频| 亚洲成人免费电影在线观看| 桃色一区二区三区在线观看| 色综合亚洲欧美另类图片| 久久久久久久久中文| 夜夜爽天天搞| 欧美+亚洲+日韩+国产| 国产激情久久老熟女| 久久精品影院6| 女性被躁到高潮视频| 国产亚洲av嫩草精品影院| 亚洲成人精品中文字幕电影| 黄片播放在线免费| а√天堂www在线а√下载| 午夜激情av网站| 女警被强在线播放| 久久99热这里只有精品18| 国产成人精品久久二区二区免费| 亚洲五月天丁香| www.自偷自拍.com| 天堂影院成人在线观看| 久久中文字幕一级| 国产一区二区三区视频了| 久久久久国产一级毛片高清牌| 亚洲熟妇中文字幕五十中出| 又黄又爽又免费观看的视频| 操出白浆在线播放| 黄网站色视频无遮挡免费观看| 人人妻人人澡人人看| 黄色女人牲交| 中文字幕久久专区| 免费人成视频x8x8入口观看| 伊人久久大香线蕉亚洲五| 成人18禁在线播放| 听说在线观看完整版免费高清| 亚洲三区欧美一区| 国产成人欧美在线观看| 亚洲精品美女久久av网站| 欧美不卡视频在线免费观看 | 9191精品国产免费久久| 黄片大片在线免费观看| 国产精品久久久av美女十八| 啦啦啦韩国在线观看视频| 欧美在线一区亚洲| 波多野结衣高清无吗| 91大片在线观看| 国产伦人伦偷精品视频| 丝袜在线中文字幕| 老司机福利观看| 亚洲成人国产一区在线观看| 中文字幕久久专区| 亚洲性夜色夜夜综合| 国产精品国产高清国产av| 精品久久久久久久末码| 国产99白浆流出| 亚洲欧美激情综合另类| 亚洲第一av免费看| 91大片在线观看| 久久久久久国产a免费观看| 亚洲欧美日韩高清在线视频| 90打野战视频偷拍视频| 欧美绝顶高潮抽搐喷水| 男女视频在线观看网站免费 | 大型黄色视频在线免费观看| 国产日本99.免费观看| 一级a爱片免费观看的视频| 午夜成年电影在线免费观看| 18禁国产床啪视频网站| 黄片播放在线免费| 变态另类丝袜制服| 精品高清国产在线一区| 一二三四社区在线视频社区8| 国产成人av激情在线播放| 99在线人妻在线中文字幕| 国产精品野战在线观看| 哪里可以看免费的av片| 精品国产乱子伦一区二区三区| 天堂动漫精品| 黄片播放在线免费| 国产成年人精品一区二区| 欧美一级a爱片免费观看看 | 亚洲第一av免费看| 亚洲第一欧美日韩一区二区三区| 国产免费av片在线观看野外av| 婷婷丁香在线五月| 亚洲精品一区av在线观看| 久99久视频精品免费| 午夜福利在线在线| 女人被狂操c到高潮| 午夜福利欧美成人| svipshipincom国产片| 国产亚洲精品第一综合不卡| 真人一进一出gif抽搐免费| 免费在线观看影片大全网站| 十八禁人妻一区二区| 手机成人av网站|