超臨界CO2/R41小通道內(nèi)的換熱特性

代寶民，李敏霞，呂佳桐，王派，馬一太

?

超臨界CO2/R41小通道內(nèi)的換熱特性

代寶民，李敏霞，呂佳桐，王派，馬一太

（天津大學機械工程學院，中低溫熱能高效利用教育部重點實驗室，天津 300072）

對R41和混合工質(zhì)CO2/R41 (20.5/79.5)、CO2/R41（51.4/48.6）在直徑為2 mm的水平光滑圓管中的超臨界冷卻流動換熱特性進行了實驗研究。質(zhì)量流速范圍為400～800 kg·m-2·s-1，壓力為6.0～8.0 MPa，熱通量為12～48 kW·m-2，流體溫度為20～80℃。3種工質(zhì)的對流傳熱系數(shù)的極值隨CO2含量的增加而增大，在相同條件下R41的傳熱系數(shù)小于CO2/R41的傳熱系數(shù)?；旌衔锏某R界傳熱系數(shù)變化規(guī)律與純R41相同。實驗條件下，3種流體的傳熱系數(shù)在2～25 kW·m-2·K-1之間，壓力的影響顯著，越接近臨界壓力對應壓力條件下的傳熱系數(shù)極值越高。在遠離準臨界點的區(qū)域傳熱系數(shù)隨熱通量變化不明顯，而在準臨界點附近對流傳熱系數(shù)的極值隨熱通量的增加而小幅減小。將實驗結(jié)果與經(jīng)驗關(guān)聯(lián)式計算結(jié)果進行了比較，有4個關(guān)聯(lián)式的預測效果較好，誤差均在±30%以內(nèi)，預測誤差隨CO2含量的增加而下降。

超臨界CO2；R41；混合物；對流；換熱；小通道

引 言

歐盟“禁氟令”規(guī)定：自2015年1月1日起，GWP（全球變暖潛能值）大于150的制冷劑停止在家用制冷和冷凍系統(tǒng)使用[1]。對于耗能較高的制冷空調(diào)領(lǐng)域，傳統(tǒng)制冷劑（HFC，HCFC）由于其較高的GWP，使用量會逐步消減，直到完全停止使用。為了尋求新型節(jié)能環(huán)保的制冷劑，以下幾種途徑在制冷領(lǐng)域受到了廣泛關(guān)注：①低GWP制冷劑；②自然工質(zhì)；③混合制冷劑[2]。

以R1234yf和R1234ze為代表的新一代制冷劑雖然GWP較低（分別為4和6），但其合成路線復雜，產(chǎn)量較低，造成現(xiàn)階段價格較昂貴[3]；自然工質(zhì)如丙烷、丁烷等烴類化合物（HC）由于其較高的可燃性限制了其在大型設(shè)備中的應用；氨由于較高的毒性以及近年來事故頻發(fā)，其安全使用也引起制冷界的廣泛討論；CO2由于其優(yōu)良的熱物性及傳輸特性、GWP1、無毒不可燃、易獲取等優(yōu)點，被認為是最具潛力的制冷劑，但其運行壓力較高、能效較低的缺點也限制了其廣泛推廣應用。Kim等[4-5]、Niu等[6]、Zhang等[7]、Sarkar等[8]、Onaka等[9]和Hakkaki-Fard等[10]的理論及實驗研究表明，將R290、R600、R600a、RE170和R32等制冷劑與CO2按一定比例混合后用于熱泵及制冷系統(tǒng)，可以顯著提升系統(tǒng)的COP（制熱能效比），并且降低系統(tǒng)的運行高壓。將CO2與普通工質(zhì)按一定比例混合，既可以降低CO2較高的運行壓力，也可以通過CO2的阻燃性提高可燃低GWP工質(zhì)的安全性，并且混合后的制冷劑依然保持較低的GWP，滿足環(huán)保的要求。

R41的物理性質(zhì)與CO2相近，臨界溫度較低（表1），僅為44.1℃，ODP為0，GWP為107，屬低GWP工質(zhì)，但有可燃性。Dai 等[11]對CO2與10種低GWP工質(zhì)的混合物用于熱泵熱水器進行了分析和篩選，結(jié)果表明R41與CO2混合后系統(tǒng)的性能最優(yōu)，在最優(yōu)配比條件下系統(tǒng)COP提高了4%，并大幅降低了運行壓力。因此，將R41與CO2混合用于熱泵熱水器是改善CO2系統(tǒng)性能的有效方法。之前的分析[11]表明，將R41與CO2以任意比例混合，其放熱過程均發(fā)生在超臨界狀態(tài)。

表1 CO2及R41的物理性質(zhì)、安全及環(huán)保特性

CO2流體在放熱過程中壓力較高（8～12 MPa），常規(guī)換熱器的承壓能力有限、體積較大，氣冷器的運行和制造成本較高。而微小通道換熱器承壓能力強、結(jié)構(gòu)緊湊、換熱高效，并可減小工質(zhì)充灌量，更加適用于類似CO2等流體的高壓運行工況。然而到目前為止R41及其與CO2混合物的換熱數(shù)據(jù)非常少，本文對純R41以及不同混合比的超臨界CO2/R41流體在水力直徑為2 mm的水平光滑圓管中的換熱特性進行了實驗研究，為CO2/R41制冷及熱泵系統(tǒng)的氣冷器設(shè)計提供理論依據(jù)。

1 實驗系統(tǒng)和數(shù)據(jù)處理

1.1 實驗系統(tǒng)介紹

實驗系統(tǒng)及實驗段如圖1所示。由圖1 (a)可以看到，工質(zhì)由齒輪泵推動，通過流量計后進入恒溫水浴加熱，之后進入電加熱段繼續(xù)加熱，通過控制電加熱量獲得不同溫度的流體，加熱量用測功儀測量。工質(zhì)加熱至預定溫度后進入實驗段，與換熱流體換熱后進入冷凝器，冷卻至過冷液，流過過濾器進入泵入口，完成制冷劑側(cè)循環(huán)。與工質(zhì)進行換熱的流體為水，同樣由泵推動，流過流量計后進入實驗段，與工質(zhì)換熱后進入恒溫水浴冷卻，完成換熱流體側(cè)循環(huán)。工質(zhì)側(cè)及水側(cè)的流量均用Coriolis質(zhì)量流量計測定。圖1 (b)為實驗段的軸向截面圖，為水平放置同軸套管換熱段，工質(zhì)走管內(nèi)側(cè)，水走管外的環(huán)狀通道，兩側(cè)流體逆向流動。與水換熱的有效長度為480 mm。12根線徑為0.13 mm的T型熱電偶焊接在6個截面上，每個截面上下外壁面各布置一根熱電偶以測量壁面溫度。用鎧裝鉑電阻測量工質(zhì)側(cè)和水側(cè)進出口以及電加熱入口的溫度。用壓力傳感器測量工質(zhì)進出口及電加熱入口的壓力。實驗段的截面圖及尺寸如圖1 (c)所示。測試物理量的不確定度見表2，可以看到測量誤差滿足實驗精度的要求。測試的工質(zhì)包括組分質(zhì)量比為20.5/79.5和51.4/58.6的兩種CO2/R41混合物以及純R41?；旌衔锏慕M分通過文獻[12]中的方法進行測定。

圖1 實驗系統(tǒng)及實驗段

表2 不確定度分析結(jié)果

1.2 實驗數(shù)據(jù)處理

冷卻的熱通量通過水側(cè)換熱量與實際換熱面積確定

式中，water為水側(cè)換熱量，通過式(2)計算

超臨界流體在管內(nèi)換熱過程中溫度發(fā)生改變，定義其整體溫度b為

工質(zhì)內(nèi)壁面溫度通過測得的外壁面溫度推算

假設(shè)工質(zhì)局部溫度b()沿實驗段工質(zhì)流動方向線性分布，工質(zhì)側(cè)的局部對流傳熱系數(shù)為

平均對流傳熱系數(shù)為

相對平均偏差（mean absolute error，MAE）通過式(7)計算

83.6%的實驗數(shù)據(jù)的熱平衡在±10%以內(nèi)，其中熱平衡較差的數(shù)據(jù)出現(xiàn)在準臨界點（在某一特定超臨界壓力條件下比定壓熱容在某一溫度處存在最大值，這一狀態(tài)點定義為準臨界點，對應的溫度為準臨界溫度）附近，這是因為在準臨界點附近工質(zhì)的進出口溫差較小，并且焓值在準臨界點附近隨溫度變化劇烈，導致工質(zhì)側(cè)換熱量的測量誤差較大。工質(zhì)及換熱流體的物性均通過REFPROP 9.0[13]計算。

2 結(jié)果及討論

2.1 壓力的影響

圖2為純R41、CO2/R41（20.5/79.5）和CO2/R41（51.4/48.6）3種工質(zhì)在壓力為6.0～8.0 MPa時對流傳熱系數(shù)隨整體溫度的變化關(guān)系。熱通量恒定為24 kW·m-2，圖2 (a)和圖2 (b)的質(zhì)量流速分別恒定為400 kg·m-2·s-1和800 kg·m-2·s-1。混合物的傳熱系數(shù)變化規(guī)律與純工質(zhì)相同?？梢钥吹綄α鱾鳠嵯禂?shù)隨溫度升高先逐漸增大，在稍高于準臨界溫度處達到極大值，圖中的虛線為對應壓力下的準臨界溫度；而后隨著溫度的進一步升高，對流傳熱系數(shù)逐漸減小。在準臨界溫度附近對流傳熱系數(shù)急劇增大，但在遠離準臨界溫度時對流傳熱系數(shù)隨溫度的變化并不明顯。

圖2 壓力對傳熱系數(shù)的影響

對于常規(guī)的管內(nèi)對流換熱，傳熱系數(shù)隨溫度的變化不明顯，這是因為流體的物性隨溫度變化較小，但在超臨界壓力條件下流體的物性隨溫度變化劇烈，尤其是在準臨界點附近。Kim等[14]的分析結(jié)果表明，對于超臨界流體，邊界層的比定壓熱容對對流傳熱系數(shù)影響顯著，較大的比熱容對應較大的對流傳熱系數(shù)。當邊界層的溫度等于準臨界溫度時，對流傳熱系數(shù)取得極大值，此時整體溫度稍高于準臨界溫度。因此，對流傳熱系數(shù)出現(xiàn)在稍高于準臨界溫度處。從圖2還可以看到，壓力越靠近臨界壓力，對應壓力下的傳熱系數(shù)的極值越高，同樣是受比定壓熱容的影響。圖3為不同工質(zhì)的比定壓熱容隨整體溫度的變化規(guī)律，可以看到對流傳熱系數(shù)的變化趨勢與比定壓熱容隨溫度的變化一致。

圖3 不同混合物組分對應的比定壓熱容

2.2 質(zhì)量流速的影響

圖4為不同質(zhì)量流速下3種工質(zhì)的對流傳熱系數(shù)。對應的壓力為7.5 MPa，熱通量恒定為24 kW·m-2?？梢钥吹?，對于3種工質(zhì)，對流傳熱系數(shù)均隨質(zhì)量流速的增加而增大，并且均在稍高于準臨界溫度時取得最大值，不受質(zhì)量流速以及工質(zhì)濃度配比的影響。質(zhì)量流速為800 kg·m-2·s-1與 400 kg·m-2·s-1相比，質(zhì)量流速增大了1倍，對流傳熱系數(shù)提高了1.9倍左右。這是由于實驗中的工況Reynolds數(shù)（）均大于2300，為湍流流動。隨著管內(nèi)的質(zhì)量流速的增加，湍流強度增強，換熱效果明顯提升。

圖4 質(zhì)量流速對傳熱系數(shù)的影響

2.3 熱通量的影響

圖5為3種工質(zhì)在壓力為7.0 MPa、質(zhì)量流速為800 kg·m-2·s-1時不同熱通量下對應的對流傳熱系數(shù)。在準臨界溫度附近，對流傳熱系數(shù)隨熱通量的增加而明顯減小，尤其在對流傳熱系數(shù)的極值處。這是因為在準臨界點附近區(qū)域比定壓熱容隨溫度的變化劇烈，熱通量越大，工質(zhì)進出口溫差越大，導致在準臨界點處的平均比熱容越小。由2.1節(jié)的分析可知，對流傳熱系數(shù)與比定壓熱容的變化趨勢一致，因此平均比熱容的減小導致對流傳熱系數(shù)的極值隨熱通量的增加而減小。但在遠離準臨界點的區(qū)域，對流傳熱系數(shù)隨熱通量的變化不明顯。這是因為在遠離臨界點的區(qū)域流體更接近常規(guī)流體，其物性隨溫度的變化幅度較小，流體的物性對熱通量的變化不敏感。

圖5 熱通量對傳熱系數(shù)的影響

2.4 混合物組分的影響

由圖2、圖4和圖5可以看到，對于純R41、CO2/R41（20.5/79.5）和CO2/R41（51.4/48.6）3種工質(zhì)，在壓力、質(zhì)量流速和熱通量均相同的工況下，CO2/R41（51.4/48.6）的對流傳熱系數(shù)極值最高，純R41的最低?；旌衔锏呐R界壓力隨CO2含量的增加而升高，3種工質(zhì)的臨界壓力分別為5.897、6.247、6.748 MPa。在壓力相同的條件下，混合物中CO2含量越高，越接近準臨界壓力，準臨界點處的定壓比熱容越高。如圖3 (a)和圖3 (b)所示，對于3種工質(zhì)，在壓力分別為7.0 MPa和7.5 MPa的條件下，CO2/R41（51.4/48.6）比定壓熱容的極值最高。由2.1節(jié)的分析可知，比定壓熱容是決定對流傳熱系數(shù)的關(guān)鍵因素，比定壓熱容越高，對流傳熱系數(shù)越大。

由圖2 (a)和圖2 (b)也可以看到，壓力為7.0 MPa時，不同混合比的工質(zhì)對應的對流傳熱系數(shù)差異明顯，CO2/R41 (51.4/48.6)的對流傳熱系數(shù)的極值比純R41提高了101%；壓力升高至7.5 MPa后，不同工質(zhì)間的差異縮小，CO2/R41 (51.4/48.6)的對流傳熱系數(shù)的極值比純R41僅提高了56%，在數(shù)量上同圖3所示的兩個壓力條件下的比定壓熱容的變化規(guī)律一致。因此，工作壓力越接近臨界壓力，物性變化得越劇烈，超臨界對流傳熱強度的提升越明顯。

2.5 與換熱關(guān)聯(lián)式比較

為了對實際制冷熱泵系統(tǒng)的氣冷器進行設(shè)計和優(yōu)化，需要精確地預測超臨界流體的對流傳熱系數(shù)。下面選取了12個經(jīng)驗關(guān)聯(lián)式進行對比，其中Gnielinski[15]的關(guān)聯(lián)式適用于2300＜＜106的過冷液或過熱氣的湍流流動換熱，其他的均適用于超臨界流體的對流換熱，其中除了Petrov等[16]的關(guān)聯(lián)式適用于CO2、水及氦以外，其余關(guān)聯(lián)式均基于CO2提出。下面對收集的關(guān)聯(lián)式對CO2/R41的混合物及純R41的超臨界對流傳熱系數(shù)預測的準確性進行驗證。

關(guān)聯(lián)式預測的平均絕對誤差見表3?？梢钥吹?，F(xiàn)ang[17]、Pitla等[18]、Dang等[19]、Fang等[20]的關(guān)聯(lián)式的預測效果最好。以上4個關(guān)聯(lián)式對純R41的預測平均絕對誤差均在±20%以內(nèi)，CO2/R41（20.5/79.5）的在±25%以內(nèi)，CO2/R41（51.4/48.6）的在±30%以內(nèi)。Gnielinski[15]、Petrov等[16]、Kuang等[21]、Petrov等[22]、Liao等[23]的預測結(jié)果偏低，而Yoon等[24]、Huai等[25]、Son等[26]的預測結(jié)果與實驗值偏離較大。其中Son等[26]對3種工質(zhì)的預測結(jié)果均偏高，對CO2/R41 (51.4/48.6)的預測結(jié)果偏離更加明顯；Yoon等[24]的預測結(jié)果對純R41及CO2/R41（20.5/79.5）的預測結(jié)果明顯偏高，但對CO2/R41（51.4/48.6）的結(jié)果偏低，嚴重偏離了±30%的誤差線。

表3 關(guān)聯(lián)式計算的平均絕對誤差

預測較準確的4個關(guān)聯(lián)式的預測誤差如圖6所示?？梢钥吹?，對于R41、CO2/R41（20.5/79.5）和CO2/R41（51.4/48.6）3種工質(zhì)的預測精度基本均隨CO2含量的增加而降低，其中CO2/R41 (51.4/48.6)的預測值較前面兩種更加分散。這是因為，隨著CO2含量的增加，對應的壓力工況越接近臨界壓力，此時物性變化越劇烈，對流傳熱過程中壁溫和流體的溫差較小，導致實驗的測試精度有所降低，造成預測精度有所下降。

圖6 換熱關(guān)聯(lián)式的預測結(jié)果

由以上計算結(jié)果可以看到，雖然兩種工質(zhì)進行了混合，但對于超臨界流動換熱基本不存在相變換熱中的傳質(zhì)阻力對換熱的惡化效應。這是因為超臨界流體呈現(xiàn)一種近似單相流動的狀態(tài)，在流動過程中兩種工質(zhì)充分混合，在傳熱過程中沒有發(fā)生相態(tài)的變化。可以將其近似看作一種特殊的單相流體：在遠離準臨界點的左側(cè)表現(xiàn)為過冷液體，在遠離準臨界點的右側(cè)表現(xiàn)為過熱氣體的特性。雖然混合后為非共沸混合物，但超臨界流動換熱不發(fā)生相變就不存在氣液相之間的濃度差，不會導致傳熱惡化。

3 結(jié) 論

對R41、CO2/R41（20.5/79.5）和CO2/R41（51.4/48.6）在水力直徑為2 mm的水平光滑圓管中的超臨界冷卻流動換熱特性進行了實驗研究。質(zhì)量流速范圍為400～800 kg·m-2·s-1，壓力為6.0～8.0 MPa，熱通量為12～48 kW·m-2，流體溫度為20～80℃。對壓力、質(zhì)量流速、熱通量和混合組分等因素對傳熱系數(shù)的影響進行了實驗研究，得出以下結(jié)論。

（1）3種工質(zhì)超臨界冷卻對流系數(shù)變化的規(guī)律與其他工質(zhì)一樣，在稍高于準臨界溫度處出現(xiàn)極大值。實驗條件下，傳熱系數(shù)范圍為2～25 kW·m-2·K-1，在遠離準臨界點的區(qū)域小幅度變化，越靠近臨界點變化越劇烈。

（2）壓力對對流傳熱系數(shù)的影響明顯，越接近臨界壓力傳熱系數(shù)極值越高。超臨界對流傳熱系數(shù)隨質(zhì)量流速的增加而增大，質(zhì)量流速由400 kg·m-2·s-1增大至800 kg·m-2·s-1，對流傳熱系數(shù)平均提高了1.9倍。

（3）在準臨界點附近，對流傳熱系數(shù)隨熱通量的增加而減??；在遠離準臨界點的區(qū)域，對流傳熱系數(shù)隨熱通量的變化不明顯。在相同的運行壓力下，3種工質(zhì)的對流傳熱系數(shù)的極值隨R41含量增加而減小。

（4）與實驗值比較，發(fā)現(xiàn)Fang[17]、Pitla等[18]、Dang等[19]、Fang等[20]的關(guān)聯(lián)式的預測效果最好，平均絕對誤差均在±30%以內(nèi)，但對3種工質(zhì)的預測誤差隨CO2含量的增加而下降。

References

[1] The European Parliament, the Council of the European Union. Regulation (EU) No 517/2014 of the European Parliament and of the Council on fluorinated greenhouse gases and repealing Regulation (EC) No 842/2006 [J]., 2014, 150: 195-230

[2] Miyara A, Onaka Y, Koyama S. Ways of next generation refrigerants and heat pump/refrigeration systems [J]., 2012, 20 (1): 1130002

[3] McLinden M O, Kazakov A F, Steven Brown J, Domanski P A. A thermodynamic analysis of refrigerants: possibilities and tradeoffs for low-GWP refrigerants [J]., 2014, 38: 80-92

[4] Kim J H, Cho J M, Lee I H, Kim M S. Circulation concentration of CO2/propane mixtures and the effect of their charge on the cooling performance in an air-conditioning system [J]., 2007, 30 (1): 43-49

[5] Kim J H, Cho J M, Kim M S. Cooling performance of several CO2/propane mixtures and glide matching with secondary heat transferfluid [J]., 2008, 31 (5): 800-806

[6] Niu B, Zhang Y. Experimental study of the refrigeration cycle performance for the R744/R290 mixtures [J]., 2007, 30 (1): 37-42

[7] Zhang X, Wang F, Fan X, Wei X, Wang F. Determination of the optimum heat rejection pressure in transcritical cycles working with R744/R290 mixture [J]., 2013, 54 (1): 176-184

[8] Sarkar J, Bhattacharyya S. Assessment of blends of CO2with butane and isobutane as working fluids for heat pump applications [J]., 2009, 48 (7): 1460-1465

[9] Onaka Y, Miyara A, Tsubaki K. Experimental study on evaporation heat transfer of CO2/DME mixture refrigerant in a horizontal smooth tube [J]., 2010, 33 (7): 1277-1291

[10] Hakkaki-Fard A, Aidoun Z, Ouzzane M. Applying refrigerant mixtures with thermal glide in cold climate air-source heat pumps [J]., 2014, 62 (2): 714-722

[11] Dai B, Dang C, Li M, Ma Y. Thermodynamic performance assessment of carbon dioxide blends with low-global warming potential (GWP) working fluids used for a heat pump water heater [J]., 2015, DOI: 10.1016/j.ijrefrig.2014.11.009

[12] Grauso S, Mastrullo R, Mauro A W,. CO2and propane blends: experiments and assessment of predictive methods for flow boiling in horizontal tubes [J]., 2011, 34 (4): 1028-1039

[13] Lemmon E, McLinden M, Huber M. NIST Standard Reference Database 23. Reference Fluid Thermodynamic and Transport Properties Database (REFPROP) [CP]. Gaithersburg, MD; National Institute of Standards and Technology (NIST), 2010

[14] Kim D E, Kim M H. Two layer heat transfer model for supercritical fluid flow in a vertical tube[J]., 2011, 58(1): 15-25

[15] Gnielinski V. New equations for heat and mass transfer in the turbulent flow in pipes and channels [J]., 1975, 16 (2): 359-368

[16] Petrov N, Popov V. Heat transfer and hydraulic resistance with turbulent flow in a tube of water at supercritical parameters of state [J]., 1988, 35 (10): 577-580

[17] Fang X. Modeling and analysis of gas coolers. Air Conditioning and Refrigeration Center [R]. College of Engineering, University of Illinois at Urbana-Champaign, 1999

[18] Pitla S S, Groll E A, Ramadhyani S. New correlation to predict the heat transfer coefficient during in-tube cooling of turbulent supercritical CO2[J]., 2002, 25 (7): 887-895

[19] Dang C, Hihara E. In-tube cooling heat transfer of supercritical carbon dioxide (Ⅰ): Experimental measurement [J]., 2004, 27 (7): 736-747

[20] Fang X, Xu Y. Modified heat transfer equation for in-tube supercriticalCO2cooling [J]., 2011, 31 (14): 3036-3042

[21] Kuang G, Ohadi M, Dessiatoun S. Semi-empirical correlation of gas cooling heat transfer of supercritical carbon dioxide in microchannels [J]., 2008, 14 (6): 861-870

[22] Petrov N, Popov V. Heat transfer and resistance of carbon dioxide being cooled in the supercritical region [J]., 1985, 32 (3): 131-134

[23] Liao S, Zhao T. Measurements of heat transfer coefficients from supercritical carbon dioxide flowing in horizontal mini/micro channels [J]., 2002, 124 (3): 413-420

[24] Yoon S H, Kim J H, Hwang Y W, Kim M S, Min K, Kim Y. Heat transfer and pressure drop characteristics during the in-tube cooling process of carbon dioxide in the supercritical region [J]., 2003, 26 (8): 857-864

[25] Huai X, Koyama S, Zhao T. An experimental study of flow and heat transfer of supercritical carbon dioxide in multi-port mini channels under cooling conditions [J]., 2005, 60 (12): 3337-3345

[26] Son C H, Park S J. An experimental study on heat transfer and pressure drop characteristics of carbon dioxide during gas cooling process in a horizontal tube [J]., 2006, 29 (4): 539-546

Heat transfer characteristics of supercritical CO2/R41 flowing in mini-channel

DAI Baomin, LI Minxia, Lü Jiatong, WANG Pai, MA Yitai

Key Laboratory of Efficient Utilization of Low and Medium Grade EnergyMOESchool of Mechanical EngineeringTianjin UniversityTianjinChina

An experiment was conducted to study the heat transfer characteristics of supercritical pure R41, mixtures of CO2/R41 (20.5/79.5), and CO2/R41 (51.4/48.6) cooled in a horizontal smooth mini-channel tube with inner-diameter of 2 mm. Mass flow rate was in the range of 400 to 800 kg·m-2·s-1, pressure changed from 6.0 to 8.0 MPa，heat flux ranged from 12 to 48 kW·m-2and bulk temperature varied from 20℃to 80℃. The maximum heat transfer coefficient () of the three fluids increased with the CO2mass fraction. Additionally,of pure R41 is smaller than that of the mixtures. However,variation tendency of the mixtures is the same as that of pure R41.of the three working fluids are ranging from 2 kW·m-2·K-1to 25 kW·m-2·K-1under the present test condition. The influence of pressure is significant. The closer of the pressure approaches the critical pressure, the higher of the maximumis. The heat flux has little influence onwhen the bulk temperature is far away from the critical temperature. However for the bulk temperature near the critical temperature,increases with the decrease of heat flux. The experimental results are compared with prediction values calculated by twelve correlations. It is concluded that four correlations predicting with good precisions, and the errors are within ±30%. Nevertheless, the prediction errors increase with CO2mass fraction.

supercritical CO2; R41; mixture; convection; heat transfer; mini-channel

2014-09-18.

LI Minxia, tjmxli@tju.edu.cn

10.11949/j.issn.0438-1157.20141404

TK 124

A

0438—1157（2015）03—0924—08

國家自然科學基金項目（50976075）；天津市科委基金項目（12JCYBJC13800）。

2014-09-18收到初稿，2014-11-19收到修改稿。

聯(lián)系人：李敏霞。第一作者：代寶民（1987—），男，博士研究生。

supported by the National Natural Science Foundation of China (50976075) and the Tianjin Municipal Science and Technology Commission Foundation (12JCYBJC13800)