劉永坤，陶于兵，唐宗斌

（西安交通大學能源與動力工程學院，熱流科學與工程教育部重點實驗室，陜西 西安 710049）

①相變蓄熱技術為促進節(jié)能減排、實現(xiàn)能量的合理利用提供了有效途徑。但由于受相變材料導熱系數(shù)低的影響，相變蓄熱系統(tǒng)的儲、放熱速率也較低，因而，急需開展相變蓄熱系統(tǒng)的性能強化研究。Adine等[1]對采用多級相變材料的相變蓄熱方案進行了數(shù)值研究，證明了多級相變材料可以強化相變蓄熱性能。Zhao和Wu[2-3]對在相變材料內(nèi)添加金屬泡沫和膨脹石墨來強化其傳熱性能的方案進行了實驗研究，結(jié)果表明，儲、放熱速率都得到了提升。Yang等[4]研究了應用金屬泡沫和翅片的相變蓄熱器，結(jié)果表明，傳熱效果得到提高，儲、放熱時間減少。Lacroix[5]數(shù)值研究了采用強化傳熱翅片的套管式相變蓄熱單元的蓄熱性能，結(jié)果表明，環(huán)形翅片對中等質(zhì)量流速、小進口溫度的情況最有效。Castell等[6]實驗研究了相變材料側(cè)采用縱向翅片對相變材料放熱過程的強化效果，發(fā)現(xiàn)采用縱向翅片可以有效降低相變材料的放熱時間。Tao等[7]數(shù)值研究了帶有內(nèi)翅片的強化傳熱管用于強化套管式相變蓄熱單元性能的可行性，證明了強化傳熱流體側(cè)的傳熱性能也是提高總體蓄熱性能的一個有效途徑。Sohif和Abduljalil等[8-9]建立了二維數(shù)值模型，對帶有內(nèi)外翅片的三螺旋管式相變蓄熱器中的熔化和凝固過程進行了模擬，研究了不同結(jié)構(gòu)參數(shù)對熔化、凝固過程的影響。

然而，現(xiàn)有的關于相變蓄熱過程的數(shù)值研究大多采用二維簡化模型，忽略了液態(tài)相變材料的自然對流，或者采用有效導熱系數(shù)方法取代自然對流。為了分析自然對流對相變儲熱單元綜合性能的影響，本文建立了考慮自然對流的三維相變蓄熱模型，開發(fā)了三維柱坐標下的數(shù)值模擬程序，探討了自然對流對相變蓄熱過程的影響；并設計了采用外縱向翅片的外強化傳熱管以及帶有內(nèi)螺紋翅片和外縱向翅片的雙側(cè)強化傳熱管的殼管式相變蓄熱單元，數(shù)值研究了強化傳熱管對蓄熱性能的強化效果，為相變蓄熱單元的設計及性能強化提供參考。翅片上底寬0.2 mm、下底寬0.48 mm。

圖1 外強化傳熱管結(jié)構(gòu)圖Fig.1 Schematic for the outside enhanced tube

1 物理及數(shù)學模型

1.1 物理模型

首先針對相變材料導熱系數(shù)小、儲熱速率低的缺點，在套管式蓄熱單元的內(nèi)管外壁面上設置沿軸線方向的翅片，形成外強化傳熱管（圖1），一定程度上解決環(huán)形空間內(nèi)相變蓄熱材料導熱系數(shù)低的缺陷。然后在外強化傳熱管的基礎上，在內(nèi)管內(nèi)壁面采用內(nèi)螺紋翅片結(jié)構(gòu)來進一步強化傳熱流體側(cè)的對流傳熱性能（圖2）。蓄熱單元的總體幾何參數(shù)及相變材料熱物性參數(shù)見表1，其中相變材料采用NaCl-Na2CO3-NaOH 的混合熔鹽（質(zhì)量比為 7.8︰6.4︰85.8）。對于強化傳熱管，外翅片選取鋁型材，沿圓周方向設置12個翅片、翅高9 mm、翅寬2 mm；內(nèi)螺紋翅片高與內(nèi)管內(nèi)徑之比為0.0244、翅片片間距與翅高比為2.577、翅片數(shù)為45個、螺旋角為48°、

圖2 雙側(cè)強化傳熱管結(jié)構(gòu)圖Fig. 2 Schematic for the bilateral enhanced tube

表1 蓄熱單元幾何參數(shù)及相變材料的熱物性參數(shù)Table 1 Geometric parameters for the LHS tube and thermophysical parameters for PCM

傳熱流體采用導熱油Syltherm800，物性擬合公式見式(1)～(4)[10]

為了簡化模型，假設：① 綜合考慮計算精度和計算時間，傳熱流體采用一維分布參數(shù)模型，傳熱流體軸向?qū)岷宛ば院纳⒑雎圆挥?；?相變材料熱物性取常數(shù)；③ 相變蓄熱單元內(nèi)管壁熱阻忽略不計，外管表面當絕熱處理；④ 相變材料的密度符合Boussinesq假設。

1.2 控制方程

依據(jù)上述假設，傳熱流體側(cè)采用一維分布參數(shù)法模型，控制方程如下

相變材料側(cè)為了考慮自然對流的影響，采用三維柱坐標模型，動量方程見式(6)～(8)

能量方程為

1.3 初始條件及邊界條件

由于對稱性選取整個蓄熱管的半個圓周作為研究對象，初始時刻相變材料及傳熱流體的溫度為150 ℃，邊界條件設置如下。

（1）傳熱流體側(cè) 采用給定入口溫度和流速的第一類邊界條件，入口溫度為300 ℃，入口質(zhì)量流量為0.0125 kg/s。

（2）相變材料側(cè)ω方向上，采用對稱邊界條件。r方向上，速度邊界在內(nèi)外表面采用固壁邊界條件；溫度邊界在外表面采用絕熱邊界，內(nèi)表面采用第三類對流傳熱邊界，其中對流換熱系數(shù)采用關聯(lián)式進行計算[7,11]。z方向上，在兩端面速度邊界采用固壁邊界條件；溫度邊界采用絕熱邊界。

2 模型驗證

本文根據(jù)建立的物理及數(shù)學模型，基于SIMPLER算法，自行開發(fā)了Fortran計算程序。經(jīng)過網(wǎng)格獨立性及時間步長驗證后，選取 40×40×40的網(wǎng)格作為計算網(wǎng)格，時間步長設定為5 s。收斂準則取為傳熱流體和蓄熱材料的最大相對熱平衡偏差小于10–3。然后，通過與文獻[5]中的實驗結(jié)果對比，來驗證模型及模擬程序的正確性。驗證中，利用本文的計算模型，采用文獻中的幾何結(jié)構(gòu)及熱物性參數(shù)，在相同工況下計算了傳熱流體相對入口溫度為20 ℃（即比相變材料的熔化溫度高20 ℃）時，蓄熱單元在r=0.0089 m、z=0.95 m以及r=0.0099 m、z=0.51 m兩位置溫度T1、T2隨時間的變化關系，由于本模型采用的是三維模型，T1、T2的溫度采用的是對應的兩圓周上的平均溫度。圖 3為模擬結(jié)果與實驗結(jié)果及文獻結(jié)果的對比圖，由圖可見，模擬結(jié)果和實驗數(shù)據(jù)吻合很好，驗證了模型的可靠性。

圖3 模型驗證Fig.3 Model validation

3 結(jié)果及分析

3.1 強化傳熱管對相變材料熔化時間的影響

圖4是采用光管、外強化傳熱管和雙側(cè)強化傳熱管三種結(jié)構(gòu)時，相變材料熔化分數(shù)隨時間的變化曲線。采用外強化傳熱管時，與光管結(jié)構(gòu)相比，由于強化了相變材料側(cè)的導熱系數(shù)，使得相變材料先開始熔化，而且相變材料的熔化速率高于采用光管時的熔化速率，因而相變材料完全熔化時間提前。采用雙側(cè)強化傳熱管結(jié)構(gòu)時，不僅強化了相變材料側(cè)的導熱系數(shù)，傳熱流體側(cè)的對流傳熱系數(shù)也大幅度增加，因而相變材料開始熔化和完全熔化的時間大幅度提前。采用光管，200 min時相變材料完全熔化；采用外強化傳熱管，164 min時相變材料完全熔化，熔化時間縮短 18.0%；而采用雙側(cè)強化傳熱管，95 min時相變材料已經(jīng)完全熔化，熔化時間縮短 52.5%?？梢姡啾扔诓捎猛鈴娀瘋鳠峁?，采用雙側(cè)強化傳熱管可以明顯提升相變材料的熔化速率，縮短熔化時間。

圖4 強化傳熱管對熔化分數(shù)的影響Fig.4 Effects of enhanced tubes on melting fraction

3.2 強化傳熱管對蓄熱速率及蓄熱量的影響

圖 5顯示了不同傳熱管結(jié)構(gòu)對儲熱速率的影響。從圖中可以看出，熔化前的顯熱蓄熱階段，由于采用外強化傳熱管的蓄熱單元相變材料側(cè)的相對導熱系數(shù)增大，因此，其儲熱速率高于光管結(jié)構(gòu)。但隨著時間推移，采用外強化傳熱管時，相變材料的溫度升高更快，因此，相變材料和傳熱流體的傳熱溫差逐漸小于光管時的傳熱溫差，造成采用外強化傳熱管方案時，儲熱速率反而低于光管。相變蓄熱階段由于外強化換熱結(jié)構(gòu)和光管結(jié)構(gòu)的傳熱溫差相同，相對導熱系數(shù)大的外強化傳熱管的儲熱速率略高。相變蓄熱完成后的顯熱蓄熱階段，與熔化前的顯熱蓄熱階段相似，隨著時間的進行，采用外強化傳熱管時相變材料的溫度比采用光管時的溫度更高，傳熱溫差減小，儲熱速率降低，所以完全熔化后外強化傳熱管蓄熱單元的蓄熱速率小于光管結(jié)構(gòu)。采用雙側(cè)強化傳熱管的蓄熱單元，一方面相變材料側(cè)的導熱系數(shù)增加；另一方面?zhèn)鳠崃黧w側(cè)的對流傳熱系數(shù)增加，因而，在顯熱階段前期的儲熱速率遠高于光管結(jié)構(gòu)。相變蓄熱階段，由于傳熱溫差相同，雙側(cè)強化傳熱管具有更高的傳熱系數(shù)，因此，其儲熱速率明顯高于光管結(jié)構(gòu)。而隨著時間的進行，雙側(cè)強化傳熱管結(jié)構(gòu)先完成相變蓄熱，進入顯熱蓄熱階段，此時采用雙側(cè)強化傳熱管的蓄熱單元的傳熱溫差進一步減小，儲熱速率也再次低于光管結(jié)構(gòu)。

圖5 強化傳熱管對儲熱速率的影響Fig.5 Effects of enhanced tubes on LHS rate

圖 6顯示了不同傳熱管結(jié)構(gòu)對總蓄熱量的影響。在熔化前的顯熱蓄熱階段，蓄熱過程剛開始時，采用外強化傳熱管的蓄熱單元的儲熱速率高于光管，這段時間內(nèi)采用外強化傳熱管的蓄熱單元的總蓄熱量增長較快。而隨著時間推移，外強化傳熱管儲熱速率低于光管，總蓄熱量增長速度減小，光管結(jié)構(gòu)的總蓄熱量逐漸大于外強化傳熱管結(jié)構(gòu)。進入相變蓄熱階段后，外強化傳熱管的蓄熱單元的儲熱速率略高于光管結(jié)構(gòu)，因此，一段時間后，總蓄熱量相比于光管也在一定程度上增大。采用光管時，相變蓄熱過程在281 min結(jié)束，總蓄熱量達到最大值；采用外強化傳熱管時，蓄熱過程在216 min結(jié)束，蓄熱時間縮短了23%。采用雙側(cè)強化傳熱管的蓄熱單元，相變材料完全熔化前的時間內(nèi)，儲熱速率均高于光管結(jié)構(gòu)，這個時間段內(nèi)采用雙側(cè)強化傳熱管的蓄熱單元的總蓄熱量也明顯高于光管結(jié)構(gòu)。完全熔化后進入顯熱蓄熱階段，隨著蓄熱過程的進行，總蓄熱量在146 min達到最大值，蓄熱過程結(jié)束，相比于光管蓄熱時間縮短了48%。另一方面，與光管結(jié)構(gòu)相比，外強化傳熱管和雙側(cè)強化傳熱管結(jié)構(gòu)，由于外翅片占據(jù)了一定空間，相變材料的質(zhì)量減小，最大總蓄熱量小于光管結(jié)構(gòu)。

圖6 強化傳熱管對總蓄熱量的影響Fig.6 Effects of enhanced tubes on thermal energy storage capacity

3.3 強化傳熱管對固–液界面分布的影響

圖7是總?cè)刍謹?shù)f=0.5時，不同軸向位置(z)截面上熔化分數(shù)場圖（固–液界面分布圖）。通過圖7(a)～(c)的綜合對比可以看出，總體上沿著流動方向，傳熱流體的溫度逐漸降低，傳熱溫差逐漸減小，因此，蓄熱速率減慢，引起截面內(nèi)的熔化分數(shù)逐漸減小。同時，由于受自然對流的影響，底部的液態(tài)相變材料在自然對流的驅(qū)動下，會向頂部流動，引起靠近頂部的相變材料的熔化速率明顯高于底部相變材料的熔化速率，從而固–液界面分布呈現(xiàn)明顯的非均勻性，頂部的相變材料完全熔化時，底部相變材料絕大部分還處于固體狀態(tài)。而采用強化傳熱管時，一方面由于傳熱效果增強，沿著流動方向，傳熱流體的溫度降低得更快，因此，沿著流動方向截面內(nèi)的熔化分數(shù)減小得更明顯，特別是采用雙側(cè)強化傳熱管時，入口截面附近相變材料已接近完全熔化，而出口界面附件相變材料幾乎尚未開始熔化。另一外翅片的存在會抑制自然對流現(xiàn)象，使得自然對流引起的固–液界面的非均勻性減弱。

可見，自然對流的存在雖然可以強化相變材料側(cè)的傳熱性能，但同時也會引起嚴重的固–液界面分布不均勻現(xiàn)象；采用相變材料側(cè)加翅片的強化方案，可以有效削弱自然對流引起的固–液界面分布不均勻現(xiàn)象，并起到強化蓄熱速率的作用；采用雙側(cè)強化傳熱管可以進一步強化蓄熱速率。

圖7 總?cè)刍謹?shù)f=0.5時，z=0.1、z=0.5、z=0.9截面內(nèi)的熔化分數(shù)場圖Fig. 7 PCM melting interface in the cross-section of z =0.1, z =0.5, z =0.9 at f=0.5

4 結(jié) 論

本文建立了考慮自然對流的三維相變蓄熱單元的計算模型，提出了采用縱向翅片的外強化傳熱管和同時采用外縱向翅片和內(nèi)螺紋翅片的雙側(cè)強化傳熱管兩種強化結(jié)構(gòu)，以此來強化蓄熱單元的蓄熱性能，并數(shù)值研究了兩種強化方案對相變蓄熱性能的前后效果。研究結(jié)果表明：

（1）采用縱向翅片的外強化傳熱管結(jié)構(gòu)可以有效提高相變材料區(qū)域的有效導熱系數(shù)，一定程度上提高儲熱速率，縮短相變蓄熱時間，但效果較弱。

（2）采用外壁加縱向翅片、內(nèi)壁加螺紋翅片的雙側(cè)強化傳熱管結(jié)構(gòu)，可以同時提高相變材料側(cè)的有效導熱系數(shù)和傳熱流體側(cè)的對流傳熱系數(shù)，相變蓄熱速率得到大幅提升，蓄熱時間明顯縮短。

（3）液態(tài)相變材料的自然對流雖然可以強化相變材料側(cè)的傳熱性能，但同時也會引起固–液界面分布不均勻現(xiàn)象，采用相變材料側(cè)添加翅片的強化方案，不僅可以有效削弱這一現(xiàn)象，而且可以強化相變蓄熱單元的蓄熱性能。

[1]Adine H A，Qarnia H E. Numerical analysis of the thermal behaviour of a shell-and-tube heat storage unit using phase change materials[J].Applied Mathematical Modelling，2009，33（4）：2134-2144.

[2]Zhao C Y，Wu Z G. Heat transfer enhancement of high temperature thermal energy storage using metal foams and expanded graphite[J].Solar Energy Materials & Solar Cells，2011，95（2）：636-643.

[3]Wu Z G，Zhao C Y. Experimental investigations of porous materials in high temperature thermal energy storage systems[J].Solar Energy，2011，85（7）：1371-1380.

[4]Yang J L，Du X Z，Yang L J，Yang Y P. Numerical analysis on the thermal behavior of high temperature latent heat thermal energy storage system[J].Solar Energy，2013，98，PartC：543-552.

[5]Lacroix M. Study of the heat transfer behavior of a latent heat thermal energy storage unit with a finned tube[J].International Journal Heat Transfer，1993，36（8）：2083-2092.

[6]Castell A，Sole C，Medrano M，et al. Natural convection heat transfer coefficients in phase change material（PCM）modules with external vertical fins[J].Applied Thermal Engineering，2008，28（13）：1676-1686.

[7]Tao Y B，He Y L，Qu Z G. Numerical study on performance of molten salt phase change thermal energy storage system with enhanced tubes[J].Solar Energy，2012，86：1155-1163.

[8]Mat S，Al-Abidi A A，Sopian K，Sulaiman M Y，Mohammad A T.Enhance heat transfer for PCM melting in triplex tube with internal-external fins[J].Energy Conversion and Management，2013，74：223-236.

[9]Al-Abidi A A，Mat S，Sopian K，Sulaiman M Y，Mohammad A T.Numerical study of PCM solidification in a triplex tube heat exchanger with internal and external fins[J].International Journal of Heat and Mess Transfer，2013，61：684-695.

[10]Dow Corning Corporation. Properties of syltherm 800 heat transfer liquid[R]. Michigan，1985.

[11]Gregory J Z，Louay M C，Walters D K. Correlating heat transfer and friction in helically-finned tubes using artificial neural networks[J].International Journal of Heat and Mass Transfer，2007，50（23-24）：4713-4723.