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

    CO2水合物漿作為空調(diào)載冷劑的流動(dòng)和傳熱特性研究進(jìn)展

    2014-07-05 16:02:45謝應(yīng)明謝振興范興龍
    化工進(jìn)展 2014年1期
    關(guān)鍵詞:冷劑漿體水合物

    謝應(yīng)明,謝振興,范興龍

    (上海理工大學(xué)能源與動(dòng)力工程學(xué)院,上海 200093)

    進(jìn)展與述評(píng)

    CO2水合物漿作為空調(diào)載冷劑的流動(dòng)和傳熱特性研究進(jìn)展

    謝應(yīng)明,謝振興,范興龍

    (上海理工大學(xué)能源與動(dòng)力工程學(xué)院,上海 200093)

    綜述了國(guó)內(nèi)外在CO2水合物漿的流動(dòng)和傳熱特性方面的最新研究進(jìn)展。在流動(dòng)特性部分,介紹了流體的五大類型(牛頓流體、假塑性流體、膨脹性流體、賓翰流體和赫-巴流體)以及CO2水合物漿的表觀黏度計(jì)算公式,研究表明CO2水合物漿是非牛頓流體,表觀黏度較小,因而具有優(yōu)秀的流動(dòng)性能;在傳熱特性部分,介紹了水合物漿在板式換熱器中的總體換熱系數(shù)計(jì)算方法,研究表明CO2水合物漿的對(duì)流換熱系數(shù)高達(dá)3658 W/(m2·K),因此有優(yōu)秀的傳熱性能,而且傳熱性能會(huì)受到晶體直徑和分布規(guī)律、流速、固體含量和管道尺寸等因素的影響。綜上所述,CO2水合物漿具有良好的流動(dòng)和傳熱特性,在空調(diào)領(lǐng)域有著廣闊的應(yīng)用前景。最后簡(jiǎn)要展望了今后的研究重點(diǎn)。

    二氧化碳;水合物漿;多相流;空調(diào)載冷劑;流動(dòng);傳熱

    《蒙特利爾議定書(shū)》對(duì)氟里昂制冷劑,包括氯氟烴(CFCs)和氫氯氟烴(HCFCs)等作了嚴(yán)格的管制規(guī)定,因?yàn)樗鼈兇嬖诔粞鯎p耗的潛能。為了避免這些制冷劑對(duì)臭氧的影響,氫氟烴(HFCs)被提出來(lái)作為CFCs和HCFCs的替代物,但是HFCs卻

    具有高的溫室效應(yīng),受到《京都議定書(shū)》的約束。隨著溫室氣體的減排進(jìn)程加快,HFCs制冷劑逐漸被淘汰是不可避免的。許多替代制冷劑,比如CO2、NH3、R290等自然工質(zhì),HFO-1234yf、HFE143m等新型環(huán)保制冷劑開(kāi)始得到重視[1-6]。然而,上述制冷劑或多或少存在缺陷,如毒性、可燃性、與潤(rùn)滑油的兼容性差和金屬腐蝕性等。在這樣的背景下,間接制冷,也稱二次制冷(secondary refrigeration),是一種比較理想的制冷替代方法[7-9]。實(shí)際上,二次制冷的研究主要在于環(huán)保載冷劑的研究。載冷劑作用是攜帶冷量并把冷量從制冷裝置運(yùn)輸?shù)叫枥淇臻g,這樣可以減小制冷劑的使用量和泄漏量[10]。常用的載冷劑有水、冰和共晶鹽。水作為載冷劑是利用顯熱蓄冷,從而導(dǎo)致體積龐大的蓄水箱和大功率的水泵;冰漿是水和冰組成的兩相載冷劑,其采取潛熱蓄冷,相比水具有更小的蓄冷體積,但是制冰過(guò)程中過(guò)低的蒸發(fā)溫度和過(guò)大的過(guò)冷度使得制冰漿效率低、能耗大;共晶鹽的相變潛熱低,時(shí)間長(zhǎng)了容易老化失效、腐蝕裝置。四丁基溴化銨(TBAB)、四氫呋喃(THF)水合物漿作為載冷劑已經(jīng)受到廣泛關(guān)注,被認(rèn)為是替代水的理想介質(zhì)選擇[11-12],但是TBAB和THF具有一定的毒性和腐蝕性,且價(jià)格沒(méi)有優(yōu)勢(shì)。顯然,尋找一種理想的載冷劑成為關(guān)鍵,而新型環(huán)保的CO2水合物漿就是一種選擇。CO2水合物的相變潛熱為500 kJ/kg,比冰的潛熱(333 kJ/kg)大,并且生成溫度可以在0~10 ℃變化,生成壓力在添加劑的作用下可降低70%~90%,低于1 MPa[8-9,13-14],非常適合用在蓄冷或空調(diào)領(lǐng)域。因此在生產(chǎn)和運(yùn)輸環(huán)節(jié),CO2水合物漿比冰漿更加高效和節(jié)能。為了保證漿體所攜帶的冷量滿足要求,CO2水合物漿中的固體含量必須適量[15-16]。固體含量過(guò)少,冷量不夠;固體含量多了,流動(dòng)性降低,泵耗功大,嚴(yán)重時(shí)會(huì)堵塞管路。嚴(yán)格來(lái)說(shuō),CO2水合物漿在管路中流動(dòng)時(shí)包含水、CO2水合物固體和微小CO2氣泡,是一種多相混合物[17]。CO2水合物漿能實(shí)際應(yīng)用的前提是具有良好的換熱和流動(dòng)性能,因此對(duì)其流動(dòng)特性和傳熱特性進(jìn)行研究和優(yōu)化至關(guān)重要。

    1 CO2水合物漿

    CO2水合物漿是一種類似冰漿、主要由細(xì)小的CO2水合物顆粒和水組成的多相懸浮液。單純的CO2+H2O體系生成的水合物為Ⅰ型水合物,水合物結(jié)構(gòu)包含8個(gè)晶穴,其中2個(gè)小晶穴、6個(gè)大晶穴[18-20]。如果所有的晶穴都被CO2分子占據(jù)的話,那么8個(gè)CO2分子對(duì)應(yīng)著46個(gè)水分子。CO2水合物的生成是結(jié)晶放熱過(guò)程,當(dāng)CO2+H2O體系生成Ⅰ型水合物時(shí),理想狀態(tài)下每消耗1 mol CO2將有46/8=5.75 mol的水參與反應(yīng)過(guò)程。但是,氣體水合物是非化學(xué)計(jì)量絡(luò)合物,水的物質(zhì)的量n與生成時(shí)的壓力和溫度相關(guān)。根據(jù)大小晶穴的填充率不同,水的物質(zhì)的量n取值也不同,見(jiàn)表1[21]。

    表1 CO2水合物晶穴占據(jù)比例

    表1中水的物質(zhì)的量為5.75 mol表示晶穴全部被占據(jù),為理想狀態(tài);水的物質(zhì)的量≥9.2 mol,說(shuō)明反應(yīng)的氣體CO2質(zhì)量少。根據(jù)文獻(xiàn)[22],在實(shí)驗(yàn)條件下,水的物質(zhì)的量的值一般近似取7.23。即CO2水合物漿生成方程可表述為式(1)。

    Hu等[23]研究飽和CO2水溶液在過(guò)冷度為1~9 ℃時(shí)水合物生成情況。他們觀察到清澈的CO2水溶液不到2 s變成渾濁溶液,這是由于水中生成許多微小的水合物晶體造成的。此現(xiàn)象很符合水合物快速成核理論,水合物成核與生長(zhǎng)類似于鹽的結(jié)晶過(guò)程。

    生成CO2水合物漿的方式有攪拌式、噴霧式、鼓泡式和噴射誘導(dǎo)式等多種方式。圖1為CO2水合物在無(wú)添加劑和有添加劑四氫呋喃(THF)的相平衡曲線圖[13]。從圖1中易看出,無(wú)添加劑時(shí),水合物生成溫度隨著壓力的升高而急劇升高。當(dāng)生成壓力過(guò)高時(shí),對(duì)于蓄冷或載冷劑等空調(diào)系統(tǒng)不利,因?yàn)檫^(guò)高的壓力導(dǎo)致成本上升、安全系數(shù)降低,而且容易泄漏。但是當(dāng)有添加劑后,平衡壓力大大降低,一般降低到原來(lái)的70%~90%。THF質(zhì)量分?jǐn)?shù)高于10%,生成壓力在12 ℃左右都能保持在1 MPa以下,因此可以根據(jù)空調(diào)系統(tǒng)的要求調(diào)節(jié)生成參數(shù),制備出滿足要求的漿體。

    圖1 CO2水合物相平衡曲線圖

    2 流動(dòng)特性

    2.1 流體的劃分

    在氣體水合物漿的研究中,漿體的流動(dòng)特性關(guān)系到其在管道中能否被高效的輸運(yùn)。一般說(shuō)來(lái),流體流動(dòng)特性可描述為剪切應(yīng)力τ與剪切率.γ的函數(shù),見(jiàn)式(2)。

    根據(jù)屈服應(yīng)力τo和流變指數(shù)n的值不同,可以把流體大致分為5類[24],如表2所示。表中k為黏度系數(shù)。

    根據(jù)表2,這5種流體的流動(dòng)特征變化趨勢(shì)可由圖2表示。

    2.2 表觀黏度計(jì)算征時(shí)采用奧氏黏度法(Ostwald viscosimeter method)對(duì)流體進(jìn)行計(jì)算,見(jiàn)式(3)~式(6)。最后,CO2水合物漿的表觀黏度就能定義為式(7)。

    表2 流體種類劃分

    圖2 流體特征曲線

    當(dāng)前,大多數(shù)研究者在研究水合物漿的流動(dòng)特

    式中,Q為體積流量,m3/s;τw為壁面剪切應(yīng)力,Pa;.

    γ為剪切率;D為圓管直徑,m;ΔP為壓降,Pa;L為壓降兩端距離,m;u為流體速度,m/s;n為流變指數(shù)。

    計(jì)算水合物漿體的流動(dòng)特征研究中以上是目前使用最多的方法。據(jù)不同研究者[25-28]研究表明,水合物漿是非牛頓流體,即使在很低的水合物質(zhì)量分?jǐn)?shù)的情況下依然保持這種特性,恰恰與呈現(xiàn)牛頓流體特征的冰漿相反[29]。因此,為了探究CO2水合物漿流體特征,國(guó)內(nèi)外研究者紛紛對(duì)其展開(kāi)了研究。

    2.3 研究現(xiàn)狀

    法國(guó)的Delahaye等[15]通過(guò)注射的方式生成CO2水合物漿,并考察了CO2水合物漿流動(dòng)特性,他們認(rèn)為漿體的黏度與其水合物體積分?jǐn)?shù)φs密切相關(guān)。當(dāng)φs<5%時(shí)為膨脹性流體;5%≤φs<10%為赫-巴流體(膨脹性趨勢(shì));φs= 10%時(shí)為賓翰流體;φs> 10%時(shí)為赫-巴流體(假塑性趨勢(shì))。此后他們還研究了有防聚劑Caflond 情況下CO2水合物漿的流變特性,實(shí)驗(yàn)表明CO2水合物漿呈現(xiàn)出牛頓流體的特征。正如作者所說(shuō)的那樣,防聚劑使?jié){體更加均勻,水合物晶體更加細(xì)小,從而有著較低的表觀黏度。隨后Jerbi等[30]通過(guò)攪拌制備CO2水合物漿,實(shí)驗(yàn)表明漿體的表觀黏度隨著水合物含量的增加而升高,并且是非牛頓流體,這些很好地與之前的研究相符合;他們還發(fā)現(xiàn)攪拌生成的水合物漿的黏度很低,這很可能是攪拌使得CO2水合物與水充分混合,降低水合物晶體凝聚,從而降低了流動(dòng)阻力。國(guó)內(nèi),陳偉軍[31]研究表明,在CO2水合物漿的固體含量為2.75%~6.34%時(shí),CO2水合物漿在TBAB作用下屬于假塑性流體,并且其表觀黏度略大于純的CO2水合物漿,而且比TBAB水合物漿的表觀黏度小得多。從文獻(xiàn)[8]可知,TBAB能使生成的CO2水合物變成半籠形型水合物,相比Ⅰ型,具有較大的晶體體積,導(dǎo)致了漿體黏度上升。戴海鳳[32]研究認(rèn)為,當(dāng)漿體在管壁面的剪切率在 610 s-1< < 640 s-1的范圍時(shí),表觀黏度較小,并且存在最優(yōu)表觀黏度,使得漿體流動(dòng)摩擦損失最低。

    然而,瑞士的Hu等[23,33]的研究結(jié)果卻與之前的研究成果有差異。他們得出的結(jié)論為:CO2飽和水溶液生成的水合物漿中的水合物固體含量對(duì)漿體的黏度影響不明顯;密度大于1021 kg/m3(對(duì)應(yīng)的固體含量為34%)時(shí),漿體黏度迅速增大,但即使固體含量高達(dá)40%,其黏度還保持較小的值,為3.4 mPa·s,比相同質(zhì)量濃度冰漿的黏度更小[34],因此表現(xiàn)出了優(yōu)秀的流動(dòng)性能和穩(wěn)定性。但與此同時(shí),他們還發(fā)現(xiàn)較高固體含量的漿體在低流速時(shí)具有較高的壓力降,并把造成這種現(xiàn)象的原因歸咎于水合物晶體與管道壁面之間的摩擦力。以上研究者的成果歸納總結(jié)如表3所示。

    從表3中可以看出,盡管體積分?jǐn)?shù)各有不同,但CO2水合物漿的表觀黏度保持在較低的水平,因此有著優(yōu)秀的流動(dòng)性能。

    再者,肖睿等[36]在研究TBAB水合物漿時(shí)發(fā)現(xiàn)流速超過(guò)一定值后,會(huì)出現(xiàn)“再層流化”現(xiàn)象,固相含量越大越有利于維持流動(dòng)的層流狀態(tài),抑制流動(dòng)進(jìn)入湍流狀態(tài),并且漿體的流動(dòng)會(huì)出現(xiàn)“阻力降低區(qū)”。根據(jù)以上論斷,漿體輸送中選取好流速之后應(yīng)該采用合適的固相含量,以盡可能使流動(dòng)進(jìn)入“阻力降低區(qū)”,最大限度降低漿體的管道輸送阻力。在載冷劑循環(huán)中,載冷劑的作用是把冷量送到需冷的場(chǎng)所,輸送距離遠(yuǎn),輸送過(guò)程中應(yīng)該盡量減少泵耗功和冷量損失。固體含量過(guò)低則輸送冷量小,而且易進(jìn)入湍流區(qū),增加了流動(dòng)阻力;過(guò)高時(shí),雖然輸送冷量大了,也容易保證層流流態(tài),但黏度迅速增大,流動(dòng)阻力也會(huì)迅速攀升。他們認(rèn)為出現(xiàn)再層流化現(xiàn)象的原因主要在于漿體中的固相粒子對(duì)液相湍流動(dòng)能具有吸收作用。這種現(xiàn)象是否也在CO2水合物漿中出現(xiàn),還需進(jìn)一步的實(shí)驗(yàn)研究。

    3 傳熱特性

    傳熱特性的研究是有效利用水合物漿的一個(gè)重要環(huán)節(jié)。研究表明在相同的雷諾數(shù)條件下,相比于單相顯熱傳熱,漿體中固相顆粒的相變潛熱能夠提高傳熱效率,并且隨著固體含量的增大而提高,比如冰漿和TBAB水合物漿的換熱特性研究[37-39]中,固相顆粒使得近壁處傳熱邊界層擾動(dòng)從而提高換熱系數(shù)。因此,漿體中的固體含量是影響傳熱性能的一個(gè)重要因素[40]。

    3.1 換熱系數(shù)計(jì)算

    表3 CO2水合物漿的流體特性

    雖然有關(guān)CO2水合物漿的傳熱特性研究國(guó)內(nèi)外都比較少,但是可以參考其他漿體的傳熱研究成果,進(jìn)行預(yù)測(cè)和對(duì)比。參照文獻(xiàn)[39],作者根據(jù)大量傳熱實(shí)驗(yàn)結(jié)果和經(jīng)驗(yàn)公式總結(jié)出了適合水合物漿傳熱計(jì)算的公式。此計(jì)算公式主要用來(lái)對(duì)水合物漿在板式換熱器中與熱水進(jìn)行換熱計(jì)算,因此對(duì)于CO2水合物漿具有一定的參考意義??倱Q熱系數(shù)h總先通過(guò)實(shí)驗(yàn)測(cè)出,最后與理論計(jì)算進(jìn)行校核??倱Q熱系數(shù)h總[W/(m2·K)]和換熱器兩側(cè)(h熱側(cè),h冷側(cè))(熱側(cè)為熱水側(cè),冷側(cè)為水合物漿側(cè))的換熱系數(shù)關(guān)系為式(8)。

    式中,δ為板式換熱器的壁厚,m;λ為板式換熱器的熱導(dǎo)率,W/(m·K)。

    換熱器兩側(cè)的局部換熱系數(shù)h熱側(cè)和h冷側(cè)可通過(guò)關(guān)系式Nu=aRebPrc分別計(jì)算出來(lái)。此關(guān)系式中,一般取普朗克指數(shù)c=1/3,因此式(8)變?yōu)槭剑?)。

    式中,D為換熱管徑,m;常數(shù)a和b的值可通過(guò)迭代的方法計(jì)算出。設(shè)a,b∈(0,1),步長(zhǎng)取為0.01代入式(9)進(jìn)行迭代計(jì)算,在誤差允許的范圍得出a和b分別為0.49和0.63,得出熱側(cè)的努賽爾方程為式(10)。

    由于水合物漿吸熱分解是相變過(guò)程,顆粒不斷融化,固體含量不斷降低,導(dǎo)致雷諾數(shù)Re和普朗特?cái)?shù)Pr等換熱相關(guān)參數(shù)變化較大,因此需要對(duì)冷側(cè)的方程進(jìn)行修正。修正關(guān)系式為式(11)。同獲得式(10)的方法,最終冷側(cè)的修正努賽爾方程為式(12)。

    式中,Δω為質(zhì)量濃度變化,%;ΔH為水合物相變潛熱,kg/kJ;Cp為熱容量,kJ/(kg·K);ΔT為對(duì)數(shù)平均溫差,K。則在平均對(duì)數(shù)溫差ΔT范圍內(nèi),ΔωΔH表示水合物總的分解熱,CpΔT表示水合物的顯熱。最后結(jié)合式(8)、式(10)和式(11)就可以算出總的換熱系數(shù)h總。研究表明,上述計(jì)算過(guò)程能很好地用于水合物漿的換熱計(jì)算,誤差保持在5%左右。

    3.2 研究現(xiàn)狀

    Hu等[23]用套管加熱質(zhì)量分?jǐn)?shù)為15%的CO2水合物漿來(lái)計(jì)算理想狀態(tài)下漿體的對(duì)流換熱系數(shù)。實(shí)驗(yàn)結(jié)果表明內(nèi)管中的CO2水合物漿與管外逆流的水發(fā)生換熱,漿體側(cè)的對(duì)流換熱系數(shù)高達(dá)3658 W/(m2·K)。戴海鳳[32]實(shí)驗(yàn)測(cè)出了在靜態(tài)水的加熱下,質(zhì)量分?jǐn)?shù)為13.2%、流速為0.45 m/s的CO2水合物漿對(duì)流換熱系數(shù)為1800 W/(m2·K)。造成如此大的差異可能是因?yàn)椴煌膶?shí)驗(yàn)條件和計(jì)算方法導(dǎo)致的,因此更多實(shí)驗(yàn)研究工作需開(kāi)展起來(lái)。

    馬志偉等[41]研究了TBAB水合物漿的換熱特性,發(fā)現(xiàn)隨著紊流程度的增大, 加熱功率和相變潛熱對(duì)漿體換熱產(chǎn)生的影響減弱,這時(shí)換熱情況主要受流動(dòng)的影響。再者,Stamatiou等[42]認(rèn)為水合物晶粒分布對(duì)于換熱特性影響非常明顯。Clarke等[43-44]采用聚焦光束反射(FBRM)的方法研究了CO2水合物粒度分布情況;Myriam等[45]發(fā)現(xiàn)TBAB水合物晶體的粒徑分布在10~100 μm之間;Delahaye 等[46]探究了冰漿和THF水合物漿的晶體粒度分布規(guī)律及其晶粒大小。張鵬等[39]研究了TBAB水合物漿在板式換熱器中的換熱性能,他們發(fā)現(xiàn)局部換熱系數(shù)主要與漿體流速有關(guān),與漿體的固體含量關(guān)系較?。划?dāng)漿體變成TBAB溶液后,由于較低的黏度和較高的雷諾數(shù),局部換熱系數(shù)明顯提高;在恒定的流速下,總換熱系數(shù)隨著固體含量的提高增大。水合物晶粒大小有可能直接影響漿體的傳熱特性,更小的晶粒尺寸使得漿體更加均勻,傳熱效率更高,黏度更小,流動(dòng)性能更好等。表4歸納總結(jié)了不同種類晶粒分布研究結(jié)果。從表4中可以看出,CO2水合物晶體的平均直徑最小,因此有著更優(yōu)越的流動(dòng)和傳熱特性。

    綜上所述,晶體分布規(guī)律及大小、漿體的流速、固體含量和管道尺寸等,甚至氣態(tài)CO2都會(huì)對(duì)CO2水合物漿的傳熱特性產(chǎn)生影響。特別是水合物溶解后產(chǎn)生CO2氣體,這時(shí)氣-固-液3種狀態(tài)同時(shí)存在,使得傳熱特性研究就變得更加復(fù)雜。

    表4 不同晶粒分布的研究

    4 結(jié) 語(yǔ)

    詳細(xì)介紹了CO2水合物漿的流動(dòng)及其傳熱特性研究概況。CO2水合物漿作為一種新型環(huán)保節(jié)能的載冷劑,有著相變潛熱大、制備簡(jiǎn)單、流動(dòng)和傳熱性能良好等特點(diǎn),用在載冷循環(huán)中,相比于其他載冷漿體更加經(jīng)濟(jì)。但是,CO2水合物漿的流動(dòng)和傳熱等基礎(chǔ)性研究還十分欠缺。

    為了使CO2水合物漿技術(shù)真正實(shí)用化,今后的主要工作有以下幾方面。

    (1)研究快速、連續(xù)、循環(huán)制備CO2水合物漿的技術(shù),進(jìn)一步提高制備效率。

    (2)在實(shí)驗(yàn)和理論的基礎(chǔ)上對(duì)漿體流動(dòng)和傳熱過(guò)程進(jìn)行數(shù)值模擬,優(yōu)化流動(dòng)和傳熱特性。

    (3)換熱器中分離氣態(tài)CO2,增強(qiáng)換熱效果,并且使分離的CO2氣體重新生成水合物進(jìn)行循環(huán)。

    (4)防止水合物顆粒集聚現(xiàn)象[16,48], CO2水合物漿的流動(dòng)時(shí)間長(zhǎng)了,水合物晶體容易附聚結(jié)塊,不利于漿體的流動(dòng)和傳熱。

    [1] 郝紅,熊國(guó)華,徐蕾. 制冷劑研究現(xiàn)狀及發(fā)展動(dòng)向[J]. 化工進(jìn)展,2001,20(3):5-6,12.

    [2] 謝應(yīng)明,梁德青,郭開(kāi)華,等. 氣體水合物蓄冷技術(shù)研究現(xiàn)狀與展望[J]. 暖通空調(diào),2004,34(9):25-29.

    [3] 李連生. 制冷劑替代技術(shù)研究進(jìn)展及發(fā)展趨勢(shì)[J]. 制冷學(xué)報(bào),2011,32(6):53-58.

    [4] 王鑫,史琳,朱明善,等. 近年來(lái)國(guó)內(nèi)外制冷劑的研究狀況[J]. 暖通空調(diào),2007,37(10):40-43.

    [5] 高鳳龍. 從全球氣候變化看制冷劑的替代[J]. 制冷與空調(diào),2010,24(3):115-118.

    [6] 李剛,謝應(yīng)明,劉道平. 蓄冷空調(diào)新工質(zhì)——?dú)怏w水合物的替代選擇概述[J]. 制冷學(xué)報(bào),2008,29(3):18-24.

    [7] Jerbi S,Delahaye A,F(xiàn)ournaison L,et al. Characterization of CO2hydrate formation and dissociation kinetics in a flow loop[J]. Int. J. Refrig.,2010,33(8):1625-1631

    [8] Lin W,Delahaye A,F(xiàn)ournaison L. Phase equilibrium and dissociation enthalpy for semi-clathrate of CO2+TBAB[J]. Fluid Phase Equilibri.,2008,264(1-2):220-227.

    [9] Marinhas S,Delahaye A,F(xiàn)ournaison L,et al. Modeling of the available latent heat of a CO2hydrate slurry in an experimental loop applied to secondary refrigeration[J]. Chemical Engineering and Processing,2006,45(3):184-192.

    [10] Martinez M C,Dalmazzone D,F(xiàn)ust W,et al. Thermodynamic properties of THF+CO2hydrates in relation with refrigeration applications[J]. AIChE J.,2008,54(4):1088-1095.

    [11] 馬志偉,張鵬,王如竹. 水合物漿體蓄冷空調(diào)系統(tǒng)的運(yùn)行特性[J].化工學(xué)報(bào),2010,61(S2):93-98.

    [12] 王武昌,李玉星,樊栓獅,等. 四氫呋喃水合物漿流動(dòng)特性[J]. 化工進(jìn)展,2010,29(8):1418-1422

    [13] Delahaye A,F(xiàn)ournaison L,Marinhas S,et al. Effect of THF on equilibrium pressure and dissociation enthalpy of CO2hydrates applied to secondary refrigeration[J]. Ind. Eng. Chem. Res.,2006,45 (1):391-397.

    [14] 姚豪,周春艷,梁德清,等. 水合物漿和冰漿高密度潛熱輸送研究進(jìn)展[J]. 化工學(xué)報(bào),2003,54(S1):57-61.

    [15] Delahaye A,F(xiàn)ournaison L,Marinhas S,et al. Rheological study of CO2hydrate slurry in a dynamic loop applied to secondary Refrigeration[J]. Chem. Eng. Sci. Res.,2008,63(13):3551-3559.

    [16] Delahaye A,F(xiàn)ournaison L,Jerbi S,et al. Rheological properties of CO2hydrate slurry flow in the presence of additives[J]. Ind. Eng. Chem. Res.,2011,50(13):8344-8353.

    [17] Sari O,Hu J,Brun F,et al. Proceedings of the 22nd International Congress of Refrigeration[C]. Beijing:International Institute of Refrigeration,2007.

    [18] 祁影霞,張華. 添加水合物促進(jìn)二氧化碳水合物生成的實(shí)驗(yàn)研究[J]. 高?;瘜W(xué)工程學(xué)報(bào),2010,24(5):842-847.

    [19] 王趙,張偉,李文強(qiáng),等. CO2水合物氣體分子籠占據(jù)狀態(tài)的第一性原理研究[J]. 四川師范大學(xué)學(xué)報(bào):自然科學(xué)版,2010,33(3):356-360.

    [20] Marinhas S,Delahaye A,F(xiàn)ournaison L. Solid fraction modeling for CO2and CO2-THF hydrate slurries used as secondary refrigerants[J]. Int.J.Refrig., 2007,30(5):758-766.

    [21] Teng H,Yamasaki A,Shindo Y. Stability of the hydrate layer formed on the surface of a CO2droplet in high-pressure,low-temperature water[J]. Chem. Eng. Sci.Res.,1996,51(22):4979-4986.

    [22] Kang S P,Lee H,Lee C S,et al. Hydrate phase equilibria of guest mixtures containing CO2,N2and tetrahydrofuran[J]. Fluid Phase Equilib.,2001,185(1-2):101-109.

    [23] Hu J,Sari O,Mahmed C,et al. 9th International Conference on Phase-Change Materials and Slurries for Refrigeration and Air Conditioning[C]. Sofia,Bulgaria:International Institute of Refrigeration,2010.

    [24] 何歸庭. 工程流體力學(xué)[M]. 北京:科學(xué)出版社,2007:17-18.

    [25] Xiao R,Wu S,Tang L,et al. Proceedings of 10th International Conference on thermal energy storage[C]. USA,New Jersey:Ecostock,2006.

    [26] Wang W,F(xiàn)an S,Liang D,et al. Experimental study on flow characters of CH3CCl2F hydrate slurry[J]. Int. J. Refrig.,2008,31 (3):371-378.

    [27] Clain P,Delahaye A,F(xiàn)ournaison L,et al. Rheological properties of tetra-n-butylphosphonium bromide hydrate slurry flow[J]. Chem. Eng. J.,2012,193:112-122.

    [28] Kumano H,Hirata T,Kudoh T. Experimental study on the flow and heat transfer characteristics of a tetra-n-butyl ammonium bromide hydrate slurry[J]. Int. J. Refrig.,2011,34(8):1963-1971.

    [29] Youssef A Z,Delahaye A,Huang L,et al. State of the art on phase change material slurries[J]. Energy Conversion and Management,2013,65(S1):120-132.

    [30] Jerbi S,Delahaye A,Oignet J,et al. Rheological properties of CO2hydrate slurry produced in a stirred tank reactor and a secondary refrigeration loop[J]. Int. J. Refrig.,2013,36(4):1294-1301.

    [31] 陳偉軍. CO2水合物漿的連續(xù)制備及流動(dòng)特性實(shí)驗(yàn)研究[D]. 上海:上海理工大學(xué),2012.

    [32] 戴海鳳. CO2水合物漿流動(dòng)傳熱特性的研究[D]. 上海:上海理工大學(xué),2013.

    [33] Sari O,Hu J,Eicher S,et al. International Refrigeration and Air Conditioning Conference[C]. Purdue:Purdue University,July 14-17,2008.

    [34] Kitanovski A,Poredos A. Concentration distribution and viscosity of ice-slurry in heterogeneous flow[J]. Int. J. Refrig.,2002,25(6):827-835.

    [35] Oyama H,Ebinuma T,Shimada W,et al. 4th Int. Conf. on Gas Hydrates[C]. Yokohama,Japan,2002.

    [36] 肖睿,宋文吉,黃沖,等. TBAB包絡(luò)化合物漿的管內(nèi)流動(dòng)再層流化現(xiàn)象[J]. 工程熱物理學(xué)報(bào),2009,30(6):971-973.

    [37] Ma Z W,Zhang P,Wang R Z,et al. Forced flow and convective melting heat transfer of clathrate hydrate slurry in tubes[J]. Int. J. Heat Mass Transf .,2010,53(19-20):3745-3757.

    [38] 宋文吉,肖睿,董凱軍,等. TBAB包絡(luò)化合物漿水平管內(nèi)的傳熱特性研究[J]. 工程熱物理學(xué)報(bào),2009,30(9):1502-1504.

    [39] Ma Z W,Zhang P. Pressure drop and heat transfer characteristics of clathrate hydrate slurry in a plate heat exchanger[J]. Int. J. Refrig.,2011,34(3):796-806.

    [40] Zhai X Q,Wang X L,Wang T,et al. A review on phase change cold storage in air-conditioning system:Materials and applications[J]. Renewable and Sustainable Energy Reviews,2013,22:108-120.

    [41] 馬志偉,張鵬,王如竹. 水合物漿體在直管中的流動(dòng)換熱特性研究[J]. 工程熱物理學(xué)報(bào),2010,31(8):1398-1402.

    [42] Stamatiou E,Meewisse J W,Kawaji M. Ice slurry generation involving moving parts[J]. Int. J. Refrig.,2005,28(1):60-72.

    [43] Clarke M A,Bishnoi P R. Determination of the intrinsic rate constant and activation energy of CO2gas hydrate decomposition using in-situ particle size analysis[J]. Chem. Eng. Sci.,2004,59(14):2983-2993.

    [44] Clarke M A,Bishnoi P R. Determination of the intrinsic kinetics of CO2gas hydrate formation using in situ particle size analysis[J]. Chem. Eng. Sci.,2005,60(3):695-709.

    [45] Myriam D,Michel C,Jean-Michel H. Rheological study of TBAB hydrate slurries as secondary two-phase refrigerants[J]. Int. J. Refrig.,2005,28(5):663-671.

    [46] Delahaye A,F(xiàn)ournaison L,Guilpart J. Characterization of ice and THF hydrate slurry crystal size distribution by microscopic observation method[J]. Int. J. Refrig.,2010,33(8):1639-1647.

    [47] Inada T,Modak P R. Growth control of ice crystals by poly(vinyl alcohol) and antifreeze protein in ice slurries[J]. Chem. Eng. Sci.,2006,61(10):3149-3158.

    [48] Huo Z,F(xiàn)reer E,Lamar M,et al. Hydrate plug prevention by anti-agglomeration[J]. Chem. Eng. Sci.,2001,56(17):4979-4991.

    Review on flow and heat transfer properties of CO2hydrate slurry as secondary refrigerant

    XIE Yingming,XIE Zhenxing,F(xiàn)AN Xinglong
    (School of Energy and Power Engineering,University of Shanghai for Science and Technology,Shanghai 200093,China)

    The recent development of flow and heat transfer research of CO2hydrate slurry was reviewed as a secondary refrigerant. In the aspect of flow research, five flow patterns(Newtonian flow, Pseudo-plastic flow, Dilatant flow, Bingham plastic flow and Herschel-Bulkley flow) were introduced. It shows that CO2hydrate slurry is a kind of non-Newtonian liquid with small apparent viscosity indicating excellent flow properties. With regard to heat transfer aspect, the calculation method of hydrate slurry overall heat transfer coefficient in the plate heat exchanger was presented, and CO2hydrate slurry has excellent heat transfer properties with a convective heat transfer coefficient of 3658W/(m2·K). Moreover, heat transfer properties can be affected by crystal diameter and distribution, flow velocity, solid content and tube dimensions. In conclusion, CO2hydrate slurry has good prospects in air conditioning application with good flow and heat transfer properties. Finally, the future research priorities on CO2hydrate slurry were discussed and summarized.

    carbon dioxide;hydrate slurry;multiphase flow;secondary refrigerant;flow;heat transfer

    TQ 025.3;TB 61+2;TU 831.6

    A

    1000-6613(2014)01-0010-07

    10.3969/j.issn.1000-6613.2014.01.002

    2013-07-15;修改稿日期:2013-08-26。

    國(guó)家自然科學(xué)基金(50806050)及上海市教委科研創(chuàng)新項(xiàng)目(14YZ097)。

    及聯(lián)系人:謝應(yīng)明(1976—),男,博士,副教授,研究方向?yàn)橹评淇照{(diào)技術(shù)、氣體水合物技術(shù)、新能源技術(shù)。E-mail xymbox@ 163.com。

    猜你喜歡
    冷劑漿體水合物
    天然氣液化工廠冷劑儲(chǔ)存區(qū)的布置淺析
    云南化工(2021年7期)2021-12-21 07:27:50
    氣井用水合物自生熱解堵劑解堵效果數(shù)值模擬
    漿體輸送中彎管磨損原因分析及預(yù)防措施
    昆鋼科技(2020年5期)2021-01-04 01:41:24
    熱水吞吐開(kāi)采水合物藏?cái)?shù)值模擬研究
    基于混料設(shè)計(jì)的MRC工藝中混合冷劑配比的優(yōu)化
    淺析吸收式制冷機(jī)組冷劑水故障問(wèn)題
    西部皮革(2018年2期)2018-02-17 06:51:01
    天然氣水合物保壓轉(zhuǎn)移的壓力特性
    我國(guó)海域天然氣水合物試采成功
    K+和Na+在C3S-納米SiO2漿體上的吸附和脫附特性
    長(zhǎng)距離漿體管道正排量泵智能協(xié)同系統(tǒng)的設(shè)計(jì)
    亚洲高清免费不卡视频| 三级国产精品片| 麻豆乱淫一区二区| 亚洲av电影不卡..在线观看| 大又大粗又爽又黄少妇毛片口| 精品久久久噜噜| 国产伦理片在线播放av一区| 国产免费一级a男人的天堂| 国产精品无大码| 久久久久久国产a免费观看| 禁无遮挡网站| 国产v大片淫在线免费观看| 天堂av国产一区二区熟女人妻| 国产白丝娇喘喷水9色精品| av播播在线观看一区| 亚洲av成人av| 午夜老司机福利剧场| 亚洲国产精品成人久久小说| 亚洲天堂国产精品一区在线| 国产激情偷乱视频一区二区| 国产精品久久久久久久电影| 久久99蜜桃精品久久| 国产欧美另类精品又又久久亚洲欧美| 国产精品久久视频播放| 免费看av在线观看网站| 久久久欧美国产精品| 成人毛片60女人毛片免费| 久久久久久伊人网av| 丝袜美腿在线中文| 麻豆精品久久久久久蜜桃| 亚洲av中文字字幕乱码综合| 1000部很黄的大片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精华霜和精华液先用哪个| 久久久精品欧美日韩精品| 麻豆成人午夜福利视频| 成人美女网站在线观看视频| 色网站视频免费| 久久精品国产自在天天线| av在线播放精品| 乱系列少妇在线播放| 国产日韩欧美在线精品| 肉色欧美久久久久久久蜜桃 | 美女内射精品一级片tv| 久久久久久久国产电影| 成人高潮视频无遮挡免费网站| 国产亚洲5aaaaa淫片| 久久久久久久午夜电影| 在线免费十八禁| 成人美女网站在线观看视频| 三级国产精品欧美在线观看| 日韩三级伦理在线观看| 久久国内精品自在自线图片| 99热这里只有精品一区| 九草在线视频观看| 少妇被粗大猛烈的视频| 69人妻影院| 成人二区视频| 99re6热这里在线精品视频| 亚洲成人精品中文字幕电影| 成人漫画全彩无遮挡| 久久精品夜色国产| 亚洲欧美成人精品一区二区| 国产黄色免费在线视频| 综合色av麻豆| 非洲黑人性xxxx精品又粗又长| 国产精品人妻久久久久久| 日韩av在线免费看完整版不卡| xxx大片免费视频| 久久国内精品自在自线图片| 夜夜爽夜夜爽视频| 人妻夜夜爽99麻豆av| 亚洲18禁久久av| 美女cb高潮喷水在线观看| 国产在视频线在精品| 久久久久久久久久人人人人人人| 国产成人91sexporn| 亚洲精品影视一区二区三区av| 国产久久久一区二区三区| 国产白丝娇喘喷水9色精品| 久久这里只有精品中国| 性插视频无遮挡在线免费观看| 免费黄网站久久成人精品| 岛国毛片在线播放| 免费播放大片免费观看视频在线观看| 卡戴珊不雅视频在线播放| 亚洲精品国产av蜜桃| 精品酒店卫生间| 久久久久精品久久久久真实原创| 午夜福利视频精品| 亚洲国产精品国产精品| 美女脱内裤让男人舔精品视频| 美女国产视频在线观看| 亚洲国产最新在线播放| 在线观看人妻少妇| 国产午夜福利久久久久久| 亚洲av免费在线观看| 亚洲四区av| 亚洲天堂国产精品一区在线| 九九在线视频观看精品| 99久久中文字幕三级久久日本| 最后的刺客免费高清国语| 午夜久久久久精精品| 午夜福利在线观看免费完整高清在| av又黄又爽大尺度在线免费看| 毛片女人毛片| 免费黄网站久久成人精品| 欧美区成人在线视频| 色吧在线观看| 最近手机中文字幕大全| 亚洲三级黄色毛片| 男女啪啪激烈高潮av片| 成人漫画全彩无遮挡| 晚上一个人看的免费电影| 黄色一级大片看看| 国产精品久久久久久精品电影| 久久久久九九精品影院| 久久这里有精品视频免费| 男人舔奶头视频| 国产精品爽爽va在线观看网站| 久久久久性生活片| 乱码一卡2卡4卡精品| 一个人看的www免费观看视频| 男人舔女人下体高潮全视频| 亚洲高清免费不卡视频| 午夜精品一区二区三区免费看| 国产女主播在线喷水免费视频网站 | 在线观看美女被高潮喷水网站| 国产精品伦人一区二区| 亚洲欧美日韩东京热| 久久99蜜桃精品久久| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品日本国产第一区| 美女被艹到高潮喷水动态| 亚洲av中文av极速乱| 我的女老师完整版在线观看| 大香蕉97超碰在线| 国产激情偷乱视频一区二区| 亚洲熟女精品中文字幕| av卡一久久| 伊人久久国产一区二区| 天堂中文最新版在线下载 | 天美传媒精品一区二区| 婷婷色av中文字幕| 国产精品.久久久| 久久久久久伊人网av| 在线免费十八禁| 国产伦在线观看视频一区| 热99在线观看视频| 国产精品久久视频播放| av专区在线播放| av在线观看视频网站免费| 天美传媒精品一区二区| 在线观看人妻少妇| 亚洲av国产av综合av卡| 国国产精品蜜臀av免费| 美女高潮的动态| 色综合色国产| 天堂√8在线中文| 亚洲精品第二区| 欧美激情久久久久久爽电影| 美女国产视频在线观看| 51国产日韩欧美| 91狼人影院| 久久久久精品久久久久真实原创| 在线 av 中文字幕| 免费人成在线观看视频色| 亚洲精品日韩在线中文字幕| 校园人妻丝袜中文字幕| 欧美xxxx黑人xx丫x性爽| 简卡轻食公司| 一级爰片在线观看| 中国美白少妇内射xxxbb| 国产精品伦人一区二区| 免费黄色在线免费观看| 成人午夜精彩视频在线观看| 三级男女做爰猛烈吃奶摸视频| 精华霜和精华液先用哪个| 国产男人的电影天堂91| 好男人在线观看高清免费视频| 亚洲成人一二三区av| 联通29元200g的流量卡| 亚洲欧美一区二区三区黑人 | 直男gayav资源| 日韩国内少妇激情av| 男女啪啪激烈高潮av片| 亚洲国产欧美人成| 毛片一级片免费看久久久久| 亚洲国产精品成人久久小说| 免费观看精品视频网站| 一级毛片黄色毛片免费观看视频| 男女下面进入的视频免费午夜| 乱人视频在线观看| 少妇熟女欧美另类| 亚洲三级黄色毛片| 亚洲欧美精品自产自拍| 久久久久精品久久久久真实原创| 嫩草影院入口| 久久久久久久国产电影| 少妇熟女aⅴ在线视频| 少妇人妻一区二区三区视频| 日韩av免费高清视频| 一级黄片播放器| 天天一区二区日本电影三级| av免费在线看不卡| 亚洲精品第二区| 波野结衣二区三区在线| 亚洲欧洲日产国产| 中文字幕免费在线视频6| 日本黄大片高清| 国产高清国产精品国产三级 | 国产单亲对白刺激| 熟女电影av网| 国产国拍精品亚洲av在线观看| 最近的中文字幕免费完整| 九草在线视频观看| 一级毛片电影观看| 欧美 日韩 精品 国产| 美女脱内裤让男人舔精品视频| 舔av片在线| 亚洲精品一区蜜桃| 亚洲精品自拍成人| 丰满少妇做爰视频| 国产亚洲精品av在线| www.色视频.com| 99热这里只有是精品在线观看| 久久久久久久久久久免费av| 美女黄网站色视频| 麻豆精品久久久久久蜜桃| 高清午夜精品一区二区三区| 中国国产av一级| 亚洲国产精品国产精品| 国产精品福利在线免费观看| 亚洲精品日本国产第一区| www.av在线官网国产| 亚洲国产欧美人成| 国产综合懂色| 亚洲成人精品中文字幕电影| 在线 av 中文字幕| 波野结衣二区三区在线| 青青草视频在线视频观看| 午夜精品一区二区三区免费看| 日韩中字成人| 国产黄色视频一区二区在线观看| 夫妻性生交免费视频一级片| 日本免费a在线| 日韩欧美国产在线观看| 日韩欧美精品v在线| 婷婷色综合www| 日韩成人伦理影院| 2018国产大陆天天弄谢| 国产精品国产三级国产专区5o| 精品熟女少妇av免费看| 在线观看免费高清a一片| 国产欧美日韩精品一区二区| 日本wwww免费看| 少妇人妻精品综合一区二区| 亚洲欧美精品自产自拍| 嫩草影院入口| 特级一级黄色大片| 日本猛色少妇xxxxx猛交久久| 成人亚洲精品av一区二区| 国产永久视频网站| 国产亚洲精品久久久com| 尾随美女入室| 男女边吃奶边做爰视频| 亚洲人成网站高清观看| 亚洲怡红院男人天堂| 麻豆av噜噜一区二区三区| 1000部很黄的大片| 日韩一区二区视频免费看| 日韩一本色道免费dvd| 免费观看无遮挡的男女| 国产激情偷乱视频一区二区| 久久99蜜桃精品久久| 国产高清有码在线观看视频| 午夜视频国产福利| 免费看a级黄色片| 一夜夜www| 美女被艹到高潮喷水动态| 国产成人免费观看mmmm| 晚上一个人看的免费电影| 午夜日本视频在线| 性插视频无遮挡在线免费观看| 国产伦在线观看视频一区| 国产成人福利小说| 十八禁国产超污无遮挡网站| 国产一区二区在线观看日韩| 狂野欧美白嫩少妇大欣赏| 欧美三级亚洲精品| 国精品久久久久久国模美| 小蜜桃在线观看免费完整版高清| 日韩欧美三级三区| 熟妇人妻不卡中文字幕| 国产永久视频网站| 日韩亚洲欧美综合| av天堂中文字幕网| 免费观看无遮挡的男女| 内射极品少妇av片p| 国产av不卡久久| 别揉我奶头 嗯啊视频| 久久国产乱子免费精品| 亚洲欧美一区二区三区国产| 免费观看av网站的网址| 男女视频在线观看网站免费| 熟妇人妻久久中文字幕3abv| 听说在线观看完整版免费高清| 成人美女网站在线观看视频| 免费电影在线观看免费观看| 九九爱精品视频在线观看| 亚洲精品自拍成人| 国产成人aa在线观看| 国产免费福利视频在线观看| 国产精品久久视频播放| 免费无遮挡裸体视频| 一级二级三级毛片免费看| 亚洲三级黄色毛片| 亚洲一级一片aⅴ在线观看| 午夜精品国产一区二区电影 | 99热网站在线观看| 成人亚洲精品av一区二区| 91aial.com中文字幕在线观看| 能在线免费观看的黄片| 两个人的视频大全免费| 老司机影院成人| 国产精品一区二区在线观看99 | 国产精品人妻久久久久久| 国产国拍精品亚洲av在线观看| 最近最新中文字幕大全电影3| 少妇丰满av| 天堂俺去俺来也www色官网 | 老司机影院成人| 国产有黄有色有爽视频| 好男人视频免费观看在线| 日韩,欧美,国产一区二区三区| 欧美一区二区亚洲| 97超碰精品成人国产| 少妇的逼好多水| 免费少妇av软件| 免费观看a级毛片全部| 久久人人爽人人片av| 两个人的视频大全免费| 成年人午夜在线观看视频 | 六月丁香七月| 色综合色国产| 卡戴珊不雅视频在线播放| 成人高潮视频无遮挡免费网站| 综合色丁香网| 国产午夜精品论理片| 色综合站精品国产| 成人欧美大片| 久久久国产一区二区| 亚洲在线自拍视频| 我的女老师完整版在线观看| 免费黄频网站在线观看国产| 免费不卡的大黄色大毛片视频在线观看 | 日本黄大片高清| 午夜福利视频1000在线观看| 精品熟女少妇av免费看| 男人和女人高潮做爰伦理| 99视频精品全部免费 在线| 国产成人a∨麻豆精品| 少妇熟女aⅴ在线视频| 一夜夜www| 不卡视频在线观看欧美| 国精品久久久久久国模美| 日本熟妇午夜| 日日摸夜夜添夜夜爱| 九草在线视频观看| av卡一久久| 亚洲成色77777| 六月丁香七月| 精品国产露脸久久av麻豆 | www.av在线官网国产| 日日干狠狠操夜夜爽| 中文资源天堂在线| or卡值多少钱| 免费看不卡的av| 搡老乐熟女国产| 波多野结衣巨乳人妻| 2018国产大陆天天弄谢| 午夜视频国产福利| 99九九线精品视频在线观看视频| 特级一级黄色大片| 女人十人毛片免费观看3o分钟| 黄色一级大片看看| 亚洲国产精品成人综合色| 亚洲最大成人中文| av卡一久久| 老司机影院毛片| 亚洲高清免费不卡视频| 在线天堂最新版资源| 高清欧美精品videossex| 永久网站在线| 日韩 亚洲 欧美在线| 男人舔奶头视频| 99热这里只有是精品在线观看| 国产精品福利在线免费观看| 国产精品蜜桃在线观看| 日本与韩国留学比较| 男人狂女人下面高潮的视频| 亚洲成人中文字幕在线播放| 建设人人有责人人尽责人人享有的 | 亚洲精品国产av蜜桃| 成人国产麻豆网| 高清在线视频一区二区三区| 男女下面进入的视频免费午夜| 九九爱精品视频在线观看| 亚洲伊人久久精品综合| videos熟女内射| 精品久久久久久电影网| 日本三级黄在线观看| 久久久久久伊人网av| 男女那种视频在线观看| 最近最新中文字幕免费大全7| 99热这里只有是精品在线观看| 亚洲成人精品中文字幕电影| 亚洲久久久久久中文字幕| 国产成人91sexporn| av在线老鸭窝| 亚洲av.av天堂| 老女人水多毛片| 亚洲精品乱久久久久久| a级一级毛片免费在线观看| 一级毛片电影观看| 国产探花极品一区二区| 18禁在线播放成人免费| 三级经典国产精品| 国语对白做爰xxxⅹ性视频网站| 永久网站在线| 亚洲精品成人av观看孕妇| 亚洲精品视频女| 麻豆av噜噜一区二区三区| 男的添女的下面高潮视频| 美女内射精品一级片tv| 亚洲av日韩在线播放| 久久精品夜色国产| 欧美日韩精品成人综合77777| 在线免费观看不下载黄p国产| 熟女电影av网| 亚洲熟妇中文字幕五十中出| 日韩欧美三级三区| 午夜爱爱视频在线播放| 国产精品一区二区在线观看99 | av卡一久久| 欧美+日韩+精品| 国产伦精品一区二区三区四那| 2018国产大陆天天弄谢| 精品国内亚洲2022精品成人| 777米奇影视久久| videos熟女内射| 婷婷六月久久综合丁香| 欧美成人精品欧美一级黄| 日日啪夜夜撸| 男人和女人高潮做爰伦理| 熟女电影av网| 欧美xxxx黑人xx丫x性爽| 麻豆成人av视频| 成人特级av手机在线观看| 国产淫语在线视频| 久久久成人免费电影| 中文字幕免费在线视频6| 日韩欧美精品v在线| 国产69精品久久久久777片| 亚洲激情五月婷婷啪啪| 在线免费十八禁| 身体一侧抽搐| 国产精品嫩草影院av在线观看| 亚洲三级黄色毛片| 极品教师在线视频| 国产女主播在线喷水免费视频网站 | 男女边吃奶边做爰视频| 成人综合一区亚洲| 99久久精品热视频| 欧美成人精品欧美一级黄| 纵有疾风起免费观看全集完整版 | 夜夜爽夜夜爽视频| 免费不卡的大黄色大毛片视频在线观看 | 看十八女毛片水多多多| 成人漫画全彩无遮挡| 中国美白少妇内射xxxbb| 成人特级av手机在线观看| 可以在线观看毛片的网站| 日韩三级伦理在线观看| 亚洲天堂国产精品一区在线| www.色视频.com| 色尼玛亚洲综合影院| 精品人妻视频免费看| 能在线免费观看的黄片| 中文天堂在线官网| 在线免费十八禁| 男的添女的下面高潮视频| 国产精品一区www在线观看| 纵有疾风起免费观看全集完整版 | 亚洲人成网站在线观看播放| 激情五月婷婷亚洲| 人人妻人人看人人澡| 亚洲精品国产成人久久av| 高清视频免费观看一区二区 | 日韩欧美 国产精品| 狂野欧美白嫩少妇大欣赏| 26uuu在线亚洲综合色| 天堂网av新在线| 噜噜噜噜噜久久久久久91| 国产成人freesex在线| 亚洲精品影视一区二区三区av| 国产精品综合久久久久久久免费| 18禁动态无遮挡网站| 色吧在线观看| 中国美白少妇内射xxxbb| 97热精品久久久久久| 91午夜精品亚洲一区二区三区| 婷婷色av中文字幕| 黑人高潮一二区| 亚洲无线观看免费| 亚洲婷婷狠狠爱综合网| 亚洲精品视频女| 丝瓜视频免费看黄片| 国产一区二区亚洲精品在线观看| 免费无遮挡裸体视频| 成人av在线播放网站| 美女大奶头视频| 中文在线观看免费www的网站| 亚洲最大成人av| av天堂中文字幕网| 少妇熟女欧美另类| 亚洲av男天堂| 91狼人影院| 婷婷色综合www| 国产黄频视频在线观看| 久久久a久久爽久久v久久| 欧美性感艳星| 国产av码专区亚洲av| a级毛色黄片| 免费电影在线观看免费观看| 国产午夜精品论理片| 亚洲欧美清纯卡通| av在线观看视频网站免费| 中文在线观看免费www的网站| 听说在线观看完整版免费高清| 精品国产一区二区三区久久久樱花 | 五月伊人婷婷丁香| 啦啦啦韩国在线观看视频| 一夜夜www| 麻豆久久精品国产亚洲av| 国产精品蜜桃在线观看| 日本猛色少妇xxxxx猛交久久| 99久久九九国产精品国产免费| 精品人妻熟女av久视频| 亚洲怡红院男人天堂| 免费看不卡的av| 一区二区三区乱码不卡18| 午夜福利在线在线| 精品久久久久久久久久久久久| 最近视频中文字幕2019在线8| 国产黄片美女视频| 免费在线观看成人毛片| 亚洲一级一片aⅴ在线观看| 欧美人与善性xxx| 日韩亚洲欧美综合| 在线观看一区二区三区| 欧美变态另类bdsm刘玥| 大片免费播放器 马上看| 汤姆久久久久久久影院中文字幕 | 26uuu在线亚洲综合色| 国产激情偷乱视频一区二区| 欧美xxxx黑人xx丫x性爽| 久久99精品国语久久久| 国产成人a区在线观看| 免费看美女性在线毛片视频| av在线蜜桃| 晚上一个人看的免费电影| 国产亚洲午夜精品一区二区久久 | 精品国产三级普通话版| 免费av观看视频| 禁无遮挡网站| 亚洲精品日韩av片在线观看| 精品久久久久久久久久久久久| 偷拍熟女少妇极品色| 日本熟妇午夜| 国产麻豆成人av免费视频| 少妇丰满av| 欧美高清成人免费视频www| 午夜福利在线观看吧| 街头女战士在线观看网站| 国产乱人偷精品视频| 久久久久久久午夜电影| 成年版毛片免费区| 韩国av在线不卡| 色吧在线观看| 久久久久国产网址| 好男人在线观看高清免费视频| 精品国产三级普通话版| 尾随美女入室| 97热精品久久久久久| 美女大奶头视频| 亚洲av中文字字幕乱码综合| 色播亚洲综合网| 一个人看视频在线观看www免费| 嫩草影院入口| 亚洲自拍偷在线| av专区在线播放| 毛片一级片免费看久久久久| 天天躁夜夜躁狠狠久久av| 精品亚洲乱码少妇综合久久| 99九九线精品视频在线观看视频| 三级国产精品欧美在线观看| 亚洲图色成人| 最近视频中文字幕2019在线8| 欧美区成人在线视频| 欧美 日韩 精品 国产| 91久久精品国产一区二区三区| 波野结衣二区三区在线| 男人狂女人下面高潮的视频| 97超视频在线观看视频|