• <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
    一级二级三级毛片免费看| 女性被躁到高潮视频| 天天躁夜夜躁狠狠久久av| 在线播放无遮挡| 七月丁香在线播放| 简卡轻食公司| 自拍偷自拍亚洲精品老妇| 亚洲国产av新网站| 人人妻人人爽人人添夜夜欢视频 | 深夜a级毛片| 国产久久久一区二区三区| 国产精品一区www在线观看| 欧美日韩在线观看h| 日本av免费视频播放| 日本wwww免费看| 看十八女毛片水多多多| 亚洲国产精品999| 精品亚洲成a人片在线观看 | 一本色道久久久久久精品综合| 免费播放大片免费观看视频在线观看| 久久精品熟女亚洲av麻豆精品| 免费观看的影片在线观看| 伦精品一区二区三区| 十八禁网站网址无遮挡 | 成年免费大片在线观看| 国产片特级美女逼逼视频| 午夜福利影视在线免费观看| 综合色丁香网| 青青草视频在线视频观看| 国产精品女同一区二区软件| xxx大片免费视频| xxx大片免费视频| 成人漫画全彩无遮挡| 国产69精品久久久久777片| 色视频www国产| 免费黄网站久久成人精品| 老师上课跳d突然被开到最大视频| 亚洲美女搞黄在线观看| 男人添女人高潮全过程视频| 另类亚洲欧美激情| 国产精品99久久久久久久久| 日韩成人伦理影院| 久久久久国产网址| 熟女电影av网| 久久精品人妻少妇| 99热国产这里只有精品6| 中文精品一卡2卡3卡4更新| 中文字幕av成人在线电影| 五月伊人婷婷丁香| 各种免费的搞黄视频| 又粗又硬又长又爽又黄的视频| 男女边吃奶边做爰视频| 国产男人的电影天堂91| 欧美亚洲 丝袜 人妻 在线| 久久99热6这里只有精品| 中文字幕免费在线视频6| 99久久精品一区二区三区| 免费看光身美女| 亚洲中文av在线| 制服丝袜香蕉在线| 亚洲精品视频女| 一区在线观看完整版| 高清日韩中文字幕在线| 国产亚洲最大av| 国产高清有码在线观看视频| 一区二区三区精品91| 亚洲国产毛片av蜜桃av| 欧美高清成人免费视频www| 欧美日韩一区二区视频在线观看视频在线| 国产精品偷伦视频观看了| 高清日韩中文字幕在线| 国产精品人妻久久久影院| 99久久中文字幕三级久久日本| 亚洲av电影在线观看一区二区三区| 色婷婷av一区二区三区视频| 国产亚洲91精品色在线| 一边亲一边摸免费视频| 男男h啪啪无遮挡| 人妻系列 视频| 乱码一卡2卡4卡精品| 国产伦理片在线播放av一区| 亚洲成人av在线免费| 边亲边吃奶的免费视频| 只有这里有精品99| 国产精品一及| 久久久久久久精品精品| 老司机影院成人| 国产深夜福利视频在线观看| 80岁老熟妇乱子伦牲交| 18+在线观看网站| 国产高清不卡午夜福利| 久久99精品国语久久久| xxx大片免费视频| 亚洲色图综合在线观看| 下体分泌物呈黄色| 国产69精品久久久久777片| 内地一区二区视频在线| 久久国产精品男人的天堂亚洲 | 99热这里只有是精品50| 少妇的逼水好多| 人人妻人人看人人澡| 麻豆成人av视频| 国产精品一区二区三区四区免费观看| 久久精品人妻少妇| 高清日韩中文字幕在线| 黑人高潮一二区| 交换朋友夫妻互换小说| 交换朋友夫妻互换小说| 看非洲黑人一级黄片| 精品午夜福利在线看| 大片电影免费在线观看免费| 一本久久精品| 中文字幕免费在线视频6| 亚洲丝袜综合中文字幕| 午夜视频国产福利| 亚洲在久久综合| 亚洲电影在线观看av| 黄片wwwwww| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品乱码久久久v下载方式| 国产又色又爽无遮挡免| 黄色日韩在线| 夫妻性生交免费视频一级片| 国内揄拍国产精品人妻在线| 国产综合精华液| 精品人妻熟女av久视频| 中文字幕亚洲精品专区| 赤兔流量卡办理| 欧美区成人在线视频| 亚洲精品日韩av片在线观看| 97超碰精品成人国产| 最后的刺客免费高清国语| 国产精品久久久久久精品电影小说 | 少妇丰满av| 午夜激情久久久久久久| 日本爱情动作片www.在线观看| 在线免费观看不下载黄p国产| 国产一级毛片在线| 小蜜桃在线观看免费完整版高清| 国产免费福利视频在线观看| 日日撸夜夜添| 久久国产亚洲av麻豆专区| 日韩欧美精品免费久久| 伊人久久国产一区二区| 女性生殖器流出的白浆| 国产黄色免费在线视频| 激情 狠狠 欧美| 国产人妻一区二区三区在| 精品一区二区三卡| 在线天堂最新版资源| 亚洲怡红院男人天堂| 中国三级夫妇交换| 97热精品久久久久久| 中文字幕制服av| 午夜免费观看性视频| 亚洲av成人精品一区久久| 王馨瑶露胸无遮挡在线观看| 秋霞在线观看毛片| 人妻 亚洲 视频| 一级毛片黄色毛片免费观看视频| 国产一区二区在线观看日韩| 综合色丁香网| 国产精品一二三区在线看| 涩涩av久久男人的天堂| 岛国毛片在线播放| 亚洲av中文av极速乱| 中国三级夫妇交换| 亚洲色图综合在线观看| 这个男人来自地球电影免费观看 | 伦精品一区二区三区| 色视频www国产| 久久精品国产自在天天线| 久久久国产一区二区| 国产高清有码在线观看视频| 精品久久久精品久久久| 国产 一区 欧美 日韩| 亚洲国产高清在线一区二区三| 各种免费的搞黄视频| 美女cb高潮喷水在线观看| 亚洲丝袜综合中文字幕| 午夜福利在线在线| 观看美女的网站| 久久久久久久国产电影| 成人亚洲欧美一区二区av| 日本vs欧美在线观看视频 | 91久久精品国产一区二区成人| 在线亚洲精品国产二区图片欧美 | 97精品久久久久久久久久精品| 人妻 亚洲 视频| 街头女战士在线观看网站| 国产午夜精品久久久久久一区二区三区| 十分钟在线观看高清视频www | 精品国产一区二区三区久久久樱花 | 国产成人aa在线观看| 熟妇人妻不卡中文字幕| 国产又色又爽无遮挡免| 亚洲一区二区三区欧美精品| 国产精品一区二区在线不卡| 亚洲va在线va天堂va国产| 纯流量卡能插随身wifi吗| 毛片女人毛片| 亚洲欧美一区二区三区国产| 99久久精品一区二区三区| h日本视频在线播放| 精品国产乱码久久久久久小说| 国产精品一及| 亚洲电影在线观看av| 在线观看一区二区三区| 日本av手机在线免费观看| 热99国产精品久久久久久7| 免费黄频网站在线观看国产| 国产精品一区二区三区四区免费观看| 久久久亚洲精品成人影院| 久久久久久久久大av| 91精品国产国语对白视频| a级毛色黄片| kizo精华| 亚洲第一av免费看| 国产片特级美女逼逼视频| 乱码一卡2卡4卡精品| 全区人妻精品视频| 色婷婷av一区二区三区视频| 毛片一级片免费看久久久久| 老熟女久久久| 老司机影院毛片| 久久久国产一区二区| 欧美丝袜亚洲另类| 国产探花极品一区二区| 免费人妻精品一区二区三区视频| 国产精品一及| 亚洲欧美中文字幕日韩二区| 亚洲人成网站在线观看播放| 国产精品成人在线| av在线观看视频网站免费| 超碰av人人做人人爽久久| 蜜桃久久精品国产亚洲av| 亚洲成人av在线免费| 大码成人一级视频| av不卡在线播放| 美女中出高潮动态图| 久久久亚洲精品成人影院| 亚洲av电影在线观看一区二区三区| 成年免费大片在线观看| 国产在线视频一区二区| av国产免费在线观看| 啦啦啦中文免费视频观看日本| 97精品久久久久久久久久精品| 搡老乐熟女国产| 丰满迷人的少妇在线观看| 久久精品夜色国产| 女的被弄到高潮叫床怎么办| 乱码一卡2卡4卡精品| 日韩av不卡免费在线播放| 成人18禁高潮啪啪吃奶动态图 | 亚洲欧美精品自产自拍| 一级a做视频免费观看| 韩国高清视频一区二区三区| 久久精品国产亚洲av涩爱| 最黄视频免费看| 国内精品宾馆在线| 欧美老熟妇乱子伦牲交| av网站免费在线观看视频| 午夜免费鲁丝| 亚洲av在线观看美女高潮| 美女脱内裤让男人舔精品视频| 一本色道久久久久久精品综合| 女性被躁到高潮视频| 人妻一区二区av| 又大又黄又爽视频免费| 亚洲va在线va天堂va国产| 人妻 亚洲 视频| 街头女战士在线观看网站| 天美传媒精品一区二区| 欧美极品一区二区三区四区| 色综合色国产| www.色视频.com| 王馨瑶露胸无遮挡在线观看| 国产91av在线免费观看| 国产又色又爽无遮挡免| 身体一侧抽搐| 中国三级夫妇交换| 一本一本综合久久| 一级av片app| 日本黄色日本黄色录像| 中文字幕久久专区| 国产男女超爽视频在线观看| 日韩中字成人| 亚洲成人手机| 99久久人妻综合| 在线天堂最新版资源| 国产欧美另类精品又又久久亚洲欧美| 日韩免费高清中文字幕av| 蜜桃久久精品国产亚洲av| 成人毛片a级毛片在线播放| 亚洲国产欧美人成| 日本黄色片子视频| 人人妻人人澡人人爽人人夜夜| 日本vs欧美在线观看视频 | 久久精品久久精品一区二区三区| 少妇人妻 视频| 日韩成人av中文字幕在线观看| 51国产日韩欧美| 欧美日本视频| 免费不卡的大黄色大毛片视频在线观看| 最近的中文字幕免费完整| 欧美区成人在线视频| 有码 亚洲区| 韩国av在线不卡| 大香蕉久久网| 精品久久久久久久久av| 欧美成人午夜免费资源| 亚洲av欧美aⅴ国产| 国产精品一区二区在线不卡| 久久久久久久大尺度免费视频| 亚洲经典国产精华液单| 麻豆国产97在线/欧美| 欧美高清性xxxxhd video| 交换朋友夫妻互换小说| 国精品久久久久久国模美| 一个人免费看片子| 国产精品福利在线免费观看| 丰满少妇做爰视频| 一级毛片aaaaaa免费看小| 国产亚洲5aaaaa淫片| 久久97久久精品| 国产在线一区二区三区精| 99视频精品全部免费 在线| 一本色道久久久久久精品综合| 欧美xxxx性猛交bbbb| 国产白丝娇喘喷水9色精品| 夫妻午夜视频| 亚洲人成网站在线播| 男女边摸边吃奶| 亚洲va在线va天堂va国产| 久久av网站| 99久久精品热视频| 免费在线观看成人毛片| 99视频精品全部免费 在线| 国产综合精华液| 成人综合一区亚洲| 男女下面进入的视频免费午夜| 色婷婷久久久亚洲欧美| 看免费成人av毛片| 国国产精品蜜臀av免费| 精品99又大又爽又粗少妇毛片| 欧美精品一区二区大全| 少妇裸体淫交视频免费看高清| 一级爰片在线观看| 亚洲av国产av综合av卡| 伦精品一区二区三区| 夜夜看夜夜爽夜夜摸| 夫妻午夜视频| 亚洲美女搞黄在线观看| 欧美精品人与动牲交sv欧美| 91在线精品国自产拍蜜月| 新久久久久国产一级毛片| 亚洲精品日韩在线中文字幕| 欧美成人a在线观看| 国产一区有黄有色的免费视频| 国产爱豆传媒在线观看| 黄色欧美视频在线观看| 成人高潮视频无遮挡免费网站| 亚洲精品乱码久久久久久按摩| 亚洲,一卡二卡三卡| 联通29元200g的流量卡| 亚洲天堂av无毛| 中文字幕精品免费在线观看视频 | 男女啪啪激烈高潮av片| 99精国产麻豆久久婷婷| 一级毛片电影观看| 伊人久久国产一区二区| 亚洲精品aⅴ在线观看| 国产黄片视频在线免费观看| av专区在线播放| 久久精品久久久久久噜噜老黄| 毛片女人毛片| 亚洲最大成人中文| 国产精品成人在线| 久久99精品国语久久久| 校园人妻丝袜中文字幕| 国产色婷婷99| 国产伦精品一区二区三区视频9| 人妻系列 视频| 亚洲av二区三区四区| 国产高清国产精品国产三级 | 久久热精品热| 王馨瑶露胸无遮挡在线观看| 高清视频免费观看一区二区| 欧美国产精品一级二级三级 | 国产淫语在线视频| 亚洲人成网站在线观看播放| 亚洲国产毛片av蜜桃av| 午夜福利在线观看免费完整高清在| 偷拍熟女少妇极品色| 成人美女网站在线观看视频| 欧美97在线视频| 天堂中文最新版在线下载| 在线观看美女被高潮喷水网站| 伊人久久国产一区二区| 观看av在线不卡| 少妇 在线观看| 久久国产精品大桥未久av | 欧美bdsm另类| 99久国产av精品国产电影| 日韩,欧美,国产一区二区三区| 久久久久久久久大av| 天堂俺去俺来也www色官网| 国产成人精品久久久久久| 国精品久久久久久国模美| 97超碰精品成人国产| 有码 亚洲区| 欧美一级a爱片免费观看看| 少妇人妻久久综合中文| 婷婷色av中文字幕| 久久精品人妻少妇| 国产淫片久久久久久久久| 久久久久久久久久人人人人人人| 欧美性感艳星| 80岁老熟妇乱子伦牲交| 一级毛片 在线播放| 国产精品久久久久久精品古装| 国产精品99久久久久久久久| 一级爰片在线观看| 汤姆久久久久久久影院中文字幕| 国产精品久久久久成人av| 色5月婷婷丁香| 国产精品免费大片| 97在线人人人人妻| 欧美精品一区二区免费开放| 国产美女午夜福利| 欧美区成人在线视频| 精品熟女少妇av免费看| 日韩人妻高清精品专区| 久久国产亚洲av麻豆专区| 日韩在线高清观看一区二区三区| 国产男人的电影天堂91| 纯流量卡能插随身wifi吗| 久久久久久九九精品二区国产| 国产精品久久久久久av不卡| 欧美+日韩+精品| 青春草亚洲视频在线观看| h视频一区二区三区| 成人国产av品久久久| 卡戴珊不雅视频在线播放| 99久久精品一区二区三区| 色吧在线观看| 久久久久久久大尺度免费视频| 国产 一区精品| 日韩 亚洲 欧美在线| 三级国产精品欧美在线观看| 校园人妻丝袜中文字幕| 国产成人精品一,二区| 国产欧美另类精品又又久久亚洲欧美| 在线亚洲精品国产二区图片欧美 | 午夜福利影视在线免费观看| h日本视频在线播放| 国产亚洲最大av| 婷婷色av中文字幕| 男女免费视频国产| 少妇猛男粗大的猛烈进出视频| 国产伦理片在线播放av一区| 91久久精品国产一区二区成人| 亚洲精品,欧美精品| 欧美xxxx黑人xx丫x性爽| 亚洲最大成人中文| 嫩草影院新地址| 91精品一卡2卡3卡4卡| 亚洲内射少妇av| 国产在线一区二区三区精| 亚洲真实伦在线观看| 久久精品夜色国产| 永久免费av网站大全| 性色avwww在线观看| 亚洲最大成人中文| 在线观看国产h片| 丝袜喷水一区| 我的女老师完整版在线观看| 精品亚洲乱码少妇综合久久| 26uuu在线亚洲综合色| 国产一级毛片在线| 欧美日韩视频精品一区| 成人二区视频| 国产午夜精品久久久久久一区二区三区| 男人爽女人下面视频在线观看| 亚洲成人av在线免费| 国产免费一区二区三区四区乱码| 免费av不卡在线播放| 18禁裸乳无遮挡免费网站照片| 国产有黄有色有爽视频| 久久韩国三级中文字幕| 日本一二三区视频观看| 亚洲成人av在线免费| 欧美xxxx黑人xx丫x性爽| 亚洲成色77777| 国产免费视频播放在线视频| 免费大片黄手机在线观看| 黄色怎么调成土黄色| 久久久久久久亚洲中文字幕| 国产在线男女| 九色成人免费人妻av| 男人添女人高潮全过程视频| 日日撸夜夜添| 天堂俺去俺来也www色官网| 午夜免费观看性视频| a级毛片免费高清观看在线播放| 久热这里只有精品99| 麻豆国产97在线/欧美| 人妻 亚洲 视频| 我要看日韩黄色一级片| 中文资源天堂在线| 天美传媒精品一区二区| 看十八女毛片水多多多| 日韩欧美精品免费久久| 在线天堂最新版资源| 毛片女人毛片| 亚洲国产高清在线一区二区三| 99久久精品国产国产毛片| 九九在线视频观看精品| 国产精品熟女久久久久浪| 纵有疾风起免费观看全集完整版| 熟女人妻精品中文字幕| 中文欧美无线码| 内地一区二区视频在线| 国产伦精品一区二区三区视频9| 99热这里只有精品一区| 一区二区av电影网| 精品亚洲成国产av| 少妇的逼好多水| 中文资源天堂在线| 777米奇影视久久| 人妻 亚洲 视频| 日韩强制内射视频| 91精品一卡2卡3卡4卡| 青青草视频在线视频观看| 五月开心婷婷网| 精品久久久久久久末码| 亚洲欧美一区二区三区黑人 | 一区二区av电影网| 在线观看一区二区三区激情| 99久久中文字幕三级久久日本| 国产一区有黄有色的免费视频| 日本vs欧美在线观看视频 | 成人美女网站在线观看视频| 亚洲一级一片aⅴ在线观看| 老女人水多毛片| av又黄又爽大尺度在线免费看| 尤物成人国产欧美一区二区三区| 欧美bdsm另类| 亚洲精品国产av蜜桃| 精品国产三级普通话版| 国产91av在线免费观看| 视频区图区小说| 精品亚洲成国产av| 少妇的逼好多水| 人人妻人人看人人澡| 亚洲精品国产色婷婷电影| 小蜜桃在线观看免费完整版高清| 日本猛色少妇xxxxx猛交久久| av专区在线播放| 国产高清国产精品国产三级 | 亚洲内射少妇av| 日韩一区二区视频免费看| 国产精品三级大全| 国产在线男女| 一级二级三级毛片免费看| 亚洲欧美日韩东京热| 97热精品久久久久久| 美女脱内裤让男人舔精品视频| 国语对白做爰xxxⅹ性视频网站| 九草在线视频观看| 青青草视频在线视频观看| 伦精品一区二区三区| 校园人妻丝袜中文字幕| 日韩国内少妇激情av| 日韩强制内射视频| 欧美老熟妇乱子伦牲交| 日本av手机在线免费观看| 黄色配什么色好看| 成人无遮挡网站| 美女脱内裤让男人舔精品视频| 一本一本综合久久| 成人毛片a级毛片在线播放| 黄色视频在线播放观看不卡| 肉色欧美久久久久久久蜜桃| 美女cb高潮喷水在线观看| 中文字幕av成人在线电影| 一级黄片播放器| 日韩av在线免费看完整版不卡| 秋霞在线观看毛片| 美女主播在线视频| 国产在线一区二区三区精| 亚洲精品色激情综合| 午夜视频国产福利| 少妇裸体淫交视频免费看高清| 国产欧美亚洲国产| 午夜免费鲁丝| 波野结衣二区三区在线| 亚洲精品国产av成人精品| 亚洲精品乱码久久久久久按摩| 嘟嘟电影网在线观看| 成人黄色视频免费在线看| 亚洲精品456在线播放app| 草草在线视频免费看| av在线老鸭窝| 国国产精品蜜臀av免费| 国产探花极品一区二区| 亚洲色图综合在线观看| 一级毛片黄色毛片免费观看视频| 亚洲国产精品成人久久小说| 男人舔奶头视频| 中文字幕亚洲精品专区| 成人国产av品久久久| 韩国av在线不卡| 少妇被粗大猛烈的视频| av在线老鸭窝|