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

    導(dǎo)流錐式超聲速旋流分離裝置流動特性

    2016-07-07 12:12:38胡大鵬王熒光任文文趙健華劉培啟大連理工大學(xué)化工機械學(xué)院遼寧大連116023
    化工學(xué)報 2016年6期

    胡大鵬,王熒光,任文文,趙健華,劉培啟(大連理工大學(xué)化工機械學(xué)院,遼寧 大連 116023)

    ?

    導(dǎo)流錐式超聲速旋流分離裝置流動特性

    胡大鵬,王熒光,任文文,趙健華,劉培啟
    (大連理工大學(xué)化工機械學(xué)院,遼寧 大連 116023)

    摘要:采用二維軸對稱模型,以理想空氣為介質(zhì)對超聲速旋流分離裝置內(nèi)的流場特性進(jìn)行了數(shù)值模擬。同時搭建實驗平臺,對小壓比條件下影響流場特性和分離性能的結(jié)構(gòu)和操作參數(shù)進(jìn)行了研究。研究結(jié)果表明:進(jìn)出口壓比為1.4同時排液通道外壁張角小于12°時,超聲速噴管擴張段內(nèi)在面積比達(dá)1.27的情形下仍不存在氣動激波,同時對擴張段內(nèi)存在激波產(chǎn)生的情況進(jìn)行分析,得出其原因為在排液口內(nèi)產(chǎn)生反向壓縮波,并向噴管上游移動。通過實驗研究得出,壓比為1.4,面積比為1.27時超聲速旋流分離器分離效率最高,達(dá)到20.5%。

    關(guān)鍵詞:超聲速流動;流體動力學(xué);數(shù)值模擬;激波;離心分離;Mach數(shù)

    2015-12-02收到初稿,2016-03-06收到修改稿。

    聯(lián)系人:劉培啟。第一作者:胡大鵬(1963—),男,教授。

    Received date: 2015-12-02.

    Foundation item: supported by the National Natural Science Foundation of China (21206013, 21476036) and the Basic Research Project of Key Laboratory of Liaoning Provincial Education Department(LZ2015019).

    引 言

    天然氣作為一種潔凈的化工原料被廣泛應(yīng)用于諸多行業(yè),然而從天然氣井直接開采出的天然氣中含有大量的水蒸氣,使得天然氣在輸運和處理過程中易形成液態(tài)水或固態(tài)冰,會大大降低設(shè)備的性能[1]。目前常用的天然氣脫水方法有[2-6]:J-T閥和透平膨脹機、三甘醇脫水技術(shù)和分子篩脫水技術(shù)等。而且天然氣脫水裝置正朝著體積小、操作方便、投資和使用成本低和無污染的方向發(fā)展。超聲速冷凝旋流分離技術(shù)[7-10]作為一種近年來出現(xiàn)的新型脫水技術(shù),采用低溫冷凝法和強大的旋流場可以快速脫除含濕天然氣中的重?zé)N組分,其最大的優(yōu)點是集膨脹制冷和旋流分離于一體,而且由于其較大的軸向速度,介質(zhì)在設(shè)備內(nèi)停留時間短,只有幾毫秒,不易生成水合物,進(jìn)而不會造成設(shè)備和管道的堵塞。因此,與傳統(tǒng)的分離方法相比,超聲速旋流分離技術(shù)具有工藝簡單、結(jié)構(gòu)緊湊和無轉(zhuǎn)動部件的顯著優(yōu)勢,具有廣闊的應(yīng)用前景[11-12]。

    目前,關(guān)于超聲速分離器的結(jié)構(gòu)研究主要集中在Laval噴管、旋流發(fā)生器和用于回收壓力的擴壓器3個方面。Laval噴管是裝置分離段產(chǎn)生均勻超聲速氣流的重要部件,其出口截面與喉部面積比保證了裝置達(dá)到設(shè)計Mach數(shù)的要求,噴管曲線線型直接決定氣流的均勻度[13]。曹學(xué)文等[14]對3種不同的噴管設(shè)計方法進(jìn)行對比分析,得出:當(dāng)噴管漸縮段采用維托辛思基曲線、喉部為一段光滑圓弧、擴張段采用福爾士法進(jìn)行設(shè)計時,噴管出口氣流均勻,達(dá)到設(shè)計要求。楊文等[15]結(jié)合雙三次曲線法、圓弧加直線方法以及邊界層黏性修正對Laval噴管的結(jié)構(gòu)設(shè)計進(jìn)行研究,認(rèn)為邊界層位移厚度沿軸向線性發(fā)展,線性修正角為0.5°。Wyslouzil等[16]在超聲速實驗中采用直線型噴管,漸擴段傾角1.8°,且在噴管喉部采用三次樣條曲線使噴管收縮和擴張段光滑過渡,這樣可以避免在噴管下游產(chǎn)生弱激波。旋流發(fā)生器作為超聲速分離器的另一個核心部件也受到了廣泛的關(guān)注,文獻(xiàn)[17-18]研究了旋流對噴管內(nèi)超聲速流動的影響,對有旋流和無旋流的幾何模型進(jìn)行了數(shù)值和實驗研究,得出旋流運動會引起噴管內(nèi)超聲速流動的不均勻,使噴管內(nèi)出現(xiàn)Mach數(shù)沿截面徑向分布的不一致,且研究結(jié)果表明隨著旋流強度的增加,這種速度沿徑向分布的不一致性會更加明顯。目前,噴管結(jié)構(gòu)和旋流發(fā)生器的研究主要集中在其對超聲速流動穩(wěn)定性和均勻性方面,而用來降低流速,回收壓力能的擴壓器的研究主要集中在能量損失方面。文獻(xiàn)[19]研究了激波前Mach數(shù)和氣體絕熱指數(shù)對擴壓器回收壓力能的影響,得出激波前Mach數(shù)越大,壓力恢復(fù)系數(shù)越小,即壓力能損失越大。因此,Laval噴管、旋流發(fā)生器和擴壓器的設(shè)計要綜合考慮。在滿足設(shè)計要求的前提下,應(yīng)保證噴管內(nèi)氣流的均勻性和更小的壓力能損失。排液口作為超聲速分離器的重要組成部分,位于噴管出口和擴壓器進(jìn)口的相交處,對噴管內(nèi)的超聲速流動和可凝組分的分離效率有著重要的影響,目前關(guān)于排液口的結(jié)構(gòu)及幾何尺寸對噴管內(nèi)流動特性影響的報道還較少。

    本文設(shè)計了一種新型導(dǎo)流錐式超聲速旋流分離器,通過調(diào)節(jié)導(dǎo)流錐的位置實現(xiàn)了超聲速噴管喉部面積可調(diào),噴管擴張段采用圓形截面,與環(huán)形噴管截面相比,壓力損失小,分離效率高,尤其適用于小壓比工況。實現(xiàn)進(jìn)口氣體含濕組分冷凝分離的首要前提是,在超聲速旋流分離擴展段建立穩(wěn)定的低溫旋流場,盡量避免排液口對上游噴管內(nèi)流場的影響。本文首先建立二維軸對稱超聲速旋流數(shù)值模型,以理想干空氣為介質(zhì)對超聲速旋流分離裝置內(nèi)的流場特性進(jìn)行了數(shù)值模擬,重點分析了排液口結(jié)構(gòu)對上游旋流超聲速流場的影響。在對超聲速旋流裝置結(jié)構(gòu)優(yōu)化的基礎(chǔ)上,以含濕空氣為實驗介質(zhì)實驗研究了裝置的分離效率。主要分析不同噴管面積比、排液通道外傾角和壓比對噴管內(nèi)流動特性和設(shè)備分離效率的影響。

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

    為了能準(zhǔn)確反映噴管內(nèi)的流動狀態(tài),同時保證控制方程推導(dǎo)的方便,假設(shè)噴管內(nèi)的流動為二維軸對稱、絕熱、可壓縮流動,建立如下的質(zhì)量、動量和能量守恒方程。

    質(zhì)量守恒方程

    動量守恒方程

    能量守恒方程

    式中,ui(i=1,2)分布為流場x、y方向的速度;p、T為流場的壓力和溫度;ρ為流體密度;k為流動介質(zhì)的傳熱系數(shù);cp為比熱容;τij為湍流應(yīng)力張力;ρgi為i方向上的體積力。由于超聲速旋流器進(jìn)口氣體在噴管內(nèi)膨脹加速到Mach數(shù)Ma > 1后,在排液口附近,受壁面幾何形狀突變的影響會在噴管漸擴段內(nèi)形成駐激波,產(chǎn)生壓力和溫度的間斷面。求解該問題時選用建立在Godunov方法基礎(chǔ)上的Roe通量差分分裂二階迎風(fēng)離散格式[20],在時間域上采用全隱式時間積分方案。根據(jù)噴管內(nèi)流體的流動特性和不同湍流模型的適用范圍而采用標(biāo)準(zhǔn)的k-ε模型。求解方程時采用所有變量聯(lián)立的耦合式解法。由于本文研究小壓比下干空氣在超聲速旋流分離裝置中的流動特性,氣體狀態(tài)方程采用理想氣體狀態(tài)方程[21]。

    2 模型驗證

    本文采用文獻(xiàn)[22]中的實驗數(shù)據(jù)對所建模型進(jìn)行驗證,實驗噴管的結(jié)構(gòu)如圖1所示,操作參數(shù)和結(jié)構(gòu)參數(shù)見表1、表2。

    圖1 實驗噴管的幾何模型Fig.1 Geometrical model of experimental nozzle

    表1 實驗噴管的操作參數(shù)Table 1 Operating parameters of experimental nozzle

    表2 實驗噴管的結(jié)構(gòu)參數(shù)Table 2 Structural parameters of experimental nozzle

    圖2 噴管內(nèi)靜壓沿軸線的分布Fig.2 Pressure distribution along axis of nozzle

    根據(jù)實驗噴管的操作參數(shù),用于模型驗證的邊界條件:進(jìn)出口均采用壓力邊界條件,數(shù)值分別為0.46和0.10 MPa;進(jìn)出口溫度分別為288和291 K。圖2為模擬所得噴管中心靜壓與實驗測量值的對比,x表示噴管任意截面與進(jìn)口的軸向距離,噴管喉部位置為x = 50 mm。從圖中可以看出,當(dāng)x在0~100 mm范圍內(nèi)時,模擬值與實驗值基本吻合;當(dāng)x > 100 mm時,則存在一定的偏差。主要原因是實驗采用的介質(zhì)是含濕氣體,在x = 100 mm位置濕空氣中水蒸氣產(chǎn)生自發(fā)凝結(jié)現(xiàn)象釋放大量凝結(jié)潛熱使氣體膨脹偏離等熵膨脹過程。本文數(shù)值模型中采用的介質(zhì)是干空氣,并沒有考慮含濕組分的凝結(jié)特性,使得噴管軸線上的壓降過程繼續(xù)沿著等熵膨脹線下降,以至于x > 100 mm之后位置的靜壓模擬值相對于實驗值略低。對比x在0~100 mm范圍內(nèi)的靜壓模擬和實驗曲線,本文所建立的數(shù)學(xué)模型能正確反映干空氣在超聲速噴管內(nèi)的流場特性。

    3 幾何模型及數(shù)值模擬結(jié)果分析

    3.1幾何模型

    當(dāng)不考慮裝置前的旋流器時,超聲速旋流分離器幾何模型具有旋轉(zhuǎn)軸對稱的特性,因此采用二維軸對稱模型對超聲速旋流分離器內(nèi)的流場特性進(jìn)行數(shù)值研究。圖3(a)、(b)分別為二維全模型和排液口結(jié)構(gòu)示意圖。影響超聲速旋流器內(nèi)流場特性的結(jié)構(gòu)參數(shù)有:噴管漸縮和漸擴長度L1、L2,排液間隙Lb,排液通道內(nèi)外傾角及其出口面積與喉部面積之比A2/A1。下面采用本文建立的數(shù)值模型分別計算以上結(jié)構(gòu)參數(shù)對噴管內(nèi)流場特性的影響。以超聲速旋流器高壓進(jìn)口邊為軸向長度原點。

    圖3 超聲速分離器幾何模型Fig.3 Geometrical model of supersonic gas separator

    采用Gambit對模型進(jìn)行網(wǎng)格劃分,如圖4所示,圖4(a)、(b)分別為噴管喉部和排液口的網(wǎng)格示意圖。為減小計算量,模型主要采用四邊形結(jié)構(gòu)化網(wǎng)格,近壁面采用邊界層網(wǎng)格,分別對噴管喉部和排液口附近區(qū)域進(jìn)行網(wǎng)格細(xì)化。為了避免網(wǎng)格疏密程度對計算結(jié)果造成影響,對表3中的不同網(wǎng)格數(shù)模型分別進(jìn)行計算,并對比和分析不同網(wǎng)格數(shù)量下Mach數(shù)沿軸向分布曲線。

    圖4 二維數(shù)值計算網(wǎng)格劃分Fig.4 Grid of two-dimensional numerical computation

    表3 不同網(wǎng)格數(shù)量Table 3 Different computational grids

    圖5為壓比1.4,面積比1.27時,不同網(wǎng)格數(shù)量下的Mach數(shù)沿中心軸線的分布曲線。從圖中可以看出,模型2和模型3的分布曲線基本一致,即當(dāng)最小網(wǎng)格尺寸取0.01 mm時,繼續(xù)細(xì)化網(wǎng)格對計算結(jié)果沒什么影響。因此,采用模型網(wǎng)格數(shù)量為46235,近壁面邊界層網(wǎng)格最小厚度為0.01 mm。

    圖5 不同網(wǎng)格數(shù)下Mach數(shù)沿軸向的變化曲線Fig.5 Curves of Mach number distribution along axis with different computational grids

    利用計算流體力學(xué)CFD軟件Fluent對建立的數(shù)值模型進(jìn)行求解,采用基于密度定常隱式求解器,在空間上采用二階迎風(fēng)有限體積格式對控制進(jìn)行離散。進(jìn)口及濕、干氣出口均采用壓力邊界條件,如圖3(a)所示,對稱軸給定軸對稱邊界條件,壁面采用無滑移邊界條件,壓力進(jìn)口邊界條件給定總壓和總溫,壓力出口邊界條件給定靜壓值,同時出口溫度采用第2類邊界條件。用于建模的結(jié)構(gòu)參數(shù):L1=100 mm,L2=500 mm,Lb=1.5 mm,inner angle=8°,outer angle=12°~24°,A2/A1=1.27~1.35;數(shù)值計算的操作參數(shù)入口總壓pin=0.12、0.13、0.14、0.15 MPa,入口總溫Tin=300 K,干、濕氣出口的靜壓pdry、pwet均為0.1 MPa。操作參數(shù)和結(jié)構(gòu)參數(shù)的確定遵循以下流程:根據(jù)實驗室所能提供的流量條件和最大壓比確定噴管的喉部直徑→根據(jù)噴管的理論計算公式確定噴管進(jìn)出口截面積→利用數(shù)值計算軟件Fluent來驗證設(shè)計的合理性→實驗研究→確定最優(yōu)結(jié)構(gòu)參數(shù)和操作參數(shù)。

    3.2計算結(jié)果及分析

    3.2.1結(jié)構(gòu)和操作參數(shù)對流動的影響超聲速旋流分離器噴管擴張段內(nèi)軸線上的Mach數(shù)分布是反映噴管內(nèi)流場特性的一個重要因素?;诒疚慕⒌亩S軸對稱模型,研究了不同排液口傾角、面積比及壓比下噴管內(nèi)Mach數(shù)沿軸向的分布情況。

    圖6(a)是在壓比為1.4,排液通道外傾角為12°,噴管面積比分別為1.27、1.29、1.34和1.35時Mach數(shù)沿軸向的分布曲線,從圖中可以看出,不同面積比下的計算結(jié)果有著相同的變化趨勢,當(dāng)氣體運動到噴管喉部時達(dá)到聲速,在噴管漸擴段內(nèi),氣體繼續(xù)膨脹加速,產(chǎn)生低溫環(huán)境,使含濕氣體在噴管內(nèi)凝結(jié);其后在排液口附近以及擴壓段內(nèi)產(chǎn)生激波,氣體速度降低,靜壓升高。顯然,隨著面積比的增大,噴管漸擴段內(nèi)相同位置處的Mach數(shù)也逐漸增大。

    為了探究在排液口附近產(chǎn)生激波的原因,對排液通道外傾角以及壓比對流動狀態(tài)的影響進(jìn)行研究,圖6(b)描述的是在壓比為1.4,面積比為1.27,排液通道外傾角分別為12°、16°、20°和24°時Mach數(shù)沿軸向的分布曲線,從圖中可以看出,隨著排液通道外傾角的增大,在排液口附近形成的激波逐漸向噴管上游移動,當(dāng)外傾角為12°時,排液口對噴管內(nèi)的流動特性的影響最小。

    圖6 不同參數(shù)下Mach數(shù)沿軸向的變化曲線Fig.6 Curves of Mach number distribution along axis with different parameters

    圖6(c)描述了面積比為1.27,排液通道外傾角為12°,壓比分別為1.2、1.3、1.4和1.5時Mach數(shù)沿軸向的分布曲線,從圖中可以看出,隨著壓比的增大,噴管內(nèi)激波產(chǎn)生的位置向下游移動,但是當(dāng)壓比大于1.4,繼續(xù)增大壓比對激波所產(chǎn)生的位置影響不大;從圖中還可以明顯看出,隨著壓比的增大,擴壓段內(nèi)激波前的Mach數(shù)越大,即總壓損失越大,所以在此面積比條件下,采用1.4的壓比即可。

    3.2.2原因分析由以上的分析可知,結(jié)構(gòu)參數(shù)和操作參數(shù)的改變并不能避免排液口附近激波的產(chǎn)生,然而激波的存在會導(dǎo)致排液口附近溫度升高,凝結(jié)的液滴可能會再次揮發(fā)與干氣混合,降低分離效率。因此,本節(jié)以壓比1.4,面積比1.27,排液通道外傾角12°為例進(jìn)行分析。圖7為噴管出口與排液口附近靠近壁面處軸向速度沿軸向的分布曲線,從圖中可以看出,在x = 598 mm處,軸向速度由正向變?yōu)榉聪?,且隨著軸向距離的增加,軸線速度的大小也顯著增大,在x = 602 mm處,軸向速度恢復(fù)成向下游移動,產(chǎn)生這種流動現(xiàn)象的原因是排液口壓力邊界條件設(shè)置的壓力為0.10 MPa,噴管內(nèi)的高速低壓氣體撞擊排液腔內(nèi)的高壓氣體,形成一道反向壓縮波,并向上游移動,引起返流現(xiàn)象。由于噴管內(nèi)是超聲速,反向壓縮波并不能向上游移動,但是在靠近壁面處流速低且速度梯度較大,反向壓縮波的存在會對靠近壁面的流速產(chǎn)生較大的影響,使靠近壁面的流速降低,增加了邊界層的厚度,邊界層的存在使噴管出口附近實際截面積減小,從而使氣體不能繼續(xù)膨脹加速,流速降低。

    圖7 噴管出口與排液口入口附近靠近壁面處軸向速度沿軸向的分布曲線Fig.7 Curves of axial velocity close to wall along axis of nozzle outlet and liquid discharge inlet

    圖8為噴管出口、排液口和漸擴器入口相接處的速度等值線分布云圖,從圖中可以看出在x = 587 mm處壁面邊界層逐漸變厚,對噴管內(nèi)流速的影響也越來越明顯,使噴管出口附近流速由400 m·s?1降低至300 m·s?1。因此,為了避免回流對分離效率的影響,在實驗的過程中可通過調(diào)節(jié)干氣出口閥門控制干氣出口背壓,在保證壓比的條件下提高噴管內(nèi)的靜壓,使?jié)駳獬隹谟袣怏w流出。

    圖8 噴管出口、排液口和漸擴器入口相接處的速度等值線分布云圖Fig.8 Velocity counter map of nozzle outlet, inlet of liquid discharge and diffuser

    4 實驗研究

    4.1SGS實驗裝置結(jié)構(gòu)

    導(dǎo)流式超聲速冷凝旋流分離器作為整個實驗流程的核心部件,其整體結(jié)構(gòu)如圖9所示,噴管和排液腔的結(jié)構(gòu)尺寸與前文用于二維數(shù)值計算的尺寸一致。采用軸流式旋流器起旋,混合氣體沿著軸線方向進(jìn)入旋流器中,沿著葉片流道流動而形成強旋流。葉片的成型原理[23]是一條直線與一個圓柱面相交成一定的角度,直線沿著圓柱面上的某條圓弧線移動所形成的曲面即葉片型面。圖10為旋流發(fā)生器的結(jié)構(gòu)示意圖,本文實驗研究采用的旋流發(fā)生器高度為160 mm,葉片出口角為55°。

    圖9 SGS整體結(jié)構(gòu)Fig.9 Whole structural diagram of SGS

    圖10 旋流發(fā)生器的結(jié)構(gòu)Fig.10 Structural diagram of swirling generator

    4.2實驗流程

    圖11為超聲速旋流分離裝置分離性能研究的實驗流程。穩(wěn)壓罐出來的高壓空氣與蒸氣發(fā)生裝置提供的酒精蒸氣混合,經(jīng)過立式氣液旋流分離器進(jìn)行預(yù)分離。經(jīng)過預(yù)分離后的混合氣體進(jìn)入過濾裝置進(jìn)行再次分離,得到的空氣和乙醇蒸氣的混合氣體作為實驗介質(zhì)?;旌蠚怏w進(jìn)入導(dǎo)流式超聲速旋流分離器,氣流經(jīng)設(shè)備處理后被分為兩股氣流,濕度較小的稱為干氣從干氣出口6排出,另一股稱為濕氣,從排液口7排出。

    圖11 超聲速旋流分離裝置性能實驗研究流程Fig.11 Schematic diagram of experimental process of supersonic gas separator1—low pressure air; 2—high pressure air; 3—liquid water (or ethanol); 4—vapor(water or ethanol); 5—mixture of air and vapor; 6—dry gas; 7—wet gas

    4.3性能評價指標(biāo)

    為了評價裝置的分離性能,研究排液裝置結(jié)構(gòu)的影響,以重組分脫除率、露點降和干氣率作為裝置分離性能評價指標(biāo)。

    重組分脫除率δsteam的計算公式為

    式中,xsteam,in和xsteam,dry分別為超聲速分離設(shè)備進(jìn)口和干氣出口混合氣體所含乙醇的質(zhì)量分?jǐn)?shù),由色譜分析儀測得。

    露點降ΔTd的計算公式如下

    式中,Td,in和Td,dry分別為超聲速分離設(shè)備進(jìn)口和干氣出口混合氣體的露點,可以由露點儀直接測得。

    4.4實驗結(jié)果及分析

    在實驗中,通過激光筆照射透明玻璃管內(nèi),可以看到垂直光柱(圖12),這一典型的丁達(dá)爾現(xiàn)象說明噴管旋流氣內(nèi)含有冷凝液滴。通過肉眼也可以觀測到整個透明管內(nèi)液滴的跡線(圖13),混合氣體以與水平方向成固定角度的螺旋線軌跡運動?;旌蠚怏w通過旋流器獲得切向速度,在進(jìn)氣與出氣壓差的作用下產(chǎn)生軸向運動速度,因此混合氣體以螺旋線的軌跡運動。在透明有機玻璃管壁面處可觀測到液體沿壁面螺旋線運動,即冷凝的液滴在離心力作用下已到達(dá)壁面處,可以從排液口處分離。

    圖12 激光照射圖Fig.12 Picture of laser irradiating nozzle

    圖13 液相跡線圖Fig.13 Track diagram of liquid phase

    圖14為壓比1.4時,不同面積比下混合氣體的露點降和乙醇脫除率隨面積比的變化曲線,顯然,不同面積比條件下,乙醇脫除率與露點降的變化趨勢基本相同,均隨著面積比的增大而減??;當(dāng)面積比為1.27時,乙醇脫除率和干氣露點降均達(dá)到最大值,分別為20.5%和4.2 K。

    圖14 乙醇露點降和脫除率隨面積比的變化Fig.14 Variation of ΔTdand δsteamwith AR

    根據(jù)前面模擬結(jié)果可知,當(dāng)壓比為1.4,排液腔外傾角為12°時,不同面積比下,噴管內(nèi)均能達(dá)到超聲速流動狀態(tài),且隨著面積比的增大,噴管內(nèi)Mach數(shù)逐漸增大,然而分離效率卻呈現(xiàn)降低的趨勢,其原因是Mach數(shù)增大,即軸向速度增加,液滴在噴管內(nèi)的停留時間降低,已經(jīng)凝結(jié)的液滴未甩到壁面上就已經(jīng)從干氣出口排出,乙醇脫除率和干氣露點降表現(xiàn)出下降的趨勢,因此,在設(shè)計超聲速噴管時應(yīng)該選擇合適的面積比,在保證噴管內(nèi)的低溫環(huán)境足以實現(xiàn)可凝蒸氣的凝結(jié),且保證液滴在噴管內(nèi)的停留時間。

    5 結(jié) 論

    通過對導(dǎo)流錐式超聲速旋流分離器進(jìn)行數(shù)值和實驗研究,得出如下結(jié)論。

    (1)減小排液通道的外傾角可以有效改善返流的影響,避免激波在噴管內(nèi)產(chǎn)生,但是適當(dāng)增加面積比和壓比并不能消除回流現(xiàn)象,在實驗過程中應(yīng)調(diào)節(jié)裝置干氣出口閥,在保證壓比的前提下,提高背壓使?jié)駳饽軓呐乓呵慌懦觥?/p>

    (2)搭建了超聲速旋流分離器的實驗平臺,在實驗中發(fā)現(xiàn),可通過加工透明玻璃管和激光照射的方法判斷噴管內(nèi)可凝組分是否已經(jīng)開始凝結(jié)。

    (3)在壓比一定條件下,適當(dāng)?shù)卦龃竺娣e比即是為含濕氣體的凝結(jié)提供了更低溫的環(huán)境,但是在面積比為1.27時,噴管內(nèi)的低溫環(huán)境已可以保證含濕氣體的凝結(jié),繼續(xù)增大面積比即加大了氣體在噴管內(nèi)膨脹的力度使軸向速度增大縮短液滴沉降時間從而影響分離效率。在壓比為1.4,面積比為1.27時,乙醇脫除率達(dá)到20.5%。

    符號說明

    A1, A2——分別為噴管喉部截面積、出口截面積,mm2

    cp——流體的比定壓熱容,J·kg?1·K?1

    L1, L2——分別為噴管漸縮段和漸擴段的長度,mm

    PRC——壓力恢復(fù)系數(shù)

    pin, pdry, pwet——分別為噴管進(jìn)口、干氣出口和濕氣出口的壓力,MPa

    Td,in, Td,dry——分別為超聲速分離設(shè)備進(jìn)口和干氣出口混合氣體的露點,K

    ΔTd——混合氣體露點降,K

    Tin, Tdry, Twet——分別為噴管進(jìn)口、干氣出口和濕氣出口的溫度,K

    u1, u2——分別為流場x、y方向的速度,m· s?1

    x——噴管任一截面與進(jìn)口的軸向距離,mm

    xsteam,in, xsteam,dry——分別為超聲速分離設(shè)備進(jìn)口和干氣出口混合氣體所含乙醇的質(zhì)量分?jǐn)?shù)

    δsteam——混合氣體重組分脫除率

    κ——流體的傳熱系數(shù),W·m?2·K?1

    ρ——流體密度,kg·m?3

    ρgi——i方向上的重力體積力,N·m?3

    References

    [1]吳小林, 熊至宜, 姬忠禮. 天然氣凈化用旋風(fēng)分離器氣液分離性能[J]. 化工學(xué)報, 2010, 61(9): 2430-2436. WU X L, XIONG Z Y, JI Z L. Gas-liquid separation performance of cyclone separator for purification of natural gas [J]. CIESC Journal, 2010, 61(9): 2430-2436.

    [2]何策, 程雁. 天然氣超音速脫水技術(shù)評析[J]. 石油機械, 2006, 34(5): 70-72. HE C, CHENG Y. Analysis of supersonic dehydration technology of natural gas [J]. China Petroleum Machinery, 2006, 34(5): 70-72.

    [3]胡耀強, 何飛, 韓建紅. 天然氣脫水技術(shù)[J]. 化學(xué)工程與裝備, 2013, (3): 151-153. HU Y Q, HE F, HAN J H. Natural gas dehydration technology [J]. Chemical Engineering & Equipment, 2013, (3): 151-153.

    [4]王旭鋒, 逯國英, 張萍, 等. 鹽穴儲氣庫乙二醇與JT閥制冷脫水工藝設(shè)計[J]. 中國新技術(shù)新產(chǎn)品, 2014, (23): 59-60. WANG X F, LU G Y, ZHANG P, et al. Design of JT valve refrigeration dehydration process and salt cavity gas storage ethylene glycol [J]. China New Technologies and New Products, 2014, (23): 59-60.

    [5]蔣洪, 楊昌平, 吳敏, 等. 天然氣三甘醇脫水裝置節(jié)能分析[J]. 石油與天然氣化工, 2010, 39(2): 122-127. JIANG H, YANG C P, WU M, et al. Analysis on energy saving of three glycol dehydration unit of natural gas [J]. Chemical Engineering of Oil & Gas, 2010, 39(2): 122-127.

    [6]KUZNICKI S M, BELL V A, NAIR S. A titanosilicate molecular sieve with adjustable pores for size-selective adsorption of molecules [J]. Nature, 2001, 412(6848): 720-724.

    [7]高曉根, 計維安, 劉薔, 等. 超音速分離技術(shù)及在氣田地面工程中的應(yīng)用[J]. 石油與天然氣化工, 2011, 40(1): 42-46. GAO X G, JI W A, LIU Q, et al. The supersonic vortex tube separation technology and its application in the gas-field [J]. Chemical Engineering of Oil & Gas, 2011, 40(1): 42-46.

    [8]BETTING M M, HOLTEN T V, et al. Supersonic separator apparatus and method: US 6776825B2 [P]. 2004-08-17.

    [9]SCHINKELSHOEK P, EPSOM H. Supersonic gas conditioning-low pressure drop twister for NGL recovery [C]// Offshore Technology Conference, 2006.

    [10]楊勇, 張琨, 宋煜. 水蒸氣超音速非平衡流動的變壓凝結(jié)特性[J].化工學(xué)報, 2010, 61(4): 820-824. YANG Y, ZHANG K, SONG Y. Condensation characteristics of supersonic non-equilibrium steam flow under different pressures [J]. CIESC Journal, 2010, 61(4): 820-824.

    [11]涂輝, 蔣洪, 劉曉強. 超音速分離在天然氣脫水中的應(yīng)用[J]. 管道技術(shù)與設(shè)備, 2008, (3): 1-3. TU H, JIANG H, LIU X Q. Application of supersonic separator to natural gas dehydration[J]. Pipeline Technique and Equipment, 2008, (3): 1-3.

    [12]Knott T, 高衛(wèi)東, 紀(jì)常杰. 新型超聲波分離裝置[J]. 國外油田工程, 2001, 17(2): 34-35. KNOTT T, GAO W D, JI C J. New-type supersonic separation device [J]. Foreign Oilfield Engineering, 2001, 17(2): 34-35.

    [13]張宏義, 張發(fā)啟, 孫權(quán). 二維超音速噴管設(shè)計仿真研究[J]. 機械設(shè)計與制造, 2012, (10): 132-134. ZHANG H Y, ZHANG F Q, SUN Q. Designing and numerical simulation on two dimensions supersonic nozzle [J]. Machinery Design & Manufacture, 2012, (10): 132-134.

    [14]曹學(xué)文, 陳麗, 林宗虎, 等. 用于超聲速旋流分離器中的超聲速噴管研究[J]. 天然氣工業(yè), 2007, 27(7): 112-114. CAO X W, CHEN L, LIN Z H, et al. The nozzle used in supersonic swirling separators [J]. Natural Gas Industry, 2007, 27(7): 112-114.

    [15]楊文, 曹雪文. Laval 噴管設(shè)計及在天然氣液化中的應(yīng)用研究[J].西安石油大學(xué)學(xué)報(自然科學(xué)版), 2015, 30(2): 75-79. YANG W, CAO X W. Design Laval nozzle and its application in natural gas liquidation [J]. Journal of Xi’an Shiyou University (Natural Science Edition), 2015, 30(2): 75-79.

    [16]WYSLOUZIL B E, HEATH C H, CHEUNG J L, et al. Binary condensation in a supersonic nozzle [J]. Journal of Chemical Physics, 2000, 113(17): 7317-7329.

    [17]WEN C, CAO X W, YANG Y, et al. Swirling effects on the performance of supersonic separators for natural gas separation [J]. Chemical Engineering & Technology, 2011, 34(9): 1575-1580.

    [18]WEN C, CAO X W, YANG Y. Swirling flow of natural gas in supersonic separators [J]. Chemical Engineering and Processing: Process Intensification, 2011, 50(7): 644-649.

    [19]YANG Y, WEN C, WANG S L, et al. Theoretical and numerical analysis on pressure recovery of supersonic separators for natural gas dehydration [J]. Applied Energy, 2014, 132: 248-253.

    [20]VAN LEER B. Upwind-difference methods for aerodynamics problems governed by the Euler equations of gas dynamics [J]. Lectures in Applied Mathematics, 1985, (22): 327-336.

    [21]YANG Y, CAO W, WANG S, et al. Numerical simulation of real gas flows in natural gas supersonic separation processing [J]. Journal of Natural Gas Science & Engineering, 2014, 21: 829-836.

    [22]蔣文明, 劉中良, 劉曉麗, 等. 濕空氣超音速凝結(jié)特性[J].化工學(xué)報, 2011, 62(S1): 97-102. JIANG W M, LIU Z L, LIU X L, et al. Experimental study on wet air in supersonic condensation flows [J]. CIESC Journal, 2011, 62(S1): 97-102.

    [23]金有海, 范超. 導(dǎo)葉式旋風(fēng)管葉片參數(shù)設(shè)計方法的研究[J]. 化工機械, 1999, 26(1): 21-24. JIN Y H, FAN C. Research on vane design method of guide vane cyclone [J]. Chemical Machinery, 1999, 26(1): 21-24.

    Flow characteristic of supersonic gas separator with diversion cone

    HU Dapeng, WANG Yingguang, REN Wenwen, ZHAO Jianhua, LIU Peiqi
    (School of Chemical Machinery, Dalian University of Technology, Dalian 116023, Liaoning, China)

    Abstract:The two-dimensional axial model is adopted to study the flow characteristic in supersonic nozzle using ideal gas as medium. Experimental platform is set up to study the influence of structural and operating parameters on flow field and separation efficiency under low pressure ratio. The numerical results indicate that there is no shock wave in the nozzle divergent section when pressure ratio is 1.4, outer angles of liquid discharge chamber is less than 12° and area ratio is 1.27. Meanwhile, the analysis of the leading shock waves indicates that the reverse compression waves appear in the liquid discharge chamber and move to nozzle upstream. The experimental results show that the highest separation efficiency of the supersonic gas separator (SGS) can reach 20.5% when the pressure ratio and area ratio are set to 1.4 and 1.27, respectively.

    Key words:supersonic flow; hydrodynamics; numerical simulation; shock wave; centrifugation; Mach number

    中圖分類號:TQ 051.8

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

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

    DOI:10.11949/j.issn.0438-1157.20151817

    基金項目:國家自然科學(xué)基金項目(21206013,21476036);遼寧省教育廳重點實驗室基礎(chǔ)研究項目(LZ2015019)。

    Corresponding author:Prof. LIU Peiqi, lpq21cn@dlut.edu.cn

    欧美精品啪啪一区二区三区| 嫩草影视91久久| 一区二区日韩欧美中文字幕| 亚洲精品中文字幕一二三四区| 欧美乱妇无乱码| 亚洲av电影不卡..在线观看| 亚洲一区中文字幕在线| 亚洲片人在线观看| 日韩av在线大香蕉| 一进一出抽搐动态| 午夜免费鲁丝| 亚洲电影在线观看av| 国产xxxxx性猛交| 午夜福利免费观看在线| 老司机午夜十八禁免费视频| 久久精品国产综合久久久| 国产成人欧美| 99在线人妻在线中文字幕| 日韩高清综合在线| 欧美久久黑人一区二区| 国产一级毛片七仙女欲春2 | 老汉色av国产亚洲站长工具| 久久久国产成人免费| 久久午夜综合久久蜜桃| 一区二区三区激情视频| 国产成年人精品一区二区| av超薄肉色丝袜交足视频| 女人被狂操c到高潮| 精品国产国语对白av| ponron亚洲| 亚洲久久久国产精品| 这个男人来自地球电影免费观看| 亚洲欧美精品综合久久99| 91九色精品人成在线观看| 午夜两性在线视频| 757午夜福利合集在线观看| 久久性视频一级片| 满18在线观看网站| 视频区欧美日本亚洲| 亚洲第一青青草原| 中文字幕人成人乱码亚洲影| 美女大奶头视频| 欧美黑人精品巨大| 亚洲av美国av| 91成年电影在线观看| 欧美成人性av电影在线观看| a在线观看视频网站| 在线观看日韩欧美| 自拍欧美九色日韩亚洲蝌蚪91| 国产av一区在线观看免费| 亚洲欧美日韩高清在线视频| 精品国产国语对白av| 欧美精品亚洲一区二区| 国产主播在线观看一区二区| 69av精品久久久久久| 亚洲一区高清亚洲精品| 色老头精品视频在线观看| 91成人精品电影| 亚洲无线在线观看| 久久精品国产99精品国产亚洲性色 | 欧美国产日韩亚洲一区| 一级黄色大片毛片| 久久精品国产清高在天天线| 国产精华一区二区三区| 国内精品久久久久精免费| 国产精品1区2区在线观看.| 成人三级黄色视频| 国产精品精品国产色婷婷| 女警被强在线播放| 色播在线永久视频| 男女午夜视频在线观看| 欧美成人性av电影在线观看| 熟女少妇亚洲综合色aaa.| 韩国av一区二区三区四区| 亚洲少妇的诱惑av| 一夜夜www| 女人被狂操c到高潮| 久久热在线av| 怎么达到女性高潮| 国产精品一区二区在线不卡| 国产精品野战在线观看| 可以在线观看毛片的网站| 精品国产乱子伦一区二区三区| 久久午夜综合久久蜜桃| 国内久久婷婷六月综合欲色啪| 久久国产精品人妻蜜桃| 中文字幕高清在线视频| 色老头精品视频在线观看| 亚洲精品国产区一区二| 人人妻人人澡欧美一区二区 | 国产精品亚洲av一区麻豆| 国产人伦9x9x在线观看| 免费一级毛片在线播放高清视频 | 少妇裸体淫交视频免费看高清 | 久久精品亚洲熟妇少妇任你| 91九色精品人成在线观看| 真人一进一出gif抽搐免费| 亚洲色图 男人天堂 中文字幕| 999精品在线视频| 国产精品综合久久久久久久免费 | 成人免费观看视频高清| 色播亚洲综合网| 中文字幕人妻熟女乱码| 熟女少妇亚洲综合色aaa.| 精品久久久久久久人妻蜜臀av | 99热只有精品国产| 男女之事视频高清在线观看| 91国产中文字幕| 国产成人欧美在线观看| 欧美一级a爱片免费观看看 | 国产精品美女特级片免费视频播放器 | 在线十欧美十亚洲十日本专区| 夜夜躁狠狠躁天天躁| 成人精品一区二区免费| 久久热在线av| 久久精品国产亚洲av高清一级| 满18在线观看网站| 一级作爱视频免费观看| 在线播放国产精品三级| 可以在线观看毛片的网站| 精品福利观看| 欧美+亚洲+日韩+国产| avwww免费| 国产亚洲欧美在线一区二区| 国产精品 国内视频| 国产成人av教育| 伦理电影免费视频| 欧美乱色亚洲激情| 侵犯人妻中文字幕一二三四区| 国产精品影院久久| 免费在线观看完整版高清| 亚洲性夜色夜夜综合| 制服诱惑二区| 老司机福利观看| 久久国产亚洲av麻豆专区| 九色亚洲精品在线播放| av天堂在线播放| 99热只有精品国产| 久久久久久久精品吃奶| 女人精品久久久久毛片| 国产精品野战在线观看| 不卡一级毛片| 一进一出好大好爽视频| 午夜亚洲福利在线播放| 无人区码免费观看不卡| 91老司机精品| 亚洲欧美精品综合一区二区三区| 亚洲自拍偷在线| 久久久久久国产a免费观看| 免费看十八禁软件| 亚洲av电影不卡..在线观看| 久久久久久久久中文| 999久久久精品免费观看国产| 男女下面进入的视频免费午夜 | www.自偷自拍.com| 午夜福利,免费看| 亚洲国产欧美一区二区综合| 十分钟在线观看高清视频www| 久久 成人 亚洲| 一区二区三区高清视频在线| 久久精品国产清高在天天线| 满18在线观看网站| 美女 人体艺术 gogo| 精品一区二区三区av网在线观看| 亚洲午夜精品一区,二区,三区| 丁香欧美五月| 成在线人永久免费视频| 波多野结衣一区麻豆| 中文字幕最新亚洲高清| 日韩国内少妇激情av| 神马国产精品三级电影在线观看 | 国产亚洲av高清不卡| avwww免费| 欧美成人性av电影在线观看| 欧美丝袜亚洲另类 | 国产97色在线日韩免费| 脱女人内裤的视频| 亚洲三区欧美一区| 在线观看免费视频网站a站| 18禁美女被吸乳视频| 女人精品久久久久毛片| 午夜日韩欧美国产| 欧美亚洲日本最大视频资源| 日韩欧美一区二区三区在线观看| 51午夜福利影视在线观看| 性少妇av在线| 精品久久久久久久人妻蜜臀av | 成人亚洲精品一区在线观看| 久9热在线精品视频| 亚洲欧美精品综合久久99| 午夜日韩欧美国产| 国产高清视频在线播放一区| 久久中文看片网| 国产精品久久视频播放| 国产亚洲欧美精品永久| 国产免费男女视频| 国产野战对白在线观看| 久久久国产精品麻豆| 麻豆国产av国片精品| 制服人妻中文乱码| 一二三四在线观看免费中文在| 97人妻精品一区二区三区麻豆 | 亚洲中文字幕一区二区三区有码在线看 | 国产欧美日韩一区二区精品| 亚洲av电影在线进入| 69精品国产乱码久久久| 国产精品一区二区免费欧美| 亚洲一区中文字幕在线| 97人妻精品一区二区三区麻豆 | 亚洲av美国av| 18禁裸乳无遮挡免费网站照片 | 老司机午夜十八禁免费视频| 国产1区2区3区精品| 国产高清视频在线播放一区| 长腿黑丝高跟| 欧美乱妇无乱码| 精品国产乱子伦一区二区三区| 国产99白浆流出| 亚洲精品粉嫩美女一区| 亚洲一卡2卡3卡4卡5卡精品中文| 在线av久久热| 成人18禁高潮啪啪吃奶动态图| 亚洲国产欧美网| 九色亚洲精品在线播放| 在线观看免费视频网站a站| 夜夜爽天天搞| 国产人伦9x9x在线观看| 九色亚洲精品在线播放| 国产精品久久久久久精品电影 | 亚洲人成电影免费在线| 日日摸夜夜添夜夜添小说| 日本免费一区二区三区高清不卡 | 国产1区2区3区精品| 久久人人97超碰香蕉20202| 国产三级黄色录像| 日韩有码中文字幕| 免费一级毛片在线播放高清视频 | 午夜免费激情av| 日韩欧美在线二视频| 国语自产精品视频在线第100页| 亚洲精品久久成人aⅴ小说| a级毛片在线看网站| 丰满人妻熟妇乱又伦精品不卡| 精品福利观看| 久久精品aⅴ一区二区三区四区| 色播亚洲综合网| 日韩欧美在线二视频| 精品不卡国产一区二区三区| 亚洲,欧美精品.| 青草久久国产| 可以在线观看毛片的网站| 亚洲精品一区av在线观看| 欧美激情极品国产一区二区三区| 可以在线观看的亚洲视频| 久久国产精品人妻蜜桃| 动漫黄色视频在线观看| 美女高潮喷水抽搐中文字幕| 好男人在线观看高清免费视频 | 亚洲自偷自拍图片 自拍| 亚洲国产精品999在线| 中文字幕另类日韩欧美亚洲嫩草| 亚洲成国产人片在线观看| 咕卡用的链子| 久久久久亚洲av毛片大全| 制服丝袜大香蕉在线| 午夜久久久在线观看| 女人被躁到高潮嗷嗷叫费观| 亚洲av第一区精品v没综合| 久久草成人影院| 久久国产亚洲av麻豆专区| 亚洲电影在线观看av| 午夜免费观看网址| 男人的好看免费观看在线视频 | 色综合站精品国产| bbb黄色大片| 男男h啪啪无遮挡| 纯流量卡能插随身wifi吗| 亚洲中文av在线| 在线永久观看黄色视频| or卡值多少钱| 欧美中文日本在线观看视频| 欧美激情久久久久久爽电影 | 91大片在线观看| 久久久久久人人人人人| 亚洲精品美女久久av网站| 动漫黄色视频在线观看| 天堂√8在线中文| 午夜精品在线福利| 91老司机精品| 悠悠久久av| 午夜福利免费观看在线| 亚洲一区高清亚洲精品| 亚洲色图av天堂| 免费少妇av软件| 91av网站免费观看| 亚洲成国产人片在线观看| 色婷婷久久久亚洲欧美| 一级a爱片免费观看的视频| 亚洲一码二码三码区别大吗| 一卡2卡三卡四卡精品乱码亚洲| 国产精品99久久99久久久不卡| 色综合亚洲欧美另类图片| 高清毛片免费观看视频网站| 国内精品久久久久精免费| 国产成人欧美| 久久香蕉国产精品| 国内精品久久久久久久电影| 亚洲av成人一区二区三| 国产精品久久久av美女十八| 一边摸一边抽搐一进一小说| 久久人妻av系列| 自线自在国产av| 视频在线观看一区二区三区| 日韩国内少妇激情av| av有码第一页| 亚洲欧美日韩高清在线视频| 亚洲成人国产一区在线观看| 午夜成年电影在线免费观看| 丁香欧美五月| 久久这里只有精品19| 日本黄色视频三级网站网址| 亚洲五月天丁香| 国产亚洲精品一区二区www| 国产成人欧美在线观看| 精品人妻1区二区| 好男人在线观看高清免费视频 | 老司机靠b影院| 日韩欧美免费精品| 看黄色毛片网站| 欧美黄色片欧美黄色片| 国产aⅴ精品一区二区三区波| 搡老岳熟女国产| 性少妇av在线| 天天躁夜夜躁狠狠躁躁| 好男人在线观看高清免费视频 | 欧美激情 高清一区二区三区| 欧美中文综合在线视频| 乱人伦中国视频| 麻豆国产av国片精品| 久久国产精品影院| 国产精品 国内视频| 啦啦啦韩国在线观看视频| 亚洲精品在线观看二区| 超碰成人久久| 免费看a级黄色片| 亚洲aⅴ乱码一区二区在线播放 | 国产免费av片在线观看野外av| 99精品久久久久人妻精品| 精品久久久精品久久久| 亚洲精品久久成人aⅴ小说| 久9热在线精品视频| 国产精品亚洲美女久久久| 亚洲成人国产一区在线观看| 免费在线观看亚洲国产| 在线观看免费视频日本深夜| 亚洲中文av在线| 成人亚洲精品av一区二区| 身体一侧抽搐| 少妇裸体淫交视频免费看高清 | 此物有八面人人有两片| 免费人成视频x8x8入口观看| 久久久久久国产a免费观看| 两人在一起打扑克的视频| 欧美不卡视频在线免费观看 | 九色国产91popny在线| 亚洲国产欧美日韩在线播放| 国产一区二区三区在线臀色熟女| 桃色一区二区三区在线观看| 国产精品98久久久久久宅男小说| 在线天堂中文资源库| 精品一区二区三区四区五区乱码| 免费无遮挡裸体视频| 19禁男女啪啪无遮挡网站| 精品午夜福利视频在线观看一区| 黑人巨大精品欧美一区二区mp4| 亚洲人成77777在线视频| 美女扒开内裤让男人捅视频| 中文字幕精品免费在线观看视频| 亚洲成人国产一区在线观看| 日韩大码丰满熟妇| 国产亚洲精品av在线| 亚洲激情在线av| 三级毛片av免费| 国产亚洲av高清不卡| 国产精品九九99| 欧美另类亚洲清纯唯美| 欧美亚洲日本最大视频资源| 免费高清视频大片| 很黄的视频免费| 女同久久另类99精品国产91| 亚洲人成77777在线视频| 黄色丝袜av网址大全| 久久人人精品亚洲av| 桃红色精品国产亚洲av| 亚洲午夜理论影院| 国产精品久久久人人做人人爽| 欧美乱妇无乱码| 亚洲一区二区三区不卡视频| 亚洲第一青青草原| 女人精品久久久久毛片| 成人特级黄色片久久久久久久| 免费看美女性在线毛片视频| 国产在线精品亚洲第一网站| 久久国产乱子伦精品免费另类| 麻豆av在线久日| 精品国产国语对白av| 黄片小视频在线播放| 亚洲免费av在线视频| 身体一侧抽搐| 午夜福利免费观看在线| 级片在线观看| 日本五十路高清| 91国产中文字幕| 麻豆一二三区av精品| 国产精品久久久久久精品电影 | 性欧美人与动物交配| 91九色精品人成在线观看| 在线免费观看的www视频| 欧美中文日本在线观看视频| 午夜福利高清视频| 色精品久久人妻99蜜桃| 美国免费a级毛片| 国产欧美日韩精品亚洲av| 婷婷丁香在线五月| 国产精品综合久久久久久久免费 | 国产单亲对白刺激| 国产精品综合久久久久久久免费 | 中文字幕av电影在线播放| 欧美老熟妇乱子伦牲交| 精品久久久久久成人av| 男女午夜视频在线观看| 一级作爱视频免费观看| 黄网站色视频无遮挡免费观看| 高清黄色对白视频在线免费看| 久久 成人 亚洲| 18禁国产床啪视频网站| 精品久久蜜臀av无| 久久欧美精品欧美久久欧美| 18禁裸乳无遮挡免费网站照片 | 日韩大尺度精品在线看网址 | 欧美丝袜亚洲另类 | 日韩精品中文字幕看吧| 99在线人妻在线中文字幕| 午夜福利成人在线免费观看| 国产99白浆流出| cao死你这个sao货| 91在线观看av| 一二三四社区在线视频社区8| 美女扒开内裤让男人捅视频| 成人国产一区最新在线观看| 18美女黄网站色大片免费观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产精品合色在线| 91成年电影在线观看| 久久久国产成人精品二区| 精品一区二区三区四区五区乱码| 最好的美女福利视频网| 无人区码免费观看不卡| 又黄又粗又硬又大视频| 国产亚洲欧美在线一区二区| 欧美一级毛片孕妇| 超碰成人久久| av有码第一页| 亚洲免费av在线视频| 日韩国内少妇激情av| 久久亚洲精品不卡| 久久久精品国产亚洲av高清涩受| 久久精品亚洲熟妇少妇任你| 国产精品久久久久久精品电影 | 88av欧美| 人人妻人人澡人人看| 欧美中文日本在线观看视频| 看黄色毛片网站| 国产真人三级小视频在线观看| 丝袜人妻中文字幕| 国产1区2区3区精品| 亚洲久久久国产精品| 欧美丝袜亚洲另类 | xxx96com| 男女做爰动态图高潮gif福利片 | 欧美精品啪啪一区二区三区| 淫秽高清视频在线观看| 99国产精品一区二区三区| 天堂动漫精品| 好男人在线观看高清免费视频 | 亚洲人成伊人成综合网2020| 精品高清国产在线一区| 日本免费一区二区三区高清不卡 | 久久久国产成人免费| 午夜成年电影在线免费观看| 在线观看免费视频日本深夜| av在线播放免费不卡| 久久精品成人免费网站| 9色porny在线观看| 纯流量卡能插随身wifi吗| 国产成人精品久久二区二区91| 午夜久久久久精精品| 此物有八面人人有两片| 亚洲av成人av| 国内精品久久久久久久电影| 在线观看66精品国产| 亚洲人成77777在线视频| 99精品在免费线老司机午夜| 免费少妇av软件| 国产97色在线日韩免费| 夜夜看夜夜爽夜夜摸| 欧美日韩中文字幕国产精品一区二区三区 | 大陆偷拍与自拍| 男男h啪啪无遮挡| 啦啦啦免费观看视频1| 高清毛片免费观看视频网站| 美国免费a级毛片| 亚洲性夜色夜夜综合| 九色国产91popny在线| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成a人片在线一区二区| 一区二区三区激情视频| 国内毛片毛片毛片毛片毛片| 中文字幕高清在线视频| 黄片小视频在线播放| 久久久久久亚洲精品国产蜜桃av| 国产精品亚洲一级av第二区| 亚洲精品国产精品久久久不卡| 日韩大尺度精品在线看网址 | 国产伦人伦偷精品视频| 9色porny在线观看| 国产亚洲精品第一综合不卡| 午夜免费激情av| 亚洲人成网站在线播放欧美日韩| 欧美日韩亚洲综合一区二区三区_| 操美女的视频在线观看| 久久香蕉精品热| 午夜福利一区二区在线看| 法律面前人人平等表现在哪些方面| www.自偷自拍.com| 在线观看午夜福利视频| 亚洲国产欧美网| or卡值多少钱| 非洲黑人性xxxx精品又粗又长| 一级黄色大片毛片| 999久久久国产精品视频| 满18在线观看网站| 亚洲欧美日韩高清在线视频| 性欧美人与动物交配| 久久欧美精品欧美久久欧美| 精品熟女少妇八av免费久了| 女人被躁到高潮嗷嗷叫费观| 母亲3免费完整高清在线观看| 欧美不卡视频在线免费观看 | 50天的宝宝边吃奶边哭怎么回事| 免费少妇av软件| 成人三级做爰电影| 亚洲av成人不卡在线观看播放网| 美女 人体艺术 gogo| 日韩中文字幕欧美一区二区| 桃红色精品国产亚洲av| 亚洲欧美日韩高清在线视频| 国产精品一区二区三区四区久久 | 精品人妻在线不人妻| 非洲黑人性xxxx精品又粗又长| 久久欧美精品欧美久久欧美| 久久性视频一级片| 欧美绝顶高潮抽搐喷水| 99国产综合亚洲精品| 最新美女视频免费是黄的| 午夜精品国产一区二区电影| 高清黄色对白视频在线免费看| 亚洲一区中文字幕在线| 女性生殖器流出的白浆| www.熟女人妻精品国产| 午夜a级毛片| 欧美中文日本在线观看视频| 色老头精品视频在线观看| 女生性感内裤真人,穿戴方法视频| 亚洲国产中文字幕在线视频| 无遮挡黄片免费观看| 真人一进一出gif抽搐免费| 免费一级毛片在线播放高清视频 | 欧美中文综合在线视频| 色哟哟哟哟哟哟| 最近最新免费中文字幕在线| 又紧又爽又黄一区二区| 国产成人精品久久二区二区91| 成人精品一区二区免费| 日韩视频一区二区在线观看| 国产成人av激情在线播放| 日韩欧美一区二区三区在线观看| 国产av精品麻豆| 精品一区二区三区av网在线观看| 在线观看一区二区三区| 亚洲成av人片免费观看| 亚洲五月色婷婷综合| 中文字幕另类日韩欧美亚洲嫩草| 人妻丰满熟妇av一区二区三区| 人人妻人人爽人人添夜夜欢视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲中文字幕一区二区三区有码在线看 | 精品国产亚洲在线| 精品欧美国产一区二区三| 一区福利在线观看| 级片在线观看| 国内毛片毛片毛片毛片毛片| 午夜激情av网站| 国产97色在线日韩免费| 日韩 欧美 亚洲 中文字幕| 中文字幕av电影在线播放| 国产精品av久久久久免费| 亚洲国产欧美日韩在线播放| 色av中文字幕| 亚洲,欧美精品.| 女性被躁到高潮视频| 免费搜索国产男女视频| 亚洲欧美激情在线| 国产熟女午夜一区二区三区| 美女免费视频网站|