石綺云,王嘉程,張江輝,劉 麗,項燕龍,何 燕,馬連湘,李 蔚,3*
(1.青島科技大學 機電工程學院,山東 青島 266061;2.寧波金田管業(yè)有限公司,浙江 寧波 315000;3.浙江大學 能源工程學院,浙江 杭州 310027)
螺旋微翅片結(jié)構(gòu)是一種典型的強化換熱管結(jié)構(gòu),它在提供較明顯的換熱提升能力同時將壓降的損失控制在理想可接受的范圍之內(nèi)[1-4]。WU 等[5]實驗研究了R22 和R410A 在5 根具有相同5.0 mm 外徑的單凹槽微翅片管內(nèi)的蒸發(fā)換熱和摩擦壓降特性。他們提出,翅片高度和液膜厚度之間的比值是影響微翅片管內(nèi)換熱特性的關鍵因素。當兩者比值接近于1 的時候,換熱系數(shù)的值更高。文獻[6-7]等分別對內(nèi)徑為2.4和3.48 mm 的小管徑微翅片管進行管內(nèi)流動沸騰換熱實驗的研究。結(jié)果發(fā)現(xiàn),在高熱流密度下,換熱系數(shù)不會被質(zhì)量流速影響,干度的影響同樣也不大,因為此時核態(tài)沸騰占據(jù)了沸騰換熱過程的主要影響因素。而在低干度的區(qū)域,換熱系數(shù)會隨著熱流密度的增加而增加,這是因為核態(tài)沸騰效果會隨著熱流密度的增加而得到增強。
近十年來,有關三維強化管的研究逐漸成熟,并引起廣泛關注。AROONRAT 等[8]同樣對R134a在內(nèi)徑為8.1 mm 的水平凹坑強化管內(nèi)的流動沸騰換熱和壓降特性進行實驗探究。實驗的蒸發(fā)溫度范圍在7~13℃,熱流密度在20~30 k W·m-2,質(zhì)量流速范圍在300~500 kg·(m2·s)-1。實驗結(jié)果證明,具有凹坑的強化換熱管能夠?qū)Q熱系數(shù)提升至光滑管的1.7倍,同時壓降損失大約在光滑管壓降損失的2.78~5.78倍。郭雨等[9]針對不銹鋼三維強化管管內(nèi)流動沸騰和冷凝傳熱進行了實驗研究。結(jié)果表明,由一系列花瓣和凹坑組成的強化換熱管1EHT 表現(xiàn)出最佳的兩相換熱性能。課題組此前探究過對單一結(jié)構(gòu)的強化換熱管外的冷凝換熱特性實驗[10],實驗中所涉及的強化結(jié)構(gòu)為單一的不規(guī)則球形凹坑或規(guī)則的矩形凹坑,而并沒有涉及到針對同時具有2種不同強化結(jié)構(gòu)的復合強化換熱管的相關研究。因此本研究針對這一研究的空白領域,實驗探究了3種具有不同外表面強化結(jié)構(gòu)換熱管的管外環(huán)形側(cè)流動沸騰換熱性能和表現(xiàn)。參考經(jīng)典的流動沸騰換熱模型并對其進行修正,提出適合本研究實驗結(jié)果的管外對流沸騰換熱關聯(lián)式。
如圖1所示為實驗過程中所使用的系統(tǒng)裝置示意圖。整個實驗系統(tǒng)由制冷劑循環(huán)、預熱段循環(huán)回路、實驗段循環(huán)回路以及過冷段循環(huán)回路組成。
圖1 實驗系統(tǒng)裝置的示意圖Fig.1 Schematic diagram of the experimental setup
制冷劑由儲液罐中流出,依次經(jīng)過齒輪泵和流量計流進預熱段裝置,在預熱段與預熱段的循環(huán)水進行換熱被加熱到合適的工況,接著流進套管式逆流換熱器組成的實驗段。制冷劑從管外的環(huán)形側(cè)流過,與管內(nèi)的去離子水進行熱交換,以達到出口所需工況。此后,制冷劑經(jīng)過冷段裝置被冷卻到一定的過冷度之后重新回到儲液罐,等待開啟下一次循環(huán)。
實驗開始之后,低溫低壓的制冷劑通過干燥設備經(jīng)由齒輪泵流過科式流量計。經(jīng)過預熱段的熱量吸收過程,通過位于試驗段進口處的壓力和溫度傳感器測量溫度和壓力值,并將數(shù)據(jù)記錄進數(shù)據(jù)采集系統(tǒng)。繼而流過實驗段完成換熱過程,從實驗段出口流出。實驗段出口處同樣設置高精度壓力和溫度傳感器對工質(zhì)的溫度、壓力等參數(shù)進行采集。整個實驗過程的數(shù)據(jù)采集由具有20通道A/D 采集卡的Agilent34970A 數(shù)據(jù)采集儀完成,并將獲得的實驗數(shù)據(jù)實時上傳至上位機檢測系統(tǒng)中。通常每20 s對實驗數(shù)據(jù)進行1次采集,當上位機檢查窗口在連續(xù)的20 min內(nèi)的溫度波動值小于0.1 K,壓力波動值小于2 kPa,流量波動值小于0.5 kg·h-1,干度波動值小于0.02的時候,認定工況點達到了穩(wěn)定狀態(tài)。選取連續(xù)的10個數(shù)據(jù)點進行數(shù)據(jù)記錄。為進一步減小實驗數(shù)據(jù)的誤差值,對各個參數(shù)取平均值后作為一個工況點并記錄。
本研究的實驗對象為3種具有不同表面結(jié)構(gòu)的強化換熱管,分別為HB、DIM 和DIM/HB 管。3根強化換熱管具有相同的外徑12.7 mm(do),外套管的內(nèi)徑是17 mm(Di)。管外側(cè)走制冷劑工質(zhì),管內(nèi)側(cè)工質(zhì)為去離子水。實驗被測試管的具體結(jié)構(gòu)參數(shù)如圖2所示。
圖2 強化管細節(jié)展示圖Fig.2 Details of the enhanced tubes
被測換熱管均為不銹鋼材質(zhì),公稱外徑12.7 mm,公稱內(nèi)徑11.5 mm,測試管長度2 m。3根強化換熱管的表面結(jié)構(gòu)具體參數(shù)如表1所示,其中DIM/HB管的表面結(jié)構(gòu)為HB和DIM 2種結(jié)構(gòu)的復合。
表1 強化換熱管的結(jié)構(gòu)參數(shù)Table 1 Geometric parameters of tested tubes
管外流動沸騰換熱實驗在飽和溫度為6℃下進行,所使用的工質(zhì)為R410A。質(zhì)量流速的范圍從75~225 kg·(m2·s)-1,進出口干度值維持在0.2/0.8。在整個實驗過程中,管內(nèi)側(cè)水的進出口平均溫度范圍是278~292 K,進出口的溫差保持在4~10 K,流速范圍是200~305 kg·h-1。
為了驗證實驗臺準確性,進行單相熱平衡實驗。實驗工況選擇飽和溫度為311 K,制冷劑側(cè)質(zhì)量流速范圍在50~125 kg·(m2·s)-1。實驗過程中測試段管內(nèi)水側(cè)的換熱量
其中:Qw為水側(cè)換熱總熱量,cp,w為水側(cè)進出口的平均溫度下水的定壓比熱容,mw為水側(cè)質(zhì)量流速,Tw,in和Tw,out分別為測試段水側(cè)進口溫度和出口溫度。管外環(huán)形側(cè)的制冷劑換熱量
其中:Qre是測試段制冷劑側(cè)的總換熱量,hre,in和hre,out分別為管外測試段制冷劑的進、出口焓值。焓值的獲取可以通過測量制冷劑進出口側(cè)的溫度Tte,re,in、Tte,re,out和測試段制冷劑側(cè)進出口的壓力值pte,re,in、pte,re,out,再結(jié)合物性查詢軟件REFPROP 9.1[24]查詢得到。
單相熱平衡實驗結(jié)果如圖3所示。由圖3可以看出,不論是光管或是具有表面強化結(jié)構(gòu)的換熱管,其熱量損失都維持在5%以內(nèi)。由此可以推斷本實驗臺的熱密封性能良好,實驗結(jié)果真實可靠。
圖3 單相熱平衡實驗結(jié)果Fig.3 Single-phase heat balance of tested tubes
實驗過程的對數(shù)平均溫差如下:
在不考慮管壁熱阻的情況下,制冷劑側(cè)的換熱系數(shù)
其中:Ao為制冷劑側(cè)換熱面積,Ai為水側(cè)換熱面積,λwall為管壁材料的導熱系數(shù),l為測試管長度,do和di分別為測試管的公稱外徑和公稱內(nèi)徑。
由文獻[11]可查證,Gnielinski公式作為預測光滑管內(nèi)水側(cè)單相換熱系數(shù)的關聯(lián)式具有相對較高的準確性。該公式的適用范圍是:0.5≤Pr≤2 000,并且3 000≤Re≤5×106。
其中:(μbulk/μwall)0.14項是物性參數(shù)比修正項,以減小壁面溫度對近壁面處制冷劑黏度的影響,μbulk和μwall 分別是水側(cè)流體的中心溫度和壁面溫度所對應的液態(tài)黏度,λw為水側(cè)進出口平均溫度下所對應的導熱系數(shù),dh是管內(nèi)側(cè)水力直徑。關聯(lián)式中的fw為范寧摩擦壓降系數(shù),由Petukhov 關聯(lián)式[12]得出
為了修正水側(cè)的換熱系數(shù)預測關聯(lián)式預測結(jié)果與實際工況下?lián)Q熱系數(shù)值之間的偏差,引入一個強化倍率的概念對預測的水側(cè)換熱系數(shù)值進行調(diào)整。因此,實驗的強化換熱管內(nèi)的換熱系數(shù)計算公式被修正為
其中:C為強化倍率。
強化倍率的計算本研究采用威爾遜圖解法[12]。威爾遜圖解法的原理在于,當管內(nèi)側(cè)水的流速以等間距增加并增大到一定值,總的換熱熱阻值會逐漸趨近于另外2個換熱熱阻值之和。此時,總換熱熱阻值減去管壁熱阻即得到管外制冷劑側(cè)的換熱熱阻值:
實驗得到的3根不同強化換熱管的強化倍率如圖4所示。
圖4 威爾遜圖解法Fig.4 Wilson plot method
采用威爾遜圖解法[13]計算得到ST、HB、DIM和DIM/HB管的管內(nèi)強化倍率C的值分別:1.06、1.01、1.74和1.77。其中光滑管ST 和HB 管的強化倍率所得結(jié)果在1.1以內(nèi),證明實驗結(jié)果與關聯(lián)式預測的偏差范圍小于10%。由此可以看出Gnielinski公式[11]適用于本研究實驗工況,進一步驗證實驗系統(tǒng)的準確性。
目前實驗參數(shù)誤差的計算方法采用誤差傳遞算法。誤差傳遞算法[14]的基本計算原理
其中:y是因變量,r(y)是因變量的相對誤差,x i是自變量,σ2(x i)是自變量的標準偏差值。計算結(jié)果顯示,管外流動沸騰換熱系數(shù)的不確定度計算值結(jié)果為12.31%。
其余各參數(shù)值的誤差計算結(jié)果如表2、3所示。
表2 直接測量參量誤差值Table 2 Uncertainties of measuring parameters
表3 間接計算參量誤差值Table 3 Uncertainties of calculation parameters
注:壓力范圍:0~40 bar;壓降范圍:0~100 k Pa;水的質(zhì)量流量范圍:0~12 L·min-1;制冷劑的質(zhì)量流量范圍:0~60 kg·h-1。
3根不同的強化換熱管在實驗過程中所展示出的蒸發(fā)換熱特性如圖5(a)、(b)所示。
圖5 蒸發(fā)換熱特性隨質(zhì)量流速變化的特性Fig.5 Evaporation heat transfer characteristics as a function of mass flux
圖5(a)展示了3種強化換熱管在實驗工況下的換熱系數(shù)隨著質(zhì)量流速的變化情況。從圖5中可以看出,換熱系數(shù)均隨質(zhì)量流速的增加而上升,且DIM 管的沸騰換熱系數(shù)高于HB管。3種強化換熱管中,DIM/HB管的沸騰換熱系數(shù)最高,為光管的1.54~2.23倍。圖5(b)展示了4種不同換熱管管外沸騰換熱壓降實驗結(jié)果??梢钥吹?HB 管的管外壓降比光滑管的管外壓降平均增加了9%,DIM管的管外壓降平均值比光滑管提高了5%,DIM/HB 管的管外壓降平均值是光滑管的1.26 倍。DIM 管相比較HB 管具有更高的換熱系數(shù)和更低的壓降損失,復合表面強化換熱管DIM/HB管在本次實驗中也表現(xiàn)出了較高的壓降損失。
針對實驗得到的光滑管管外蒸發(fā)換熱系數(shù)實驗數(shù)據(jù),選取GUNGOR 和WINTERTON[15]、LIU 和WINTERTON[16]、KATTAN 等[17],以 及KANDLIKAR等[18]所提出的關聯(lián)式進行預測對比。其中文獻[15-16,18]關聯(lián)式的建立都是在CHEN[19]所提出的關于飽和狀態(tài)流體對流沸騰換熱模型的基礎之上。文獻[15]則對關聯(lián)式中的對流沸騰強化因子E和核態(tài)沸騰抑制因子S分別修正。而文獻[17]在計算方程中將熱壁面劃分成了2部分:一部分是被工質(zhì)潤濕的周長,稱為潤濕周長;另一部分是未被工質(zhì)潤濕的周長,稱為干涸周長。
圖6展示了采用3種不同水平光滑管內(nèi)沸騰換熱系數(shù)的關聯(lián)式計算模型所得到的預測值和實際值的比較情況。
從圖6可以看到,光滑管管外蒸發(fā)換熱的數(shù)據(jù)可以被上述的4種不同的蒸發(fā)換熱關聯(lián)式模型很好地預測。但質(zhì)量流速在75 kg·(m2·s)-1時,文獻[15]關聯(lián)式的預測誤差在20% 以上,達到了34.7%。此外,文獻[18]關聯(lián)式的預測誤差最小最為精準,文獻[15]關聯(lián)式的預測值整體偏高,而文獻[16]關聯(lián)式的預測值整體偏低。
本研究針對3根強化換熱管的管外蒸發(fā)數(shù)據(jù)的預測分別在文獻[18]和文獻[20]關聯(lián)式基礎上提出2種新的修正預測模型。
早在1990年,KANDLIKAR 在文獻[18]中提出了基于CHEN[19]關于飽和流體水平方向流動沸騰換熱加和機理的光滑管內(nèi)換熱系數(shù)預測關聯(lián)式:
其中,換熱系數(shù)hl表示當制冷工質(zhì)為純液態(tài)時所產(chǎn)生的換熱系數(shù)預測值,其計算采用經(jīng)典的Dittus-Boelter[21]關聯(lián)式:
本實驗所對應的工況為水平流動且Frl>0.04,因此相對應的C5*=0。此外對于近共沸混合制冷劑R410A,關聯(lián)式中的Ffl=1。本工作實驗數(shù)據(jù)點全部符合原文獻[21]中所規(guī)定的對流換熱區(qū)域情況,因此式(2)中 的C1、C2、C3、C4、C5系數(shù)值分別為1.136、-0.9、667.2、0.7、0。
從圖6中可以看到,Kandlikar[18]關聯(lián)式在預測光滑管外換熱系數(shù)時精度較高,誤差在5%以內(nèi)。因此將該關聯(lián)式拓展應用到管外流動沸騰換熱工況下,以達到預期的擬合效果。Kandlikar為研究微尺度下的流動沸騰換熱,分別提出了K1、K2、K33個無量綱量[22]。在新型強化換熱管中,起到強化換熱的結(jié)構(gòu)分別是微翅片和凹坑結(jié)構(gòu)。分析管外流動沸騰的強化機理如下:一是微翅片和凹坑結(jié)構(gòu)的存在會使得對流沸騰過程中的湍流度增強,二是流動沸騰過程中,較薄的液膜會在微翅片及凹坑結(jié)構(gòu)的影響下表面張力發(fā)生變化,從而對核態(tài)沸騰產(chǎn)生影響。因此在對流沸騰項添加一個修飾因子Ecv(方程(12)),在核態(tài)沸騰項中添加另一個修飾因子Ebo(方程(13))。通過f/D的形式從幾何結(jié)構(gòu)參數(shù)的角度進一步量化湍流度的增強效應。
需要注意的是,方程(14)中的f在HB管中表示翅片高度,在DIM 管中表示凹坑深度,而在DIM/HB管中則采用等效值
添加修正因子之后的關聯(lián)式最終形式為
使用MATLAB軟件對實驗數(shù)據(jù)進行回歸分析得出參考系數(shù)a1、a2、b1、b2分別如表4所示。
表4 參考系數(shù)值Table 4 Reference values of coefficients
關聯(lián)式修正之前的預測結(jié)果與實驗數(shù)據(jù)對比結(jié)果如圖7所示。可以看出原始關聯(lián)式在預測強化結(jié)構(gòu)換熱管外的流動沸騰換熱系數(shù)時候整體的預測值都是偏低的。且預測值大部分與實驗值之間的誤差范圍在15%~50%,在某些工況下誤差范圍甚至達到50%以上。在原始關聯(lián)式的基礎之上添加了修飾因子之后,所得到的換熱系數(shù)關聯(lián)式預測值與實際值之間的誤差除個別數(shù)據(jù)點之外,絕大部分都保持在±15%以內(nèi),對比結(jié)果如圖8所示。由此可見所提出的新關聯(lián)式針對本研究實驗數(shù)據(jù)擬合結(jié)果較為精準。
圖7 原始關聯(lián)式預測值與實驗值的比較Fig.7 Comparison between the predicted values of the original correlation and the experimental values
圖8 新關聯(lián)式預測值與實驗值的比較Fig.8 Comparison between the predicted values of the new correlation and the experimental values
THOME[19]對微翅片管內(nèi)換熱系數(shù)的預測提出一個通用的換熱模型,關聯(lián)式中涉及到2個關鍵的強化因子分別是:ERB和Emf。2個因子分別是針對微翅片結(jié)構(gòu)本身和膜狀流動狀態(tài)所提出的修正因子。因此推測該關聯(lián)式同樣適用于本研究中所涉及的管外強化結(jié)構(gòu)及對流沸騰換熱狀態(tài)。針對強化結(jié)構(gòu)換熱管的不同結(jié)構(gòu)特性和流動狀態(tài),另提出修正之后的關聯(lián)式見表5。
表5 修正之后的關聯(lián)式方程Table 5 Modified correlations for enhanced tubes
在使用修正之后的新關聯(lián)式進行換熱系數(shù)預測之前,先利用原始關聯(lián)式對本課題的實驗數(shù)據(jù)進行預測。預測所得的結(jié)果如圖9所示。
圖9 原始Thome關聯(lián)式預測值與實驗數(shù)據(jù)的比較情況Fig.9 Comparison of the test results and the original correlations
從圖9 中可以看出,無論是HB、DIM 還是DIM/HB管,所有的預測值與實驗結(jié)果相比都偏低。此外,圖9中還標注了原始Thome[20]關聯(lián)式在預測3種強化換熱管外環(huán)形側(cè)的流動沸騰換熱系數(shù)的誤差范圍。對于HB、DIM、DIM/HB 管,預測值誤差范圍的最大值分別是36%、45%和29%。
修正后的新關聯(lián)式與本實驗結(jié)果之間的擬合情況如圖10所示。由圖10看到,修正之后的誤差極小。這是由于本研究所選取的數(shù)據(jù)點總量較小,因此針對本研究中強化換熱管外蒸發(fā)換熱情況進行擬合所得到的關聯(lián)式預測精度相對較高。此外,根據(jù)本研究對于管外復合結(jié)構(gòu)蒸發(fā)換熱關聯(lián)式中強化換熱因子的擬合過程可以尋找到其與各組成結(jié)構(gòu)的強化換熱因子之間的聯(lián)系。
圖10 修正后的關聯(lián)式與實驗數(shù)據(jù)的比較情況Fig.10 Comparison of the test results and the modified correlations
圖10展示了基于Thome[20]的修正關聯(lián)式最終的預測結(jié)果與實驗數(shù)據(jù)之間的誤差可以達到5%以內(nèi),相較于上一小節(jié)探討的基于Kandlikar[17]的修正預測模型,其準確度有所提升。由此可以看出,在本實驗條件下,基于Thome[20]的修正關聯(lián)式更適用于對強化換熱管換熱系數(shù)的預測。
采用蒸發(fā)/冷凝高精度實驗臺對3種具有不同管外強化結(jié)構(gòu)的換熱管和光滑管進行了環(huán)形側(cè)的流動沸騰實驗研究。所得出的結(jié)論如下:
1) 同光滑管相比,HB 管、DIM 管和DIM/HB管均在不同程度上提高了管外環(huán)形側(cè)流動沸騰的換熱系數(shù)。其中,DIM/HB管對流動沸騰換熱系數(shù)的提升效果最為明顯,是光管的1.54~2.23倍。由此可見具有凹坑與微翅片的復合結(jié)構(gòu)強化換熱管在管外側(cè)流動沸騰換熱工況下的表現(xiàn)明顯優(yōu)于光滑管以及僅具有單一結(jié)構(gòu)的強化換熱管。
2) 針對3根強化換熱管,分別在Kandlikar[18]和Thome[20]關聯(lián)式的基礎之上進行一定程度的修正。修正之后所得到的兩種不同形式的新關聯(lián)式與原有的關聯(lián)式相比較,預測精度有了較大的提升。其中,基于Thome[20]關聯(lián)式的修正模型預測精度更高。需要注意的是,由于本研究所提出的相關換熱系數(shù)預測關聯(lián)式的數(shù)據(jù)基礎為文章實驗所取得的數(shù)據(jù)點。故其適用工況條件為近共沸混合制冷劑R410A 在飽和溫度6℃、質(zhì)量流速為75~225 kg·(m2·s)-1下的水平管外環(huán)形側(cè)的流動沸騰換熱。