• <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
    国产1区2区3区精品| 欧美日韩黄片免| 精品久久久久久久久久久久久| 亚洲五月婷婷丁香| 日日干狠狠操夜夜爽| 男女床上黄色一级片免费看| 亚洲国产看品久久| 天天躁夜夜躁狠狠躁躁| 午夜日韩欧美国产| 亚洲在线自拍视频| 国产亚洲欧美98| 俄罗斯特黄特色一大片| 无遮挡黄片免费观看| 白带黄色成豆腐渣| 精品第一国产精品| 成人三级做爰电影| 欧美乱妇无乱码| 亚洲一码二码三码区别大吗| 国产av不卡久久| 巨乳人妻的诱惑在线观看| 成人av一区二区三区在线看| 无限看片的www在线观看| 久久久久久久久免费视频了| 99久久精品国产亚洲精品| 亚洲精品中文字幕一二三四区| 久久99热这里只有精品18| 国产精品影院久久| 国产精品1区2区在线观看.| 国产激情欧美一区二区| 最新美女视频免费是黄的| 国内精品久久久久久久电影| 国产精品久久视频播放| 听说在线观看完整版免费高清| 天堂√8在线中文| 美女黄网站色视频| 99re在线观看精品视频| 成人av一区二区三区在线看| 免费无遮挡裸体视频| 欧美日韩国产亚洲二区| 亚洲国产欧美网| 韩国av一区二区三区四区| 亚洲人成伊人成综合网2020| 久久久久久久精品吃奶| 亚洲熟妇中文字幕五十中出| 国产精品久久久人人做人人爽| 男人舔奶头视频| 亚洲精品在线美女| 亚洲一区二区三区不卡视频| 看黄色毛片网站| 一级毛片精品| 亚洲在线自拍视频| 国产精品一区二区免费欧美| 国产人伦9x9x在线观看| 欧美成人免费av一区二区三区| 18禁黄网站禁片午夜丰满| 精品欧美国产一区二区三| а√天堂www在线а√下载| 国产黄色小视频在线观看| 国产私拍福利视频在线观看| 午夜福利在线在线| 精品人妻1区二区| 91字幕亚洲| 日韩免费av在线播放| 99久久久亚洲精品蜜臀av| 婷婷亚洲欧美| 白带黄色成豆腐渣| av福利片在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲午夜理论影院| 中文字幕av在线有码专区| 国产av麻豆久久久久久久| 一级毛片高清免费大全| 免费在线观看日本一区| 久久久久久久久免费视频了| 50天的宝宝边吃奶边哭怎么回事| 国产高清videossex| 中文在线观看免费www的网站 | 国产精品一区二区三区四区免费观看 | 视频区欧美日本亚洲| 亚洲av成人av| 麻豆国产97在线/欧美 | 亚洲精华国产精华精| 99国产综合亚洲精品| 精品国产乱码久久久久久男人| 欧美日韩中文字幕国产精品一区二区三区| 亚洲熟女毛片儿| 日本熟妇午夜| 国产精品一区二区免费欧美| 精品一区二区三区av网在线观看| 欧美久久黑人一区二区| 法律面前人人平等表现在哪些方面| 国产单亲对白刺激| 特大巨黑吊av在线直播| 欧美日韩亚洲国产一区二区在线观看| 狂野欧美激情性xxxx| 欧洲精品卡2卡3卡4卡5卡区| 国产精品久久电影中文字幕| 成人三级做爰电影| 欧美中文综合在线视频| 亚洲avbb在线观看| www.精华液| www.精华液| 国产精品久久电影中文字幕| 国产精品永久免费网站| 成人18禁在线播放| 国产精品 欧美亚洲| 日韩欧美在线二视频| 中文在线观看免费www的网站 | 精华霜和精华液先用哪个| 久久欧美精品欧美久久欧美| 亚洲成人精品中文字幕电影| 亚洲美女黄片视频| 国产成人系列免费观看| 一a级毛片在线观看| 男人舔女人下体高潮全视频| 欧美在线黄色| 国产精品av视频在线免费观看| 深夜精品福利| 国产三级黄色录像| 黄色视频,在线免费观看| 中亚洲国语对白在线视频| 国产成人影院久久av| 色播亚洲综合网| 日日干狠狠操夜夜爽| 亚洲五月天丁香| 国产精华一区二区三区| 免费在线观看黄色视频的| 午夜激情福利司机影院| 久久久久性生活片| 亚洲国产精品久久男人天堂| 国产aⅴ精品一区二区三区波| 50天的宝宝边吃奶边哭怎么回事| 亚洲av成人一区二区三| 在线观看舔阴道视频| 搞女人的毛片| 亚洲欧美精品综合一区二区三区| 不卡一级毛片| 波多野结衣高清无吗| 亚洲美女黄片视频| 国产在线观看jvid| 国产私拍福利视频在线观看| 亚洲精品美女久久久久99蜜臀| 国产一区二区在线av高清观看| 变态另类丝袜制服| 久久久水蜜桃国产精品网| 久久久久久久午夜电影| 亚洲人成77777在线视频| 亚洲狠狠婷婷综合久久图片| 51午夜福利影视在线观看| 夜夜爽天天搞| 久久久久国内视频| 在线观看免费午夜福利视频| 国产69精品久久久久777片 | 亚洲一码二码三码区别大吗| 午夜免费激情av| 欧美色视频一区免费| 欧美3d第一页| 亚洲第一电影网av| 国产激情欧美一区二区| 又黄又爽又免费观看的视频| 又紧又爽又黄一区二区| 18禁观看日本| 五月伊人婷婷丁香| 一二三四社区在线视频社区8| 国内毛片毛片毛片毛片毛片| 亚洲一码二码三码区别大吗| 俺也久久电影网| 亚洲欧美日韩高清专用| 99精品久久久久人妻精品| 别揉我奶头~嗯~啊~动态视频| 在线观看日韩欧美| 少妇的丰满在线观看| 观看免费一级毛片| 亚洲中文日韩欧美视频| 757午夜福利合集在线观看| 不卡一级毛片| 久久精品夜夜夜夜夜久久蜜豆 | 18美女黄网站色大片免费观看| bbb黄色大片| 蜜桃久久精品国产亚洲av| 中文字幕高清在线视频| 岛国在线观看网站| 精品国产亚洲在线| 啦啦啦免费观看视频1| 五月玫瑰六月丁香| 香蕉丝袜av| 在线观看午夜福利视频| 亚洲成av人片免费观看| 免费在线观看成人毛片| 亚洲avbb在线观看| 日日夜夜操网爽| 成年版毛片免费区| 波多野结衣高清无吗| 丰满人妻一区二区三区视频av | 这个男人来自地球电影免费观看| 在线观看免费午夜福利视频| 成人三级做爰电影| 国产av麻豆久久久久久久| 欧美精品亚洲一区二区| 亚洲美女黄片视频| 国产成年人精品一区二区| 精品国产乱码久久久久久男人| 午夜免费激情av| 成人国语在线视频| 脱女人内裤的视频| 老熟妇乱子伦视频在线观看| 婷婷精品国产亚洲av在线| 国产91精品成人一区二区三区| 亚洲狠狠婷婷综合久久图片| www.www免费av| 精品无人区乱码1区二区| 国产片内射在线| 午夜日韩欧美国产| 中文亚洲av片在线观看爽| 又粗又爽又猛毛片免费看| av在线播放免费不卡| 一二三四社区在线视频社区8| 岛国在线免费视频观看| 成年人黄色毛片网站| 两人在一起打扑克的视频| 国产片内射在线| 熟妇人妻久久中文字幕3abv| 婷婷精品国产亚洲av| 日韩欧美国产一区二区入口| av在线播放免费不卡| 免费高清视频大片| 国产成人精品无人区| 国产成人影院久久av| 老司机福利观看| 一级毛片精品| 激情在线观看视频在线高清| 老司机靠b影院| 99久久精品热视频| 18禁美女被吸乳视频| 高清毛片免费观看视频网站| 成人三级做爰电影| 麻豆一二三区av精品| 亚洲成人免费电影在线观看| 日韩 欧美 亚洲 中文字幕| 99国产综合亚洲精品| 老司机午夜十八禁免费视频| 国产精品永久免费网站| 91九色精品人成在线观看| 不卡一级毛片| av在线播放免费不卡| 嫩草影院精品99| 黄片大片在线免费观看| 成人三级做爰电影| 黑人操中国人逼视频| 日韩中文字幕欧美一区二区| 五月伊人婷婷丁香| 亚洲狠狠婷婷综合久久图片| 香蕉国产在线看| 欧洲精品卡2卡3卡4卡5卡区| 精品高清国产在线一区| 动漫黄色视频在线观看| 美女免费视频网站| 90打野战视频偷拍视频| 国产真实乱freesex| 国产激情久久老熟女| 久久久久久久久中文| 12—13女人毛片做爰片一| 在线永久观看黄色视频| 久久亚洲真实| 99久久精品国产亚洲精品| 一本精品99久久精品77| 国产亚洲精品一区二区www| 国产一区二区三区视频了| 国产免费av片在线观看野外av| 日本免费a在线| 狂野欧美白嫩少妇大欣赏| 亚洲人成网站高清观看| 一二三四在线观看免费中文在| 一级a爱片免费观看的视频| 色综合站精品国产| 午夜免费激情av| 精品久久久久久成人av| 中文资源天堂在线| 日本免费一区二区三区高清不卡| 国产av不卡久久| 两性夫妻黄色片| 国产aⅴ精品一区二区三区波| 老熟妇乱子伦视频在线观看| 精华霜和精华液先用哪个| 91av网站免费观看| 99久久精品国产亚洲精品| 国产成人精品久久二区二区免费| 久久婷婷成人综合色麻豆| 国产91精品成人一区二区三区| 免费人成视频x8x8入口观看| www国产在线视频色| 欧美日本视频| 国产男靠女视频免费网站| 丝袜美腿诱惑在线| 午夜老司机福利片| 免费观看精品视频网站| 色哟哟哟哟哟哟| 99精品欧美一区二区三区四区| 桃红色精品国产亚洲av| 久久人妻av系列| 欧美成人午夜精品| 免费电影在线观看免费观看| 黄色视频,在线免费观看| 18禁美女被吸乳视频| 精品国产美女av久久久久小说| 国内久久婷婷六月综合欲色啪| 国产精品久久久久久亚洲av鲁大| 国产成人av教育| 日日干狠狠操夜夜爽| 三级毛片av免费| 亚洲av中文字字幕乱码综合| 色哟哟哟哟哟哟| www.999成人在线观看| tocl精华| 一级毛片女人18水好多| 国产高清视频在线播放一区| 琪琪午夜伦伦电影理论片6080| 久久天躁狠狠躁夜夜2o2o| 国产又黄又爽又无遮挡在线| 亚洲欧美精品综合一区二区三区| 男女做爰动态图高潮gif福利片| 欧美日韩亚洲综合一区二区三区_| 18禁黄网站禁片午夜丰满| 日日夜夜操网爽| 大型黄色视频在线免费观看| 女生性感内裤真人,穿戴方法视频| 午夜精品久久久久久毛片777| 精品一区二区三区av网在线观看| 麻豆成人av在线观看| 国产精品一及| 日本一二三区视频观看| 久久国产精品人妻蜜桃| 久久国产精品影院| 国产精品,欧美在线| 亚洲五月婷婷丁香| 国产精品电影一区二区三区| 国产爱豆传媒在线观看 | 老司机福利观看| 亚洲av五月六月丁香网| 婷婷六月久久综合丁香| 啪啪无遮挡十八禁网站| 欧美日韩亚洲国产一区二区在线观看| 欧美精品亚洲一区二区| 91麻豆精品激情在线观看国产| 丝袜人妻中文字幕| 久久热在线av| 久久香蕉精品热| www.999成人在线观看| 成人av在线播放网站| 国产1区2区3区精品| 99国产极品粉嫩在线观看| 欧美成人性av电影在线观看| 欧美中文综合在线视频| 国产精品免费视频内射| 国产精品一及| 午夜精品在线福利| 国内精品久久久久精免费| 18禁观看日本| 国产精品自产拍在线观看55亚洲| 国产区一区二久久| 国内精品一区二区在线观看| 在线观看一区二区三区| 又大又爽又粗| 亚洲av成人精品一区久久| 亚洲乱码一区二区免费版| 国产精品av久久久久免费| 99精品久久久久人妻精品| 一区二区三区国产精品乱码| 无遮挡黄片免费观看| 日本 欧美在线| 久久久精品大字幕| 亚洲全国av大片| 成在线人永久免费视频| 欧美成人免费av一区二区三区| 一级黄色大片毛片| 亚洲熟妇中文字幕五十中出| 欧美色欧美亚洲另类二区| 久久久久国产精品人妻aⅴ院| 香蕉av资源在线| 毛片女人毛片| 天堂av国产一区二区熟女人妻 | 亚洲精品粉嫩美女一区| 久久久久久亚洲精品国产蜜桃av| 欧美精品亚洲一区二区| 99热只有精品国产| 欧美色视频一区免费| 日韩欧美在线乱码| 中文字幕人成人乱码亚洲影| 在线观看免费日韩欧美大片| 欧美最黄视频在线播放免费| 久久国产乱子伦精品免费另类| 搡老妇女老女人老熟妇| 99精品在免费线老司机午夜| 国产午夜福利久久久久久| 19禁男女啪啪无遮挡网站| 国产成人欧美在线观看| 又紧又爽又黄一区二区| 九色成人免费人妻av| 精品无人区乱码1区二区| 欧美日韩精品网址| 亚洲av成人av| 中出人妻视频一区二区| 一区二区三区激情视频| 久久中文字幕一级| 日本在线视频免费播放| 久久久久久人人人人人| av视频在线观看入口| 国产av一区在线观看免费| 精品国产亚洲在线| videosex国产| 黄色丝袜av网址大全| 国产精品久久久久久亚洲av鲁大| 两个人视频免费观看高清| 精品一区二区三区av网在线观看| 欧美又色又爽又黄视频| 一个人观看的视频www高清免费观看 | 国产高清视频在线播放一区| 欧美日韩亚洲综合一区二区三区_| 精品久久久久久久久久久久久| 两个人的视频大全免费| 久久久久久久午夜电影| 亚洲精品国产一区二区精华液| 一进一出抽搐gif免费好疼| 国产一区二区三区在线臀色熟女| 日韩欧美在线二视频| 国产高清视频在线播放一区| 两个人视频免费观看高清| 精品一区二区三区av网在线观看| 亚洲自拍偷在线| 久久午夜亚洲精品久久| 欧美日韩国产亚洲二区| 亚洲精品色激情综合| 可以免费在线观看a视频的电影网站| 在线观看免费视频日本深夜| 国产一区二区在线av高清观看| 99久久综合精品五月天人人| 每晚都被弄得嗷嗷叫到高潮| 午夜亚洲福利在线播放| 最新美女视频免费是黄的| 欧美在线黄色| 欧美人与性动交α欧美精品济南到| 日韩国内少妇激情av| 国产午夜精品久久久久久| √禁漫天堂资源中文www| 亚洲专区中文字幕在线| 亚洲五月天丁香| 在线国产一区二区在线| 欧美黑人精品巨大| 中文字幕精品亚洲无线码一区| av免费在线观看网站| 亚洲九九香蕉| 黄色a级毛片大全视频| 日韩成人在线观看一区二区三区| 一本久久中文字幕| 超碰成人久久| 亚洲精品粉嫩美女一区| 日本熟妇午夜| 亚洲中文字幕日韩| 18禁美女被吸乳视频| 欧美久久黑人一区二区| 欧美日本亚洲视频在线播放| 久久久国产成人免费| 又黄又爽又免费观看的视频| 国产一区二区在线av高清观看| 亚洲欧洲精品一区二区精品久久久| 美女免费视频网站| 午夜老司机福利片| 日韩中文字幕欧美一区二区| 国产亚洲精品久久久久5区| 夜夜躁狠狠躁天天躁| 99久久精品热视频| aaaaa片日本免费| 精品不卡国产一区二区三区| 最近视频中文字幕2019在线8| 两个人看的免费小视频| 国产av在哪里看| 国产精华一区二区三区| 国产精品国产高清国产av| 国产亚洲欧美在线一区二区| 黑人欧美特级aaaaaa片| 国产又色又爽无遮挡免费看| 在线看三级毛片| 12—13女人毛片做爰片一| av超薄肉色丝袜交足视频| 国产1区2区3区精品| 精品乱码久久久久久99久播| 精品久久蜜臀av无| 黄色 视频免费看| 亚洲五月婷婷丁香| 国产精品99久久99久久久不卡| 在线观看美女被高潮喷水网站 | 可以免费在线观看a视频的电影网站| 91字幕亚洲| 男人的好看免费观看在线视频 | 久久亚洲真实| 午夜日韩欧美国产| 久久精品国产清高在天天线| 久久精品国产99精品国产亚洲性色| 国产爱豆传媒在线观看 | 777久久人妻少妇嫩草av网站| 午夜福利免费观看在线| 亚洲在线自拍视频| 亚洲美女视频黄频| 亚洲第一欧美日韩一区二区三区| 久久欧美精品欧美久久欧美| 国产av又大| 黄色成人免费大全| 久久精品国产亚洲av高清一级| 亚洲av中文字字幕乱码综合| 宅男免费午夜| 久久精品国产99精品国产亚洲性色| 亚洲精品中文字幕在线视频| 亚洲五月天丁香| 亚洲av五月六月丁香网| 黄频高清免费视频| 亚洲av成人一区二区三| 1024香蕉在线观看| 亚洲avbb在线观看| 搡老岳熟女国产| 中文字幕精品亚洲无线码一区| 欧美日韩黄片免| 18美女黄网站色大片免费观看| 又黄又粗又硬又大视频| 老汉色av国产亚洲站长工具| 日韩欧美在线二视频| 国产野战对白在线观看| 国产真实乱freesex| 全区人妻精品视频| 少妇粗大呻吟视频| 这个男人来自地球电影免费观看| www日本在线高清视频| 久久久精品欧美日韩精品| 成人欧美大片| 男女床上黄色一级片免费看| www.精华液| 国产不卡一卡二| 又爽又黄无遮挡网站| 人人妻,人人澡人人爽秒播| 后天国语完整版免费观看| 亚洲电影在线观看av| 久久天躁狠狠躁夜夜2o2o| 白带黄色成豆腐渣| 国产成人精品久久二区二区91| 又黄又粗又硬又大视频| 亚洲av电影在线进入| 99国产精品99久久久久| 天堂动漫精品| 国内精品久久久久精免费| 欧美黄色片欧美黄色片| 国内久久婷婷六月综合欲色啪| 一区福利在线观看| 亚洲精品国产一区二区精华液| 亚洲国产精品999在线| 久久精品国产99精品国产亚洲性色| 18禁观看日本| 黄色片一级片一级黄色片| 大型黄色视频在线免费观看| 国产成人aa在线观看| 一进一出抽搐gif免费好疼| 成人av在线播放网站| 国产精品综合久久久久久久免费| 麻豆国产av国片精品| 老熟妇仑乱视频hdxx| 午夜精品久久久久久毛片777| 老司机靠b影院| 久久久国产成人免费| 午夜视频精品福利| 别揉我奶头~嗯~啊~动态视频| 12—13女人毛片做爰片一| 女警被强在线播放| 好男人电影高清在线观看| 一级作爱视频免费观看| 12—13女人毛片做爰片一| 成人av在线播放网站| 欧美+亚洲+日韩+国产| 久久久精品国产亚洲av高清涩受| 男人舔奶头视频| 老司机午夜十八禁免费视频| 国产av在哪里看| 久久久国产成人精品二区| 搡老岳熟女国产| 精品久久蜜臀av无| 在线观看一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 中文字幕人成人乱码亚洲影| 一卡2卡三卡四卡精品乱码亚洲| 国产主播在线观看一区二区| 亚洲第一欧美日韩一区二区三区| 日韩 欧美 亚洲 中文字幕| 99久久综合精品五月天人人| 日韩 欧美 亚洲 中文字幕| 国产精品久久电影中文字幕| 少妇被粗大的猛进出69影院| 午夜久久久久精精品| 中文字幕人妻丝袜一区二区| xxx96com| 国产三级黄色录像| 亚洲成av人片在线播放无| 搡老妇女老女人老熟妇| 国产精品国产高清国产av| 亚洲精华国产精华精| 欧美三级亚洲精品| 亚洲一卡2卡3卡4卡5卡精品中文| 大型黄色视频在线免费观看| 午夜老司机福利片| 日韩成人在线观看一区二区三区| 精品久久久久久成人av| 香蕉av资源在线| 窝窝影院91人妻| videosex国产| 欧美乱码精品一区二区三区| 精品免费久久久久久久清纯| 老司机在亚洲福利影院| 丝袜人妻中文字幕| 日韩欧美 国产精品|