CO2冷風(fēng)機(jī)換熱性能仿真及實(shí)驗(yàn)研究

申江 邊煜竣 黃冰

(天津商業(yè)大學(xué)天津市制冷技術(shù)重點(diǎn)實(shí)驗(yàn)室 天津 300134)

CO2冷風(fēng)機(jī)換熱性能仿真及實(shí)驗(yàn)研究

申江 邊煜竣 黃冰

(天津商業(yè)大學(xué)天津市制冷技術(shù)重點(diǎn)實(shí)驗(yàn)室 天津 300134)

本文利用穩(wěn)態(tài)分布參數(shù)法對(duì)冷風(fēng)機(jī)建立仿真模型，并利用冷風(fēng)機(jī)性能實(shí)驗(yàn)臺(tái)對(duì)冷風(fēng)機(jī)樣機(jī)進(jìn)行實(shí)驗(yàn)研究，利用實(shí)驗(yàn)研究與數(shù)值模擬相結(jié)合的方法，對(duì)冷風(fēng)機(jī)換熱性能進(jìn)行分析研究。在校準(zhǔn)箱內(nèi)溫度為-25～0℃范圍內(nèi)，循環(huán)倍率在2～5范圍內(nèi)變化時(shí)，冷風(fēng)機(jī)總換熱系數(shù)隨著校準(zhǔn)箱溫度的升高而增大;制冷工質(zhì)為CO2時(shí)冷風(fēng)機(jī)的制冷量明顯高于制冷工質(zhì)為NH3時(shí)，在校準(zhǔn)箱內(nèi)溫度為0℃時(shí)高42%，-20℃時(shí)高26%;管內(nèi)側(cè)壓降隨著循環(huán)倍率的增大而增大;換熱系數(shù)隨著循環(huán)倍率的增大先增大后逐漸減小，在循環(huán)倍率為3左右時(shí)，換熱系數(shù)達(dá)到最大。仿真結(jié)果與測(cè)試結(jié)果趨勢(shì)相同，但存在一定誤差。模擬計(jì)算得出NH3換熱系數(shù)值與測(cè)試結(jié)果的誤差約為16%，CO2換熱系數(shù)值與測(cè)試結(jié)果的誤差約為8%。

換熱系數(shù);管內(nèi)壓降;校準(zhǔn)箱溫度;循環(huán)倍率

AbstractIn this study a cooling-fan simulation model was established using the steady-state distributed parameter method，and the performance of a cooling-fan prototype was tested in an air-blower performance test rig.Further，a cooler heat transfer performance analysis was conducted based on the experimental data and numerical simulation.For a temperature calibration of-25-0℃and a circulation rate varying in the range of 2-5，it was found that the total heat transfer coefficient increases with increasing calibration-box temperature.When the refrigerant is CO2， the cooling capacity of the cooling fan is obviously higher than that for an NH3refrigerant.Further， the cooling capacity is 42%higher when the temperature in the tank is 0℃and 26%higher for a tank temperature of-20℃.The pressure drop of inner side of the tube increases with increases in the circulation rate.Further，the heat transfer coefficient first increases and then decreases with increases in the circulation rate.The circulation ratio is approximately 3，corresponding to the maximum change in the thermal coefficient.The simulation results exhibit the same trend as the test results，but some errors exist.The error of heat transfer coefficient between the numerical and test results is approximately 16%for the NH3system，and approximately 8%for the CO2system.

Keywordsheat transfer coefficient;pressure drop inside the tube;calibration temperatures;circulation rate

冷風(fēng)機(jī)在制冷系統(tǒng)中作為蒸發(fā)器，依靠軸流風(fēng)機(jī)使冷庫內(nèi)空氣與排管翅片發(fā)生強(qiáng)制對(duì)流換熱，實(shí)現(xiàn)冷卻降溫的效果。與傳統(tǒng)頂排管相比，具有1)傳熱效率高:由于采用銅管(或者鋼管)與翅片緊密接觸，管片成為一體，大大增加了換熱面積，從而使傳熱效率增高;2)防腐性能優(yōu)異:冷風(fēng)機(jī)采用熱浸鋅防腐，是電鍍防腐性能的10倍至15倍;3)能耗小、噪音低;4)結(jié)構(gòu)型式多:冷風(fēng)機(jī)型式有落地頂吹式、落地上側(cè)吹風(fēng)式、落地平吹式、吊頂式、拼裝式等;5)適用范圍廣:可廣泛應(yīng)用于肉、蛋、瓜果、蔬菜、海產(chǎn)品等冷庫內(nèi)。冷風(fēng)機(jī)的性能對(duì)制冷系統(tǒng)的運(yùn)行與冷庫內(nèi)食品品質(zhì)都會(huì)造成十分重要的影響[1-2]。

然而，目前我國(guó)對(duì)CO2冷風(fēng)機(jī)性能的研究還相對(duì)較少。當(dāng)前在工程上選型計(jì)算所使用的方法均是按照冷風(fēng)機(jī)的標(biāo)牌名義制冷量上乘以一定系數(shù)，這種選型方式會(huì)導(dǎo)致CO2冷風(fēng)機(jī)的選型誤差，也不能保證CO2冷風(fēng)機(jī)能一直處于高效的運(yùn)行工況[3]。

隨著現(xiàn)代計(jì)算機(jī)科學(xué)與技術(shù)的發(fā)展，以冷風(fēng)機(jī)內(nèi)部傳熱傳質(zhì)機(jī)理作為理論依據(jù)，用計(jì)算機(jī)仿真計(jì)算方法來實(shí)現(xiàn)冷風(fēng)機(jī)性能預(yù)測(cè)和產(chǎn)品開發(fā)，不僅降低對(duì)實(shí)際樣機(jī)測(cè)試的依賴程度，還可大大提高冷風(fēng)機(jī)性能預(yù)測(cè)的效率和準(zhǔn)確性，因此冷風(fēng)機(jī)性能的仿真研究已經(jīng)成為現(xiàn)代產(chǎn)品設(shè)計(jì)的主要發(fā)展方向。

P.Jostein等[4-7]對(duì)CO2在微通道內(nèi)的蒸發(fā)過程進(jìn)行了理論分析和實(shí)驗(yàn)研究，并對(duì)蒸發(fā)過程流型變化進(jìn)行了可視化研究與分析，其實(shí)驗(yàn)段由25個(gè)圓形的通道組成，其中每個(gè)通道的內(nèi)徑為0.79 mm。主要研究?jī)?nèi)容為質(zhì)量流速、蒸發(fā)溫度及熱流密度等因素對(duì)換熱系數(shù)的影響。

J.Patino等[8]對(duì)多篇文獻(xiàn)中CO2兩相區(qū)管內(nèi)流動(dòng)沸騰換熱計(jì)算公式進(jìn)行總結(jié)，并建立數(shù)學(xué)模型進(jìn)行仿真計(jì)算，讓其結(jié)果與實(shí)驗(yàn)值進(jìn)行對(duì)比分析，發(fā)現(xiàn)用Cheng Lixin等[9-10]的流型圖法計(jì)算結(jié)果誤差最小，對(duì)CO2制冷能力的計(jì)算誤差為8.4%，對(duì)制冷工質(zhì)出口溫度的計(jì)算誤差為0.2%，制冷工質(zhì)出口壓力的計(jì)算誤差為0.4%。

Sun Z.等[11-13]對(duì)在水平管內(nèi)CO2蒸發(fā)過程中流動(dòng)的沸騰換熱特性、流型變化及換熱關(guān)聯(lián)式進(jìn)行研究，對(duì)CO2蒸發(fā)時(shí)干涸的預(yù)測(cè)以及在干涸發(fā)生前、后及干涸區(qū)三個(gè)流動(dòng)區(qū)域的換熱計(jì)算均分別建立了模型。

1 仿真模型及實(shí)驗(yàn)裝置

考慮到制冷工質(zhì)在管內(nèi)的蒸發(fā)換熱過程為相變換熱，同時(shí)為了更好的反映冷風(fēng)機(jī)的實(shí)際換熱情況，根據(jù)制冷工質(zhì)在蒸發(fā)器內(nèi)的實(shí)際流動(dòng)及換熱狀態(tài)對(duì)冷風(fēng)機(jī)建立仿真計(jì)算模型，并分析確定制冷劑狀態(tài)參數(shù)，當(dāng)制冷劑在蒸發(fā)器內(nèi)處于過熱時(shí)，將蒸發(fā)器模型進(jìn)行分區(qū)，并在各區(qū)內(nèi)采用穩(wěn)態(tài)分布參數(shù)法進(jìn)行計(jì)算。為了簡(jiǎn)化模型計(jì)算并方便模型建立，進(jìn)行以下假設(shè):1)冷風(fēng)機(jī)內(nèi)制冷劑為穩(wěn)態(tài)運(yùn)行;2)管內(nèi)側(cè)制冷工質(zhì)及管外濕空氣的流動(dòng)均為一維穩(wěn)態(tài)流動(dòng);3)管內(nèi)工質(zhì)和管外濕空氣在各點(diǎn)的流量不隨時(shí)間發(fā)生變化，流量均勻分布;4)換熱管管壁徑向溫度近似一致，即忽略換熱管管壁的傳熱熱阻;5)所有熱量均在制冷工質(zhì)跟空氣間進(jìn)行，忽略軸向傳熱;6)對(duì)于沿管長(zhǎng)方向的每個(gè)控制單元，其濕空氣側(cè)、制冷工質(zhì)側(cè)以及管壁上的物性參數(shù)均近似一致。

根據(jù)冷風(fēng)機(jī)實(shí)際的結(jié)構(gòu)參數(shù)，每一個(gè)管程的換熱管及其翅片均作為一個(gè)單獨(dú)的控制單元進(jìn)行仿真計(jì)算。將每個(gè)管程的管長(zhǎng)均分成八排，前一排出口狀態(tài)參數(shù)作為下一排進(jìn)口狀態(tài)參數(shù)。在兩相區(qū)內(nèi)，單個(gè)控制單元根據(jù)制冷劑CO2的流動(dòng)方向進(jìn)行劃分，由于兩相區(qū)內(nèi)制冷劑焓值和壓降變化幅度較大，因此在兩相區(qū)內(nèi)按照干度將管長(zhǎng)分成N個(gè)控制單元。在過熱區(qū)，制冷劑的溫度變化較大，在此階段將按照剩余的管長(zhǎng)劃分成M個(gè)控制單位，采用同樣的方法進(jìn)行計(jì)算。最后將各個(gè)微元控制體的傳熱模型疊加，得到完整的傳熱模型。即可得到局部傳熱系數(shù)、傳熱系數(shù)、管內(nèi)壓降、制冷量等參數(shù)的變化趨勢(shì)?？刂茊卧膭澐秩鐖D1所示。

圖1 控制單元的劃分Fig.1 Division of the control-unit

在每個(gè)控制單元內(nèi)，制冷劑側(cè)的沸騰換熱過程與空氣側(cè)的冷卻過程遵循能量守恒定律。其中，制冷劑側(cè)的換熱方程為:

相對(duì)應(yīng)的，空氣側(cè)的換熱方程為:

制冷劑側(cè)的換熱量與空氣側(cè)的換熱量的轉(zhuǎn)換關(guān)系為:在冷風(fēng)機(jī)性能測(cè)試實(shí)驗(yàn)過程中，制冷劑在換熱管內(nèi)沸騰換熱，不只是吸收冷庫內(nèi)的熱負(fù)荷，同時(shí)還吸收換熱管管壁等其他裝置的熱負(fù)荷，從而造成制冷劑側(cè)的換熱量與空氣側(cè)的換熱量不可能完全相等。根據(jù)相關(guān)的文獻(xiàn)和大量的實(shí)驗(yàn)數(shù)據(jù)，可得到兩者的關(guān)聯(lián)式，如式(3)所示:

式中:r根據(jù)運(yùn)行工況確定，本文取值0.9。

冷風(fēng)機(jī)的總傳熱系數(shù)的計(jì)算公式與空氣側(cè)工況有關(guān)?？諝鈧?cè)工況大致分為三類:干工況(露點(diǎn)溫度以上的工況)、凝露的工況以及結(jié)霜的工況。因此應(yīng)用于制冷系統(tǒng)中的冷風(fēng)機(jī)，在仿真計(jì)算和實(shí)驗(yàn)研究的過程中為干工況，總傳熱系數(shù)的計(jì)算公式如下:

干工況:

實(shí)驗(yàn)所用制冷系統(tǒng)為泵供液NH3制冷系統(tǒng)中和NH3做制冷劑CO2做載冷劑的載冷循環(huán)系統(tǒng)。系統(tǒng)原理如圖2所示。

圖2 系統(tǒng)原理Fig.2 The principle of refrigeration system

當(dāng)采用NH3泵供液制冷循環(huán)時(shí)，打開截止閥9、12、15，關(guān)閉截止閥 8、17、20，經(jīng)過液體調(diào)節(jié)站后NH3液體分為兩路，一路供環(huán)境間制冷系統(tǒng)，另一路供被測(cè)冷風(fēng)機(jī)22。低壓循環(huán)桶11中的液NH3通過循環(huán)泵13后向被測(cè)冷風(fēng)機(jī)供液，換熱后氣液兩相的NH3回到氣液分離器進(jìn)行氣液分離，分離后的氣體部分經(jīng)過熱器過熱后進(jìn)入壓縮機(jī)，開始下一次循環(huán)。

當(dāng)采用CO2載冷循環(huán)制冷系統(tǒng)時(shí)，打開截止閥8、17、20，關(guān)閉截止閥 9、12、15，氨高壓儲(chǔ)液器 4 出來的液NH3分成兩路，分別向環(huán)境間和冷凝蒸發(fā)器6供液。在冷凝蒸發(fā)器中，節(jié)流后的低溫NH3液體與CO2進(jìn)行換熱。冷凝后的CO2進(jìn)入儲(chǔ)液器16，然后經(jīng)CO2泵進(jìn)入被測(cè)冷風(fēng)機(jī)22進(jìn)行蒸發(fā)換熱。從被測(cè)冷風(fēng)機(jī)出來的CO2進(jìn)入冷凝蒸發(fā)器6被冷凝后進(jìn)入儲(chǔ)液器開始下一循環(huán)。仿真計(jì)算及實(shí)驗(yàn)工況如表1和表2所示，所用冷風(fēng)機(jī)結(jié)構(gòu)參數(shù)見表3。

表1 不同校準(zhǔn)箱溫度時(shí)的實(shí)驗(yàn)工況Tab.1 Experimental operating conditions under different calibration temperature

表2 不同循環(huán)倍率時(shí)的實(shí)驗(yàn)工況Tab.2 Experimental operating conditions under different circulating ratios

實(shí)驗(yàn)測(cè)試所用冷風(fēng)機(jī)為煙臺(tái)冰輪股份有限公司的實(shí)驗(yàn)樣機(jī)，以通過實(shí)驗(yàn)方式準(zhǔn)確測(cè)量冷風(fēng)機(jī)實(shí)際性能。實(shí)驗(yàn)用冷風(fēng)機(jī)材料為鋁管鋁翅片，換熱管布置為叉排，主要結(jié)構(gòu)尺寸的具體參數(shù)如表3所示。

表3 冷風(fēng)機(jī)結(jié)構(gòu)參數(shù)Tab.3 The structural parameters of air coolers

2 實(shí)驗(yàn)結(jié)果分析及與仿真結(jié)果對(duì)比

實(shí)驗(yàn)過程中，在不同工況下對(duì)冷風(fēng)機(jī)風(fēng)量和風(fēng)速進(jìn)行測(cè)定，測(cè)量結(jié)果顯示風(fēng)量和風(fēng)速分別在12 000 m3/h、2.4 m/s左右變化，波動(dòng)范圍≤5%，故冷風(fēng)機(jī)風(fēng)量和風(fēng)速按以上測(cè)定值計(jì)算。

圖3給出了制冷工質(zhì)分別為NH3和CO2時(shí)，冷風(fēng)機(jī)制冷量隨校準(zhǔn)箱溫度的變化情況。由圖可知，在校準(zhǔn)箱溫度為-25～0℃范圍內(nèi)，冷風(fēng)機(jī)制冷量隨著校準(zhǔn)箱內(nèi)溫度的降低而減小，符合制冷系統(tǒng)的經(jīng)典規(guī)律。制冷工質(zhì)為NH3時(shí)，變化范圍在11.83～12.88 kW，制冷量的下降梯度較小;制冷工質(zhì)為CO2時(shí)，變化范圍在15.3～18.29 kW，制冷量的下降梯度較大。制冷工質(zhì)為CO2時(shí)冷風(fēng)機(jī)的制冷量明顯高于制冷工質(zhì)為NH3時(shí)，在校準(zhǔn)箱內(nèi)溫度為0℃時(shí)高42%，-20℃時(shí)高26%，實(shí)驗(yàn)溫度范圍內(nèi)換熱系數(shù)平均高35%左右。分析其原因，在其他條件相同的情況下，管內(nèi)CO2流動(dòng)沸騰換熱系數(shù)明顯高于NH3，且CO2的蒸發(fā)潛熱也明顯大于NH3的值。

圖3 制冷量隨校準(zhǔn)箱溫度的變化Fig.3 The refrigerating capacity change over the calibration temperature

對(duì)不同校準(zhǔn)箱溫度下冷風(fēng)機(jī)換熱性能進(jìn)行仿真模擬計(jì)算，并將仿真計(jì)算結(jié)果與測(cè)試數(shù)據(jù)進(jìn)行對(duì)比，結(jié)果如圖4所示。由圖可知，仿真與實(shí)驗(yàn)中換熱系數(shù)均隨著校準(zhǔn)箱內(nèi)溫度的升高而增加，變化趨勢(shì)相似且仿真模擬結(jié)果均比實(shí)驗(yàn)結(jié)果偏大。制冷工質(zhì)為NH3時(shí)，仿真計(jì)算結(jié)果顯示換熱系數(shù)在30.08～31.85 W/(m2·K)范圍內(nèi)變化，實(shí)驗(yàn)結(jié)果顯示換熱系數(shù)在25.32～27.25 W/(m2·K)范圍內(nèi)變化，仿真計(jì)算結(jié)果與實(shí)驗(yàn)測(cè)試誤差為18%左右。制冷工質(zhì)為CO2時(shí)，換熱系數(shù)在28.73～30.82 W/(m2·K)范圍內(nèi)變化，不同計(jì)算式的計(jì)算結(jié)果誤差并不相同，其中D.S.Jung等[14-16]計(jì)算公式的仿真結(jié)果最為相近，與實(shí)驗(yàn)結(jié)果誤差小于16%，用S.H.Yoon等[17]的計(jì)算公式的仿真結(jié)果與實(shí)驗(yàn)結(jié)果相比平均誤差為15%左右，用Cheng Lixin等[9-10]的流型圖法的仿真模型計(jì)算結(jié)果與測(cè)試結(jié)果誤差為9%左右。由幾種CO2計(jì)算公式模擬結(jié)果對(duì)比可知，在CO2管內(nèi)流動(dòng)沸騰換熱關(guān)聯(lián)式中，流型圖法計(jì)算結(jié)果誤差最小。經(jīng)過以上分析，本文選用精度較高的流型圖法作為CO2管內(nèi)流動(dòng)換熱系數(shù)計(jì)算方法。

分析所有仿真模擬的結(jié)果均比實(shí)驗(yàn)結(jié)果偏大的原因可能是在工程實(shí)踐中，由于冷風(fēng)機(jī)內(nèi)換熱管的脹接和焊接中操作不符合工藝等導(dǎo)致?lián)Q熱管與翅片間不完全接觸，從而產(chǎn)生接觸熱阻;在冷風(fēng)機(jī)長(zhǎng)期使用過程中換熱管內(nèi)可能存在污垢，產(chǎn)生污垢熱阻;在換熱器實(shí)際運(yùn)行過程中，制冷工質(zhì)中夾帶油從而降低了制冷工質(zhì)管內(nèi)換熱。而以上這些原因在仿真模型計(jì)算中均沒有充分考慮，從而導(dǎo)致仿真模擬計(jì)算出的冷風(fēng)機(jī)換熱系數(shù)大于實(shí)驗(yàn)測(cè)量計(jì)算得出的冷風(fēng)機(jī)的換熱系數(shù)。

管內(nèi)側(cè)壓降隨循環(huán)倍率的變化趨勢(shì)如圖5可知。由圖可知，仿真和實(shí)驗(yàn)結(jié)果相似，載冷循環(huán)的循環(huán)倍率在2～5范圍內(nèi)變化時(shí)，隨著循環(huán)倍率的增加，管內(nèi)制冷工質(zhì)質(zhì)量流量和壓降均逐漸增加，兩者誤差約為30%。誤差較大的原因:一方面仿真模擬過程中，為方便計(jì)算對(duì)實(shí)際情況做了理想化的前提假設(shè)，與實(shí)際情況有所不同;另一方面，實(shí)驗(yàn)用冷風(fēng)機(jī)中換熱管內(nèi)表面不光滑、有一定粗糙度，管內(nèi)工質(zhì)不純混入潤(rùn)滑油等因素。

圖4 不同校準(zhǔn)箱溫度時(shí)仿真模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)對(duì)比Fig.4 Comparison of simulation results and experiment data under different calibration temperatures

圖5 管內(nèi)側(cè)壓降隨循環(huán)倍率的變化Fig.5 Change of pressure drop inside the tube over the circulate ratio

圖6所示為不同循環(huán)倍率時(shí)冷風(fēng)機(jī)換熱系數(shù)的模擬與實(shí)驗(yàn)結(jié)果。制冷工質(zhì)為NH3及CO2時(shí)的變化趨勢(shì)相似，當(dāng)循環(huán)倍率在2～5范圍內(nèi)變化時(shí)，實(shí)驗(yàn)測(cè)試結(jié)果顯示:隨著循環(huán)倍率的增加，換熱系數(shù)先升高后下降，在循環(huán)倍率為3有最大值;模擬結(jié)果顯示隨著循環(huán)倍率的增加冷風(fēng)機(jī)換熱系數(shù)先增大后趨于穩(wěn)定，在循環(huán)倍率超過3時(shí)，換熱系數(shù)增加梯度相對(duì)較小，其絕對(duì)值趨于穩(wěn)定。模擬結(jié)果與實(shí)驗(yàn)結(jié)果的差異主要是因?yàn)樵诜抡婺Ｐ椭杏?jì)算換熱系數(shù)時(shí)忽略了因?yàn)楣軆?nèi)壓降引起的管內(nèi)制冷工質(zhì)溫度的降低。

圖6 換熱系數(shù)隨循環(huán)倍率的變化Fig.6 Change of heat transfer coefficient over the circulate ratio

3 結(jié)論

本文對(duì)冷風(fēng)機(jī)的換熱性能進(jìn)行了建模分析與實(shí)驗(yàn)測(cè)試得出以下結(jié)論:

在校準(zhǔn)箱內(nèi)溫度為-25～0℃范圍內(nèi)，載冷系統(tǒng)多倍供液，工質(zhì)分別為NH3與CO2時(shí)均顯示:1)冷風(fēng)機(jī)總換熱系數(shù)隨著校準(zhǔn)箱溫度的升高而增大;2)載冷循環(huán)的循環(huán)倍率在2～5范圍內(nèi)變化時(shí)，管內(nèi)側(cè)壓降隨著循環(huán)倍率的增大而增大;3)當(dāng)循環(huán)倍率在2～5范圍內(nèi)變化時(shí)，換熱系數(shù)隨著循環(huán)倍率的增大先增大后逐漸減小，在循環(huán)倍率為3左右時(shí)，換熱系數(shù)達(dá)到最大。

通過對(duì)仿真模擬計(jì)算結(jié)果與測(cè)試結(jié)果的對(duì)比分析可知，仿真結(jié)果與測(cè)試結(jié)果趨勢(shì)相同，但存在一定誤差。模擬計(jì)算得出NH3換熱系數(shù)值與測(cè)試結(jié)果的誤差為16%左右，CO2換熱系數(shù)值與測(cè)試結(jié)果的誤差為8%左右?？傮w來說，計(jì)算值與實(shí)驗(yàn)值的變化趨勢(shì)基本相同，二者吻合程度較高，表明仿真計(jì)算模型對(duì)實(shí)驗(yàn)工況變化引起的冷風(fēng)機(jī)換熱性能變化具有較好的仿真能力，可以滿足工程上的應(yīng)用。

[1]申江，楊萌.食品冷鏈的技術(shù)發(fā)展[J].包裝工程，2015，36(15):1-9.(SHEN Jiang，YANG Meng.Technological development of food cold chain[J].Packaging Engeering，2015，36(15):1-9.)

[2]國(guó)家發(fā)展和改革委員會(huì).農(nóng)產(chǎn)品冷鏈物流發(fā)展規(guī)劃[R].北京:國(guó)家發(fā)展和改革委員會(huì)，2010.(National Development and Reform Commission.Agricultural products cold chain logistics development planning[R].Beijing: NationalDevelopmentand Reform Commission，2010.)

[3]申江，黃冰.CO2冷風(fēng)機(jī)性能研究[D].天津:天津商業(yè)大學(xué)，2015.(SHEN Jiang， HUANG Bing.Performance study of the CO2air-cooler[D].Tianjin:Tianjin University of Commerce，2015.)

[4]JOSTEIN P.Flow vaporization of CO2in micro-channel tubes[J].Experimental Thermal and Fluid Science，2004，28(5):111-121.

[5]PETTERSEN J， RIEBERER R， MUNKEJORD S T.Heat transfer and pressure drop characteristics of evaporating carbon dioxide in micro-channel tubes[C]//Preliminary Proceedings of the 4th IIR-Gustav Lorentzen Conference on Natural Working Fluids.Purdue:Purdue University， 2000，22(7):279-284.

[6]PETTERSEN J.Flow vaporization of CO2in micro-channel tubes，Part 1:Experimental method and two-phase flow pattern[C]//Preliminary Proceedings of the 5th IIR-Gustav Lorentzen Conference on Natural Working Fluids.Guangzhou，2002.

[7]PETTERSEN J.Flow vaporization of CO2in micro-channel tubes，Part 2:Heat transfer， pressure drop and correlation[C]//PreliminaryProceedingsofthe 5th IIR-Gustav Lorentzen Conference on Natural Working Fluids.Guangzhou，2002.

[8]PATINO J，LLOPIS R，SANCHEZ D，et al.A comparative analysis of a CO2evaporator model using experimental heat transfer correlations and a flow pattern map[J].International Journal of Heat and Mass Transfer，2014，71:361-375.

[9]CHENG Lixin，RIBATSKI G，QUIBE′N J M，et al.New prediction methods for CO2evaporation inside tubes:Part I—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/2):111-124.

[10]CHENG Lixin，RIBATSKI G，THOME J R.New prediction methods for CO2evaporation inside tubes:Part II—An updated general flow boiling heat transfer model based on flow patterns[J].International Journal of Heat and Mass Transfer，2008，51(1/2):125-135.

[11]SUN Z，GROLL E A.CO2flow boiling heat transfer in horizontal tubes，Part 1:Flow regime and prediction of dry-out[C]//PreliminaryProceedingsofthe 5th IIR-Gustav Lorentzen Conference on Natural Working Fluids.Guangzhou，2002:116-125.

[12]SUN Z，GROLL E A.CO2flow boiling heat transfer in horizontal tubes，Part 2:Experimental[C]//Preliminary Proceedings of the 5th IIR-Gustav Lorentzen Conference onNatural Working Fluids.Guangzhou，2002:126-132.

[13]SUN Z，GROLL E A.CO2flow boiling heat transfer in horizontal tubes，Part 1:Prediction of heat transfer coefficient[C]//PreliminaryProceedingsofthe 5th IIR-Gustav Lorentzen Conference on Natural Working Fluids.Guangzhou，2002:133-142.

[14]JUNG D S，MCLINDEN M，RADERMACHER R，et al.A study of flow boiling heat transfer with refrigerant mixtures[J].International Journal of Heat and Mass Transfer，1989，32(9):1751-1764.

[15]GUNGOR K E，WINTERTON R H S.A general correlation for flow boiling in tubes and annuli[J].Heat Mass Transfer， 1986，29(3):351-358.

[16]KANDLIKAR S G.A general correlation for saturated twophase flow boiling heat transfer inside horizontal and vertical tubes[J].Journal of Heat Transfer， 1990，112(1):219-228.

[17]YOON S H，CHO E S， HWANG Y W，et al.Characteristics of evaporative heat transfer and pressure drop of carbon dioxide and correlation development[J].International Journal of Refrigeration，2004，27(2):111-119.

Simulation and Experimental Study on Heat Transfer Performance of CO2Air Cooler

Shen Jiang Bian Yujun Huang Bing
(Refrigeration Key Laboratory of Tianjin， Tianjin University of Commerce， Tianjin， 300134， China)

TB657;TP391.9

A

2016年11月2日

0253-4339(2017)05-0051-06

10.3969/j.issn.0253-4339.2017.05.051

邊煜竣，男，碩士研究生，天津商業(yè)大學(xué)機(jī)械工程學(xué)院，13682126797，E-mail:424381844@qq.com。 研究方向:二氧化碳制冷。

About the corresponding authorBian Yujun， male， graduate， School of Mechanical Engineering，Tianjin University of Commerce， +86 13682126797， E-mail:424381844@qq.com.Research fields:energy saving and optimization of refrigeration system.