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

    規(guī)整填料塔中離子液體吸收CO2的傳質(zhì)與流體力學(xué)性能

    2015-11-26 03:20:09代成娜項銀雷志剛
    化工學(xué)報 2015年8期
    關(guān)鍵詞:規(guī)整傳質(zhì)陽離子

    代成娜,項銀,雷志剛

    ?

    規(guī)整填料塔中離子液體吸收CO2的傳質(zhì)與流體力學(xué)性能

    代成娜,項銀,雷志剛

    (北京化工大學(xué)化工資源有效利用國家重點實驗室,北京100029)

    在規(guī)整填料塔中采用離子液體吸收二氧化碳氣體,利用計算流體力學(xué)(CFD)軟件建立可靠的數(shù)學(xué)模型,系統(tǒng)考察了離子液體結(jié)構(gòu)及規(guī)整填料幾何參數(shù)對吸收過程的傳質(zhì)特性和流體力學(xué)性能的影響規(guī)律。結(jié)果表明,床層壓降隨氣體流速增大而增大,液相傳質(zhì)系數(shù)隨液體流速的增大而增大。相同陰離子時,隨著陽離子碳鏈長度的增長,吸收過程壓降增大,同時液相傳質(zhì)系數(shù)減小。相同陽離子時,不同陰離子的離子液體壓降大的同時傳質(zhì)系數(shù)也大。但離子液體的結(jié)構(gòu)對壓降影響不明顯。離子液體篩選主要考慮傳質(zhì)系數(shù)和溶解度因素,但二者與離子液體結(jié)構(gòu)表現(xiàn)出相反的規(guī)律。兩種折線結(jié)構(gòu)的規(guī)整填料傳質(zhì)性能優(yōu)于傳統(tǒng)的X型和Y型結(jié)構(gòu)。

    離子液體;傳質(zhì);二氧化碳;吸收;規(guī)整填料;壓降;計算流體力學(xué)

    引 言

    溫室效應(yīng)引起的全球變暖是當(dāng)前最為顯著的環(huán)境問題之一,而在總的溫室效應(yīng)中二氧化碳(CO2)的作用約占一半。從美國能源部二氧化碳信息分析中心(CDIAC)收集的數(shù)據(jù)可見,至2010年CO2的全球年排放量達到335億噸,而我國的CO2排放量達到82.4億噸,占世界總排放量的25%,居于榜首,比美國日本兩國之和還要多。因此,我國面臨著嚴峻的CO2減排和捕集挑戰(zhàn)[1]。我國CO2的工業(yè)排放主要來源于兩類:一類是煤等化石燃料燃燒所產(chǎn)生的煙道氣;另一類是化學(xué)工業(yè)中產(chǎn)生的合成氣。合成氣化工在現(xiàn)代化學(xué)工業(yè)中占有十分重要的地位,如從合成氣出發(fā)合成氨、合成烯烴、費托合成與制氫等。根據(jù)后續(xù)工藝的不同,工業(yè)上對合成氣的分離與凈化要求亦各不相同,有時非常嚴格(須脫除至20 ml·m-3以下),例如對于合成氨過程,CO2既是氨合成催化劑的毒物,又會導(dǎo)致后續(xù)過程產(chǎn)生NH4HCO3或(NH4)2CO3結(jié)晶進而造成管道和設(shè)備的堵塞。

    傳統(tǒng)的從工業(yè)氣體中吸收分離CO2的方法主要包括物理吸收法、化學(xué)吸收法、吸附分離法、膜分離法等,其中吸收法應(yīng)用最為廣泛[2-4]。化學(xué)吸收法主要發(fā)生在常溫常壓條件下,其優(yōu)點是吸收能力強、選擇性高,特別適宜常壓煙道氣中CO2的脫除,但其缺點是溶劑回收能耗高,且一般化學(xué)分離劑對反應(yīng)器設(shè)備腐蝕性大。物理吸收法通常通過改變溫度和壓力對氣體進行吸收或解吸,該方法一般適用于中高壓氣體的吸收,其優(yōu)點是易于實現(xiàn)溶劑的回收和循環(huán)利用。甲醇因其價廉易得而廣泛應(yīng)用于工業(yè)生產(chǎn)中,特別地,對于分離要求十分苛刻的CO2捕集過程,目前工業(yè)上廣泛采用的是低溫甲醇洗工藝(亦稱Rectisol法)[5-6],但由于甲醇具有很高的揮發(fā)性,在吸收過程中會造成大量的甲醇揮發(fā)損失,進而導(dǎo)致后續(xù)一系列的水洗及甲醇稀水溶液精餾回收等復(fù)雜工藝。因而迫切需要研究新的分離劑或新的分離技術(shù)以強化CO2吸收過程。

    離子液體(ionic liquids,ILs)是一種由特定陰、陽離子組成的100℃以下呈液體形態(tài)的有機鹽,由于具有蒸氣壓低、熱穩(wěn)定性好、溶解能力優(yōu)良、結(jié)構(gòu)與性質(zhì)可調(diào)節(jié)、可循環(huán)使用等優(yōu)勢而被認為是一種環(huán)境友好的綠色溶劑,廣泛應(yīng)用于氣體的吸收和分離過程[7-12]。大量研究發(fā)現(xiàn),二氧化碳與離子液體之間有著強烈的相互作用,使其具有較高溶解度。近年來大量研究人員及資源經(jīng)費投入離子液體吸收CO2的研究,研究成果頗豐。但目前研究主要集中于離子液體吸收CO2的相平衡研究,對其動力學(xué)特性研究鮮見報道。

    相對散堆填料來說,規(guī)整填料形狀較為規(guī)范,且在塔內(nèi)是均勻整齊堆砌的。規(guī)整填料是在20世紀60年代后出現(xiàn),因其分離效率高,阻力小,壓降低,且通量和操作彈性大,在石油化工、精細化工、香料、化肥等領(lǐng)域有著不可比擬的地位[13-17]。因此本文主要采用計算流體力學(xué)軟件CFD,研究在規(guī)整填料塔中離子液體吸收CO2過程的傳質(zhì)及其流體力學(xué)特性,探索其構(gòu)效關(guān)系,為后續(xù)進一步實驗或工業(yè)化提供基礎(chǔ)數(shù)據(jù)和理論支撐。

    1 填料及計算方法介紹

    1.1 規(guī)整填料

    本文的研究對象是比表面積為500 m2·m-3的金屬板波紋規(guī)整填料,峰高約6.3 mm,波紋板傾斜角分別為常規(guī)的30°(稱為500X型填料)和45°(稱為500Y型填料)填料,以及兩種新型的過渡結(jié)構(gòu)(折線結(jié)構(gòu))(稱為45°-30°-45°和30°-45°-30°),相鄰兩塊波紋板的傾斜方向相反,規(guī)整填料結(jié)構(gòu)如圖1所示。此種填料比表面積大,傳質(zhì)效率較高,適用于分離要求較高的操作。

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

    1.2.1 模型假設(shè)與計算域 本文采用CFD軟件,利用Fluent單元模塊對金屬波紋板規(guī)整填料中離子液體吸收CO2過程的壓降和傳質(zhì)性能進行研究。整體式結(jié)構(gòu)化填料具有規(guī)律的結(jié)構(gòu),因此在本文中截取其中一個單元進行數(shù)學(xué)模型的建立,計算域如圖2所示,包括方向兩個完整的三角形通道的寬度(25.2 mm)、方向兩塊金屬波紋板(傾斜方向相反)、方向(一塊整體式結(jié)構(gòu)化催化劑填料的高度,100 mm)。

    在CFD模擬過程中,離子液體吸收劑是第一相,CO2氣體作為模擬中的第二相。為考察規(guī)整填料內(nèi)氣液兩相流體力學(xué)性能,進行如下假設(shè):

    (1)所有的計算均在穩(wěn)態(tài)下進行;

    (2)為了區(qū)分出計算域中液相和氣相進出口區(qū)域,在模型中假設(shè)液相的出口和入口為很多分布均勻的小孔,如圖2所示;

    (3)忽略熱量傳遞;

    (4)設(shè)定均勻的氣液進口的有效速度、溫度和濃度,液體進口CO2濃度為零,氣液進出口溫度適中為298.15 K;

    (5)氣液出口為壓力出口。

    1.2.2 控制方程 本文沒有考慮能量的變化,采用了歐拉-歐拉兩相模型,考慮了相間相互作用,將氣液兩相看作是相互穿插的連續(xù)流動的流體,并且對每一相均做了連續(xù)性方程和動量平衡方程的求解。氣液兩相的連續(xù)性方程、動量和質(zhì)量平衡方程的偏微分形式如下:

    連續(xù)性方程:

    動量平衡方程:

    式中,為黏度,mPa·s;LG為氣液兩相相互作用力,在此數(shù)學(xué)模型中通過Schiller-Naumann方程[18]求解。

    質(zhì)量平衡方程:

    式中,(kg·m-3)、u(m·s-1)和φ分別代表密度、平均速度和體積分數(shù),表示氣相或者液相;D(m·s-2)是第相(氣相或液相)中組分的擴散系數(shù);w表示第相(氣相或液相)中組分的質(zhì)量分數(shù)。在計算域中,氣相體積分數(shù)和液相體積分數(shù)加和為1。

    吸收過程中液相傳質(zhì)系數(shù)按照低含量氣體吸收計算得出。

    式中,OL(m) 和OL是液相總傳質(zhì)單元高度和傳質(zhì)單元數(shù);Ka(kmol·s-1·m-3)和是傳質(zhì)系數(shù)和液相中CO2的摩爾分數(shù);(kmol·s-1·m-2)為液相流量。在計算過程設(shè)定液相入口不含CO2,即in=0;出口處濃度可以從模擬計算的結(jié)果中得出,而e為CO2的平衡濃度,可通過文獻獲得。

    1.2.3 邊界條件 規(guī)整填料模擬中邊界條件的設(shè)定如下:

    (1)均勻的氣液入口速度;

    (2)壓力出口作為液相和氣相的出口條件;

    (3)所有的外周壁視作對稱結(jié)構(gòu);

    (4)利用標準的壁面函數(shù)作為近壁區(qū)域和壁面的處理,均為無滑移壁面。

    1.2.4 數(shù)值方法 首先采用2.3.16版本的Gambit軟件構(gòu)建了金屬波紋板規(guī)整填料計算域,并進行了網(wǎng)格劃分,如圖3所示。氣液兩相的入口和出口面均使用非結(jié)構(gòu)化的三角網(wǎng)格加密,體網(wǎng)格采用非結(jié)構(gòu)化的Tet/Hybrid網(wǎng)格。網(wǎng)格劃分后將其導(dǎo)入6.3.26版本的Fluent軟件中進行求解,在動量方程的壓力-速度耦合中采用SIMPLE算法,離散化方法采用一階上風(fēng)方案。連續(xù)方程的殘差值在10-3,而其他方程的殘差值均在10-4以下。

    1.2.5 吸收劑及CO2氣體的物性參數(shù) CO2在離子液體中擴散系數(shù)(cm2·s-1)采用如下方法得到[19]

    式中,m(CO2)為CO2摩爾體積(34 cm3·mol-1);(mPa·s)和(g·cm-3)分別為溶劑的黏度和密度,從文獻中查閱得到[20-29]。

    常壓下298 K時,吸收劑和CO2的物性參數(shù)列于表1中。

    表1 離子液體和CO2的物性參數(shù)

    2 結(jié)果與討論

    2.1 網(wǎng)格無關(guān)性與模型可靠性驗證

    在CFD模擬過程中,連續(xù)的計算域被離散化,因此為保證網(wǎng)格劃分數(shù)量及質(zhì)量對模擬結(jié)果的準確性,首先對網(wǎng)格數(shù)量無關(guān)性進行了驗證。選擇500X型的規(guī)整填料、[EMIM][Tf2N]作為溶劑、液速為0.1 m·s-1、氣速為2.0 m·s-1的條件下進行驗證,結(jié)果如圖4(a)所示。在較低的網(wǎng)格數(shù)量下氣體壓降隨著網(wǎng)格數(shù)的增加而增大,當(dāng)網(wǎng)格數(shù)量大于52.6萬后壓降趨于固定值,對模擬結(jié)果影響不明顯。本文的后續(xù)計算過程中所劃分網(wǎng)格數(shù)均大于60萬。

    對比了500X型填料空氣-水體系的實驗和模擬結(jié)果。如圖4(b)所示,模擬結(jié)果和實驗結(jié)果基本相符,說明所建立的數(shù)學(xué)模型是基本可靠的。其中,實驗數(shù)據(jù)來自于文獻值[30]。模擬結(jié)果與實驗值的偏差可能與模型假設(shè)中液相出口和入口的小孔大小、數(shù)量及其分布有關(guān)。為進一步提高模擬的準確性,在后續(xù)的研究中應(yīng)具體加以考慮。

    2.2 離子液體結(jié)構(gòu)對吸收過程壓降及傳質(zhì)性能的影響

    為考察離子液體結(jié)構(gòu)對吸收過程壓降及傳質(zhì)性能的影響,選擇了6種常見的咪唑型離子液體([EMIM][Tf2N]、[BMIM][Tf2N]、[HMIM][Tf2N]、[EMIM][BF4]、[EMIM][TfO]、[EMIM][TFA]),離子液體名稱列于符號說明中。

    2.2.1 陽離子對壓降和傳質(zhì)系數(shù)的影響 首先考察了298.15 K、常壓下,在500X型金屬波紋板規(guī)整填料中,陰離子為雙(三氟甲基磺酰基)酰亞胺([Tf2N])、陽離子為不同烷基鏈長度的咪唑離子([EMIM]、[BMIM]和[HMIM])時,陽離子結(jié)構(gòu)對壓降和傳質(zhì)系數(shù)的影響。

    圖5給出了在相同液速(L=0.05 m·s-1)情況下,壓降和傳質(zhì)系數(shù)隨氣速的變化規(guī)律,可見隨氣速的增大,所有離子液體壓降均增大,但液相傳質(zhì)系數(shù)基本保持不變,可見傳質(zhì)過程主要是受液膜控制。

    圖6給出了在相同氣速(G=1.8 m·s-1)情況下壓降和傳質(zhì)系數(shù)隨液速的變化規(guī)律,可見增大液速,傳質(zhì)系數(shù)會快速增大。對于液膜控制吸收過程,增大液體流速,能夠使液體湍動增強,減小液膜厚度,同時增大氣液相的接觸面積,從而提高液相傳質(zhì)系數(shù)。相同氣速和液速的情況下,壓降從大到小[HMIM][Tf2N]>[BMIM][Tf2N]>[EMIM][Tf2N],液相傳質(zhì)系數(shù)從大到小依次是[EMIM][Tf2N]> [BMIM][Tf2N]>[HMIM][Tf2N]??梢?,陰離子相同時增加離子液體陽離子烷基鏈長度,壓降增大,同時傳質(zhì)系數(shù)減小,但CO2氣體溶解度會增加。隨著陽離子碳鏈的增長,離子液體中的自由體積增大,故溶解度增加;但同時離子液體黏度增大,CO2氣體在離子液體中的擴散系數(shù)變小,且傳質(zhì)阻力增大,故傳質(zhì)系數(shù)變小。因此篩選離子液體陽離子時需權(quán)衡CO2的溶解度(熱力學(xué)角度)和傳質(zhì)及動力消耗(動力學(xué)角度)。

    2.2.2 陰離子對壓降和傳質(zhì)系數(shù)的影響 隨后考察了298.15 K、常壓下,在500X型金屬波紋板規(guī)整填料中,相同陽離子([EMIM]),陰離子不同([Tf2N]、[BF4]、[TfO]和[TFA])的4種離子液體吸收CO2過程的壓降和傳質(zhì)性能。

    圖7和圖8分別給出了在相同液速(L=0.05 m·s-1)和相同氣速(G=1.8 m·s-1)情況下,壓降和傳質(zhì)系數(shù)隨氣速的變化規(guī)律,可見隨氣速的增大,所有離子液體壓降均增大;增大液速,傳質(zhì)系數(shù)會快速增大。相同氣速和液速的情況下,壓降從大到小依次為[EMIM][BF4] > [EMIM][TFA] ≈ [EMIM][TfO] > [EMIM][Tf2N],液相傳質(zhì)系數(shù)從大到小依次是[EMIM][BF4] > [EMIM][TFA] ≈ [EMIM][TfO] > [EMIM][Tf2N]??梢?,陽離子相同時,不同陰離子的離子液體壓降大的同時傳質(zhì)系數(shù)也增大,即傳質(zhì)系數(shù)的增大是依靠更多的動力消耗為代價的,這與陽離子對動力學(xué)影響表現(xiàn)出相反的規(guī)律。但不論陰離子還是陽離子對壓降的影響不是很明顯,所以離子液體篩選過程中應(yīng)以傳質(zhì)系數(shù)和溶解度為關(guān)鍵因素。對比表1中的溶解度數(shù)據(jù)可以發(fā)現(xiàn),不管陰離子還是陽離子結(jié)構(gòu)的離子液體,傳質(zhì)系數(shù)大的氣體的溶解度小。

    2.3 填料結(jié)構(gòu)對傳遞性能的影響

    圖9和圖10給出了波紋板傾斜角度對離子液體吸收CO2過程中壓降和傳質(zhì)系數(shù)的影響,包括4種不同的傾斜角結(jié)構(gòu)(500X: 30°、500Y: 45°和兩種折線結(jié)構(gòu): 30°-45°-30°和45°-30°-45°),在相同氣速和液速下,4種結(jié)構(gòu)壓降大小順序為30°<30°-45°-30°<45°-30°-45°<45°,傳質(zhì)系數(shù)大小順序為30°<45°≈30°-45°-30°<45°-30°-45°。可見兩種折線結(jié)構(gòu)盡管沒有降低氣體的壓降,但傳質(zhì)性能相對于傳統(tǒng)的X型和Y型有所提高,主要是因為兩次傾斜角的變化促進了氣液表面的更新,從而提高了氣液間的流動傳質(zhì)效率。因此,權(quán)衡壓降和傳質(zhì)性能的折線結(jié)構(gòu)填料是一個不錯的選擇。

    3 結(jié) 論

    本文主要對金屬波紋板規(guī)整填料中離子液體吸收CO2過程的傳質(zhì)和流體力學(xué)性能進行了研究,采用CFD軟件建立了可靠的數(shù)學(xué)模型,系統(tǒng)考察了離子液體陽離子和陰離子結(jié)構(gòu)以及規(guī)整填料波紋板的傾斜角度對吸收過程壓降和傳質(zhì)性能的影響。結(jié)果表明:

    (1)氣體通過單位高度床層壓降隨氣體流速增大而增大,液相傳質(zhì)系數(shù)隨液體流速的增大而增大。

    (2)相同陰離子的離子液體,隨著陽離子碳鏈長度的增長,吸收過程壓降增大([HMIM][Tf2N] > [BMIM][Tf2N] > [EMIM][Tf2N]),同時液相傳質(zhì)系數(shù)減?。╗EMIM][Tf2N] > [BMIM][Tf2N] > [HMIM][Tf2N])。

    (3)相同陽離子的離子液體,不同陰離子的離子液體壓降大的同時傳質(zhì)系數(shù)也增大,即傳質(zhì)系數(shù)的增大是依靠更多的動力消耗為代價的。壓降和傳質(zhì)系數(shù)從大到小均為[EMIM][BF4]>[EMIM][TFA]≈ [EMIM][TfO] > [EMIM][Tf2N]。

    (4)離子液體的結(jié)構(gòu)對吸收過程的壓降影響不明顯,所以離子液體篩選過程中主要考慮傳質(zhì)系數(shù)和溶解度,但二者與離子液體結(jié)構(gòu)表現(xiàn)出相反的規(guī)律。

    (5)4種不同的波紋板傾斜角規(guī)整填料壓降大小順序為30°<30°-45°-30°<45°-30°-45°<45°,傳質(zhì)系數(shù)大小順序為30°<45°≈30°-45°-30°<45°- 30°-45°。可見兩種折線結(jié)構(gòu)盡管沒有降低氣體的壓降,但傳質(zhì)性能優(yōu)于傳統(tǒng)的X型和Y型規(guī)整填料。

    符 號 說 明

    D——擴散系數(shù),m·s-2 F——相互作用力,N HOL——液相總傳質(zhì)單元高度,m Kxa——液相傳質(zhì)系數(shù),kmol·s-1·m-3 L——液相流量,kmol·s-1·m-2 M——摩爾質(zhì)量,g·mol-1 NOL——傳質(zhì)單元數(shù) Δp——壓降,Pa·m-1 u——速度,m·s-1 Vm——摩爾體積,cm3·mol-1 w——質(zhì)量分數(shù) x——摩爾分數(shù) α——傾斜角,(°) m——黏度,mPa·s r——密度,kg·m-3 下角標 e——平衡狀態(tài) G——氣體 in——入口 L——液體 out——出口 離子液體名稱 [EMIM][Tf2N]——1-ethyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide [BMIM][Tf2N]——1-butyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide [HMIM][Tf2N]——1-hexyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide [EMIM][BF4]——1-ethyl-3-methylimidazolium tetrafluoroborate [EMIM][TfO]——1-ethyl-3-methylimidazolium trifluoromethanesulfonate [EMIM][TFA]——1-ethyl-3-methylimidazolium trifluoroacetate

    [1] Fei Weiyang(費維揚), Ai Ning(艾寧), Chen Jian(陳健). Capture and separation of greenhouse gases CO2— the challenge and opportunity for separation technology [J].(化工進展), 2005, 24(1): 1-4.

    [2] Li Sifang(黎四芳), Ren Zhengwei(任錚偉), Li Pansheng(李盤生), Lu Qionghua(路瓊花). Absorption of CO2by aqueous MDEA-MEA solution [J].() (化工學(xué)報), 1994, 45(6): 698-703.

    [3] Krishnamurthy S, Rao V R, Gunkuka S, Sharratt P, Haghpanah R, Rajendran A, Amanullah M, Karimi I A, Farooq S. CO2capture from dry flue gas by vacuum swing adsorption: a pilot plant study [J].., 2014, 60(5): 1830-1842.

    [4] Aaron D, Tsouris C. Separation of CO2from flue gas: a review [J]., 2005, 40(1): 321-348.

    [5] Mei Lihua(梅麗華). Analysis of present production of low-temperature methanol absorbing device in China [J].(煤炭技術(shù)), 2008, 27(4): 118-120.

    [6] Qin Li(秦莉), Shi Xianlin(史賢林), Li Yu’an(李玉安), Zhao Haiquan(趙海泉). Process simulation and parameters optimization of Rectisol [J].(安徽化工), 2010, 36(3): 51-53.

    [7] Lei Z G, Dai C N, Chen B H. Gas solubility in ionic liquids [J]., 2014, 114(2): 1289-1326.

    [8] Dai C N, Lei Z G, Xiao L, Wang W, Chen B H. Group contribution lattice fluid equation of state (GCLF EOS) for CO2-ionic liquid systems [J].2013, 59(11): 4399-4412.

    [9] Lei Z G, Dai C N, Wang W, Chen B H. UNIFAC model for ionic liquid-CO2systems [J].2014, 60(2): 716-729.

    [10] Anthony J L, Anderson J L, Maginn E J, Brennecke J F. Anion effects on gas solubility in ionic liquids [J]., 2005, 109(13): 6366-6374.

    [11] Li Tiancheng(李天成), Kong Qinghua(孔慶華), Chen Chen(陳郴), Piao Xianglan(樸香蘭). CO2capture with complex absorbent of ionic liquid and water [J].(化工學(xué)報), 2013, 64(2): 600-608.

    [12] Wu Yanhui(伍艷輝), Zhang Tianhao(張?zhí)旌?, Duan Yongchao(段永超), Yu Shikun(于世昆). Supported ionic liquid membranes for CO2/CH4separation [J].(化工學(xué)報), 2013, 64(S1): 111-117.

    [13] Ratheesh S, Kannan A. Holdup and pressure drop studies in structured packings with catalysts [J]., 2004, 104(1-3): 45-54.

    [14] Hosseini S H, Shojaee S, Ahmadi G, Zivdar M. Computational fluid dynamics studies of dry and wet pressure drops in structured packings [J]., 2012, 18(4): 1465-1473.

    [15] Oluji? ?, Seibert A F, Kaibel B, Jansen H, Rietfort T, Zich E. Performance characteristics of a new high capacity structured packing [J]., 2003, 42(1): 55-60.

    [16] Zhou Yuan(周媛), Li Qunsheng(李群生), Zhang Zeting(張澤廷). The hydrodynamic performance of novel corrugated web packing [J].(北京化工大學(xué)學(xué)報), 2005, 32(3): 13-15.

    [17] Luo S J, Fei W Y, Song X Y, Li H Z. Effect of channel opening angle on the performance of structured packings [J]., 2008, 144(2): 227-234.

    [18] Schiller L, Naumann L. A drag coefficient correlation [J]., 1935, 77(1): 318-320

    [19] Morgan D, Ferguson L, Scovazzo P. Diffusivities of gases in room-temperature ionic liquids: data and correlations obtained using a lag-time technique [J].., 2005, 44(13): 4815-4823.

    [20] Won Y S, Chung D K, Mills A F. Density, viscosity, surface tension, and carbon dioxide solubility and diffusivity of methanol, ethanol, aqueous propanol, and aqueous ethylene glycol at 25℃[J]., 1981, 26(2): 140-141.

    [21] Yokozek A, Shiflett M B, Junk C P, Grieco L M, Foo T. Physical and chemical absorptions of carbon dioxide in room-temperature ionic liquids [J]., 2008, 112(51): 16654-16663.

    [22] Valderrama J O, Robles P A. Critical properties, normal boiling temperatures, and acentric factors of fifty ionic liquids [J]., 2007, 46(4): 1338-1344.

    [23] Yu G R, Zhao D C, Wen L, Yang S D, Chen X C. Viscosity of ionic liquids: database, observation, and quantitative structure-property relationship analysis [J]., 2012, 58(9): 2885-2899.

    [24] Anthony J L, Anderson J L, Maginn E J, Brennecke J F. Anion effects on gas solubility in ionic liquids [J]., 2005, 109(13): 6366-6374.

    [25] Muldoon M J, Aki S N V K, Anderson J L, Dixon J K, Brennecke J F. Improving carbon dioxide solubility in ionic liquids [J]., 2007, 111(30): 9001-9009.

    [26] Lei Z G, Zhang J G, Li Q S, Chen B H. UNIFAC model for ionic liquids [J]., 2009, 48(5): 2697-2704.

    [27] Valderrama J O, Sanga W W, Lazzús J A. Critical properties, normal boiling temperatures, and acentric factors of another 200 ionic liquids [J]., 2008, 47(4): 1318-1330.

    [28] Arce A, Rodil E, Soto A. Volumetric and viscosity study for the mixtures of 2-ethoxy-2-methylpropane, ethanol, and 1-ethyl-3- methylimidazolium ethyl sulfate ionic liquid [J]., 2006, 51(4): 1453-1457.

    [29] Shiflett M B, Yokozeki A. Phase behavior of carbon dioxide in ionic liquids: [emim][acetate], [emim][trifluoroacetate], and [emim][acetate] + [emim][trifluoroacetate] mixtures [J]., 2009, 54(1): 108-114.

    [30] Dai C N, Lei Z G, Li Q S, Chen B H. Pressure drop and mass transfer study in structured catalytic packings [J].2012, 98(1): 78-87.

    Mass transfer and hydraulic performance of CO2absorptionby ionic liquids over structured packings

    DAI Chengna, XIANG Yin, LEI Zhigang

    State Key Laboratory of Chemical Resource EngineeringBeijing University of Chemical TechnologyBeijingChina

    In this work, ionic liquids (ILs) were proposed for the capture of CO2over the structured packings. The pressure drop and mass transfer coefficient were investigated using computational fluid dynamics (CFD). The relationship among the structure of ILs, performance of pressure drop and mass transfer coefficient, as well as the relationship among the corrugation angles of structured packings, performance of pressure drop and mass transfer coefficient were identified. The results showed that under the same operating conditions with same anion, the increase of carbon number in the alkyl chain on the cation of ILs leads to the increase of the pressure drop (., [HMIM][Tf2N] > [BMIM][Tf2N] > [EMIM][Tf2N]), and decrease of the mass transfer coefficient (., [EMIM][Tf2N] > [BMIM][Tf2N] > [HMIM][Tf2N]). Thus, the cation [EMIM] is the optimum choice with regard to the pressure drop and mass transfer performance. Moreover, for ILs with different anions the pressure drop and mass transfer coefficient show the same trend,., [EMIM][BF4] > [EMIM][TFA] ≈ [EMIM][TfO] > [EMIM][Tf2N]. There is an inverse relationship between the mass transfer coefficient and solubility of CO2for different ILs. Thus, there is a tradeoff when choosing IL for the capture of CO2. In addition, four types of structured packings with different corrugation angles (X type, Y type, and two transition types) were compared. From the viewpoint of mass transfer performance the two transition types structured packings are superior to the traditional X or Y types.

    ionic liquid; mass transfer; carbon dioxide; absorption; structured packing; pressure drop; computational fluid dynamics (CFD)

    2015-06-01.

    Prof. LEI Zhigang, leizhg@mail.buct.edu.cn

    10.11949/j.issn.0438-1157.20150753

    TQ 028.8

    A

    0438—1157(2015)08—2953—09

    雷志剛。

    代成娜(1984—),女,博士,講師。

    國家自然科學(xué)基金項目(21476009,21406007,U1462104)。

    2015-06-01收到初稿,2015-06-09收到修改稿。

    supported by the National Natural Science Foundation of China (21476009, 21406007, U1462104).

    猜你喜歡
    規(guī)整傳質(zhì)陽離子
    300kt/a硫酸系統(tǒng)規(guī)整填料使用情況簡介
    烷基胺插層蒙脫土的陽離子交換容量研究
    提高日用玻璃陶瓷規(guī)整度和表面光滑度的處理方法
    佛山陶瓷(2016年11期)2016-12-23 08:50:27
    電梯的建筑化藝術(shù)探索
    大觀(2016年9期)2016-11-16 10:31:30
    氨基酸鹽吸收二氧化碳過程的傳質(zhì)特性
    基于發(fā)音機制的貪婪自適應(yīng)語音時長規(guī)整算法
    計算機工程(2015年8期)2015-07-03 12:20:34
    PTFE膜吸收CO2的工藝參數(shù)對傳質(zhì)性能的影響
    陽離子Gemini表面活性劑的應(yīng)用研究進展
    清潔轉(zhuǎn)向酸H+表面?zhèn)髻|(zhì)行為實驗研究
    煤顆粒熱解的傳熱傳質(zhì)分析
    計算物理(2014年1期)2014-03-11 17:00:36
    亚洲精品久久久久久婷婷小说| 久久精品国产自在天天线| 五月伊人婷婷丁香| 深夜a级毛片| 大片电影免费在线观看免费| 一级毛片我不卡| 久久久久久久大尺度免费视频| av专区在线播放| 性高湖久久久久久久久免费观看| 91aial.com中文字幕在线观看| 国产男女内射视频| 亚洲国产色片| a级一级毛片免费在线观看| 一级毛片黄色毛片免费观看视频| 十分钟在线观看高清视频www | 熟妇人妻不卡中文字幕| 亚洲欧美日韩东京热| 亚洲精品亚洲一区二区| 亚洲精品国产色婷婷电影| 大香蕉97超碰在线| 精品亚洲成国产av| 国产有黄有色有爽视频| 一本色道久久久久久精品综合| 国产成人a∨麻豆精品| 亚洲av综合色区一区| 人人妻人人看人人澡| 我的老师免费观看完整版| 桃花免费在线播放| 午夜91福利影院| 在线亚洲精品国产二区图片欧美 | 国产极品天堂在线| 一级a做视频免费观看| 一本大道久久a久久精品| 狠狠精品人妻久久久久久综合| 国产毛片在线视频| 精华霜和精华液先用哪个| 亚洲av电影在线观看一区二区三区| 亚洲国产精品成人久久小说| 一级毛片aaaaaa免费看小| 免费观看在线日韩| 欧美精品国产亚洲| 国产精品久久久久久精品电影小说| 国产色爽女视频免费观看| 夫妻性生交免费视频一级片| 女人精品久久久久毛片| 99久久精品国产国产毛片| 久久精品国产鲁丝片午夜精品| 女人久久www免费人成看片| 一区在线观看完整版| 久久午夜综合久久蜜桃| 免费看日本二区| 亚洲性久久影院| av在线观看视频网站免费| 男女国产视频网站| 丰满饥渴人妻一区二区三| 国产精品国产av在线观看| 成人毛片60女人毛片免费| 国产一区亚洲一区在线观看| 在线 av 中文字幕| 中文字幕久久专区| 亚洲一区二区三区欧美精品| 亚洲精品国产成人久久av| 国产 精品1| 久久精品熟女亚洲av麻豆精品| 亚洲欧美一区二区三区黑人 | av视频免费观看在线观看| 午夜福利网站1000一区二区三区| 欧美少妇被猛烈插入视频| 亚洲欧美中文字幕日韩二区| 91精品一卡2卡3卡4卡| 亚洲精品视频女| 在线观看三级黄色| 91在线精品国自产拍蜜月| 在线天堂最新版资源| 欧美激情国产日韩精品一区| 亚洲欧美成人精品一区二区| 日韩精品免费视频一区二区三区 | 又黄又爽又刺激的免费视频.| 成人美女网站在线观看视频| 日韩一区二区视频免费看| 人人妻人人看人人澡| 尾随美女入室| 亚洲经典国产精华液单| 大陆偷拍与自拍| 肉色欧美久久久久久久蜜桃| 婷婷色麻豆天堂久久| 深夜a级毛片| 国产精品一二三区在线看| 国内精品宾馆在线| 午夜激情福利司机影院| 99热这里只有精品一区| 熟妇人妻不卡中文字幕| 亚洲精品一区蜜桃| 国产色婷婷99| 国产成人一区二区在线| 亚洲经典国产精华液单| 亚洲av日韩在线播放| 亚洲真实伦在线观看| 极品人妻少妇av视频| 亚洲国产欧美在线一区| 亚洲国产精品一区三区| av国产久精品久网站免费入址| 狂野欧美白嫩少妇大欣赏| 一级毛片aaaaaa免费看小| 久久久亚洲精品成人影院| av一本久久久久| 麻豆精品久久久久久蜜桃| 午夜福利,免费看| 国产午夜精品久久久久久一区二区三区| 18禁在线播放成人免费| 国产爽快片一区二区三区| 亚洲经典国产精华液单| 色哟哟·www| 少妇人妻精品综合一区二区| 久久av网站| 日本色播在线视频| 免费观看的影片在线观看| 久久国产精品男人的天堂亚洲 | tube8黄色片| 高清午夜精品一区二区三区| 亚洲精品国产色婷婷电影| 欧美丝袜亚洲另类| 99精国产麻豆久久婷婷| a级毛片免费高清观看在线播放| 久久人人爽人人片av| 午夜福利,免费看| 最黄视频免费看| 亚洲婷婷狠狠爱综合网| 99热这里只有是精品50| 在线观看免费日韩欧美大片 | 欧美成人精品欧美一级黄| 精品酒店卫生间| 日本vs欧美在线观看视频 | 在线观看免费日韩欧美大片 | 蜜桃久久精品国产亚洲av| 人人妻人人澡人人爽人人夜夜| 最近最新中文字幕免费大全7| 中文乱码字字幕精品一区二区三区| 国产欧美日韩精品一区二区| 日本色播在线视频| 大香蕉97超碰在线| 国产精品人妻久久久久久| 欧美精品一区二区免费开放| 亚洲精品国产av成人精品| 成年女人在线观看亚洲视频| 最黄视频免费看| 午夜老司机福利剧场| 中文字幕久久专区| 国产一区二区在线观看日韩| 女性被躁到高潮视频| 久久毛片免费看一区二区三区| 国产精品偷伦视频观看了| 啦啦啦视频在线资源免费观看| 成年女人在线观看亚洲视频| 欧美精品一区二区大全| 欧美日韩视频高清一区二区三区二| 18禁动态无遮挡网站| 人人妻人人澡人人爽人人夜夜| 青青草视频在线视频观看| 成年av动漫网址| 肉色欧美久久久久久久蜜桃| 校园人妻丝袜中文字幕| 青春草视频在线免费观看| 免费高清在线观看视频在线观看| 午夜老司机福利剧场| 秋霞伦理黄片| 精品国产一区二区久久| 99热全是精品| 国产av码专区亚洲av| 2018国产大陆天天弄谢| 亚洲国产成人一精品久久久| 伊人亚洲综合成人网| 少妇裸体淫交视频免费看高清| 亚洲av在线观看美女高潮| 在线看a的网站| 成人18禁高潮啪啪吃奶动态图 | 国产乱人偷精品视频| 亚洲国产精品一区二区三区在线| 久久6这里有精品| 丝袜喷水一区| 深夜a级毛片| 国语对白做爰xxxⅹ性视频网站| 日韩欧美一区视频在线观看 | 日日啪夜夜爽| 国产亚洲5aaaaa淫片| 91久久精品国产一区二区成人| 在线观看一区二区三区激情| 3wmmmm亚洲av在线观看| 亚洲人与动物交配视频| 精品国产一区二区三区久久久樱花| 国产精品三级大全| 乱人伦中国视频| av网站免费在线观看视频| 亚洲精品久久午夜乱码| 国产欧美日韩精品一区二区| 色婷婷av一区二区三区视频| 狂野欧美白嫩少妇大欣赏| 国产在线免费精品| 搡老乐熟女国产| 黑人巨大精品欧美一区二区蜜桃 | av天堂中文字幕网| 亚洲av电影在线观看一区二区三区| 高清毛片免费看| 边亲边吃奶的免费视频| 午夜91福利影院| 午夜免费男女啪啪视频观看| 丰满人妻一区二区三区视频av| 一级黄片播放器| h日本视频在线播放| 有码 亚洲区| 嘟嘟电影网在线观看| 夜夜看夜夜爽夜夜摸| av一本久久久久| av线在线观看网站| 如日韩欧美国产精品一区二区三区 | 这个男人来自地球电影免费观看 | a级毛片免费高清观看在线播放| 欧美国产精品一级二级三级 | 岛国毛片在线播放| 国产午夜精品久久久久久一区二区三区| 国产成人免费观看mmmm| 亚洲情色 制服丝袜| 好男人视频免费观看在线| 亚洲欧美一区二区三区国产| 高清毛片免费看| 在线天堂最新版资源| 亚洲精品乱久久久久久| 婷婷色综合www| 又爽又黄a免费视频| av有码第一页| 精品国产乱码久久久久久小说| 国产精品女同一区二区软件| 欧美精品亚洲一区二区| 美女中出高潮动态图| 全区人妻精品视频| 七月丁香在线播放| 国产一级毛片在线| 亚洲精品视频女| 日韩欧美精品免费久久| 国产在线视频一区二区| av天堂久久9| 在线观看人妻少妇| 日韩伦理黄色片| 国产美女午夜福利| 日本-黄色视频高清免费观看| 在线观看美女被高潮喷水网站| 一级爰片在线观看| 日韩制服骚丝袜av| 亚洲国产成人一精品久久久| 国产视频首页在线观看| 蜜桃久久精品国产亚洲av| 男女边摸边吃奶| 18禁在线无遮挡免费观看视频| 日韩在线高清观看一区二区三区| 2018国产大陆天天弄谢| 久久久a久久爽久久v久久| 亚洲欧美精品专区久久| av不卡在线播放| 麻豆成人av视频| 国产精品人妻久久久影院| 亚洲国产欧美日韩在线播放 | 国产69精品久久久久777片| 成年美女黄网站色视频大全免费 | 国产精品蜜桃在线观看| 亚洲精品,欧美精品| 2021少妇久久久久久久久久久| 久久久久久久久久成人| 亚洲精品国产成人久久av| 极品少妇高潮喷水抽搐| 国产欧美日韩综合在线一区二区 | 永久网站在线| 欧美日韩国产mv在线观看视频| 国产伦理片在线播放av一区| a级毛片免费高清观看在线播放| 国产日韩欧美亚洲二区| 一区二区av电影网| 中文字幕制服av| freevideosex欧美| 乱系列少妇在线播放| 精品国产一区二区三区久久久樱花| 青青草视频在线视频观看| 精品视频人人做人人爽| 日韩三级伦理在线观看| 精品亚洲成a人片在线观看| 免费观看a级毛片全部| 国产美女午夜福利| 午夜福利网站1000一区二区三区| 国产精品久久久久成人av| 黄色一级大片看看| 人人澡人人妻人| 精品视频人人做人人爽| 中文精品一卡2卡3卡4更新| 国产欧美日韩一区二区三区在线 | 国国产精品蜜臀av免费| 亚洲三级黄色毛片| 中国三级夫妇交换| 日日啪夜夜爽| 国产一区二区在线观看av| 国产精品久久久久久精品古装| av网站免费在线观看视频| 亚洲va在线va天堂va国产| www.av在线官网国产| 亚洲av电影在线观看一区二区三区| 欧美97在线视频| 久久久久久久久大av| 特大巨黑吊av在线直播| 2022亚洲国产成人精品| av女优亚洲男人天堂| 插逼视频在线观看| 91久久精品国产一区二区成人| 亚洲国产精品999| 亚洲内射少妇av| videos熟女内射| 欧美bdsm另类| 亚洲精品国产av蜜桃| 最新的欧美精品一区二区| 五月天丁香电影| 99视频精品全部免费 在线| 免费观看无遮挡的男女| 岛国毛片在线播放| 亚洲色图综合在线观看| 人妻制服诱惑在线中文字幕| 亚洲情色 制服丝袜| 亚洲人成网站在线观看播放| 黑人猛操日本美女一级片| 亚洲第一av免费看| 日日撸夜夜添| 日本91视频免费播放| 少妇丰满av| 男人狂女人下面高潮的视频| 久久国产乱子免费精品| 免费看日本二区| 国产亚洲5aaaaa淫片| 亚洲人与动物交配视频| 男女边摸边吃奶| 久久久久网色| 欧美精品一区二区大全| 国产精品99久久99久久久不卡 | 国产国拍精品亚洲av在线观看| 乱人伦中国视频| 精品国产乱码久久久久久小说| 欧美3d第一页| 亚洲精品日韩av片在线观看| 极品教师在线视频| 3wmmmm亚洲av在线观看| 精品99又大又爽又粗少妇毛片| 如日韩欧美国产精品一区二区三区 | 日本91视频免费播放| 国产精品久久久久久久电影| 青春草亚洲视频在线观看| 国产一区有黄有色的免费视频| 最黄视频免费看| 成人二区视频| 日本色播在线视频| 最近的中文字幕免费完整| 久久久久国产网址| 国产黄片视频在线免费观看| 亚洲人成网站在线观看播放| 亚洲av成人精品一二三区| 99热全是精品| 久久99蜜桃精品久久| 韩国av在线不卡| 日韩av不卡免费在线播放| av女优亚洲男人天堂| 亚洲国产精品一区二区三区在线| 国产片特级美女逼逼视频| 亚洲欧美精品自产自拍| 色婷婷久久久亚洲欧美| 人妻夜夜爽99麻豆av| 一本大道久久a久久精品| 精品国产一区二区三区久久久樱花| 六月丁香七月| 黄色日韩在线| 人妻人人澡人人爽人人| 欧美三级亚洲精品| 麻豆成人午夜福利视频| 久久免费观看电影| 一级黄片播放器| 2022亚洲国产成人精品| 国产无遮挡羞羞视频在线观看| av在线观看视频网站免费| 日韩一区二区三区影片| 如日韩欧美国产精品一区二区三区 | 午夜老司机福利剧场| 国产高清不卡午夜福利| 国产亚洲精品久久久com| 免费观看在线日韩| 麻豆成人午夜福利视频| 日韩av在线免费看完整版不卡| 简卡轻食公司| 午夜日本视频在线| 热re99久久精品国产66热6| 国产午夜精品久久久久久一区二区三区| 性色av一级| 国产一区二区三区av在线| 免费看av在线观看网站| 少妇被粗大猛烈的视频| av线在线观看网站| 免费不卡的大黄色大毛片视频在线观看| 一个人免费看片子| 国产亚洲91精品色在线| 赤兔流量卡办理| 搡老乐熟女国产| 亚洲国产精品999| 少妇的逼好多水| 十八禁高潮呻吟视频 | 久久午夜福利片| 亚洲欧美成人精品一区二区| 亚州av有码| 亚洲成人一二三区av| 欧美日韩综合久久久久久| 亚洲情色 制服丝袜| 精品99又大又爽又粗少妇毛片| 在线看a的网站| 色吧在线观看| 久久精品久久久久久久性| 午夜福利在线观看免费完整高清在| 久久韩国三级中文字幕| 日韩欧美一区视频在线观看 | 香蕉精品网在线| 国产永久视频网站| 人妻人人澡人人爽人人| 国产在线视频一区二区| 欧美另类一区| av在线老鸭窝| 69精品国产乱码久久久| 国产黄色视频一区二区在线观看| 最近中文字幕2019免费版| 18禁在线播放成人免费| 免费少妇av软件| av福利片在线| 日韩一区二区三区影片| 美女主播在线视频| 新久久久久国产一级毛片| 爱豆传媒免费全集在线观看| 国产伦精品一区二区三区四那| 亚洲精品国产成人久久av| 国产毛片在线视频| 午夜福利视频精品| 成年美女黄网站色视频大全免费 | 男女无遮挡免费网站观看| 国产亚洲午夜精品一区二区久久| a级片在线免费高清观看视频| 在线观看美女被高潮喷水网站| 老熟女久久久| av免费在线看不卡| 国产伦精品一区二区三区四那| 久久久久人妻精品一区果冻| 看免费成人av毛片| 最近手机中文字幕大全| 美女大奶头黄色视频| 三级国产精品片| 美女主播在线视频| 丝袜喷水一区| av一本久久久久| 久久精品国产亚洲av涩爱| 男男h啪啪无遮挡| 久久国产精品大桥未久av | 2022亚洲国产成人精品| 中文字幕制服av| 我的女老师完整版在线观看| 久久这里有精品视频免费| 搡老乐熟女国产| 久久鲁丝午夜福利片| 99九九在线精品视频 | 欧美3d第一页| 极品人妻少妇av视频| 国产免费又黄又爽又色| 久久午夜综合久久蜜桃| 久久久久久久久久久免费av| 日本欧美国产在线视频| 男人舔奶头视频| 亚洲国产欧美日韩在线播放 | 丝瓜视频免费看黄片| 26uuu在线亚洲综合色| 香蕉精品网在线| 国产一区二区在线观看日韩| 男人狂女人下面高潮的视频| 中文欧美无线码| 亚洲精品中文字幕在线视频 | 人人澡人人妻人| 欧美人与善性xxx| 观看美女的网站| 国产毛片在线视频| 综合色丁香网| 少妇的逼水好多| 多毛熟女@视频| 视频中文字幕在线观看| 精品少妇黑人巨大在线播放| 色婷婷av一区二区三区视频| 自线自在国产av| 少妇人妻 视频| 中文字幕亚洲精品专区| av福利片在线| 99热这里只有精品一区| 少妇丰满av| 亚洲人成网站在线播| 欧美bdsm另类| 中文字幕免费在线视频6| 欧美bdsm另类| 黑人巨大精品欧美一区二区蜜桃 | 老女人水多毛片| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久久大尺度免费视频| 成人黄色视频免费在线看| 丰满饥渴人妻一区二区三| 久久久亚洲精品成人影院| 熟女人妻精品中文字幕| 男的添女的下面高潮视频| 精品一区二区三区视频在线| 午夜福利在线观看免费完整高清在| 精品久久久精品久久久| 欧美日韩综合久久久久久| 99九九在线精品视频 | 大片免费播放器 马上看| 国产成人精品福利久久| 久久精品久久精品一区二区三区| 精品人妻熟女毛片av久久网站| 成人漫画全彩无遮挡| 尾随美女入室| 在线观看免费高清a一片| 在线天堂最新版资源| 久久精品国产亚洲网站| 一区二区av电影网| av视频免费观看在线观看| 卡戴珊不雅视频在线播放| 国产成人精品久久久久久| 啦啦啦啦在线视频资源| 中国三级夫妇交换| 高清视频免费观看一区二区| 最近手机中文字幕大全| 国产精品福利在线免费观看| 久久久国产欧美日韩av| 久久人人爽av亚洲精品天堂| 国产熟女欧美一区二区| 国产乱来视频区| 亚洲国产精品专区欧美| 亚洲欧美成人精品一区二区| 日韩欧美精品免费久久| 一区二区三区精品91| 乱码一卡2卡4卡精品| 成人漫画全彩无遮挡| 精品国产露脸久久av麻豆| 国产精品国产三级国产专区5o| 亚洲av男天堂| 国产毛片在线视频| 国内精品宾馆在线| 日本午夜av视频| 精品视频人人做人人爽| 91精品国产国语对白视频| 亚洲国产精品一区二区三区在线| 特大巨黑吊av在线直播| 精品久久久久久久久av| 少妇高潮的动态图| 国产黄色视频一区二区在线观看| 国产精品秋霞免费鲁丝片| 久久这里有精品视频免费| 国产成人a∨麻豆精品| 亚洲激情五月婷婷啪啪| 日日啪夜夜爽| 偷拍熟女少妇极品色| 免费观看在线日韩| 久久精品久久精品一区二区三区| 亚洲国产精品一区二区三区在线| 18禁动态无遮挡网站| 午夜久久久在线观看| h日本视频在线播放| 亚洲精品成人av观看孕妇| 久久久久久久久大av| 色5月婷婷丁香| 日本欧美国产在线视频| 婷婷色综合大香蕉| 午夜免费观看性视频| 精品卡一卡二卡四卡免费| 一级毛片黄色毛片免费观看视频| 午夜激情福利司机影院| 99久久中文字幕三级久久日本| 人人妻人人爽人人添夜夜欢视频 | 十八禁网站网址无遮挡 | 亚洲天堂av无毛| 亚洲三级黄色毛片| 国产成人a∨麻豆精品| 婷婷色麻豆天堂久久| 一个人免费看片子| av国产精品久久久久影院| 22中文网久久字幕| 日日啪夜夜爽| 综合色丁香网| 边亲边吃奶的免费视频| 国产深夜福利视频在线观看| 亚洲国产精品一区三区| 在线看a的网站| 内地一区二区视频在线| 在线观看国产h片| 制服丝袜香蕉在线| 99热这里只有精品一区| 搡老乐熟女国产| 国产成人精品久久久久久| 亚洲人成网站在线播| 噜噜噜噜噜久久久久久91| 国产精品久久久久久久电影| 晚上一个人看的免费电影| 丰满乱子伦码专区| 欧美日韩一区二区视频在线观看视频在线| 免费看av在线观看网站| 国产免费又黄又爽又色| 久久国产亚洲av麻豆专区| 亚洲av欧美aⅴ国产| av黄色大香蕉| 亚洲精品久久午夜乱码| 久久免费观看电影| 丁香六月天网| 日日爽夜夜爽网站| 免费高清在线观看视频在线观看| 久久国产精品男人的天堂亚洲 |