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

    水平微細(xì)管內(nèi)CO2流動(dòng)沸騰壓降特性

    2017-12-22 05:36:54姜林林柳建華張良趙越
    化工學(xué)報(bào) 2017年12期
    關(guān)鍵詞:干度流率流態(tài)

    姜林林,柳建華,張良,趙越

    (1上海理工大學(xué)制冷技術(shù)研究所,上海 200093;2上海市計(jì)量測(cè)試技術(shù)研究院,上海 200093)

    水平微細(xì)管內(nèi)CO2流動(dòng)沸騰壓降特性

    姜林林1,柳建華1,張良1,趙越2

    (1上海理工大學(xué)制冷技術(shù)研究所,上海 200093;2上海市計(jì)量測(cè)試技術(shù)研究院,上海 200093)

    對(duì)CO2在內(nèi)徑1.5 mm水平微細(xì)管內(nèi)流動(dòng)沸騰換熱摩擦壓降特性進(jìn)行了實(shí)驗(yàn)研究。實(shí)驗(yàn)工況:熱通量(7.5~30 kW·m-2)、質(zhì)量流率(300~600 kg·m-2·s-1)、飽和溫度(-40~0℃)。實(shí)驗(yàn)結(jié)果表明:熱通量的增加對(duì)摩擦壓降影響很小,幾乎為零;質(zhì)量流率是影響摩擦壓降的最主要因素;隨著飽和溫度的升高摩擦壓降減小;干度對(duì)摩擦壓降影響主要由管內(nèi)流型變化導(dǎo)致。將實(shí)測(cè)摩擦壓降變化趨勢(shì)繪制于CO2流態(tài)圖中,比較發(fā)現(xiàn)理論預(yù)測(cè)摩擦壓降最大值落在環(huán)狀流末端區(qū)域。實(shí)驗(yàn)過程中對(duì)各個(gè)工況管內(nèi)流態(tài)進(jìn)行可視化研究,理論分析所采用的流態(tài)形式與實(shí)際CO2在微細(xì)通道內(nèi)所具有的流態(tài)類型基本一致。

    CO2;微細(xì)管;流動(dòng)沸騰;摩擦壓降;流態(tài)

    引 言

    微細(xì)通道換熱器具有結(jié)構(gòu)緊湊、較大單位體積接觸面積與較小制冷劑充注量等優(yōu)點(diǎn),對(duì)于目前廣泛使用的HFC、HCFC制冷劑而言,雖然隨著換熱器管徑的減小可以顯著提高其傳熱系數(shù),但隨之而來的問題是換熱過程中摩擦壓降大幅上升。如果采用CO2作為制冷劑用于微細(xì)通道換熱器內(nèi)其摩擦壓降要小得多,主要由于CO2飽和溫度所對(duì)應(yīng)的蒸發(fā)壓力要比傳統(tǒng)制冷劑高得多,且具有較高的蒸氣密度、較低的液體黏度與表面張力。Zhao等[1]研究表明飽和溫度10℃時(shí),CO2在0.86 mm水力直徑矩形通道內(nèi)流動(dòng)沸騰換熱的摩擦壓降僅為R134a摩擦壓降的60%,且在不同工況下類似;Pamitran等[2]的研究數(shù)據(jù)顯示同樣在飽和溫度 10℃且相同熱通量工況下,CO2在直徑1.5 mm管內(nèi)流動(dòng)沸騰換熱的摩擦壓降僅為R22摩擦壓降的16%;Yun等[3]測(cè)得1.53 mm矩形管內(nèi) 5℃飽和溫度時(shí),在相同熱通量與質(zhì)量流率工況下,R134a的流動(dòng)沸騰換熱摩擦壓降是CO2的 4.5倍。流動(dòng)沸騰換熱過程導(dǎo)致的摩擦壓降會(huì)引起制冷過程中產(chǎn)生溫度滑移,由于CO2蒸發(fā)壓力較高,因此摩擦壓降變化所引起蒸發(fā)壓力變化對(duì)于飽和溫度的影響較小,對(duì)于給定摩擦壓降CO2所產(chǎn)生的溫度滑移要遠(yuǎn)小于傳統(tǒng)制冷劑。

    綜合國內(nèi)外近年來關(guān)于CO2微細(xì)通道中流動(dòng)沸騰換熱過程中摩擦壓降特性實(shí)驗(yàn)研究:CO2流動(dòng)沸騰換熱摩擦壓降隨著熱通量和質(zhì)量流率的增加而增加,隨著飽和溫度的上升而減小,CO2具有較低的壓力損失且小于傳統(tǒng)制冷劑[4-13]。但是對(duì)于 CO2摩擦壓降定量研究在現(xiàn)有文獻(xiàn)中還沒有理想的理論預(yù)測(cè)模型或經(jīng)驗(yàn)關(guān)聯(lián)式。針對(duì)上述現(xiàn)象,本文對(duì)CO2在低溫工況下水平微細(xì)管內(nèi)流動(dòng)沸騰摩擦壓降特性進(jìn)行實(shí)驗(yàn)研究,結(jié)合理論預(yù)測(cè)流態(tài)圖分別分析熱通量、質(zhì)量流率、飽和溫度、干度、流態(tài)對(duì)摩擦壓降的影響,通過實(shí)際觀測(cè)的流態(tài)驗(yàn)證理論流態(tài)圖可靠性,為獲得一個(gè)可靠、精確的CO2沸騰換熱過程壓降預(yù)測(cè)方法提供理論基礎(chǔ)與數(shù)據(jù)支持。

    1 實(shí)驗(yàn)裝置與數(shù)據(jù)處理

    1.1 實(shí)驗(yàn)系統(tǒng)介紹

    本實(shí)驗(yàn)系統(tǒng)由3部分組成:CO2循環(huán)系統(tǒng)、脈管制冷機(jī)系統(tǒng)和數(shù)據(jù)采集系統(tǒng)。實(shí)驗(yàn)系統(tǒng)原理圖如圖1所示,實(shí)驗(yàn)系統(tǒng)主要設(shè)備配置如表1所示。

    表1 實(shí)驗(yàn)系統(tǒng)主要設(shè)備配置Table 1 Main equipment of experimental system

    圖1 實(shí)驗(yàn)系統(tǒng)原理圖Fig.1 Schematic of test system

    1.2 可視化與干度測(cè)量

    為了驗(yàn)證理論預(yù)測(cè)流態(tài)圖的可靠性,本實(shí)驗(yàn)系統(tǒng)在測(cè)試段進(jìn)出口位置設(shè)計(jì)了用于觀測(cè)測(cè)試段流態(tài)的可視段。在測(cè)試段流態(tài)觀測(cè)時(shí),為了避免外界環(huán)境與可視段熱量傳遞對(duì)可視段內(nèi)流態(tài)產(chǎn)生影響,實(shí)驗(yàn)系統(tǒng)將可視段置于封閉空間內(nèi),實(shí)驗(yàn)過程中采用真空泵抽空該封閉空間的空氣,維持有機(jī)玻璃管外部處于真空環(huán)境,避免可視段與外界環(huán)境的熱交換。

    測(cè)試段總長(zhǎng)300 mm,外徑6 mm,內(nèi)徑1.5 mm,加熱段 150 mm。實(shí)驗(yàn)系統(tǒng)在測(cè)試段處采用加熱絲均勻纏繞于測(cè)試管上,對(duì)測(cè)試段所施加電功率由功率計(jì)直接測(cè)得,在加熱段前后布置引壓管測(cè)試總壓差,測(cè)試段布置如圖2所示。

    圖2 實(shí)驗(yàn)段布置Fig.2 Experiment section layout

    1.3 數(shù)據(jù)處理及不確定度

    制冷劑進(jìn)入預(yù)熱器時(shí)為過冷狀態(tài),制冷劑從過冷直至飽和時(shí)溫度由能量平衡方程得出,飽和后制冷劑溫度由飽和壓力換算而得[式(1)]。加熱段焓值的變化由加熱段中間位置計(jì)算獲得,局部干度計(jì)算見式(2)、式(3)。

    根據(jù)相關(guān)文獻(xiàn)研究[14]選擇 Thome等[15]提出的空隙率關(guān)聯(lián)式,該關(guān)聯(lián)式在Rouhani等[16]提出的漂移流模型基礎(chǔ)上進(jìn)行了更新,在較大實(shí)驗(yàn)條件跨度下仍具有較高的精確度和適用性,且參數(shù)少、易于計(jì)算,如式(5)所示。

    將測(cè)試段進(jìn)出口干度值代入,即可得到加速壓降,最終的摩擦壓降計(jì)算見式(6)。

    數(shù)據(jù)采集系統(tǒng)對(duì)所有測(cè)量數(shù)據(jù)均采用多次測(cè)量取平均值的方式獲得,為了對(duì)實(shí)驗(yàn)結(jié)果分析獲得較高概率的置信區(qū)間,依據(jù)技術(shù)規(guī)范《JJF1059—1999 測(cè)量不確定度評(píng)定與表示》對(duì)測(cè)試數(shù)據(jù)進(jìn)行不確度A類評(píng)定,合成不確定度評(píng)定以及擴(kuò)展不確定度評(píng)定。依據(jù)不確定度計(jì)算公式可得摩擦壓降不確定度為 6.7%~11.3%,干度不確定度為 5.2%~9.5%,含氣率不確定度為8.6%~14.3%。

    2 實(shí)驗(yàn)結(jié)果與分析

    2.1 熱通量對(duì)壓降的影響

    如圖3所示,在熱通量7.5~30 kW·m-2范圍內(nèi)進(jìn)行了摩擦壓降實(shí)驗(yàn)對(duì)比研究。實(shí)驗(yàn)結(jié)果表明:不同熱通量下摩擦壓降對(duì)干度的變化趨勢(shì)是一致的,在低干度區(qū)域摩擦壓降隨著干度的增加快速增加,在高干度區(qū)域摩擦壓降的增速放緩,在高干度區(qū)域出現(xiàn)極值,隨后摩擦壓降回落,在整個(gè)干度范圍內(nèi)熱通量大小對(duì)摩擦壓降影響較弱,幾乎為零。對(duì)此解釋為:無論熱通量如何變化,在相同質(zhì)量流率相同干度相同飽和溫度下,氣液兩相相對(duì)速度沒有變化,CO2物性也沒有變化,故熱通量對(duì)摩擦壓降幾乎沒有影響。通過理論分析隨著熱通量的增加,管內(nèi)液體蒸發(fā)越發(fā)劇烈,流態(tài)也提前進(jìn)入高速的環(huán)狀流,故達(dá)到摩擦壓降極值所對(duì)應(yīng)的時(shí)間會(huì)變小,這是熱通量對(duì)摩擦壓降最大影響之處。

    圖3 飽和溫度-10℃、質(zhì)量流率 400 kg·m-2·s-1、不同熱通量時(shí)管內(nèi)CO2實(shí)驗(yàn)?zāi)Σ翂航岛蛡鳠嵯禂?shù)Fig.3 Experimental frictional pressure drop and heat transfer coefficient for CO2 at Tsat=-10℃ with initial mass flow rate G=400 kg·m-2·s-1 at different heat flux

    與CO2流動(dòng)沸騰換熱特性研究實(shí)驗(yàn)結(jié)果對(duì)比來看,摩擦壓降極值點(diǎn)一般都出現(xiàn)在環(huán)狀流區(qū)域,而傳熱系數(shù)最大值大都出現(xiàn)在環(huán)狀流向霧狀流轉(zhuǎn)變的區(qū)域[17-18],如圖3所示,在同一工況下壓降極值點(diǎn)所對(duì)應(yīng)的干度比傳熱系數(shù)極值點(diǎn)所對(duì)應(yīng)的干度要小。對(duì)此解釋為:隨著干度的增加,工質(zhì)流型在逐漸發(fā)生變化,在中高干度區(qū)域大部分工質(zhì)處于環(huán)狀流流態(tài),環(huán)狀流的形成增強(qiáng)了環(huán)形液膜與管壁面以及氣液兩相分界面的擾動(dòng),摩擦壓降值達(dá)到最大,隨著干度繼續(xù)增大,環(huán)形液膜逐漸變薄直至蒸干,在此過程中管內(nèi)工質(zhì)開始部分干涸,此時(shí)液膜與壁面的應(yīng)力以及氣液兩相分界面的擾動(dòng)雙雙下降,因此摩擦壓降值逐漸開始下降;而對(duì)于傳熱系數(shù)而言,隨著工質(zhì)干度的增大甚至是部分干涸,中心高速氣流使得液滴夾帶現(xiàn)象增加,部分干涸區(qū)域仍然有較高的概率重新濕潤(rùn),因此傳熱系數(shù)還有繼續(xù)上升的空間。理論分析和實(shí)驗(yàn)結(jié)果表明在同樣工況下壓降極值點(diǎn)的出現(xiàn)優(yōu)先于傳熱系數(shù)極值點(diǎn)。

    2.2 質(zhì)量流率對(duì)壓降的影響

    如圖4所示,對(duì)質(zhì)量流率300~600 kg·m-2·s-1范圍內(nèi)進(jìn)行了摩擦壓降實(shí)驗(yàn)對(duì)比研究。實(shí)驗(yàn)結(jié)果顯示:微細(xì)通道內(nèi)摩擦壓降隨質(zhì)量流率的增加而增大,且增幅較為明顯;隨著干度逐漸增加,質(zhì)量流率對(duì)摩擦壓降提升效果更加明顯,這與高干度區(qū)流型轉(zhuǎn)變?yōu)榄h(huán)狀流有很大的關(guān)系;當(dāng)質(zhì)量流率較小時(shí),摩擦壓降隨著干度的增加略微增加,質(zhì)量流率越大,摩擦壓降值隨干度變化速率越快。實(shí)驗(yàn)結(jié)果表明:隨著質(zhì)量流率的增大,摩擦壓降最大值有左移的趨勢(shì)。對(duì)此解釋為:隨著換熱的進(jìn)行中心氣流速度加快,質(zhì)量流率增大使得流態(tài)提前向高速環(huán)狀流轉(zhuǎn)變,摩擦壓降極大值提前出現(xiàn);與此同時(shí)由于過高的流速造成液滴夾帶現(xiàn)象增加,對(duì)管壁面液膜造成了破壞,受液滴沖擊形成的干涸現(xiàn)象提前發(fā)生,摩擦壓降值下降。通過數(shù)據(jù)分析可得:在相同干度下,隨著質(zhì)量流率的增大,摩擦壓降明顯增大,其原因在于質(zhì)量流率的增加使管內(nèi)兩相平均流速的增加,摩擦壓降數(shù)值顯著增大;而同樣質(zhì)量流率情況下,隨著干度的增大,壓降的變化主要受兩相流流型轉(zhuǎn)變所帶來的影響。

    圖4 飽和溫度-10℃、熱通量30 kW·m-2、不同質(zhì)量流率時(shí)管內(nèi)CO2實(shí)驗(yàn)?zāi)Σ翂航礔ig.4 Experimental frictional pressure drop for CO2 at Tsat=-10℃ with initial heat flux q=30 kW·m-2 at different mass flow rate

    對(duì)于兩相流而言,氣相和液相各自的速度是影響摩擦壓降的另一個(gè)重要因素。氣相與液相速度分別由式(7)、式(8)計(jì)算獲得。

    則兩相相對(duì)速度由式(9)計(jì)算獲得

    如圖5所示,CO2兩相相對(duì)速度隨著質(zhì)量流率增大而增大,且增幅隨著干度增大而增大。兩相相對(duì)速度與摩擦壓降梯度的關(guān)系與質(zhì)量流率的趨勢(shì)是一致的,因此可以判斷質(zhì)量流率和干度對(duì)摩擦壓降的影響還體現(xiàn)在氣液兩相相對(duì)速度的改變,相對(duì)速度的改變使得兩相流壓降中重要的組成部分-兩相流氣液分界面摩擦壓降改變,最終導(dǎo)致摩擦壓降值的變化。

    圖5 飽和溫度-10℃、熱通量30 kW·m-2、不同質(zhì)量流率時(shí)管內(nèi)CO2實(shí)驗(yàn)兩相相對(duì)速度Fig.5 Experimental gas-liquid phase relative speed for CO2 at Tsat=-10℃ with initial heat flux q=30 kW·m-2 at different mass flow rate

    結(jié)合CO2流動(dòng)沸騰換熱特性研究實(shí)驗(yàn)結(jié)果,質(zhì)量流率對(duì)傳熱系數(shù)的影響不大,但質(zhì)量流率的增加導(dǎo)致干涸出現(xiàn)干度有所降低,減小高效換熱區(qū)域,影響整體平均傳熱系數(shù)[19-20]?;谏鲜隼碚撗芯酷槍?duì)CO2微細(xì)通道換熱器設(shè)計(jì)適用于低質(zhì)量流率運(yùn)行工況且能維持較高的傳熱系數(shù),降低壓降損失,有利于換熱器均勻補(bǔ)液。

    2.3 飽和溫度對(duì)壓降的影響

    如圖6所示,對(duì)飽和溫度-40~0℃范圍內(nèi)進(jìn)行了摩擦壓降實(shí)驗(yàn)對(duì)比研究。實(shí)驗(yàn)結(jié)果表明:飽和溫度對(duì)CO2摩擦壓降特性具有顯著的影響,隨著飽和溫度的升高摩擦壓降變小,且較高飽和溫度下摩擦壓降隨干度上升的斜率明顯低于低飽和溫度下上升的斜率。對(duì)此解釋為:飽和溫度的變化改變了CO2熱物性,隨著飽和溫度的上升,一方面氣相密度升高,使得兩相流內(nèi)氣液相流速和平均流速雙雙下降,最終導(dǎo)致摩擦壓降下降,另一方面飽和壓力的上升使得液體黏度降低、在中干度區(qū)域強(qiáng)烈的核態(tài)沸騰使管內(nèi)壁面液膜變薄更易破裂,流體表面張力作用減弱,使得流體與壁面之間的摩擦壓降減小,故總摩擦壓降減小。

    圖6 熱通量 30 kW·m-2、質(zhì)量流率 300 kg·m-2·s-1、不同飽和溫度時(shí)管內(nèi)CO2實(shí)驗(yàn)?zāi)Σ翂航堤荻菷ig.6 Experimental frictional pressure drop for CO2 at G=300 kg·m-2·s-1 with initial heat flux q=30 kW·m-2 at different saturation temperature

    2.4 干度對(duì)壓降的影響

    實(shí)驗(yàn)結(jié)果表明:CO2水平微細(xì)通道內(nèi)流通沸騰摩擦壓降隨干度變化趨勢(shì)基本一致,在中低干度區(qū)域摩擦壓降梯度隨著干度增大而增大,在高干度區(qū)域增速放緩,摩擦壓降在某個(gè)干度存在極值,越過極值后摩擦壓降隨干度增大而緩慢減?。荒Σ翂航翟谫|(zhì)量流率較高的條件下增速較快,在對(duì)于低質(zhì)量流率工況下,當(dāng)干度增大時(shí),摩擦壓降的增大速度趨于平緩;在飽和壓力較低的條件下摩擦壓降隨著干度呈正相關(guān)增長(zhǎng),當(dāng)達(dá)到極大值后摩擦壓降出現(xiàn)明顯下降,高壓力條件下摩擦壓降梯度增長(zhǎng)平緩。究其原因主要是管內(nèi)流型的改變引起的,在核態(tài)沸騰開始階段,管內(nèi)流型大多為塞狀流,在管內(nèi)壁面不斷生成的氣泡擾動(dòng)下,兩相流動(dòng)摩擦壓降顯著增加;此外隨著內(nèi)部液態(tài)工質(zhì)不斷蒸發(fā),生成的氣泡聚合作用使得對(duì)液膜的擠壓作用力增大,易形成較為穩(wěn)定的環(huán)狀流,促使摩擦壓降增加,直至部分干涸現(xiàn)象開始出現(xiàn),引起摩擦壓降增大的擾動(dòng)作用達(dá)到最大,此后管內(nèi)壁面生成的氣泡擾動(dòng)作用有限,壓降趨于平穩(wěn)。當(dāng)流型由環(huán)狀流向霧狀流轉(zhuǎn)變時(shí),氣相變成連續(xù)相,液態(tài)CO2以滴狀的形式彌散在中心高速氣流中,此時(shí)流體與壁面的摩擦阻力相對(duì)于環(huán)狀流時(shí)會(huì)有所下降。

    3 流態(tài)圖對(duì)比驗(yàn)證

    實(shí)驗(yàn)結(jié)果表明摩擦壓降值的變化與管內(nèi)流態(tài)具有密不可分的聯(lián)系。圖7所示為實(shí)驗(yàn)中觀察到穩(wěn)態(tài)條件下氣液兩相流流型匯總,分別為泡狀流(bubbly flow)、塞狀流(slug flow)、間歇流(intermittent flow)、層流(stratified-wavy flow)、波狀流(wavy flow)、環(huán)狀流(annular flow)、霧狀流(mist flow)共7種流型。通過對(duì)同一工況不同干度下流態(tài)的分析,推演出流態(tài)發(fā)展過程:流體受熱在管壁面生成氣泡,隨著氣泡逐漸生長(zhǎng)脫離壁面形成一個(gè)個(gè)直徑小于管截面的圓形或細(xì)長(zhǎng)形氣泡懸浮于管內(nèi),形成泡狀流,氣泡的形狀隨質(zhì)量流率與熱通量而變化。隨著換熱持續(xù)進(jìn)行,氣泡生成速度加快,大量的氣泡聚并在一起氣泡直徑充滿管內(nèi)徑,形成間歇流或塞狀流,當(dāng)換熱使大量液體氣化導(dǎo)致管內(nèi)流速迅速增加時(shí)管內(nèi)流態(tài)演變?yōu)橹虚g具有快速流動(dòng)的蒸汽核心,管壁附著一層薄液膜,液膜與蒸汽界面表現(xiàn)為波狀,形成波狀流。在這其中分層流出現(xiàn)過,由于換熱過程極快,所以出現(xiàn)時(shí)間極為短暫。隨著干度繼續(xù)增大管內(nèi)中心氣流速度繼續(xù)加快,液體被擠壓到管內(nèi)壁兩側(cè)形成環(huán)狀流。環(huán)狀流過程中蒸汽流流動(dòng)時(shí)夾帶著不規(guī)則的液滴,同時(shí)薄液膜與管壁間的核態(tài)沸騰受到抑制,當(dāng)管壁液膜隨換熱過程由薄變干直至消失后整個(gè)換熱流態(tài)進(jìn)入霧狀流。實(shí)驗(yàn)結(jié)果表明理論分析所采用的流態(tài)形式與實(shí)際 CO2在微細(xì)通道內(nèi)所具有的流態(tài)類型基本一致。

    圖7 管內(nèi)CO2流動(dòng)沸騰換熱流態(tài)Fig.7 Flow pattern of CO2 during flow boiling heat transfer

    研究流態(tài)圖的目的就是揭示熱通量、質(zhì)量流率、飽和溫度等參數(shù)對(duì)流態(tài)的影響,用于解釋換熱特性和壓降特性變化趨勢(shì)。Cheng等[21-22]通過對(duì)流態(tài)的研究獲得了CO2流態(tài)預(yù)測(cè)模型,該流態(tài)圖反映了隨流體干度和質(zhì)量流率變化所對(duì)應(yīng)的流態(tài)轉(zhuǎn)變曲線,具有重要的實(shí)際意義,通過確定CO2流動(dòng)沸騰換熱過程中不同質(zhì)量流率所對(duì)應(yīng)的轉(zhuǎn)變干度即可獲得其流態(tài)圖。圖8~圖11為不同換熱工況下實(shí)測(cè)摩擦壓降與對(duì)應(yīng)理論預(yù)測(cè)流態(tài),圖中流態(tài)轉(zhuǎn)變曲線為泡狀流-間歇流轉(zhuǎn)變曲線(G bubble)、波狀流-環(huán)狀流轉(zhuǎn)變曲線(G wavy)、層流-波狀流轉(zhuǎn)變曲線(G strat)、間歇流-環(huán)狀流轉(zhuǎn)變曲線(X IA)、環(huán)狀流-干涸轉(zhuǎn)變曲線(G dryout)、干涸-霧狀流轉(zhuǎn)變曲線(G mist)。

    圖8 飽和溫度-10℃、質(zhì)量流率400 kg·m-2·s-1、熱通量7.5 kW·m-2時(shí)CO2實(shí)測(cè)摩擦壓降梯度與對(duì)應(yīng)流態(tài)Fig.8 Experimental frictional pressure drop data for CO2 at T=-10℃,G=400 kg·m-2·s-1,q=7.5 kW·m-2 and corresponding flow pattern map

    從理論流型圖與實(shí)驗(yàn)值對(duì)比可以看出,熱通量對(duì)干涸區(qū)和霧狀區(qū)有決定性作用,熱通量越大,代表環(huán)狀流區(qū)域的面積越小,包括橫掠干度和縱掠質(zhì)量流率都對(duì)應(yīng)減小。通常來說隨著熱通量逐漸增大,干涸區(qū)的起始干度減小,但干涸區(qū)和霧狀區(qū)的干度跨度增大,這嚴(yán)重影響CO2流動(dòng)沸騰換熱特性。雖然熱通量對(duì)高干度區(qū)域的流型轉(zhuǎn)變影響重大,但是熱通量對(duì)摩擦壓降影響較小,摩擦壓降值不會(huì)發(fā)生突變,且間歇流、環(huán)狀流對(duì)摩擦壓降變化趨勢(shì)起最主要影響。實(shí)驗(yàn)中不同質(zhì)量流率對(duì)流型影響非常大,直接影響管內(nèi)摩擦壓降。隨著質(zhì)量流率增大干涸干度明顯減小,不同質(zhì)量流率直接決定了其換熱過程中所經(jīng)歷的流態(tài)。在低干度區(qū)域大多出現(xiàn)間歇流,此時(shí)摩擦壓降絕對(duì)值和壓降波動(dòng)都相對(duì)較小,摩擦壓降小幅度上升;隨著干度增大流型向環(huán)狀流轉(zhuǎn)變,過度階段壓降發(fā)生一定的波動(dòng),出現(xiàn)時(shí)大時(shí)小的現(xiàn)象,絕大多數(shù)壓降最大值都落在環(huán)狀流末端區(qū)域。對(duì)此解釋為:環(huán)狀流流型時(shí)氣液兩相相對(duì)流速差值非常大,氣液兩相之間的摩擦壓降值很大,工質(zhì)與壁面之間摩擦壓降值也一直上升,這樣無論是工質(zhì)與壁面之間的作用力還是工質(zhì)氣液兩相之間作用力都達(dá)到最大;在高干度區(qū)域管內(nèi)壁面附著的液膜越來越薄,中心氣流速度加快,液滴夾帶現(xiàn)象增強(qiáng),管道上部開始出現(xiàn)部分干涸現(xiàn)象,此時(shí)摩擦壓降雖無大幅度波動(dòng),但絕對(duì)值已開始下降;隨著干度繼續(xù)上升,摩擦壓降一直減小。研究發(fā)現(xiàn)當(dāng)流型從環(huán)狀流向霧狀流轉(zhuǎn)變過程中,并沒有明顯的壓降突變。這說明即使環(huán)狀流液膜變薄并發(fā)生干涸,壓降也不會(huì)有大的波動(dòng)。飽和溫度對(duì)流型轉(zhuǎn)變有重要影響,隨著飽和溫度的降低,代表干涸區(qū)域的面積越小,包括橫掠干度和縱掠質(zhì)量流率都相應(yīng)減小。這主要是因?yàn)殡S著飽和溫度的降低,液氣密度比、表面張力、液體黏度、液體熱導(dǎo)率增大,使得管內(nèi)氣流速度加快,液相流體能夠更好地附著在管內(nèi)壁面形成環(huán)狀流,在高熱通量該質(zhì)量流率工況下,管內(nèi)液滴夾帶現(xiàn)象比較明顯,管內(nèi)流態(tài)直接由環(huán)狀流轉(zhuǎn)向霧狀流。

    圖9 飽和溫度-10℃、質(zhì)量流率400 kg·m-2·s-1、熱通量30 kW·m-2時(shí)CO2實(shí)測(cè)摩擦壓降梯度與對(duì)應(yīng)流態(tài)Fig.9 Experimental frictional pressure drop data for CO2 at T=-10℃,G=400 kg·m-2·s-1,q=30 kW·m-2 and corresponding flow pattern map

    圖10 飽和溫度-10℃、質(zhì)量流率600 kg·m-2·s-1、熱通量30 kW·m-2時(shí)CO2實(shí)測(cè)摩擦壓降梯度與對(duì)應(yīng)流態(tài)Fig.10 Experimental frictional pressure drop data for CO2 at T=-10℃,G=600 kg·m-2·s-1,q=30 kW·m-2 and corresponding flow pattern map

    圖11 飽和溫度-35℃、質(zhì)量流率300 kg·m-2·s-1、熱通量30 kW·m-2時(shí)CO2實(shí)測(cè)摩擦壓降梯度與對(duì)應(yīng)流態(tài)Fig.11 Experimental frictional pressure drop data for CO2 at T=-35℃,G=300 kg·m-2·s-1,q=30 kW·m-2 and corresponding flow pattern map

    通常情況下摩擦壓降隨著氣體和液體流速增加而增加。在相同的氣速下,表觀液速越大,流體與管壁的相對(duì)速度也越大,摩擦壓力降也就越大;而液相流速相當(dāng)時(shí)氣速越高,液層越容易被排擠到兩側(cè),與管壁接觸更加密切,液層與管壁的作用力也更強(qiáng)烈,從而摩擦壓力降隨氣速增大而增大。實(shí)驗(yàn)結(jié)果顯示在間歇流、環(huán)狀流、干涸區(qū)、霧狀流流態(tài)轉(zhuǎn)變時(shí)摩擦壓降變化較為平緩,沒有發(fā)生突變,尤其是在干涸區(qū)和霧狀流區(qū)域,雖然干涸區(qū)的摩擦壓降值開始降低,但變化平緩,并沒有像傳熱系數(shù)一樣有突降[23-24]。

    4 結(jié) 論

    本文針對(duì)CO2在水平微細(xì)通道管內(nèi)流動(dòng)沸騰換熱過程進(jìn)行了實(shí)驗(yàn)研究,分別分析熱通量、質(zhì)量流率、飽和溫度、干度對(duì)壓降特性和流態(tài)轉(zhuǎn)變的影響,并將實(shí)驗(yàn)結(jié)果與理論流態(tài)預(yù)測(cè)模型進(jìn)行對(duì)比分析,得到如下結(jié)論。

    (1)熱通量對(duì)摩擦壓降影響很小,幾乎為零,隨著熱通量的增加,達(dá)到摩擦壓降最大值所對(duì)應(yīng)的時(shí)間會(huì)變短;質(zhì)量流率是影響摩擦壓降最主要的因素,摩擦壓降隨質(zhì)量流率增加大幅增加,與此同時(shí)摩擦壓降最大值所對(duì)應(yīng)的干度值減??;CO2不同飽和溫度物性的改變是造成其不同溫度時(shí)壓降特性差異的主要原因,隨著飽和溫度的升高,摩擦壓降變小;摩擦壓降在低干度區(qū)域隨干度增速較快,在高干度區(qū)域增速放緩,在某個(gè)干度值存在極值,越過極值后摩擦壓降緩慢下降。

    (2)通過理論計(jì)算CO2微細(xì)通道內(nèi)流動(dòng)沸騰流態(tài)圖,分析得熱通量對(duì)高干度區(qū)域流態(tài)轉(zhuǎn)變有顯著影響,熱通量越高,代表環(huán)狀流與干涸區(qū)面積越??;質(zhì)量流率大小直接決定了換熱過程所經(jīng)歷的流態(tài),飽和溫度對(duì)流型轉(zhuǎn)變有重要影響,主要是由于不同飽和溫度對(duì)應(yīng)的熱物性不同,改變了流型轉(zhuǎn)變特性。

    (3)實(shí)驗(yàn)拍攝到不同工況不同干度下CO2微細(xì)通道內(nèi)流動(dòng)沸騰換熱過程中流態(tài),并與理論計(jì)算流態(tài)圖進(jìn)行對(duì)比分析,實(shí)驗(yàn)結(jié)果顯示觀測(cè)到的流態(tài)與理論預(yù)測(cè)流態(tài)較為吻合,理論模型預(yù)測(cè)摩擦壓降最大值大都落在環(huán)狀流末端區(qū)域。

    符 號(hào) 說 明

    cp——比定壓熱容,J·kg-1·K-1

    D——水力直徑,m

    G——質(zhì)量流率,kg·m-2·s-1

    h——焓,kJ·kg-1

    k——熱導(dǎo)率,W·m-2·K-1

    L——加熱段長(zhǎng)度,m

    M——摩爾質(zhì)量,kg·kmol-1

    m——質(zhì)量流量,kg·s-1

    q——熱通量,W·m-2

    x——干度

    α——空隙率

    δ——液膜厚度,m

    ρ——密度,kg·m-3

    下角標(biāo)

    in——進(jìn)口

    l——液體

    out——出口

    pre——預(yù)熱

    r——制冷劑

    sul——過冷

    tp——兩相

    v——?dú)怏w

    [1]ZHAO Y,MOLKI M,OHADI M,et al.Flow boiling of CO2in micro-channels[J].ASHRAE Trans.,2000,106:437-445.

    [2]PAMITRAN A S,CHOI K I,OH J T,et al.Two-phase pressure drop during CO2vaporization in horizontal smooth mini-channels[J].International Journal of Refrigeration,2008,31(8):1375-1383.

    [3]YUN R,KIM Y.Two-phase pressure drop of CO2in mini tubes and micro-channels[J].Microscale Thermophysical Engineering,2004,8(3):259-270.

    [4]PETTERSEN J.Flow vaporization of CO2in microchannel tubes[J].Experimental Thermal and Fluid Science,2004,28(2):111-121.

    [5]PETTERSEN J.Two-phase flow pattern,heat transfer,and pressure drop in microchannel vaporization of CO2[J].ASHRAE Trans.,2003,109:523-532.

    [6]JEONG S,CHO E,KIM H K.Evaporative heat transfer and pressure drop of CO2in a microchannel tube[C]//Proceedings of the 3rd International Conference on Micro-channels and Mini-channels.Toronto,2005:103-108.

    [7]YAMAMOTO T,UEDA Y,ISHIHARA I,et al.Flow boiling heat transfer of carbon dioxide at high pressure in horizontal mini-channels[C]//Proceedings of the 6th International Conference on Multiphase Flow.Leipzig,2007.

    [8]JEONG S,PARK D.Evaporating heat transfer and pressure drop of CO2in a multichannel micro-tube[C]//Proceedings of the International Congress of Refrigeration.2007:1-7.

    [9]CHENG L X,THOME J R.Cooling of microprocessors using flow boiling of CO2in a micro-evaporator:preliminary analysis and performance comparison[J].Applied Thermal Engineering,2009,29(11):2426-2432.

    [10]OH H K,SON C H.Flow boiling heat transfer and pressure drop characteristics of CO2in horizontal tube of 4.57 mm inner diameter[J].Applied Thermal Engineering,2011,31(2):163-172.

    [11]WU J,KOETTIG T,FRANKE C,et al.Investigation of heat transfer and pressure drop of CO2two-phase flow in a horizontal mini-channel[J].International Journal of Heat and Mass Transfer,2011,54(9):2154-2162.

    [12]DUCOULOMBIE M,COLASSON S,BONJOUR J,et al.Carbon dioxide flow boiling in a single microchannel(Ⅰ):Pressure drops[J].Experimental Thermal and Fluid Science,2011,35(4):81-596.

    [13]吳曉敏,趙紅藝,王維城,等.CO2在細(xì)徑管內(nèi)蒸發(fā)換熱的實(shí)驗(yàn)研究[J].工程熱物理學(xué)報(bào),2005,26(5):823-825.WU X M,ZHAO H Y,WANG W C,et al.Experimental study on evaporating heat transfer of CO2in thin tube[J].Journal of Engineering Thermophysics,2005,26(5):823-825.

    [14]FANG X D.Review of correlations of flow boiling heat transfer coefficients for carbon dioxide[J].International Journal of Refrigeration,2013,36(8):2017-2039.

    [15]THOME J R,EI H.Two-phase flow pattern map for evaporation in horizontal tubes[C]// Proceeding of the First International Conference on Heat Transfer,Fluid Mechanics and Thermodynamics.Kruger Park,2002:182-188.

    [16]ROUHANI S Z,AXELSSON E.Calculation of void volume fraction in the subcooled and quality boiling regions[J].International Journal of Heat and Mass Transfer,1970,13(2):383-393.

    [17]JIANG L L,LIU J H,ZHANG L,et al.Characteristics of heat transfer for CO2flow boiling at low temperature in mini-channel[J].International Journal of Heat and Mass Transfer,2017,108:2120-2129.

    [18]陸至羚,柳建華,張良,等.微細(xì)通道內(nèi)CO2沸騰換熱與干涸特性[J].化工進(jìn)展,2015,34(8):2961-2966.LU Z L,LIU J H,ZHANG L,et al.Heat transfer and dry-out characteristics of CO2in mini-channel [J].Chemical Industry and Engineering Progress,2015,34(8):2961-2966.

    [19]姜林林,柳建華,張良,等.微細(xì)管內(nèi)CO2流動(dòng)沸騰換熱特性研究[J].制冷學(xué)報(bào),2014,35(6):58-63.JIANG L L,LIU J H,ZHANG L,et al.Characteristics of heat transfer for CO2flow boiling in mini-channel[J].Journal of Refrigeration,2014,35(6):58-63.

    [20]吳昊,柳建華,張良,等.CO2微細(xì)通道流動(dòng)沸騰換熱干涸特性[J].化工學(xué)報(bào),2015,66(5):1676-1682.WU H,LIU J H,ZHANG L,et al.Characteristics of dryout for CO2flow boiling heat transfer in mini-channels[J].CIESC Journal,2015,66(5):1676-1682.

    [21]CHENG L X,RIBATSKI G,QUIBEN J M.New prediction methods for CO2evaporation inside tubes(Ⅰ):A two-phase flow pattern map and a flow pattern based phenomenological model for two-phase flow frictional pressure drops[J].International Journal of Heat and Mass Transfer,2008,51(1):111-124.

    [22]CHENG L X,RIBATSKI G,THOME J R.New prediction methods for CO2evaporation inside tubes(Ⅱ):An updated general flow boiling heat transfer model based on flow patterns[J].International Journal of Heat and Mass Transfer,2008,51(1):125-135.

    [23]OH H K,SON C H.Flowing boiling heat transfer and pressure drop characteristics of CO2in horizontal tube of 4.57 mm inner diameter[J].Applied Thermal Engineering,2011,31(2):163-172.

    [24]ALBERTO C,DAVIDE D C,LUISA R.Heat transfer and pressure drop of natural refrigerants in mini-channels (low charge equipment)[J].International Journal of Refrigeration,2013,36(2):287-300.

    date:2017-05-15.

    Prof.LIU Jianhua,lwnlwn_liu@163.com

    supported by the Capacity Building Plan for Some Non-military Universities and Colleges of Shanghai Scientific Committee(16060502600).

    Flow boiling pressure drop characteristics of CO2in horizontal micro tube

    JIANG Linlin1,LIU Jianhua1,ZHANG Liang1,ZHAO Yue2
    (1Institute of Refrigeration Technology,University of Shanghai for Science and Technology,Shanghai200093,China;2Shanghai Institute of Metrology and Measurement Technology,Shanghai200093,China)

    The study on flow boiling frictional pressure drop characteristics of CO2in horizontal micro tube which internal diameter is 1.5mm was made.Experimental conditions:heat flux(7.5—30 kW·m-2),mass flow rate(300—600 kg·m-2·s-1),and saturation temperature(-40~0℃).Experimental results show:The increase of heat flux makes little effect on frictional pressure drop,almost zero; mass flow rate is the main factor that affects frictional pressure drop; the frictional pressure drop decreases as the saturation temperature increases; the effect of vapor quality on the frictional pressure drop is mainly caused by the change of the flow pattern in the tube.The trend of the experimental measured frictional pressure drop is plotted in CO2flow regime chart and it is found that the theoretical prediction of the maximum value of the frictional pressure drop falls in the end of the annular flow region.Making visible study on flow regime change in each working condition during the experiment process and the theoretical pattern is consistent with the flow pattern of the actual CO2in the micro tube in general.

    carbon dioxide; micro tube; flow boiling; frictional pressure drop; flow regime

    TK 121; TB 61+2

    A

    0438—1157(2017)12—4576—09

    10.11949/j.issn.0438-1157.20170613

    2017-05-15收到初稿,2017-09-19收到修改稿。

    聯(lián)系人:柳建華。

    姜林林(1990—),男,博士研究生。

    上海市部分地方院校能力建設(shè)專項(xiàng)計(jì)劃(16060502600)。

    猜你喜歡
    干度流率流態(tài)
    側(cè)邊機(jī)組故障對(duì)泵站前池流態(tài)的影響
    CDN 緩存空間調(diào)整優(yōu)化方法
    帶交換和跳躍的一維雙向自驅(qū)動(dòng)系統(tǒng)的仿真研究
    多束光在皮膚組織中傳輸?shù)拿商乜_模擬
    改進(jìn)邊界條件的非恒定流模型在城市河流橡膠壩流態(tài)模擬中的應(yīng)用
    聚合物氣輔共擠中熔體體積流率對(duì)擠出脹大的影響
    注汽鍋爐應(yīng)用干度自控系統(tǒng)提高稠油采收率
    動(dòng)態(tài)流態(tài)冰蓄冷系統(tǒng)在千級(jí)凈化廠房的應(yīng)用
    基于TM遙感影像的河口流態(tài)信息半定量化研究
    稠油注汽鍋爐蒸汽干度測(cè)控裝置的應(yīng)用
    成人三级做爰电影| 亚洲国产高清在线一区二区三 | 色在线成人网| 久久久精品欧美日韩精品| 免费看十八禁软件| 国产片内射在线| 国产免费男女视频| 18禁黄网站禁片免费观看直播| 90打野战视频偷拍视频| xxxwww97欧美| 中文资源天堂在线| 好男人电影高清在线观看| 亚洲熟妇中文字幕五十中出| 男人的好看免费观看在线视频 | 中文资源天堂在线| 亚洲国产精品sss在线观看| 国产三级黄色录像| 亚洲人成电影免费在线| 国产精品电影一区二区三区| 日韩精品青青久久久久久| 亚洲熟妇熟女久久| or卡值多少钱| 首页视频小说图片口味搜索| 亚洲一区二区三区不卡视频| 波多野结衣高清无吗| АⅤ资源中文在线天堂| 丝袜人妻中文字幕| 我的亚洲天堂| 一区二区三区高清视频在线| 色综合站精品国产| 中文字幕精品亚洲无线码一区 | 在线永久观看黄色视频| 久久久久九九精品影院| 夜夜躁狠狠躁天天躁| 日本黄色视频三级网站网址| www.999成人在线观看| 久久久久九九精品影院| 亚洲第一电影网av| 精品国产超薄肉色丝袜足j| 黄色女人牲交| 一区二区三区精品91| 级片在线观看| 久久婷婷成人综合色麻豆| 亚洲男人天堂网一区| 高清毛片免费观看视频网站| 啦啦啦韩国在线观看视频| 亚洲精品国产区一区二| 色综合亚洲欧美另类图片| 免费搜索国产男女视频| 91大片在线观看| 亚洲五月色婷婷综合| 久热这里只有精品99| 人成视频在线观看免费观看| 天堂影院成人在线观看| 一区二区三区高清视频在线| 亚洲精品美女久久久久99蜜臀| 99国产综合亚洲精品| 国产一卡二卡三卡精品| 国产欧美日韩一区二区三| 亚洲精品粉嫩美女一区| 国产精华一区二区三区| 夜夜夜夜夜久久久久| 久久伊人香网站| 成人一区二区视频在线观看| 亚洲国产精品合色在线| 在线观看www视频免费| 男人舔奶头视频| 国产主播在线观看一区二区| 十八禁网站免费在线| 国产三级黄色录像| 极品教师在线免费播放| 精品久久久久久久末码| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人一区二区三区免费视频网站| 成年女人毛片免费观看观看9| 在线天堂中文资源库| 欧美在线一区亚洲| a在线观看视频网站| 精品免费久久久久久久清纯| 天堂动漫精品| 欧美性猛交╳xxx乱大交人| 亚洲一码二码三码区别大吗| 男人操女人黄网站| 中文字幕精品亚洲无线码一区 | 法律面前人人平等表现在哪些方面| 观看免费一级毛片| 在线播放国产精品三级| 久久精品国产99精品国产亚洲性色| 日韩欧美国产在线观看| 日日爽夜夜爽网站| 国产精品乱码一区二三区的特点| 亚洲欧美精品综合一区二区三区| 亚洲 国产 在线| 97人妻精品一区二区三区麻豆 | 男女之事视频高清在线观看| 日本a在线网址| 国产真实乱freesex| 久久久国产成人精品二区| 国产成人av激情在线播放| 亚洲国产欧美网| 亚洲精品在线美女| 国产精品综合久久久久久久免费| 国产成人系列免费观看| 亚洲精品在线观看二区| 村上凉子中文字幕在线| 12—13女人毛片做爰片一| 久久香蕉精品热| 老司机在亚洲福利影院| 色婷婷久久久亚洲欧美| 高潮久久久久久久久久久不卡| 女同久久另类99精品国产91| 色综合亚洲欧美另类图片| 久久人妻福利社区极品人妻图片| 后天国语完整版免费观看| 很黄的视频免费| 99在线视频只有这里精品首页| 日韩高清综合在线| 男女视频在线观看网站免费 | 亚洲av电影在线进入| 亚洲av电影在线进入| 日韩高清综合在线| 亚洲人成电影免费在线| 级片在线观看| 高清在线国产一区| 精品国产乱子伦一区二区三区| 久久青草综合色| 欧美乱色亚洲激情| 亚洲一区二区三区色噜噜| 国产精品久久久久久亚洲av鲁大| 啦啦啦免费观看视频1| 久久久久久久久免费视频了| 91大片在线观看| 看免费av毛片| 又黄又粗又硬又大视频| 国产精品乱码一区二三区的特点| 女性被躁到高潮视频| 天天一区二区日本电影三级| 老司机靠b影院| 在线观看免费视频日本深夜| 欧美精品亚洲一区二区| 欧美亚洲日本最大视频资源| 国产主播在线观看一区二区| 黄频高清免费视频| 亚洲欧美激情综合另类| 中文字幕久久专区| 韩国av一区二区三区四区| 一本大道久久a久久精品| 国产熟女xx| 超碰成人久久| 老熟妇乱子伦视频在线观看| 亚洲 欧美 日韩 在线 免费| 久久香蕉精品热| 人人妻人人澡欧美一区二区| 一级a爱片免费观看的视频| 精品人妻1区二区| 18禁国产床啪视频网站| 久久久久久久精品吃奶| 一进一出抽搐动态| 亚洲第一电影网av| 精品一区二区三区av网在线观看| 麻豆成人av在线观看| 久久人人精品亚洲av| 久久久久久大精品| 啪啪无遮挡十八禁网站| 欧美中文日本在线观看视频| 亚洲欧美日韩高清在线视频| 日韩中文字幕欧美一区二区| 18禁黄网站禁片免费观看直播| 国产精品爽爽va在线观看网站 | 国产成人精品无人区| 久久精品91无色码中文字幕| 99精品久久久久人妻精品| 特大巨黑吊av在线直播 | 亚洲av美国av| 久久亚洲真实| 女警被强在线播放| 日本一本二区三区精品| 精品卡一卡二卡四卡免费| 欧美三级亚洲精品| 国产高清视频在线播放一区| 99热这里只有精品一区 | 成人特级黄色片久久久久久久| 伊人久久大香线蕉亚洲五| 国产91精品成人一区二区三区| av在线播放免费不卡| 国产精品九九99| 欧美色视频一区免费| www.精华液| 成人一区二区视频在线观看| 99久久99久久久精品蜜桃| 欧美 亚洲 国产 日韩一| 免费在线观看影片大全网站| 天天躁夜夜躁狠狠躁躁| 免费av毛片视频| 一夜夜www| 亚洲av第一区精品v没综合| 男女之事视频高清在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色片一级片一级黄色片| 欧美一级毛片孕妇| 国产片内射在线| 中文资源天堂在线| 国产精品99久久99久久久不卡| 亚洲中文字幕日韩| 国产精品影院久久| 一区福利在线观看| 黄色视频,在线免费观看| 男人操女人黄网站| 久久人妻av系列| 久久国产亚洲av麻豆专区| 精品久久久久久久久久免费视频| 国产亚洲精品一区二区www| 亚洲成人免费电影在线观看| 亚洲精品中文字幕在线视频| 亚洲自拍偷在线| 人人妻,人人澡人人爽秒播| 国产私拍福利视频在线观看| 欧美激情极品国产一区二区三区| 日韩av在线大香蕉| 性欧美人与动物交配| 午夜日韩欧美国产| 亚洲成人免费电影在线观看| 在线视频色国产色| 午夜久久久在线观看| 国产高清videossex| 欧美最黄视频在线播放免费| 欧美日韩瑟瑟在线播放| 此物有八面人人有两片| 18禁国产床啪视频网站| 国产欧美日韩精品亚洲av| 国产亚洲欧美在线一区二区| 免费人成视频x8x8入口观看| 满18在线观看网站| 两个人视频免费观看高清| 日韩有码中文字幕| 99久久国产精品久久久| 欧美zozozo另类| 久久久久久国产a免费观看| 18禁美女被吸乳视频| 18美女黄网站色大片免费观看| 精品日产1卡2卡| 可以在线观看毛片的网站| 久久精品91蜜桃| 日本 av在线| 国内毛片毛片毛片毛片毛片| 麻豆av在线久日| 美女 人体艺术 gogo| 黑人欧美特级aaaaaa片| 中文资源天堂在线| 18禁观看日本| 宅男免费午夜| 亚洲黑人精品在线| 麻豆国产av国片精品| 国产成人啪精品午夜网站| 久久 成人 亚洲| 亚洲 欧美 日韩 在线 免费| 亚洲片人在线观看| 欧美性猛交黑人性爽| 欧美在线一区亚洲| av福利片在线| www.精华液| 麻豆一二三区av精品| 一级a爱视频在线免费观看| 神马国产精品三级电影在线观看 | 两性午夜刺激爽爽歪歪视频在线观看 | 婷婷六月久久综合丁香| 亚洲精品粉嫩美女一区| 一个人观看的视频www高清免费观看 | 午夜免费激情av| 免费看日本二区| 高潮久久久久久久久久久不卡| 91九色精品人成在线观看| 婷婷精品国产亚洲av| 免费看a级黄色片| 国产91精品成人一区二区三区| 亚洲中文av在线| av超薄肉色丝袜交足视频| 国产麻豆成人av免费视频| 级片在线观看| 美国免费a级毛片| 999精品在线视频| 成人三级黄色视频| 露出奶头的视频| 亚洲一码二码三码区别大吗| 亚洲人成电影免费在线| 熟女少妇亚洲综合色aaa.| 黄片播放在线免费| 国产视频一区二区在线看| 亚洲精品国产精品久久久不卡| 精品一区二区三区四区五区乱码| 国产亚洲精品一区二区www| 国产一区在线观看成人免费| 欧美黑人巨大hd| 日韩一卡2卡3卡4卡2021年| www.熟女人妻精品国产| 欧美日本视频| 悠悠久久av| 香蕉av资源在线| 亚洲人成77777在线视频| 久久久久国内视频| 久久精品国产99精品国产亚洲性色| 国产精品野战在线观看| 一边摸一边做爽爽视频免费| 久久精品夜夜夜夜夜久久蜜豆 | 男女午夜视频在线观看| 99精品欧美一区二区三区四区| 亚洲av第一区精品v没综合| 嫁个100分男人电影在线观看| 特大巨黑吊av在线直播 | 男女做爰动态图高潮gif福利片| 黄色视频不卡| 男女视频在线观看网站免费 | 午夜精品久久久久久毛片777| 亚洲国产精品久久男人天堂| 欧美绝顶高潮抽搐喷水| 大型黄色视频在线免费观看| 亚洲熟妇熟女久久| 中国美女看黄片| 精品无人区乱码1区二区| 丰满人妻熟妇乱又伦精品不卡| 少妇 在线观看| 在线免费观看的www视频| 成年人黄色毛片网站| 国产在线精品亚洲第一网站| 国产亚洲精品久久久久久毛片| av在线播放免费不卡| 1024视频免费在线观看| 高潮久久久久久久久久久不卡| 亚洲五月色婷婷综合| 久久99热这里只有精品18| 久9热在线精品视频| 美国免费a级毛片| avwww免费| 久久天躁狠狠躁夜夜2o2o| 50天的宝宝边吃奶边哭怎么回事| 嫩草影视91久久| 免费人成视频x8x8入口观看| 国产在线精品亚洲第一网站| 欧美一区二区精品小视频在线| 一夜夜www| 一级毛片精品| 男女下面进入的视频免费午夜 | 黄片播放在线免费| 亚洲 国产 在线| 亚洲国产精品成人综合色| 免费在线观看日本一区| 久久国产乱子伦精品免费另类| 久久久久国内视频| 国产成人精品久久二区二区免费| avwww免费| 中文字幕人成人乱码亚洲影| 国内精品久久久久久久电影| 韩国精品一区二区三区| x7x7x7水蜜桃| 无遮挡黄片免费观看| 在线看三级毛片| 韩国精品一区二区三区| 夜夜夜夜夜久久久久| 亚洲精品一区av在线观看| xxxwww97欧美| 精品国产超薄肉色丝袜足j| 18禁美女被吸乳视频| 亚洲人成网站高清观看| 成人av一区二区三区在线看| 中文字幕另类日韩欧美亚洲嫩草| 欧美国产日韩亚洲一区| 国产av在哪里看| 又大又爽又粗| 麻豆一二三区av精品| 国产av一区二区精品久久| 黄片大片在线免费观看| 久久青草综合色| 日韩 欧美 亚洲 中文字幕| 母亲3免费完整高清在线观看| 91老司机精品| 91成人精品电影| 青草久久国产| 在线十欧美十亚洲十日本专区| 国产一卡二卡三卡精品| e午夜精品久久久久久久| 国产精品亚洲一级av第二区| 琪琪午夜伦伦电影理论片6080| 欧美中文日本在线观看视频| 亚洲人成网站在线播放欧美日韩| 19禁男女啪啪无遮挡网站| 国产单亲对白刺激| 亚洲熟女毛片儿| 丝袜美腿诱惑在线| 亚洲天堂国产精品一区在线| 亚洲全国av大片| 1024视频免费在线观看| 国产成人精品久久二区二区免费| 午夜福利一区二区在线看| 免费在线观看黄色视频的| 欧美性长视频在线观看| 亚洲国产欧洲综合997久久, | 在线观看66精品国产| 99国产综合亚洲精品| 最新美女视频免费是黄的| 最好的美女福利视频网| 欧美成狂野欧美在线观看| 白带黄色成豆腐渣| 欧美性猛交黑人性爽| 制服人妻中文乱码| av中文乱码字幕在线| 日韩av在线大香蕉| 亚洲国产欧美网| 久久亚洲精品不卡| 欧美色欧美亚洲另类二区| 亚洲精品国产区一区二| 亚洲av熟女| 看黄色毛片网站| 99热只有精品国产| 日韩欧美在线二视频| 日韩成人在线观看一区二区三区| 午夜精品在线福利| 在线观看66精品国产| 欧美在线一区亚洲| 超碰成人久久| 丝袜美腿诱惑在线| 国产亚洲精品综合一区在线观看 | 51午夜福利影视在线观看| 亚洲av成人一区二区三| 他把我摸到了高潮在线观看| 精品久久蜜臀av无| 久久久久亚洲av毛片大全| 青草久久国产| 欧美精品亚洲一区二区| 制服丝袜大香蕉在线| 麻豆一二三区av精品| 国产成年人精品一区二区| 中出人妻视频一区二区| 韩国精品一区二区三区| 禁无遮挡网站| 美国免费a级毛片| 亚洲精品色激情综合| 国产区一区二久久| 亚洲中文字幕日韩| 国产亚洲精品久久久久5区| 日本一本二区三区精品| 最好的美女福利视频网| 国产av又大| 国产精品1区2区在线观看.| 琪琪午夜伦伦电影理论片6080| 亚洲天堂国产精品一区在线| 自线自在国产av| 夜夜爽天天搞| 人人澡人人妻人| 亚洲自拍偷在线| 大香蕉久久成人网| 成人永久免费在线观看视频| 久久久久久久久久黄片| 国产成人啪精品午夜网站| 国产97色在线日韩免费| 亚洲一区二区三区色噜噜| 亚洲国产精品成人综合色| 超碰成人久久| 男人舔女人的私密视频| 91成年电影在线观看| 窝窝影院91人妻| 在线观看一区二区三区| 国产单亲对白刺激| 2021天堂中文幕一二区在线观 | 9191精品国产免费久久| 一二三四社区在线视频社区8| 制服人妻中文乱码| 他把我摸到了高潮在线观看| 校园春色视频在线观看| 成人三级黄色视频| 亚洲成a人片在线一区二区| 久久久久亚洲av毛片大全| 欧美日韩福利视频一区二区| 亚洲精品中文字幕在线视频| 一本久久中文字幕| 一边摸一边做爽爽视频免费| 亚洲国产欧美一区二区综合| 手机成人av网站| 亚洲av熟女| 两性午夜刺激爽爽歪歪视频在线观看 | 天天躁狠狠躁夜夜躁狠狠躁| 美女 人体艺术 gogo| 91字幕亚洲| 欧美日韩一级在线毛片| 少妇粗大呻吟视频| 久99久视频精品免费| 一区二区三区高清视频在线| 亚洲成人久久爱视频| 美女大奶头视频| 制服人妻中文乱码| 黑丝袜美女国产一区| 国语自产精品视频在线第100页| 久久午夜亚洲精品久久| 亚洲自拍偷在线| 亚洲天堂国产精品一区在线| 成人亚洲精品一区在线观看| 亚洲午夜精品一区,二区,三区| 给我免费播放毛片高清在线观看| 99国产极品粉嫩在线观看| www.精华液| 亚洲精品色激情综合| 在线国产一区二区在线| 人成视频在线观看免费观看| 丝袜在线中文字幕| 性欧美人与动物交配| 宅男免费午夜| 亚洲午夜理论影院| 精品久久久久久久末码| 999久久久国产精品视频| 成熟少妇高潮喷水视频| 亚洲av美国av| 精品国产国语对白av| 母亲3免费完整高清在线观看| 国产成人欧美| 国产单亲对白刺激| 久久精品国产亚洲av香蕉五月| www.精华液| 久久久久久人人人人人| 亚洲熟妇中文字幕五十中出| 50天的宝宝边吃奶边哭怎么回事| 男女床上黄色一级片免费看| 国产精品野战在线观看| 免费高清在线观看日韩| 天天一区二区日本电影三级| 久久天躁狠狠躁夜夜2o2o| 亚洲国产高清在线一区二区三 | 欧洲精品卡2卡3卡4卡5卡区| 精品国产美女av久久久久小说| 中文在线观看免费www的网站 | 免费高清视频大片| 午夜福利成人在线免费观看| 啦啦啦免费观看视频1| 久久久久久亚洲精品国产蜜桃av| 亚洲专区中文字幕在线| 禁无遮挡网站| 精品国产亚洲在线| 亚洲五月婷婷丁香| 身体一侧抽搐| 国产av又大| 久久精品91蜜桃| 一本久久中文字幕| 又紧又爽又黄一区二区| 欧美黑人巨大hd| 777久久人妻少妇嫩草av网站| 一个人观看的视频www高清免费观看 | 亚洲国产欧美日韩在线播放| 精品熟女少妇八av免费久了| 欧美一区二区精品小视频在线| 黄色片一级片一级黄色片| av中文乱码字幕在线| 狂野欧美激情性xxxx| 久久狼人影院| 亚洲欧洲精品一区二区精品久久久| 欧美激情 高清一区二区三区| 高潮久久久久久久久久久不卡| 久久香蕉国产精品| 亚洲午夜理论影院| 51午夜福利影视在线观看| 777久久人妻少妇嫩草av网站| 桃红色精品国产亚洲av| 悠悠久久av| 精品一区二区三区四区五区乱码| 亚洲国产毛片av蜜桃av| 亚洲精品国产一区二区精华液| 午夜免费观看网址| 91大片在线观看| 亚洲狠狠婷婷综合久久图片| 侵犯人妻中文字幕一二三四区| 一进一出抽搐动态| 久久精品aⅴ一区二区三区四区| 黄色毛片三级朝国网站| 成人午夜高清在线视频 | 中文字幕精品亚洲无线码一区 | 中文字幕精品亚洲无线码一区 | 一进一出抽搐动态| 一区二区日韩欧美中文字幕| 看片在线看免费视频| 国产精品久久久久久人妻精品电影| 久久久久久九九精品二区国产 | 精品国产国语对白av| 99riav亚洲国产免费| 最好的美女福利视频网| 91大片在线观看| av在线天堂中文字幕| 日本黄色视频三级网站网址| 97超级碰碰碰精品色视频在线观看| 亚洲第一青青草原| 91在线观看av| 在线观看免费日韩欧美大片| 久久国产乱子伦精品免费另类| 黄片播放在线免费| 狠狠狠狠99中文字幕| 国产高清视频在线播放一区| 熟女少妇亚洲综合色aaa.| 国产午夜精品久久久久久| 曰老女人黄片| 国产国语露脸激情在线看| 日韩欧美一区二区三区在线观看| 午夜日韩欧美国产| 嫁个100分男人电影在线观看| 午夜影院日韩av| 我的亚洲天堂| 欧美在线一区亚洲| 日韩精品青青久久久久久| 精品久久久久久久久久久久久 | 亚洲第一电影网av| 久久久精品欧美日韩精品| 一区二区三区激情视频| 久久久久国产一级毛片高清牌| 久久久久久国产a免费观看| 国产精品亚洲一级av第二区| 首页视频小说图片口味搜索| 天堂动漫精品| 白带黄色成豆腐渣| 非洲黑人性xxxx精品又粗又长| 黄片小视频在线播放|