相變蓄熱材料導(dǎo)熱系數(shù)對太陽能通風(fēng)墻性能的影響

穆林， 雷勇剛， 宋翀芳， 程遠(yuǎn)達(dá)

(太原理工大學(xué) 環(huán)境科學(xué)與工程學(xué)院， 山西 太原 030024)

相變蓄熱材料導(dǎo)熱系數(shù)對太陽能通風(fēng)墻性能的影響

穆林， 雷勇剛， 宋翀芳， 程遠(yuǎn)達(dá)

(太原理工大學(xué) 環(huán)境科學(xué)與工程學(xué)院， 山西 太原 030024)

通過對相變蓄熱型太陽能煙囪模型通風(fēng)蓄放熱變化過程的計算，分析比較不同相變蓄熱墻導(dǎo)熱系數(shù)對太陽能煙囪性能的影響.計算結(jié)果表明：吸熱板最大表面溫度隨著相變蓄熱墻導(dǎo)熱系數(shù)的增大越接近相變蓄熱墻的相變溫度；蓄熱階段，入口平均風(fēng)速隨著相變蓄熱墻導(dǎo)熱系數(shù)的增大而減小；放熱階段，入口平均風(fēng)速隨著相變蓄熱墻導(dǎo)熱系數(shù)的增大反而越大；相變蓄熱墻導(dǎo)熱系數(shù)越大，蓄熱型太陽能煙囪系統(tǒng)16 h的累計通風(fēng)量越高，但在導(dǎo)熱系數(shù)增大到0.66 W·(m·K)-1后，再增大材料導(dǎo)熱系數(shù)，累計通風(fēng)量幾乎不再增加.

太陽能煙囪； 相變材料； 蓄熱墻； 通風(fēng)； 導(dǎo)熱系數(shù)

太陽能煙囪是利用太陽能加熱空氣產(chǎn)生煙囪效應(yīng)，從而促進通風(fēng)效果的裝置，它具有節(jié)能效果顯著、與建筑結(jié)合方便、初投資相對較少等優(yōu)點.然而，太陽能易受氣候環(huán)境影響，限制了太陽能煙囪的發(fā)展應(yīng)用.國內(nèi)外學(xué)者對太陽能煙囪相變蓄熱的特性進行了許多研究[1-9].Arkar等[1]提出一種相變蓄熱模塊(LHTES)進行建筑自然冷卻通風(fēng)，對不同氣候地點自然冷卻潛力進行分析.Stritih等[2]設(shè)計結(jié)合黑色石蠟相變材料的太陽能煙囪，相變材料儲存的熱量被用來加熱和通風(fēng).Amori等[3]對具有相變溫度為65.6 ℃石蠟的太陽能煙囪開展了實驗和數(shù)值研究.David等[4]建立具有相變材料的物理模型，對自然、強制和混合對流條件下進行計算.Liu等[5]對具有不同厚度的雙層Na2CO3·H2O蓄熱墻的太陽能煙囪系統(tǒng)進行恒熱流實驗.Giacia等[6-7]對雙層通道結(jié)合SP-22相變材料進行了實驗研究，結(jié)果表明,應(yīng)用PCM夾層對于建筑的熱效果提高明顯.李風(fēng)雷等[8]利用TRNSYS軟件,模擬計算以R141b為制冷劑的蓄熱型太陽能噴射制冷系統(tǒng)在夏季典型日的性能.Li等[9-10]對具有42 ℃相變溫度石蠟的太陽能煙囪系統(tǒng)進行實驗測試，研究穩(wěn)定條件下,熱流密度和石蠟層導(dǎo)熱系數(shù)對太陽能煙囪系統(tǒng)熱效率和通風(fēng)傳熱特性的影響.吳彥廷等[11]建立了二維相變蓄熱太陽能煙囪數(shù)理模型，分析空氣通道寬度對太陽能相變蓄熱集熱墻系統(tǒng)熱特性的影響.李永財[12]研究了結(jié)合相變蓄熱材料的太陽能通風(fēng)屋頂系統(tǒng)，結(jié)果顯示9月份60 ℃的日間通風(fēng)量較大，晚間的40 ℃的通風(fēng)量較大.盧軍等[13]對具有相變溫度的相變材料的太陽能通風(fēng)屋頂進行數(shù)值模擬，獲得不同相變溫度下,系統(tǒng)的夜間通風(fēng)量的變化規(guī)律.王元明等[14]以石蠟為材料,制造出針對于太陽能煙囪應(yīng)用復(fù)合相變材料，結(jié)果表明,碳納米管的添加可以有效提高相變材料的熱導(dǎo)率和相變過程的傳熱速率.雖然對太陽能煙囪結(jié)合相變蓄熱的研究較多，但是蓄熱墻導(dǎo)熱系數(shù)對太陽能煙囪系統(tǒng)的性能影響的成果比較匱乏.因此，本文建立了相變蓄熱型太陽能煙囪系統(tǒng)的三維非穩(wěn)態(tài)模型，并用文獻[3]的實驗結(jié)果驗證了其準(zhǔn)確性，研究不同蓄熱墻導(dǎo)熱系數(shù)對太陽能煙囪系統(tǒng)通風(fēng)放熱性能的影響.

1 物理模型

(a) 物理模型 (b) 計算區(qū)域圖1 相變蓄熱型太陽能煙囪系統(tǒng)(單位：mm)Fig.1 Phase change and thermal storage type of solar chimney system (unit:mm)

相變蓄熱型太陽能煙囪系統(tǒng)由半透明玻璃蓋板、吸熱板、相變蓄熱材料板、絕熱保溫板組成，其物理模型如圖1(a)所示.其工作原理：半透明玻璃蓋板和吸熱板之間形成空氣通道，空氣由入口進入通道，被吸熱板加熱后從出口排出，強化了煙囪系統(tǒng)自然通風(fēng).相對于傳統(tǒng)的太陽能煙囪系統(tǒng)，本模型采取了相變蓄熱墻，通過蓄熱材料的相變蓄熱可有效提高系統(tǒng)的蓄熱效果，延長系統(tǒng)放熱時間.該結(jié)構(gòu)高度為2 m，寬度為1 m，深度為0.25 m，高深比為10∶1[15]，其中相變材料層厚度為0.05 m，玻璃蓋板和吸熱板厚度均為0.004 m.本模型X為正南方向，計算區(qū)域如圖1(b)所示.

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

數(shù)值模擬假設(shè)如下:1) 工作流體為不可壓縮牛頓流體；2) 流體密度只與溫度有關(guān)，滿足Boussinesq假設(shè)；3) 沿著吸熱板邊緣的熱傳導(dǎo)忽略不計；4) 相變材料相變溫度為定值，固液相物性參數(shù)一致，無固液相分離；5) 忽略接觸熱阻.使用有限容積法求解控制方程，模擬采用RNGk-ε湍流模型、融化/凝固模型和離散坐標(biāo)(DO)模型，速度與壓力之間的耦合采用SIMPLE算法，對流項的離散格式使用二階迎風(fēng)格式，離散方程的求解應(yīng)用亞松弛TDMA算法.

1) 連續(xù)性方程為

2) 動量方程為

3) 能量方程為

4)k-方程為

5)ε-方程為

式(5)～(7)中：ui為xi方向上的平均速度矢量；xi為坐標(biāo)；p為平均壓力；t為非穩(wěn)態(tài)時間；vt，v分別為絮流和層流黏滯系數(shù)；gi為i方向的重力加速度；β為空氣膨脹系數(shù)；T，Tin為平均溫度和室內(nèi)溫度；Γ為廣義擴散系數(shù)；k為湍流脈動動能；ε為流體脈動動能的耗散率；C1，C2為經(jīng)驗系數(shù)；Gk為湍流動能的產(chǎn)生項；Pr為Prandel數(shù)；σk，σε，σt取經(jīng)驗值.

6) 相變方程為

7) 發(fā)射和散射方程為

式(8),(9)中：H為焓值；r材質(zhì)所在位置；s為直射方向矢量；s′為散射方向矢量；δ為散射系數(shù)；n為折射系數(shù)；α為吸收系數(shù)；Ф為相位函數(shù)；Ω為空間立體角.

3 邊界條件和物性參數(shù)

3.1邊界條件

1) 太陽輻射得熱.太陽輻射強度是影響太陽能煙囪相變蓄熱的主要因素.通過給定時間、日期和位置的太陽光束方向和輻射照度，對法向直接輻照應(yīng)用ASHRAE手冊的晴朗天氣條件的方程為Edn=A/exp(B/sinβ)，對漫射輻射的表達(dá)式為Ed=CYEdn，計算得到的參數(shù)作為太陽射線跟蹤算法半透明邊界條件的輸入條件.玻璃蓋板邊界條件：玻璃蓋板外邊界采取第三類邊界條件，對流換熱系數(shù)h=4 W·(m2·K)-1，外界空氣溫度取夏季通風(fēng)溫度為300 K.

2) 其他壁面面邊界條件.吸熱面采取無內(nèi)熱源耦合邊界，絕熱木板為絕熱邊界.

3) 出入口邊界條件.入口采取壓力入口邊界條件，入口空氣溫度取太原夏季通風(fēng)溫度300 K,而出口采取壓力出口邊界條件.

4) 初始條件.初始溫度為300 K，初始相對壓力為0，初始速度為0.

3.2物性參數(shù)

玻璃蓋板(4 mm玻璃)、吸熱板(4 mm鋼板)、空氣和相變石蠟的物性參數(shù)，如表1所示.表1中：ρ為密度；λ為導(dǎo)熱系數(shù)；cp為比熱容；α為吸收率；θ為發(fā)射率；τ為透射率；t為相變溫度；Q為相變潛熱.

表1 材料的物性參數(shù)Tab.1 Physical properties of material

4 網(wǎng)格獨立性考核與模型有效性驗證

由于模型結(jié)構(gòu)形式比較規(guī)則，計算區(qū)域采用結(jié)構(gòu)化正六面體網(wǎng)格.圖2為計算網(wǎng)格(Z軸方向截面圖)，上面為相變材料，下面為空氣通道.為了使計算結(jié)果準(zhǔn)確且節(jié)省計算資源，進行了網(wǎng)格獨立性考核.對5套不同的網(wǎng)格下10：00的入口風(fēng)速進行計算，網(wǎng)格7 608，13 639，56 724，107 642，210 378時的入口平均風(fēng)速分別為2.210，1.320，0.253，0.224，0.223 m·s-1.可以看出，第4套和第5套網(wǎng)格的計算偏差小于1%，因此，得到了網(wǎng)格無關(guān)的獨立解.考慮到計算精度和計算速度要求，以第4套網(wǎng)格作為所有數(shù)值計算網(wǎng)格，其網(wǎng)格數(shù)為107 642.

圖2 太陽能煙囪相變蓄熱系統(tǒng)網(wǎng)格劃分Fig.2 Grid generation of phase change and thermal storage type solar chimney system

圖3 模擬結(jié)果與實驗數(shù)據(jù)比較Fig.3 Simulation results compared with experimental data

為了驗證該模型，根據(jù)文獻[3]實驗參數(shù)，建立970 mm×150 mm×2 000 mm的相變蓄熱型太陽能煙囪模型.其相變材料層厚度為40 mm，室外溫度和室外太陽輻射強度與文獻[3]一致.吸熱板平均溫度模擬結(jié)果與實驗結(jié)果比較,如圖3所示.

由圖3可知：模擬計算與實驗所測溫度隨時間變化趨勢基本吻合，最大溫度都是在13:00；放熱階段的偏差較大，主要是因為模擬計算時忽略了接觸熱阻和能量損耗.由圖3還可以看出：本模擬計算結(jié)果與文獻[3]的實驗結(jié)果誤差在±5%以內(nèi)，表明本模型的可靠性.

5 實驗結(jié)果與分析

5.1吸熱板表面平均溫度與時間的變化關(guān)系

圖4 吸熱板表面溫度隨時間變化Fig.4 Endothermic board surface temperature over different times

不同蓄熱墻導(dǎo)熱系數(shù)下，吸熱板表面平均溫度隨時間的變化關(guān)系，如圖4所示.由圖4可知：由于太陽的輻射，吸熱面的溫度均從初始時刻的300 K開始增大，到13：00時達(dá)到最大值，然后開始下降.這與太陽入射輻射13：00時出現(xiàn)最大值相對應(yīng).由圖4還可知：吸熱板表面最大溫度隨著相變蓄熱墻導(dǎo)熱系數(shù)的減小反而越大，采用1.10 W·(m·K)-1相變蓄熱墻時，吸熱板表面的最大溫度為317.62 K，而采用0.33 W·(m·K)-1相變蓄熱墻時，吸熱板表面最大溫度卻為322.2 K.這是由于相變材料達(dá)到相變溫度后開始相變蓄熱，導(dǎo)熱系數(shù)大時,吸熱板與相變材料之間的傳熱熱阻比較小，較多的熱量會被相變材料儲存，對同樣的入射輻射熱流，吸熱板表面溫度就比較低.

另外，相變蓄熱墻導(dǎo)熱系數(shù)越低,吸熱板蓄熱和放熱時的溫差就越大，隨著相變材料導(dǎo)熱系數(shù)的變大，這種差別逐漸降低直至消失，最終都會等于相變溫度.此時,蓄熱墻的儲存熱量達(dá)到最大值，導(dǎo)熱系數(shù)為0.33 W·(m·K)-1的相變材料吸熱板最大溫度與相變溫度的差值為6.2 K，而導(dǎo)熱系數(shù)為1.10 W·(m·K)-1時僅為1.62 K.由此可以看出，圖4中20：00時刻都有一個明顯的轉(zhuǎn)折點.這是因為太陽入射輻射在20：00降為0，系統(tǒng)進入單純放熱階段.

5.2入口平均風(fēng)速與時間的變化關(guān)系

不同蓄熱墻導(dǎo)熱系數(shù)下,入口平均風(fēng)速(v)隨時刻的變化關(guān)系，如圖5所示.由圖5可知：不同蓄熱墻導(dǎo)熱系數(shù)下，系統(tǒng)入口平均風(fēng)速在到達(dá)13：00以前均是不斷增大的，到13：00時均達(dá)到最大值，且相變蓄熱墻導(dǎo)熱系數(shù)越小，入口平均風(fēng)速越大，這與吸熱面溫度相對應(yīng)的；在18：00以后，相變蓄熱墻進入放熱階段.由圖5還可知：導(dǎo)熱系數(shù)越大，入口平均風(fēng)速越大，反之則越小.這是由于導(dǎo)熱系數(shù)小，相變蓄熱墻傳遞到空氣中的熱流就會少；而導(dǎo)熱系數(shù)大，相變蓄熱墻傳遞到空氣中的熱流就會越大.但是，當(dāng)導(dǎo)熱系數(shù)增加到0.66 W·(m·K)-1時，再增加，效果就會越來越不明顯；當(dāng)導(dǎo)熱系數(shù)為1.10 W·(m·K)-1時，幾乎和0.66 W·(m·K)-1相同.

不同通風(fēng)時間內(nèi)的累計通風(fēng)量(Q)，如圖6所示.由圖6可知：在通風(fēng)12 h以前，導(dǎo)熱系數(shù)越大,累計通風(fēng)量越大;通風(fēng)12 h以后，導(dǎo)熱系數(shù)越大,累計通風(fēng)量反而越小;而且在通風(fēng)12 h時，不同導(dǎo)熱系數(shù)下累計通風(fēng)量幾乎相等.但是,當(dāng)導(dǎo)熱系數(shù)超過0.66 W·(m·K)-1以后,再增大導(dǎo)熱系數(shù)對累計通風(fēng)量幾乎不再有影響.

垂直Y軸中心面局部速度流圖，如圖7所示.由圖7可知：在出口處有回流現(xiàn)象.

圖5 入口平均風(fēng)速隨時間變化 圖6 不同通風(fēng)時間內(nèi)累計通風(fēng)量 Fig.5 Entrance average wind speed change over time Fig.6 Total ventilation rate over different times

5.3吸熱板中心溫度沿高度方向的變化關(guān)系

吸熱板中心溫度沿高度方向的變化關(guān)系，如圖8所示.由圖8可知：吸熱板中心溫度在13:00沿高度方向均不斷增大，且高度在0.2～0.8 m和1.2～1.8 m時變化趨勢比較大，高度在0.8～1.2 m時變化趨勢變得比較平緩.另外，導(dǎo)熱系數(shù)越小，吸熱板溫度越高，導(dǎo)熱系數(shù)越大，吸熱板溫度越低，且隨著高度的增加，這種差別越來越大.由圖8還可知：在高度為0.2 m時，導(dǎo)熱系數(shù)0.33,1.10 W·(m·K)-1對應(yīng)的溫差幾乎為0 K;在高度為1.0 m時，導(dǎo)熱系數(shù)0.33和1.10 W·(m·K)-1對應(yīng)的溫差為6.67 K;而在高度為1.8 m時，導(dǎo)熱系數(shù)0.33和1.10 W·(m·K)-1對應(yīng)的溫差則增大為9.8 K.

此外，隨著導(dǎo)熱系數(shù)的增大，吸熱板中心溫度沿高度方向的變化越來越平緩，而且當(dāng)導(dǎo)熱系數(shù)增大到0.66 W·(m·K)-1后,再增大導(dǎo)熱系數(shù)，吸熱板中心溫度沿高度方向的變化越來越接近相變溫度.這和圖4結(jié)論相同.

圖7 中心面局部速度流場圖 圖8 吸熱板中心溫度沿高度方向的變化關(guān)系 Fig.7 Local flow field at core surface Fig.8 Heat absorbing plate core temperature along height direction

6 結(jié)論

建立太陽能煙囪相變蓄熱通風(fēng)系統(tǒng)的三維數(shù)理模型，并對太原市7月15日8：00-24:00的蓄放熱和通風(fēng)性能隨時間變化的關(guān)系開展數(shù)值模擬研究，得到以下3個主要結(jié)論.

1) 相變蓄熱墻導(dǎo)熱系數(shù)越低，吸熱板蓄熱和放熱時,溫差就越大；隨著相變蓄熱墻導(dǎo)熱系數(shù)的變大，這種差別逐漸降低直至消失，最終等于相變溫度，此時,蓄熱墻儲存熱量達(dá)到最大值；當(dāng)相變蓄熱墻導(dǎo)熱系數(shù)增大到0.66 W·(m·K)-1以后，再增大導(dǎo)熱系數(shù)，吸熱板表面溫度變化變得越來越弱.

2) 蓄熱階段入口平均風(fēng)速隨相變蓄熱墻導(dǎo)熱系數(shù)增大而增大，且均在13：00出現(xiàn)最大風(fēng)速；放熱階段入口平均風(fēng)速隨相變蓄熱墻導(dǎo)熱系數(shù)增大反而減小.導(dǎo)熱系數(shù)低的相變蓄熱墻熱阻比較大，白天蓄熱階段存儲的熱量較少，更多的熱量用于通風(fēng)使入口風(fēng)速較大；晚上放熱階段放出的熱量較少，使入口風(fēng)速變得較小.在通風(fēng)12 h以前，導(dǎo)熱系數(shù)越大累計通風(fēng)量越大;通風(fēng)12 h以后，導(dǎo)熱系數(shù)越大累計通風(fēng)量反而越小，而且在通風(fēng)12 h時，不同導(dǎo)熱系數(shù)下累計通風(fēng)量幾乎相等.但是,當(dāng)導(dǎo)熱系數(shù)超過0.66 W·(m·K)-1以后,再增大導(dǎo)熱系數(shù)后,對累計通風(fēng)量幾乎不再有影響.

3) 吸熱板中心溫度沿高度方向不斷增大，而且吸熱板中心溫度沿高度方向的變化隨著導(dǎo)熱系數(shù)越大越來越平緩.當(dāng)導(dǎo)熱系數(shù)增大到0.66 W·(m·K)-1后,再增大導(dǎo)熱系數(shù)，吸熱板中心溫度沿高度方向的變化越來越接近相變溫度.

[1] ARKAR C，MEDVED S.Free cooling of a building using PCM heat storage integrated into the ventilation system[J].Solar Energy，2007，81(9):1078-1087.DOI:10.1016/j.solener.2007.01.010.

[2] STRITIH U，NOVAK P.Solar heat storage wall for building ventilation[J].Renewable Energy，1996，8(1/2/3/4):268-271.DOI:10.1016/0960-1481(96)88860-4.

[3] AMORI K E，MOHAMMED S W.Experimental and numerical studies of solar chimney for natural ventilation in Iraq[J].Energy and Buildings，2012，47(4):450-457.DOI:10.1016/j.enbuild.2011.12.014.

[4] DAVID D，KUZNIK F,ROUX J J.Numerical study of the influence of the convective heat transfer on the dynamical behaviour of a phase change material wall[J].Applied Thermal Engineering，2011，31(16):3117-3124.DOI:10.1016/j.applthermaleng.2011.04.

[5] LIU Bin,LI Peng,MA Xiaoyan,etal.Chimney effect of solar hybrid-double wall with different thickness PCM of Na2CO3·10H2O[J].Energy Procedia,2015,75:521-526.DOI:10.1016/j.egypro.2015.07.443.

[6] GRACIA A D,NAVARRO L,CASTELL A,etal.Experimental study of a ventilated facade with PCM during winter period[J].Energy and Buildings,2012,58(2):324-332.

[7] GRACIA A D,NAVARRO L,etal.Solar absorption in a ventilated facade with PCM: Experimental results[J].Energy Procedia,2012,30(1):986-994.DOI:10.1016/j.enbuild.2012.10.026.

[8] 李風(fēng)雷,任艷玲.典型氣象日蓄熱型太陽能噴射制冷系統(tǒng)性能分析[J].華僑大學(xué)學(xué)報(自然科學(xué)版),2014,35(2):185-190.DOI:10.11830/ISSN.1000-5013.2014.02.0185.

[9] LI Yongcai,LIU Shuli.Experimental study on thermal performance of a solar chimney combined with PCM[J].Applied Energy,2014,114:172-178.DOI:10.1016/j.apenergy.2013.09.022.

[10] LI Yongcai,LIU Shuli,SHUKLA A.Experimental analysis on use of thermal conductivity enhancers (TCEs) for solar chimney applications with energy storage layer[J].Energy and Buildings,2016,116:35-44.DOI:10.1016/j.enbuild.2015.12.049.

[11] 吳彥廷，周國兵.太陽能相變蓄熱集熱墻二維非穩(wěn)態(tài)模型及分析[J].太陽能學(xué)報,2010，33(6):949-952.

[12] 李永財.太陽能通風(fēng)技術(shù)結(jié)合相變蓄熱材料自然通風(fēng)效果實驗研究[D].重慶:重慶大學(xué),2009:1-101.

[13] 盧軍，趙娟,黃光勤,等.太陽能蓄能通風(fēng)系統(tǒng)理論模型[J].土木建筑與環(huán)境工程,2012，34(3):64-74.

[14] 王元明.石蠟/碳納米管復(fù)合相變材料的性能研究[J].現(xiàn)代化工，2015，35(1):130-134.

[15] RYAN D，BUREK S A M.Experimental study of the influence of collector height on the steady state performance of a passive solar air heater[J].Solar Energy，2010，84(9):1676-1684.DOI:10.1016/J.SOLENER.2010.05.018.

(責(zé)任編輯： 黃曉楠英文審校： 劉源崗)

EffectofThermalConductivityonPerformanceofSolarChimneyWithPhaseChangeAccumulator

MU Lin， LEI Yonggang， SONG Chongfang， CHENG Yuanda

(School of Environment Science and Engineering， Taiyuan University of Technology， Taiyuan 030024， China)

The influence of different thermal conductivities of phase change and thermal storage wall on the performance of solar chimney was analyzed and compared through the calculation of the ventilation and heat storage or release change process of the phase change and thermal storage type of solar chimney model. The results show that the maximum temperature of the heat absorber plate is more close to the phase change temperature with the increase of the thermal conductivity of the phase change and heat storage wall. The average inlet wind velocity reduces with the increase of the thermal conductivity for the process of heat storage， while increases for the process of heat release. The higher the thermal conductivity of the phase change and thermal storage wall is, the higher the cumulative ventilation in 16 hours of the regenerative solar chimney system will be.But when the thermal conductivity is increased to 0.66 W·(m·K)-1, the cumulative ventilation has almost no increase even if increasing the thermal conductivity of the material.

solar chimney； phase change materials； thermal storage wall； ventilation; thermal conductivity

10.11830/ISSN.1000-5013.201612040

TU 834

A

1000-5013(2017)06-0830-06

2016-12-20

雷勇剛(1976-)，男，教授，博士，主要從事低功耗強化傳熱技術(shù)、供熱技術(shù)和建筑節(jié)能技術(shù)的研究.E-mail:yongganglei@gmail.com.

國家自然科學(xué)基金資助項目(51106103, 51408391)