劉建昌，李隆鍵，申憲文，馮 雅

(1.重慶大學(xué) 動(dòng)力工程學(xué)院, 重慶 400030； 2.中國建筑西南設(shè)計(jì)研究院， 成都 610042)

EPS板薄抹灰外墻外保溫系統(tǒng)熱濕耦合傳遞的數(shù)值模擬

劉建昌1，李隆鍵1，申憲文1，馮 雅2

(1.重慶大學(xué) 動(dòng)力工程學(xué)院, 重慶 400030； 2.中國建筑西南設(shè)計(jì)研究院， 成都 610042)

EPS板薄抹灰外墻外保溫系統(tǒng)是我國應(yīng)用最廣的一種外保溫系統(tǒng)。為了分析該保溫系統(tǒng)在哈爾濱地區(qū)使用過程中的保溫性能和熱濕傳遞變化情況，對其使用過程中的熱濕傳遞進(jìn)行模擬?？紤]墻體內(nèi)部水分的氣、固、液三相相變，建立以溫度梯度和體積含濕量為驅(qū)動(dòng)勢的熱濕耦合傳遞模型，并對模型進(jìn)行驗(yàn)證。結(jié)果表明：在模擬建筑材料的結(jié)融冰過程中，不能忽略氣液相變，這樣才能更好地模擬出墻體濕度的變化情況；混凝土墻體在使用過程中前3個(gè)月的干燥速度最快，之后變緩，冬季干燥速度最慢；在整個(gè)墻體結(jié)構(gòu)中，EPS保溫板的干燥速度最慢；隨著保溫系統(tǒng)內(nèi)部水分減少，其傳熱系數(shù)逐漸降低，而在冬季由于保溫層結(jié)冰，使得傳熱系數(shù)明顯上升。

多層建筑墻體；三相相變潛熱；熱濕耦合傳遞

Abstract: EPS plates thin plaster external wall insulation systems is one of the most widely used systems in our country. In order to analyze the thermal insulation properties and heat and moisture transfer of use process in Harbin, the heat and moisture coupling was simulated. Considering the three-phase moisture in the wall, the mathematical model of coupled heat and moisture transfer using temperature and moisture gradient as driving forces was build.The results show that gas-liquid phase transition can’t be ignored in order to simulate the changes of humidity when mathematical model is built. Drying speed is the fastest in the first three months and then slowed down.The drying speed in winter is slowest. In the whole structure of the wall, the drying rate of EPS board is the slowest. The heat transfer coefficient decreases gradually during the drying process, while it increases obviously due to the freezing of the insulation layer.

Keywords: multilayer wall; three-phase moisture latent; coupled heat and moisture transfer

研究建筑墻體的熱濕傳輸特性對于提高建筑保溫效果、發(fā)展節(jié)能建筑有著重要的意義。建筑墻體是一個(gè)典型的多層多孔結(jié)構(gòu)，其內(nèi)部孔洞中充滿濕空氣、液態(tài)水或冰，結(jié)構(gòu)內(nèi)的熱傳遞、濕傳遞及空氣滲透是一個(gè)典型的熱濕耦合過程[1]。幾十年來，國內(nèi)外學(xué)者對墻體熱濕耦合傳遞做了大量的研究，主要集中在墻體內(nèi)部水分氣液兩相相變情況下的熱濕特性研究[2-7]，而針對墻體內(nèi)部水分氣、固、液三相相變的熱濕耦合特性研究較少?？追布t等[8-9]針對新建圍護(hù)結(jié)構(gòu)熱濕耦合傳遞進(jìn)行模擬，建立凍融條件下液、固兩相的一維的熱濕傳遞模型，并結(jié)合實(shí)驗(yàn)進(jìn)行了研究，但是其采用建筑材料的毛細(xì)壓力擬合曲線與實(shí)際情況誤差較大，且在能量方程中忽略了相變潛熱對墻體內(nèi)部濕傳遞的影響。Wolfgang等[10]研究了凍融條件下多孔介質(zhì)的傳熱傳輸特性，著重研究了毛細(xì)壓力對濕傳遞的影響，但其研究的多孔介質(zhì)內(nèi)部含水量較高，有一定的局限性。呂恒林等[11]建立了以溫度和水蒸氣分壓力為驅(qū)動(dòng)勢的三相水分共存的多層墻體的熱濕耦合傳輸模型，其層與層交界處忽略接觸熱阻和濕傳遞阻，以水蒸氣分壓力作為水分的驅(qū)動(dòng)勢，得到的結(jié)果具有滯后性。文獻(xiàn)[8,11-12]在建立能量傳遞模型時(shí)認(rèn)為，濕傳遞引起的能量變化與相變產(chǎn)生的熱量相比很小，往往在能量方程中忽略濕傳遞項(xiàng)。但是由于EPS板薄抹灰外墻外保溫系統(tǒng)在使用初期，其內(nèi)部濕含量大，液水傳導(dǎo)率較大，故濕傳遞引起的能量變化不能忽略。

本文采用Whitaker體積平均理論對建筑墻體進(jìn)行處理，以溫度梯度和體積含濕量梯度為熱濕傳遞的驅(qū)動(dòng)勢，建立多層建筑墻體氣、固、液三相熱濕耦合傳遞模型，并在能量方程中保留濕傳遞項(xiàng)。對EPS板薄抹灰外墻外保溫系統(tǒng)在哈爾濱地區(qū)的使用過程進(jìn)行模擬，并分析模擬結(jié)果。

1 氣、固、液三相熱濕耦合傳遞模型

1.1 假設(shè)條件

① 建筑墻體為剛性多孔材料，結(jié)冰時(shí)不發(fā)生形變，各向同性；② 墻體內(nèi)部介質(zhì)連續(xù)，各相不被阻斷，局部熱力學(xué)平衡；③ 忽略濕分傳遞過程中的壓縮功和粘性耗散；④ 建筑材料內(nèi)部介質(zhì)的物性參數(shù)(比熱、導(dǎo)熱系數(shù)等)不隨溫度的變化而變化，空氣壓力恒定為大氣壓；⑤ 建筑材料內(nèi)部氣、固、液三相滿足蒸發(fā)平衡，并且水蒸氣分壓滿足開爾文方程；⑥ 固相不移動(dòng)，不蒸發(fā)，為剛性結(jié)構(gòu)；⑦ 空氣和水蒸氣組成的混合氣相滿足理想氣體方程；⑧ 由于建筑圍護(hù)結(jié)構(gòu)高度和寬度遠(yuǎn)大于厚度，所以本文的模型為一維模型。

1.2 質(zhì)量方程

根據(jù)質(zhì)量守恒定律，控制體內(nèi)各相組分質(zhì)量增量等于通過控制體邊界進(jìn)入的質(zhì)量與相變產(chǎn)生的質(zhì)量之和。

(1)

1.3 能量方程

相較于液固相變，氣固相變很小，因此本文忽略冰相的升華和水蒸氣的凝華。根據(jù)能量守恒定律，控制體內(nèi)總能量的變化量等于導(dǎo)熱傳入熱流、濕傳遞遷移帶入控制體的熱流以及內(nèi)熱源放熱熱流的總和，將相變吸放熱看作內(nèi)熱源的吸放熱。

(2)

ρtut(T)-ρsθsus(T)=ρlθlul(T)+

ρgθgug(T)+ρiθiui(T)

(3)

1.4 結(jié)冰點(diǎn)確定

徐宇工等[13]根據(jù)非均相系統(tǒng)中多相平衡的吉布斯相律推導(dǎo)出了冰水平衡時(shí)建筑材料內(nèi)部毛細(xì)壓力與結(jié)冰溫度的關(guān)系：

(4)

由此通過水分保持特性曲線可以得某一溫度下對應(yīng)的最大不結(jié)冰含濕量θl,max，并以此作為結(jié)冰發(fā)生點(diǎn)[14]。

1.5 邊界條件

1) 墻體內(nèi)邊界條件

質(zhì)量邊界為：

βv,in(Pvin-Pvsurf)=(Jl+Jv)|surf

(5)

能量邊界為：

αin(Tin-Tsurf)+hv-l(T0)βv,in(Pvin-Pvsurf)=

(6)

2) 墻體外邊界條件

質(zhì)量邊界：

βv,out(Pvsurf-Pvout)=[Jl+Jv]|surf

(7)

能量邊界：

αout(Tsurf-Tout)+hv-l(T0)βv,out(Pvsurf-Pvout)=

(8)

式中：αin、αout分別為室內(nèi)、外熱交換系數(shù)(W/m2K)；βv,in、βv,out分別為室內(nèi)、外水汽交換系數(shù)(W/m2K)；Tin、Tout分別為室內(nèi)、外環(huán)境溫度(K)；Tsurf為墻體表面溫度(K)；Pvin、Pvout、Pvsurf分別為室內(nèi)環(huán)境、室外環(huán)境、墻體表面的水蒸氣分壓力(Pa)。

2 求解方法

2.1 網(wǎng)格離散

本文采用內(nèi)節(jié)點(diǎn)法劃分網(wǎng)格。與外節(jié)點(diǎn)法相比，層與層界面處的節(jié)點(diǎn)只包含一種材料，在處理時(shí)不需要將兩種材料納入同一個(gè)節(jié)點(diǎn)中去，易于網(wǎng)格離散和程序編寫。然后通過控制容積法對質(zhì)量方程和能量方程進(jìn)行離散，時(shí)間格式為全隱格式。采用附加源項(xiàng)法對邊界處節(jié)點(diǎn)進(jìn)行處理，并使用耦合TDMA算法[15]求解。

2.2 層與層間界面處理

層與層邊界處節(jié)點(diǎn)分布如圖1所示。

圖1 層與層交界面處節(jié)點(diǎn)分布

在界面e處：建筑材料a通過界面e向建筑材料b進(jìn)行熱濕傳遞，這時(shí)有：

(9)

(10)

(11)

其中：

(12)

(13)

(14)

(15)

(16)

(17)

式中：Ple、Pve、Te分別為界面處液水壓力、水蒸氣分壓力、溫度。聯(lián)立上面的方程(9)-(17)可以得到界面e的導(dǎo)熱熱流量、水蒸氣流量、液水流量：

(18)

(19)

(20)

然后采用附加源項(xiàng)法對節(jié)點(diǎn)i和節(jié)點(diǎn)i+1處理即可。

2.3 熱濕物性參數(shù)選取及擬合

文獻(xiàn)[8]在計(jì)算模擬過程中選用文獻(xiàn)[16]提供的經(jīng)驗(yàn)公式:

(21)

式中：θSat為建筑材料毛細(xì)飽和體積含濕量;ΘHyg為相對濕度為95%時(shí)對應(yīng)的體積含濕量。顯然通過經(jīng)驗(yàn)公式計(jì)算得到的毛細(xì)壓力與體積含濕量是一個(gè)線性關(guān)系，這與實(shí)際情況不相符。Champs軟件數(shù)據(jù)庫中包含利用壓盤法實(shí)驗(yàn)獲得的EPS保溫板的水分保持特性曲線，如圖2所示。

從EPS保溫板經(jīng)驗(yàn)公式和實(shí)驗(yàn)數(shù)據(jù)對比發(fā)現(xiàn)：在較低含濕量的情況下，Champs軟件數(shù)據(jù)庫提供的實(shí)驗(yàn)數(shù)據(jù)可以更好地體現(xiàn)建筑材料毛細(xì)壓力的變化趨勢，所以本文在模擬計(jì)算過程中采用Champs軟件庫提供的實(shí)驗(yàn)數(shù)據(jù)，并利用3次樣條插值對其擬合。

圖2 EPS保溫板水分保持特性曲線

3 模型驗(yàn)證

為了驗(yàn)證模型的正確性，將模型計(jì)算結(jié)果與文獻(xiàn)[8]的實(shí)驗(yàn)結(jié)果進(jìn)行對比。文獻(xiàn)[8]中，測試了單層爐渣混凝土墻體在哈爾濱地區(qū)一年內(nèi)的溫濕度變化情況，將測點(diǎn)分別布置在墻體兩側(cè)表面和墻體中間。對比結(jié)果如圖3～4所示。

圖3 溫度對比

圖4 內(nèi)表面與中間切面相對濕度對比

圖5 相對濕度對比

通過對比結(jié)果發(fā)現(xiàn)：本文模擬得到的溫度結(jié)果與實(shí)驗(yàn)結(jié)果吻合度高；而內(nèi)表面處模擬得到的冬季(210～270 d)相對濕度與實(shí)驗(yàn)結(jié)果相比有較大的差別，同時(shí)模擬值始終高于實(shí)驗(yàn)值。這是由于冬季室內(nèi)外溫差大，無法保證實(shí)驗(yàn)墻體完全是一維傳熱傳濕，導(dǎo)致模擬結(jié)果與測量結(jié)果誤差較大[8]。從圖5可以看出：與文獻(xiàn)[8]相比，本文模擬得到的冬季外表面墻體相對濕度與實(shí)驗(yàn)值吻合較好。這是因?yàn)槎緣w右側(cè)結(jié)冰，而本文建立的結(jié)冰模型中考慮了氣液相變，與實(shí)際情況相符合，可以很好地模擬得到相對濕度的變化情況。

4 保溫系統(tǒng)熱濕耦合傳遞模擬

由于新建EPS板薄抹灰外墻外保溫系統(tǒng)內(nèi)部含濕量較高，墻體保溫效果較差，因此了解其內(nèi)部熱濕傳遞情況對改善其保溫性能有著重要意義。本文將對新建保溫墻體前3年的熱濕傳遞情況進(jìn)行模擬，保溫墻體結(jié)構(gòu)剖面如圖6所示。

圖6 墻體結(jié)構(gòu)剖面

4.1 初始條件

首先將內(nèi)側(cè)石灰水泥砂漿與混凝土層干燥1個(gè)月，之后加上EPS保溫層和外側(cè)石灰水泥砂漿抹面?？紤]到墻體內(nèi)含濕量較大，將各層初始條件設(shè)置為高含濕量。

內(nèi)側(cè)石灰水泥砂漿(10 mm)，t=0, 0≤x≤10 mm時(shí)：

T(x,0)=295 K，θl(x,0)=0.25 m3/m3

θi(x,0)=0 m3/m3

混凝土層(240 mm)，t=0, 10≤x≤250 mm 時(shí)：

T(x,0)=295 K，θl(x,0)=0.2 m3/m3

θi(x,0)=0 m3/m3

EPS保溫層(100 mm)，t=30 d, 250≤x≤350 mm 時(shí)：

T(x,0)=295 K，θl(x,0)=0.022 5 m3/m3

θi(x,0)=0 m3/m3

外側(cè)石灰水泥砂漿(10 mm)，t=30 d, 350≤x≤360 mm 時(shí)：

T(x,0)=295 K，θl(x,0)=0.25 m3/m3

θi(x,0)=0 m3/m3

4.2 室內(nèi)外環(huán)境參數(shù)

室內(nèi)環(huán)境溫度分別取Tin=295 K，濕度φin=50%，水汽交換系數(shù)βv,in=8.8×10-8s/m，熱交換系數(shù)αin=8.7 W/(m2K)；水汽交換系數(shù)βv,out=2.3×10-7s/m，熱交換系數(shù)αout=23.3 W/(m2K)。室外溫濕度條件為哈爾濱3年間的月平均溫濕度變化值，如圖7所示，模擬起始時(shí)刻從6月份開始。

圖7 哈爾濱全年月平均溫濕度變化

4.3 結(jié)果分析

由于EPS保溫板水液水傳導(dǎo)率和蒸汽傳導(dǎo)率都很小，混凝土內(nèi)的水分主要由室內(nèi)一側(cè)散出。圖8為混凝土層平均含濕量以及室外溫度隨時(shí)間變化圖。如圖所示，混凝土墻體在前3個(gè)月干燥速度較快，之后墻體含濕量變少，干燥速度變緩，冬季干燥速度最慢。這是因?yàn)楦稍锍跗冢簯B(tài)水含量較多，液水傳導(dǎo)率大，水分向液態(tài)水較少的室內(nèi)一側(cè)傳輸較快，干燥速度變快。隨著墻體含濕量的減少，水分傳輸主要以蒸汽的形式傳遞，在冬季室內(nèi)外溫差大，蒸汽向室內(nèi)一側(cè)的傳遞速度變慢，干燥速度降低。

圖8 混凝土層平均含濕量以及室外溫度隨時(shí)間的變化

EPS保溫板為不親水材料，液水傳導(dǎo)率很小，其內(nèi)部的濕份傳輸方式主要為蒸汽滲透?；炷翆咏?jīng)過1個(gè)月的干燥，與保溫板接觸一側(cè)相對濕度低，同時(shí)保溫板右側(cè)砂漿干燥速度快，使得整個(gè)墻體中保溫板內(nèi)部相對濕度最大。保溫板內(nèi)部濕份將透過相鄰墻體結(jié)構(gòu)散濕到環(huán)境中，這一過程中保溫板內(nèi)部的含濕量不斷降低，如圖9所示。另外，在冬季保溫層內(nèi)部結(jié)冰，使得液水含量突降。

圖10為墻體換熱系數(shù)和保溫板結(jié)冰量隨時(shí)間的變化圖。如圖所示，隨著墻體的不斷干燥，墻體傳熱系數(shù)不斷減小。但是在冬季保溫層結(jié)冰，使得墻體傳熱系數(shù)升高，結(jié)冰量越多則墻體保溫性能越差，因此避免保溫層結(jié)冰是提高墻體保溫措施的重要措施。

圖9 EPS保溫板平均含濕量隨時(shí)間的變化

圖10 墻體換熱系數(shù)以及保溫層結(jié)冰量隨時(shí)間變化

5 結(jié)論

1) 在模擬建筑材料的結(jié)融冰過程中，不能忽略氣液相變，這樣才能更好地模擬出墻體濕度的變化情況。

2) 混凝土墻體在使用過程中前3個(gè)月的干燥速度最快，之后變緩，冬季干燥速度最慢。

3) 在整個(gè)墻體結(jié)構(gòu)中，EPS保溫板的干燥速度最慢，在使用前要防止保溫層受潮。

4) 隨著保溫系統(tǒng)內(nèi)部水分減少，其傳熱系數(shù)逐漸降低，而在冬季由于保溫層結(jié)冰，使得傳熱系數(shù)明顯上升。

[1] 孔凡紅.節(jié)能建筑圍護(hù)結(jié)構(gòu)熱質(zhì)耦合傳遞的影響綜述[J].太陽能學(xué)報(bào),2012(S1):91-97.

[2] LIU Y,WANG Y,WANG D,et al.Effect of moisture transfer on internal surface temperature[J].Energy & Buildings,2013,60(6):83-91.

[4] BART M,MOISSETTE S,OUMEZIANE Y A.A hybrid modeling for transient heat and moisture transfer in porous hygroscopic building materials[J].Special Topics & Reviews in Porous Media,2013,4(3):273-286.

[5] 張景欣,郭興國,陳友明,等.墻體內(nèi)熱濕耦合過程的時(shí)域遞歸展開算法[J].土木建筑與環(huán)境工程,2015(6):147-152.

[6] 黃建恩,呂恒林,馮偉.多層墻體熱濕耦合傳遞模型及驗(yàn)證[J].土木建筑與環(huán)境工程,2015,37(1):18-22.

[7] VASILYEV G P,LICHMAN V A,PESKOV N V,et al.Simulation of heat and moisture transfer in a multiplex structure[J].Energy & Buildings,2015,86:803-807.

[8] 孔凡紅.新建建筑圍護(hù)結(jié)構(gòu)干燥特性及其影響研究[D].哈爾濱：哈爾濱工業(yè)大學(xué),2008.

[9] KONG F,ZHANG Q.Effect of heat and mass coupled transfer combined with freezing process on building exterior envelope[J].Energy & Buildings,2013,62(62):486-495.

[10] BLO?FELD W M,BLUHM J,RICKEN T.Simulation of Freezing and Thawing Processes with Capillary Effects in fluid filled porous media[J].Pamm,2012,12(1):365-366.

[11] 呂恒林,黃建恩,馮偉,等.三相水分共存的多層墻體熱濕耦合傳遞模型[J].建筑材料學(xué)報(bào),2016,19(2):310-316.

[12] BLUHM J,RICKEN T,BLO?FELD W M.Simulation of freeze-thaw-cycles in liquid and gas saturated porous media[J].Pamm,2015,10(1):359-360.

[13] 徐宇工,李笑.多孔建筑材料凍融現(xiàn)象研究中孔隙水凍結(jié)點(diǎn)的確定方法[J].哈爾濱建筑大學(xué)學(xué)報(bào),2001,34(3):91-95.

[14] KURYLYK B L,WATANABE K.The mathematical representation of freezing and thawing processes in variably-saturated,non-deformable soils[J].Advances in Water Resources,2013,60(60):160-177.

[15] MENDES N,PHILIPPI P C,LAMBERTS R.A new mathematical method to solve highly coupled equations of heat and mass transfer in porous media[J].International Journal of Heat & Mass Transfer,2002,45(3):509-518.

(責(zé)任編輯林 芳)

SimulationonCoupledHeatandMoistureTransferforEPSPlatesThinPlasterExternalWallSystems

LIU Jianchang1, LI Longjian1, SHEN Xianwen1, FENG Ya2

(1.College of Power Engineering, Chongqing University, Chongqing 400030, China; 2.China Southwest Architectural Design and Research Institute, Chengdu 610042, China)

2017-03-27

國家自然科學(xué)基金資助項(xiàng)目(51268478)

劉建昌(1990—)，男，碩士研究生，主要從事多層建筑熱濕特性研究，E-mail：1142890262@qq.com； 通訊作者 李隆鍵(1946—)，男，教授，博士生導(dǎo)師，E-mail：longjian@cqu.edu.cn。

劉建昌，李隆鍵，申憲文，等.EPS板薄抹灰外墻外保溫系統(tǒng)熱濕耦合傳遞的數(shù)值模擬[J].重慶理工大學(xué)學(xué)報(bào)(自然科學(xué))，2017(9):90-96.

formatLIU Jianchang, LI Longjian, SHEN Xianwen, et al.Simulation on Coupled Heat and Moisture Transfer for EPS Plates Thin Plaster External Wall Systems[J].Journal of Chongqing University of Technology(Natural Science)，2017(9):90-96.

10.3969/j.issn.1674-8425(z).2017.09.015

TU111.4

A

1674-8425(2017)09-0090-07