丁 楊 柳建華 葉方平 姜林林 鄂曉雪 吳 昊

(1上海理工大學(xué)能源與動力工程學(xué)院 上海 200093; 2 浙江新勁空調(diào)設(shè)備有限公司 龍泉 323700)

R134a水平微細(xì)管內(nèi)流動沸騰換熱的實驗研究

丁 楊1柳建華1葉方平2姜林林1鄂曉雪1吳 昊1

(1上海理工大學(xué)能源與動力工程學(xué)院 上海 200093; 2 浙江新勁空調(diào)設(shè)備有限公司 龍泉 323700)

本文對R134a在水平微細(xì)管內(nèi)的流動沸騰進(jìn)行了實驗研究。實驗測試段選用了內(nèi)徑為1 mm、2 mm、3 mm共3種不同的水平光滑不銹鋼管，實驗的飽和溫度為5～30 ℃，熱流密度為2～70 kW/m2，流量范圍為200～1500 kg/(m25s)。實驗結(jié)果表明：相同條件下，干涸前2 mm管較3 mm管換熱系數(shù)平均增幅為11.6%，1 mm管較2 mm管換熱增幅為26.3%，1 mm管徑換熱系數(shù)比3 mm管徑平均增大40.8%。隨著管徑的減小，換熱系數(shù)在更低的干度開始減小，質(zhì)量流速和強制對流蒸發(fā)作用對換熱系數(shù)的影響變小，熱流密度的影響依然顯著；塞狀流和彈狀流區(qū)域減小，泡狀流和環(huán)狀流區(qū)域增大。

R134a；微細(xì)管；流動沸騰；干涸；換熱系數(shù)

在能源日益緊缺的今天，高效、緊湊的換熱器由于在制冷、空調(diào)系統(tǒng)中表現(xiàn)出良好的節(jié)能效應(yīng)而日益受到大家的關(guān)注，微小型換熱器的研制也成為研究熱點。因此，研究微小尺度管徑下制冷劑的流動沸騰換熱特性，具有非常重要的實際意義。盡管從最早Lazarek G M & Black S H[1]對管內(nèi)流動沸騰的實驗到現(xiàn)在已經(jīng)超過30年，已經(jīng)有大量關(guān)于這一方面的研究，但是在微細(xì)管徑(管徑小于3 mm)下的流動沸騰換熱中，各個文獻(xiàn)之間的換熱系數(shù)趨勢與熱流密度、飽和溫度、質(zhì)量流速、干度等等因素的變化關(guān)系存在較大差異，普遍的換熱機理以及管徑、流型等因素對換熱的影響尚不明確[2-4]。

首先，微細(xì)通道內(nèi)流動沸騰換熱機理存在較大的爭議，根據(jù)Ong C L & Thome J R[5]總結(jié)主要的換熱機理主要分為四種：1)核態(tài)沸騰，換熱與熱流密度有著很強的關(guān)聯(lián)性，如Lazarek G M & Black S H[1]、Bao Z Y等[6]所述；2)對流沸騰，換熱與質(zhì)量流速和干度有很強的依賴性，與熱流密度關(guān)聯(lián)不大，如Qu W & Mudawar I[2]、Kuwahara K等[7]所述；3)核態(tài)沸騰和對流沸騰共存，即低干度區(qū)核態(tài)沸騰主導(dǎo)，高干度區(qū)對流沸騰主導(dǎo)，兩種機理此消彼長的混合換熱機理，如Lin S等[8-9]、Yan Y Y & Lin T F[10]所述；4)薄膜蒸發(fā)，即管內(nèi)加長氣泡與內(nèi)壁之間薄膜的蒸發(fā)換熱機理，與氣泡的產(chǎn)生頻率有關(guān)，進(jìn)而與熱流密度相關(guān)，如Consolini L & Thome J R[11]所述。其次，微細(xì)通道流動沸騰換熱的實驗結(jié)果有很大的差異[12-14]，甚至存在規(guī)律完全相反的情況。如Lin S等[8-9]研究R141b在不同管徑下的流動沸騰換熱結(jié)果顯示，在干涸點以前，換熱系數(shù)一直隨著干度的增加而增加；而Huo X[15]研究R134a水平圓管內(nèi)流動沸騰換熱時發(fā)現(xiàn)在較大的干度范圍內(nèi)，換熱系數(shù)隨著干度的增加而減小。

本文利用實驗方法，研究R134a在3種不同管徑水平管內(nèi)的流動沸騰換熱特性，重點在于系統(tǒng)地分析管道尺寸不斷減小的條件下，質(zhì)量流速、熱流密度以及流型對管內(nèi)流動沸騰換熱的影響，為緊湊型換熱器的設(shè)計提供有效的數(shù)據(jù)支持。同時將實驗結(jié)果與文獻(xiàn)關(guān)聯(lián)式進(jìn)行比較，找出預(yù)測效果較好的關(guān)聯(lián)式類型。

1 實驗裝置與原理

1.1 實驗裝置

以R134a為實驗工質(zhì)，R134a管內(nèi)流動沸騰換熱實驗系統(tǒng)由兩個循環(huán)系統(tǒng)(實驗工質(zhì)循環(huán)系統(tǒng)和水系統(tǒng))和一套測控系統(tǒng)組成。實驗系統(tǒng)原理圖如圖1所示。

實驗工質(zhì)循環(huán)系統(tǒng)由柱塞計量泵、質(zhì)量流量計、預(yù)熱器、實驗段、節(jié)流閥、冷凝器、儲液器和過冷器組成。實驗時，R134a液體由柱塞計量泵壓出，流過質(zhì)量流量計，在預(yù)熱器中被加熱，再經(jīng)過實驗段，然后在冷凝器中被冷凝成飽和液體，流進(jìn)過冷器，最后進(jìn)入柱塞計量泵完成一個循環(huán)。預(yù)熱器對R134a進(jìn)行預(yù)熱以達(dá)到實驗所要求的干度，根據(jù)預(yù)熱器入口溫度、壓力傳感器分別計算出其對應(yīng)焓值與實驗段入口制冷劑焓值之差，再根據(jù)制冷劑流量可計算預(yù)熱器進(jìn)出口加熱量，調(diào)節(jié)預(yù)熱器加熱量從而可實現(xiàn)實驗段入口干度的控制，由此來看，預(yù)熱段前工質(zhì)的溫度是實驗段入口處干度計算的重要參數(shù)。冷凝器可調(diào)節(jié)R134a的冷凝溫度，以維持實驗段的壓力，過冷器使進(jìn)入柱塞計量泵的R134a保持過冷，防止R134a通過柱塞計量泵的進(jìn)出口閥組時產(chǎn)生閃發(fā)蒸氣。

1.2 測量方法

實驗段水平放置在預(yù)熱器之后，節(jié)流閥之前，實驗段的進(jìn)出口分別裝有相同內(nèi)徑的石英玻璃管，可利用高速攝像機記錄進(jìn)出測試管的流型。實驗過程中，實驗段兩端固定電極與直流穩(wěn)壓電源連接，通過低電壓、大電流的方法直接加熱，這樣實驗段的加熱均勻且調(diào)節(jié)方便。

實驗段外壁溫Tw,o是一個很重要的參數(shù)，分別在實驗段外壁的相隔均勻的12個位置的上、下、左、右四個方向上各布置4個T型熱電偶，對實驗段外壁溫進(jìn)行測量，通過求平均值獲得外壁溫。具體實驗條件見表1，實驗中各測量裝置與參數(shù)的不確定度見表2。

為了整個實驗系統(tǒng)測試的可靠性和穩(wěn)定性，首先對單相對流換熱系數(shù)進(jìn)行了測試，在所有實驗工況下，熱平衡偏差<2.5%。在此熱平衡測試的基礎(chǔ)上，進(jìn)行了兩相流動沸騰換熱實驗。換熱系數(shù)由公式(1)得出。

(1)

式中：Tw,i為實驗段管內(nèi)壁溫，測試過程中制冷劑在圓管換熱可假設(shè)為常物性、具有內(nèi)熱源的一維穩(wěn)態(tài)導(dǎo)熱模型，由下式(2)和式(3)計算得出：

(2)

(3)

式中：Tr為制冷劑溫度，制冷劑進(jìn)入實驗段時為過冷或飽和狀態(tài)，制冷劑從過冷直至飽和時的溫度通過能量平衡方程得出，見式(4)；飽和后制冷劑的溫度通過飽和壓力計算而得，假定L>Lsat時的壓降是線性的，實驗段位置z處飽和壓力和溫度如公式(5)和(6)所示：

(4)

(5)

Tr(z)=Tsat(psat(z))

(6)

當(dāng)某測試工況下系統(tǒng)運行穩(wěn)定后，所有的溫度、壓力和流量測量信號均由數(shù)據(jù)采集儀采集并輸入計算機，利用VB采集軟件處理和存儲。數(shù)據(jù)分析時用到的流體熱物性通過NIST的Refprop v8.0[16]查得。

2 實驗結(jié)果與分析

2.1 熱流密度(q)和質(zhì)量流速(G)對換熱系數(shù)的影響

圖2(a)～(c)分別顯示了飽和溫度10 ℃時熱流密度q和質(zhì)量流速G對3種管徑條件下流體沸騰換熱的影響, 圖中的橫坐標(biāo)和縱坐標(biāo)分別為工質(zhì)的干度和換熱系數(shù)。以圖2(a)所示的3 mm管徑為例，在低干度區(qū)(x<0.6)，當(dāng)質(zhì)量流速G=300 kg/(m2·s)，熱流密度q=24 kW/m2條件下的換熱系數(shù)要高于q=15 kW/m2，而當(dāng)熱流密度q保持在15 kW/m2，質(zhì)量流速G=200 kg/(m2·s)條件下的換熱系數(shù)與G=300 kg/(m2·s)時相近；相比之下，在高干度區(qū)(x>0.6)，當(dāng)熱流密度q保持在15 kW/m2，質(zhì)量流速G=300 kg/(m2·s)條件下的換熱系數(shù)要高于G=200 kg/(m2·s)，而當(dāng)質(zhì)量流速G固定為300 kg/(m2·s)，熱流密度q=24 kW/m2條件下的換熱系數(shù)與q=15 kW/m2時相近。類似的規(guī)律在其他曲線也比較明顯，這些實驗結(jié)果表明：在低干度區(qū)R134a流動沸騰換熱系數(shù)隨熱流密度的增大而增大，與質(zhì)量流速的關(guān)系很小，可推斷此時的流動沸騰換熱由核態(tài)沸騰換熱占主導(dǎo)；在高干度區(qū)R134a流動沸騰換熱系數(shù)隨質(zhì)量流速的增大而顯著增大，而與熱流密度的關(guān)系很小，可推斷此時的流動沸騰換熱由強制對流換熱占主導(dǎo)。3種管徑對比來看，隨著管徑的減小，在低干度區(qū)熱流密度對R134a流動沸騰換熱呈明顯的正相關(guān)關(guān)系，單獨的熱流密度q和質(zhì)量流速G的增大對換熱的加強都是有利的，且當(dāng)管徑減小時，對流沸騰換熱的主導(dǎo)作用在削弱，核態(tài)沸騰換熱作用依舊顯著，尤其是1 mm管徑下熱流密度對換熱的提升十分巨大，說明此時核態(tài)沸騰占據(jù)絕對的主導(dǎo)地位。

總體來看：相同條件下干涸前1 mm管徑換熱系數(shù)要高出2 mm管徑和3 mm管徑，2 mm管徑較3 mm管換熱平均增幅較小為11.6%，1 mm管較2 mm管徑換熱平均增幅較大為26.3%，1 mm管比3 mm管平均增加40.8%；隨著干度的增加，3 mm管徑在x<0.5區(qū)域換熱不隨干度而改變，2 mm管徑這個區(qū)域變小至x<0.4，而1 mm管徑的換熱系數(shù)在整個干度區(qū)域內(nèi)已消失平緩區(qū)，整個換熱曲線全部呈上升趨勢，且明顯陡于2 mm和3 mm管徑下的換熱曲線，說明一定范圍內(nèi)管徑的減小確實有利于換熱系數(shù)的提升。另外管徑減小時，換熱系數(shù)hexp的減小發(fā)生在更低的干度條件下，3 mm管在干度x=0.9處換熱驟減，2 mm管和1 mm管則分別在x=0.82和x=0.6處，即管徑的減小，導(dǎo)致表面張力成為主導(dǎo)力，此時微尺度效應(yīng)顯現(xiàn)，干涸在更低的干度下發(fā)生。

2.2 流型對換熱系數(shù)的影響

換熱系數(shù)的變化與流型的轉(zhuǎn)變密不可分，從流型的角度研究換熱以及干涸，可以更好地理解換熱規(guī)律和干涸現(xiàn)象[4,17]。

隨著管道直徑逐漸減小至微細(xì)尺度，管道內(nèi)的流型與常規(guī)尺度的流型存在很大的不同，這是由于管道內(nèi)的重力和慣性力的作用慢慢變小，取而代之的是表面張力和剪切力的作用逐漸變大。Revellin R & Thome J R[18]對微尺度流型研究表明：在常規(guī)大尺度管中的層流或者分層波狀流受到抑制，彈狀流和柱塞流的聚合形成了長氣泡流型，長氣泡流態(tài)在較低的干度下轉(zhuǎn)化成混狀環(huán)流或環(huán)狀流。

本實驗通過德國MIKROTRONCube4高速攝像機實現(xiàn)可視化研究，觀察到泡狀流、塞狀流、彈狀流和環(huán)狀流，而常規(guī)管徑下的光滑分層流和波狀分層流現(xiàn)象在微尺度管徑的流型研究實驗中都沒有觀察到。在微尺度管徑的流型實驗中觀察到的流型轉(zhuǎn)變的界線非常明顯，易于捕捉和記錄。下圖3(a)～(c)分別顯示了飽和溫度為30 ℃時1～3 mm管徑下R134a的流型圖以及轉(zhuǎn)變界線，流型圖的橫坐標(biāo)和縱坐標(biāo)分別為流體的氣相速度和液相速度，圖中的粗實線為流型轉(zhuǎn)變的邊界線，虛線為一定質(zhì)量流速下流型隨著干度的變化曲線。

圖3(a)是管徑3 mm下的流型圖，從圖中可以看出：在質(zhì)量流速為300 kg/(m2·s)和600 kg/(m2·s)時，隨著干度的增大，流型由塞狀流轉(zhuǎn)變?yōu)閺棤盍?，再轉(zhuǎn)變?yōu)榄h(huán)狀流；質(zhì)量流速為1500 kg/(m2·s)時，流型保持為泡狀流；從各個圖之間的對比來看，隨著管徑的減小，泡狀流、塞狀流和彈狀流的區(qū)域減小，環(huán)狀流的區(qū)域增大，尤其是1 mm和3 mm對比，泡狀流和彈狀流的差距十分明顯；一定質(zhì)量流速下，如當(dāng)G=300 kg/(m2·s)時，塞狀流和彈狀流的臨界干度、彈狀流和環(huán)狀流的臨界干度從3 mm管徑下的0.098、0.4下降為1 mm管徑下的0.084、0.15，即一定質(zhì)量流速時，隨著管徑的下降，塞狀流和彈狀流的區(qū)域范圍減小，環(huán)狀流區(qū)域增大；隨著管徑的減小，泡狀流的區(qū)域也在減小，在實驗條件質(zhì)量流速G=1500 kg/(m2·s)以內(nèi)，3 mm管徑下G≥1000 kg/(m2·s)時會出現(xiàn)泡狀流，1 mm管徑下G≥1350 kg/(m2·s)時才會出現(xiàn)泡狀流，即隨著管徑的下降，管內(nèi)產(chǎn)生泡狀流的臨界質(zhì)量流速G增大。從以上規(guī)律可以看出：在一定質(zhì)量流速范圍內(nèi)，微尺度條件下管徑的下降縮小了泡狀流的區(qū)域，致使產(chǎn)生泡狀流的臨界質(zhì)量流速增大，說明縮小管徑有利于氣泡的成長、變大，有利于氣泡從孤立小氣泡轉(zhuǎn)變?yōu)榫酆洗髿馀?，對塞狀流和彈狀流的產(chǎn)生有促進(jìn)作用；同時，管徑的減小也縮小了塞狀流和彈狀流的區(qū)域，促使核態(tài)沸騰的主導(dǎo)機制在更低的干度結(jié)束，環(huán)狀流區(qū)域的增大，促使對流沸騰的主導(dǎo)機制在更低干度開始；這樣的結(jié)果是隨著管徑的減小，換熱系數(shù)在干涸點以前有所增大，但同時干涸點的提前造成的換熱系數(shù)下降對管內(nèi)流動沸騰換熱存在不利影響。

2.3 與關(guān)聯(lián)式比較

根據(jù)對核態(tài)沸騰傳熱和強制對流傳熱兩種傳熱機理的合成方式可將關(guān)聯(lián)式分為三種類型：加和模型、增強模型和漸進(jìn)模型。從適用工質(zhì)來看，可以分為通用關(guān)聯(lián)式和專用關(guān)聯(lián)式，鑒于通用關(guān)聯(lián)式對于R134a的換熱預(yù)測普遍不如專用關(guān)聯(lián)式，本文選取了R134a專用關(guān)聯(lián)式：Saitoh公式[19]和Fang公式[20]，旨在找出預(yù)測R134a換熱精度較高的關(guān)聯(lián)式，為換熱器的設(shè)計運算提供數(shù)據(jù)基礎(chǔ)。Saitoh公式在原有Chen公式[21]基礎(chǔ)上考慮到管徑的影響，通過加入韋伯?dāng)?shù)對公式進(jìn)行了修正；Fang公式以19篇文獻(xiàn)的2286組R134a數(shù)據(jù)為基礎(chǔ)，運用了與流動沸騰換熱緊密相關(guān)的無量綱數(shù)Fa，并體現(xiàn)出較高的精確度。

Saitoh公式如下：

hTP=Fhl+Shn

(7)

(8)

(9)

(10)

(11)

(12)

ReTP=RelF1.25

(13)

Fang公式如下：

(14)

(15)

(16)

(17)

(18)

圖4和圖5分別顯示了實驗值與Saitoh關(guān)聯(lián)式和Fang關(guān)聯(lián)式的對比情況，圖中的橫坐標(biāo)為換熱實驗值hexp，縱坐標(biāo)為關(guān)聯(lián)式計算值hcal。從對比可以看出，Saitoh公式對3 mm管的預(yù)測效果較好，90.9%點都處于±30%以內(nèi)，干涸區(qū)域有部分偏差；對于2 mm管和1 mm管的預(yù)測精度要差一些，尤其在干度小于0.2和大于0.8時誤差最大，2 mm管的實驗值偏大，而1 mm管的實驗值偏小。Fang公式對三種管徑的換熱預(yù)測都要好于Saitoh公式，尤其是3 mm管和2 mm管的預(yù)測精度很高，而1 mm管的較多實驗值偏小，盡管在干度很小和很大的情況誤差也增大，但也好于Saitoh公式。共同點方面來看，兩個公式與實驗值的趨勢是吻合的，都是換熱隨干度先增大、后減小，F(xiàn)ang公式的跟隨性更好?；谝陨蠈Ρ确治觯梢酝茖?dǎo)流動沸騰換熱在低干度區(qū)為核態(tài)沸騰主導(dǎo)，此時換熱與干度無關(guān)，隨著干度增加至干涸或部分干涸前，核態(tài)沸騰和對流換熱共同影響換熱，對流換熱作用增強，當(dāng)換熱到達(dá)峰值后，出現(xiàn)部分干涸致使換熱系數(shù)下降。雖然以上對比是在飽和溫度為10 ℃的情況下，但其他工況下數(shù)據(jù)的規(guī)律和適用性是一致的，F(xiàn)ang公式無論是整體的精確度還是單一管徑的適用度都要好于Saitoh公式。

3 結(jié)論

通過實驗對內(nèi)徑分別為1 mm、2 mm和3 mm的3種水平不銹鋼光滑管內(nèi)R134a的流動沸騰換熱進(jìn)行了研究，重點在于研究管徑對于流動沸騰換熱的影響，實驗的飽和溫度為5～30 ℃，熱流密度為2～70 kW/m2，流量范圍為200～1500 kg/(m2·s)。通過實驗得到以下結(jié)論：

1)3種管徑下，單獨的熱流密度q和質(zhì)量流速G的增大對換熱的加強都是有利的，當(dāng)管徑減小時，對流沸騰換熱的主導(dǎo)作用在削弱，核態(tài)沸騰換熱作用依舊顯著。

2)一定條件下，管徑的減小確實有利于換熱系數(shù)的提升，但同時導(dǎo)致表面張力成為主導(dǎo)力，此時微尺度效應(yīng)顯現(xiàn)，干涸在更低的干度下發(fā)生。

3)流型與換熱系數(shù)的變化密不可分，隨著管道直徑逐漸減小至微細(xì)尺度，常規(guī)通道下的光滑分層流和波狀分層流都沒有觀察到。小管徑有利于氣泡的成長、變大，有利于氣泡從孤立小氣泡轉(zhuǎn)變?yōu)榫酆洗髿馀?，對塞狀流和彈狀流的產(chǎn)生有促進(jìn)作用；管徑的減小也促使核態(tài)沸騰的主導(dǎo)機制在更低的干度結(jié)束，流沸騰的主導(dǎo)機制在更低干度開始；最終的結(jié)果是隨著管徑的減小，換熱系數(shù)在干涸點以前有所增大，但同時干涸點的提前造成的換熱系數(shù)下降對管內(nèi)流動沸騰換熱存在不利影響。

4)Saitoh公式對3 mm管的換熱預(yù)測效果較好，對于2 mm管和1 mm管的預(yù)測精度要差一些；Fang公式對三種管徑的換熱預(yù)測都要好于Saitoh公式，尤其是3 mm管和2 mm管的預(yù)測精度很高。

[1] Lazarek G M, Black S H. Evaporative heat transfer, pressure drop and critical heat flux in small vertical tube with R113[J]. International Journal of Heat and Mass Transfer, 1982, 25(7): 945-960.

[2] Qu W, Mudawar I. Flow boiling heat transfer in two-phase micro-channel heat sinks—I. Experimental investigation and assessment of correlation methods[J]. International Journal of Heat and Mass Transfer, 2003, 46(15): 2755-2771.

[3] Jong-Taek Oh A S，Pamitran, Kwang-Il Choi, et al. Experimental investigation on two-phase flow boiling heat transfer of five refrigerants in horizontal small tubes of 0.5, 1.5 and 3.0 mm inner diameters[J]. International Journal of Heat and Mass Transfer, 2011, 54(9/10): 2080-2088.

[4] Ong C L, Thome J R. Macro-to-microchannel transition in two-phase flow: Part 1—Two-phase flow patterns and film thickness measurements[J]. Experimental Thermal and Fluid Science, 2011, 35(1): 37-47.

[5] Ong C L, Thome J R. Macro-to-microchannel transition in two-phase flow: Part 2— Flow boiling heat transfer and critical heat flux[J]. Experimental Thermal and Fluid Science, 2011, 35(6): 873-886.

[6] Bao Z Y, Fletcher D F, Haynes B S. Flow boiling heat transfer of Freon R11 and HCFC123 in narrow passages[J]. International Journal of Heat and Mass Transfer, 2000, 43(18): 3347-3358.

[7] Kuwahara K, Koyama S, Hashimoto Y. Characteristics of evaporation heat transfer and flow pattern of pure refrigerant HFC134a in a horizontal capillary tube[C]//R134a flow boiling heat transfer in small diameter tubes: Proceedings of the Fourth JSME-KSME Thermal Engineering Conference. 2000: 385-390.

[8] Lin S, Kew P A, Cornwell K. Two-phase heat transfer to a refrigerant in a 1 mm tube[J]. International Journal of Refrigeration, 2001, 24(1): 51-56.

[9] Lin S, Kew P A, Cornwell K. Flow boiling of refrigerant R141b in small tubes [J]. Chemical Engineering Research and Design Chemical Engineering Research and Design, 2001, 79(4): 417-424.

[10] Yan Y Y, Lin T F. Evaporation heat transfer and pressure drop of refrigerant R134a in a small pipe[J]. International Journal of Heat and Mass Transfer, 1998, 41: 4183-4194.

[11] Consolini L, Thome J R. A heat transfer model for evaporation of coalescing bubbles in micro-channel flow[J]. International Journal of Heat and Fluid Flow, 2010, 31(1): 115-125.

[12] Sira Saisorn, Jatuporn Kaew-On, Somchai Wongwises. Two-phase flow of R-134a refrigerant during flow boiling through a horizontal circular mini-channel[J]. Experimental Thermal and Fluid Science, 2011, 35(6): 887-895.

[13] Grauso S, Mastrullo R, Mauro A, et al. Flow pattern map, heat transfer and pressure drops during evaporation of R-1234ze(E) and R134a in a horizontal, circular smooth tube: Experiments and assessment of predictive methods[J]. International Journal of Refrigeration, 2013, 36(2): 478-491.

[14] Cristiano Bigonha Tibiri?á, Ribatski G, Thome J R. Saturated flow boiling heat transfer and critical heat flux in small horizontal flattened tubes[J]. International Journal of Heat and Mass Transfer, 2012, 55(25/26): 7873-7883.

[15] Huo X. Experimental study of boiling heat transfer in small diameter tube[D]. UK: London South Bank University, 2006.

[16] Lemmon E W, Huber M L, McLinden M O. NIST Standard Reference Databases 23, Version 8.0 [DB/CD]. Gaithersburg: MD, USA, US Department of Commerce, Technology Administration, National Institute of Standards and Technology, 2007.

[17] Charnay R, Bonjour J, Rémi Revellin. Experimental investigation of R-245fa flow boiling in minichannels at high saturation temperatures: Flow patterns and flow pattern maps[J]. International Journal of Heat and Fluid Flow, 2014, 46:1-16.

[18] Revellin R, Thome J R. A new type of diabatic flow pattern map for boiling heat transfer in microchannels[J].Journal of Micromechanics and Microengineering, 2007, 17(4): 788-796.

[19] Saitoh S, Daiguji H, Hihara E. Correlation for boiling heat transfer of R-134a in horizontal tubes including effect of tube diameter[J]. International Journal of Heat and Mass Transfer, 2007, 50(25/26): 5215-5225.

[20] Fang X D. A new correlation of flow boiling heat transfer coefficients based on R134a data[J]. International Journal of Heat and Mass Transfer, 2013, 66: 279-283.

[21] Chen J C. Correlation for boiling heat transfer to saturated fluids inconvective flow[J]. Ind. Eng. Chem. Proc. Des. Dev., 1966,5(3): 322-329.

About the author

Ding Yang, male, master candidate, School of Energy and Power Engineering, University of Shanghai for Science and Technology,+86 21-55271153,E-mail: dingyang19900102@163.com.Research fields: refrigeration testing technology and automation.

Experimental Studies on Flow Boiling of R134a in Horizontal Small Tubes

Ding Yang1Liu Jianhua1Ye Fangping2Jiang Linlin1E Xiaoxue1Wu Hao1

(1.School of Energy and Power Engineering, University of Shanghai for Science and Technology, Shanghai, 200093, China; 2.Zhejiang Xinjing Air Conditioning Equipment Co., Ltd., Longquan, 323700, China)

Experimental investigation on the flow boiling of R134a in horizontal mini-channels was carried out in this study. The test sections are three different horizontal smooth stainless steel tubes with inner diameter of 1mm, 2mm and 3mm. The experiment was performed on the saturated temperatures ranging from 5 to 30 ℃, heat fluxes from 2 to 70 kW/m2and mass fluxes from 200 to 1500 kg/(m25s). The results indicated that on the same conditions the pre-dryout heat transfer coefficients in 2-mm-ID tube increased 11.6% averagely than the 3-mm-ID tube, the coefficients in 1-mm-ID tube increased 26.3% averagely than the 2-mm-ID tube and the heat transfer coefficients in the 1-mm-ID tube were 40.8% greater in average than those in the 3-mm-ID tube. With decreasing tube diameter the coefficients started decreasing at lower vapor quality, the effect of mass flux and the contribution of forced convective evaporation on the local heat transfer coefficient decreased while the effect of heat flux was still strong in all three tubes, the zone of plug flow and slug flow decreased while the zone of bubbly flow and annular flow increased.

R134a; mini-channel; flow boiling; dryout; heat transfer coefficient

0253- 4339(2015) 01- 0090- 08

10.3969/j.issn.0253- 4339.2015.01.090

2014年3月17日

TB61+2; TQ051.5; TB61+1

A

丁楊，男，碩士研究生，上海理工大學(xué)能源與動力工程學(xué)院，(021)55271153，E-mail：dingyang19900102@163.com。研究方向：制冷測試技術(shù)與自動化。