相變材料的融化/凝固對建筑墻體壁面溫度的影響分析

陳占秀，馬秀琴，馬 赫

（河北工業(yè)大學 能源與環(huán)境工程學院，天津 300401）

相變材料的融化/凝固對建筑墻體壁面溫度的影響分析

陳占秀，馬秀琴，馬 赫

（河北工業(yè)大學 能源與環(huán)境工程學院，天津 300401）

采用enthalpy-porosity法建立建筑墻體內(nèi)置相變板層融化與凝固的傳熱模型，室內(nèi)外空氣溫度采用正弦周期溫度波邊界條件，計算了相變板層內(nèi)融化與凝固引起室內(nèi)壁面溫度的響應．結(jié)果表明，在相同的環(huán)境條件下，合適的相變溫度能使相變材料（PhaseChangeMaterial，PCM）完全融化/凝固，室內(nèi)壁面溫度波動變化最??；相變材料的相變溫度過高或過低，會導致部分PCM融化/凝固，過高或過低的相變溫度將導致室內(nèi)壁面溫度波動增大；相變過程溫度范圍越窄，相變材料融化與凝固的份額越高，引起壁面溫度波動也越?。惠^大的相變潛熱量可以使相變材料儲存較多的熱量，但相變潛熱量超過一定量時，再增大潛熱量對室內(nèi)壁溫影響較?。^大的導熱系數(shù)有利于強化傳熱，但也導致熱阻減小，不利于保溫，導熱系數(shù)的綜合效果對室內(nèi)壁溫影響較?。?/p>

相變材料；相變溫度；溫度波；融化/凝固；相變潛熱

0 引言

相變儲能材料在相變過程中吸收或釋放的熱量，將白天儲存的熱量在晚上釋放，在時間和空間上可以解決能量供求不均衡的矛盾．將相變材料與建筑圍護結(jié)構(gòu)結(jié)合可以增加建筑圍護結(jié)構(gòu)的等效熱容，提高可再生能源的利用率，減少室內(nèi)溫度波動同時增加室內(nèi)環(huán)境的舒適度，相變材料儲熱技術(shù)的研究與應用正在受到更多的關(guān)注[1-5]．

Alawadhi等[6]采用數(shù)值和實驗方法分析了在混凝土屋頂上內(nèi)填充PCM，發(fā)現(xiàn)通過室內(nèi)空間的熱通量可以減少39%，室內(nèi)環(huán)境溫度波動減?。甂alousck等[7]在夏熱冬暖的氣候條件下，在屋頂上加入相變材料板層，與普通房間的熱舒適性對比，使用相變材料的屋頂，夏季房間內(nèi)表面的溫度和空氣溫度分別降低3.5℃和2.1℃．Bogdan M等[8]設(shè)計由二層不同PCM層和傳統(tǒng)混凝土層組成的墻體，采用潛熱法模擬計算了全年氣象條件下室內(nèi)溫度波動情況，發(fā)現(xiàn)此墻可以有效減少冷卻和加熱負載的峰值．Peippo等[9]在夏暖冬冷氣候條件下，把相變材料加入房間的墻體內(nèi)，平均每年可以節(jié)省15%的能耗，并發(fā)現(xiàn)相變材料蓄熱一般發(fā)生在高出房間平均溫度1～3℃的相變過程溫度范圍，節(jié)能量的多少與不同地區(qū)氣候條件有關(guān)．

周國兵[10]根據(jù)北京地區(qū)的氣象條件，計算內(nèi)墻為定型相變材料板在外界周期溫度波作用下的熱特性分析，結(jié)果表明室內(nèi)舒適度增加；肖偉等[11]研究了輕質(zhì)建筑中相變蓄能石膏板作為內(nèi)隔墻體在冬季使用的效果，并給出了北京地區(qū)相變內(nèi)隔墻房間的設(shè)計原則．Zhang Yinping等[12]研究了相變墻板的儲熱和非線性傳熱熱性，對地板和墻體內(nèi)選擇和應用相變材料具有指導作用．

一般對相變墻體的相變潛熱采用等效熱容的方法計算，對相變材料的融化與凝固特性不考慮，而實際的相變材料都有一定的相變溫度和相變過程溫度范圍，相變材料是在其相變溫度范圍內(nèi)實現(xiàn)儲熱和放熱的．本文根據(jù)天津地區(qū)的氣象特點，從相變材料的融化與凝固特性分析，采用enthalpy-porosity[13]法建立建筑墻體中相變材料的融化與凝固模型，對相變材料的相變溫度、相變過程溫度范圍變化、相變潛熱和導熱系數(shù)等分析相變材料的儲與放熱對墻體壁溫的影響．

1 模型的建立

1.1 采用焓-孔隙率法建立融化與凝固模型

為了簡化計算，采用的墻體模型由2層材料組成，傳統(tǒng)墻體為紅磚墻和水泥砂漿構(gòu)成，相變墻體為紅磚墻和定形復合相變材料構(gòu)成，如圖1所示．

圖1 建筑墻體的物理模型Fig.1 Physicalmodelofwall

采用enthalpy-porosity法[13]建立建筑墻體伴有融化與凝固的非穩(wěn)態(tài)一維熱傳導模型，以焓和溫度為函數(shù)，能量控制方程如式 (1)和式2所示

式 (1)中的焓h和溫度T作為待求函數(shù)，k為導熱系數(shù)．由于相變材料有凝固狀態(tài)、固液兩相共存的糊狀區(qū)及液相狀態(tài)，計算焓值時，需要采用不同區(qū)域的比熱值，如式 (2)所示，其中 cp,s、cp,m和cp,l分別表示固體、固液兩相共存及液相的比熱．

采用enthalpy-porosity法求解導熱方程式 (1)，墻體內(nèi)的焓值由式 (2)確定，其中的融化與凝固狀態(tài)計算不直接追蹤融化、凝固交界面，而是引入液相分數(shù)來表征控制單元體內(nèi)PCM中液相所占據(jù)的份額．液相分數(shù)被定義為式 (3)[13]

式中：Ts和Tl分別為凝固點和融化點溫度；液相分數(shù) 表示融化的物質(zhì)占相變材料的體積百分數(shù)．假設(shè)PCM在固液兩相共存糊狀區(qū)內(nèi)的比熱容和導熱系數(shù)均呈線性變化，用式 (4)、式 (5)[13]表示

式中：cp,s、cp,l分別表示相變材料固、液相的比熱容；kp,s、kp,l分別表示固、液相導熱系數(shù)； 為液相分數(shù)．

為了簡化能量方程的求解過程，做以下3個假設(shè)：

1）固-液兩相共有區(qū)的物性參數(shù)（如導熱系數(shù)、比熱及密度）按常量計算；

2）忽略PCM液相區(qū)可能存在的自然對流；

3）忽略PCM凝固產(chǎn)生的過冷效應．

1.2 設(shè)置邊界和初始條件

假設(shè)室內(nèi)外空氣溫度為正弦函數(shù)規(guī)律變化的溫度波，為便于分析，輻射換熱當量成等效溫度的形式；導熱僅沿墻體厚度方向，其它表面絕熱，以天津地區(qū)冬季氣候條件為基礎(chǔ)，室外溫度變化為261～273 K，室內(nèi)溫度變化為289～297 K，室內(nèi)、室外對流換熱系數(shù)分別為8.7W/m2K和23W/m2K．作為對比計算的相變材料的物性參數(shù)為：相變溫度287.3～287.8 K、相變潛熱量60 J/g，相變板層厚度50mm、導熱系數(shù)0.8W/m2K．

墻體與室內(nèi)外空氣換熱的邊界條件為

1.3 數(shù)值計算

建筑墻體模型采用控制容積有限差分法對方程進行離散求解，整個計算區(qū)域采用六面體網(wǎng)格，相變板層局部區(qū)域進行網(wǎng)格加密，在每個時間計算節(jié)點處殘差均小于107時認為迭代收斂．

1.4 模型驗證

根據(jù)本文建立的模型，邊界條件與吳梁玉等[13]等計算交變溫度作用下建筑墻體的邊界條件相同，計算內(nèi)壁溫的分布，結(jié)果如圖2所示，計算結(jié)果與文獻 [13]結(jié)果非常接近，誤差可能因為計算精度不同導致的，因此證實所采用的模擬計算方法是正確的．

圖2 模型驗證Fig.2 Indoor temperature validation data

2 計算結(jié)果及分析

2.1 相變溫度的影響

在室內(nèi)外溫度波作用下，計算室內(nèi)壁面溫度和紅磚墻與相變材料貼緊一側(cè)的液相分數(shù)．

復合相變材料中的PCM具有相同的相變過程溫度范圍（0.5 K），3種PCM的融化與凝固溫度Tm范圍分別為285.3～285.8K、287.3～287.8K和289.3～289.8K，計算結(jié)果如圖3、圖4所示．由圖3可知，當相變溫度Tm范圍為287.3～287.8K時，液相分數(shù)曲線剛好經(jīng)歷從0到1完整的正弦函數(shù)曲線變化；PCM在1 d溫度變化范圍內(nèi)的高溫段內(nèi)可以完全融化儲存最多的熱量，在低溫段凝固可以釋放全部儲存的熱量．在圖4中表現(xiàn)為室內(nèi)壁面溫度波動最小．

當相變溫度Tm為285.3～285.8 K，低于287.3～287.8 K，PCM較長時間處于高溫段，液相分數(shù)曲線在較長時間等于1，最小的液相分數(shù)值為0.4，PCM能全部融化但不能全部凝固，相變材料在1 d的高溫段所蓄熱量在低溫段時間內(nèi)不能得到有效釋放；當相變溫度Tm為289.3～289.8 K，PCM較長時間內(nèi)處于較低溫度環(huán)境，液相分數(shù)在較長時間等于0，不能充分蓄熱，也不能有效地利用PCM的蓄放熱性能．這2種情況在圖4中均表現(xiàn)為室內(nèi)壁面溫度波動較大．說明當相變溫度Tm范圍過低（285.3～285.8 K）或過高（289.3～289.8 K）時，相變材料有一個時間段不放熱或是不蓄熱，只能部分轉(zhuǎn)移能量，從而使得室內(nèi)壁面溫度波動增加，只有相變溫度Tm范圍能實現(xiàn)完全融化和凝固時，蓄熱和放熱量最大，使得室內(nèi)壁面溫度幅值變化最小，如圖4所示．

圖3 相變溫度不同時液相分數(shù)分布Fig.3 Liquid fraction varied with phase change temperature

圖4 相變溫度不同時相變板內(nèi)層溫度分布Fig.4 The distribution of indoorwall temperature varied with phase change temperature

2.2 相變過程溫度范圍對液相分數(shù)的影響

相變過程溫度范圍選為：287.3～287.8 K、287.3～288.8 K、287.3～290.8 K，相變過程溫度范圍依次為0.5 K、1 K和3 K．計算結(jié)果如圖5、圖6所示．由圖5可知，在1 d的高溫段，相變過程溫度范圍為287.3～287.8K時，相變材料可以完全融化和凝固，表現(xiàn)其液相分數(shù)剛好經(jīng)歷從0到1正弦函數(shù)曲線變化；在相變過程溫度范圍為1 K和3 K時，相變材料的液相分數(shù)的最大值減小，不能到達1，說明相變材料不能全部融化，相變過程溫度范圍變大，蓄熱量變小，在圖6中，相變過程溫度范圍為0.5 K時，室內(nèi)壁溫波動最小，而當相變過程溫度范圍增大到1K和3K時，室內(nèi)壁溫波動隨之增加；也說明PCM轉(zhuǎn)移的熱量隨著相變過程溫度范圍的增加而減少，導致室內(nèi)壁溫波動增加．

圖5 相變過程溫度范圍不同時液相分數(shù)分布Fig.5 Liquid fraction varieswith phase change temperature range

圖6 相變過程溫度范圍不同時相變板層內(nèi)壁溫分布Fig.6 Thedistributionof indoorwallsurface temperature varieswith phase change temperature range

2.3 PCM相變潛熱對融化與凝固過程的影響

相變溫度287.3～287.8K范圍內(nèi)，計算相變潛熱量分別為30 J/g、60 J/g和90 J/g的液相分數(shù)和室內(nèi)壁面溫度的變化，結(jié)果如圖7、圖8．當相變潛熱值為30 J/g和60 J/g，PCM能夠完全融化和完全凝固；但是在潛熱值為30 J/g時，PCM較長時間處于等于1和0不變化的狀態(tài)，處于既不吸熱也不放熱的狀態(tài)，而不是剛好等于1和0的正弦函數(shù)曲線變化規(guī)律．從圖8室內(nèi)壁面溫度波動分析，潛熱量為30 J/g時的室內(nèi)壁面溫度波動大于潛熱量為60 J/g時的室內(nèi)壁面溫度波動，說明相變潛熱量為30 J/g時所儲存的熱量明顯小于相變潛熱量為60 J/g時所儲存的熱量．當相變潛熱量為90 J/g時，液相分數(shù)不能達到1也不能達到0，PCM既不能完全融化也不能完全凝固，但是與相變潛熱量為60 J/g的壁溫比較發(fā)現(xiàn)，兩者相變板層內(nèi)外溫度曲線基本重合，說明2種相變材料的儲熱量基本相同，即潛熱量達到一定值后，增加PCM潛熱量也不會吸收更多的能量．與周國兵[10]計算的相變潛熱量有極值情況吻合．

圖7 相變潛熱量不同時液相分數(shù)分布Fig.7 Liquid fraction varieswith latentheat

圖8 相變潛熱量不同時相變板內(nèi)層溫度分布Fig.8 The distribution of indoor surface temperature varieswith latentheat

2.4 相變材料的導熱系數(shù)對融化與凝固的影響

導熱系數(shù)分別選取0.4W/m2K 、0.8W/m2K 、1.2W/m2K，計算結(jié)果如圖9、圖10所示．隨著PCM 導熱系數(shù)的增加，液相分數(shù)變化幅值逐漸由小變大，說明相變材料由部分融化到完全融化，蓄放熱量增加，當導熱系數(shù)較大時能較多的轉(zhuǎn)移能量，但是室內(nèi)壁面溫度波動沒有隨著液相分數(shù)的變化幅度而產(chǎn)生較大的波動．從圖10中發(fā)現(xiàn)，隨著導熱系數(shù)的增加室內(nèi)壁溫波動并不大，表明室內(nèi)壁面溫度不僅與相變材料蓄放熱量的大小有關(guān)，而且還和導熱系數(shù)相關(guān)，較大的導熱系數(shù)雖然能使傳熱效果增加，使得PCM可以吸收較多的熱量，但同時也使墻體的換熱增加，墻體的熱阻減小了，表現(xiàn)在室內(nèi)壁溫仍然較低．

圖9 導熱系數(shù)不同時液相分數(shù)分布Fig.9 Liquid fraction varieswith thermalconductivity

圖10 導熱系數(shù)不同時相變板內(nèi)層溫度分布Fig.10 Indoor surface temperature varieswith thermal conductivity

3 結(jié)論

本文計算了不同相變溫度、相變過程溫度范圍、相變潛熱和導熱系數(shù)對相變材料融化與凝固及其對室內(nèi)壁面溫度的影響，得出以下結(jié)論：

1）PCM的物性參數(shù)不同時，其融化與凝固量是不同的，從而造成吸熱量和放熱量的不同．其中相變溫度是影響相變板層重要因素，合適的相變溫度能使相變材料完全融化和凝固，過高的相變溫度和過低的相變溫度，只能使PCM部分融化與凝固，不能起到充分蓄放熱的目的．

2）相變過程溫度范圍較大時，融化和凝固的溫度區(qū)間都將變窄，使得PCM不能完全融化或凝固，相對蓄放熱量也減小，因此選擇相變材料應盡量采用較小的相變過程溫度范圍的材料．

3）相變潛熱量較大可以吸收較多的能量，當潛熱量增加到一定值時，相變材料只能達到最大吸熱量，潛熱量再增加不能增加吸熱量．

4）相變材料板導熱系數(shù)較小時，相變材料融化與凝固量減小，但是相變材料的熱阻較大；導熱系數(shù)較大時有利于吸熱放熱，但PCM板層的熱阻較小，兩者綜合起來對相變板層室內(nèi)壁面溫度影響不大．

[1]SoaresN，Costa JJ，GasparA R，etal．Review ofpassive PCM latentheatthermalenergy storagesystems towardsbuildings'energy efficiency[J]．Energy and Buildings，2013，59：82-103．

[2]Kuznik F，Virgone J．Experimental investigation of wallboard containing phase changematerial：data for validation of numericalmodeling[J]．Energy and buildings，2009，41（5）：561-570．

[3]Mettawee EB S，Assassa GM R．Thermal conductivity enhancement in a latentheatstorage system[J]．Solar Energy，2007，81：839-845．

[4]張仁元．相變材料與相變儲能技術(shù) [M]．北京：科學出版社，2009．

[5]呂石磊，馮國會，朱能．脂酸類相變材料在節(jié)能建筑中應用的可行性研究 [J]．沈陽建筑大學學報（自然科學版），2006，22（1）：129-132．

[6]A lawadhiEM，AlqallafH J．Building roofwith conicalholescontaining PCM to reduce thecooling load：numericalstudy[J]．Energyand Buildings，2008，40（3）：193-203．

[7]Kalousck M，Hirs J．Simulation of the summer indoor thermal com fort by using wallboard with phase changematerial[J]．Eurosun，Bologna，2002，19：263-270．

[8]Bogdan M Diaconu，Mihao Cruceru．Novel concept of composite phase changematerialwall system for year-round thermal energy savings[J]．Energy and buildings，2010，42：1759-1772．

[9]Peippo K，Kauranen P，Lund PD．Amulti-componentPCM walloptim ized for passive solarheating[J]．Energy Build，1991，17：259-270．

[10]周國兵．溫度波作用下建筑內(nèi)墻定型相變材料板的熱特性 [J]．太陽能學報，2011，32（8）：1211-1216．

[11]肖偉，王馨，張寅平，等．輕質(zhì)建筑中相變蓄能石膏板熱性能研究 [J]．建設(shè)科技，2008，123（10）：84-88．

[12]Zhang Yinping，Lin Kunping，Jiang Yi，et al．Thermal storage and nonlinear heat-transfer characteristics of PCM wallboard[J]．Energy and buildings，2008，40：1771-1779．

[13]吳梁玉，陳永平，張程賓，等．交變溫度作用下填充PCM磚墻的動態(tài)傳熱特 [J]．太陽能學報，2012，33（9）：1593-1599．

[責任編輯 田 豐]

Analysisof indoor temperature influenced bymelting and solidification of phase changematerial inwallboard

CHEN Zhanxiu，MA Xiuqin，MA He

(Schoolof Energy and EnvironmentEngineering,HebeiUniversity of Technology,Tianjin 300401,China)

The thermalconductionmodelaccompaniedwith solidification andmelting processeswasestablished inwallboard incorporating phase changematerial layerusing enthalpy-porositymethod.The dynam ic oscillation of the outside and insidewall in solidification andmelting processesof PCM layerunder differentsinusoidal temperaturewave cycle of outdoor and indoorwas calculated.The results showed that phase change temperature plays an important role in phase transition under the same boundary conditions.The suitable phase transition temperaturemake PCM meltand solidify completely;however,if thehigheror lowerphase transition temperaturemake the PCM solidify ormeltpartly,itwould lead to greatervolatility in indoorwall temperature.When the phase transition temperature rangebetween themelting and solidification point tend to have smaller value,and the interval is too large,which lead to processofmelting and solidification happen in a narrow temperature range,the fluctuation of indoorwall temperature is still increased.The larger latentheatcaused smaller fluctuation of indoor temperature,thehigher thermalconductivity isconducive to theheattransfer,which also leads to heat resistance reduced,and isnotconducive to thermal insulation.

phase changematerial;phase transition temperature;temperaturewave;melting/solidification;latentheat of phase change

TK124；TB34

A

1007-2373(2016)01-0045-06

10.14081/j.cnki.hgdxb.2016.01.009

2015-04-10

河北省科技支撐計劃（14214304D，13274517，14243701D）

陳占秀（1969-），女（漢族），副教授，博士，chenzhanxiu@126.com．

數(shù)字出版日期：2015-12-09數(shù)字出版網(wǎng)址：http：//www.cnki.net/kcms/detail/13.1208.T.20151209.0854.002.htm l