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

    提升管出口SVQS自然旋風(fēng)長度的數(shù)值模擬

    2016-09-18 09:56:40程兆龍鄂承林盧春喜
    化工學(xué)報 2016年8期
    關(guān)鍵詞:切向速度汽提氣速

    程兆龍,鄂承林,盧春喜

    (中國石油大學(xué)(北京)重質(zhì)油國家重點(diǎn)實(shí)驗室,北京 102249)

    提升管出口SVQS自然旋風(fēng)長度的數(shù)值模擬

    程兆龍,鄂承林,盧春喜

    (中國石油大學(xué)(北京)重質(zhì)油國家重點(diǎn)實(shí)驗室,北京 102249)

    采用RNG k-ε湍流模型對一套?600 mm×4150 mm的SVQS旋流快分裝置的氣相流場及自然旋風(fēng)長度進(jìn)行了模擬研究。結(jié)果表明,SVQS旋風(fēng)尾渦的截面位置與隔流筒下端面最大切向速度衰減達(dá)88%時的截面位置相吻合,并據(jù)此定義了SVQS的自然旋風(fēng)長度?;赟VQS旋流快分切向速度的分布規(guī)律發(fā)現(xiàn),SVQS的自然旋風(fēng)長度隨噴口氣速的增加和汽提氣速的減小而逐漸增大。參照旋風(fēng)分離器自然旋風(fēng)長度的計算方法,基于SVQS旋流快分的結(jié)構(gòu)特點(diǎn)及模擬計算結(jié)果,提出了提升管出口SVQS旋流快分自然旋風(fēng)長度的計算關(guān)系式。

    提升管;分離;計算流體力學(xué)

    引 言

    催化裂化提升管出口旋流快分是實(shí)現(xiàn)油劑高效分離及油氣快速引出的重要設(shè)備,其原理是利用提升管出口油氣的高速旋轉(zhuǎn)流動所產(chǎn)生的離心力場進(jìn)行氣固分離的。針對現(xiàn)有旋流快分噴口處存在部分上行短路流而使分離效率降低的問題,盧春喜等[1]開發(fā)了一種帶隔流筒和隔流蓋板的提升管出口旋流快分系統(tǒng)(super vortex quick separation system,SVQS)。孫鳳俠等[2-4]通過研究發(fā)現(xiàn),隔流筒的直徑存在著一個最佳適中值,同時隔流筒的長度應(yīng)以稍大為宜。胡艷華等[5]則將原有的圓筒狀隔流筒改為折邊隔流筒和圓錐形隔流筒,發(fā)現(xiàn)圓錐形隔流筒的氣固分離效果更佳。除了隔流筒的結(jié)構(gòu)和尺寸外,SVQS旋流快分的自然旋風(fēng)長度也是影響其氣固分離效率及油氣停留時間分布的重要參數(shù)。在旋流快分分離效率模型建立的過程中,其高度采用的是自然旋風(fēng)長度而不是物理高度。不僅如此,如果在工業(yè)設(shè)計SVQS旋流快分時能夠知曉SVQS的自然旋風(fēng)長度,這對確定SVQS封閉罩的安裝高度以降低整套工業(yè)裝置的高度和設(shè)備投資費(fèi)用都將大有裨益。因此,自然旋風(fēng)長度的確定對SVQS旋流快分的優(yōu)化設(shè)計具有重要意義。

    目前關(guān)于自然旋風(fēng)長度的研究還主要集中在旋風(fēng)分離器上,而對提升管出口旋流快分自然旋風(fēng)長度的研究卻還未見報道。在旋風(fēng)分離器自然旋風(fēng)長度的研究方面,Alexander[6]提出旋風(fēng)分離器自然旋風(fēng)長度的概念,并通過實(shí)驗發(fā)現(xiàn)長筒型旋風(fēng)分離器內(nèi)的氣相流場存在著一個穩(wěn)定的旋轉(zhuǎn)流動轉(zhuǎn)折點(diǎn),此轉(zhuǎn)折點(diǎn)(尾渦位置)與升氣管下口截面之間的軸向高度即為旋風(fēng)分離器的自然旋風(fēng)長度。在Leith等[7]以及 Dietz[8]的計算旋風(fēng)分離器分離效率的公式中,也都應(yīng)用了Alexander提出的自然旋風(fēng)長度的計算公式。Hoffmann等[9]根據(jù)實(shí)驗,分析了入口顆粒速度和顆粒濃度對旋風(fēng)分離器自然旋風(fēng)長度的影響規(guī)律。Bryant等[10]認(rèn)為影響旋風(fēng)分離器自然旋風(fēng)長度的重要因素為旋風(fēng)分離器的入口截面比和升氣管的下口直徑比。姬忠禮等[11]通過實(shí)驗發(fā)現(xiàn),在旋風(fēng)分離器的尾渦處存在著一個灰環(huán),灰環(huán)截面處的切向速度分布較為平坦,其所在截面的最大切向速度約為升氣管下口截面的12%,并將這一截面與升氣管下口截面之間的軸向高度定義為旋風(fēng)分離器的自然旋風(fēng)長度。杜德喜等[12]通過數(shù)值計算發(fā)現(xiàn),旋風(fēng)分離器的自然旋風(fēng)長度是多個自然轉(zhuǎn)折旋風(fēng)串聯(lián)的結(jié)果。高翠芝等[13]認(rèn)為當(dāng)旋風(fēng)分離器內(nèi)外旋流之間的能量傳遞達(dá)到穩(wěn)定狀態(tài)時,旋轉(zhuǎn)氣流才到達(dá)漩渦的尾端位置,并據(jù)此推導(dǎo)出了旋風(fēng)長度的計算公式。魏耀東等[14]采用激光多普勒測速技術(shù)對旋風(fēng)分離器內(nèi)的氣相流場進(jìn)行了測定,并基于測量結(jié)果定義了升氣管下口到準(zhǔn)自由渦消失位置處的距離為旋風(fēng)分離器的自然旋風(fēng)長度。錢付平等[15]則基于曲面響應(yīng)法分析了旋風(fēng)分離器自然旋風(fēng)長度的影響因素,并得到了旋風(fēng)長度的預(yù)測模型。

    鑒于目前還沒有對提升管出口旋流快分自然旋風(fēng)長度方面的相關(guān)研究,本文通過數(shù)值計算的方法,在旋風(fēng)分離器自然旋風(fēng)長度研究的基礎(chǔ)上,對一套?600 mm×4150 mm的提升管出口SVQS旋流快分在不同噴口氣速及汽提氣速下的自然旋風(fēng)長度進(jìn)行了模擬研究,以期為SVQS旋流快分的優(yōu)化設(shè)計提供基礎(chǔ)數(shù)據(jù)。

    1 數(shù)值計算方法及其驗證

    本文采用可以良好處理流線彎曲程度較大流動情況的RNG k-ε湍流模型來模擬SVQS旋流快分系統(tǒng)內(nèi)的氣相流場,其通式見式(1),控制方程見表1[16-17]。

    所有方程均采用有限體積法進(jìn)行建立,壓力-速度耦合方程采用壓力耦合方程半隱法(sem i-implicit method for pressure linked equations,SIMPLE)進(jìn)行求解。動量方程、湍動能方程及湍流擴(kuò)散率方程均采用對流項二次迎風(fēng)差值格式(quadratic upw ind interpolation of convective Kinematics,QUICK)進(jìn)行離散。本文基于文獻(xiàn)[2]的實(shí)驗結(jié)果對RNG k-ε湍流模型進(jìn)行了驗證,模擬值與測量值的驗證結(jié)果及網(wǎng)格無關(guān)性檢驗詳見文獻(xiàn)[18]。

    表1 控制方程Table 1 Control equation

    圖1 SVQS實(shí)驗裝置結(jié)構(gòu)及尺寸示意圖Fig.1 Diagram of experimental apparatus and size

    2 SVQS旋流快分網(wǎng)格劃分及邊界條件

    2.1網(wǎng)格劃分

    圖1(a)、(b)分別為本文模擬計算的提升管出口SVQS旋流快分實(shí)驗裝置的3D和2D簡圖。實(shí)驗裝置主要由旋流頭、隔流筒、隔流蓋板、封閉罩及提升管等組成。裝置總高約5360 mm,封閉罩直徑為?600 mm、高度為 4150 mm,提升管直徑為 ?100 mm,隔流筒直徑為?380 mm、高度為475 mm。旋流快分頭為 3旋臂型SVQS旋流快分裝置,其噴口截面尺寸為88 mm×29 mm。由于頂部旋風(fēng)分離器對SVQS旋流快分自然旋風(fēng)長度的影響很小,所以在計算過程中忽略了旋風(fēng)分離器區(qū)域,取計算區(qū)域如圖1(c)所示。圖2為3旋臂SVQS旋流快分的網(wǎng)格結(jié)構(gòu)示意圖,SVQS旋流快分均采用六面體結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分,整個SVQS旋流快分實(shí)驗裝置共劃分了4886046個網(wǎng)格。

    2.2邊界條件

    圖2 SVQS實(shí)驗裝置網(wǎng)格結(jié)構(gòu)示意圖Fig.2 Mesh of SVQS experimental apparatus

    表2 提升管入口邊界條件Table 2 Boundary conditions of riser inlet

    (2)實(shí)驗裝置的出口設(shè)定為壓力出口,出口壓力為標(biāo)準(zhǔn)大氣壓,水力學(xué)直徑Dh為?100 mm,出口面積Ao為0.00786 m2。

    (3)實(shí)驗裝置內(nèi)的所有壁面均采用標(biāo)準(zhǔn)壁面函數(shù)進(jìn)行處理。

    3 結(jié)果分析與討論

    3.1自然旋風(fēng)長度的確定

    為了對提升管出口SVQS旋流快分的自然旋風(fēng)長度進(jìn)行預(yù)測,沿圖1(b)中的封閉罩軸向自下而上取8個截面作為計算截面,其距封閉罩下部錐形擋板上端的軸向距離分別為:0.97、1.0、1.1、1.5、1.9、2.3、2.7及3.1 m,8個計算截面的具體位置如圖3所示。

    圖4 SVQS旋流頭噴口中心處水平截面上速度分布(Vs=22 m·s-1)Fig.4 Velocity distribution on horizontal section across center of SVQS vortex spouts (Vs=22 m·s-1)

    圖5 SVQS氣體跡線(Vs=22 m·s-1)Fig.5 Stream line of SVQS (Vs=22 m·s-1)

    圖4為SVQS旋流快分噴出口中心處橫截面上的速度分布。由圖4可見,旋流快分噴出口中心處的速度分布呈中心對稱分布,氣體從旋流臂噴口噴出后,在隔流筒和隔流蓋板的約束下緊貼封閉罩內(nèi)壁向下做強(qiáng)旋流動。圖5為在旋流快分噴口氣速為22 m·s-1時,噴口氣體在封閉罩內(nèi)的旋流跡線圖。由圖5可見,進(jìn)入3旋臂型SVQS旋流快分的提升管氣體,在旋流臂的約束下,從旋流臂出口沿封閉罩內(nèi)壁切向噴出,同時沿封閉罩軸向向下做旋轉(zhuǎn)流動。向下旋轉(zhuǎn)流動的氣體并沒有在到達(dá)封閉罩下部錐形擋板后才向上折返,而是在距封閉罩下部錐形擋板上端約0.97 m時便開始向上折返流動。隔流筒下端面距該折返截面(旋渦尾端)的軸向高度2.18 m即為SVQS旋流快分在噴口氣速為22 m·s-1時的自然旋風(fēng)長度。旋風(fēng)分離器的有效分離空間為旋渦尾端的上部空間,而在旋渦尾端的下部,顆粒沿壁面到出口的傳遞是低效的。如果直筒體的高度過小使旋風(fēng)尾渦進(jìn)入到底部的錐形擋板中,則急速旋轉(zhuǎn)的氣流就會把已經(jīng)沉降到錐體段內(nèi)表面的顆粒重新卷起,并被旋轉(zhuǎn)向上的內(nèi)旋流帶出,而使SVQS旋流快分的分離效率大為降低,同時錐形擋板上隨氣流一起旋轉(zhuǎn)的顆粒由于受到離心力的作用會被強(qiáng)烈地擠壓到錐形擋板上而對汽提擋板的內(nèi)表面產(chǎn)生磨蝕。

    圖6為在噴出口氣速為22 m·s-1時,量綱1切向速度Vt/Vo在封閉罩8個軸向截面上的徑向分布曲線,圖中Vo為封閉罩與提升管環(huán)形截面上的表觀氣速。由圖6可見,切向速度沿徑向向外不斷增大并在靠近封閉罩壁面處達(dá)到最大值。在軸向方向上,隨著軸向高度的降低,最大切向速度不斷減小,當(dāng)軸向高度由z=3.1 m下降至z=1.1 m時,最大切向速度下降得較為緩慢,共下降了60%左右;繼續(xù)降低軸向高度,其值迅速減小,在下降至 z=0.97 m截面處,最大切向速度減小至不足隔流筒下口截面(z=3.1 m)處最大切向速度的12%。通過與圖5中相同噴口氣速條件下的旋流跡線圖進(jìn)行對比可以發(fā)現(xiàn),最大切向速度衰減達(dá)88%時的軸向位置恰好位于自然旋風(fēng)長度的尾部,這與文獻(xiàn)[11]在旋風(fēng)分離器中的測量結(jié)果一致。為此,本文取隔流筒下端面最大切向速度衰減達(dá)88%時的截面與隔流筒下端面間的軸向高度為SVQS旋流快分的自然旋風(fēng)長度。

    圖6 SVQS量綱1切向速度分布(Vs=22 m·s-1)Fig.6 Dimensionless tangential velocity distribution in SVQS (Vs=22 m·s-1)

    3.2噴口氣速對SVQS自然旋風(fēng)長度的影響

    圖7為不同噴口氣速下SVQS旋流快分系統(tǒng)內(nèi)量綱1切向速度在不同軸向截面上的分布曲線。由圖7可知,隨著旋流頭噴出口氣速的增加,同一截面高度上的切向速度也逐漸增大,并且隨著軸向位置的降低,切向速度的徑向梯度也在逐漸減小且分布更加平緩。當(dāng)截面的軸向高度一定時,隨著徑向位置的增大,切向速度在r/R<0.5的區(qū)域內(nèi)變化較小,在r/R>0.5的徑向區(qū)域則呈現(xiàn)先增加后降低的變化趨勢,并在邊壁附近處達(dá)到最大值。隨著噴口氣速的增加,最大切向速度衰減達(dá)88%時的截面位置逐漸降低,當(dāng)噴口氣速由 8 m·s-1增加到 22 m·s-1時,最大切向速度衰減達(dá)88%時的截面位置逐漸由z=2.36 m降低至z=0.97 m,自然旋風(fēng)長度則由0.79 m逐漸增加到2.18 m。由此可知,SVQS旋流快分的自然旋風(fēng)長度隨噴口氣速的增大而明顯增加。

    圖7 噴出口氣速對SVQS切向速度的影響Fig.7 Dimensionless tangential velocity w ith various SVQS spout velocities

    3.3汽提氣速對SVQS自然旋風(fēng)長度的影響

    在SVQS旋流快分系統(tǒng)的下部設(shè)有汽提段,以回收催化劑夾帶的油氣。為了考察汽提氣速對SVQS旋流快分自然旋風(fēng)長度的影響規(guī)律,本文在噴口氣速Vs=22 m·s-1,不同汽提氣速下模擬計算了SVQS系統(tǒng)量綱1切向速度沿軸向高度的變化規(guī)律,如圖8所示。

    由圖8可知,汽提氣的引入對SVQS旋流快分系統(tǒng)上部空間氣相流場的影響較小,而對下部空間氣相流場的影響較大,在汽提氣速Vq由0.15 m·s-1增加到0.90 m·s-1的過程中,最大切向速度衰減達(dá)88%時的截面位置從截面z=0.905 m不斷上升至截面z=1.36 m處,自然旋風(fēng)長度由2.245 m縮短為1.790 m。由此可知,汽提氣的引入縮短了SVQS旋流快分的自然旋風(fēng)長度。

    3.4SVQS自然旋風(fēng)長度的計算關(guān)系式

    SVQS旋流快分旋風(fēng)尾渦的軸向位置在不同操作條件下具有較大的變化范圍,因此采用實(shí)驗的方法較難測量出不同條件下SVQS旋流尾渦的軸向位置,而采用模擬計算的方法卻可以相對容易地預(yù)測不同操作條件下SVQS旋流尾渦的截面位置,進(jìn)而獲得SVQS旋流快分的自然旋風(fēng)長度。本文通過參照旋風(fēng)分離器自然旋風(fēng)長度的計算方法,基于SVQS旋流快分的結(jié)構(gòu)特點(diǎn)及計算模擬結(jié)果,確定了提升管出口SVQS旋流快分自然旋風(fēng)長度的計算關(guān)系式。

    圖8 不同汽提氣速的切向速度分布(Vs=22 m·s-1)Fig.8 Dimensionless velocity distribution w ith various stripping velocities (Vs=22 m·s-1)

    圖9 SVQS旋流快分的自然旋風(fēng)長度Fig.9 Natural vortex length of SVQS

    圖10 優(yōu)化高度后的SVQS量綱1切向速度分布及流線圖(Vs=24 m·s-1)Fig.10 Dimensionless tangential velocity distribution and stream line of optim ized SVQS (Vs=24 m·s-1)

    式(2)表明,在確定的流量條件下,SVQS的自然旋風(fēng)長度隨隔流筒直徑的減小而增大,這與魏耀東等[14]所測實(shí)驗結(jié)果相一致;同時自然旋風(fēng)長度隨封閉罩筒體直徑的增加而增大,這與Alexander[6]所測實(shí)驗結(jié)果相一致。因此,本文給出的SVQS自然旋風(fēng)長度的計算公式能夠較準(zhǔn)確地描述SVQS旋流快分系統(tǒng)的旋轉(zhuǎn)流動特性,對工程設(shè)計有一定的指導(dǎo)意義。

    3.5SVQS旋流快分的結(jié)構(gòu)優(yōu)化

    通過上述分析可知,SVQS旋流快分的自然旋風(fēng)長度隨旋流頭噴出口氣速的增加而增大,隨汽提氣速的增加而減小。當(dāng)噴口氣速Vs為22 m·s-1、汽提氣速Vq分別為0和0.90 m·s-1時,旋風(fēng)尾渦距錐形擋板上沿的軸向高度分別為 970和 1360 mm。因此,模擬過程中設(shè)定的 SVQS旋流快分系統(tǒng)的封閉罩筒體有相當(dāng)一部分區(qū)域處于無旋流區(qū)域。為此,本文在封閉罩直筒段高度降低 800 mm(隔流筒下端面距錐形擋板上沿的高度為 2350mm)、旋流頭噴口氣速為24 m·s-1的條件下,模擬計算了SVQS旋流快分的氣相流場,其流動軌跡如圖10(b)所示。由圖10可知,當(dāng)噴出口氣速Vs=24 m·s-1時,最大切向速度衰減達(dá)88%時的截面為錐形擋板上沿以下40 mm的截面位置處,此時的封閉罩直筒段全部處于旋流穩(wěn)定區(qū),SVQS的自然旋風(fēng)長度約為 2390 mm。在該條件下根據(jù)式(2)計算得到的自然旋風(fēng)長度為2474 mm,其與模擬計算值間的相對偏差約為3.50%,說明采用式(2)能夠相對簡單、準(zhǔn)確地計算提升管出口SVQS旋流快分的自然旋風(fēng)長度。

    4 結(jié) 論

    本文采用RNG k-ε湍流模型模擬了SVQS旋流快分系統(tǒng)內(nèi)的氣相流場并得到了SVQS旋流快分自然旋風(fēng)長度的影響因素,結(jié)論如下。

    (1)通過模擬發(fā)現(xiàn),隔流筒下端面最大切向速度衰減達(dá)88%時的截面與SVQS旋風(fēng)尾渦截面相吻合,這與實(shí)驗測量結(jié)果一致,說明采用RNG k-ε湍流模型對 SVQS旋流快分系統(tǒng)進(jìn)行模擬是有效的,并據(jù)此定義了提升管出口SVQS的自然旋風(fēng)長度。

    (2)通過考察不同旋流頭噴口氣速及不同汽提線速條件下SVQS不同軸向高度截面上軸切向速度徑向分布的規(guī)律發(fā)現(xiàn),SVQS旋流快分的自然旋風(fēng)長度隨旋流頭噴出口氣速的增加和汽提線速的減小而增加。

    符號說明

    Ao——出口面積,m2

    Ar——提升管入口截面積,m2

    a ——SVQS旋流頭噴出口高度,mm

    b ——SVQS旋流頭噴出口寬度,mm

    C1ε,C2ε,Cμ——常數(shù)

    D ——封閉罩直徑,m

    De——隔流筒直徑,m

    Dh——水力學(xué)直徑,m

    d ——平衡管截面邊長,mm

    Eij——應(yīng)變率張量

    Gk——平均速度梯度引起的產(chǎn)生項,J

    I ——湍動強(qiáng)度

    k——湍動動能,m2·s-2

    L——自然旋風(fēng)長度,m

    n——旋流臂總數(shù)

    Po——出口壓力,Pa

    Pt——?dú)怏w總壓,Pa

    p——流體壓力,Pa

    Q——入口氣量,m3·s-1

    R——封閉罩半徑,mm

    r——測點(diǎn)距中心軸線半徑,mm

    S——源項

    u——x方向速度,m·s-1

    Ve——隔流筒與提升管環(huán)形空間表觀氣速,m·s-1

    Vo——提升管與封閉罩環(huán)形空間表觀氣速,m·s-1

    Vq——汽提氣速,m·s-1

    Vr——提升管入口氣速,m·s-1

    Vs——旋流頭噴出口氣速,m·s-1

    Vt——切向速度,m·s-1

    V

    z——軸向速度,m·s-1

    v——y方向速度,m·s-1

    w——z方向速度,m·s-1

    z——截面高度,m

    αk,αε,β——常數(shù)

    Γ——擴(kuò)散系數(shù)

    ε——湍流耗散率,m2·s-3

    η0——常數(shù)

    μ——動力黏度,Pa·s

    μeff——有效黏度,Pa·s

    μt——湍動黏度,Pa·s

    ρ——密度,kg·m-3

    References

    [1] 盧春喜, 徐文清, 魏耀東, 等. 新型緊湊式催化裂化沉降系統(tǒng)的實(shí)驗研究 [J]. 石油學(xué)報(石油加工), 2007, 23 (6): 6-12. DOI: 10.3969/j.issn.1001-8719.2007.06.002.

    LU C X, XUN W Q, WEI Y D, et al. Experimental studies of a novel compact FCC disengage [J]. Acta Petrolei Sinica (Petroleum Processing Section), 2007, 23 (6): 6-12. DOI: 10.3969/j.issn.1001-8719.2007.06.002.

    [2]孫鳳俠, 盧春喜, 時銘顯. 旋流快分器內(nèi)氣相流場的實(shí)驗與數(shù)值模擬研究 [J]. 石油大學(xué)學(xué)報(自然科學(xué)版), 2005, 29 (3): 106-111. DOI: 10.3321/j.issn:1000-5870.2005.03.024.

    SUN F X, LU C X, SHI M X. Experiment and numerical simulation of gas flow field in new vortex quick separation system [J]. Journal of the University of Petroleum, 2005, 29 (3): 106-111. DOI: 10.3321/j.issn:1000-5870.2005.03.024.

    [3]孫鳳俠, 盧春喜, 時銘顯. 催化裂化沉降器新型高效旋流快分器內(nèi)氣固兩相流動 [J]. 化工學(xué)報, 2005, 56 (12): 2280-2287.

    SUN F X, LU C X, SHI M X. Numerical simulation of gas particles flow field in new vortex quick separation system of FCC disengage [J]. Journal of Chemical Industry and Engineering (China), 2005, 56 (12): 2280-2287.

    [4]孫鳳俠, 盧春喜, 時銘顯. 催化裂化沉降器新型高效旋流快分器的結(jié)構(gòu)優(yōu)化與分析 [J]. 中國石油大學(xué)學(xué)報(自然科學(xué)版), 2007, 31 (5): 109-113. DOI: 10.3321/j.issn:1000-5870.2007.05.024.

    SUN F X, LU C X, SHI M X. Configuration optimization and analysis of new ly developed vortex quick separator for FCC disengage [J]. Journal of China University of Petroleum (Natural Science Edition), 2007, 31 (5): 109-113. DOI: 10.3321/j.issn:1000-5870.2007.05.024.

    [5]胡艷華, 汪洋, 盧春喜, 等. 催化裂化提升管出口旋流快分系統(tǒng)內(nèi)隔流筒結(jié)構(gòu)的優(yōu)化改進(jìn) [J]. 石油學(xué)報(石油加工), 2008, 24 (2): 177-183. DOI: 10.3969/j.issn.1001-8719.2008.02.010.

    HU Y H, WANG Y, LU C X, et al. Improvement of the isolate cylinder in vortex quick separation system at FCCU riser outlet [J]. Acta Petrolei Sinica (Petroleum Processing Section), 2008, 24 (2): 177-183. DOI: 10.3969/j.issn.1001-8719.2008.02.010.

    [6]ALEXANDER R M. Fundamentals of cyclone design and operation [J]. Proceedings of the Australia Institute of Minerals and Metals (New Series), 1949, 152 (3): 203-228.

    [7]LEITH D, LICTH W. The collection efficiency of cyclone type particle collectors: a new theoretical approach [J]. AIChE Symp. Series, 1972, 68 (126): 196-206.

    [8]DIETZ P W. Collection efficiency of cyclone separators [J]. AIChE Journal, 1981, 27 (27): 888-892.

    [9]HOFFMANN A C, JONEGE R D, ARENDS H, et al. Evidence of the natural vortex length and its effect on the separation efficiency of gas cyclones [J]. Filtration & Separation, 1995, 32 (8): 799-804.

    [10]BRYANT H S, SILVERMAN R M, ZENE F A. How dust in gas affects cyclone pressure drop [J]. HydrocarbonProcessing, 1983, 62 (6): 87-90.

    [11]姬忠禮, 吳小林, 時銘顯. 旋風(fēng)分離器自然旋風(fēng)長的實(shí)驗研究 [J].石油學(xué)報(石油加工), 1993, 9 (4): 86-91.

    JI Z L, WU X L, SHI M X. Experimental research on the natural turning length in the cyclone [J]. Acta Petrolei Sinica (Petroleum Processing Section), 1993, 9 (4): 86-91.

    [12]杜德喜, 朱麗娟, 董秀奇. CFD模擬旋風(fēng)分離器自然旋風(fēng)長 [J].礦山機(jī)械, 2007, 35 (4): 94-96.

    DU D X, ZHU L J, DONG X Q. Simulation to the length of the natural vortex in cyclone separator with CFD [J]. M ining & Processing Equipment, 2007, 35 (4): 94-96.

    [13]高翠芝, 孫國剛, 董瑞倩. 旋風(fēng)分離器旋風(fēng)長度的分析計算 [J].石油學(xué)報(石油加工), 2012, 28 (1): 94-98. DOI: 10.3969/j.issn.1001-8719.2012.01.017.

    GAO C Z, SUN G G, DONG R Q. Analysis calculation of the vortex length in a gas cyclone [J]. Acta Petrolei Sinica (Petroleum Processing Section), 2012, 28 (1): 94-98. DOI: 10.3969/j.issn.1001-8719.2012.01.017.

    [14]魏耀東, 張靜, 宋健斐, 等. 旋風(fēng)分離器自然旋風(fēng)長的實(shí)驗研究[J]. 熱能動力工程, 2010, 25 (2): 206-210.

    WEI Y D, ZHANG J, SONG J F, et al. Experimental study of the natural cyclone length of a cyclone separator [J]. Journal of Engineering for Thermal Energy and Power, 2010, 25 (2): 206-210.

    [15]錢付平, 章名耀. 基于響應(yīng)曲面法旋風(fēng)分離器的自然旋風(fēng)長 [J].東南大學(xué)學(xué)報(自然科學(xué)版), 2006, 36 (2): 247-251. DOI: 10.3321/j.issn:1001-0505.2006.02.014.

    QIAN F P, ZHANG M Y. Natural vortex lengths of cyclone separators based on response surface methodology [J]. Journal of Southeast University (Natural Science Edition), 2006, 36 (2): 247-251. DOI: 10.3321/j.issn:1001-0505.2006.02.014.

    [16]王福軍. 計算流體動力學(xué)分析——CFD軟件原理與應(yīng)用 [M]. 北京:清華大學(xué)出版社, 2009: 122-125.

    WANG F J. Computational Hydrodynamics Analysis—CFD Software Theory and Application [M]. Beijing: Tsinghua University Press,2009: 122-125.

    [17]王江云, 毛羽, 劉美麗, 等. 用改進(jìn)的RNG k-ε模型模擬旋風(fēng)分離器內(nèi)的強(qiáng)旋流動 [J]. 石油學(xué)報(石油加工), 2010, 26 (1): 8-13. DOI: 10.3969/j.issn.1001-8719.2010.01.002.

    WANG J Y, MAO Y, LIU M L, et al. Numerical simulation of strongly sw irling flow in cyclone separator by using an advanced RNG k-ε model [J]. Acta Petrolei Sinica (Petroleum Processing Section),2010, 26 (1): 8-13. DOI: 10.3969/j.issn.1001-8719.2010.01.002.

    [18]程兆龍, 姚秀穎, 鄂承林, 等. 帶隔流筒旋流快分系統(tǒng)新型工業(yè)催化裂化沉降器內(nèi)氣相流場的數(shù)值模擬 [J]. 過程工程學(xué)報, 2016,16 (1): 86-95.

    CHENG Z L, YAO X Y, E C L, et al. Numerical simulation of gas phase flow field in a novel FCC disengager w ith flow partition cylinder quick separation system [J]. The Chinese Journal of Process Engineering, 2016, 16 (1): 86-95.

    [19]宋健斐, 魏耀東, 時銘顯. 蝸殼式旋風(fēng)分離器內(nèi)氣相流場非軸對稱特性分析 [J]. 化工學(xué)報, 2007, 58 (5): 1091-1096.

    SONG J F, WEI Y D, SHI M X. Analysis of asymmetry of gas phase flow field in volute cyclone [J]. Journal of Chemical Industry and Engineering (China), 2007, 58 (5): 1091-1096.

    [20]WEI Y D, SONG J F, SHI M X, et al. Numerical simulation of the asymmetry of gas-phase flow field in a volute cyclone separator [J]. Progress in Natural Science, 2005, 15: 98-104.

    Numerical simulation of natural vortex length in riser term ination w ith SVQS system

    CHENG Zhaolong, E Chenglin, LU Chunxi
    (State Key Laboratory of Heavy Oil Processing, China University of Petroleum, Beijing 102249, China)

    To optim ize the cylinder height of the closing section, the gas phase flow field and the natural vortex length in a ?600 mm×4150 mm super vortex quick separation (SVQS) system were simulated by using the RNG k-ε turbulent model. The simulation results show that the cross-section location of the vortex end was agreeable to that of 88% attenuation of the maximum tangential velocity in the entrance section of flow partition column,which hence defined the natural vortex length. From the tangential velocity distribution of gas phase in SVQS, the natural vortex length of the quick separation was found to increase w ith the increase of spout velocity and the decrease of stripping velocity. Based on the numerical simulation and structural characteristics of SVQS, an equation for calculating the natural vortex length in SVQS was developed.

    riser; separation; computational fluid dynam ics

    date: 2016-03-30.

    Prof. LU Chunxi, lcx725@sina.com

    supported by the National Basic Research Program of China (2012CB215000).

    O 351; TE 624; TQ 021.1

    A

    0438—1157(2016)08—3347—10

    10.11949/j.issn.0438-1157.20160381

    2016-03-30收到初稿,2016-07-16收到修改稿。

    聯(lián)系人:盧春喜。第一作者:程兆龍(1990—),男,碩士研究生。

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計劃項目(2012CB215000)。

    猜你喜歡
    切向速度汽提氣速
    旋風(fēng)分離器內(nèi)氣相旋轉(zhuǎn)流不穩(wěn)定性的實(shí)驗研究
    提高含硫氣田水汽提效率技術(shù)及其應(yīng)用
    元壩氣田含硫污水負(fù)壓汽提脫硫工藝影響因素分析與優(yōu)化
    雙旋流耦合式旋流反應(yīng)器內(nèi)切向速度分布研究
    濕法煙氣脫硫吸收塔阻力特性實(shí)驗研究
    汽提條件對兩段提升管催化裂化待生劑汽提過程的影響
    新型折板除霧器的流場和壓降數(shù)值模擬
    D120 mm流化床冷模實(shí)驗研究
    化工科技(2014年5期)2014-06-09 05:17:22
    重介質(zhì)微型旋流器內(nèi)切向速度的數(shù)值模擬
    化學(xué)汽提新工藝
    欧美激情久久久久久爽电影| 精品欧美国产一区二区三| 国产精品,欧美在线| 精品人妻1区二区| 国产亚洲精品久久久久久毛片| 又黄又爽又免费观看的视频| 久久精品国产亚洲av涩爱 | 伦精品一区二区三区| 老司机福利观看| 给我免费播放毛片高清在线观看| 亚洲精品粉嫩美女一区| 亚洲熟妇中文字幕五十中出| 97碰自拍视频| 精品一区二区三区人妻视频| 日韩中文字幕欧美一区二区| 天堂√8在线中文| 国产精品一区二区性色av| 伊人久久精品亚洲午夜| 久久久久免费精品人妻一区二区| 联通29元200g的流量卡| 日本三级黄在线观看| 精品一区二区三区av网在线观看| 日本熟妇午夜| 又紧又爽又黄一区二区| 亚洲av二区三区四区| 亚洲国产精品sss在线观看| 亚洲自偷自拍三级| 桃红色精品国产亚洲av| 一卡2卡三卡四卡精品乱码亚洲| 99精品在免费线老司机午夜| 人妻夜夜爽99麻豆av| 国产精品福利在线免费观看| 丰满人妻一区二区三区视频av| 国产一区二区三区在线臀色熟女| 国产成人av教育| 精品久久久久久久久久免费视频| 精品不卡国产一区二区三区| 3wmmmm亚洲av在线观看| 国产精品爽爽va在线观看网站| 给我免费播放毛片高清在线观看| 亚洲人成网站高清观看| 欧美区成人在线视频| 在线观看一区二区三区| 两个人视频免费观看高清| 亚洲国产色片| 99精品在免费线老司机午夜| 亚洲av中文av极速乱 | 十八禁国产超污无遮挡网站| 又紧又爽又黄一区二区| 亚洲精品色激情综合| 午夜福利在线观看吧| 成人二区视频| 久久久久久久久中文| 国产精品98久久久久久宅男小说| 久久精品国产自在天天线| 日本爱情动作片www.在线观看 | 韩国av一区二区三区四区| 99久久久亚洲精品蜜臀av| 在线免费观看的www视频| 色精品久久人妻99蜜桃| 亚洲美女搞黄在线观看 | 国产一区二区在线观看日韩| 哪里可以看免费的av片| 99久久精品国产国产毛片| 成人鲁丝片一二三区免费| 亚洲在线观看片| 日韩国内少妇激情av| 日本-黄色视频高清免费观看| 日日撸夜夜添| 在线免费十八禁| 国产精品永久免费网站| 在线a可以看的网站| 久久久久久九九精品二区国产| 国产高清视频在线观看网站| 成年人黄色毛片网站| 成人美女网站在线观看视频| 国内精品美女久久久久久| 深夜a级毛片| 18禁黄网站禁片免费观看直播| 色综合亚洲欧美另类图片| 色综合站精品国产| 最近在线观看免费完整版| 性欧美人与动物交配| 亚洲无线观看免费| 深夜a级毛片| 乱人视频在线观看| 两人在一起打扑克的视频| 九色国产91popny在线| 999久久久精品免费观看国产| eeuss影院久久| 俄罗斯特黄特色一大片| 极品教师在线视频| 亚洲经典国产精华液单| 99久久精品热视频| 俺也久久电影网| 国产成人a区在线观看| 1000部很黄的大片| 女生性感内裤真人,穿戴方法视频| 国产精品久久视频播放| 国产精品一区二区免费欧美| 两个人的视频大全免费| 五月玫瑰六月丁香| 麻豆av噜噜一区二区三区| 99精品在免费线老司机午夜| 国产老妇女一区| 男人舔女人下体高潮全视频| 69av精品久久久久久| 精品日产1卡2卡| 中文字幕高清在线视频| 亚洲精品在线观看二区| 日本爱情动作片www.在线观看 | 少妇熟女aⅴ在线视频| 精品人妻1区二区| 我要搜黄色片| 欧美高清性xxxxhd video| 国产毛片a区久久久久| av天堂在线播放| 久久精品国产亚洲av香蕉五月| 美女xxoo啪啪120秒动态图| 51国产日韩欧美| 国产高清三级在线| 国产三级在线视频| 欧美成人性av电影在线观看| 亚洲男人的天堂狠狠| 欧美日韩国产亚洲二区| 中文字幕av成人在线电影| 天堂动漫精品| 国产日本99.免费观看| 99热精品在线国产| 2021天堂中文幕一二区在线观| 九色成人免费人妻av| videossex国产| 极品教师在线视频| av专区在线播放| 中文字幕av成人在线电影| 亚洲av免费在线观看| 成年人黄色毛片网站| 亚洲在线自拍视频| 国产在线男女| 男女啪啪激烈高潮av片| 日韩欧美免费精品| 精品人妻熟女av久视频| 人人妻人人看人人澡| 22中文网久久字幕| 色综合站精品国产| 人妻久久中文字幕网| 精品福利观看| 天堂动漫精品| 能在线免费观看的黄片| 中国美白少妇内射xxxbb| 亚洲欧美日韩无卡精品| АⅤ资源中文在线天堂| 国产人妻一区二区三区在| 成人美女网站在线观看视频| 久久久久久九九精品二区国产| 日韩强制内射视频| aaaaa片日本免费| 黄色配什么色好看| 国产精品一区二区三区四区久久| 一级黄片播放器| 亚洲,欧美,日韩| 国产一区二区三区av在线 | 欧美性猛交╳xxx乱大交人| 热99re8久久精品国产| 成人欧美大片| 亚洲最大成人中文| 精品久久久久久久末码| 久9热在线精品视频| 欧美日韩亚洲国产一区二区在线观看| 18+在线观看网站| av.在线天堂| 一区福利在线观看| 国产精品免费一区二区三区在线| 91久久精品国产一区二区三区| 给我免费播放毛片高清在线观看| 天堂影院成人在线观看| 国产精品不卡视频一区二区| 国产激情偷乱视频一区二区| 看免费成人av毛片| 嫩草影视91久久| 久久精品国产自在天天线| www日本黄色视频网| 有码 亚洲区| 男人的好看免费观看在线视频| 如何舔出高潮| 99久久九九国产精品国产免费| 久久亚洲精品不卡| 国模一区二区三区四区视频| 级片在线观看| 黄色配什么色好看| 一区二区三区高清视频在线| 男女啪啪激烈高潮av片| 美女高潮喷水抽搐中文字幕| 亚洲精华国产精华精| 久久精品国产亚洲av天美| 在线观看免费视频日本深夜| 国产精品一及| 亚洲av免费高清在线观看| 啦啦啦观看免费观看视频高清| 国产精品一区二区三区四区免费观看 | 国产欧美日韩精品一区二区| 日本-黄色视频高清免费观看| 天堂动漫精品| 给我免费播放毛片高清在线观看| 免费黄网站久久成人精品| 亚洲美女搞黄在线观看 | 国产精品1区2区在线观看.| 一区福利在线观看| 久久久久久久精品吃奶| 一卡2卡三卡四卡精品乱码亚洲| 国产精品1区2区在线观看.| 国产三级中文精品| 3wmmmm亚洲av在线观看| 丰满乱子伦码专区| 午夜精品久久久久久毛片777| 亚洲人成伊人成综合网2020| 又爽又黄无遮挡网站| 99热精品在线国产| 婷婷亚洲欧美| 精品人妻视频免费看| 国产精品三级大全| 亚洲欧美清纯卡通| 免费看日本二区| 啦啦啦韩国在线观看视频| 熟女电影av网| 国产 一区精品| 丰满乱子伦码专区| 美女cb高潮喷水在线观看| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美国产在线观看| 在线观看66精品国产| 18禁黄网站禁片免费观看直播| 少妇的逼好多水| 久久久久九九精品影院| 精品人妻一区二区三区麻豆 | 亚洲中文日韩欧美视频| 国产不卡一卡二| 国产精华一区二区三区| 欧美最黄视频在线播放免费| av国产免费在线观看| 精品久久国产蜜桃| 级片在线观看| 一级a爱片免费观看的视频| 国产成人影院久久av| 精品久久久噜噜| 自拍偷自拍亚洲精品老妇| 一边摸一边抽搐一进一小说| 伦精品一区二区三区| 国产精品,欧美在线| 无人区码免费观看不卡| 国产亚洲精品久久久久久毛片| 淫秽高清视频在线观看| 直男gayav资源| 99精品在免费线老司机午夜| 不卡一级毛片| 大型黄色视频在线免费观看| 伦精品一区二区三区| 国产伦一二天堂av在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 精品一区二区三区人妻视频| 国产精品人妻久久久影院| videossex国产| 麻豆成人午夜福利视频| 亚洲久久久久久中文字幕| 免费av观看视频| 国产精品1区2区在线观看.| 久久久成人免费电影| 久久久国产成人免费| 日本熟妇午夜| 日韩精品有码人妻一区| 麻豆成人av在线观看| 1000部很黄的大片| 黄色女人牲交| 国产黄色小视频在线观看| 一级av片app| 一个人看视频在线观看www免费| 亚洲精品一区av在线观看| 18禁在线播放成人免费| 老司机午夜福利在线观看视频| 婷婷亚洲欧美| 精品乱码久久久久久99久播| 麻豆国产av国片精品| 国产 一区 欧美 日韩| 欧美日韩瑟瑟在线播放| av中文乱码字幕在线| 国产黄色小视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产精品98久久久久久宅男小说| 美女高潮喷水抽搐中文字幕| 国产美女午夜福利| 日本一本二区三区精品| 亚洲av不卡在线观看| 欧美日韩瑟瑟在线播放| 99热网站在线观看| 真人一进一出gif抽搐免费| 波多野结衣巨乳人妻| 国产午夜精品久久久久久一区二区三区 | av在线老鸭窝| 国产爱豆传媒在线观看| 俄罗斯特黄特色一大片| 亚洲图色成人| 动漫黄色视频在线观看| 国产一区二区在线av高清观看| 久久精品国产亚洲网站| 国产午夜精品论理片| 国产三级在线视频| 日本 欧美在线| 可以在线观看毛片的网站| 97碰自拍视频| av视频在线观看入口| 老女人水多毛片| 麻豆久久精品国产亚洲av| 亚洲精品国产成人久久av| 国产日本99.免费观看| 色综合色国产| 国产精品免费一区二区三区在线| 日韩亚洲欧美综合| 国产午夜精品论理片| 亚洲美女视频黄频| 1024手机看黄色片| 国产伦在线观看视频一区| 久久精品国产清高在天天线| 国产精品一区www在线观看 | 欧美性猛交黑人性爽| 国产视频内射| 亚洲av五月六月丁香网| 亚洲熟妇中文字幕五十中出| 少妇人妻精品综合一区二区 | 丝袜美腿在线中文| 久久久精品大字幕| av中文乱码字幕在线| 国产精品久久久久久av不卡| 亚洲三级黄色毛片| 深夜a级毛片| 亚洲va在线va天堂va国产| 成人综合一区亚洲| 日本与韩国留学比较| 亚洲精品影视一区二区三区av| 日韩大尺度精品在线看网址| 欧美日本亚洲视频在线播放| 免费人成视频x8x8入口观看| 色哟哟·www| 午夜爱爱视频在线播放| 国产免费男女视频| 日韩欧美一区二区三区在线观看| 欧美色视频一区免费| 91狼人影院| 搡老岳熟女国产| 嫁个100分男人电影在线观看| 人人妻人人看人人澡| 日本 av在线| 乱系列少妇在线播放| 伊人久久精品亚洲午夜| 亚洲综合色惰| 国产中年淑女户外野战色| 少妇熟女aⅴ在线视频| 天堂动漫精品| 久久草成人影院| 永久网站在线| 国产伦人伦偷精品视频| 欧美3d第一页| 成年版毛片免费区| 少妇丰满av| 日本欧美国产在线视频| netflix在线观看网站| 美女免费视频网站| 国产精华一区二区三区| 少妇的逼水好多| 久久精品夜夜夜夜夜久久蜜豆| 成年人黄色毛片网站| 欧美日韩瑟瑟在线播放| 五月玫瑰六月丁香| xxxwww97欧美| 亚洲经典国产精华液单| 国产成人福利小说| 亚洲男人的天堂狠狠| 免费在线观看影片大全网站| 国产伦人伦偷精品视频| 中文在线观看免费www的网站| 久久精品国产自在天天线| 国产成人影院久久av| 国产精品野战在线观看| 欧美+亚洲+日韩+国产| 最近最新免费中文字幕在线| 国产探花在线观看一区二区| 无遮挡黄片免费观看| 在线观看午夜福利视频| 精品国内亚洲2022精品成人| 国产精品久久久久久亚洲av鲁大| 欧美黑人欧美精品刺激| 久久久色成人| 可以在线观看毛片的网站| av福利片在线观看| 久久久久久大精品| 色尼玛亚洲综合影院| 一区二区三区高清视频在线| 天天躁日日操中文字幕| 少妇的逼水好多| 99热网站在线观看| 亚洲成人中文字幕在线播放| 村上凉子中文字幕在线| 国产精品久久视频播放| 免费av观看视频| 婷婷六月久久综合丁香| 免费在线观看成人毛片| 麻豆成人av在线观看| 狠狠狠狠99中文字幕| 国产伦精品一区二区三区四那| 亚洲成人中文字幕在线播放| 高清毛片免费观看视频网站| 国产精品一区二区三区四区免费观看 | 男女那种视频在线观看| 午夜免费成人在线视频| 午夜精品久久久久久毛片777| 日韩,欧美,国产一区二区三区 | 禁无遮挡网站| 在线观看舔阴道视频| 亚洲18禁久久av| 在线播放国产精品三级| 一个人看的www免费观看视频| 亚洲精品色激情综合| 国产v大片淫在线免费观看| 欧美三级亚洲精品| 亚洲熟妇中文字幕五十中出| 欧美中文日本在线观看视频| 午夜视频国产福利| 男女做爰动态图高潮gif福利片| 毛片一级片免费看久久久久 | 人人妻人人澡欧美一区二区| 色播亚洲综合网| 看片在线看免费视频| 精品一区二区免费观看| 午夜激情福利司机影院| 亚洲中文字幕日韩| 国产精品一区二区三区四区久久| 国产单亲对白刺激| 欧美区成人在线视频| 综合色av麻豆| 人妻久久中文字幕网| 国产精品不卡视频一区二区| eeuss影院久久| 熟女人妻精品中文字幕| 在线观看舔阴道视频| 久久精品国产清高在天天线| 好男人在线观看高清免费视频| 亚洲国产精品合色在线| 搡女人真爽免费视频火全软件 | 国产一区二区在线av高清观看| 欧洲精品卡2卡3卡4卡5卡区| 1000部很黄的大片| 男女之事视频高清在线观看| 十八禁国产超污无遮挡网站| 自拍偷自拍亚洲精品老妇| 在线观看午夜福利视频| 亚洲va日本ⅴa欧美va伊人久久| 88av欧美| 校园春色视频在线观看| 国产一区二区三区av在线 | 欧美性猛交╳xxx乱大交人| 天堂网av新在线| 88av欧美| 国产av一区在线观看免费| 国产精品国产高清国产av| 最近最新中文字幕大全电影3| www日本黄色视频网| 日本撒尿小便嘘嘘汇集6| 精品福利观看| 亚洲内射少妇av| 免费在线观看日本一区| 国产精品一及| 麻豆av噜噜一区二区三区| 国产午夜福利久久久久久| 性插视频无遮挡在线免费观看| 在线观看舔阴道视频| av在线天堂中文字幕| 国产精品一区二区免费欧美| 久久久久久久久大av| 成年女人毛片免费观看观看9| 亚洲精品亚洲一区二区| 美女cb高潮喷水在线观看| 日本成人三级电影网站| 国产精品自产拍在线观看55亚洲| 少妇的逼水好多| 精品久久久久久久久av| 日韩欧美国产在线观看| 男人舔奶头视频| 日韩高清综合在线| 国产精品一区二区三区四区免费观看 | 身体一侧抽搐| 搡老妇女老女人老熟妇| 国内精品久久久久精免费| 亚洲综合色惰| 亚洲美女搞黄在线观看 | 久久99热这里只有精品18| 成人av在线播放网站| 亚洲av二区三区四区| 99riav亚洲国产免费| 日本黄色视频三级网站网址| 日本欧美国产在线视频| 日韩一本色道免费dvd| 日韩 亚洲 欧美在线| 国产白丝娇喘喷水9色精品| 精品久久久久久久久av| 精品午夜福利视频在线观看一区| 熟女电影av网| 男女视频在线观看网站免费| 日韩一区二区视频免费看| 亚洲国产精品sss在线观看| 国产美女午夜福利| 又黄又爽又免费观看的视频| 久久久久久久久大av| 如何舔出高潮| 亚洲美女黄片视频| 97超级碰碰碰精品色视频在线观看| 亚洲自拍偷在线| 亚洲一区高清亚洲精品| 国产乱人视频| 毛片女人毛片| 麻豆成人av在线观看| 一边摸一边抽搐一进一小说| av福利片在线观看| 嫩草影院新地址| 色5月婷婷丁香| 深爱激情五月婷婷| 久久久久久久亚洲中文字幕| 亚洲av.av天堂| 女的被弄到高潮叫床怎么办 | 三级毛片av免费| 永久网站在线| 嫩草影院入口| 国产精品不卡视频一区二区| 国产精品一区二区三区四区久久| 99热这里只有是精品50| 久久久久久伊人网av| 日本三级黄在线观看| 国内揄拍国产精品人妻在线| 亚洲一区二区三区色噜噜| 欧美zozozo另类| 欧美区成人在线视频| 午夜老司机福利剧场| 亚洲国产精品合色在线| 九九久久精品国产亚洲av麻豆| 国产视频内射| 国产av在哪里看| 成年免费大片在线观看| 久久精品91蜜桃| 日本欧美国产在线视频| 伦理电影大哥的女人| 久久久久久久精品吃奶| 亚洲色图av天堂| 一级毛片久久久久久久久女| 国产伦一二天堂av在线观看| 观看美女的网站| 久久精品国产鲁丝片午夜精品 | 人妻少妇偷人精品九色| 嫩草影视91久久| 久久久久九九精品影院| 国内精品美女久久久久久| 国产真实乱freesex| 国产白丝娇喘喷水9色精品| 一边摸一边抽搐一进一小说| 国内精品一区二区在线观看| 日韩av在线大香蕉| 国产又黄又爽又无遮挡在线| 久久久久久伊人网av| 精品久久久久久成人av| 黄片wwwwww| 婷婷丁香在线五月| 国内精品美女久久久久久| 毛片女人毛片| 午夜精品一区二区三区免费看| 国产精品免费一区二区三区在线| 观看免费一级毛片| 欧美黑人欧美精品刺激| 久久草成人影院| 黄色女人牲交| 免费观看精品视频网站| 日韩大尺度精品在线看网址| 自拍偷自拍亚洲精品老妇| 99在线视频只有这里精品首页| 熟女电影av网| 人妻少妇偷人精品九色| 色哟哟哟哟哟哟| 精品人妻熟女av久视频| 免费av不卡在线播放| 亚洲国产精品成人综合色| 熟女电影av网| 国产91精品成人一区二区三区| 亚洲国产精品久久男人天堂| 3wmmmm亚洲av在线观看| 国产中年淑女户外野战色| 22中文网久久字幕| 亚洲性夜色夜夜综合| av在线亚洲专区| 色噜噜av男人的天堂激情| 亚洲最大成人av| 免费av观看视频| 欧美激情在线99| 国产亚洲精品av在线| 一级毛片久久久久久久久女| 国产精品野战在线观看| 少妇裸体淫交视频免费看高清| 日韩,欧美,国产一区二区三区 | 国产伦精品一区二区三区四那| 他把我摸到了高潮在线观看| 亚洲人成网站高清观看| 亚洲av.av天堂| 男人狂女人下面高潮的视频| 黄色日韩在线| 九色成人免费人妻av| 亚洲avbb在线观看| 他把我摸到了高潮在线观看| 国产在线精品亚洲第一网站| 国产一区二区亚洲精品在线观看| 国产精品野战在线观看|