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

    新型文丘里空化器及其空化特性的數(shù)值模擬

    2023-10-25 11:41:50陳偉政
    船舶力學(xué) 2023年10期
    關(guān)鍵詞:文丘里喉部空泡

    海 青,陳偉政,顏 開

    (1.中國(guó)船舶科學(xué)研究中心,江蘇無(wú)錫 214082;2.深海技術(shù)科學(xué)太湖實(shí)驗(yàn)室,江蘇無(wú)錫 214082)

    0 引 言

    隨著人類社會(huì)的快速發(fā)展,工農(nóng)業(yè)生產(chǎn)活動(dòng)和日常生活產(chǎn)生的廢水逐漸呈現(xiàn)出成份多樣、危害較大的特點(diǎn),傳統(tǒng)的水處理方法已難以有效降解日趨復(fù)雜的廢水??栈夹g(shù)作為一種高級(jí)氧化工藝(advance oxidation process,AOP)可有效降解廢水中的芳香烴、雜環(huán)類、酚類等有機(jī)污染物[1],同時(shí)也可以應(yīng)用于消毒[2]、處理藍(lán)藻等微生物和過(guò)程強(qiáng)化[3]等領(lǐng)域??栈F(xiàn)象降解廢水的機(jī)理可以簡(jiǎn)要概括為:利用空泡潰滅時(shí)產(chǎn)生的局部極端物理?xiàng)l件與因此發(fā)生的化學(xué)反應(yīng)[4]來(lái)達(dá)到凈水目的。在空泡潰滅瞬間,會(huì)在爆點(diǎn)周圍的極小范圍內(nèi)產(chǎn)生高溫(103~104K)、高壓(102~103MPa)和高速微射流,在高溫高壓環(huán)境下,水分子分解產(chǎn)生具有強(qiáng)氧化性的羥基自由基(·OH),·OH 與污染物發(fā)生反應(yīng)達(dá)到降解的效果;同時(shí),高速微射流的強(qiáng)剪切效應(yīng)可以打斷大分子物質(zhì)的主鏈化學(xué)鍵,也可以破壞微生物細(xì)胞壁,從而達(dá)到降解和使微生物失活的效果[5]。

    水力空化技術(shù)因其設(shè)備結(jié)構(gòu)簡(jiǎn)單、能耗較低、產(chǎn)生空泡數(shù)目多以及易于工業(yè)放大的特點(diǎn)而更加具有應(yīng)用價(jià)值,相較于超聲空化等方法已被視為最具前景的空化降解方式[6]。水力空化發(fā)生器的類型和結(jié)構(gòu)是空化降解系統(tǒng)的關(guān)鍵,其流場(chǎng)特性直接影響降解效率。水力空化發(fā)生器類型一般通過(guò)幾何結(jié)構(gòu)變化使液體流速增加或發(fā)生流動(dòng)分離,在有限區(qū)域內(nèi)形成低壓,進(jìn)而發(fā)生空化。常見的水力空化發(fā)生器有孔板、文丘里管、漩渦型和對(duì)沖射流型[7]。其中,文丘里管作為一種典型且有效的水力空化發(fā)生器受到了廣泛關(guān)注。

    目前,針對(duì)文丘里管空化的數(shù)值模擬方法已經(jīng)相對(duì)成熟,基于穩(wěn)態(tài)RANS 方法可以對(duì)空化流場(chǎng)進(jìn)行計(jì)算,并以此來(lái)評(píng)估文丘里管內(nèi)的空化區(qū)域、壓力場(chǎng)和速度場(chǎng)等分布情況,為水力空化設(shè)備提供了一種有效的分析和優(yōu)化方法。在幾何結(jié)構(gòu)上,文丘里管中喉部的周長(zhǎng)與流動(dòng)面積的比值α反映了剪切區(qū)域的大小。Bashir等[8]對(duì)5種不同喉部截面的文丘里管(α=2~10)進(jìn)行了數(shù)值模擬,如圖1所示,研究發(fā)現(xiàn)在流動(dòng)面積相同且入口壓力一致的條件下,圓形和矩形截面的文丘里管均在環(huán)形狀態(tài)下獲得了更大的空化區(qū)域,雙矩環(huán)形截面(圖1(e))的文丘里管的空化區(qū)域最大;Abbasi 等[9]對(duì)多組文丘里管空化流場(chǎng)進(jìn)行了CFD計(jì)算,并利用響應(yīng)面方法(RSM)進(jìn)行優(yōu)化設(shè)計(jì),以此來(lái)確定最優(yōu)的幾何結(jié)構(gòu)與運(yùn)行參數(shù)組合;Jain等[10]研究了喉部截面分別為圓形和窄縫的兩種文丘里管空化面積的大小隨著運(yùn)行及幾何參數(shù)的變化規(guī)律,包括不同的入口壓力、擴(kuò)張角度和長(zhǎng)徑比。數(shù)值模擬結(jié)果顯示:對(duì)于圓截面文丘里管,當(dāng)入口壓力為0.8 MPa時(shí)空化面積最大,而窄縫截面文丘里管可以在入口壓力為0.5MPa時(shí)達(dá)到最大值,可見窄縫截面的文丘里管在較低入口壓力時(shí)就可以獲得最大空化效率,節(jié)省能源;Dutta等[11]針對(duì)文丘里管內(nèi)壁粗糙度和擴(kuò)張段線形進(jìn)行了數(shù)值模擬研究,研究發(fā)現(xiàn),增大內(nèi)壁粗糙度和將擴(kuò)張段線形由直線改為曲線(類喇叭狀)也可以達(dá)到延長(zhǎng)下游低壓區(qū)的效果,以此增大空化區(qū)域;王澤鵬等[12]利用CFD 方法分析了湍流作用下入口壓力及不同喉部直徑對(duì)空泡運(yùn)動(dòng)特性及其所形成壓力脈動(dòng)的影響規(guī)律,發(fā)現(xiàn)在潰滅區(qū)流體湍動(dòng)能達(dá)到峰值,說(shuō)明湍動(dòng)能大小可以在一定程度上反映潰滅強(qiáng)弱;董波等[13]探究了文丘里管的喉徑與管徑比、入口角度、喉徑與喉長(zhǎng)比以及出口角度等結(jié)構(gòu)參數(shù)對(duì)空化效果的影響規(guī)律;宋永興等[14]通過(guò)數(shù)值方法研究了不同入口錐角和不同喉部直徑的文丘里管的空化特性,發(fā)現(xiàn)喉部直徑變化對(duì)下游含氣率影響較大,當(dāng)喉部直徑不變時(shí),入口收縮角的大小對(duì)初生空化數(shù)的影響較大。

    圖1 文獻(xiàn)[8]中文丘里管不同喉部結(jié)構(gòu)示意圖Fig.1 Schematic of different throat structures of Venturi tube in Ref.[8]

    綜合已報(bào)道的研究結(jié)果,發(fā)現(xiàn)大多數(shù)用于降解有機(jī)廢水的文丘里型空化器均以傳統(tǒng)文丘里管為主,主要的優(yōu)化內(nèi)容包括喉部長(zhǎng)徑比、收縮、擴(kuò)張角度大小和擴(kuò)張段線形等幾何參數(shù),空化器的組成結(jié)構(gòu)改變較少。同時(shí),設(shè)計(jì)與優(yōu)化文丘里管應(yīng)根據(jù)空化降解廢水的機(jī)理展開,即提高·OH產(chǎn)量和增加有機(jī)物與·OH 的摻混程度來(lái)提升降解效率。空化區(qū)域的大小反映了空泡數(shù)目的多少,所以空化體積大小可以作為評(píng)價(jià)文丘里管空化性能的重要指標(biāo),即空化體積越大,表示產(chǎn)生·OH的潛在能力越強(qiáng)。基于上述認(rèn)識(shí),本文提出一種新型文丘里管,為了探究其空化流場(chǎng)特性,采用Fluent 軟件對(duì)不同幾何參數(shù)的新型文丘里管在不同入口壓力下進(jìn)行數(shù)值模擬和分析。

    1 新型文丘里管的結(jié)構(gòu)設(shè)計(jì)

    一般的文丘里管由收縮段、喉部和擴(kuò)張段構(gòu)成。收縮角度對(duì)于空化產(chǎn)生的影響相對(duì)較小,研究發(fā)現(xiàn)收縮角度在22.6°~90°之間變化時(shí),空化體積分?jǐn)?shù)僅變化2%左右[14]。業(yè)內(nèi)多數(shù)學(xué)者將收縮段角度選定為40°~50°之間;擴(kuò)張角度的大小影響喉部下游壓力分布與空泡發(fā)展,從降解試驗(yàn)結(jié)果來(lái)看,文丘里型空化器的最優(yōu)擴(kuò)張角度一般為10°~14°;喉部的長(zhǎng)徑比影響空化發(fā)展,前人研究表明長(zhǎng)徑比對(duì)于文丘里型空化器的影響不及孔板型顯著,一般最優(yōu)值為1:1[7]。綜上,本文的優(yōu)化出發(fā)點(diǎn)是收縮角為45°、擴(kuò)張段角度β為8°~14°的矩形截面文丘里管。喉道長(zhǎng)、寬、高分別為12 mm、10 mm、10 mm;來(lái)流主管道邊長(zhǎng)為50 mm。

    本研究提出的新型文丘里管在其喉部固定一個(gè)分流部件,如圖2(a)所示。喉部分流部件的縱剖面為梯形結(jié)構(gòu),位于喉部的中心位置,如圖2(b)~(d)所示,其兩端與文丘里管喉部的內(nèi)壁固連,將喉部分隔為上下兩個(gè)通道,其總長(zhǎng)與喉道長(zhǎng)度Lth一致,為12 mm。一方面,喉部加裝的部件具有增加剪切周長(zhǎng)的作用,可以增加發(fā)生空化現(xiàn)象區(qū)域,通過(guò)增加空化汽泡數(shù)的方法增大·OH產(chǎn)量;另一方面,梯形截面可以相對(duì)減少流動(dòng)損失,增加被處理溶液在單位時(shí)間內(nèi)通過(guò)空化發(fā)生器的次數(shù)。為保證喉部流速具有可比性,需要流通面積一致,故積梯形部件的厚度H應(yīng)為固定值。考慮到過(guò)小的喉部流通面積會(huì)導(dǎo)致下游壓力恢復(fù)過(guò)快,則選定H=4 mm,此時(shí)喉部流動(dòng)面積與來(lái)流管道面積比為0.024。分流部件的特征長(zhǎng)度L影響喉部的壓力分布,為使空化有足夠的區(qū)域發(fā)展,且不超出喉部總長(zhǎng),在本研究中L介于4 mm到10 mm之間,喉部控制參數(shù)γ的范圍介于0.33~0.83。定義為

    圖2 新型文丘里管幾何結(jié)構(gòu)Fig.2 Geometric structure of novel Venturi tube

    2 數(shù)值方法

    2.1 多相流與空化模型

    空化流場(chǎng)是由液態(tài)水和汽態(tài)水蒸氣組成的多相流場(chǎng),流動(dòng)過(guò)程中存在空化氣泡的生產(chǎn)、發(fā)展與潰滅。在數(shù)值模擬中需要多相流模型和描述相間傳質(zhì)的空化模型。本研究選擇均相流模型,汽液兩相混合物的連續(xù)性方程和動(dòng)量方程為

    式中,下標(biāo)m 代表混合相,l 與v 分別代表液相與汽相,下標(biāo)i與j表示坐標(biāo)方向,ρ、μ、α、p和u分別表示流體密度、粘性系數(shù)、體積分?jǐn)?shù)、壓力和速度。

    描述質(zhì)量輸運(yùn)的空化模型選擇基于R-P方程發(fā)展而來(lái)的Schnerr-Sauer模型:

    上式中,式(5)表示汽化過(guò)程,式(6)表示液化過(guò)程。Pv表示液體的飽和蒸汽壓,它隨溶液溫度變化;RB表示汽核半徑,由式(7)確定。

    式中,n為汽核密度,一般取n=1×1013。

    2.2 湍流模型

    目前計(jì)算空化流場(chǎng)常用的湍流模型為大渦模擬(LES)和雷諾平均模型(RANS)。為了提高計(jì)算速度,滿足多種工況的計(jì)算要求,本研究采用RANS 方法在穩(wěn)態(tài)條件下進(jìn)行數(shù)值計(jì)算。根據(jù)文丘里管空化流場(chǎng)剪切強(qiáng)的特點(diǎn),本文選擇Realizablek-ε模型,該模型對(duì)于模擬射流、存在流動(dòng)分離、邊界層存在強(qiáng)逆壓梯度的流動(dòng)具有良好表現(xiàn)。湍動(dòng)能k和湍流耗散率ε的輸運(yùn)方程為

    式中,ρ是液體密度,μ是液體動(dòng)力粘度,μt是渦粘系數(shù),S是平均應(yīng)變率張量的模。

    在式(9)中,系數(shù)C1受下式控制:

    上式中,模型常數(shù)C2=1.9,σk=1,σε=1.2。

    2.3 求解器設(shè)置

    求解器選擇Ansys-Fluent,根據(jù)上文選擇的多相流和湍流模型進(jìn)行相應(yīng)設(shè)置,在空化模型中將蒸汽壓設(shè)置為20℃時(shí)水的飽和蒸汽壓值,即Pv=2338 Pa。在湍流模型設(shè)置中開啟增強(qiáng)壁面函數(shù)(Enhanced Wall Treatment);在求解算法中選擇的壓力-速度耦合求解算法為Coupled,由于空化現(xiàn)象是一個(gè)非定常過(guò)程,為了避免出現(xiàn)出口回流現(xiàn)象并提高計(jì)算精度,選擇開啟偽瞬態(tài)算法(Pseudo Transient);壓力離散方法選擇PRESTO,體積分?jǐn)?shù)離散方法選擇QUICK,其他物理量選擇二階迎風(fēng)格式。邊界條件設(shè)置:入口、出口選擇壓力入口和壓力出口,入口壓力范圍在0.2~0.6 MPa,出口壓力為0.1 MPa。壁面選擇無(wú)滑移邊界條件,對(duì)稱面設(shè)置為對(duì)稱邊界。

    2.4 網(wǎng)格劃分

    本文所設(shè)計(jì)的新型文丘里管因其不具備軸對(duì)稱特點(diǎn),需要進(jìn)行三維模擬,流動(dòng)域具有兩個(gè)對(duì)稱平面,則計(jì)算域選取1/4 流動(dòng)域,可節(jié)約計(jì)算資源、提高計(jì)算效率。文丘里管結(jié)構(gòu)相對(duì)規(guī)則,選用結(jié)構(gòu)化網(wǎng)格,在壁面附近和下游空化區(qū)域進(jìn)行加密,如圖3所示。

    圖3 計(jì)算域網(wǎng)格劃分Fig.3 Computational domain mesh

    為保證計(jì)算精度要求,選擇了最大空化體積、喉部速度最大值和喉部平均速度進(jìn)行網(wǎng)格無(wú)關(guān)性檢驗(yàn),在0.4 MPa 的入口壓力下進(jìn)行計(jì)算。結(jié)果如圖4 所示,可見網(wǎng)格總數(shù)在80萬(wàn)左右滿足計(jì)算要求。

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

    3.1 新型文丘里管空化流場(chǎng)特點(diǎn)

    文丘里管對(duì)稱面上的空化區(qū)域示意圖與空化體積分?jǐn)?shù)云圖如圖5~6所示,可見空化起始于喉部分流體平臺(tái)區(qū)前緣處和擴(kuò)張段與喉部的交界處,空化區(qū)域也因此分為兩個(gè)區(qū)域,即文丘里管擴(kuò)展段壁面附近的空化區(qū)域①,以及喉部分流體上方與下游的空化區(qū)域②。這正是新型文丘里管可以提升空化產(chǎn)量的原因。圖6 中,觀察距離入口160~250 mm 的截面體積分?jǐn)?shù)云圖,可以發(fā)現(xiàn)在喉部射流沖出的核心區(qū)域不產(chǎn)生空化(紅色)。在高速射流的周圍形成環(huán)形空化區(qū)域,這是由于該區(qū)域的流速快、慣性強(qiáng),喉部的分流體后方形成類似于后臺(tái)階流動(dòng)的低壓區(qū),綜合兩方面因素,空化區(qū)域被限制在分流體正后方;之后隨著流速變慢,汽相在下游融合,上下分離不明顯。

    圖4 網(wǎng)格無(wú)關(guān)性驗(yàn)證Fig.4 Mesh independence verification

    圖6 汽相體積分?jǐn)?shù)云圖Fig.6 Volume fraction contour of vapor phase

    從能量的觀點(diǎn)出發(fā),在文丘里水力空化發(fā)生器中,液體通過(guò)喉部后的永久壓降是空化能Eca的來(lái)源??张葑缘蛪簠^(qū)產(chǎn)生后,隨流動(dòng)在壓力恢復(fù)區(qū)潰滅。潰滅過(guò)程瞬間發(fā)生,伴隨著極端物理現(xiàn)象,主要包含:機(jī)械能Em(噪聲、高速微射流、激波等)、高溫產(chǎn)生的熱能Eth、電磁能Eel、產(chǎn)生化學(xué)反應(yīng)的化學(xué)能Ech。可近似認(rèn)為潰滅過(guò)程中空化能量的轉(zhuǎn)化過(guò)程為

    在特定的空化條件下,流場(chǎng)中單個(gè)空泡的潰滅能量是有限且大致固定的,因此,空泡數(shù)量就顯得更為重要。從圖6 可以觀察到,新型文丘里管的空化區(qū)域②增加了下游的空化體積,一方面,空泡數(shù)目因此增加,說(shuō)明有更多的空化能Eca參與潰滅瞬間的能量轉(zhuǎn)化,另一方面,增強(qiáng)了潰滅產(chǎn)生的羥基自由基在液體中的擴(kuò)散效果,有利于提高羥基自由基與污染物的接觸概率,提升降解效率。雖然潰滅產(chǎn)生的各種形式的能量難以分別量化計(jì)算,但高溫、微射流、化學(xué)反應(yīng)都參與了有機(jī)物降解和微生物滅活的過(guò)程,所以,總的空化能Eca越大,文丘里管的空化能力也就越強(qiáng)。

    3.2 不同入口壓力的影響

    入口壓力是水力空化器在使用過(guò)程中的重要運(yùn)行參數(shù)之一,空泡的發(fā)展與潰滅行為對(duì)入口壓力變化較為敏感??张莸臐鐝?qiáng)度主要取決于喉道下游部分湍流脈動(dòng)強(qiáng)度和出口恢復(fù)壓力的最終大小。入口壓力的增加會(huì)導(dǎo)致液體流速和系統(tǒng)中能量耗散率的增加,從而導(dǎo)致通過(guò)喉道的永久壓降更大,這會(huì)使系統(tǒng)向單位質(zhì)量液體輸入的功率增加[15]。同時(shí),入口壓力的大小會(huì)改變空化數(shù)大小??栈瘮?shù)Cv是描述空化發(fā)生劇烈程度的無(wú)量綱參數(shù),定義為

    式中,Pout是出口壓力,本文設(shè)為1 個(gè)標(biāo)準(zhǔn)大氣壓;Pv是液體飽和蒸汽壓;Vthroat是喉部流速,在計(jì)算中,Vthroat取喉部截面的平均速度。

    文丘里管的幾何結(jié)構(gòu)為L(zhǎng)=8 mm(γ=0.67),β=10°,對(duì)其在工況為出口壓力0.1 MPa,入口壓力0.2~0.6 MPa之間進(jìn)行數(shù)值模擬。空化區(qū)域的體積大小Vc定義為

    其中,n為計(jì)算域控制體網(wǎng)格的數(shù)目,αi是控制體內(nèi)蒸汽相體積分?jǐn)?shù),Vi是控制體的體積。

    如圖7所示,空化數(shù)隨入口壓力的增大而逐漸減小,說(shuō)明空化程度逐漸增強(qiáng)。文丘里管喉部的最大速度和平均速度也因壓力增大而出現(xiàn)近似線性增加??栈w積Vc變化趨勢(shì)如圖8 所示,入口壓力的增大使得空化體積迅速增加,當(dāng)入口壓力為0.2 MPa 時(shí),Vc僅為19.56 mm3,此時(shí)空化汽泡僅存在于文丘里管喉部,產(chǎn)率極低,故后續(xù)分析以0.3 MPa 作為起始入口壓力;而當(dāng)入口壓力提升至0.5 MPa 以上時(shí),空化體積較0.2 MPa 的工況已有了3 個(gè)量級(jí)的提升。可見,入口壓力大小對(duì)空泡產(chǎn)量的影響是顯著的。

    圖7 不同入口壓力下空化數(shù)和速度變化Fig.7 Variation of cavitation number and velocity at different inlet pressures

    圖8不同入口壓力下空化體積變化Fig.8 Variation of cavitation volume at different inlet pressures

    為說(shuō)明新型文丘里管與一般文丘里管空化效果的差異,選取普通圓形截面文丘里管空化器作為比較對(duì)象,二者的幾何參數(shù)如表1。兩型空化器的喉道過(guò)流面積一致,以確保在相同入口壓力時(shí)喉部的速度大小和空化數(shù)相同。將兩型空化器在相同工況下進(jìn)行數(shù)值計(jì)算,采用的數(shù)值計(jì)算方法與前文一致。計(jì)算結(jié)果如圖8 所示,除了在入口壓力較低的0.2 MPa 之外,新型文丘里管的空化區(qū)域體積均大于普通文丘里管,在入口壓力大于0.3 MPa后,兩者的空化區(qū)域體積相差越來(lái)越顯著,體積增加率約為31.4%~72%。由此可見,新型文丘里管在空化效率方面有所提升,這得益于其喉部剪切區(qū)域的增加和喉部下游低壓區(qū)的增大。

    表1 兩型文丘里管幾何參數(shù)Tab.1 Geometric parameters of two-type venturies

    圖9 顯示了文丘里管內(nèi)部距離入口100~500 mm 區(qū)間的壓力分布情況,在收縮段液體的流速陡然增加,壓力因此快速降低,在喉部降到最小。發(fā)生空化現(xiàn)象后,壓力在空化區(qū)域保持為近似液體飽和蒸汽壓。當(dāng)空化潰滅后,在汽液兩相交接區(qū)域出現(xiàn)壓力波動(dòng),并在擴(kuò)張段逐漸恢復(fù)到約等于出口壓力。隨著入口壓力的增大,喉部下游的壓力恢復(fù)速率減慢,使得空化區(qū)域得以延伸。

    圖9 壓力分布曲線(100~500mm)Fig.9 Pressure distribution(100~500 mm)

    一方面,較大的入口壓力提升了空化體積;另一方面,文丘里管下游擴(kuò)張段內(nèi)液體的湍動(dòng)能大小也隨之增加。如圖10所示,喉部?jī)?nèi)嵌的分流體下游空化區(qū)域湍動(dòng)能大于擴(kuò)張段壁面處的空化區(qū)湍動(dòng)能,湍動(dòng)能極值出現(xiàn)在液體含汽率變化較快的潰滅區(qū)。湍流越劇烈潰滅強(qiáng)度越強(qiáng),越有利于提高·OH產(chǎn)率,并強(qiáng)化·OH擴(kuò)散到周圍的液體中參與降解反應(yīng)。在較高的入口壓力下,空泡的潰滅變得更加劇烈,從而在空泡潰滅時(shí)產(chǎn)生更高的壓力脈沖,水分子的解離程度增強(qiáng),產(chǎn)生更多的羥基自由基。

    圖10 對(duì)稱面(XY平面)汽相體積分?jǐn)?shù)與湍動(dòng)能分布云圖Fig.10 Vapor volume fraction and turbulence kinetic energy contours on symmetric plane (XY plane)

    3.3 不同擴(kuò)張角度的影響

    空化汽泡的潰滅發(fā)生在文丘里管喉部的下游,空化的總體效果取決于產(chǎn)生的空泡數(shù)量及其潰滅強(qiáng)度。擴(kuò)張段角度大小直接影響了喉部下游的沿流向壓力恢復(fù)速率和流動(dòng)分離程度,從而影響空化的發(fā)展與潰滅,故擴(kuò)張段幾何結(jié)構(gòu)是文丘里管設(shè)計(jì)中所關(guān)注的重點(diǎn)。一般來(lái)說(shuō),壓力恢復(fù)率隨著擴(kuò)張角的增加而增加,較大的擴(kuò)張角會(huì)導(dǎo)致空泡迅速潰滅,不利于空泡的發(fā)展;而對(duì)于較小的擴(kuò)張角,壓力恢復(fù)相對(duì)平穩(wěn),有利于空泡的生長(zhǎng)。但是擴(kuò)張角度過(guò)小時(shí),流過(guò)喉部的液體流動(dòng)分離減弱,導(dǎo)致下游湍動(dòng)能變小。而潰滅強(qiáng)度與湍動(dòng)能大小正相關(guān),因此過(guò)小的擴(kuò)張角度會(huì)削弱潰滅強(qiáng)度,也不利于·OH在液體中的擴(kuò)散。

    圖11 不同擴(kuò)張角度下的空化體積變化Fig.11 Variation of cavitation volume at different expansion angles

    圖12 不同擴(kuò)張角度下的汽相體積分?jǐn)?shù)云圖(入口壓力0.5 MPa)Fig.12 Vapor volume fraction contour at different expansion angles

    本文選擇幾何結(jié)構(gòu)為L(zhǎng)=8 mm,擴(kuò)張角度β=8°、10°、12°和14°的文丘里管在入口壓力為0.3~0.6 MPa區(qū)間內(nèi)進(jìn)行數(shù)值模擬,結(jié)果如圖11~12所示。由圖11可見,當(dāng)文丘里管擴(kuò)張角度為10°時(shí),在多數(shù)情況下空化體積均為最大值,故選擇10°為最佳擴(kuò)張角度。圖12以入口壓力為0.5 MPa為例展示了不同擴(kuò)張角度下的空化體積變化情況。擴(kuò)張角度的變化同時(shí)影響擴(kuò)張段壁面附近和梯形部件下游兩個(gè)空化區(qū)域的體積,主要是因?yàn)樵搮^(qū)域的壓力梯度主要體現(xiàn)在流動(dòng)方向上,而在徑向的壓力梯度變化并不顯著。由于在本文的數(shù)值模型中不考慮汽體的可壓縮性,根據(jù)連續(xù)性方程,擴(kuò)張角度越小,速度在流動(dòng)方向的衰減越緩慢。從機(jī)械能守恒的角度看,動(dòng)能大的高速的流動(dòng)區(qū)域環(huán)境壓力較低,維持了空化區(qū)域的發(fā)展。

    3.4 喉部分流體長(zhǎng)度(L)對(duì)于空化的影響

    在3.2 節(jié)和3.3 節(jié)中已經(jīng)可以清晰發(fā)現(xiàn)在喉部分流體的平臺(tái)前緣開始發(fā)生空化,并在其下游得到發(fā)展。分流體的幾何結(jié)構(gòu)可以影響喉部與下游的流動(dòng)結(jié)構(gòu),進(jìn)一步影響空化體積。平臺(tái)長(zhǎng)度影響空化在喉部的發(fā)展,本文選擇擴(kuò)張角度β=10°,L=4 mm、6 mm、7 mm、8 mm、9 mm、10 mm 進(jìn)行數(shù)值模擬,與之對(duì)應(yīng)的γ值在0.33~0.83之間,如圖13所示。

    圖13 不同喉部結(jié)構(gòu)(L=4~10 mm)Fig.13 Different throat structures(L=4-10 mm)

    不同喉部特征長(zhǎng)度下空化體積變化趨勢(shì)如圖14 所示,當(dāng)入口壓力小于0.3 MPa 時(shí),L的大小對(duì)空化體積的影響較小;反之,空化體積隨L的變化而發(fā)生顯著波動(dòng),且入口壓力越大空化體積變波動(dòng)越大,說(shuō)明在較低空化數(shù)條件下,空化體積的大小隨特征長(zhǎng)度L的變化更明顯。對(duì)于所有0.3 MPa 以上的工況,在L大于6 mm(γ>0.5)時(shí),空化體積先增大后減小,在達(dá)到8 mm 時(shí),空化體積達(dá)到峰值,即γ=0.67 時(shí),空化效果最佳;當(dāng)L=10 mm 時(shí)空化體積顯著下降。

    總的來(lái)說(shuō),存在一個(gè)最優(yōu)值L0使得空化體積在特定入口壓力下達(dá)到最大。當(dāng)L比L0小時(shí),空化發(fā)展并不充分,流動(dòng)在喉部發(fā)生的第一次分離位置過(guò)于靠前,較短的空化發(fā)展區(qū)間可能是降低總空化體積的原因;L比L0大時(shí),喉部與擴(kuò)張段的交匯處流動(dòng)面積變化更為劇烈,不利于空化氣泡在擴(kuò)張段的發(fā)展,該情況類似于孔板型空化器,L=9 mm、10 mm的分流體后緣存在近似直角的突變,雖然這種結(jié)構(gòu)可以形成更強(qiáng)的壓力突變,但其下游的壓力恢復(fù)速率較快,不利于空化繼續(xù)發(fā)展。在本文所計(jì)算的參數(shù)中,L0=8 mm。

    圖14 不同喉部特征長(zhǎng)度(L)下的空化體積變化趨勢(shì)Fig.14 Variation trend of cavitation volume at different throat characteristic lengths(L)

    4 結(jié) 論

    本文圍繞新型文丘里管的空化特性,基于RANS 方法在穩(wěn)態(tài)條件下對(duì)其進(jìn)行了數(shù)值模擬研究。以空泡潰滅時(shí)的能量轉(zhuǎn)化過(guò)程為切入點(diǎn),并以空化體積作為衡量該型文丘里管空化能力的表征參數(shù),圍繞入口壓力(0.2~0.6 MPa)、擴(kuò)張段擴(kuò)張角度(8°~14°)和喉部幾何特征(L=4~10mm,γ=0.33~0.83)進(jìn)行了比對(duì)分析,得到以下結(jié)論:

    (1)新型文丘里管與普通文丘里管相比,空化區(qū)域體積有所增大。除了在擴(kuò)張段壁面附近產(chǎn)生空化外,還在擴(kuò)張段的中心位置存在一個(gè)因分流部件而產(chǎn)生的空化區(qū),說(shuō)明下游產(chǎn)生了額外的低壓區(qū),這主要?dú)w因于分流部件使流體產(chǎn)生了額外的流動(dòng)分離。兩個(gè)空化區(qū)在流動(dòng)方向上的延伸長(zhǎng)度基本一致??栈w積與入口壓力大小呈正相關(guān),與擴(kuò)張角度大小呈負(fù)相關(guān),其特征長(zhǎng)度L與喉道長(zhǎng)度Lth的比值γ對(duì)空化產(chǎn)率具有較大影響,同一空化數(shù)下,空化體積隨γ的增大先增加后減小,且空化數(shù)越小,該參數(shù)的影響越顯著。

    (2)根據(jù)數(shù)值模擬結(jié)果,為達(dá)到最大空化體積,新型文丘里管的優(yōu)選幾何參數(shù)γ應(yīng)該為0.67,擴(kuò)張角度為10°,來(lái)流入口壓力為0.6 MPa。當(dāng)入口壓力足夠高時(shí),空化數(shù)進(jìn)一步下降,可能會(huì)出現(xiàn)超空化現(xiàn)象,將不利于降解效果的提升[16]。最佳入口壓力的大小還需要通過(guò)試驗(yàn)進(jìn)一步確定。

    (3)本文提出的新型文丘里管在一定程度上提升了傳統(tǒng)文丘里型水力空化發(fā)生器的空化產(chǎn)率,因此在污水降解方面該水力空化發(fā)生器將更具應(yīng)用價(jià)值。

    猜你喜歡
    文丘里喉部空泡
    文丘里洗滌器在凈化粗煤氣中的應(yīng)用及優(yōu)化
    云南化工(2023年6期)2023-07-04 01:35:20
    軸排凝汽器喉部設(shè)計(jì)
    水下航行體雙空泡相互作用數(shù)值模擬研究
    電子喉鏡聯(lián)合窄帶成像技術(shù)對(duì)喉部早期惡性病變的診斷價(jià)值研討
    高溫高壓飛灰過(guò)濾器文丘里管故障原因分析及應(yīng)對(duì)措施
    化工管理(2020年10期)2020-04-30 10:20:36
    文丘里管在國(guó)六排放標(biāo)準(zhǔn)中的應(yīng)用分析
    北京汽車(2019年3期)2019-07-19 01:44:08
    凝汽器喉部流場(chǎng)數(shù)值模擬
    基于LPV的超空泡航行體H∞抗飽和控制
    基于CFD的對(duì)轉(zhuǎn)槳無(wú)空泡噪聲的仿真預(yù)報(bào)
    船海工程(2015年4期)2016-01-05 15:53:28
    咽及喉部鱗癌放射治療技術(shù)研究進(jìn)展
    免费高清在线观看视频在线观看| 国产黄色视频一区二区在线观看| 中文精品一卡2卡3卡4更新| 国产熟女欧美一区二区| h视频一区二区三区| 国产欧美日韩一区二区三区在线 | 最后的刺客免费高清国语| 久久ye,这里只有精品| 日韩精品有码人妻一区| 男的添女的下面高潮视频| 人人妻人人澡人人爽人人夜夜| av黄色大香蕉| 在线免费观看不下载黄p国产| 国产高清有码在线观看视频| 少妇人妻久久综合中文| 黑丝袜美女国产一区| 成人一区二区视频在线观看| a级毛色黄片| 啦啦啦中文免费视频观看日本| 色婷婷av一区二区三区视频| 永久免费av网站大全| 午夜免费男女啪啪视频观看| 蜜桃亚洲精品一区二区三区| 老司机影院成人| 黑人高潮一二区| 久久精品国产自在天天线| 久久久午夜欧美精品| 成人免费观看视频高清| 国产91av在线免费观看| 一本久久精品| 色哟哟·www| 国产午夜精品久久久久久一区二区三区| 少妇高潮的动态图| 一级爰片在线观看| 免费黄色在线免费观看| 91久久精品国产一区二区成人| 免费久久久久久久精品成人欧美视频 | 大话2 男鬼变身卡| 国产黄片美女视频| 国产在线男女| 五月天丁香电影| 国产在线免费精品| 亚洲av在线观看美女高潮| 色婷婷久久久亚洲欧美| 亚洲,一卡二卡三卡| 亚洲av电影在线观看一区二区三区| 欧美精品一区二区免费开放| 人体艺术视频欧美日本| 欧美变态另类bdsm刘玥| 熟妇人妻不卡中文字幕| 麻豆成人午夜福利视频| 国产精品精品国产色婷婷| 国产精品av视频在线免费观看| 亚洲av日韩在线播放| 一个人看视频在线观看www免费| 蜜桃久久精品国产亚洲av| 婷婷色综合www| 在线观看免费高清a一片| 日日摸夜夜添夜夜添av毛片| 精品人妻视频免费看| 国产真实伦视频高清在线观看| 精品人妻偷拍中文字幕| 日本免费在线观看一区| 久久av网站| 97在线人人人人妻| 一级黄片播放器| 国产免费一区二区三区四区乱码| 久久久精品94久久精品| 精品久久久久久电影网| 亚洲av国产av综合av卡| 成人亚洲欧美一区二区av| 少妇的逼水好多| 在线观看人妻少妇| 天堂俺去俺来也www色官网| 美女cb高潮喷水在线观看| 在线观看美女被高潮喷水网站| 亚洲综合色惰| 街头女战士在线观看网站| 国产熟女欧美一区二区| 不卡视频在线观看欧美| 麻豆成人午夜福利视频| 国产精品一及| 日韩不卡一区二区三区视频在线| 少妇丰满av| 观看av在线不卡| 国产精品国产三级专区第一集| 久久国产精品大桥未久av | 国产精品国产三级国产专区5o| 亚洲欧美日韩另类电影网站 | 中文字幕免费在线视频6| 久久99热6这里只有精品| a 毛片基地| av国产久精品久网站免费入址| 校园人妻丝袜中文字幕| 国产综合精华液| 国产亚洲午夜精品一区二区久久| 少妇人妻久久综合中文| 日本午夜av视频| 久久国产乱子免费精品| 亚洲怡红院男人天堂| 免费不卡的大黄色大毛片视频在线观看| 少妇丰满av| 国产一区二区在线观看日韩| 亚洲欧洲国产日韩| 亚洲熟女精品中文字幕| 国产午夜精品久久久久久一区二区三区| 91久久精品电影网| 九九在线视频观看精品| 黑丝袜美女国产一区| 国产精品精品国产色婷婷| 天美传媒精品一区二区| 国产高清有码在线观看视频| 大片免费播放器 马上看| 国国产精品蜜臀av免费| 国产亚洲91精品色在线| 成人国产av品久久久| 中文精品一卡2卡3卡4更新| 国产精品秋霞免费鲁丝片| 少妇人妻一区二区三区视频| 男女边摸边吃奶| 两个人的视频大全免费| 九九在线视频观看精品| 女人久久www免费人成看片| 久久亚洲国产成人精品v| 男人添女人高潮全过程视频| 97在线人人人人妻| 国产久久久一区二区三区| 干丝袜人妻中文字幕| 麻豆乱淫一区二区| 国产精品99久久99久久久不卡 | 国产亚洲精品久久久com| 一区二区三区免费毛片| 亚洲aⅴ乱码一区二区在线播放| 欧美一区二区亚洲| 国精品久久久久久国模美| 亚洲精品久久午夜乱码| 日日摸夜夜添夜夜爱| 一区二区三区四区激情视频| av卡一久久| av女优亚洲男人天堂| 亚洲精品日本国产第一区| 欧美xxxx黑人xx丫x性爽| 国产成人免费无遮挡视频| 亚洲精品久久午夜乱码| 亚洲国产最新在线播放| 王馨瑶露胸无遮挡在线观看| 高清不卡的av网站| 久久97久久精品| 天美传媒精品一区二区| 亚洲国产精品一区三区| 多毛熟女@视频| 国产精品欧美亚洲77777| 少妇的逼好多水| 国产男女内射视频| 高清av免费在线| 午夜激情久久久久久久| 国产亚洲91精品色在线| 欧美精品国产亚洲| 久久人人爽人人爽人人片va| 国产欧美另类精品又又久久亚洲欧美| 亚洲伊人久久精品综合| 欧美bdsm另类| 最黄视频免费看| 91午夜精品亚洲一区二区三区| 在线观看av片永久免费下载| 中文字幕av成人在线电影| 男女边摸边吃奶| 亚洲精品久久久久久婷婷小说| 2022亚洲国产成人精品| 国产极品天堂在线| 自拍偷自拍亚洲精品老妇| 久久99蜜桃精品久久| 91aial.com中文字幕在线观看| 久久99蜜桃精品久久| 成人黄色视频免费在线看| 亚洲成人中文字幕在线播放| 亚洲怡红院男人天堂| 日日撸夜夜添| av在线播放精品| 好男人视频免费观看在线| 麻豆成人午夜福利视频| 国产精品熟女久久久久浪| 王馨瑶露胸无遮挡在线观看| 啦啦啦视频在线资源免费观看| 卡戴珊不雅视频在线播放| 国产成人免费无遮挡视频| 精品久久久久久久末码| 午夜老司机福利剧场| 99久国产av精品国产电影| 国产老妇伦熟女老妇高清| 丝袜脚勾引网站| 日韩 亚洲 欧美在线| 青春草国产在线视频| 国产精品久久久久久久久免| 亚洲图色成人| 在线观看一区二区三区激情| 国产欧美日韩精品一区二区| 日韩在线高清观看一区二区三区| 成人高潮视频无遮挡免费网站| 肉色欧美久久久久久久蜜桃| 91午夜精品亚洲一区二区三区| 下体分泌物呈黄色| 欧美丝袜亚洲另类| 国产色爽女视频免费观看| 18+在线观看网站| .国产精品久久| 国产爱豆传媒在线观看| 国产精品偷伦视频观看了| 国产伦精品一区二区三区视频9| 2018国产大陆天天弄谢| 久久99热这里只频精品6学生| 欧美日韩国产mv在线观看视频 | 国产毛片在线视频| 亚洲成色77777| 午夜福利高清视频| 亚洲av国产av综合av卡| www.色视频.com| 国产精品精品国产色婷婷| 一区二区三区免费毛片| 久久精品夜色国产| 欧美 日韩 精品 国产| 国产日韩欧美在线精品| 亚洲最大成人中文| 成人亚洲欧美一区二区av| 我的老师免费观看完整版| 久久久午夜欧美精品| av免费观看日本| av免费在线看不卡| 一级片'在线观看视频| 久久影院123| av卡一久久| 国产淫片久久久久久久久| 国产乱来视频区| 亚洲精品乱码久久久v下载方式| 国产一区有黄有色的免费视频| 99久久精品一区二区三区| 久久久久久久久久成人| 国产男女超爽视频在线观看| 久久精品人妻少妇| 精品久久久久久电影网| 深爱激情五月婷婷| 午夜福利影视在线免费观看| 啦啦啦啦在线视频资源| 中文在线观看免费www的网站| 日本欧美视频一区| 人人妻人人看人人澡| 亚洲久久久国产精品| 欧美最新免费一区二区三区| 精品人妻偷拍中文字幕| 哪个播放器可以免费观看大片| av免费在线看不卡| 丰满迷人的少妇在线观看| 免费观看无遮挡的男女| 国产乱来视频区| 老熟女久久久| 一个人看视频在线观看www免费| 国产爱豆传媒在线观看| 亚洲怡红院男人天堂| 久久人人爽人人爽人人片va| 色视频在线一区二区三区| 日本欧美国产在线视频| 日本黄色片子视频| 亚洲三级黄色毛片| 国产高清不卡午夜福利| 哪个播放器可以免费观看大片| 色吧在线观看| 深爱激情五月婷婷| 我要看黄色一级片免费的| 青春草国产在线视频| 国产日韩欧美亚洲二区| kizo精华| 在线天堂最新版资源| 亚洲精品日韩av片在线观看| 国产精品成人在线| 国产男女内射视频| 一本久久精品| 中文字幕免费在线视频6| 久久鲁丝午夜福利片| 久久精品人妻少妇| 亚洲综合色惰| 久久久精品94久久精品| 丰满少妇做爰视频| 成年免费大片在线观看| 中文在线观看免费www的网站| 人人妻人人爽人人添夜夜欢视频 | 国产精品国产三级国产专区5o| 国产高清三级在线| 国产精品一区二区三区四区免费观看| 中文字幕av成人在线电影| 丰满人妻一区二区三区视频av| 久久 成人 亚洲| 日韩 亚洲 欧美在线| 国产日韩欧美在线精品| 黄色一级大片看看| 亚洲精品视频女| 99久久精品热视频| 18禁在线播放成人免费| 亚洲欧美清纯卡通| 欧美日韩亚洲高清精品| 国产在线男女| 肉色欧美久久久久久久蜜桃| 观看av在线不卡| 国产视频首页在线观看| 午夜福利影视在线免费观看| 蜜桃亚洲精品一区二区三区| 精品久久久噜噜| 三级经典国产精品| 国产午夜精品久久久久久一区二区三区| 大码成人一级视频| 日日摸夜夜添夜夜爱| 99热全是精品| 精品一区二区三卡| 久久久久久久久久久免费av| 午夜福利在线在线| 久久人人爽人人爽人人片va| 只有这里有精品99| 久久久久人妻精品一区果冻| 国产精品熟女久久久久浪| 男女无遮挡免费网站观看| 你懂的网址亚洲精品在线观看| 精品国产乱码久久久久久小说| 毛片一级片免费看久久久久| 国产精品伦人一区二区| 成人午夜精彩视频在线观看| 看非洲黑人一级黄片| 多毛熟女@视频| 五月天丁香电影| 国产在线一区二区三区精| 国产精品一区二区三区四区免费观看| 精品人妻视频免费看| 久久久久视频综合| 成人无遮挡网站| 高清毛片免费看| 免费大片黄手机在线观看| 婷婷色麻豆天堂久久| 黑人猛操日本美女一级片| 久久久色成人| 老女人水多毛片| 美女脱内裤让男人舔精品视频| 少妇人妻 视频| 亚洲婷婷狠狠爱综合网| 尤物成人国产欧美一区二区三区| 夫妻性生交免费视频一级片| 91aial.com中文字幕在线观看| 2018国产大陆天天弄谢| 国产精品久久久久成人av| 99热这里只有是精品50| freevideosex欧美| 搡老乐熟女国产| 欧美精品亚洲一区二区| 国产精品国产av在线观看| 午夜日本视频在线| 中国国产av一级| 建设人人有责人人尽责人人享有的 | videossex国产| 99热这里只有精品一区| 国产在线免费精品| 亚洲欧美成人综合另类久久久| 99久久综合免费| 中文字幕人妻熟人妻熟丝袜美| 久久国产精品大桥未久av | 国产精品一区二区性色av| 国产精品一及| 老司机影院成人| 在线观看免费视频网站a站| 99精国产麻豆久久婷婷| 国产精品国产三级国产av玫瑰| 我要看黄色一级片免费的| 又爽又黄a免费视频| 夫妻性生交免费视频一级片| 少妇被粗大猛烈的视频| 少妇裸体淫交视频免费看高清| 97精品久久久久久久久久精品| 成人综合一区亚洲| 欧美区成人在线视频| 亚洲怡红院男人天堂| 免费观看性生交大片5| 最近2019中文字幕mv第一页| 中文在线观看免费www的网站| 人人妻人人澡人人爽人人夜夜| 日韩伦理黄色片| 国精品久久久久久国模美| 亚洲精品中文字幕在线视频 | 亚洲精品日本国产第一区| 又黄又爽又刺激的免费视频.| 亚洲成人av在线免费| 精品久久久久久电影网| 美女福利国产在线 | 国产欧美亚洲国产| 国产精品av视频在线免费观看| 99久久综合免费| av在线老鸭窝| 国内少妇人妻偷人精品xxx网站| 成人特级av手机在线观看| 久热久热在线精品观看| 色综合色国产| 26uuu在线亚洲综合色| 高清不卡的av网站| 欧美高清性xxxxhd video| 久久国产精品大桥未久av | 亚洲精品自拍成人| 欧美另类一区| 午夜激情福利司机影院| 蜜桃在线观看..| 2022亚洲国产成人精品| 久久久久人妻精品一区果冻| 建设人人有责人人尽责人人享有的 | 新久久久久国产一级毛片| 天堂中文最新版在线下载| 国产亚洲午夜精品一区二区久久| 韩国av在线不卡| 日本黄色日本黄色录像| 少妇被粗大猛烈的视频| 91精品国产国语对白视频| 久久久久久人妻| 熟女电影av网| 中文精品一卡2卡3卡4更新| 亚洲欧美中文字幕日韩二区| 亚洲av中文字字幕乱码综合| 汤姆久久久久久久影院中文字幕| 国产精品精品国产色婷婷| 嘟嘟电影网在线观看| 亚洲中文av在线| 国产成人免费观看mmmm| 色婷婷久久久亚洲欧美| kizo精华| 精品国产三级普通话版| 日韩免费高清中文字幕av| 午夜免费观看性视频| 国产真实伦视频高清在线观看| av国产久精品久网站免费入址| 最近最新中文字幕大全电影3| 男的添女的下面高潮视频| 久久人人爽av亚洲精品天堂 | 亚洲精品乱码久久久久久按摩| 久久久久久久精品精品| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产高清在线一区二区三| 亚洲自偷自拍三级| 偷拍熟女少妇极品色| 国产乱人偷精品视频| 国内少妇人妻偷人精品xxx网站| 久久影院123| 少妇精品久久久久久久| 国产黄频视频在线观看| 精品久久久久久久末码| 精品视频人人做人人爽| 观看av在线不卡| 毛片女人毛片| 欧美激情国产日韩精品一区| 美女国产视频在线观看| 亚洲欧美中文字幕日韩二区| 久久av网站| 少妇人妻久久综合中文| 欧美精品亚洲一区二区| 夫妻性生交免费视频一级片| 国产精品成人在线| 毛片一级片免费看久久久久| 免费黄频网站在线观看国产| 国产高潮美女av| 18禁在线播放成人免费| 九九爱精品视频在线观看| 欧美xxⅹ黑人| 国产精品一区二区在线不卡| 国产欧美日韩一区二区三区在线 | 成人漫画全彩无遮挡| 国产 一区 欧美 日韩| 国产精品偷伦视频观看了| videossex国产| 国产欧美日韩精品一区二区| 少妇丰满av| 在线免费十八禁| 国产精品一区二区三区四区免费观看| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲精品成人av观看孕妇| 精品酒店卫生间| 九色成人免费人妻av| 一级av片app| 久久人妻熟女aⅴ| 天堂8中文在线网| 国产高清三级在线| 91aial.com中文字幕在线观看| 最近手机中文字幕大全| 中国美白少妇内射xxxbb| 我的女老师完整版在线观看| 久久精品夜色国产| 一个人看的www免费观看视频| 有码 亚洲区| 国产91av在线免费观看| 男人舔奶头视频| 国产免费视频播放在线视频| 日日啪夜夜爽| 午夜福利在线在线| 日韩成人av中文字幕在线观看| 性色av一级| 亚洲av欧美aⅴ国产| 成人无遮挡网站| 亚洲,欧美,日韩| 国产一区亚洲一区在线观看| 欧美成人精品欧美一级黄| 精品人妻偷拍中文字幕| 美女福利国产在线 | 亚洲在久久综合| 亚洲国产欧美人成| 欧美高清性xxxxhd video| 日产精品乱码卡一卡2卡三| 一区二区三区乱码不卡18| 欧美日本视频| 久久久久视频综合| 亚洲av不卡在线观看| 在线看a的网站| 天天躁日日操中文字幕| 久久久欧美国产精品| 97超视频在线观看视频| 国产精品久久久久久精品古装| 美女脱内裤让男人舔精品视频| 在线观看国产h片| 久久国产精品男人的天堂亚洲 | av播播在线观看一区| 五月天丁香电影| av天堂中文字幕网| 视频中文字幕在线观看| av专区在线播放| 成人国产av品久久久| 免费观看在线日韩| 91久久精品国产一区二区成人| 成人午夜精彩视频在线观看| 色吧在线观看| 精品99又大又爽又粗少妇毛片| 日韩av在线免费看完整版不卡| 春色校园在线视频观看| 女人久久www免费人成看片| 伦理电影免费视频| 日本欧美视频一区| 亚洲伊人久久精品综合| 欧美成人精品欧美一级黄| 日韩av不卡免费在线播放| 一本久久精品| av天堂中文字幕网| 91在线精品国自产拍蜜月| 国产在线男女| 欧美精品国产亚洲| 99热这里只有是精品50| 九九久久精品国产亚洲av麻豆| 国产亚洲一区二区精品| 亚洲精品aⅴ在线观看| 欧美日韩国产mv在线观看视频 | 菩萨蛮人人尽说江南好唐韦庄| 香蕉精品网在线| 人妻夜夜爽99麻豆av| 一级毛片我不卡| 国产白丝娇喘喷水9色精品| 久久人人爽av亚洲精品天堂 | 亚洲欧洲国产日韩| 欧美成人精品欧美一级黄| 街头女战士在线观看网站| 日韩 亚洲 欧美在线| 亚洲欧美一区二区三区黑人 | 黑丝袜美女国产一区| 国产黄片视频在线免费观看| 又爽又黄a免费视频| 亚州av有码| 国产在线男女| 男女免费视频国产| 久久久久网色| 欧美精品人与动牲交sv欧美| 色婷婷久久久亚洲欧美| 日本-黄色视频高清免费观看| 高清欧美精品videossex| 亚洲精品aⅴ在线观看| 美女高潮的动态| 国产黄色免费在线视频| 国产亚洲一区二区精品| freevideosex欧美| 一级av片app| 欧美日韩精品成人综合77777| 丝袜脚勾引网站| 2022亚洲国产成人精品| 热re99久久精品国产66热6| 草草在线视频免费看| 久久久久久久大尺度免费视频| 最黄视频免费看| 少妇精品久久久久久久| 2018国产大陆天天弄谢| 国产免费又黄又爽又色| 丰满人妻一区二区三区视频av| 老女人水多毛片| 日本猛色少妇xxxxx猛交久久| 99久久精品一区二区三区| 校园人妻丝袜中文字幕| 欧美区成人在线视频| 免费观看的影片在线观看| 男女啪啪激烈高潮av片| 婷婷色av中文字幕| 五月伊人婷婷丁香| 99久久人妻综合| 各种免费的搞黄视频| 国产精品国产av在线观看| 久久精品久久精品一区二区三区| 美女脱内裤让男人舔精品视频| 日韩 亚洲 欧美在线| 最黄视频免费看| 久久国产精品男人的天堂亚洲 | 我的女老师完整版在线观看| 国产亚洲av片在线观看秒播厂| 黄片wwwwww| 精品亚洲成a人片在线观看 | 久久久久久人妻| 少妇人妻 视频| av黄色大香蕉| 爱豆传媒免费全集在线观看| 日本欧美视频一区| 麻豆国产97在线/欧美| kizo精华| 久久久久人妻精品一区果冻|