劉顯茜，趙振超，鄒三全，張雪波

（昆明理工大學(xué)機(jī)電工程學(xué)院，云南昆明 650500）

0 引言

2017 年建筑能源消費(fèi)占全國(guó)總能源消費(fèi)約21.10%，換算為9.47 億標(biāo)準(zhǔn)煤，建筑領(lǐng)域節(jié)能減排任務(wù)十分艱巨。建筑墻體的耗熱量約占建筑能耗的25%，新型節(jié)能墻體材料的開發(fā)與推廣應(yīng)用是實(shí)現(xiàn)建筑節(jié)能與降低建筑能源消耗的一項(xiàng)重要措施。

新型墻體材料開發(fā)與應(yīng)用發(fā)展很快，相繼有蒸壓加氣混凝土砌塊、再生混凝土砌塊、復(fù)合混凝土砌塊等在工程中應(yīng)用，復(fù)合混凝土砌塊研究在原材料、受力性能、結(jié)構(gòu)及熱工性能方面取得一些成果。于秋波等對(duì)煤矸石燒結(jié)的復(fù)合保溫砌塊進(jìn)行砌體受壓性能研究，研究結(jié)果表明煤矸石燒結(jié)的復(fù)合混凝土砌塊的局部破壞會(huì)導(dǎo)致砌體試件整體失穩(wěn)；白玲對(duì)由新型陶粒復(fù)合自保溫砌塊砌成的墻體進(jìn)行整體結(jié)構(gòu)與熱力學(xué)性能研究，結(jié)果證明新型陶粒自保溫砌體結(jié)構(gòu)穩(wěn)定，保溫性能符合國(guó)家標(biāo)準(zhǔn)，選材低碳環(huán)保；楊召通等采用有限元數(shù)值模擬和加權(quán)求和法對(duì)混凝土砌塊進(jìn)行多目標(biāo)優(yōu)化，研究表明相同的砌塊隨著孔排的交錯(cuò)程度增加其保溫隔熱性能也增加；馮偉運(yùn)用Kelvins 定律將水蒸氣和液態(tài)水傳遞量轉(zhuǎn)變?yōu)橐韵鄬?duì)濕度梯度為驅(qū)動(dòng)勢(shì)的統(tǒng)一函數(shù)，并利用多物理場(chǎng)耦合軟件對(duì)復(fù)合混凝土砌塊內(nèi)部的熱濕耦合傳遞進(jìn)行數(shù)值模擬計(jì)算，通過實(shí)驗(yàn)驗(yàn)證數(shù)值模擬計(jì)算的準(zhǔn)確性與適用性。

前人對(duì)復(fù)合砌塊的研究成果很多，但對(duì)復(fù)合混凝土砌塊內(nèi)部熱濕耦合傳熱規(guī)律研究尚不深入，特別是外部非穩(wěn)態(tài)濕度因素對(duì)砌塊內(nèi)部各位置傳熱的影響研究較少。

采用數(shù)值模擬仿真的方法對(duì)建筑墻體進(jìn)行研究越來越多。與數(shù)值模擬仿真相比，實(shí)驗(yàn)研究的方法存在材料昂貴、時(shí)間成本高、容易受外部環(huán)境影響等缺點(diǎn)。鑒于此，本文基于多孔介質(zhì)瞬態(tài)熱濕耦合傳遞基本理論，結(jié)合COMSOL Multiphysics多物理場(chǎng)耦合軟件，采用數(shù)值模擬仿真方法模擬復(fù)合混凝土砌塊的瞬態(tài)熱濕耦合傳遞過程，研究復(fù)合混凝土砌塊內(nèi)部的傳熱規(guī)律。

1 控制方程

復(fù)合混凝土砌塊的主體材料、填充保溫材料、固化后砂漿均為多孔介質(zhì)，可以運(yùn)用多孔介質(zhì)內(nèi)部瞬態(tài)熱濕耦合模型來處理砌塊內(nèi)部的熱量和水分傳遞，采用最基本的Fick 和Darcy 定律推導(dǎo)出熱量和水分在多孔砌塊內(nèi)部的瞬態(tài)耦合傳遞方程。常用的熱濕耦合傳遞驅(qū)動(dòng)勢(shì)有溫度梯度、相對(duì)濕度梯度、水蒸氣分壓力等。劉向偉、呂恒林等采用的驅(qū)動(dòng)勢(shì)為溫度梯度與相對(duì)濕度梯度共同作用，本文采用與之相同的方法并對(duì)方程進(jìn)行簡(jiǎn)化，以提高方程運(yùn)用場(chǎng)景。

1.1 水分傳遞方程

在孔隙中，水分會(huì)以水蒸氣、液態(tài)水和冰的形態(tài)存在。本文采用的氣候數(shù)據(jù)為夏熱冬冷地區(qū)，年平均氣溫在0℃以上，水分只能以氣態(tài)和液態(tài)形式傳遞，主要傳遞形式是蒸汽擴(kuò)散、毛細(xì)吸附或兩種形式共存。

由質(zhì)量守恒定律得水分傳遞方程：

式中，

w

是含水量函數(shù)，單位kg∕kg；

g

是蒸汽擴(kuò)散通量，單位kg∕（m·s）；

g

是液態(tài)水通量，單位kg∕（m·s）。

通過Fick 定律和Darcy 定律來描述，則公式（1）變?yōu)椋?/p>

式中，

δ

是水蒸氣滲透率，單位kg∕（m·Pa·s）；

p

是部分水蒸氣壓力，單位Pa；

K

是液態(tài)水滲透率，單位kg∕（m·Pa·s）；

p

是毛細(xì)管壓力，單位Pa；

p

=

φP

，

P

是飽和水蒸氣壓力，是一個(gè)關(guān)于溫度

T

（

k

）的函數(shù)。由于含水量

w

在兩個(gè)材料的界面處是不連續(xù)的，所以要轉(zhuǎn)化為溫度與相對(duì)濕度

φ

（%）的雙驅(qū)動(dòng)函數(shù)：

1.2 熱傳遞方程

因流入控制單元體的能量減去流出控制單元體的能量為流入控制單元體的凈能量，等于控制單元體內(nèi)焓的變化，所以能量守恒可表達(dá)為：

式中，

ρ

是干材料密度，單位kg∕m；

c

是干材料比熱，單位J∕（kg·K）；

c

是液態(tài)水比熱，單位J∕（kg·K）；

q

是熱傳導(dǎo)熱流密度，單位W∕m；

q

是熱對(duì)流熱流密度，單位W∕m。

式中，

k

是材料導(dǎo)熱系數(shù)，單位W∕（m·K）；

c

是干空氣的比熱，單位J∕（kg·K）；

h

是水蒸氣汽化潛熱，單位J∕kg。

1.3 空氣流動(dòng)方程

根據(jù)Poiseuille 定律，通過多孔介質(zhì)砌塊的空氣通量可表示為：

式中，

k

是空氣通過多孔介質(zhì)砌塊的滲透率，單位kg∕（m·Pa·s），是空氣流動(dòng)速率與壓力梯度的比值。

在建筑物理應(yīng)用中，空氣常被認(rèn)為是不可壓縮流體，在實(shí)際應(yīng)用中，空氣通過多孔介質(zhì)的孔隙時(shí)會(huì)呈現(xiàn)出極低的流速、較低的壓力以及溫度變化不明顯特性。因此，

1.4 邊界條件方程

傳質(zhì)邊界條件：

傳熱邊界條件：

式中，

β

是水分遷移系數(shù)，單位s∕m；

φ

是相對(duì)濕度；

h

是傳熱系數(shù)，單位W∕（m·K）；

p

是飽和蒸汽壓力，單位Pa；

T

是溫度，單位K；下標(biāo)中：

surf

表示墻體表面，

e

表示外部，

i

表示內(nèi)部。

2 模擬對(duì)象

通過查詢國(guó)家相關(guān)標(biāo)準(zhǔn)及調(diào)研市場(chǎng)相關(guān)產(chǎn)品，以市場(chǎng)上常見390mm×240mm×190mm 尺寸的復(fù)合混凝土砌塊為研究對(duì)象，圓倒角設(shè)為5 mm，研究復(fù)合混凝土砌塊在穩(wěn)定條件及非穩(wěn)態(tài)環(huán)境下的熱濕耦合傳遞，并通過數(shù)值模擬顯示計(jì)算結(jié)果。

根據(jù)熱濕空氣耦合控制方程，進(jìn)行簡(jiǎn)化后便可在非穩(wěn)態(tài)環(huán)境下對(duì)砌塊的熱濕耦合傳遞進(jìn)行數(shù)值模擬，得出研究對(duì)象內(nèi)部溫度的分布情況。對(duì)模擬做出如下假設(shè)：砌塊砂漿、保溫材料均與砌塊實(shí)體緊密接觸，不考慮接觸熱阻，不考慮熱輻射，熱量沿砌塊的厚度方向傳遞，其余面為熱濕絕緣狀態(tài)，不考慮風(fēng)速及太陽(yáng)光輻射。

2.1 COMSOL Multiphysics

多物理耦合仿真軟件COMSOL Multiphysics 以有限元法作為求解偏微分問題基礎(chǔ)，當(dāng)某種物理現(xiàn)象可用偏微分方程描述便可使用該軟件進(jìn)行模擬仿真。COMSOL Multiphysics 集成了流體流動(dòng)、熱傳導(dǎo)、結(jié)構(gòu)力學(xué)、電磁分析等多個(gè)物理場(chǎng)模塊，并以高效的計(jì)算能力和杰出的多場(chǎng)雙向耦合分析能力實(shí)現(xiàn)高精度的模擬仿真。本文采用該軟件對(duì)砌塊的熱濕傳遞流體流動(dòng)過程進(jìn)行仿真。

2.2 參數(shù)設(shè)置

模擬采用的幾何模型俯視圖如圖1（彩圖掃OSID 碼可見，下同）所示。砌塊有6個(gè)大小不一的排孔，中間填充保溫材料以提高砌塊整體的保溫性能，但不會(huì)影響砌塊的實(shí)用性以及安全性。兩端分別有一個(gè)開放的排孔，方便砌墻時(shí)與其他砌塊連接，使整個(gè)墻體更加牢固。

Fig.1 Top view of block（Unit：mm）圖1 砌塊俯視圖（單位：mm）

幾何模型建立后，因計(jì)算域包含多個(gè)不同的材料，因此需要設(shè)定不同材料傳熱所需要的屬性。復(fù)合混凝土砌塊主要材料成分配比為：煤矸石：頁(yè)巖：稻殼=0.85：0.1：0.05，保溫材料材料選用聚苯乙烯泡沫塑料（EPS）。

材料的基本物性參數(shù)如表1 所示。

其中，砌塊主體的蒸汽滲透系數(shù)與含水量相關(guān)，函數(shù)公式如下：

Table 1 Basic physical property parameters表1 基本物性參數(shù)

2.3 網(wǎng)格劃分及網(wǎng)格無關(guān)性驗(yàn)證

采用多物理場(chǎng)耦合軟件自帶的網(wǎng)格生成器進(jìn)行自由三角形網(wǎng)格劃分，進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證，不同網(wǎng)格大小的同一模型計(jì)算如圖2 所示。當(dāng)最小網(wǎng)格大小為4 mm 時(shí)，域網(wǎng)格大小不會(huì)影響仿真結(jié)果，因此域網(wǎng)格數(shù)量最終確定為68 534個(gè)。

Fig.2 Temperature changes with time for different grid size圖2 不同網(wǎng)格大小溫度隨時(shí)間變化

2.4 仿真模型設(shè)置

為驗(yàn)證對(duì)混凝土砌塊進(jìn)行熱濕耦合傳熱的必要性以及減少不必要的工作，本文進(jìn)行兩組數(shù)值模擬：第一組進(jìn)行穩(wěn)定邊界條件下熱濕耦合傳熱分析，第二組進(jìn)行非穩(wěn)態(tài)邊界條件下的熱濕耦合傳熱分析?；炷疗鰤K結(jié)構(gòu)可簡(jiǎn)化為二維模型以減少計(jì)算量。

2.5 邊界條件與初始條件

（1）根據(jù)GB 50176-2016《民用建筑熱工設(shè)計(jì)規(guī)范》規(guī)定，穩(wěn)定邊界條件下的砌塊內(nèi)外表面換熱系數(shù)分別為8.7W∕（m·K）和23W∕（m·K）。為使砌塊兩端熱量進(jìn)行傳遞，將兩端設(shè)為溫度邊界條件，無相對(duì)濕度梯度，溫度差為40K，所以將外溫度設(shè)為263.15K，內(nèi)溫度設(shè)置為303.15K；內(nèi)外邊界相對(duì)濕度設(shè)置與砌塊初始相對(duì)濕度相同；其余邊界設(shè)置為第二邊界條件，熱流密度設(shè)置為零。用于模擬的砌塊初始溫度與冷邊界溫度相等。用相對(duì)濕度描述復(fù)合混凝土砌塊內(nèi)部的水分含量，復(fù)合混凝土砌塊內(nèi)部的初始相對(duì)濕度分別為0、10%、30%、50%、70%、90%。驅(qū)動(dòng)勢(shì)采用溫度梯度和相對(duì)濕度梯度共同作用的雙驅(qū)動(dòng)模型，模擬穩(wěn)定邊界條件下砌塊內(nèi)部的熱濕耦合傳遞。

（2）根據(jù)《民用建筑熱工設(shè)計(jì)規(guī)范》與《建筑用標(biāo)準(zhǔn)氣象數(shù)據(jù)手冊(cè)》中的規(guī)定設(shè)置非穩(wěn)態(tài)環(huán)境下的邊界條件，外部邊界條件選用《建筑用標(biāo)準(zhǔn)氣象數(shù)據(jù)手冊(cè)》長(zhǎng)沙地區(qū)累年1 月份的環(huán)境干球溫度平均值和環(huán)境相對(duì)濕度；室內(nèi)環(huán)境因選用為非取暖房，所以室內(nèi)邊界條件設(shè)為285.15 K、60%和285.15 K、30%。

3 模擬結(jié)果分析

3.1 穩(wěn)定邊界條件下的結(jié)果分析

在砌塊內(nèi)部取兩點(diǎn)，如圖1 所示，坐標(biāo)分別為（195，52.5）、（195，230）。砌塊內(nèi)部?jī)牲c(diǎn)處的溫度隨時(shí)間變化如圖3、圖4 所示。C 代表熱濕耦合傳熱；0-5 代表進(jìn)行熱濕耦合傳遞數(shù)值模擬的復(fù)合混凝土砌塊相對(duì)濕度分別為0、10%、30%、50%、70%、90%。

Fig.3 Temperature change at point 1inside the block圖3 砌塊內(nèi)部點(diǎn)1 處的溫度變化

Fig.4 Temperature change at point 2 inside the block圖4 砌塊內(nèi)部點(diǎn)2 處的溫度變化

由圖3 知，砌塊內(nèi)部點(diǎn)1 處6個(gè)小組模型48h 溫度都不隨時(shí)間變化，說明達(dá)到了平衡。C0、C1、C2、C3 小組分別在45h 達(dá)到平衡，達(dá)到平衡時(shí)的溫度是相近的，為271.6 K。隨著砌塊相對(duì)濕度的增加，砌塊達(dá)到平衡的時(shí)間不斷縮短；C4、C5 小組在25h 達(dá)到平衡狀態(tài)；砌塊相對(duì)濕度為90%時(shí)，砌塊在點(diǎn)1 處平衡時(shí)溫度為275 K。

通過圖4 可以看出，相較于點(diǎn)1，點(diǎn)2 距離高溫邊界較近，溫度以一個(gè)非?？斓乃俣壬仙黠@高于點(diǎn)1 處溫度上升的速度，在很短的時(shí)間內(nèi)點(diǎn)2 處的溫度就達(dá)到平衡狀態(tài)。C0、C1、C2、C3 小組在達(dá)到平衡時(shí)的溫度都為300K；C4、C5小組達(dá)到平衡時(shí)溫度分別為297.8 K 和292.9K。

所以在相對(duì)濕度低于50%時(shí)砌塊傳熱速度不會(huì)有太大影響，而當(dāng)相對(duì)濕度高于50%時(shí)，砌塊在進(jìn)行熱濕耦合達(dá)到平衡狀態(tài)的速度明顯高于砌塊相對(duì)濕度低于50%進(jìn)行熱濕耦合達(dá)到平衡狀態(tài)的速度，說明砌塊中的含水量有利于溫度傳遞，砌塊中的含水量越高，傳熱的速率越快。

3.2 非穩(wěn)態(tài)邊界條件結(jié)果分析

在砌塊長(zhǎng)度上取距離橫截面195 mm 位置、厚度方向上取距離外墻面10 mm、120 mm、230 mm 位置，初始相對(duì)濕度30%時(shí)溫度隨時(shí)間變化趨勢(shì)如圖5 所示，Amth1 代表外部非穩(wěn)態(tài)環(huán)境氣象數(shù)據(jù)。

Fig.5 Temperature changes with time at different positions along the thickness direction圖5 沿厚度方向不同位置溫度隨時(shí)間變化曲線

Fig.6 Temperature changes with time of different humidity along the thickness direction圖6 沿厚度方向不同濕度的溫度隨時(shí)間變化曲線

圖5 顯示，同一種材料的同一截面不同點(diǎn)溫度有一定差別，而且這種差別不僅在數(shù)值模擬的不同時(shí)刻有所體現(xiàn)，在不同濕度時(shí)也呈現(xiàn)出一定的差別。如圖6 所示，實(shí)心代表相對(duì)濕度為30%，空心代表相對(duì)濕度為60%。因此，混凝土砌塊的熱濕耦合傳遞受到不同材料物理性質(zhì)差異的影響，呈現(xiàn)出不同的影響程度。砌塊沿厚度方向距離外墻面10 mm 的位置受室外溫度起伏的影響最大，所呈現(xiàn)的波形與室外溫度波形相近。隨著離外墻面距離的增加，受室外溫度起伏的影響越小。距離內(nèi)墻面10 mm 的位置幾乎不受室外溫度影響，溫度變化曲線較為平緩。隨著距外墻面距離的增加，砌塊內(nèi)部溫度波動(dòng)有一定的滯后性。

4 結(jié)語(yǔ)

由多孔介質(zhì)內(nèi)部瞬態(tài)熱濕耦合傳遞拓展到復(fù)合混凝土砌塊內(nèi)部熱濕耦合傳遞，根據(jù)數(shù)值模擬所需計(jì)算參數(shù)的設(shè)置方法，本文運(yùn)用多物理場(chǎng)耦合仿真軟件對(duì)穩(wěn)定邊界及非穩(wěn)態(tài)邊界條件下砌塊的熱濕耦合傳熱規(guī)律進(jìn)行數(shù)值模擬研究，得出以下結(jié)論：

（1）將運(yùn)用在多孔介質(zhì)內(nèi)部瞬態(tài)熱濕耦合的傳遞方程進(jìn)行推導(dǎo)和簡(jiǎn)化，使之可以適用于砌塊內(nèi)部熱濕耦合傳遞；對(duì)二維模型進(jìn)行計(jì)算，節(jié)約了更多的計(jì)算時(shí)間與成本。

（2）在穩(wěn)定邊界條件下，砌塊內(nèi)部的水分含量會(huì)影響砌塊的傳熱。隨著水分含量的增加，砌塊的傳熱速度也會(huì)增加；與不考慮傳濕相比，熱濕耦合傳遞對(duì)砌塊內(nèi)部傳熱有一定影響。

（3）在室外非穩(wěn)態(tài)邊界條件下，越靠近室外的位置受室外環(huán)境波動(dòng)的影響越大；內(nèi)部各位置的溫度變化趨勢(shì)存在滯后，越靠近室內(nèi)滯后時(shí)間越長(zhǎng)。

墻體含水量增加會(huì)導(dǎo)致墻體傳熱性能增強(qiáng)，降低了墻體保溫性能。因此，在實(shí)際工程中要采取一定措施減少墻體內(nèi)部的含水量。不僅要使用保溫性能、絕濕性能好的墻體材料，也要輔以相應(yīng)的墻面材料。而本文只針對(duì)墻體砌塊展開研究，具有一定的局限性，后續(xù)將會(huì)以墻體構(gòu)件作為參照物，構(gòu)建更加符合實(shí)際的模型，以增加研究的應(yīng)用價(jià)值。