• <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)用
    亚洲av成人精品一二三区| 老女人水多毛片| 亚洲,欧美,日韩| 高清黄色对白视频在线免费看 | 国产精品一区www在线观看| 熟妇人妻不卡中文字幕| 国产一区二区三区综合在线观看 | 三级经典国产精品| 精品国产一区二区久久| 亚洲成人手机| 晚上一个人看的免费电影| 欧美高清成人免费视频www| 熟女电影av网| 亚洲精品乱久久久久久| 中文字幕精品免费在线观看视频 | 啦啦啦啦在线视频资源| 国产伦理片在线播放av一区| 亚洲精品自拍成人| 九色成人免费人妻av| 夜夜骑夜夜射夜夜干| 熟女av电影| 91久久精品电影网| 国国产精品蜜臀av免费| 久久人人爽人人片av| 国产淫语在线视频| 黑丝袜美女国产一区| 欧美日韩精品成人综合77777| 青春草亚洲视频在线观看| 水蜜桃什么品种好| 人人妻人人看人人澡| 欧美精品人与动牲交sv欧美| 在线观看三级黄色| 亚洲欧美中文字幕日韩二区| 国产精品人妻久久久久久| 91久久精品电影网| 国产精品一区二区在线不卡| 性色av一级| 欧美高清成人免费视频www| 日韩av不卡免费在线播放| 国产精品免费大片| 99热全是精品| 欧美激情国产日韩精品一区| 欧美日韩精品成人综合77777| 少妇裸体淫交视频免费看高清| 午夜激情福利司机影院| 国产探花极品一区二区| 久久精品久久久久久噜噜老黄| 午夜免费鲁丝| 亚洲美女搞黄在线观看| 欧美区成人在线视频| 精品国产一区二区久久| 美女cb高潮喷水在线观看| 日日啪夜夜爽| 欧美xxⅹ黑人| 一级毛片电影观看| 亚洲av综合色区一区| 久久av网站| 国产精品成人在线| 偷拍熟女少妇极品色| 秋霞在线观看毛片| 免费看光身美女| 中文欧美无线码| 观看美女的网站| 日本黄色片子视频| 亚洲经典国产精华液单| 成人国产麻豆网| 久久 成人 亚洲| 晚上一个人看的免费电影| 看免费成人av毛片| 丰满人妻一区二区三区视频av| 午夜老司机福利剧场| 又爽又黄a免费视频| 夫妻性生交免费视频一级片| 亚洲自偷自拍三级| 国产日韩欧美视频二区| 日韩,欧美,国产一区二区三区| 国产真实伦视频高清在线观看| 久久人人爽av亚洲精品天堂| 老熟女久久久| 美女脱内裤让男人舔精品视频| 我的女老师完整版在线观看| 黄色视频在线播放观看不卡| 亚洲,欧美,日韩| 国产熟女午夜一区二区三区 | 亚洲国产欧美日韩在线播放 | 国产亚洲精品久久久com| av在线老鸭窝| 蜜臀久久99精品久久宅男| 国产精品一区二区在线不卡| 免费黄网站久久成人精品| 亚洲中文av在线| 亚洲国产av新网站| 久久久久网色| 国产午夜精品一二区理论片| 精华霜和精华液先用哪个| 久久青草综合色| 久久久久久久国产电影| 国产av码专区亚洲av| 亚州av有码| 日本欧美国产在线视频| 久久av网站| 观看av在线不卡| 卡戴珊不雅视频在线播放| 成年av动漫网址| 国产高清不卡午夜福利| 啦啦啦在线观看免费高清www| 国产精品久久久久久精品电影小说| 在线观看免费视频网站a站| 男的添女的下面高潮视频| 免费av中文字幕在线| 久久精品国产亚洲网站| 老司机影院毛片| 一本—道久久a久久精品蜜桃钙片| 久久国内精品自在自线图片| 国产精品人妻久久久久久| 欧美日韩视频精品一区| 久久韩国三级中文字幕| 99热这里只有是精品50| av播播在线观看一区| 黄色毛片三级朝国网站 | 国产精品三级大全| 国产精品久久久久久久久免| 全区人妻精品视频| 久久影院123| 国产无遮挡羞羞视频在线观看| 黄色一级大片看看| 大片电影免费在线观看免费| 国产精品久久久久久精品古装| 国产欧美亚洲国产| 自拍欧美九色日韩亚洲蝌蚪91 | 日本色播在线视频| 免费观看的影片在线观看| 成人综合一区亚洲| 69精品国产乱码久久久| 热99国产精品久久久久久7| 精品一品国产午夜福利视频| 亚洲精品亚洲一区二区| 精品少妇内射三级| 久久久久久久亚洲中文字幕| 欧美亚洲 丝袜 人妻 在线| 色视频在线一区二区三区| 中文字幕免费在线视频6| 国产高清不卡午夜福利| 在线观看国产h片| 亚洲精品国产av蜜桃| 夫妻性生交免费视频一级片| 国产精品一区二区在线不卡| 日本91视频免费播放| 精品人妻偷拍中文字幕| 99久久综合免费| 日本欧美视频一区| 国产亚洲5aaaaa淫片| 久久97久久精品| 久久精品国产鲁丝片午夜精品| 国产成人精品久久久久久| a级毛片免费高清观看在线播放| 亚洲精品中文字幕在线视频 | 成人漫画全彩无遮挡| 嘟嘟电影网在线观看| 色94色欧美一区二区| 亚洲成人av在线免费| 欧美精品国产亚洲| 建设人人有责人人尽责人人享有的| 亚洲精品一二三| 国产一区二区在线观看日韩| 亚洲欧美清纯卡通| 亚洲性久久影院| 久久久久精品久久久久真实原创| 丝袜脚勾引网站| 伦精品一区二区三区| 日本vs欧美在线观看视频 | 国产一区二区三区av在线| 99九九在线精品视频 | 久久97久久精品| 高清毛片免费看| 天天操日日干夜夜撸| 日韩成人av中文字幕在线观看| 久久久久久久久久久丰满| 交换朋友夫妻互换小说| 在线观看免费高清a一片| 丝袜喷水一区| 精品国产乱码久久久久久小说| 久久国内精品自在自线图片| 久久鲁丝午夜福利片| 五月玫瑰六月丁香| 久久久久久人妻| 三级国产精品欧美在线观看| 久久97久久精品| 国产成人精品一,二区| 久久久午夜欧美精品| 99久久综合免费| 久热这里只有精品99| 国产欧美日韩综合在线一区二区 | 女性生殖器流出的白浆| 搡老乐熟女国产| 自线自在国产av| 国产女主播在线喷水免费视频网站| 亚洲精品久久午夜乱码| 日韩精品免费视频一区二区三区 | 十八禁网站网址无遮挡 | 中国美白少妇内射xxxbb| 国产黄色免费在线视频| 国产精品三级大全| 国产日韩欧美视频二区| 亚洲欧美一区二区三区黑人 | 国产av码专区亚洲av| 国产高清国产精品国产三级| 精品亚洲成国产av| 一区在线观看完整版| 一本一本综合久久| 大香蕉久久网| 男女边吃奶边做爰视频| 中文字幕制服av| 久久99一区二区三区| 一级二级三级毛片免费看| 乱人伦中国视频| 日韩人妻高清精品专区| 日韩免费高清中文字幕av| 少妇被粗大的猛进出69影院 | 夫妻性生交免费视频一级片| 国产高清国产精品国产三级| 男人添女人高潮全过程视频| www.色视频.com| 麻豆乱淫一区二区| 国产精品久久久久久久电影| 亚洲色图综合在线观看| 亚洲人成网站在线观看播放| 亚洲一级一片aⅴ在线观看| 妹子高潮喷水视频| 亚洲激情五月婷婷啪啪| 亚洲国产最新在线播放| 一个人看视频在线观看www免费| 国产亚洲91精品色在线| 自线自在国产av| 亚洲精品乱码久久久久久按摩| 久久国产亚洲av麻豆专区| 成年女人在线观看亚洲视频| 国产精品国产三级国产av玫瑰| 性色av一级| 在线天堂最新版资源| 国产综合精华液| 天堂俺去俺来也www色官网| 国产在线视频一区二区| av一本久久久久| 丝袜脚勾引网站| 一个人看视频在线观看www免费| 午夜激情久久久久久久| 乱人伦中国视频| 国产黄频视频在线观看| 国产精品一区www在线观看| 亚洲av二区三区四区| 欧美少妇被猛烈插入视频| 亚洲精品一区蜜桃| 国产一区二区在线观看日韩| 久久精品久久久久久久性| 寂寞人妻少妇视频99o| 国产高清有码在线观看视频| 亚洲精品一二三| 久久久久视频综合| 久久精品国产a三级三级三级| 汤姆久久久久久久影院中文字幕| 国产淫片久久久久久久久| 亚洲成色77777| 激情五月婷婷亚洲| 男人添女人高潮全过程视频| 久久久久视频综合| 新久久久久国产一级毛片| 成人黄色视频免费在线看| 国产成人aa在线观看| 九九在线视频观看精品| 久久午夜福利片| 一个人看视频在线观看www免费| 国产精品三级大全| 日韩电影二区| av又黄又爽大尺度在线免费看| 亚洲国产欧美日韩在线播放 | 自线自在国产av| 亚洲欧美精品专区久久| 麻豆精品久久久久久蜜桃| 制服丝袜香蕉在线| 免费观看的影片在线观看| 国产精品.久久久| 国产日韩一区二区三区精品不卡 | 久久av网站| 九色成人免费人妻av| 中文在线观看免费www的网站| 欧美日韩精品成人综合77777| 国产中年淑女户外野战色| 全区人妻精品视频| 一区二区三区乱码不卡18| 九草在线视频观看| 校园人妻丝袜中文字幕| 亚洲人成网站在线观看播放| 欧美另类一区| 国产av一区二区精品久久| 桃花免费在线播放| 乱人伦中国视频| 国产熟女午夜一区二区三区 | 高清视频免费观看一区二区| 91久久精品国产一区二区三区| 在线观看人妻少妇| 七月丁香在线播放| 国产日韩欧美亚洲二区| 国产精品秋霞免费鲁丝片| 嘟嘟电影网在线观看| 2022亚洲国产成人精品| 午夜激情福利司机影院| 亚洲av中文av极速乱| 色婷婷久久久亚洲欧美| 亚洲av日韩在线播放| 亚洲成人手机| av福利片在线观看| 毛片一级片免费看久久久久| www.av在线官网国产| 午夜福利在线观看免费完整高清在| 啦啦啦视频在线资源免费观看| 国产毛片在线视频| 国产综合精华液| 蜜桃在线观看..| 大片电影免费在线观看免费| 大又大粗又爽又黄少妇毛片口| 国产成人免费观看mmmm| 精品亚洲成国产av| 一本色道久久久久久精品综合| 亚洲精品乱码久久久v下载方式| 少妇人妻一区二区三区视频| 国产精品久久久久久精品古装| 国产极品天堂在线| 免费在线观看成人毛片| 午夜日本视频在线| 亚洲精品第二区| 在线观看av片永久免费下载| 狂野欧美白嫩少妇大欣赏| 成年美女黄网站色视频大全免费 | 狠狠精品人妻久久久久久综合| 久久99精品国语久久久| 韩国av在线不卡| 国产女主播在线喷水免费视频网站| 一本大道久久a久久精品| 大陆偷拍与自拍| 美女内射精品一级片tv| 欧美区成人在线视频| 六月丁香七月| 精品国产露脸久久av麻豆| 最近中文字幕高清免费大全6| 国产伦在线观看视频一区| 在线观看国产h片| 午夜免费男女啪啪视频观看| 国产日韩一区二区三区精品不卡 | 国产黄片美女视频| 成人无遮挡网站| 色94色欧美一区二区| 在线观看免费高清a一片| 国产在线免费精品| 国产成人免费无遮挡视频| 黄色配什么色好看| 国产免费视频播放在线视频| 国产精品一区二区在线观看99| 久久青草综合色| 精品少妇黑人巨大在线播放| 日日摸夜夜添夜夜爱| 精品少妇黑人巨大在线播放| 中文字幕av电影在线播放| 十八禁网站网址无遮挡 | 日本-黄色视频高清免费观看| 精品久久久精品久久久| 欧美xxⅹ黑人| 精品人妻熟女av久视频| 亚洲精品视频女| 国产精品99久久99久久久不卡 | 精品一区在线观看国产| 国产又色又爽无遮挡免| 丰满迷人的少妇在线观看| 久久亚洲国产成人精品v| 久久影院123| 日日摸夜夜添夜夜添av毛片| 日本av手机在线免费观看| 免费看光身美女| 亚洲精品久久午夜乱码| 性高湖久久久久久久久免费观看| 久久午夜福利片| 99精国产麻豆久久婷婷| 日本欧美国产在线视频| 久久热精品热| 丰满饥渴人妻一区二区三| 国产亚洲5aaaaa淫片| 国产亚洲一区二区精品| 观看av在线不卡| 国产爽快片一区二区三区| 亚洲美女视频黄频| 又爽又黄a免费视频| 桃花免费在线播放| 大又大粗又爽又黄少妇毛片口| 日韩av免费高清视频| 日韩电影二区| 中国国产av一级| a级片在线免费高清观看视频| 久久午夜综合久久蜜桃| 日本黄色日本黄色录像| 美女福利国产在线| 最近中文字幕2019免费版| 亚洲精品国产av成人精品| 亚洲第一区二区三区不卡| 色婷婷久久久亚洲欧美| 亚洲精品亚洲一区二区| 一级av片app| 久久精品久久久久久噜噜老黄| 99久久人妻综合| 久久久久久伊人网av| 国产一区有黄有色的免费视频| 爱豆传媒免费全集在线观看| 日韩欧美一区视频在线观看 | 国产精品一区二区三区四区免费观看| 亚洲高清免费不卡视频| 免费观看在线日韩| 狂野欧美白嫩少妇大欣赏| 免费高清在线观看视频在线观看| 另类精品久久| 亚洲欧洲国产日韩| 丁香六月天网| 国产国拍精品亚洲av在线观看| 久久婷婷青草| 中国国产av一级| 国产成人freesex在线| 美女福利国产在线| 永久网站在线| 午夜激情福利司机影院| 欧美日韩精品成人综合77777| 老熟女久久久| 国产男人的电影天堂91| 人妻 亚洲 视频| 美女xxoo啪啪120秒动态图| 亚洲国产日韩一区二区| 国产免费一区二区三区四区乱码| 三级国产精品欧美在线观看| 日韩一区二区三区影片| 久久久久精品久久久久真实原创| 国产精品麻豆人妻色哟哟久久| 欧美精品国产亚洲| 久久亚洲国产成人精品v| 国产高清有码在线观看视频| 国产有黄有色有爽视频| 国产美女午夜福利| 婷婷色麻豆天堂久久| 丝袜在线中文字幕| 极品人妻少妇av视频| 成人美女网站在线观看视频| 在线观看国产h片| 成年人免费黄色播放视频 | 国产欧美亚洲国产| 亚洲精品国产av成人精品| av网站免费在线观看视频| 亚洲欧美日韩另类电影网站| 成人黄色视频免费在线看| 久久精品国产a三级三级三级| 少妇猛男粗大的猛烈进出视频| 久久久国产精品麻豆| 大香蕉97超碰在线| 汤姆久久久久久久影院中文字幕| 午夜精品国产一区二区电影| 黑人高潮一二区| 亚洲人成网站在线观看播放| 亚洲,一卡二卡三卡| 夫妻性生交免费视频一级片| 亚洲欧美中文字幕日韩二区| 少妇人妻一区二区三区视频| 国产日韩欧美亚洲二区| 一级毛片电影观看| 国产在线免费精品| 女性被躁到高潮视频| 只有这里有精品99| 久久精品国产鲁丝片午夜精品| 一级片'在线观看视频| 中文资源天堂在线| 18禁在线无遮挡免费观看视频| 在线看a的网站| 国产黄色视频一区二区在线观看| 国产精品一区二区性色av| 色婷婷久久久亚洲欧美| 国产欧美亚洲国产| freevideosex欧美| 国产精品久久久久久精品电影小说| 天天操日日干夜夜撸| 美女cb高潮喷水在线观看| 不卡视频在线观看欧美| 丰满乱子伦码专区| 成人午夜精彩视频在线观看| 亚洲av成人精品一二三区| 久久精品国产a三级三级三级| 99久久精品热视频| 久久久久久久国产电影| 国产欧美日韩综合在线一区二区 | 国产av一区二区精品久久| 亚洲欧洲日产国产| 亚洲三级黄色毛片| 久久久久久伊人网av| 亚洲精品乱码久久久v下载方式| 观看av在线不卡| 国产精品国产三级国产专区5o| 日本-黄色视频高清免费观看| 午夜老司机福利剧场| 久久久国产精品麻豆| 视频区图区小说| 日韩亚洲欧美综合| av福利片在线观看| 能在线免费看毛片的网站| 欧美最新免费一区二区三区| 热re99久久国产66热| 大香蕉久久网| 亚洲精品久久午夜乱码| 日韩成人伦理影院| √禁漫天堂资源中文www| 一级毛片aaaaaa免费看小| 中国三级夫妇交换| 午夜精品国产一区二区电影| 欧美精品一区二区大全| 日韩亚洲欧美综合| 你懂的网址亚洲精品在线观看| 国产成人免费无遮挡视频| 国产成人freesex在线| 一级二级三级毛片免费看| 99久久精品一区二区三区| 丰满人妻一区二区三区视频av| av线在线观看网站| 一区二区av电影网| 成年av动漫网址| av福利片在线| 妹子高潮喷水视频| 一区二区三区精品91| 狂野欧美白嫩少妇大欣赏| 高清黄色对白视频在线免费看 | 国产毛片在线视频| 97在线视频观看| 日韩欧美一区视频在线观看 | 亚洲精品一区蜜桃| 国产精品秋霞免费鲁丝片| 亚洲成人一二三区av| 国内少妇人妻偷人精品xxx网站| 国内精品宾馆在线| 国产男女超爽视频在线观看| 热re99久久国产66热| 看免费成人av毛片| 国产精品女同一区二区软件| 国产色婷婷99| 亚洲怡红院男人天堂| 欧美日韩亚洲高清精品| 91aial.com中文字幕在线观看| 国产黄色视频一区二区在线观看| 欧美日韩一区二区视频在线观看视频在线| 91久久精品国产一区二区成人| 亚洲精品一二三| 欧美国产精品一级二级三级 | 又大又黄又爽视频免费| 欧美丝袜亚洲另类| 亚洲精品自拍成人| 亚洲欧美一区二区三区国产| 99精国产麻豆久久婷婷| 久久韩国三级中文字幕| 国产一区二区三区综合在线观看 | 午夜视频国产福利| 五月玫瑰六月丁香| videos熟女内射| 搡女人真爽免费视频火全软件| 热re99久久精品国产66热6| 亚洲精品自拍成人| 日韩 亚洲 欧美在线| 永久免费av网站大全| 91精品国产国语对白视频| 亚洲中文av在线| av国产精品久久久久影院| 99国产精品免费福利视频| 在线观看美女被高潮喷水网站| 两个人免费观看高清视频 | 精品一品国产午夜福利视频| 日韩三级伦理在线观看| 久久久国产欧美日韩av| 国内揄拍国产精品人妻在线| 婷婷色av中文字幕| 这个男人来自地球电影免费观看 | 九九在线视频观看精品| 男男h啪啪无遮挡| 色94色欧美一区二区| 久久国产亚洲av麻豆专区| 亚洲中文av在线| 成人漫画全彩无遮挡| 日韩一区二区视频免费看| 精品久久久久久久久av| 欧美一级a爱片免费观看看| 精品国产一区二区久久| 国产在线一区二区三区精| 美女主播在线视频| 熟妇人妻不卡中文字幕| 欧美日韩精品成人综合77777| 啦啦啦啦在线视频资源| 深夜a级毛片| 午夜精品国产一区二区电影| 亚洲欧美成人综合另类久久久| 久久精品国产a三级三级三级| 五月玫瑰六月丁香| 少妇的逼水好多| 亚洲情色 制服丝袜| 国产欧美亚洲国产| 国产免费一区二区三区四区乱码| 人体艺术视频欧美日本| 2021少妇久久久久久久久久久| av网站免费在线观看视频| 麻豆乱淫一区二区| 69精品国产乱码久久久| 国产日韩欧美亚洲二区| 免费久久久久久久精品成人欧美视频 | 国产亚洲一区二区精品| 国产熟女午夜一区二区三区 | a级毛色黄片| 国产极品粉嫩免费观看在线 | 大香蕉97超碰在线|