微通道內流動沸騰特性研究

姜圣列　馬虎根　徐法堯

摘 要： 對國內外微通道流動和換熱的研究實驗作了總結，闡述了影響微通道換熱系數的因素，如熱流密度、過熱度和干度等.對去離子水在內徑為0.65 mm、長為102 mm的圓形管道內流動沸騰換熱進行了實驗研究，得到了局部換熱系數隨干度的變化關系，進而根據換熱系數的變化趨勢討論了飽和流動沸騰區微通道內主導的換熱機制.結果表明：從換熱系數隨干度的變化關系很難判定主導的換熱機制；將實驗數據與已發表的預測關聯式進行了比較，發現大多關聯式都失效，說明基于常規理論的模型不再適用于微通道.

關鍵詞：微通道； 流動沸騰換熱； 干度； 換熱關聯式； 模型

中圖分類號： TK 124 文獻標志碼： A

隨著傳統機電系統微型化的發展及其應用領域的不斷擴大，有關流體在微通道中流動和換熱的研究也越來越引起人們的重視，傳統流動和傳熱理論在微尺度下是否依然適用是目前研究的重要課題之一[1].但微細通道內的流動沸騰特點還沒有被清晰地闡明.因此國際上逐步形成了微細尺度傳熱這一新的分支學科.Lazarek等[2]對R113在小通道和大管道的流動沸騰特征作了對比，結果與Kandlikar[3]的關聯式吻合較好，但高干度區存在一定差別.Kamidis等[4]對1.59 mm小管內R113的流動沸騰換熱實驗數據作了整理，并與Kandlikar關聯式進行了比較，結果發現，液相雷諾數Re1=5 720時，兩者吻合較好，但Re1=2 370時，關聯式結果偏低.Wambsganss[5]的實驗顯示出熱流系數對熱流密度和質量流率的依賴性，這表明核態沸騰和對流沸騰機理對換熱有重要影響.

Tran等[6]實驗研究了R12在管徑為2.92 mm圓管內的流動沸騰換熱特性，實驗結果表明，當壁面過熱度大于2.75 K時，局部換熱系數受質量流量（50～695 kg·m-2·s-1）和平衡干度（0.20～0.75）的影響不大，但會隨著熱流密度的變化而改變.根據常規通道內流動沸騰換熱理論，Tran認為微通道內流動沸騰總是核態沸騰換熱模式，并將實驗數據擬合成核態沸騰主導關聯式形式，即q=aΔTn，其中：q為熱流密度；ΔT為過熱度；a為常數；n=2.7.這種形式是典型的池沸騰關聯式.實驗結果表明，當壁面過熱度小于2.75 K，沸騰曲線的斜率發生很大的改變.Tran認為換熱模式轉變為強制對流主導，換熱系數只是隨著熱流密度的變化而改變.

Zhao等[7]的實驗結果表明，在一定熱流密度區間，換熱系數不隨熱流密度變化而變化.而Hwang等[8]的實驗表明，在很寬泛的干度區間熱流密度及質量流量都對換熱系數有很大的影響.

Jatuporn等[9]研究了R134a在水力直徑分別為1.1 mm和1.2 mm微通道內的流動沸騰換熱特性.結果表明，換熱系數變化趨勢與Jatuporn關聯式[10]符合較好.Jatuporn關聯式是在單相強制對流基礎上加上兩相修正因子S得到，在該關聯式基礎上又考慮了熱流密度和質量流量的影響，分別以沸騰數Bo（代表加熱量與蒸發量關系）和韋伯數We（代表慣性力和表面張力效應之比）表示，最終得到換熱系數關聯式h=f（S，Bo，We，h1），其中h1為單相時換熱系數.

Li等[11]認為與沸騰換熱有關的力有：重力、慣性力、黏性力和表面張力，用無量綱數Bond（浮力和表面張力的比值）和雷諾數Re表示.Li從力學角度擬合實驗數據，將換熱系數表示為這些無量綱數的關聯式.此外，考慮到微通道內的沸騰實驗表明了換熱系數與熱流密度和質量流量的關系，引入Bo.最終換熱系數可表示為Bo、Bond和Re的函數，從而建立適用于微通道飽和流動沸騰換熱關聯式.

Thome等[12]提出了微管道內流動沸騰的三區模型.建立該模型的思路為：假定氣泡在微通道內快速成核并生長到通道直徑大小，氣泡的生長空間受到管壁的限制，從而形成一系列拉長的氣泡.與此同時，一薄層液膜在拉長的氣泡和管壁間形成，此液膜將擔負起熱量傳遞的主要任務.周期性產生的氣泡將下游的流動區域劃分為液塞區和拉長氣泡區.隨著傳熱過程的進行，如果拉長氣泡區部分液膜被蒸干還將產生汽塞區.在通道某一橫截面位置液塞、拉長泡狀流和/或汽塞隨時間周期性通過.流體和管壁間的熱量傳遞主要通過液塞的對流換熱、拉長氣泡液膜的導熱和干涸區的蒸汽對流換熱三部分進行.根據一個周期內三種換熱模式在該截面處分別所占據的時間權重，計算一個周期內該截面處的局部時均換熱系數.該模型可定性解釋實驗得到的換熱系數的變化趨勢，定量分析結果也較合理（67%的實驗數據在模型預測值的±30%之內）.

綜上可知，對于微通道流動沸騰還沒有通用的換熱關聯式.本文對去離子水在直徑為0.65 mm不銹鋼圓形截面微通道內流動沸騰換熱進行實驗，得到換熱系數隨干度的變化趨勢，然后將得到的飽和流動沸騰局部換熱系數與已發表的基于常/微通道建立的預測關聯式進行比較.

1 實驗過程

1.1 實驗裝置

設計了標準大氣壓力下去離子水在微細通道內流動和換熱研究實驗裝置，如圖1所示.

該實驗裝置包括兩部分：流體流動與傳熱部分和數據采集部分.實驗運行流程為：儲液箱內的去離子水在恒流泵的驅動下，以某一恒定流量進入預熱段吸收熱量達到某一恒定的出口溫度.之后去離子水進入實驗段，去離子水在實驗段流動過程中吸收由直流加熱電源提供的加熱量，同時Agilent數據采集儀采集由壓力傳感器a、b分別測得的流體進、出口壓力，由熱電阻a、b分別測得的流體進、出口溫度，熱電偶a、b、c、d測得的實驗段壁面溫度，并將采集數據輸入計算機.去離子水通過冷凝裝置后流入燒杯，用電子天平稱得液體的質量，用秒表記錄稱重時間，從而得出去離子水的流量.進、出口溫度采用4線制熱電阻測量，熱電阻Pt100精度為A級；進、出口壓力采用U.S.Setra 209系列壓力傳感器測量，精度為0.25%；壁面溫度采用E型熱電偶測量，標定精度為0.1℃；質量采用FA2004電子天平測量，精度為0.1 mg.微通道內徑為0.65 mm，長度為102 mm.

1.2 實驗步驟

實驗過程中設定泵的轉速為一定值，調節預熱段加熱器使實驗段入口溫度恒定在40℃，然后調節實驗段的直流輸出電源.隨著實驗段電源輸出功率的增加實驗段出口溫度逐漸升高（此時入口溫度一般也會發生小的波動，這主要是由熱量沿微通道管壁和流體的軸向導熱引起的，該部分本文并沒有考慮）.在每個實驗段輸出功率對應的工況下，待流動和傳熱達到穩定狀態約需30 min采集壓力和溫度數據，測量10 min的流量完成一個實驗工況的工作.調節實驗段直流輸出電源的輸出功率，過冷度很大時（本實驗中壁面溫度為70～98℃）一般以出口溫度升高5℃為一個實驗工況.完成這一流量下的實驗工況重新設定泵的轉速按照上述的步驟進行實驗，直至完成實驗.

2 實驗數據處理

過冷段長度Lsp可根據熱平衡計算得到，即式中：M為質量流量；cp為去離子水在進、出口平均壓力時的定壓比熱；Tsat為去離子水在進、出口平均壓力時的飽和溫度；Tfi為去離子水在微通道入口處的溫度；L為實驗段的長度；Qeff為內表面有效加熱量.

微通道內流體溫度Tfz分為兩部分：過冷段流體溫度和飽和段流體溫度.為此，要確定微通道內流體溫度時首先要判定測量點是處于過冷段還是飽和段.

Lz為測量點距離微通道入口的距離.當Lsp>Lz時，測量點處的管道截面流體處于過冷段，則

式中：hi為微通道入口處焓值；hsat為測量點處壓力對應的飽和焓值；hfg為進、出口平均壓力值對應的汽化潛熱；x的最大不確定度為6.8%.

3 實驗結果及分析

實驗過程中維持去離子水進口溫度為65℃，質量流率G=339.1 kg·m-2·s-1.實驗過程中質量流率產生了波動，干度較小時質量流率較大，干度較大時質量流率較小，此處取平均值.在求解局部換熱系數時采用質量流率的實際測量值.逐漸增大熱流密度，研究測量點處（距入口100 mm）局部換熱系數隨干度的變化.在整個實驗工況內各參數變化范圍為：G=296.4～379.5 kg·m-2·s-1；qeff=15.1～452.4 kW·m-2；x=-0.05～0.35.

3.1 微通道內流動沸騰換熱系數

將文獻[6]提出的劃分換熱主導機制的準則與本文實驗結果進行比較.在飽和流動沸騰區（0.02

圖3為測量點處局部換熱系數隨干度的變化關系，圖中顯示：當-0.10

數隨干度增加而增大；當0.02

圖3中的負干度區域（-0.100）.根據常規通道飽和流動沸騰理論，低干度區核態沸騰主導著換熱機制，換熱系數隨著熱流密度增加而增大.隨著熱流密度的繼續增加，管內氣泡聚合形成彈狀和環狀，從而使管內流型由泡狀流演變為環狀流，塞狀流和環狀流的出現有可能抑制氣泡的產生.核態沸騰受到抑制，環狀流區域控制的強制對流蒸發成為主導的換熱機制.圖3中換熱系數變化趨勢與這些常規理論有很大的不同.換熱系數的峰值出現在0.02

3.2 實驗數據與常規關聯式的比較

圖4為實驗數據與文獻[3]和[14]中的常規關聯式的比較，從圖中可看出，這兩種基于常規通道的預測關聯式均

不能很好地預測0

圖5為實驗數據與基于微通道建立的關聯式的比較結果.從圖中可看出，除了Jatuporn[10]和Thome[12]關聯式能基本預測00.15時，換熱系數幾乎不隨干度變化而改變.這是因為在此區域時，根據Thome模型提供的公式得到的本實驗工況參數時的起始液膜厚度δ0，小于最小液膜厚度δmin=0.3×10-6 m.在實驗數據處理過程中取δ0=δmin意味著該工況下拉長氣泡周圍的液膜不再參與換熱過程，這顯然與事實不符.Thome模型中也明確指出該模型對最小液膜厚度和液膜形成過程的強烈依賴.根據本文研究可發現Thome模型中最小液膜厚度或采用液膜形成過程模型存在不合理性.

4 結 論

通過研究去離子水在直徑為0.65 mm的圓形管道內的流動沸騰換熱，得到了局部換熱系數隨干度的變化關系，進而根據換熱系數的變化趨勢討論了飽和流動沸騰區微通道內主導的換熱機制.將實驗數據與已發表的預測關聯式進行比較，主要結論為：

（1） 已發表的飽和流動沸騰區換熱主導機制的判定準則均不能有效地判別本實驗的主導換熱機制.在本實驗中飽和流動沸騰區的整個工況范圍內，從換熱系數隨著干度的變化關系很難判定主導的換熱機制.

（2） 兩種基于常規通道的預測關聯式均不能很好地預測0

參考文獻：

[1] 過增元.國際傳熱研究前沿——微細尺度傳熱[J].力學進展，2000，30（1）：1-6.

[2] LAZAREK G M，BLACK S H.Evaporative heat transfer，pressure drop and critical heat flux in a small vertical tube with R113[J].International Journal of Heat and Mass Transfer，1982，25（7）：945-960.

[3] KANDLIKAR S G.Development of a flow boiling map for subcooled and saturated flow boiling of different fluids in circular tubes[ J].Journal of Heat Transfer，1991，113（1）：190-200.

[4] KAMIDIS D E，RAVIGURAJAN T S.Single and twophase refrigerant flow in minichannels[C].Proceeding of the 33rd National Heat Transfer Conference，Albuquerque，NM，1999：1-8.

[5] WAMBSGANSS M W，FRANCE D M，JENDRZEJCZYK J A，et al.Boiling heat transfer in a horizontal smalldiameter tube[J].Journal of Heat Transfer，1993，115（4）：963-972.

[6] TRAN T N，WAMBSGANSS M W，FRANCE D M.Small circular and rectangular channel boiling with two refrigerants[J].International Journal of Multiphase Flow，1996，22（3）：485-498.

[7] ZHAO Y，MOLKI M，OHADI M M，et al.Flow boiling of CO2 in microchannels[J].ASHRAE Conference，2000，106（1）：437-445.

[8] HWANG Y W，KIM M S，RO S T.Experimental investigation of evaporative heat transfer characteristics in a smalldiameter tube using R134a[C].Proceedings of Symposium on Energy Engineering in the 21st Century，New York：Begell House，2000：965-971.

[9] JATUPORN K O，SAKAMATAPAN K，WONGWISES S.Flow boiling heat transfer of R134a in the multiport minichannel heat exchangers[J].Experimental Thermal and Fluid Science，2011，35（2）：364-374.

[10] JATUPORN K O，WONGWISES S.Experimental investigation of evaporation heat transfer coefcient and pressure of R410A in multiport minichannel[J].International Journal of Refrigeration，2009，32（1）：124-137.

[11] LI W，WU Z.A general correlation for evaporative heat transfer in micro/minichannels[J].International Journal of Heat and Mass Transfer，2010，53（9/10）：1778-1787.

[12] THOME J R，DUPONT V，JACOBI A M.Heat transfer model for evaporation in microchannels.Part Ⅰ：prosentation of the model[J].International Journal of Heat and Mass Transfer，2004，47（14-16）：3375-3385.

[13] WOJTAN L，URSENBACHER T U，THOME J R.Investigation of flow boiling in horizontal tubes：Part Ⅱ：development of a new heat transfer model for stratifiedwavy，dryout and mist flow regimes[J].International Journal of Heat and Mass Transfer，2005，48（14）：2970-2985.

[14] CHEN J C.Correlation for boiling heat transfer to saturated fluids in convective flow[J].Industrial and Engineering Chemistry Process Design and Development，1966，5（3）：322-329.

[15] LIE Y M，SU F Q，LAI R L，et al.Experimental study of evaporation heat transfer characteristics of refrigerants R134a and R407C in horizontal small tubes[J].International Journal of Heat and Mass Transfer，2006，49（1/2）：207-218.

[16] CHOI K I，PAMITRANA A S，OH C Y，et al.Boiling heat transfer of R22，R134a，and CO2 in horizontal smooth minichannels[J].International Journal of Refrigeration，2007，30（8）：1336-1346.

[17] BASU S，NDAO S，MICHNA G J，et al.Flow boiling of R134a in circular microtubes—Part I：study of heat transfer characteristics[J].Journal of Heat Transfer（ASME），2011，133（5）：0515021- 0515029.

[18] WARRIER G R，DHIR V K，MOMODA L A.Heat transfer and pressure drop in narrow rectangular channels[J].Experimental Thermal and Fluid Science，2002，26（1）：53-64.

[19] SUN L，MISHIMA K.An evaluation of prediction methods for saturated flow boiling heat transfer in minichannels[J].International Journal of Heat and Mass Transfer，2009，52（23/24）：5323-5329.

[20] BERTSCH S S，GROLL E A，GARIMELLA S V.A composite heat transfer correlation for saturated ow boiling in small channels[J].International Journal of Heat and Mass Transfer，2009，52（7/8）：2110-2118.

[21] QI S L，ZHANG P，WANG R Z，et al.Flow boiling of liquid nitrogen in microtubes：Part ⅡHeat transfer characteristics and critical heat flux[J].International Journal of Heat and Mass Transfer，2007，50（25/26）：5017-5030.