姬 旭，楊俊蘭，李金芮，白 楊

（天津城建大學(xué)能源與安全工程學(xué)院，天津 300384）

進(jìn)入21世紀(jì)以來，節(jié)能和環(huán)保是世界各國(guó)政府和學(xué)者普遍關(guān)注的問題[1-2].我國(guó)將節(jié)約資源與環(huán)境保護(hù)列為基本國(guó)策[3]，為走可持續(xù)發(fā)展道路，在空調(diào)領(lǐng)域回歸自然工質(zhì)作為制冷劑是大勢(shì)所趨[4-6].CO2是自然工質(zhì)的代表，其性質(zhì)無毒不燃，ODP為0，GWP為1，容積制冷量大，導(dǎo)熱性良好，并且價(jià)格低廉，但是CO2跨臨界系統(tǒng)工作壓力遠(yuǎn)高于一般制冷系統(tǒng)，這就對(duì)換熱器的結(jié)構(gòu)有較高的要求，將CO2跨臨界系統(tǒng)與微小通道換熱器有機(jī)結(jié)合，既可以解決微小通道內(nèi)壓降大，堵塞通道或流體流動(dòng)分配不均等問題[7]，又可以滿足CO2空調(diào)制冷系統(tǒng)的高效、緊湊等要求.

影響CO2在管內(nèi)沸騰換熱的因素較多，沸騰換熱彼此間存在較大差異.Bredesen等[8]研究了CO2在直管內(nèi)的沸騰換熱和壓降特性，表明蒸發(fā)溫度為5℃時(shí)出現(xiàn)干涸現(xiàn)象，并指出與其他制冷劑相比，CO2的核態(tài)沸騰換熱作用明顯.Jang和Hrnjak[9]通過實(shí)驗(yàn)對(duì)CO2常規(guī)通道圓管內(nèi)和微小通道內(nèi)的換熱機(jī)理進(jìn)行了深入研究，表明隨著干度的增加，流型從彈狀流向塞狀流最后向環(huán)狀流的轉(zhuǎn)變.Yoon等[10]對(duì)常規(guī)通道水平管內(nèi)的沸騰換熱進(jìn)行了實(shí)驗(yàn)研究，表明CO2在中低干度區(qū)，沸騰換熱系數(shù)隨熱流密度的增大而增大.Yun等[11-14]通過實(shí)驗(yàn)研究表明CO2換熱系數(shù)比相同工況條件下的R134a高了47%，并將干涸區(qū)域分成過渡區(qū)域和完全干涸區(qū)域兩部分.張良等[15]研究表明，在中低干度區(qū)，CO2沸騰換熱系數(shù)隨質(zhì)量流速、熱流密度的增加而增加，在高干度區(qū)，隨熱流密度增加發(fā)生干涸時(shí)的干度降低，換熱系數(shù)下降.

近年來，國(guó)內(nèi)外研究學(xué)者對(duì)大量的CO2流動(dòng)沸騰換熱預(yù)測(cè)模型和CO2微通道換熱器進(jìn)行了針對(duì)性的分析.對(duì)于常規(guī)通道來說，Yoon等[10]和Oh等[16]研究了CO2在水平圓管下的流動(dòng)沸騰換熱情況，并與Gungor[17]和Jung[18]等的預(yù)測(cè)模型進(jìn)行分析對(duì)比，Yoon計(jì)算結(jié)果表明Gungor模型精度最高為34.8%，而Oh等人的研究表明Jung的預(yù)測(cè)關(guān)聯(lián)式誤差最小為21.6%.Ducoulombier等[19]研究了CO2在水平圓管的流動(dòng)沸騰換熱情況，表明對(duì)于微通道kandlikar[20]模型誤差最高為43.04%.Cheng[21]和Fang[22-23]都是通過搜集大量的CO2管內(nèi)流動(dòng)沸騰換熱實(shí)驗(yàn)數(shù)據(jù)點(diǎn)，并對(duì)多個(gè)沸騰換熱關(guān)聯(lián)式進(jìn)行誤差分析，建立了流動(dòng)沸騰換熱實(shí)驗(yàn)數(shù)據(jù)庫.結(jié)果表明，Cheng的模型對(duì)數(shù)據(jù)庫中的霧狀流區(qū)預(yù)測(cè)結(jié)果較為準(zhǔn)確，而Fang預(yù)測(cè)的沸騰換熱關(guān)聯(lián)式的誤差僅為15.5%.

本文通過模擬研究微細(xì)通道內(nèi)CO2沸騰換熱的特性，為后續(xù)微細(xì)通道換熱器的研制提供技術(shù)參考.

1 模型與數(shù)值方程

1.1 幾何模型的建立

圖1為研究對(duì)象的幾何模擬圖，水平圓管管內(nèi)徑d為2 mm，管長(zhǎng)為1 m.CO2制冷劑在通道內(nèi)自左向右流動(dòng).

圖1 研究對(duì)象的幾何模型圖

1.2 獨(dú)立性驗(yàn)證

為了簡(jiǎn)化模型，選取網(wǎng)格數(shù)較少的進(jìn)行數(shù)值模擬，同時(shí)也能滿足模擬的精確性.網(wǎng)格無關(guān)性驗(yàn)證選取工況為飽和溫度Tsat=283 K，熱流密度q=6 kW/m2，質(zhì)量流速G=200 kg/（m2·s），干度x0為0~1時(shí)的沸騰換熱系數(shù)h.結(jié)果表明在不同干度時(shí)，隨著網(wǎng)格數(shù)量的增加，沸騰換熱系數(shù)變化不大，相鄰網(wǎng)格密度之間的數(shù)值結(jié)果誤差在5%以內(nèi)，可以認(rèn)為此時(shí)計(jì)算結(jié)果對(duì)網(wǎng)格數(shù)量不敏感，計(jì)算結(jié)果可信度高，足夠準(zhǔn)確（見圖2）.

圖2 沸騰換熱系數(shù)隨網(wǎng)格數(shù)的變化

1.3 控制方程

本研究選用多相流VOF模型[24]，對(duì)于局部的流動(dòng)區(qū)域，可以如下定義汽相（v）和液相（L）之間的質(zhì)量轉(zhuǎn)移

式（1）-（2）中：T為溫度，K；Tsat為飽和溫度，K；α為體積分?jǐn)?shù)；mvl為冷凝時(shí)質(zhì)量傳遞速率，kg/（m3·s）；mlv為蒸發(fā)時(shí)質(zhì)量傳遞速率，kg/（m3·s）；ρ為密度，kg/m3.

1.4 求解條件設(shè)定

如圖1所示，本次模擬中，壁面采用均勻熱流邊界條件，工質(zhì)自左向右流過水平圓管，邊界條件入口選用速度入口，出口選用壓力出口.入口的質(zhì)量流速G為300~600 kg/（m2·s），操作壓力為常壓，入口溫度T為253~273 K，壁面加熱的熱流密度q為7.5~30 kW/m2，水平圓管內(nèi)徑d為2 mm.

在求解問題時(shí)，模擬選用的離散化方法為有限容積法（FVM），其中選用二階迎風(fēng)格式.二階迎風(fēng)格式相對(duì)于一階格式來說，精度更高，使用更小的截?cái)嗾`差，適用于流動(dòng)與網(wǎng)格不在同一直線上的情況.

流動(dòng)沸騰換熱是較為復(fù)雜的過程，此過程可以按照湍流計(jì)算.本研究采用標(biāo)準(zhǔn)湍流模型.

1.5 驗(yàn)證模型選擇

本文選取Fang的CO2管內(nèi)流動(dòng)沸騰換熱模型作為驗(yàn)證模型[22].關(guān)聯(lián)式如下

式（3）-（5）中：Nu為努謝爾數(shù)；F為單位體積所受的表面張力，N；Re1為雷諾數(shù)；Fa為無因次數(shù)；Pr為普朗克數(shù)；μ1，f為流體溫度下的液體黏度，kg/（m·s）；μ1，w為內(nèi)壁表面溫度下的液體黏度，kg/（m·s）；Bo為沸點(diǎn)；x1為蒸汽質(zhì)量占比.

式（5）中a的取值由Re1·Fa決定

關(guān)聯(lián)式數(shù)據(jù)參數(shù)區(qū)間：d=0.529~7.750 mm，q=3.39~40 kW/m2，G=97.5~1 400 kg/（m2·s），Tsat=-40~26.8℃，x1=0.004 6~0.998.

1.6 驗(yàn)證模型

1.6.1 熱流密度的影響

通過選用不同的工況來多角度地驗(yàn)證模型，所選CO2質(zhì)量流速G為200 kg/（m2·s），飽和溫度283 K，操作壓力在4.4 MPa，探究熱流密度q在3，6 kW/m2的條件下，沸騰換熱系數(shù)隨干度x的變化如圖3所示.

圖3 不同熱流密度下沸騰換熱系數(shù)隨干度的變化

由圖3a可知，當(dāng)q=3 kW/m2時(shí)在0~0.1的干度區(qū)間內(nèi)，模擬值與預(yù)測(cè)值的相對(duì)誤差不大于18%；由圖3b可知，當(dāng)q=6 kW/m2，干度小于0.12時(shí)，模擬值與預(yù)測(cè)值的相對(duì)誤差不大于20%，干度大于0.12時(shí)誤差逐漸增大.

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

選CO2的飽和溫度283 K，操作壓力在4.4 MPa，熱流密度q為6 kW/m2，探究質(zhì)量流速G在200，400 kg/（m2·s）的條件下，表面?zhèn)鳠嵯禂?shù)隨干度變化.

由圖4a可知，當(dāng)G=400 kg/（m2·s）時(shí)在0~0.1的干度區(qū)間內(nèi)，模擬值與預(yù)測(cè)值的相對(duì)誤差不大于19%；由圖4b可知，當(dāng)G=200 kg/（m2·s）且干度x小于0.10時(shí)，模擬值與預(yù)測(cè)值的相對(duì)誤差不大于16%，干度大于0.10時(shí)誤差逐漸增大.

通過對(duì)模型的不斷修正，得到在較低干度情況下模型的模擬值與預(yù)測(cè)值的相對(duì)誤差在15%左右，總體相對(duì)誤差控制在20.25%以內(nèi)，驗(yàn)證了模型的準(zhǔn)確性.

圖4 不同質(zhì)量流速下沸騰換熱系數(shù)隨干度的變化

2 結(jié)果分析

2.1 飽和溫度對(duì)傳熱系數(shù)的影響

在熱流密度q=7.5 kW/m2，質(zhì)量流速G=300 kg/（m2·s），飽和溫度Tsat在253~273 K的工況下，圖5為CO2工質(zhì)在不同飽和溫度下水平管進(jìn)口和出口的液相容積分配圖.由圖5a可知，在飽和溫度為253 K時(shí)，水平圓管入口處的氣泡數(shù)量很少，且都是小氣泡，到出口位置附近干度在0.1左右時(shí)，氣泡的數(shù)量增多，尺寸變得大而扁平，沸騰現(xiàn)象開始明顯，換熱強(qiáng)烈.當(dāng)飽和溫度為263 K時(shí)，可以觀察到在同一位置出入口附近，CO2工質(zhì)所產(chǎn)生的氣泡數(shù)量相比于飽和溫度為253 K時(shí)有所增加，表明CO2工質(zhì)的傳熱系數(shù)也有明顯的增加.當(dāng)飽和溫度繼續(xù)增大到273 K時(shí)，水平圓管進(jìn)出口的氣泡變得更多，可見在低干度條件下隨著飽和溫度的不斷升高，換熱現(xiàn)象逐漸強(qiáng)烈，沸騰換熱系數(shù)也隨之增大.

圖5 不同飽和溫度下的液相容積分配圖

CO2沸騰換熱系數(shù)隨干度x的變化如圖6所示，隨著飽和溫度不斷升高，CO2工質(zhì)的傳熱系數(shù)也隨之增大，在低干度泡狀流區(qū)域，Tsat為273 K的CO2工質(zhì)沸騰換熱系數(shù)要比Tsat為253 K時(shí)平均高出2.1 kW/（m2·K）左右.沸騰換熱系數(shù)隨飽和溫度的升高而不斷增大，并且在Tsat較低時(shí)，沸騰換熱增長(zhǎng)趨勢(shì)更明顯.

2.2 熱流密度對(duì)傳熱系數(shù)的影響

圖7為在飽和溫度Tsat=263 K，質(zhì)量流速G=300 kg/（m2·s），熱流密度q在7.5，15，30 kW/m2的工況下，CO2工質(zhì)在不同熱流密度下水平管進(jìn)口和出口的液相容積分配圖.由圖7a可知，在熱流密度q為7.5 kW/m2時(shí)，水平圓管入口處的氣泡數(shù)量很少且體積小，到出口位置附近干度在0.1左右時(shí)，氣泡的數(shù)量明顯增多，大氣泡生成頻率增大，流場(chǎng)被加速攪混，沸騰現(xiàn)象加劇.當(dāng)增大熱流密度為30 kW/m2時(shí)，可以觀察到在同一位置入口附近，CO2工質(zhì)所產(chǎn)生的氣泡數(shù)量相比于熱流密度為7.5 kW/m2時(shí)明顯增加，氣泡尺寸也明顯增大，表明CO2工質(zhì)的傳熱系數(shù)有明顯的增加.這證明在低干度區(qū)，CO2工質(zhì)的沸騰換熱系數(shù)受熱流密度影響較大.

圖6 不同飽和溫度下沸騰換熱系數(shù)隨干度變化

圖7 不同熱流密度下的液相容積分配圖

CO2沸騰換熱系數(shù)隨干度x的變化如圖8所示，隨著熱流密度不斷增大，CO2工質(zhì)的沸騰換熱系數(shù)也隨之增大，在低干度泡狀流區(qū)域，熱流密度為30 kW/m2的CO2工質(zhì)沸騰換熱系數(shù)比15 kW/m2和7.5 kW/m2的沸騰換熱系數(shù)分別高出2.5 kW/（m2·K）和8 kW/（m2·K）左右.在較高熱流密度下CO2工質(zhì)的沸騰換熱系數(shù)遠(yuǎn)高于較低熱流密度下的沸騰換熱系數(shù)，這是因?yàn)樵诘透啥葏^(qū)域主要是泡狀流，核態(tài)沸騰占主導(dǎo)，管內(nèi)沸騰換熱加快，沸騰換熱系數(shù)提高.

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

在飽和溫度Tsat=263 K，熱流密度q=7.5 kW/m2，質(zhì)量流速G在300，400，500，600 kg/（m2·s）的工況下，圖9為CO2工質(zhì)在不同質(zhì)量流速下水平管進(jìn)口和出口的液相容積分配圖.由圖9a可知，質(zhì)量流速為300 kg/（m2·s）時(shí)，水平圓管入口處，由于核態(tài)沸騰在初始階段，氣泡數(shù)量很少，幾乎沒有，到出口位置附近干度在0.1左右，氣泡的數(shù)量增多，有很多小氣泡合成大氣泡，隨著干度x的增大，流場(chǎng)擾動(dòng)加大，沸騰換熱提高.當(dāng)質(zhì)量流速增大到400 kg/（m2·s）時(shí)，可以觀察到在同一位置出入口附近，CO2工質(zhì)所產(chǎn)生的氣泡數(shù)量相比于300 kg/（m2·s）時(shí)有所增加，表明CO2工質(zhì)的沸騰換熱系數(shù)也增加，但并不明顯.當(dāng)質(zhì)量流速繼續(xù)增大到600 kg/（m2·s）時(shí)，水平圓管進(jìn)口的小氣泡明顯增多，出口的氣泡多為扁平的大氣泡，在數(shù)量上也多于低質(zhì)量流速的情況.但隨著質(zhì)量流速的增大，一方面氣泡脫離速度加快，換熱加強(qiáng)，CO2工質(zhì)沸騰換熱系數(shù)升高；另一方面氣泡過多會(huì)增加管內(nèi)氣相工質(zhì)與壁面的接觸，使換熱發(fā)生惡化，導(dǎo)致CO2工質(zhì)沸騰換熱系數(shù)急劇下降，兩種因素都存在就導(dǎo)致兩種影響抵消，所以在低干度區(qū)的CO2工質(zhì)沸騰換熱系數(shù)受質(zhì)量流速的影響較小.

圖8 不同熱流密度下傳熱系數(shù)隨干度變化圖

圖9 不同質(zhì)量流速下的液相容積分配圖

CO2沸騰換熱系數(shù)隨干度的變化如圖10所示，隨著質(zhì)量流速不斷升高，CO2工質(zhì)的沸騰換熱系數(shù)沒有明顯的大幅度變化，在低干度泡狀流區(qū)域，質(zhì)量流速為600 kg/（m2·s）的CO2工質(zhì)沸騰換熱系數(shù)比質(zhì)量流速為300 kg/（m2·s）的CO2工質(zhì)沸騰換熱系數(shù)僅高出17.02%.說明與其他的影響因素相比，質(zhì)量流速對(duì)核態(tài)沸騰的影響較小.

圖10 不同質(zhì)量流速下傳熱系數(shù)隨干度變化

3 結(jié)論

通過Fluent軟件模擬CO2在低干度下微細(xì)通道內(nèi)的沸騰換熱特性，結(jié)論如下：

（1）隨著飽和溫度和熱流密度的提高，以及管徑的降低，液相容積分配圖出口位置的氣泡增多且體積變大，說明換熱效果增強(qiáng)，CO2沸騰換熱系數(shù)增大.

（2）隨著質(zhì)量流速的增大，氣泡脫離速度加快，換熱增強(qiáng)，但氣相工質(zhì)與壁面的接觸會(huì)惡化換熱，兩種作用相互抵消，導(dǎo)致液相容積分配圖出口位置氣泡沒有明顯變化，質(zhì)量流速對(duì)CO2沸騰換熱系數(shù)影響較小.也說明了在低干度區(qū)，CO2核態(tài)沸騰機(jī)理占主導(dǎo)地位.