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

    氣升式陶瓷膜過濾過程的氣液兩相流模擬

    2016-07-07 12:12:04林進(jìn)沈浩景文珩南京工業(yè)大學(xué)材料化學(xué)工程國家重點(diǎn)實(shí)驗(yàn)室江蘇南京210009
    化工學(xué)報(bào) 2016年6期

    林進(jìn),沈浩,景文珩(南京工業(yè)大學(xué)材料化學(xué)工程國家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210009)

    ?

    氣升式陶瓷膜過濾過程的氣液兩相流模擬

    林進(jìn),沈浩,景文珩
    (南京工業(yè)大學(xué)材料化學(xué)工程國家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210009)

    摘要:采用VOF雙流體模型對19通道氣升式陶瓷膜過濾裝置進(jìn)行氣液兩相流的流體動(dòng)力學(xué)模擬,研究了曝氣孔直徑和曝氣量對氣升式陶瓷膜過濾裝置的氣含率、環(huán)流液速、膜面剪切力及膜管內(nèi)湍流強(qiáng)度的影響,模擬結(jié)果與實(shí)驗(yàn)結(jié)果的誤差在5%~10%之間。結(jié)果表明,氣升管與降液管的氣含率都隨曝氣量增大而增大,隨曝氣孔直徑減小而增大;環(huán)流液速、膜面剪切力及膜管內(nèi)的湍流強(qiáng)度都隨曝氣量增大先增大,當(dāng)曝氣量達(dá)到400 L·h?1時(shí)其增大趨勢變緩。通過實(shí)驗(yàn)和模擬比較了3種不同孔徑的曝氣頭,環(huán)流液速與曝氣孔的直徑關(guān)系不大,僅與曝氣量相關(guān),但曝氣孔直徑越小,其膜面剪切力越大,越有利于過濾過程的進(jìn)行。

    關(guān)鍵詞:氣升式;陶瓷膜;計(jì)算流體力學(xué);VOF方法;兩相流;模擬

    2015-10-15收到初稿, 2016-02-24收到修改稿。

    聯(lián)系人:景文珩。第一作者:林進(jìn)(1989—),男,碩士研究生。

    Received date: 2015-10-15.

    Foundation item: supported by the National Natural Science Foundation of China (21176116), the Higher Education Natural Science Foundation of Jiangsu Province (15KJA530001) and the Project of Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

    引 言

    氣升式陶瓷膜過濾裝置是基于氣升式環(huán)流反應(yīng)器開發(fā)出來的新型過濾裝置,其典型的特點(diǎn)是采用氣液兩相流流動(dòng)替代傳統(tǒng)的液流過濾系統(tǒng)。通過向裝置內(nèi)通入壓縮氣體,使氣升管與降液管之間形成密度差異,推動(dòng)裝置內(nèi)的流體循環(huán)流動(dòng),并在膜管內(nèi)形成氣液兩相流,提高膜表面的湍流程度和剪切力,抑制膜表面濃差極化,降低膜面污染[1]。由于氣體密度遠(yuǎn)小于液體密度,能大幅降低過程能耗,而且具有結(jié)構(gòu)簡單、氣液傳質(zhì)效果好等優(yōu)點(diǎn),擁有重要的應(yīng)用前景與實(shí)用價(jià)值[2]。

    計(jì)算流體力學(xué)(computational fluid dynamics, CFD)是流體力學(xué)的一個(gè)分支。與傳統(tǒng)的實(shí)驗(yàn)方法相比,CFD 技術(shù)具有花費(fèi)少、設(shè)計(jì)周期短等優(yōu)點(diǎn),而且可以完整地獲取速度場等詳細(xì)的信息,對裝備的設(shè)計(jì)、選型、優(yōu)化等基礎(chǔ)研究都具有重要意義[3]。

    氣液兩相流由于具有清晰的相界面,對其進(jìn)行的模擬研究一般采用VOF(volume of fluid,流體體積)模型。賀蕭等[4]利用VOF模型對垂直上升管內(nèi)的氣液兩相流動(dòng)進(jìn)行模擬研究,預(yù)測了垂直及傾斜上升管內(nèi)Taylor氣泡形狀和壁面切應(yīng)力等流動(dòng)參數(shù)。Kulpers等[5]使用二維VOF模型與三維界面追蹤模型相比較,研究了密度比等物理特性對不同直徑下的氣泡上升運(yùn)動(dòng)及其最終上升速度的影響。Abid等[6]采用VOF模型模擬圓柱形反應(yīng)器內(nèi)連續(xù)氣泡流的流動(dòng)情況,重點(diǎn)研究氣泡大小與表觀氣速和曝氣孔直徑的關(guān)系,通過模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)的對比驗(yàn)證了該模型可準(zhǔn)確預(yù)測反應(yīng)器內(nèi)兩相流的分布規(guī)律。

    本工作采用VOF方法模擬氣升式陶瓷膜過濾裝置內(nèi)的氣液兩相流動(dòng),考察不同曝氣量及不同曝氣孔直徑對裝置內(nèi)部的氣含率與環(huán)流液速的影響,重點(diǎn)研究了膜管內(nèi)膜面剪切力與膜管內(nèi)湍流強(qiáng)度的變化。本研究結(jié)果可為氣升式陶瓷膜過濾裝置曝氣頭的選型與曝氣量的選擇提供依據(jù)。

    1 實(shí)驗(yàn)部分

    實(shí)驗(yàn)裝置如圖1所示,氣升管直徑40 mm、高470 mm,降液管直徑45 mm、高480 mm,氣液分離區(qū)直徑60 mm、高310 mm。膜組件為19通道陶瓷膜管,長250 mm、管徑30 mm,單個(gè)通道孔徑3.5 mm。曝氣頭長50 mm、寬40 mm。裝置內(nèi)液相為水,氣相為壓縮空氣。

    實(shí)驗(yàn)流程:打開閥門F1~F3,關(guān)閉閥門F4~F6,使用計(jì)量泵將水樣打入裝置內(nèi)部,待水位到達(dá)指定高度,關(guān)閉泵和F1;打開閥門F7和F8;用閥門F7調(diào)節(jié)曝氣量到一定值,壓縮空氣由氣升管底部的曝氣頭進(jìn)入。隨著氣體的進(jìn)入,氣升管與降液管之間形成密度差,驅(qū)使過濾器底部的液體由氣升管依次流經(jīng)上部連接管-氣液分離區(qū)-降液管,最后經(jīng)下部連接管回到氣升管底部,形成環(huán)流。氣體則通過氣液分離區(qū),最后經(jīng)過流量計(jì)流出,并讀出示數(shù)。

    傳統(tǒng)的氣升式反應(yīng)器大都采用1~2 mm的常規(guī)曝氣孔,在氣體徑向分布性等方面暴露出不足,而新型的微孔曝氣器采用微米級曝氣孔,產(chǎn)生的氣泡直徑小,氣泡停留時(shí)間長,而且分布均勻[7]。故本研究選取直徑為1、0.5、0.2 mm的曝氣頭,比較常規(guī)曝氣與微孔曝氣對氣升式陶瓷膜過濾裝置相關(guān)性能的影響。每個(gè)曝氣頭根據(jù)曝氣量的大小分為100、200、300、400、500、600 L·h?16組進(jìn)行模擬。

    2 計(jì)算模型及參數(shù)設(shè)定

    2.1裝置網(wǎng)格

    裝置網(wǎng)格如圖2所示,對實(shí)驗(yàn)裝置進(jìn)行1:1尺寸的建模。由于實(shí)驗(yàn)裝置結(jié)構(gòu)復(fù)雜,難以采用結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分,故采用非結(jié)構(gòu)化網(wǎng)格對其進(jìn)行劃分。考慮到本實(shí)驗(yàn)重點(diǎn)關(guān)注膜管及氣升管內(nèi)的氣液兩相流動(dòng)情況,故對膜管及氣升管內(nèi)的網(wǎng)格進(jìn)行加密,以求更好地對流場進(jìn)行分析。

    圖2 實(shí)驗(yàn)裝置網(wǎng)格Fig.2 Grid diagram of experimental equipment

    2.2多相流模型設(shè)置

    本研究選取VOF模型為多相流模型進(jìn)行模擬。其主要原理為:在整個(gè)流場中定義一種流體體積與網(wǎng)格體積比值為α的函數(shù),其中α=1處的網(wǎng)格充滿該流體,α=0處的網(wǎng)格則不含有該流體,在0<α<1的網(wǎng)格區(qū)域含有自由面。在不同的時(shí)間下對流場的體積分?jǐn)?shù)函數(shù)進(jìn)行求解,研究α的分布情況,并使用相關(guān)的重構(gòu)運(yùn)動(dòng)界面方法對流體界面進(jìn)行重構(gòu),從而得到準(zhǔn)確的界面位置,進(jìn)而追蹤流體界面的變化[8]。

    2.2.1連續(xù)性方程和動(dòng)量方程連續(xù)性方程和動(dòng)量方程為

    式中,F(xiàn)i是動(dòng)量方程源項(xiàng),τij是層流剪應(yīng)力張量,? ρu′iu′j是湍流雷諾應(yīng)力張量。

    本數(shù)值模擬采用RNG k-ε模型(重整化群k-ε模型)對動(dòng)量方程進(jìn)行封閉。

    湍流動(dòng)能k的輸運(yùn)方程

    湍流耗散率ε的輸運(yùn)方程

    式中,Sk和Sε是用戶定義的源項(xiàng),Gk表示由平均速度梯度產(chǎn)生的湍流動(dòng)能,Gb表示由浮力產(chǎn)生的湍流動(dòng)能。

    2.2.2體積分?jǐn)?shù)方程在VOF模型中,相與相之間沒有相互穿插,跟蹤相與相之間的界面是通過求解一相或多相的體積分?jǐn)?shù)的連續(xù)方程實(shí)現(xiàn)的[9]。在此

    設(shè)α為計(jì)算單元內(nèi)第i項(xiàng)的體積分?jǐn)?shù)。氣液相界面的體積分?jǐn)?shù)方程如下

    液相體積分?jǐn)?shù)通過式(6)求得

    2.2.3表面張力項(xiàng)Fluent中的表面張力模型是由Brackbill等[10]提出的連續(xù)表面張力模型,通過這個(gè)模型,VOF計(jì)算中附加的表面張力導(dǎo)致動(dòng)量方程式(2)中動(dòng)量源項(xiàng)Fi的變化。

    表面張力的動(dòng)量源項(xiàng)表示為

    式中,κ是界面曲率,可用自由表面處單位曲面法向量的散度表示。

    2.3求解參數(shù)及邊界條件的設(shè)置

    采用速度進(jìn)口條件及壓力出口條件。壓力-速度耦合方程求解采用PISO方法;方程離散化時(shí),時(shí)間項(xiàng)采用隱式格式,壓力項(xiàng)采用PRESTO!算法,湍流動(dòng)能與湍流耗散率均采用二級迎風(fēng)格式;對氣液相界面的追蹤采用具有較高精度的幾何重構(gòu)模型;非穩(wěn)態(tài)方程采用時(shí)間步長0.0005 s,直到獲得穩(wěn)定流場。計(jì)算前對流場進(jìn)行初始化,裝置中水位高度為550 mm。

    3 數(shù)值模擬結(jié)果討論

    3.1網(wǎng)格無關(guān)性驗(yàn)證

    對于非穩(wěn)態(tài)問題的數(shù)值模擬,有必要進(jìn)行計(jì)算所用的網(wǎng)格數(shù)量與計(jì)算結(jié)果之間關(guān)聯(lián)性的驗(yàn)證,即網(wǎng)格無關(guān)性驗(yàn)證。對1 mm曝氣孔直徑的過濾裝置模型劃分不同的網(wǎng)格數(shù)量,分別為19000、68000、120000、220000個(gè),考察氣升管內(nèi)氣含率隨時(shí)間變化的曲線在不同的網(wǎng)格數(shù)下是否有明顯差別。

    由圖3可知,網(wǎng)格數(shù)從19000個(gè)增加至220000個(gè)過程中氣升管內(nèi)氣含率隨時(shí)間變化不大;由表1可知,穩(wěn)定流場內(nèi)氣升管的平均氣含率差值變化很小。從以上結(jié)果可以看出此時(shí)網(wǎng)格數(shù)對計(jì)算結(jié)果的影響很小,可以認(rèn)為19000的網(wǎng)格數(shù)已經(jīng)達(dá)到網(wǎng)格無關(guān),故取19000作為計(jì)算網(wǎng)格。

    圖3 網(wǎng)格數(shù)對氣升管氣含率的影響Fig.3 Effect of grid number on gas hold-up of riser

    表1 網(wǎng)格數(shù)對氣升管平均氣含率的影響Table 1 Effect of grid number on average gas hold-up of riser

    3.2模擬結(jié)果與實(shí)驗(yàn)結(jié)果對比

    圖4 實(shí)驗(yàn)值與模擬值的對比Fig.4 Comparison of simulated and experimental data

    對數(shù)值模擬計(jì)算而言,需將模擬值與實(shí)驗(yàn)值進(jìn)行比較,以確定模擬實(shí)驗(yàn)的準(zhǔn)確性。如圖4所示,本實(shí)驗(yàn)的模擬值與實(shí)驗(yàn)值的誤差在5%~10%之間,可以認(rèn)為模擬結(jié)果可信。對實(shí)驗(yàn)值與模擬值進(jìn)行分析,回歸出一元線性方程,相關(guān)度分別為0.9760和0.9970,表明該方程可以較好地預(yù)測氣含率與曝氣量之間的關(guān)系,并為估算和預(yù)測氣含率提供了一種更為方便簡捷的方法。

    3.3膜管內(nèi)流場分析

    由圖5可知,氣升管內(nèi)的氣含率高于降液管及氣液分離區(qū)內(nèi)的氣含率。膜管內(nèi)的連續(xù)彈狀流由氣彈和液彈組成,而且氣彈與液彈的長度不一。膜管內(nèi)單個(gè)通道的直徑為3.5 mm,根據(jù)Mehendale等[11]的劃分標(biāo)準(zhǔn)屬于小通道。小通道內(nèi)表面效應(yīng)占主導(dǎo)地位,與常規(guī)通道相比有很大的不同:由于表面張力的作用,氣泡的首尾呈類似彈頭的形狀,氣泡直徑約等于管徑,長度大于管徑。在彈狀流流型下,氣彈的上升速度大于液彈,導(dǎo)致氣彈前端的液體以液膜的形式向下方流動(dòng),液膜厚度約為幾十微米[12],氣彈的尾端和液滴的射流作用形成尾渦。在氣彈區(qū)域剪切力與液彈區(qū)域相反,氣彈區(qū)域的剪切力較大,由于尾渦處于湍流狀態(tài),剪切力波動(dòng)較大[13]。

    圖5 裝置整體氣含率Fig.5 Diagram of gas hold-up in equipment

    由圖6可知,氣升管及氣液分離區(qū)的流體流速較快,在降液管內(nèi)流速較慢。下部連接管回流的液體對上升氣泡產(chǎn)生橫向的沖擊,并在膜管下端產(chǎn)生強(qiáng)烈湍流,使氣泡發(fā)生強(qiáng)烈的破碎與聚并。上部連接管進(jìn)入降液管的氣體大部分進(jìn)入氣液分離區(qū),進(jìn)而與液相分離,另一部分隨液相環(huán)流進(jìn)入降液管。膜管內(nèi)流體流速較快,而且速度分布均勻,不存在偏流與局部回流的現(xiàn)象。

    圖6 裝置整體速度矢量圖Fig.6 Diagram of velocity vector in equipment

    3.4氣含率與曝氣量及曝氣孔直徑的關(guān)系

    氣含率是指空氣在整個(gè)裝置內(nèi)或局部區(qū)域中占有的體積,它是氣升式陶瓷膜過濾裝置的一項(xiàng)重要特征參數(shù),與氣液相的混合傳質(zhì)緊密相關(guān),反映裝置內(nèi)氣泡的分散性能與空氣的利用率,而且影響環(huán)流液體的速度和在裝置內(nèi)的停留時(shí)間[14]。

    由圖7可知,氣升管與降液管的氣含率都隨曝氣量增加而增大。當(dāng)氣泡到達(dá)液面時(shí),氣泡破裂,氣體向上沖出。氣體沖出時(shí)總會(huì)把部分拉成薄膜的液體向上拋起,被拋起的液體呈大小不一的液滴狀。液滴在上升過程中相互碰撞,其中較大的液滴上升到一定的高度會(huì)沉降下來,返回氣液分離區(qū),而較小的液滴則容易被氣相夾帶而繼續(xù)上升。氣相夾帶的小液滴進(jìn)入氣體排出管,將對后續(xù)的裝置產(chǎn)生不利影響。

    圖7 不同曝氣量下裝置整體氣相體積分?jǐn)?shù)Fig.7 Gas hold-up of equipment at different aeration rate

    由圖8可知,氣升管氣含率隨曝氣量增大而增大。隨著時(shí)間的變化,氣含率呈先增大、后減小、最后穩(wěn)定的趨勢。曝氣開始時(shí)裝置內(nèi)液體呈靜止?fàn)顟B(tài),氣泡受浮力與液相阻力的共同作用,上升速度較慢,氣泡停留時(shí)間長,所以氣含率逐漸增大。隨著曝氣的進(jìn)行,氣升管與降液管之間產(chǎn)生密度差,形成環(huán)流,氣泡上升速度加快,氣含率降低,最終趨于穩(wěn)定。

    圖8 曝氣量對氣升管氣含率的影響Fig.8 Effect of aeration rate on gas hold-up of riser

    由圖9可知,氣含率隨曝氣孔直徑減小而略微增大。這是由于氣泡的大小與曝氣孔直徑呈正相關(guān)[15]。氣泡在上升過程中受到液相曳力與浮力的共同作用,體積越小的氣泡受到的浮力越小,上升速度越慢,氣泡在氣升管中停留時(shí)間越長,故氣含率越高。但由于曝氣區(qū)域存在湍流,加劇氣泡的聚并,所以氣含率增大得并不明顯。

    圖9 曝氣孔直徑對氣升管氣含率的影響Fig.9 Effect of aerator aperture size on gas hold-up of riser

    由圖10可知,降液管氣含率隨曝氣量增大而增大,隨曝氣孔直徑增大而減小。由于液相環(huán)流的存在,氣泡不可避免地隨著液相進(jìn)入降液管。降液管中的氣泡受到的浮力的方向與液相曳力的方向相反,液相環(huán)流液速越快,氣泡受到的曳力越大;氣泡的直徑越小,受到的浮力越小,氣泡就越容易進(jìn)入降液管。

    圖10 曝氣孔直徑對降液管氣含率的影響Fig.10 Effect of aerator aperture size on gas hold-up of down-comer

    3.5環(huán)流液速與曝氣量及曝氣孔直徑的關(guān)系

    環(huán)流液速是表征氣升式陶瓷膜過濾裝置性能的關(guān)鍵特征參數(shù),它關(guān)系到裝置內(nèi)混合與傳質(zhì)的進(jìn)行[16]。液相環(huán)流形成的主要原因是氣升管與降液管之間的密度差導(dǎo)致的自然環(huán)流。一般認(rèn)為環(huán)流的推動(dòng)力即為氣升管氣含率與降液管氣含率的差值。

    由圖11可知,隨著環(huán)流推動(dòng)力(εr–εd)的增大,環(huán)流液速隨之增大。

    圖11 環(huán)流推動(dòng)力對環(huán)流液速的影響Fig.11 Effect of driving force on liquid circulation velocity

    由表2可知,隨著曝氣量的增大,升/降液管的密度差增大,進(jìn)而環(huán)流推動(dòng)力增大,導(dǎo)致環(huán)流液速增大。當(dāng)曝氣量從100 L·h?1增加到300 L·h?1時(shí),環(huán)流推動(dòng)力從9%左右增加到14%左右,環(huán)流液速也由0.23左右,增至0.29左右。當(dāng)曝氣量達(dá)到400 L·h?1時(shí),環(huán)流推動(dòng)力達(dá)到最大值16%左右,環(huán)流液速也相應(yīng)地達(dá)到最大值0.3左右。之后雖繼續(xù)增大曝氣量,但環(huán)流推動(dòng)力并未繼續(xù)增大,而是穩(wěn)定在16%左右,環(huán)流液速也隨之穩(wěn)定在0.3左右。結(jié)果表明,當(dāng)曝氣量達(dá)到400 L·h?1之后繼續(xù)增大曝氣量對增大環(huán)流液速已無明顯作用,綜合能耗的考慮,曝氣量應(yīng)為400 L·h?1為宜。

    表2 環(huán)流液速與環(huán)流推動(dòng)力Table 2 Effect of driving force on liquid circulation velocity

    由圖12可知不同直徑曝氣孔下的環(huán)流液速并無明顯的差距,可以看出環(huán)流液速與曝氣孔直徑關(guān)系不大,僅與曝氣量相關(guān)。

    3.6膜管內(nèi)湍流強(qiáng)度與曝氣量及曝氣孔直徑的關(guān)系

    隨著過濾過程的進(jìn)行,膜表面逐漸形成濃差極化,增大了膜表面阻力,增加了膜過濾裝置的維護(hù)成本,縮短了陶瓷膜的使用壽命,所以如何更好地解決濃差極化問題成為膜過濾裝置開發(fā)的關(guān)鍵[17]。研究表明,提高膜管內(nèi)流體的湍流強(qiáng)度可以有效地減輕濃差極化,防止膜面污染,減少濾餅層厚度,提高膜通量[18]。

    圖12 曝氣量對環(huán)流液速的影響Fig.12 Effect of aeration rate on liquid circulation velocity

    當(dāng)氣升式陶瓷膜過濾裝置運(yùn)行時(shí),由于微氣泡的擾動(dòng)作用,流體內(nèi)部層與層之間的平衡被打破,難以形成穩(wěn)定的層流狀態(tài),而且氣彈在上升過程中在彈頭與尾渦處因氣彈與周圍液體的速度差異引起旋渦,導(dǎo)致膜管內(nèi)氣液兩相流的湍流強(qiáng)度增加[19],有利于過濾過程的順利進(jìn)行。

    由圖13可知,膜管內(nèi)湍流強(qiáng)度隨曝氣量增大而增大。當(dāng)曝氣量達(dá)到400 L·h?1之后,膜管內(nèi)湍流強(qiáng)度的增加趨勢減緩,繼續(xù)增大曝氣量對膜管內(nèi)流體湍流強(qiáng)度的增加貢獻(xiàn)不大。這是由于曝氣量較小時(shí)膜面的流體呈現(xiàn)層流狀態(tài),氣液兩相流能夠有效地去除膜表面的濃差極化;曝氣量增至400 L·h?1后膜表面的穩(wěn)態(tài)區(qū)域消失,此時(shí)增大曝氣量只能進(jìn)一步加強(qiáng)膜表面的湍動(dòng)程度,對膜面濃差極化的抑制作用影響不大[20]。故曝氣量應(yīng)選400 L·h?1為宜。

    圖13 曝氣孔直徑對膜管內(nèi)湍流強(qiáng)度的影響Fig.13 Effect of aerator aperture size on turbulence intensity in membrane tube

    由圖14可知,當(dāng)曝氣量在100~400 L·h?1時(shí),膜管內(nèi)湍流強(qiáng)度隨曝氣孔直徑減小而增大。在低曝氣量下,較小的曝氣孔直徑產(chǎn)生的氣泡直徑小、數(shù)量多,膜管內(nèi)氣含率高,氣液兩相運(yùn)動(dòng)引發(fā)的湍流運(yùn)動(dòng)更劇烈[21]。為了更好地減輕膜表面濃差極化,防止膜面污染,應(yīng)選擇直徑0.2 mm的曝氣頭為宜。

    圖14 曝氣孔直徑對膜管內(nèi)液相壁面剪切力的影響Fig.14 Effect of aerator aperture size on wall shear stress of liquid phase

    3.7膜管液相壁面剪切力與曝氣量及曝氣孔直徑的關(guān)系

    彈狀流條件下,膜管內(nèi)的氣液兩相流由氣彈與液彈組成。液彈區(qū)域又分為尾渦區(qū)與穩(wěn)態(tài)區(qū),尾渦區(qū)的壁面剪切力波動(dòng)較大,而穩(wěn)態(tài)區(qū)的壁面剪切力接近單相流時(shí)的剪切力。膜面和氣彈之間有一層液膜,由于氣泡上升速度快于液體流速,導(dǎo)致液膜沿膜管向下流動(dòng),有利于破壞膜表面濃差極化層,消除濃差極化,減輕膜表面污染[22]。

    由圖14可知,曝氣孔直徑不變,當(dāng)曝氣量從100 L·h?1增至400 L·h?1時(shí)液彈區(qū)膜面剪切力隨曝氣量增大而增大,當(dāng)曝氣量從400 L·h?1增至600 L·h?1時(shí)膜面剪切力保持穩(wěn)定。由于曝氣量越大裝置的能耗越高,從節(jié)能與減輕膜表面污染的角度綜合考慮,曝氣量應(yīng)選400 L·h?1為宜。

    膜管內(nèi)液相壁面剪切力隨曝氣孔直徑減小而增大。結(jié)合圖15進(jìn)行分析,曝氣孔直徑減小使膜管內(nèi)氣彈聚并程度減弱,氣泡的數(shù)量增多,氣彈與液彈的交替也更加頻繁,尾渦區(qū)長度增加,穩(wěn)態(tài)區(qū)長度減小[23],膜管內(nèi)的壁面剪切力隨之增大。所以,氣升式陶瓷膜過濾裝置的曝氣頭應(yīng)選擇0.2 mm直徑為宜。

    4 結(jié) 論

    (1)裝置模型通過了網(wǎng)格無關(guān)性驗(yàn)證,而且實(shí)驗(yàn)數(shù)據(jù)與模擬數(shù)據(jù)誤差在5%~10%之間,表明本模型可較好地預(yù)測裝置內(nèi)氣液兩相流的流動(dòng)情況。模擬結(jié)果可為氣升式陶瓷膜過濾裝置曝氣頭的選型與曝氣量的選擇提供依據(jù),并對裝置將來的優(yōu)化與放大提供參考。

    圖15 不同曝氣孔直徑下的膜管內(nèi)氣泡聚并程度Fig.15 Coalescence of bubble in membrane tube under different aerator aperture size

    (2)氣升式陶瓷膜過濾裝置內(nèi)的氣含率隨曝氣量增大而增大,隨曝氣孔直徑減小而增大。氣含率的增大導(dǎo)致氣液分離區(qū)液位的升高,容易引起液滴夾帶,對后續(xù)裝置產(chǎn)生不利影響,故曝氣量應(yīng)選擇400 L·h?1左右為宜。

    (3)氣升式陶瓷膜過濾裝置內(nèi)的環(huán)流液速隨曝氣量增大而增大,當(dāng)曝氣量達(dá)到400 L·h?1時(shí)環(huán)流液速達(dá)到峰值0.3左右,之后繼續(xù)增大曝氣量并不能相應(yīng)增大環(huán)流液速。環(huán)流液速大小與曝氣孔直徑關(guān)系不大。

    (4)當(dāng)曝氣量在100~400 L·h?1時(shí),膜管內(nèi)湍流強(qiáng)度與膜管內(nèi)液相剪切力都隨曝氣孔直徑減小而增大,隨曝氣量增大而增大。曝氣量大于400 L·h?1后,繼續(xù)增大曝氣量雖然可以增大膜面的湍流強(qiáng)度,但是對膜面剪切力的影響較小,而且對膜面濃差極化的抑制作用影響不大。故曝氣量應(yīng)選擇400 L·h?1,曝氣孔直徑應(yīng)選擇0.2 mm為宜。

    符號說明

    F——?jiǎng)恿糠匠淘错?xiàng),N·m?3

    g——重力加速度,m·s?2

    I——湍流強(qiáng)度,%

    k——湍流動(dòng)能,N·m

    p——壓強(qiáng),Pa

    Q——曝氣量,L·h?1

    r——直徑,mm

    t——時(shí)間,s

    U——表觀速度,m·s?1

    u——流體速度,m·s?1

    v——速率,m·s?1

    α——體積分?jǐn)?shù),%

    ε——?dú)夂剩?

    μeff——有效黏度,m2·s?1

    ρ——密度,kg·m?3

    σ——表面張力系數(shù),N·m?1

    τ——剪切力,Pa

    下角標(biāo)

    a——曝氣頭

    d——降液管

    G——?dú)怏w

    i,j——空間坐標(biāo)

    L——液體

    r——?dú)馍?/p>

    VOL——表面張力項(xiàng)

    References

    [1]CABASSUD C, LABORIE S, LAINE J M. How slug flow can improve ultrafiltration flux in organic hollow fibres [J]. Journal of Membrane Science, 1997, 128(1): 93-101.

    [2]張文暉, 李鑫鋼. 氣升式內(nèi)環(huán)流反應(yīng)器內(nèi)局部氣含率徑向分布 [J].化工學(xué)報(bào), 2010, 61(5): 1118-1122. ZHANG W H, LI X G. Local gas holdup profiles in an internal-loop airlift reactor [J]. CIESC Journal, 2010, 61(5): 1118-1122.

    [3]ROY S, DHOTRE M T, JOSHI J B. CFD simulation of flow and axial dispersion in external loop airlift reactor [J]. Chemical Engineering Research and Design, 2006, 84(8): 677-690.

    [4]賀蕭, 車德福. 垂直及傾斜上升管內(nèi)氣液兩相彈狀流壁面剪切力的模擬 [J]. 化工學(xué)報(bào), 2008, 59(6): 1391-1395. HE X, CHE D F. CFD simulation of wall shear stress in gas-liquid vertical and inclined upward slug flow [J]. Journal of Chemical Industry and Engineering (China), 2008, 59(6): 1391-1395.

    [5]VAN SINT A M, DEEN N G, KUIPERS J. Numerical simulation of gas bubbles behavior using a three dimensional volume of fluid method [J]. Chemical Engineering Science, 2005, 60(11): 2999-3011.

    [6]ABID A, VISHNU P. CFD simulations for continuous flow of bubbles through gas-liquid columns: application of VOF method [J]. Chemical Product and Process Modeling, 2007, 9(1): 1-19.

    [7]LIN J, HAN M H, WANG T F, et al. Influence of the gas distributor on the local hydrodynamic behavior of an external loop airlift reactor [J]. Chemical Engineering Journal, 2004, 102(1): 51-59.

    [8]AURO A S, SUSHANTA K M. Effect of dynamic contact angle in a volume of fluid (VOF) model for a microfluidic capillary flow [J]. Journal of Colloid and Interface Science, 2009, 339: 461-480.

    [9]谷芳, 劉春江, 袁希鋼, 等. 傾斜波紋板上液膜流動(dòng)的CFD模擬[J]. 化工學(xué)報(bào), 2005, 56(3): 462-467. GU F, LIU C J, YUAN X G, et al. CFD simulation of liquid film flow on inclined wavy plates surface [J]. Journal of Chemical Industry and Engineering (China), 2005, 56(3): 462-467.

    [10]BRACKBILL J U, KOTHE D B, ZEMACH C A. Continuum method for modeling surface tension [J]. Computational Physics, 1992, 100:335-354.

    [11]MEHENDALE S S, JACOBI A M, AHAH R K. Fluid flow and heat transfer at micro and meso-scales with application to heat exchanger design [J]. Applied Mechanics Reviews, 2000, 53(7): 175-193.

    [12]YOUNGBAE H, NAOKI S, NOBUHIDE K. Measurement of liquid film thickness in a micro parallel channel with interferometer and laser focus displacement meter [J]. International Journal of Multiphase Flow, 2011, 37: 36-45.

    [13]TAHA T, CUI Z F. CFD modelling of gas-sparged ultrafiltration in tubular membrane [J]. Journal of Membrane Science, 2002, 210(1): 13-27.

    [14]CHAFAA B, CATHERINE X, JEAN P R. Effect of gas liquid separator and liquid height on the global hydrodynamic parameters of an external loop airlift contactor [J]. Chemical Engineering Journal, 1997, 66(2): 91-95.

    [15]宋慶唐, 金家琪, 吳桂英, 等. 二維鼓泡床內(nèi)氣泡尺寸分布的實(shí)驗(yàn)與CFD模擬 [J]. 化工學(xué)報(bào), 2008, 59(2): 336-340. SONG Q T, JIN J Q, WU G Y, et al. Experiment and CFD simulation of bubble size distribution in 2D gas-liquid bubble column[J]. Journal of Chemical Industry and Engineering (China), 2008, 59(2): 336-340.

    [16]RENZO D F. Liquid circulation rates in two and three-phase external airlift reactors [J]. Chemical Engineering Journal, 2005, 109(1/2/3): 49-55.

    [17]RICE G, BARBER A, STEVENS G, et al. Fouling of NF membranes by dairy ultrafiltration permeates [J]. Journal of Membrane Science, 2009, 330(1/2): 117-126.

    [18]MERCIER B, MARANGES M, LAFFORGUE M, et al. Hydrodynamics of slug flow applied to cross-flow filtration in narrow tubes [J]. AIChE Journal, 2000, 46: 476-488.

    [19]WIBISONO Y, CORNELISSEN E R, KEMPERMAN A J B, et al. Two-phase flow in membrane processes: a technology with a future [J]. Journal of Membrane Science, 2014, 453: 566-602.

    [20]LI Q Y, CUI Z F, PEPPER D S. Effect of bubble size and frequency on the permeate flux of gas sparged ultrafiltration with tubular membranes [J]. Chemical Engineering Journal, 1997, 67(1): 71-75.

    [21]MEI H W, XU H, JING W H, et al. Application of airlift ceramic ultrafiltration membrane ozonation reactor in the degradation of humic acids [J]. Desalination and Water Treatment, 2015, 56: 285-294.

    [22]石風(fēng)強(qiáng), 景文珩, 邢衛(wèi)紅. 陶瓷膜氣升式納濾裝置過濾桿菌肽水溶液的研究 [J].高?;瘜W(xué)工程學(xué)報(bào), 2012, 26(2):189-195. SHI F Q, JING W H, XING W H. Study on filtration of bacitracin aqueous by ceramic membrane airlift nanofiltration equipment [J]. Journal of Chemical Engineering of Chinese Universities, 2012, 26(2): 189-195.

    [23]張峰, 景文珩, 邢衛(wèi)紅. 彈狀流對陶瓷膜超濾葡聚糖水溶液的影響 [J]. 化工學(xué)報(bào), 2009, 60(11): 2792-2797. ZHANG F, JING W H, XING W H. Effect of slug flow on dextran ultrafiltration in ceramic membrane process [J]. CIESC Journal, 2009, 60(11): 2792-2797.

    Simulation of gas-liquid two-phase flow for airlift ceramic membrane filtration process

    LIN Jin, SHEN Hao, JING Wenheng
    (State Key Laboratory of Materials-Oriented Chemical Engineering, Nanjing Tech University, Nanjing 210009, Jiangsu, China)

    Abstract:The influences of aerator aperture size and aeration rate on the gas hold-up, liquid circulation velocity, wall shear stress and turbulence intensity in a novel airlift ceramic membrane filtration equipment containing a 19-channel ceramic membrane was investigated by the CFD simulation. Further, the VOF model was adopted to simulate the flow state of the gas-liquid two-phase fluid and the mean error between experiment date and simulation date was 5%—10%. The results demonstrated that the gas hold-up increased with increasing aeration rate and decreasing aerator aperture size. The liquid circulation velocity, wall shear stress and turbulence intensity increased with increasing aeration rate, and the increasing trend became slow after the aeration rate reached 400 L·h?1. For the three aperture aerator with different size, the comparison of experiment and simulation results have indicated that the liquid circulation velocity was only related to the aeration rate, and had little to do with aerator aperture size. Decreasing of the aerator aperture size was conductive to the increasing of the wall shear stress and the filtration process.

    Key words:airlift; ceramic membrane; computational fluid dynamics; VOF; two-phase flow; simulation

    中圖分類號:TQ 052

    文獻(xiàn)標(biāo)志碼:A

    文章編號:0438—1157(2016)06—2246—09

    DOI:10.11949/j.issn.0438-1157.20151560

    基金項(xiàng)目:國家自然科學(xué)基金項(xiàng)目(21176116);江蘇省高校自然科學(xué)基金重大項(xiàng)目(15KJA530001);江蘇高校優(yōu)勢學(xué)科建設(shè)工程資助項(xiàng)目。

    Corresponding author:Prof. JING Wenheng, jingwh@njtech.edu.cn

    亚洲熟妇中文字幕五十中出| 99精品在免费线老司机午夜| 国产亚洲av高清不卡| 日本五十路高清| 19禁男女啪啪无遮挡网站| 午夜福利在线观看吧| 两人在一起打扑克的视频| 99热精品在线国产| 国产精品99久久久久久久久| 久久香蕉精品热| av在线蜜桃| av中文乱码字幕在线| 国产一区二区在线av高清观看| 18禁裸乳无遮挡免费网站照片| 欧美日韩亚洲国产一区二区在线观看| 夜夜躁狠狠躁天天躁| 熟女人妻精品中文字幕| 国产野战对白在线观看| 18禁观看日本| 国产精品一区二区三区四区免费观看 | 久久午夜综合久久蜜桃| 亚洲国产中文字幕在线视频| 国产av一区在线观看免费| 免费在线观看日本一区| 亚洲国产欧美网| 老熟妇乱子伦视频在线观看| 在线观看66精品国产| 国产亚洲精品久久久久久毛片| 久久精品人妻少妇| 日韩三级视频一区二区三区| 18禁黄网站禁片午夜丰满| 亚洲人成网站高清观看| 国产高潮美女av| 欧美黄色片欧美黄色片| 毛片女人毛片| 久久这里只有精品19| 精品欧美国产一区二区三| 九九在线视频观看精品| 国产在线精品亚洲第一网站| 国产欧美日韩一区二区三| 女人高潮潮喷娇喘18禁视频| 久久欧美精品欧美久久欧美| 搡老岳熟女国产| 久久精品夜夜夜夜夜久久蜜豆| 宅男免费午夜| АⅤ资源中文在线天堂| 久久中文字幕人妻熟女| 在线观看午夜福利视频| 给我免费播放毛片高清在线观看| 99久久精品国产亚洲精品| 丰满人妻熟妇乱又伦精品不卡| 成人精品一区二区免费| 超碰成人久久| 国产精品久久久久久精品电影| 男女视频在线观看网站免费| 久久热在线av| 成人国产一区最新在线观看| 91av网站免费观看| 国产久久久一区二区三区| 一个人看的www免费观看视频| 国产精品 欧美亚洲| 好看av亚洲va欧美ⅴa在| 午夜免费成人在线视频| 欧美性猛交╳xxx乱大交人| 在线免费观看不下载黄p国产 | 草草在线视频免费看| 久久中文字幕一级| 99热这里只有精品一区 | 一个人免费在线观看的高清视频| 非洲黑人性xxxx精品又粗又长| 在线国产一区二区在线| 99re在线观看精品视频| 亚洲精品久久国产高清桃花| 色尼玛亚洲综合影院| 婷婷精品国产亚洲av| 丰满的人妻完整版| 男插女下体视频免费在线播放| 中文字幕最新亚洲高清| 亚洲成人久久性| 亚洲精品粉嫩美女一区| 国产成人一区二区三区免费视频网站| 在线十欧美十亚洲十日本专区| 丁香欧美五月| 久久久久九九精品影院| 亚洲在线观看片| 999精品在线视频| 床上黄色一级片| 99久久国产精品久久久| 中文字幕高清在线视频| 美女免费视频网站| 国模一区二区三区四区视频 | 国产亚洲欧美在线一区二区| 老司机午夜福利在线观看视频| 亚洲专区字幕在线| 五月伊人婷婷丁香| 两性夫妻黄色片| 最近最新中文字幕大全电影3| 午夜两性在线视频| 久久精品综合一区二区三区| 91麻豆av在线| 亚洲国产高清在线一区二区三| 日韩欧美三级三区| 国产成人一区二区三区免费视频网站| 日本三级黄在线观看| 窝窝影院91人妻| 亚洲精品在线美女| 午夜免费观看网址| 亚洲av第一区精品v没综合| 18禁黄网站禁片午夜丰满| 窝窝影院91人妻| 每晚都被弄得嗷嗷叫到高潮| 国产精品一区二区三区四区久久| 99国产精品一区二区蜜桃av| 国产麻豆成人av免费视频| 中文字幕久久专区| 黄色成人免费大全| 午夜福利欧美成人| 在线永久观看黄色视频| 亚洲国产精品成人综合色| 最新在线观看一区二区三区| 成人欧美大片| 午夜福利高清视频| 天堂动漫精品| 毛片女人毛片| 欧美一级a爱片免费观看看| 亚洲欧美精品综合一区二区三区| 久久国产乱子伦精品免费另类| 色吧在线观看| 18美女黄网站色大片免费观看| 在线观看免费午夜福利视频| 法律面前人人平等表现在哪些方面| 免费在线观看日本一区| 搞女人的毛片| 国产午夜福利久久久久久| 国产欧美日韩精品亚洲av| 国产精华一区二区三区| 99热精品在线国产| 美女午夜性视频免费| 国产精品久久视频播放| 人妻丰满熟妇av一区二区三区| 一区二区三区高清视频在线| 香蕉av资源在线| 99精品欧美一区二区三区四区| 超碰成人久久| 国产极品精品免费视频能看的| 18禁观看日本| 99re在线观看精品视频| 国产成人福利小说| 亚洲精品久久国产高清桃花| 色综合欧美亚洲国产小说| 一区二区三区国产精品乱码| 一区二区三区激情视频| 男女下面进入的视频免费午夜| 国产精品自产拍在线观看55亚洲| 亚洲欧美日韩高清专用| 老司机午夜福利在线观看视频| 人人妻人人看人人澡| 日本黄大片高清| 99热只有精品国产| 老鸭窝网址在线观看| 亚洲一区高清亚洲精品| 色综合欧美亚洲国产小说| 亚洲欧美一区二区三区黑人| 亚洲成av人片免费观看| 最新在线观看一区二区三区| 最新美女视频免费是黄的| 好看av亚洲va欧美ⅴa在| av在线蜜桃| 国产私拍福利视频在线观看| 国产成人精品无人区| 可以在线观看的亚洲视频| 久久国产精品人妻蜜桃| 一个人看的www免费观看视频| 亚洲美女黄片视频| 天天添夜夜摸| 女警被强在线播放| 九色国产91popny在线| 国产精品一区二区精品视频观看| a级毛片a级免费在线| 色av中文字幕| 老汉色∧v一级毛片| 精品国产乱码久久久久久男人| av黄色大香蕉| 搡老岳熟女国产| 毛片女人毛片| 国产av麻豆久久久久久久| 国产黄片美女视频| 一a级毛片在线观看| 中国美女看黄片| 国产97色在线日韩免费| 真人做人爱边吃奶动态| 99久久精品热视频| 亚洲一区高清亚洲精品| 国产一区二区三区视频了| 嫩草影院精品99| 亚洲欧洲精品一区二区精品久久久| 免费看光身美女| 一进一出抽搐gif免费好疼| 日本熟妇午夜| 搡老岳熟女国产| 午夜福利免费观看在线| 一级作爱视频免费观看| 久久久精品大字幕| 午夜视频精品福利| 叶爱在线成人免费视频播放| 美女扒开内裤让男人捅视频| 亚洲中文字幕日韩| 国产一区二区激情短视频| 性欧美人与动物交配| 麻豆成人av在线观看| 亚洲av免费在线观看| 亚洲国产色片| 欧美乱色亚洲激情| 亚洲美女视频黄频| 18禁裸乳无遮挡免费网站照片| 国产探花在线观看一区二区| 97碰自拍视频| а√天堂www在线а√下载| 欧美另类亚洲清纯唯美| 757午夜福利合集在线观看| 国产又色又爽无遮挡免费看| 一边摸一边抽搐一进一小说| 国产成人av激情在线播放| 亚洲国产高清在线一区二区三| 免费电影在线观看免费观看| 日本免费a在线| 国产日本99.免费观看| 精品久久久久久久人妻蜜臀av| 午夜福利在线观看免费完整高清在 | 午夜福利在线在线| 五月玫瑰六月丁香| 国产精华一区二区三区| 天堂网av新在线| 亚洲熟女毛片儿| 中文字幕人妻丝袜一区二区| 久久午夜亚洲精品久久| 亚洲欧洲精品一区二区精品久久久| 亚洲精品在线观看二区| 亚洲欧美日韩高清专用| 精品不卡国产一区二区三区| 亚洲欧美日韩卡通动漫| 一本综合久久免费| 亚洲无线在线观看| 夜夜躁狠狠躁天天躁| 国产精品香港三级国产av潘金莲| 国产97色在线日韩免费| 一区二区三区高清视频在线| 久久性视频一级片| 欧美日韩精品网址| 性色av乱码一区二区三区2| 国产探花在线观看一区二区| 久久久久久久久久黄片| 精品久久久久久久久久久久久| 最好的美女福利视频网| 熟女人妻精品中文字幕| 亚洲熟女毛片儿| 亚洲av第一区精品v没综合| 国产精品野战在线观看| 欧美黑人欧美精品刺激| 国产 一区 欧美 日韩| 在线视频色国产色| 在线看三级毛片| 国内精品美女久久久久久| 日韩av在线大香蕉| 老司机深夜福利视频在线观看| 亚洲国产欧美人成| 脱女人内裤的视频| 99热只有精品国产| 最近最新中文字幕大全电影3| 18禁黄网站禁片免费观看直播| 成人av在线播放网站| 国产视频一区二区在线看| 91在线精品国自产拍蜜月 | 又粗又爽又猛毛片免费看| 黄片大片在线免费观看| 国产精品98久久久久久宅男小说| 亚洲aⅴ乱码一区二区在线播放| 日韩欧美在线乱码| 国产亚洲欧美在线一区二区| 久久久久久久久免费视频了| 他把我摸到了高潮在线观看| 久久伊人香网站| 亚洲国产精品sss在线观看| 国产美女午夜福利| 亚洲专区中文字幕在线| 好男人在线观看高清免费视频| 变态另类成人亚洲欧美熟女| 搡老妇女老女人老熟妇| 久久久久免费精品人妻一区二区| 亚洲狠狠婷婷综合久久图片| 99国产精品99久久久久| 久久精品综合一区二区三区| 国产精品永久免费网站| 欧美日韩黄片免| 亚洲av五月六月丁香网| 亚洲 欧美 日韩 在线 免费| 亚洲精品美女久久久久99蜜臀| 欧美乱妇无乱码| 精品国产乱子伦一区二区三区| 动漫黄色视频在线观看| 亚洲无线在线观看| 日韩欧美国产在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品久久国产高清桃花| 一级毛片女人18水好多| 91在线观看av| 亚洲无线观看免费| 色在线成人网| 麻豆av在线久日| 久久人妻av系列| 1000部很黄的大片| 黄色视频,在线免费观看| 色av中文字幕| 亚洲国产高清在线一区二区三| 90打野战视频偷拍视频| av片东京热男人的天堂| 欧美不卡视频在线免费观看| 亚洲国产欧洲综合997久久,| 国产精品 欧美亚洲| 日韩欧美国产一区二区入口| 99国产精品一区二区蜜桃av| 窝窝影院91人妻| 国产精品av视频在线免费观看| 欧美日韩一级在线毛片| 欧美乱色亚洲激情| 免费av毛片视频| 两性夫妻黄色片| 精品一区二区三区av网在线观看| 欧美+亚洲+日韩+国产| 日本免费a在线| 天堂影院成人在线观看| 久久精品人妻少妇| 校园春色视频在线观看| 国产亚洲av高清不卡| 久久久国产精品麻豆| 99久国产av精品| 亚洲中文字幕一区二区三区有码在线看 | 偷拍熟女少妇极品色| www.熟女人妻精品国产| 国产精品综合久久久久久久免费| 国产成人精品久久二区二区91| 国产av麻豆久久久久久久| 美女大奶头视频| 高清毛片免费观看视频网站| 久久精品国产综合久久久| 久久久久久久久久黄片| 人妻丰满熟妇av一区二区三区| 亚洲精品色激情综合| 久久午夜综合久久蜜桃| 亚洲18禁久久av| 中文在线观看免费www的网站| 热99在线观看视频| 天堂影院成人在线观看| 88av欧美| 99精品在免费线老司机午夜| 亚洲无线在线观看| 亚洲无线观看免费| 狠狠狠狠99中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 性色av乱码一区二区三区2| 国产精品99久久久久久久久| 久久久久久久精品吃奶| 色老头精品视频在线观看| 国产成人精品久久二区二区91| 又黄又爽又免费观看的视频| 免费人成视频x8x8入口观看| 久久国产精品影院| 毛片女人毛片| 国内少妇人妻偷人精品xxx网站 | 精华霜和精华液先用哪个| 成人鲁丝片一二三区免费| 日本与韩国留学比较| 亚洲成a人片在线一区二区| 人人妻,人人澡人人爽秒播| aaaaa片日本免费| 国产精品久久久久久亚洲av鲁大| 99久久精品国产亚洲精品| 欧美日本亚洲视频在线播放| 91麻豆av在线| 国产精品久久久久久精品电影| 成人18禁在线播放| 免费电影在线观看免费观看| 中文字幕久久专区| 久久精品综合一区二区三区| 欧美日韩国产亚洲二区| 国产成人精品无人区| 成年版毛片免费区| 这个男人来自地球电影免费观看| 亚洲电影在线观看av| tocl精华| 欧美又色又爽又黄视频| 两人在一起打扑克的视频| 2021天堂中文幕一二区在线观| 成人鲁丝片一二三区免费| 搡老岳熟女国产| 亚洲avbb在线观看| 色综合亚洲欧美另类图片| 最新美女视频免费是黄的| 日本精品一区二区三区蜜桃| 一区二区三区国产精品乱码| 最近视频中文字幕2019在线8| 久久精品国产综合久久久| 在线观看美女被高潮喷水网站 | 亚洲avbb在线观看| 国产精品免费一区二区三区在线| 99热精品在线国产| 国产精品精品国产色婷婷| 国产伦在线观看视频一区| 狠狠狠狠99中文字幕| 国产黄片美女视频| 久久九九热精品免费| 成人av一区二区三区在线看| 99热6这里只有精品| 无遮挡黄片免费观看| 亚洲av电影不卡..在线观看| 1024香蕉在线观看| 日本在线视频免费播放| 18禁黄网站禁片免费观看直播| 后天国语完整版免费观看| 国产伦在线观看视频一区| 天堂动漫精品| 国产成人av激情在线播放| 久久香蕉国产精品| 好看av亚洲va欧美ⅴa在| 岛国在线免费视频观看| 亚洲成人免费电影在线观看| 日韩有码中文字幕| 亚洲国产精品999在线| 国产激情久久老熟女| 亚洲激情在线av| 国产成人av激情在线播放| 国产亚洲欧美在线一区二区| 亚洲专区字幕在线| 久久精品91无色码中文字幕| 欧美三级亚洲精品| 亚洲人成伊人成综合网2020| 男人舔女人下体高潮全视频| 久久久久久九九精品二区国产| 欧美乱妇无乱码| 成人国产一区最新在线观看| 91麻豆av在线| 在线免费观看的www视频| 国产欧美日韩精品一区二区| 俄罗斯特黄特色一大片| 国产精品一区二区精品视频观看| 精品欧美国产一区二区三| 此物有八面人人有两片| 午夜福利成人在线免费观看| 99热只有精品国产| 亚洲av片天天在线观看| 91久久精品国产一区二区成人 | 欧美又色又爽又黄视频| 亚洲av电影不卡..在线观看| 88av欧美| 久久久久免费精品人妻一区二区| 亚洲国产欧美一区二区综合| 国内精品美女久久久久久| 精品一区二区三区四区五区乱码| 在线观看一区二区三区| 国内精品久久久久精免费| 一个人免费在线观看的高清视频| 老熟妇乱子伦视频在线观看| 日本一二三区视频观看| 午夜福利在线在线| 日本三级黄在线观看| 亚洲av成人av| 特级一级黄色大片| 嫁个100分男人电影在线观看| 日韩欧美 国产精品| 欧美黄色淫秽网站| 亚洲激情在线av| 国产1区2区3区精品| 日本a在线网址| 在线a可以看的网站| 成人性生交大片免费视频hd| 久久国产乱子伦精品免费另类| 日本三级黄在线观看| 久久精品影院6| 成人永久免费在线观看视频| 久久精品91无色码中文字幕| 男人舔女人的私密视频| 成人三级做爰电影| 久久久久久久久久黄片| 搡老熟女国产l中国老女人| 亚洲五月天丁香| 成人av一区二区三区在线看| 国产v大片淫在线免费观看| 色吧在线观看| 韩国av一区二区三区四区| 嫩草影院入口| 久久精品人妻少妇| 蜜桃久久精品国产亚洲av| 亚洲av成人av| 亚洲最大成人中文| 国产亚洲av嫩草精品影院| 国产三级中文精品| 欧美国产日韩亚洲一区| 国产午夜精品久久久久久| 国产精品亚洲美女久久久| a级毛片在线看网站| 12—13女人毛片做爰片一| 狂野欧美激情性xxxx| 国产91精品成人一区二区三区| 色综合婷婷激情| 男女视频在线观看网站免费| 18禁观看日本| 国产成人啪精品午夜网站| 亚洲熟妇熟女久久| 日韩精品中文字幕看吧| 后天国语完整版免费观看| 亚洲欧美日韩无卡精品| 最近最新中文字幕大全电影3| 精品久久久久久久毛片微露脸| 国语自产精品视频在线第100页| 一级作爱视频免费观看| 亚洲国产看品久久| 人妻丰满熟妇av一区二区三区| 观看美女的网站| 亚洲电影在线观看av| 最新在线观看一区二区三区| 在线免费观看的www视频| 亚洲无线观看免费| 免费观看的影片在线观看| 男女午夜视频在线观看| 久久热在线av| 国产伦人伦偷精品视频| 999精品在线视频| 亚洲第一欧美日韩一区二区三区| 美女高潮喷水抽搐中文字幕| 欧美绝顶高潮抽搐喷水| 国产精品电影一区二区三区| 免费大片18禁| 亚洲中文字幕一区二区三区有码在线看 | 老司机深夜福利视频在线观看| 搡老妇女老女人老熟妇| 老司机在亚洲福利影院| 国产美女午夜福利| 亚洲美女视频黄频| 黑人欧美特级aaaaaa片| 哪里可以看免费的av片| 国产亚洲精品综合一区在线观看| www.www免费av| 男女午夜视频在线观看| 黑人欧美特级aaaaaa片| 一夜夜www| 两个人的视频大全免费| 亚洲男人的天堂狠狠| 久久精品综合一区二区三区| 变态另类成人亚洲欧美熟女| 99久久精品国产亚洲精品| 免费电影在线观看免费观看| 特级一级黄色大片| 99在线视频只有这里精品首页| 亚洲国产看品久久| 麻豆一二三区av精品| 天天躁日日操中文字幕| 嫩草影视91久久| 天堂影院成人在线观看| 在线观看美女被高潮喷水网站 | 国产伦精品一区二区三区四那| 啪啪无遮挡十八禁网站| 欧美又色又爽又黄视频| 真实男女啪啪啪动态图| 精品免费久久久久久久清纯| 久久久久久人人人人人| 亚洲第一电影网av| 天天添夜夜摸| 亚洲中文字幕一区二区三区有码在线看 | 日韩欧美 国产精品| 午夜成年电影在线免费观看| 亚洲性夜色夜夜综合| 精品一区二区三区视频在线观看免费| 亚洲欧美精品综合久久99| 国产三级黄色录像| 精品不卡国产一区二区三区| 久久久国产成人精品二区| 国产av不卡久久| 国产精品1区2区在线观看.| xxxwww97欧美| 国产精品一区二区三区四区久久| 婷婷亚洲欧美| 国产日本99.免费观看| 男女床上黄色一级片免费看| 亚洲午夜精品一区,二区,三区| 在线免费观看不下载黄p国产 | 好看av亚洲va欧美ⅴa在| 久久久色成人| 男女下面进入的视频免费午夜| 色综合婷婷激情| 久久九九热精品免费| 在线观看一区二区三区| 9191精品国产免费久久| 国产精华一区二区三区| 国产不卡一卡二| 一区二区三区国产精品乱码| 午夜福利视频1000在线观看| 色噜噜av男人的天堂激情| netflix在线观看网站| 亚洲成av人片免费观看| 久久亚洲精品不卡| 欧美高清成人免费视频www| 亚洲成av人片免费观看| avwww免费| 久久久久免费精品人妻一区二区| 国产欧美日韩一区二区精品| 黄色女人牲交| 超碰成人久久| 美女 人体艺术 gogo| 国产精品免费一区二区三区在线| 欧美日韩瑟瑟在线播放| 老司机午夜十八禁免费视频| www.999成人在线观看| 免费av不卡在线播放| 亚洲精品456在线播放app | 一个人看的www免费观看视频|