李超 馬良棟* 張吉禮 趙天怡

大連理工大學(xué)建設(shè)工程學(xué)部

多相變材料蓄熱器蓄熱特性數(shù)值模擬研究

李超 馬良棟* 張吉禮 趙天怡

大連理工大學(xué)建設(shè)工程學(xué)部

本文基于單蓄熱介質(zhì)套管式蓄熱器蓄熱過程中換熱流體沿流向的溫度分布特點(diǎn)，提出了多相變材料套管式蓄熱器，并采用數(shù)值模擬的方法研究了該蓄熱器的蓄熱過程特性，分析了蓄熱過程中換熱流體沿流向的溫度變化特性，蓄熱材料的熔化及溫度分布特性。與單蓄熱介質(zhì)（切片石臘）相比，多相變材料（切片石臘和石蠟C16）蓄熱器換熱流體出口溫度降低了38%，且達(dá)到相同液相率時其蓄熱時間縮短了15%。結(jié)果表明，在蓄熱時，相變材料的熔化溫度隨換熱流體溫度降低而相應(yīng)降低將有利于提高蓄熱速率、增大單位時間內(nèi)的蓄熱能力。

多相變材料 蓄熱器 蓄熱特性 數(shù)值模擬

0 引言

相變蓄熱具有蓄能密度高、溫度波動小的特點(diǎn)在太陽能熱利用及廢熱回收利用領(lǐng)域得到廣泛利用[1]。對于大多數(shù)相變蓄熱材料來說，由于導(dǎo)熱系數(shù)低，相變蓄熱側(cè)熱阻大，不利于相變傳熱，因此很多學(xué)者在如何強(qiáng)化蓄熱傳熱上做了大量工作[2～5]。另外，還有一些學(xué)者通過利用不同熔點(diǎn)蓄熱材料的組合來提高相變蓄熱效率。Farid[6]最先提出在同一相變蓄熱容器中利用熔點(diǎn)不同的蓄熱材料進(jìn)行蓄熱，發(fā)明了蓄、放熱過程可同時進(jìn)行的蓄熱裝置，裝置均分為兩段，放入等體積不同熔點(diǎn)的蓄熱材料，實(shí)驗(yàn)研究發(fā)現(xiàn)其蓄熱效率得到了很大提高。后來，F(xiàn)arid[7～9]進(jìn)一步提出由三種不同熔點(diǎn)PCM組成的蓄熱容器。PCM封裝在外表面絕熱的長方體蓄熱床內(nèi)，其潛熱儲、放熱過程傳熱速率提高了15%。王劍鋒、陳光明等[10]也提出了一種組合式相變材料蓄熱設(shè)備，相變材料封裝在水平圓管內(nèi)，不同種類的相變材料按相變溫度的高低順序自下而上排列，傳熱流體外掠叉排圓柱管簇流動，模擬結(jié)果表明相變速率可提高15%～25%左右。

上述研究沒有考慮換熱流體流向溫度變化對蓄熱容器蓄熱過程的影響。當(dāng)換熱流體熱容較小時，尤其為空氣時，換熱流體入口與出口溫差較大，此時換熱流體的溫降必然會對相變蓄熱過程造成影響。本文基于單蓄熱介質(zhì)套管式蓄熱器蓄熱過程中換熱流體沿流向的溫度分布特點(diǎn)，提出了多相變材料套管式蓄熱器的概念，在盡可能保證高品位熱能的條件下，提高相變蓄熱過程傳熱速率。通過Fluent軟件研究了蓄熱過程流體溫度變化、蓄熱材料溫度分布等特點(diǎn)，為太陽能熱利用蓄熱系統(tǒng)設(shè)計(jì)提供基礎(chǔ)數(shù)據(jù)。

1 物理問題及數(shù)學(xué)模型

1.1 物理模型

多相變材料蓄熱器結(jié)構(gòu)如圖1（a）所示，蓄熱器為圓柱形，蓄熱材料被封裝在換熱流體管道與殼體之間，換熱流體在管道內(nèi)自下向上流動，通過換熱管與蓄熱材料進(jìn)行熱交換。沿?fù)Q熱體的流動方向，相變材料蓄熱腔分為前、后兩部分，前端部分的蓄熱材料為切片石蠟[12]，其蓄熱腔體積占總體積的四分之三，后端部分的蓄熱材料為石蠟C16[13]，其蓄熱腔體積占總體積的四分之一。兩種蓄熱材料的物性參數(shù)見表1[13]、表2[12]所示。為便于分析，對物理模型做出如下假設(shè)：

1）各相變材料各向同性；

2）換熱管熱阻忽略不計(jì)；

3）除密度外，各相變材料各物性參數(shù)為常數(shù)，且固相與液相參數(shù)彼此獨(dú)立；

4）各相變材料密度變化僅考慮動量方程中與體積有關(guān)的項(xiàng)，其余各項(xiàng)中的密度為常數(shù)，且滿足Boussinesq假設(shè)。

根據(jù)上述假設(shè)，選取蓄熱系統(tǒng)中的某一蓄熱單元作為研究對象，如圖1（b）所示。蓄熱單元長度為600mm，半徑為16mm，換熱管半徑為4mm。

表1 PCM-2石蠟C16物性參數(shù)

表2 PCM-1切片石蠟物性參數(shù)

圖1 多相變材料蓄熱器結(jié)構(gòu)示意圖

1.2 數(shù)學(xué)模型

本文采用Solidification/Melting模型[9]建立相變材料的融化模型，其計(jì)算區(qū)域?yàn)閳D1（b）中的ABCD區(qū)域，考慮相變材料溶化后自然對流的影響，基于上述假設(shè)，其二維軸對稱連續(xù)性方程、動量方程及能量方程如下：

1）連續(xù)性方程

2）動量方程

3）能量方程

式中：x為換熱流體主流流動方向；r為徑向；u為主流速度；v為徑向速度；H為焓值；ρ為密度；cp為比熱；μ為流體粘滯系數(shù)；β為熱膨脹系數(shù)；ε為固液共存區(qū)的連續(xù)數(shù)（在104～107之間，本文取105）；f為液相率。

式中：TS為相變蓄熱材料凝固溫度；TL相變蓄熱材料融化溫度。

相變材料的焓H可以用其顯熱焓h和潛熱ΔH來計(jì)算，即

顯熱焓h可表示為

式中：h0為參考焓；T0為參考溫度。

潛熱ΔH可表示為

式中：r為相變潛熱。

1.3 邊界條件與初始條件

傳熱流體與相變材料的換熱面（換熱管表面）為耦合壁面，其耦合邊界為：qw|fluid=qw|phase。

對于初始溫度，傳熱流體與相變材料區(qū)域應(yīng)該是相等的，即

式中：T0為蓄熱單元初始溫度，為了利用相變材料的顯熱儲存能力，其初始溫度一般低于相變溫度。

2 網(wǎng)格化分與參數(shù)設(shè)置

本文采用Gambit軟件建立物理模型并對其進(jìn)行網(wǎng)格劃分，網(wǎng)格采用四邊形單位，其網(wǎng)格數(shù)為9600。能量方程、動量方程的對流項(xiàng)離散采用二階迎風(fēng)差分格式，壓力梯度項(xiàng)采用PRESTO格式，壓力和速度的耦合采用SIMPLE算法。相變材料的密度、比熱容、導(dǎo)熱系數(shù)以及粘滯系數(shù)采用分段多項(xiàng)式表示（溫度的函數(shù)）。采用連續(xù)性方程和能量方程的殘差控制方程的迭代次數(shù)，計(jì)算中設(shè)置連續(xù)性方程為10-4，能量方程的殘差為10-8，時間步長選取0.2s。

3 結(jié)果分析

圖2給出了多相變材料蓄熱器與單相變材料蓄熱器液相率的比較，由圖2可知，與普通套管式相比，采用兩種蓄熱材料后融化時間縮短約46min，熔化時間縮短了近15%，并且，從圖2中可以看出，兩種蓄熱材料的液相率變化接近一致，可以說二者近似同步熔化，石蠟C16熔化稍稍快些。選取150min與300min兩時間點(diǎn)，分析換熱流體斷面平均溫度沿流向的變化，如圖3所示，可以看出，在150min時，對于單相變材料蓄熱容器，第一段溫降最大，以后依次逐漸減小，越接近出口處，溫降越來越小，沿?fù)Q熱流體流向，換熱流體溫度的降低是非線性的，即：越靠近出口端，溫度變化越小；而使用多種蓄熱材料后，在x=480mm處斷面平均溫度陡降，直接改變了曲線原來的走勢。并且，x=420mm之前的各斷面平均溫度并沒有發(fā)生很大變化，說明分段后前四分之三的換熱效果并沒有受到很大影響，并且，在相變材料改變處，溫度梯度大幅增加，在最后四分之一段，溫差增大了2倍多。當(dāng)蓄熱過程進(jìn)行到300min時，使用多種相變材料的蓄熱期內(nèi)蓄熱材料熔化基本完全，此時顯熱蓄熱成為主導(dǎo)，可以看到斷面平均溫度大大升高，說明了顯熱蓄熱效果遠(yuǎn)不及潛熱蓄熱效果。

圖2 相變材料液相率對比

圖3 多相變材料換熱流體軸向溫度分布

圖4給出了不同截面處斷面平均溫度隨時間的變化曲線，可以看出，熔化剛開始的一小段時間內(nèi)，溫度變化較快，這是因?yàn)橄嘧儾牧吓c內(nèi)導(dǎo)熱管直接接觸的部分只存在金屬導(dǎo)熱熱阻，熱量很快傳遞到相變材料中，使得這部分相變材料熔化速度很快。然后，傳熱方式變?yōu)閷?dǎo)熱與對流換熱，并且由于相變材料導(dǎo)熱率很低，致使金屬導(dǎo)熱熱阻可以忽略不計(jì)。對比圖4（d）、圖（e）、圖（f）中可以看出，在選用另一種蓄熱材料以后，后四分之一段斷面平均溫度明顯有所降低，并且兩條曲線存在交點(diǎn)。將兩種情況下?lián)Q熱流體的出口溫度進(jìn)行比較，如圖4（f）所示?？梢钥闯?，換熱流體入口與出口溫差明顯增大，換熱流體出口溫度降低5℃左右，溫差增大近38%。在流速相同的情況下，溫差的增大說明新型套管式蓄熱容器的換熱效率明顯增強(qiáng)，蓄熱過程得到了很好的強(qiáng)化。同時，圖4（f）中兩條曲線存在交點(diǎn)，大約在18500s左右，此時，多相變材料蓄熱器內(nèi)相變材料已近乎全部熔化，出口溫升迅速加大。對比圖4（d）、圖4（e）、圖4（f）發(fā)現(xiàn)兩曲線交點(diǎn)逐漸向右偏移，這種現(xiàn)象產(chǎn)生的主要原因是：單一相變材料出口四分之一段溫差很?。ㄈ鐖D3所示），而使用兩種相變材料后出口四分之一段溫差明顯增強(qiáng)（如圖3所示），而且熔化加快，斷面溫升提前，所以出現(xiàn)了交點(diǎn)向右偏移的現(xiàn)象。對比圖4（a）、圖4（b）、圖4（c）發(fā)現(xiàn)在熔化開始后大部分時間內(nèi)，斷面平均溫度基本相同，這是因?yàn)榍八姆种蜗嘧儾牧舷嗤?，但不同點(diǎn)是一段時間后出現(xiàn)分離，分離原因?yàn)槭褂枚嘞嘧儾牧虾笕刍瘯r間縮短，相變材料熔化到一定程度后率先開始顯熱蓄熱。

選取30min、60min、100min和200min四個時間段，繪出液相率分布圖，如圖5所示。圖5（a）反映出在蓄熱過程中，入口段蓄熱材料在0～60min熔化快于其他處蓄熱材料，這是由于入口段換熱流體溫度最高所致，與圖5（b）中前四分之三段熔化規(guī)律相同，伴隨著蓄熱材料中自然對流作用的加強(qiáng)，這種優(yōu)勢逐漸減弱，對于圖5（b）中的后四分之一段，并沒有出現(xiàn)這種規(guī)律，是因?yàn)楹笏姆种欢屋^短，溫差不大所致。同時圖5（a）、圖5（b）兩圖都反映出蓄熱過程中左下角與右下角蓄熱材料最難熔化，產(chǎn)生這種現(xiàn)象的原因?yàn)椋簱Q熱表面蓄熱材料受熱熔化后，溫度較高的液態(tài)蓄熱材料在浮升力的作用下向上流動，所以上部蓄熱材料最先完全熔化，例如圖5（a）中200min液相率圖所示；同時溫度相對較低的液態(tài)蓄熱材料向下流動，致使右下角與左下角蓄熱材料熔化最慢，如圖5（b）中200min液相率圖所示。本文中自然對流對相變材料熔化影響與文獻(xiàn)[14]中實(shí)驗(yàn)規(guī)律一致。

圖4 軸向斷面平均溫度對比圖

圖5 相變材料液相率分布

4 結(jié)論

本文提出了針對換熱流體溫度降選擇不同熔點(diǎn)的蓄熱材料這一概念，首先研究了普通光管套管蓄熱器換熱流體的溫度降，針對換熱流體在后半段溫度降低很慢這一現(xiàn)象，將后四分之一段的蓄熱材料換成熔點(diǎn)較低的其他蓄熱材料。計(jì)算了兩種情況下，相變蓄熱材料的融化時間、換熱流體沿軸向的溫度分布并進(jìn)行了對比分析。模擬結(jié)果表明：當(dāng)換熱流體入口與出口溫差較大時，合理采用不同熔點(diǎn)的蓄熱材料能夠大大加快蓄熱速率。沿?fù)Q熱流體溫度降，合理選擇熔點(diǎn)不同的蓄熱材料，對于相變蓄熱有極大的幫助。

[1]崔海亭,楊峰.蓄熱技術(shù)及其應(yīng)用[M].北京:化學(xué)工業(yè)出版社, 2004

[2]Ismail K A R,Alves C L F,Modesto M S.Numerical and experimental study on the solidification of PCM around a vertical axially finned isothermal cylinder[J].Applied Thermal Engineering, 2001,21(1):53-77

[3]Lamberg P.Approximate analytical model for two-phase solidification problem in a finned phase-change material storage[J].App -lied Energy,2004,77(2):131-152

[4]Hendra R,Mahlia T M I,Masjuki H H.Thermal and melting heat transfer characteristics in a latent heat storage system using mikro [J].Applied thermal engineering,2005,25(10):1503-1515

[5]Mettawee E B S,Assassa G M R.Thermal conductivity enhancement in a latent heat storage system[J].Solar Energy,2007,81(7): 839-845

[6]Farid MM.Solar energy storage with phase change[J].Solar Energy Res,1986,(4):11-29

[7]Farid M M,Kanzawa A.Thermal performance of a heat storage module using PCM’s with different melting temperatures:mathematical modeling[J].Journal of solar energy engineering,1989, 111(2):152-157

[8]Farid M M,Husian R M.An electrical storage heater using the phase-change method of heat storage[J].Energy Conversion and Management,1990,30(3):219-230

[9]Farid M M,Khudhair A M,Razack S A K,et al.A review on phase change energy storage:materials and applications[J].Energy conversion and management,2004,45(9):1597-1615

[10]王劍鋒,陳光明,陳國邦.組合相變材料儲熱系統(tǒng)的儲熱速率研究[J].太陽能學(xué)報,2000,21(3):258-264

[11]張仁元.相變材料與相變儲能技術(shù)[M].北京:科學(xué)出版社,2009

[12]晉瑞芳.太陽能相變蓄熱裝置蓄熱性能實(shí)驗(yàn)及模擬研究[D].上海:東華大學(xué),2009

[13]房叢叢.相變蓄熱裝置的數(shù)值模擬與優(yōu)化[D].濟(jì)南:山東建筑大學(xué),2012

[14]李偉,趙軍,李新國.方形槽內(nèi)水平圓管外相變蓄熱過程的數(shù)值模擬[J].熱能動力工程,2012,27(2):181-186

Num e ric a l Sim ula tion of The rm a l Ene rgy Stora ge Cha ra c te ris tic s of Multiple Pha s e Cha nge Ma te ria ls

LI Chao,MA Liang-dong*,ZHANG Ji-li,ZHAO Tian-yi
Faculty of Infrastructure Engineering,Dalian University of Technology

Based on the temperature distribution of heat transfer fluid along the direction of the fluid flow in the thermal energy storage container with a single Phase Change Material(PCM),a thermal energy storage container of multiple phase change material is proposed.The characteristics of heat transfer during the thermal energy storage process are investigated by using numerical method.The variation of some physical parameters with time are analyzed during thermal energy storage process,such as the temperature of heat transfer fluid,the temperature of multiple phase change material,and the percentage of liquid formation.The outlet temperature of heat transfer fluid of multiple phase change material decreases 38%compared with single phase change material,and the thermal energy storage time decreases 15% at the same percentage of liquid formation.The results show that the decrease of melting temperature of phase change material along the direction of heat transfer fluid flow is helpful to decrease the time of thermal energy storage and increase the thermal storage capacity per unit time.

multiple phase change materials,thermal energy storage container,characteristics of thermal energy storage,numerical simulation

1003-0344（2015）05-018-5

2014-4-7

馬良棟（1976～），男，博士，副教授；大連理工大學(xué)建設(shè)工程學(xué)部（116024）；E-mail:liangdma@dlut.edu.cn

“十二五”國家科技支撐計(jì)劃（2013BAJ10B02）