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

    虹吸管路內(nèi)氣泡動(dòng)力學(xué)行為數(shù)值模擬

    2020-12-30 12:55:30李興雨李琳譚義海
    人民黃河 2020年12期
    關(guān)鍵詞:虹吸管數(shù)值模擬

    李興雨 李琳 譚義海

    摘 要:為了研究虹吸管中氣泡運(yùn)動(dòng)行為與變形特征,對(duì)虹吸管中氣泡運(yùn)動(dòng)過程進(jìn)行數(shù)值模擬,結(jié)果表明:上行管內(nèi)氣泡縱橫比越小,氣泡速度越大,流量的變化對(duì)氣泡縱橫比影響很小,且氣泡形狀與靜水中氣泡相圖相吻合,氣泡在上升過程中會(huì)發(fā)生橫向振蕩,振蕩幅度隨著氣泡直徑的增大先增大后減小;下行管內(nèi)流量不同,氣泡運(yùn)動(dòng)方向也不同,氣泡變形情況與上行管內(nèi)一致;中行管內(nèi)氣泡為貼壁狀態(tài)流動(dòng),隨著氣泡直徑增大,氣泡速度減小,當(dāng)氣泡直徑增加到一定數(shù)值時(shí),氣泡會(huì)發(fā)生破裂,出現(xiàn)較小的氣泡,氣泡速度又出現(xiàn)增大的現(xiàn)象。

    關(guān)鍵詞:虹吸管;數(shù)值模擬;氣泡速度;縱橫比;振蕩幅度

    中圖分類號(hào):TV672.5 文獻(xiàn)標(biāo)志碼:A

    doi:10.3969/j.issn.1000-1379.2020.12.032

    Abstract:The aim of this study was to detect the movement behavior and deformation characteristics of the air bubbles in siphon pipeline. It was essential to use numerical methods to figure out the movement of the air bubbles in the siphon pipeline. The results show that the bubble aspect ratio in the ascending tube is smaller and the bubble velocity is bigger. The change of flow has little effect on the aspect ratio and the bubble shape is consistent with its phase diagram in still water. During the progress of rising, the bubble oscillates laterally and the amplitude of oscillation first increases and then decreases with the increase of bubble diameter. The flow shows a wide difference in the downpipe, the movement direction of the bubble is also different and the bubble deformation is in agreement with the ascending tube. The bubbles in the middle pipe flow in the way with attaching the sidewall. Meanwhile, the bubble diameter increases with the decreasing of the bubble velocity. The bubble breaks at a certain diameter and smaller bubbles appear, then bubble velocity increases again. The results can provide a reference for the study of gas-liquid two-phase flow in siphon pipeline.

    Key words: siphon; numerical methods; bubble velocity; aspect ratio; oscillation amplitude

    水平管段距離長(zhǎng)、真空度大的虹吸管輸水管道是坎兒井式地下水庫(kù)的重要組成部分。李琳等[1]、許史等[2-3]在新疆臺(tái)蘭河某地下水庫(kù)長(zhǎng)距離虹吸管道水力學(xué)模型試驗(yàn)中發(fā)現(xiàn),安裝高度小于7 m時(shí),虹吸管內(nèi)出現(xiàn)偽空化現(xiàn)象,流動(dòng)介質(zhì)由單一液相流轉(zhuǎn)變?yōu)闅馀萘鳌⑦^渡流和氣團(tuán)流。王夢(mèng)婷[4]、張小瑩等[5]的試驗(yàn)結(jié)果表明,虹吸管安裝高度為6~8 m時(shí),各水頭運(yùn)行時(shí)上行管均以氣泡流為主,進(jìn)入水平管路后,低水頭運(yùn)行時(shí)氣泡在運(yùn)動(dòng)過程中不斷聚合形成氣團(tuán)和長(zhǎng)度為30~60 cm的大氣囊,氣囊體積隨運(yùn)行時(shí)間增加而增大,當(dāng)其貫穿整個(gè)斷面時(shí)管路斷流。為了保證虹吸管正常運(yùn)行,應(yīng)盡量避免氣泡在一定條件下聚合成氣囊,而這一問題與氣泡在動(dòng)水和負(fù)壓運(yùn)行的管道內(nèi)的動(dòng)力學(xué)行為密切相關(guān)。

    氣泡以分散相的形式在水中的運(yùn)動(dòng)過程是一種復(fù)雜、不穩(wěn)定和非線性的水動(dòng)力學(xué)現(xiàn)象。Sanada等[6]、Wang等[7]通過試驗(yàn)研究和數(shù)值模擬對(duì)不同黏度液體中的氣泡運(yùn)動(dòng)進(jìn)行了研究,得到了不同直徑下氣泡的形狀。程軍明等[8]對(duì)氣泡在靜水中上升和破裂過程進(jìn)行了數(shù)值模擬,獲得了氣泡上升速度和氣泡形狀之間的關(guān)系。陳如艷等[9]通過試驗(yàn)研究得到水平管內(nèi)氣液兩相流中液相流速、氣相流速和氣泡尺寸的相互影響規(guī)律。Bhaga等[10]研究了浮力氣泡在黏性液體中的上升運(yùn)動(dòng)過程,得到了氣泡形狀及終端上升速度與雷諾數(shù)Re、莫頓數(shù)Mo和埃奧特沃斯數(shù)Eo之間關(guān)系的氣泡相圖。Ellingsen等[11]使用高速攝像機(jī)對(duì)氣泡在各種局部流動(dòng)條件及不同水質(zhì)情況下的動(dòng)力學(xué)特性進(jìn)行了研究,并得到了在不同條件下氣泡的運(yùn)動(dòng)及形變過程。Ohta等[12-13]先后運(yùn)用VOF數(shù)值模型,揭示了當(dāng)流動(dòng)條件為低Mo和Eo時(shí),氣泡上升行為與其初始形狀密切相關(guān)。前人的研究主要集中于靜水中氣泡的運(yùn)動(dòng)特性,而對(duì)于動(dòng)水和負(fù)壓運(yùn)行下虹吸管內(nèi)的氣泡運(yùn)動(dòng)特性尚未見報(bào)道。本文應(yīng)用標(biāo)準(zhǔn)k-ε模型,同時(shí)結(jié)合VOF方法對(duì)虹吸管上、中、下行管路中單個(gè)氣泡在不同流量條件下的運(yùn)動(dòng)行為進(jìn)行模擬,為探明虹吸管內(nèi)氣泡聚合形成氣囊和氣泡獨(dú)立運(yùn)動(dòng)產(chǎn)生條件奠定基礎(chǔ),同時(shí)為實(shí)際工程中合理布置虹吸管路排氣裝置、保證管路正常運(yùn)行提供設(shè)計(jì)依據(jù)。

    1 湍流數(shù)學(xué)模型

    1.1 控制方程

    1.2 標(biāo)準(zhǔn)k-ε模型和VOF法

    在整個(gè)虹吸管內(nèi),除了彎頭處以外,其他地方流線基本相互平行,因此選擇計(jì)算量小、計(jì)算精度和收斂性均較好的標(biāo)準(zhǔn)k-ε模型進(jìn)行計(jì)算。湍動(dòng)能k和耗散率ε方程為式中:σk、σε分別為湍動(dòng)能k和湍動(dòng)能耗散率ε對(duì)應(yīng)的Prandtl數(shù),均為1.39;t為時(shí)間;Gk為湍動(dòng)能k的生成項(xiàng),Gk=μt(uixj+ujxj)uixj,其中μt=ρCμk2ε,Cμ=0.084 5;C1ε=1.44;C2ε=1.92;C3ε=0.09;Gb是因浮力影響引起的湍動(dòng)能而產(chǎn)生的,Gb=βgiμtσtuixi,β=0.012,σt為紊動(dòng)普朗特?cái)?shù),其他參數(shù)取值見文獻(xiàn)[15];YM為有壓縮流脈動(dòng)膨脹對(duì)總耗散率的影響值,在本文中模擬的液相認(rèn)為是不可壓縮的,因此不考慮該項(xiàng)。

    精確描述氣液兩相的運(yùn)動(dòng)界面是研究氣泡行為的關(guān)鍵。VOF方法追蹤的是網(wǎng)格中流體體積,具有容易實(shí)現(xiàn)、計(jì)算量小和精度高等優(yōu)點(diǎn),因此本文采用VOF法追蹤氣液交界面。當(dāng)αw=0時(shí),表示管道內(nèi)沒有水,被氣充滿;當(dāng)水的體積分?jǐn)?shù)αw=1時(shí),表示管內(nèi)被水充滿;當(dāng)0≤αw≤1時(shí),表示管內(nèi)被水和氣充滿。αw的控制方程為

    2 計(jì)算區(qū)域的離散及邊界條件

    由于本文模擬的管道為圓形截面,氣泡為球形,管道與氣泡都具有軸對(duì)稱性,因此建立二維模型進(jìn)行數(shù)值模擬。模型由上行管、中行管和下行管三部分構(gòu)成,其中上行管和下行管高度均為6 m,中行管長(zhǎng)度為18 m,管徑為2 cm。分別對(duì)單個(gè)氣泡在上行管、中行管、下行管內(nèi)運(yùn)動(dòng)過程進(jìn)行數(shù)值模擬研究。虹吸管模型邊界采用速度進(jìn)口邊界條件,液相初始速度分別為0.06、0.15、0.2 m/s且垂直于管路的橫斷面,氣相的初始速度為0。管道出口為壓力邊界,相對(duì)壓強(qiáng)為0,管道壁面設(shè)為無滑移固壁,采用標(biāo)準(zhǔn)壁面函數(shù)進(jìn)行修正。當(dāng)t=0時(shí),氣泡位于上行管距離管口50 cm中心處。計(jì)算區(qū)域網(wǎng)格劃分采用四邊形網(wǎng)格和三角形網(wǎng)格,由于整體劃分網(wǎng)格數(shù)目較多,計(jì)算時(shí)間過長(zhǎng),因此對(duì)上行管、中行管和下行管網(wǎng)格局部加密。氣泡在不同管段運(yùn)動(dòng)時(shí),需運(yùn)動(dòng)一段距離后速度才能達(dá)到穩(wěn)定。在上行管中,對(duì)氣泡起點(diǎn)以上3 m范圍內(nèi)的管道進(jìn)行加密;當(dāng)氣泡在中行管運(yùn)動(dòng)時(shí),對(duì)距離氣泡起點(diǎn)0.5 m處之后的3 m范圍內(nèi)的管道進(jìn)行加密;下行管中,由于流量不同,氣泡運(yùn)動(dòng)方向也不同,因此針對(duì)不同流量分別做了加密方案,流量較小時(shí),對(duì)氣泡起點(diǎn)以上3 m范圍內(nèi)的管道進(jìn)行加密,流量較大時(shí),對(duì)氣泡起點(diǎn)以下3 m范圍內(nèi)的管道進(jìn)行加密;網(wǎng)格單元最小尺寸為4×10-4m??刂品匠痰碾x散采用有限體積法,從穩(wěn)定性、精確性、適用性方面考慮,控制方程中的對(duì)流項(xiàng)和擴(kuò)散項(xiàng)均采用二階迎風(fēng)格式進(jìn)行離散。離散后的線性代數(shù)方程組采用SIMPLEC算法迭代求解,計(jì)算時(shí)間步長(zhǎng)為1×10-4 s。

    使用標(biāo)準(zhǔn)的k-ε模型和VOF法,通過FLUENT軟件進(jìn)行模擬,設(shè)置表面張力系數(shù)σ=0.072 8 N/m,模擬虹吸管輸水流量分別為0.018 8、0.047 1 L/s,氣泡直徑分別為1、4、5、6、7 mm的管路內(nèi)氣泡運(yùn)動(dòng)行為特征。

    3 計(jì)算結(jié)果與分析

    3.1 數(shù)學(xué)模型的驗(yàn)證

    采用文獻(xiàn)[16]和文獻(xiàn)[17]的試驗(yàn)結(jié)果對(duì)本文的數(shù)學(xué)模型及經(jīng)驗(yàn)參數(shù)進(jìn)行驗(yàn)證。

    3.1.1 物理模型試驗(yàn)裝置

    文獻(xiàn)[16]試驗(yàn)裝置見圖1。試驗(yàn)裝置高度為1 m,長(zhǎng)、寬分別為20.32、7.56 cm,液面高度為0.8 m。使用空氣注射器注入所需大小的氣泡,通過水箱與管道的壓力差使得氣泡進(jìn)入試驗(yàn)裝置。氣泡出口位置為半球形倒置的杯罩,手動(dòng)旋轉(zhuǎn)杯罩釋放氣泡,使用高速攝像機(jī)在測(cè)量區(qū)域進(jìn)行測(cè)量。

    3.1.2 數(shù)學(xué)模型的建立與對(duì)比

    文獻(xiàn)[18]中指出計(jì)算區(qū)域的寬度大于4倍氣泡直徑時(shí),可忽略邊界對(duì)氣泡運(yùn)動(dòng)特性的影響。為了驗(yàn)證擬采用的數(shù)學(xué)模型,先建立高度為1 m、直徑為50 mm的一端封閉一端開口的豎直管道模型,出口邊界設(shè)為壓力出口(相對(duì)壓強(qiáng)為0),網(wǎng)格尺寸分別為6×10-4、5×10-4、4×10-4、3×10-4 m,分別模擬了直徑為4、6、8 mm的氣泡在靜水中的上升和變形過程,氣泡直徑為6 mm時(shí)不同網(wǎng)格尺寸下氣泡速度隨時(shí)間變化情況見圖2。由圖2可知,網(wǎng)格尺寸越小,氣泡速度越快,當(dāng)網(wǎng)格尺寸小于4×10-4 m時(shí),氣泡速度不再隨著網(wǎng)格尺寸的變化而發(fā)生變化。網(wǎng)格尺寸越小,需要的網(wǎng)格數(shù)量越大,計(jì)算耗時(shí)就越長(zhǎng),為了減少計(jì)算時(shí)間,選擇網(wǎng)格尺寸為4×10-4m進(jìn)行計(jì)算。由圖2還可得到,網(wǎng)格尺寸為4×10-4m時(shí)模擬的氣泡速度與實(shí)測(cè)的結(jié)果吻合較好。當(dāng)t=0.15 s時(shí),氣泡速度模擬結(jié)果與試驗(yàn)結(jié)果相差最小,相對(duì)誤差為1.7%;當(dāng)t=0.25 s時(shí),氣泡速度模擬結(jié)果與試驗(yàn)結(jié)果相差最大,相對(duì)誤差為4.9%。

    將氣泡直徑為6 mm時(shí)的計(jì)算結(jié)果與文獻(xiàn)[17]中實(shí)測(cè)的氣泡相圖及文獻(xiàn)[16]的試驗(yàn)結(jié)果進(jìn)行對(duì)比,結(jié)果如表1和圖2所示。表1為不同直徑氣泡在上升過程中不同時(shí)刻所對(duì)應(yīng)的Re和Eo,并將氣泡形狀模擬結(jié)果與文獻(xiàn)[17]結(jié)果進(jìn)行對(duì)比。從表1中可以看出3種直徑的氣泡運(yùn)動(dòng)過程與文獻(xiàn)[17]中氣泡相圖的結(jié)果一致。文獻(xiàn)[19]中提到,氣泡在流體中上升時(shí),氣泡形狀主要取決于氣泡雷諾數(shù)Re、埃奧特沃斯數(shù)Eo和莫頓數(shù)Mo三個(gè)無因次準(zhǔn)數(shù),它們的表達(dá)式分別為

    3.2 氣泡豎直運(yùn)動(dòng)和水平運(yùn)動(dòng)過程及形狀變化

    在負(fù)壓條件下,不同流量時(shí)虹吸管上、中、下行管內(nèi)直徑d分別為5、6、7 mm的氣泡運(yùn)動(dòng)過程及形狀變化見圖3~圖5。氣泡形變過程用不同時(shí)刻的縱橫比(即氣泡短軸與長(zhǎng)軸的比值)變化來表示,t=0時(shí),縱橫比為1,氣泡縱橫比越接近1表明氣泡變形量越小。

    圖3為上行管內(nèi)氣泡運(yùn)動(dòng)變形過程,t=0時(shí)氣泡為球形。當(dāng)流量為0.018 8 L/s時(shí),隨著氣泡上升,氣泡變成橢球狀,隨著直徑的增大,變形越來越明顯,如d=5 mm的氣泡經(jīng)過4 s后,縱橫比由t=0時(shí)的1變?yōu)?.64,而d=7 mm的氣泡經(jīng)過相同時(shí)間縱橫比由t=0時(shí)的1變?yōu)?.55。當(dāng)流量由0.018 8 L/s增大至0.047 1 L/s時(shí),氣泡仍然保持橢球狀上升,縱橫比隨氣泡直徑增大而減小,但是流量的增大對(duì)氣泡形狀的變化影響很小,如d=5 mm和d=7 mm的氣泡經(jīng)過6 s后縱橫比分別由1變?yōu)?.63和0.56,與流量為0.018 8 L/s時(shí)的縱橫比基本相同。

    圖4為球形氣泡在下行管隨水流運(yùn)動(dòng)時(shí)氣泡變形過程。氣泡在下行或上行過程中主要經(jīng)歷了球形和橢球形的變形過程。計(jì)算不同時(shí)刻的Re、Eo和Mo,與文獻(xiàn)[17]的氣泡相圖進(jìn)行對(duì)比,結(jié)果表明:動(dòng)水條件下虹吸管上、下行管內(nèi)氣泡在不同時(shí)刻的形狀變化過程模擬結(jié)果與文獻(xiàn)[17]的靜水條件下氣泡運(yùn)動(dòng)變形過程一致。當(dāng)流量為0.018 8 L/s時(shí),因流速較小,流體曳力(物體在流體中有相對(duì)運(yùn)動(dòng)時(shí),會(huì)受到流體的阻力,阻力大小由相對(duì)速度差決定)較小,在浮力作用下d=5 mm的氣泡在下行管內(nèi)向上運(yùn)動(dòng),氣泡運(yùn)動(dòng)速度小于上行管中氣泡運(yùn)動(dòng)速度,此時(shí)下行管氣泡縱橫比較大,氣泡運(yùn)動(dòng)過程中變形量略小于上行管內(nèi)的。如t=4 s時(shí),d=5 mm的氣泡在下行管內(nèi)氣泡縱橫比比上行管內(nèi)增加2%。當(dāng)流量為0.047 1 L/s時(shí),流體曳力大于浮力,氣泡隨水流一起向下運(yùn)動(dòng),隨著氣泡速度的增大,相對(duì)速度差減小,流體曳力減小,下行管氣泡運(yùn)動(dòng)速度小于上行管內(nèi)的,因此氣泡變形量仍小于上行管內(nèi)的氣泡變形量,經(jīng)過4 s后,d=7 mm的氣泡縱橫比比上行管中增大了9%。

    圖5為氣泡在中行管運(yùn)動(dòng)變形過程,不同流量下氣泡變形過程基本一致。當(dāng)t=0時(shí),氣泡為球形。當(dāng)t>0時(shí),在浮力和壁面黏附力作用下氣泡變?yōu)榘霗E球狀,自管道中心向上管壁運(yùn)動(dòng),保持貼壁流運(yùn)動(dòng)。當(dāng)Q=0.018 8 L/s時(shí),經(jīng)過相同時(shí)間,氣泡縱橫比隨氣泡直徑的增大而減小。如d=5 mm和7 mm的氣泡經(jīng)過4 s后縱橫比分別由1變?yōu)?.24和0.17,對(duì)同一直徑氣泡而言,經(jīng)過相同時(shí)間,氣泡縱橫比隨流量增大而減小。如Q=0.018 8 L/s時(shí),d=5 mm的氣泡經(jīng)過4 s后縱橫比較Q=0.047 1 L/s時(shí)增大了11%。

    3.3 虹吸管內(nèi)不同直徑氣泡速度與縱橫比

    圖6為Q=0.018 8 L/s時(shí)不同直徑氣泡的速度、縱橫比隨時(shí)間的變化情況。其中:v、v1、v2分別為氣泡直徑為5、6、7 mm時(shí)的氣泡運(yùn)動(dòng)速度,e、e1、e2分別為氣泡直徑為5、6、7 mm時(shí)的氣泡縱橫比。

    氣泡直徑一定時(shí),縱橫比越大,氣泡速度越小,這一規(guī)律與文獻(xiàn)[16]中提到的靜水條件下氣泡運(yùn)動(dòng)過程中縱橫比與其速度成反比的規(guī)律一致。流量增大時(shí),氣泡縱橫比與速度的關(guān)系和小流量時(shí)一致(限于篇幅,文中未給出其他流量時(shí)的圖形)。從圖6(a)可以看出,氣泡直徑越大,速度越小,但不同直徑氣泡的速度相差不大,如d=5 mm和d=7 mm時(shí),最大差為t=1.0 s時(shí)的12%,最小差為t=2 s時(shí)的2%。其原因是氣泡直徑增大,所受浮力增大,在浮力作用下,氣泡與液體之間的相對(duì)速度差增大,導(dǎo)致流體曳力增大,抑制氣泡上升,隨著氣泡直徑的增大,管壁黏附力將會(huì)阻礙氣泡運(yùn)動(dòng),使得大、小氣泡的運(yùn)動(dòng)速度差別不明顯。

    從圖6(b)可以看出,不同直徑的氣泡在下行管內(nèi)的流速均小于上行管內(nèi)的,且同一時(shí)刻氣泡縱橫比越大,速度越小。Q=0.018 8 L/s時(shí),在下行管中,氣泡向上運(yùn)動(dòng),由于水流方向與上行管內(nèi)水流方向相反,因此速度比上行管內(nèi)小。氣泡速度隨著氣泡直徑的增大略有減小,如氣泡直徑由5 mm增加到7 mm時(shí),氣泡速度由0.16 m/s減小19%,其原因是較大的氣泡所受管壁黏附作用影響較大,使得阻力增大。當(dāng)流量增加到0.047 1 L/s時(shí),流體曳力增大,氣泡向下運(yùn)動(dòng),氣泡的縱橫比與速度的關(guān)系也與上述一致;氣泡直徑對(duì)氣泡速度的影響很小。

    從圖6(c)可以看出,在中行管內(nèi)不同直徑的氣泡先進(jìn)行短暫的減速運(yùn)動(dòng)后達(dá)到終速度,與上行管相比速度有所減小,其原因是在中行管內(nèi)氣泡先向上運(yùn)動(dòng),直到與管壁貼合,壁面黏附力增大,浮力不再做功,氣泡速度減小。當(dāng)直徑d=5~6 mm時(shí),氣泡速度會(huì)隨著氣泡直徑的增大而減小。當(dāng)氣泡直徑進(jìn)一步增大,如d=7 mm時(shí),氣泡速度增大,其原因是氣泡直徑過大,運(yùn)動(dòng)過程中發(fā)生了破裂,形成了較小的氣泡,使得速度又出現(xiàn)增大的現(xiàn)象。由于中行管內(nèi)是貼壁流,因此氣泡縱橫比與上行管相比相差很大,如Q=0.018 8 L/s時(shí),d=5 mm的氣泡經(jīng)過4 s后,中行管內(nèi)氣泡縱橫比較上行管內(nèi)減小65%。

    3.4 氣泡運(yùn)動(dòng)軌跡

    氣泡在上升過程中會(huì)發(fā)生橫向位移,其原因是氣泡在上升的過程中會(huì)發(fā)生形變使得受力不平衡,導(dǎo)致了橫向振蕩。將氣泡中心點(diǎn)的位置坐標(biāo)(0,0.6)作為t=0時(shí)刻氣泡的運(yùn)動(dòng)位置,流量為0.047 1 L/s時(shí)不同直徑的氣泡在上升過程的運(yùn)動(dòng)軌跡見圖7,圖中x/d表示氣泡的中心點(diǎn)位置沿管徑方向的相對(duì)位置,y/hs表示氣泡的中心點(diǎn)位置與虹吸管高度的相對(duì)距離。由圖7可知,氣泡在上行管中運(yùn)動(dòng)時(shí),隨著氣泡直徑的增大,振蕩幅度發(fā)生變化。當(dāng)上升氣泡直徑較小時(shí),如d=1 mm時(shí),氣泡上升中的振蕩幅度較小,基本沿直線上升。隨著氣泡直徑的增大,如d=4 mm時(shí),氣泡呈“之”字形上升。當(dāng)氣泡直徑增大到7 mm時(shí),氣泡的振蕩幅度卻開始減小,其原因是直徑較大的氣泡在上升過程中受到壁面黏滯作用的影響增大,對(duì)氣泡的振蕩起阻礙作用,使得振蕩幅度變小。在不同流量的情況下,氣泡的運(yùn)動(dòng)軌跡呈現(xiàn)相似的變化。

    4 結(jié) 論

    本研究模擬了單個(gè)氣泡在靜水中的上升過程,并與文獻(xiàn)[17]的氣泡相圖結(jié)果進(jìn)行比較,驗(yàn)證了模擬結(jié)果的正確性,在此基礎(chǔ)上對(duì)虹吸管道水流進(jìn)行數(shù)值模擬研究,并對(duì)模擬結(jié)果進(jìn)行分析。

    (1)氣泡在上行管上升過程中,由球狀變?yōu)闄E球狀,縱橫比逐漸減小,氣泡直徑越大,縱橫比變化越劇烈,而流量的變化對(duì)氣泡縱橫比的影響很小。氣泡速度與縱橫比有關(guān),縱橫比越大,氣泡速度越小。氣泡直徑較小時(shí),氣泡沿直線上升,隨著氣泡直徑增大,氣泡發(fā)生振蕩,呈“之”字形上升。當(dāng)氣泡直徑增大到一定程度時(shí),氣泡受到壁面的黏滯作用影響增大,對(duì)氣泡的振蕩起阻礙作用,使得振蕩幅度變小。

    (2)氣泡在下行管運(yùn)動(dòng)過程中,氣泡變形情況與上行管相似,但氣泡受到浮力和曳力作用,下行管內(nèi)氣泡速度和變形程度小于上行管內(nèi)的。

    (3)氣泡在中行管運(yùn)動(dòng)過程中,氣泡經(jīng)過短暫時(shí)間運(yùn)動(dòng)至管道頂部,在此之后為貼壁運(yùn)動(dòng),氣泡由球狀變?yōu)榘霗E球狀。隨著氣泡直徑的增大,氣泡縱橫比逐漸減小;縱橫比大小與流量有關(guān),流量越大,縱橫比越小。氣泡速度隨氣泡直徑增大而逐漸減小,氣泡發(fā)生破裂后的氣泡速度大于原來氣泡的速度。

    參考文獻(xiàn):

    [1] 李琳,邱秀云,許史,等.長(zhǎng)距離虹吸管道輸水水力學(xué)模型試驗(yàn)研究[J].南水北調(diào)與水利科技,2010,8(3):106-109.

    [2] 許史,李琳,邱秀云,等.長(zhǎng)距離虹吸管輸水試驗(yàn)研究初探[J].中國(guó)農(nóng)村水利水電,2010(3):70-72.

    [3] 許史.長(zhǎng)距離虹吸管輸水試驗(yàn)研究[D].烏魯木齊:新疆農(nóng)業(yè)大學(xué),2010:13-24.

    [4] 王夢(mèng)婷.真空管道氣液兩相流動(dòng)水力特性試驗(yàn)研究[D].烏魯木齊:新疆農(nóng)業(yè)大學(xué),2015:19-69.

    [5] 張小瑩,李琳,譚義海,等.虹吸管道坡度對(duì)氣液兩相流動(dòng)特性影響的試驗(yàn)研究[J].農(nóng)業(yè)工程學(xué)報(bào),2017,33(14):122-129.

    [6] SANADA T,SUGIHARA K,SHIROTA M,et al. Motion and Drag of a Single Bubble in Super-Purified Water [J]. Fluid Dynamics Research,2008,40(7/8):534-545.

    [7] WANG Huanran,LI Yang,YANG Dong,et al. On the Shape Feature of a Single Bubble Rising in Viscous Liquids [J]. Engineering Thermophysics,2009,30(9):1492-1494.

    [8] 程軍明,吳偉烽,聶娟,等.氣泡在靜水中上升破裂產(chǎn)生射流特性的數(shù)值模擬[J].江蘇大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,35(5):513-517.

    [9] 陳如艷,仇性啟,王曉琦.水平管內(nèi)氣液兩相流研究[J].石油化工設(shè)備,2007,36(6):37-38.

    [10] BHAGA D,WEBER M E. Bubbles in Viscous Liquids: Shapes,Wake and Velocities [J]. Fluid Mech, 1981,105:61-85.

    [11] ELLINGSEN K,RISSO F. On the Rise of an Ellipsoidal Bubble in Water:Oscillatory Paths and Liquid Induced Velocity [J]. Computer Phys, 2008,227:3358-3382.

    [12] OHTA M,IMURA T, YOSHIDA Y, et al. A Computational Study of the Effect of Initial Bubble Conditions on the Motion of a Gas Bubble Rising in Viscous Liquids [J]. Multiphase Flow,2005,31(2):223-237.

    [13] BONOMETTI T,MAGNAUDET J. Transition from Spherical Cap to Toroidal Bubbles [J]. Phys Fluids,2006,18(5):1-12.

    [14] BRACKBILL J,KOTHE D. A Continuum Method for Modeling Surface Tension [J]. Journal of Computational Phvsics,1992,100:335-345.

    [15] ISSA R I,OLIVEIRA. Numerical Prediction of Phase Separation in TwoPhase Flow Through TJunctions [J]. Computers Fluids,1994,23(2):347-372.

    [16] 劉柳.垂直上升管中氣泡動(dòng)力學(xué)特性實(shí)驗(yàn)研究[D].長(zhǎng)沙:中南大學(xué),2013:13-17.

    [17] 李仲春,宋小明,姜?jiǎng)僖?靜水中較大氣泡運(yùn)動(dòng)特性實(shí)驗(yàn)研究[J].核動(dòng)力工程,2015,36(1):161-164.

    [18] 何丹,李彥鵬,劉艷艷.初始形狀對(duì)浮升氣泡動(dòng)力特性的影響[J].西安交通大學(xué)學(xué)報(bào),2011,45(1):43-47.

    [19] 劉靜如.非牛頓流體中多氣泡相互作用、聚并與破裂過程的數(shù)值模擬[D].天津:天津大學(xué),2014:1-10.

    【責(zé)任編輯 許立新】

    猜你喜歡
    虹吸管數(shù)值模擬
    虹吸管在楊溝水庫(kù)中的應(yīng)用分析
    陜西水利(2022年12期)2022-11-30 06:09:58
    虹吸管安裝高度問題淺析
    虹吸管在水庫(kù)高水頭小流量泄水狀態(tài)下的應(yīng)用
    熱虹吸管換熱器應(yīng)用于除濕系統(tǒng)的節(jié)能性分析
    太陽(yáng)能動(dòng)力虹吸管在小型水利工程中的應(yīng)用
    張家灣煤礦巷道無支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場(chǎng)數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    亚洲天堂国产精品一区在线| 淫秽高清视频在线观看| 亚洲专区中文字幕在线| 超碰成人久久| 色视频www国产| 这个男人来自地球电影免费观看| 免费av毛片视频| 18禁观看日本| 9191精品国产免费久久| 久久这里只有精品19| 国产黄a三级三级三级人| 亚洲一区二区三区色噜噜| 国产免费男女视频| 亚洲av成人精品一区久久| 国产97色在线日韩免费| 在线观看午夜福利视频| 俄罗斯特黄特色一大片| 一区福利在线观看| 成人国产综合亚洲| 99热精品在线国产| 校园春色视频在线观看| 热99re8久久精品国产| 国产亚洲精品久久久com| 免费人成视频x8x8入口观看| 色综合站精品国产| av在线蜜桃| 欧美日韩精品网址| 国产探花在线观看一区二区| 校园春色视频在线观看| 操出白浆在线播放| 久9热在线精品视频| 狠狠狠狠99中文字幕| 极品教师在线免费播放| 午夜福利免费观看在线| 亚洲av成人av| 国产野战对白在线观看| 欧美日韩黄片免| 亚洲av五月六月丁香网| 欧美日本视频| 久久久久久久久中文| 精品午夜福利视频在线观看一区| 欧美日韩一级在线毛片| 99国产综合亚洲精品| 中文字幕人成人乱码亚洲影| 国产伦精品一区二区三区视频9 | 色噜噜av男人的天堂激情| 999久久久国产精品视频| 亚洲 欧美 日韩 在线 免费| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕最新亚洲高清| 国产v大片淫在线免费观看| 成熟少妇高潮喷水视频| 久久精品国产99精品国产亚洲性色| 日本 欧美在线| 中文字幕熟女人妻在线| 最近视频中文字幕2019在线8| 99在线人妻在线中文字幕| 成年女人永久免费观看视频| 免费高清视频大片| 男女午夜视频在线观看| 久久久久国产精品人妻aⅴ院| 免费观看精品视频网站| 女人被狂操c到高潮| 久久久国产精品麻豆| 久久久久久九九精品二区国产| 两人在一起打扑克的视频| 日日干狠狠操夜夜爽| 久久久久久大精品| 美女午夜性视频免费| 激情在线观看视频在线高清| 精品久久久久久久久久免费视频| 最近在线观看免费完整版| 国产av不卡久久| 亚洲精品一区av在线观看| 国产精品久久久久久亚洲av鲁大| 国产精品亚洲美女久久久| 久久久久精品国产欧美久久久| 一进一出抽搐动态| 夜夜爽天天搞| 国产精品香港三级国产av潘金莲| 两个人的视频大全免费| 亚洲 欧美 日韩 在线 免费| 曰老女人黄片| 中文字幕av在线有码专区| 夜夜躁狠狠躁天天躁| 欧美中文综合在线视频| 国产精品亚洲av一区麻豆| 很黄的视频免费| 别揉我奶头~嗯~啊~动态视频| www.www免费av| 免费观看人在逋| 91麻豆精品激情在线观看国产| 在线观看一区二区三区| 91字幕亚洲| 久久久久久九九精品二区国产| 国产一区二区三区视频了| 99视频精品全部免费 在线 | 在线视频色国产色| 波多野结衣高清无吗| 日韩免费av在线播放| 日韩人妻高清精品专区| 久久国产精品影院| 亚洲av五月六月丁香网| 亚洲中文日韩欧美视频| 真实男女啪啪啪动态图| 一本精品99久久精品77| 青草久久国产| 国产精品爽爽va在线观看网站| 啦啦啦免费观看视频1| 一级a爱片免费观看的视频| 我的老师免费观看完整版| 悠悠久久av| 无遮挡黄片免费观看| 悠悠久久av| 啦啦啦免费观看视频1| 亚洲国产欧洲综合997久久,| 18禁美女被吸乳视频| 亚洲av成人av| 国产免费av片在线观看野外av| 国产精品久久久av美女十八| 麻豆久久精品国产亚洲av| 超碰成人久久| 中文字幕av在线有码专区| 黑人欧美特级aaaaaa片| 成人三级做爰电影| 亚洲成人久久爱视频| 桃红色精品国产亚洲av| 欧美黑人巨大hd| 91久久精品国产一区二区成人 | 99久久精品一区二区三区| 蜜桃久久精品国产亚洲av| 免费在线观看成人毛片| 香蕉国产在线看| 日本 欧美在线| 欧美激情在线99| 听说在线观看完整版免费高清| АⅤ资源中文在线天堂| 精品国产美女av久久久久小说| 少妇丰满av| 精品久久久久久,| 亚洲av成人不卡在线观看播放网| 我的老师免费观看完整版| 人人妻人人澡欧美一区二区| 欧美在线黄色| 又粗又爽又猛毛片免费看| 亚洲一区二区三区色噜噜| 一区二区三区激情视频| 国产亚洲av高清不卡| 三级国产精品欧美在线观看 | 精品久久久久久久毛片微露脸| 88av欧美| 国产午夜精品久久久久久| 国产主播在线观看一区二区| 男女午夜视频在线观看| 美女黄网站色视频| 97人妻精品一区二区三区麻豆| 男人舔女人的私密视频| 精品国产超薄肉色丝袜足j| 我的老师免费观看完整版| 国产精品乱码一区二三区的特点| 免费一级毛片在线播放高清视频| 狂野欧美激情性xxxx| 老司机午夜十八禁免费视频| 狠狠狠狠99中文字幕| 久久久国产精品麻豆| 亚洲18禁久久av| 97超视频在线观看视频| 在线永久观看黄色视频| 美女午夜性视频免费| 婷婷精品国产亚洲av在线| 女同久久另类99精品国产91| 日韩精品中文字幕看吧| 国产午夜精品论理片| 国产精品99久久久久久久久| 日韩国内少妇激情av| 一级黄色大片毛片| 欧美日韩黄片免| 在线视频色国产色| 91麻豆精品激情在线观看国产| 国产精品久久久久久精品电影| 欧美av亚洲av综合av国产av| 日本五十路高清| 久久人人精品亚洲av| 免费看a级黄色片| 最好的美女福利视频网| 此物有八面人人有两片| 男人舔女人的私密视频| 久久人人精品亚洲av| 国内揄拍国产精品人妻在线| 99久久精品热视频| 久久久成人免费电影| 国产免费av片在线观看野外av| 亚洲精品在线观看二区| 每晚都被弄得嗷嗷叫到高潮| 亚洲午夜理论影院| 久久久久久久久中文| 色综合婷婷激情| 国产伦一二天堂av在线观看| 免费看光身美女| 国产精品一区二区三区四区免费观看 | 国产精品久久久久久人妻精品电影| 亚洲精品中文字幕一二三四区| 亚洲国产欧美一区二区综合| 欧美黑人巨大hd| 欧美成人一区二区免费高清观看 | 日韩欧美三级三区| 亚洲18禁久久av| 床上黄色一级片| 国产一区二区三区在线臀色熟女| 嫩草影院精品99| 国产精品女同一区二区软件 | 亚洲成人精品中文字幕电影| 亚洲精品在线观看二区| 国产不卡一卡二| 欧美午夜高清在线| 97超视频在线观看视频| 成人特级av手机在线观看| 偷拍熟女少妇极品色| 国产精品 国内视频| 99热6这里只有精品| 中文字幕精品亚洲无线码一区| 亚洲人与动物交配视频| 国产单亲对白刺激| 亚洲欧美精品综合久久99| 丰满人妻熟妇乱又伦精品不卡| 午夜久久久久精精品| 国产av麻豆久久久久久久| 欧美乱色亚洲激情| 好看av亚洲va欧美ⅴa在| 日韩国内少妇激情av| 国产精品一区二区精品视频观看| 亚洲欧美精品综合一区二区三区| 久久国产精品影院| 999久久久精品免费观看国产| 成人av一区二区三区在线看| 国产三级中文精品| 国产精品久久久人人做人人爽| 欧美日韩中文字幕国产精品一区二区三区| 成年女人永久免费观看视频| 夜夜爽天天搞| 亚洲色图 男人天堂 中文字幕| 欧美黄色淫秽网站| 蜜桃久久精品国产亚洲av| 两个人看的免费小视频| 精品一区二区三区av网在线观看| 啦啦啦观看免费观看视频高清| 久久久久久久久中文| 国产精品99久久久久久久久| 免费无遮挡裸体视频| 黄色日韩在线| 精品国产亚洲在线| 成人国产一区最新在线观看| 久久国产乱子伦精品免费另类| 国产真人三级小视频在线观看| 黄色女人牲交| 亚洲,欧美精品.| or卡值多少钱| 日韩高清综合在线| 操出白浆在线播放| 精品久久蜜臀av无| 两个人看的免费小视频| 看片在线看免费视频| 超碰成人久久| 国产亚洲av高清不卡| 久久精品91蜜桃| 18禁观看日本| 啦啦啦免费观看视频1| 国产精品一及| 极品教师在线免费播放| 欧美xxxx黑人xx丫x性爽| 99精品久久久久人妻精品| 国内精品久久久久久久电影| 国产成人av教育| x7x7x7水蜜桃| 国内精品久久久久久久电影| 99热这里只有精品一区 | 亚洲精品乱码久久久v下载方式 | 国产97色在线日韩免费| 一级毛片女人18水好多| 老熟妇乱子伦视频在线观看| av中文乱码字幕在线| 天堂√8在线中文| 欧美一区二区国产精品久久精品| or卡值多少钱| 国产精品 国内视频| 亚洲成人久久爱视频| 欧美成人免费av一区二区三区| 蜜桃久久精品国产亚洲av| 可以在线观看毛片的网站| 国产伦在线观看视频一区| 亚洲色图av天堂| 日韩欧美 国产精品| 此物有八面人人有两片| 男人舔女人下体高潮全视频| 免费av毛片视频| 波多野结衣高清作品| 我要搜黄色片| 欧美乱妇无乱码| 99久久精品国产亚洲精品| 久久伊人香网站| 色综合婷婷激情| 成人精品一区二区免费| 国产精品国产高清国产av| 亚洲专区国产一区二区| 日韩欧美 国产精品| 免费电影在线观看免费观看| 制服丝袜大香蕉在线| 日韩有码中文字幕| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美国产在线观看| 日韩精品中文字幕看吧| 欧美又色又爽又黄视频| 色吧在线观看| 亚洲黑人精品在线| aaaaa片日本免费| 日本一本二区三区精品| 麻豆成人av在线观看| 欧美日本亚洲视频在线播放| 亚洲精品粉嫩美女一区| 国产精品久久视频播放| 18禁裸乳无遮挡免费网站照片| 免费av毛片视频| 国产高清视频在线观看网站| 免费一级毛片在线播放高清视频| 亚洲天堂国产精品一区在线| 国产精品乱码一区二三区的特点| 日韩中文字幕欧美一区二区| 午夜影院日韩av| 国产成年人精品一区二区| 一本久久中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 天堂av国产一区二区熟女人妻| 精品一区二区三区视频在线观看免费| 丰满人妻一区二区三区视频av | 中文字幕久久专区| 免费av不卡在线播放| 亚洲欧美日韩无卡精品| 天天添夜夜摸| 99久久精品国产亚洲精品| 亚洲国产日韩欧美精品在线观看 | 久久国产精品影院| 日本 av在线| 亚洲精品粉嫩美女一区| 麻豆av在线久日| 99久久综合精品五月天人人| 日日夜夜操网爽| 日韩中文字幕欧美一区二区| av国产免费在线观看| 亚洲中文av在线| 亚洲五月天丁香| 久久久国产成人精品二区| 亚洲午夜理论影院| 日本黄色视频三级网站网址| www日本黄色视频网| 女生性感内裤真人,穿戴方法视频| 一边摸一边抽搐一进一小说| 久久午夜综合久久蜜桃| 老司机深夜福利视频在线观看| 精品国产亚洲在线| 日本 av在线| 91久久精品国产一区二区成人 | 亚洲中文字幕日韩| 国产成人系列免费观看| 淫秽高清视频在线观看| 俺也久久电影网| 欧美中文日本在线观看视频| 国产精品野战在线观看| 女同久久另类99精品国产91| 国产真人三级小视频在线观看| 日本a在线网址| 一进一出抽搐动态| 怎么达到女性高潮| 在线观看免费视频日本深夜| 国产又色又爽无遮挡免费看| 狠狠狠狠99中文字幕| 国产精品美女特级片免费视频播放器 | 欧美乱妇无乱码| 美女扒开内裤让男人捅视频| 欧美不卡视频在线免费观看| 免费看a级黄色片| 亚洲熟妇中文字幕五十中出| 久久天躁狠狠躁夜夜2o2o| 夜夜夜夜夜久久久久| 亚洲av片天天在线观看| 女同久久另类99精品国产91| 国产成人系列免费观看| 91在线精品国自产拍蜜月 | 极品教师在线免费播放| 老司机午夜十八禁免费视频| 人人妻人人看人人澡| 日韩欧美国产一区二区入口| 国产午夜精品久久久久久| 国产成人aa在线观看| 国产伦精品一区二区三区四那| 伊人久久大香线蕉亚洲五| 久久久国产精品麻豆| 99久久国产精品久久久| 国产亚洲精品av在线| 婷婷精品国产亚洲av在线| 国产精品女同一区二区软件 | 欧美性猛交╳xxx乱大交人| www国产在线视频色| 人人妻人人看人人澡| 国产精品久久久久久亚洲av鲁大| 听说在线观看完整版免费高清| xxxwww97欧美| 国产亚洲精品av在线| 日韩欧美在线二视频| 757午夜福利合集在线观看| 国产欧美日韩一区二区三| 亚洲五月婷婷丁香| 最近最新中文字幕大全电影3| 精品国产超薄肉色丝袜足j| 美女高潮喷水抽搐中文字幕| www.www免费av| 高潮久久久久久久久久久不卡| 婷婷亚洲欧美| 黄色女人牲交| 免费观看的影片在线观看| 亚洲午夜理论影院| av天堂在线播放| 最近视频中文字幕2019在线8| 亚洲欧美激情综合另类| 免费看光身美女| 国内揄拍国产精品人妻在线| 欧美xxxx黑人xx丫x性爽| 精品久久久久久久久久久久久| 午夜两性在线视频| 中文字幕人妻丝袜一区二区| 日本一二三区视频观看| 国产1区2区3区精品| 91字幕亚洲| 听说在线观看完整版免费高清| 97超级碰碰碰精品色视频在线观看| 国产精品亚洲美女久久久| 少妇丰满av| av视频在线观看入口| 一二三四在线观看免费中文在| 97碰自拍视频| 嫩草影院精品99| 欧美极品一区二区三区四区| 国内精品美女久久久久久| 色综合站精品国产| 精品熟女少妇八av免费久了| 久久国产精品人妻蜜桃| 91字幕亚洲| 日韩欧美免费精品| 高清在线国产一区| 久久久久久人人人人人| 亚洲av熟女| 亚洲男人的天堂狠狠| 十八禁人妻一区二区| 香蕉国产在线看| 99久久99久久久精品蜜桃| 亚洲美女视频黄频| 国产麻豆成人av免费视频| 亚洲va日本ⅴa欧美va伊人久久| 欧美黄色片欧美黄色片| 日本黄色片子视频| 夜夜看夜夜爽夜夜摸| 久久精品影院6| 午夜激情福利司机影院| 国产精品99久久久久久久久| 午夜a级毛片| 成人无遮挡网站| 亚洲人成网站高清观看| 又粗又爽又猛毛片免费看| 黄色视频,在线免费观看| 亚洲第一欧美日韩一区二区三区| 老司机在亚洲福利影院| 国产精品久久久人人做人人爽| 好男人在线观看高清免费视频| 免费在线观看日本一区| 久久欧美精品欧美久久欧美| 我要搜黄色片| 成人三级做爰电影| 亚洲专区国产一区二区| 国产精品一区二区三区四区久久| 久久久久性生活片| 亚洲一区二区三区不卡视频| 草草在线视频免费看| 一本久久中文字幕| 国产欧美日韩精品一区二区| 女人高潮潮喷娇喘18禁视频| 国产伦人伦偷精品视频| 国产激情久久老熟女| 精品一区二区三区视频在线观看免费| 亚洲精品美女久久久久99蜜臀| 欧美日韩一级在线毛片| 中文资源天堂在线| 国产真人三级小视频在线观看| 日本 欧美在线| 精品国产乱子伦一区二区三区| 国产三级中文精品| 免费人成视频x8x8入口观看| 国产aⅴ精品一区二区三区波| 亚洲一区二区三区色噜噜| 天堂√8在线中文| 首页视频小说图片口味搜索| 黄色成人免费大全| 日本精品一区二区三区蜜桃| 国产69精品久久久久777片 | www.自偷自拍.com| 一夜夜www| 91老司机精品| 在线视频色国产色| 国模一区二区三区四区视频 | 日本 av在线| 国产日本99.免费观看| 亚洲人成网站在线播放欧美日韩| 久久久久久久午夜电影| 亚洲av电影不卡..在线观看| 日韩欧美 国产精品| 亚洲五月天丁香| 欧美日本视频| 岛国视频午夜一区免费看| av片东京热男人的天堂| 国产亚洲av高清不卡| 嫩草影院精品99| 99热这里只有是精品50| av欧美777| 免费在线观看亚洲国产| 在线a可以看的网站| 在线观看美女被高潮喷水网站 | 丁香六月欧美| 亚洲成人中文字幕在线播放| а√天堂www在线а√下载| 成人18禁在线播放| 免费看a级黄色片| 两个人视频免费观看高清| www.熟女人妻精品国产| 亚洲欧美日韩无卡精品| 69av精品久久久久久| 特级一级黄色大片| 亚洲精品粉嫩美女一区| 国产 一区 欧美 日韩| 草草在线视频免费看| 亚洲第一欧美日韩一区二区三区| 国产三级中文精品| 午夜福利在线观看吧| 亚洲国产欧美网| 久久天堂一区二区三区四区| 成人鲁丝片一二三区免费| 黄色 视频免费看| 欧美在线一区亚洲| 人人妻人人看人人澡| 日本精品一区二区三区蜜桃| 婷婷丁香在线五月| 亚洲专区中文字幕在线| 午夜免费激情av| 天天一区二区日本电影三级| 亚洲精品美女久久av网站| 精品99又大又爽又粗少妇毛片 | 国产精品一区二区免费欧美| av中文乱码字幕在线| 1024手机看黄色片| 白带黄色成豆腐渣| 免费观看人在逋| 免费看光身美女| 天天添夜夜摸| 日日摸夜夜添夜夜添小说| 久久香蕉国产精品| 成人欧美大片| 麻豆国产97在线/欧美| 丁香六月欧美| 久久精品综合一区二区三区| 一区二区三区国产精品乱码| 久久精品影院6| 亚洲欧美日韩高清在线视频| www.999成人在线观看| 人妻丰满熟妇av一区二区三区| 亚洲熟妇中文字幕五十中出| 日本一本二区三区精品| 1024香蕉在线观看| 国产人伦9x9x在线观看| 日韩 欧美 亚洲 中文字幕| 免费看日本二区| 搡老岳熟女国产| 亚洲成a人片在线一区二区| 又粗又爽又猛毛片免费看| 一进一出抽搐动态| 亚洲 国产 在线| 在线免费观看不下载黄p国产 | 国产精品爽爽va在线观看网站| 免费在线观看日本一区| 99在线人妻在线中文字幕| 免费一级毛片在线播放高清视频| 一二三四社区在线视频社区8| 久久热在线av| 亚洲精品色激情综合| 91老司机精品| 美女高潮喷水抽搐中文字幕| 一本精品99久久精品77| 1024香蕉在线观看| av黄色大香蕉| 美女黄网站色视频| 夜夜躁狠狠躁天天躁| 欧美黑人欧美精品刺激| 可以在线观看毛片的网站| 99国产综合亚洲精品| 国产伦在线观看视频一区| 日本免费a在线| 欧美日韩瑟瑟在线播放| 啪啪无遮挡十八禁网站| 一个人看的www免费观看视频| 男女视频在线观看网站免费| 亚洲av五月六月丁香网| 老鸭窝网址在线观看| 一级黄色大片毛片| 99riav亚洲国产免费| 午夜两性在线视频| 久久性视频一级片| 麻豆久久精品国产亚洲av| 男插女下体视频免费在线播放| 精品国产亚洲在线|