墻體熱、濕及空氣耦合傳遞非穩(wěn)態(tài)模型及驗證

劉向偉　陳國杰　陳友明

摘要：為了更準(zhǔn)確地預(yù)測墻體內(nèi)的溫濕度分布，研究多孔介質(zhì)墻體內(nèi)的熱、濕及空氣耦合非穩(wěn)態(tài)傳遞規(guī)律，以溫度、相對濕度和空氣壓力為驅(qū)動勢，考慮熱傳遞、濕傳遞、空氣滲透及其相互作用，建立了建筑多孔介質(zhì)墻體熱、濕及空氣耦合傳遞非穩(wěn)態(tài)模型，并采用有限元方法設(shè)計了相應(yīng)的模擬計算程序.通過對比新建模型模擬結(jié)果與國際公認(rèn)的HAMSTAD標(biāo)準(zhǔn)驗證實例，驗證了模型的正確性.

關(guān)鍵詞：建筑墻體；非穩(wěn)態(tài)；熱、空氣、濕耦合傳遞；含濕量；相對濕度

中圖分類號：TU111.4 文獻標(biāo)識碼：A

Abstract：A coupled heat， air and moisture transfer model， which takes into consideration the heat transfer， moisture transfer and air convection and their coupled effect， was developed to predict the distribution of the temperature and humidity and to investigate the rule of the coupled heat， air and moisture transfer in walls. The temperature， relative humidity and air pressure were chosen as the driving potentials. A program based on the finite element method was developed to calculate the governing equations. And the numerical results of this model were compared with the internationally accepted HAMSTAD benchmarks， and the results agree well with each other.

Key words：building walls； transient； heat， air and moisture transfer； moisture content； relative humidity

建筑墻體多為多孔介質(zhì)材料，建筑墻體熱濕傳遞研究以多孔介質(zhì)傳熱傳質(zhì)學(xué)為理論基礎(chǔ).Philip 與Devries[1] （1957年） 和Luikov[2] （1966年） 以溫度和含濕量為驅(qū)動勢，考慮多孔介質(zhì)內(nèi)熱傳遞、濕遷移及其相互作用建立了多孔介質(zhì)熱濕耦合傳遞模型.在Philip與Devries和Luikov模型的基礎(chǔ)上Kunzel[3]， Kong[4]， Chu[5]， Abahri[6]， Leskovsek[7]， Zhong[8]， Belarbi[9]， Qin[10]，李魁山[11]和郭興國[12]等各自建立了多孔介質(zhì)材料的熱濕耦合傳遞模型.空氣滲透對熱濕傳遞過程的影響在上述研究中均未考慮.

建筑墻體長期暴露在非穩(wěn)態(tài)氣候條件下，由于室內(nèi)外環(huán)境存在溫度、濕度及空氣壓力梯度，這將導(dǎo)致墻體內(nèi)的熱傳遞、濕遷移及空氣滲透.建筑墻體熱濕耦合傳遞研究中忽略空氣滲透對建筑墻體熱濕傳遞過程的影響將不利于準(zhǔn)確地分析建筑墻體內(nèi)的溫度和濕度分布.為了更準(zhǔn)確地研究建筑墻體內(nèi)的熱濕耦合傳遞規(guī)律，劉曉燕等[13]建立了建筑墻體熱、濕及空氣耦合傳遞模型.該模型從微觀的角度，通過定義氣體所占體積百分比與液態(tài)水所占體積百分比來分別計算水蒸汽含量和液態(tài)水含量，但是在自然條件下很難將水的物相分開測量，模型參數(shù)難以確定.另外，該模型采用室外空氣日平均溫度作為室外邊界條件，不能充分反應(yīng)室外氣候的逐時非穩(wěn)態(tài)變化.本文在KUNZEL的研究基礎(chǔ)上，通過考慮空氣滲透，以室外氣象條件作為邊界條件，建立了一個以溫度、相對濕度和空氣壓力為驅(qū)動勢的建筑墻體熱、濕及空氣耦合傳遞非穩(wěn)態(tài)模型.通過對比該模型模擬結(jié)果與HAMSTAD驗證實例[14]，驗證了模型的正確性.

1模型的建立

本文以相對濕度、溫度和空氣壓力為驅(qū)動勢，根據(jù)能量守恒和質(zhì)量守恒定律建立墻體熱、濕及空氣耦合傳遞非穩(wěn)態(tài)模型，并作如下假設(shè)：1）材料骨架是一個固定的、不變形的惰性骨架，不與液相、氣相發(fā)生化學(xué)反應(yīng)；2）墻體材料為各向同性的連續(xù)多孔介質(zhì)；3） 溫度對材料平衡含濕量的影響忽略不計；4）墻體材料中的水只有汽、液兩相，材料中始終存在局部濕平衡；5）孔隙內(nèi)的混合氣體（濕空氣）按理想氣體處理； 6）多層墻體層與層之間的濕傳遞主要受邊界濕傳遞阻的影響，如果不同層之間接觸十分緊密，則濕阻較小，可認(rèn)為這兩種材料的邊界表面處于濕平衡狀態(tài).

1.1濕傳遞

2模型求解

建筑墻體內(nèi)熱、空氣、濕傳遞過程相互耦合，為了獲得墻體內(nèi)的溫度和濕度分布，控制方程組需同時求解.本文采用有限元方法對控制方程組進行求解.時間步長可根據(jù)實際邊界條件確定，本文將時間步長設(shè)定為1 h.

3模型驗證

HAMSTAD驗證實例是為了評價建筑物理領(lǐng)域熱、空氣、濕傳遞機理模型而建立的一個開放性平臺.對比新建模型模擬結(jié)果與HAMSTAD驗證實例是目前國際公認(rèn)的建筑墻體熱、空氣、濕耦合傳遞模型驗證方法.

HAMSTAD包含了5個驗證實例，每個驗證實例至少包含熱傳遞、濕傳遞及空氣滲透機理中的兩項.本文通過對比新模型模擬結(jié)果與HAMSTAD驗證實例3和5來驗證模型的正確性.

3.1HAMSTAD驗證實例3

如圖 1～圖 6所示，新模型模擬結(jié)果與HAMSTAD驗證實例吻合良好，該模型能準(zhǔn)確地預(yù)測建筑墻體內(nèi)的溫度和濕度分布.圖中CTH， TUD， Technion和NRC表示參與HAMSTAD項目的其它研究機構(gòu)的模擬結(jié)果，新建模型模擬結(jié)果表示本文所提出模型的模擬結(jié)果.

如圖1，3，5所示，在20～21 d之間，由于壓力梯度方向的改變，室外低溫空氣向室內(nèi)滲透，墻體內(nèi)溫度迅速下降，快速接近室外空氣溫度.同樣，如圖2，4，6所示，由于室外空氣的含濕量低于墻體內(nèi)的含濕量，室外空氣向室內(nèi)滲透的過程中帶走大量的濕，墻體內(nèi)的含濕量快速降低.由此可見，建筑材料孔隙內(nèi)的空氣對流，對建筑墻體的熱濕性能有重要的影響.

4結(jié)論

本文在Kunzel的研究基礎(chǔ)上，通過考慮空氣滲透，建立了一個以溫度、相對濕度和空氣壓力為驅(qū)動勢的建筑墻體熱、濕及空氣耦合傳遞非穩(wěn)態(tài)模型，并采用有限元方法對該模型進行求解.該模型模擬結(jié)果與HAMSTAD驗證實例吻合良好，結(jié)果表明該模型能準(zhǔn)確地預(yù)測熱傳遞、濕傳遞及空氣滲透機理作用下建筑墻體內(nèi)的溫度和濕度分布.

參考文獻

[1]PHILIP J R，DEVRIES D A. Moisture movement in porous materials under temperature gradients[J]. Trans Am Geophys Union， 1957， 38（2）：222-232.

[2]LUIKOV A V. Heat and mass transfer in capillaryporous bodies[M]. Oxford： Pergamon， 1966：75-99.

[3]KUNZEL H M. Simultaneous heat and moisture transport in building components[R]. Verlag Suttgart： Fraunhofer IRB， 1995：38-40.

[4]KONG F， ZHENG M. Effects of combined heat and mass transfer on heating load in building drying period[J]. Energy and Buildings， 2008， 40（8）：1614-1622.

[5]CHU S S， FANG T H， CHANG W J. Modelling of coupled heat and moisture transfer in porous construction materials[J]. Mathematical and Computer Modelling， 2009， 50（7/8）：1195-1204.

[6]ABAHRI K， BELARBI R， TRABELSI A. Contribution to analytical and numerical study of combined heat and moisture transfers in porous building materials[J]. Building and Environment， 2011， 46（7）：1354-1360.

[7]LESKOVSEK U， MEDVED S. Heat and moisture transfer in fibrous thermal insulation with tight boundaries and a dynamical boundary temperature[J]. International Journal of Heat and Mass Transfer， 2011， 54（19/20）：4333-4340.

[8]ZHONG Z P. Combined heat and moisture transport modeling for residential buildings[D]. Indiana： Purdue University，2008：18-24.

[9]BELARBI R， QIN M， ABDELKARIM A M， et al. Experimental and theoretical investigation of nonisothermal transfer in hygroscopic building materials[J]. Building and Environment， 2008， 43（12）：2154-2162.

[10]QIN M， BELARBI R， ABDELKARIM A M， et al. Coupled heat and moisture transfer in multilayer building materials[J]. Construction and Building Materials， 2009， 23（2）：967-975.

[11]李魁山， 張旭， 高軍. 周期邊界條件下多層墻體內(nèi)熱濕耦合遷移[J]. 同濟大學(xué)學(xué)報：自然科學(xué)版， 2009， 37（6）： 814-818.

LI Kuishan， ZHANG Xu， GAO Jun. Coupled heat and moisture transfer in multiwall under periodic boundary conditions[J]. Journal of Tongji University： Natural Sciences， 2009， 37（6）：814-818. （In Chinese）

[12]郭興國， 陳友明. 一種新型木結(jié)構(gòu)墻體的熱濕性能分析[J]. 湖南大學(xué)學(xué)報：自然科學(xué)版， 2009， 36（12）：18-21.

GUO Xingguo， CHEN Youming. Analysis of the hygrothermal performance of a new timber structure wall[J].Journal of Hunan University： Natural Sciences， 2009， 36（12）：18-21.（In Chinese）

[13]劉曉燕， 賈永英， 王志國. 建筑墻體熱、濕及空氣耦合傳遞[J]. 太陽能學(xué)報， 2004， 25（1）： 13-18.

LIU Xiaoyan， JIA Yongying， WANG Zhiguo. The study of coupled transfer of heat， moisture and air in building walls[J]. Acta Energiae Solatis Sinica， 2004， 25（1）：13-18.（In Chinese）

[14]HAGENTOFT C E. Hamstadfinal report： methodology of HAMmodeling， Report R02：8[R]. Gothenburg： Department of Building Physics， Chalmers university of Technology， 2002：5-51.

[15]HUGO H. Building physicsheat， air and moisture： fundamentals and engineering methods with examples and exercises[M]. Berlin： Ernst & Sohn A Wiley， 2012：136.