羅戴維 劉加平 劉大龍

摘 要：極端熱濕地區(qū)常年高溫多雨，為了準(zhǔn)確預(yù)測圍護(hù)結(jié)構(gòu)內(nèi)的溫度和濕度分布，模擬高溫、高濕和高太陽輻射對圍護(hù)結(jié)構(gòu)的影響，建立了適用于極端熱濕地區(qū)的圍護(hù)結(jié)構(gòu)熱濕耦合傳遞模型，在邊界條件中加入太陽輻射和雨水負(fù)荷對傳熱傳濕的影響，考慮了隨材料含濕量不斷變化的材料物性參數(shù)。為了求解模型，提出了基于計(jì)算機(jī)軟件COMSOL的模型求解方法，利用MATLAB計(jì)算方程組不斷變化的系數(shù)，并確定了軟件容差和網(wǎng)格的設(shè)置條件：相對容差推薦設(shè)定在0.001，絕對容差設(shè)定為0.000 1，網(wǎng)格按極細(xì)化劃分。通過對比HAMSTAD標(biāo)準(zhǔn)實(shí)例驗(yàn)證了模型的準(zhǔn)確性。

關(guān)鍵詞：極端熱濕地區(qū)；熱濕耦合傳遞；太陽輻射；雨水負(fù)荷

中圖分類號：TU111.4

文獻(xiàn)標(biāo)志碼：A文章編號：1674-4764（2018）04-0036-06

Abstract：Extreme heat-moisture climate area are always high temperature and rainy. In order to accurately predict the temperature and humidity distribution in the envelope， and simulate the effects of high temperature， high humidity and high solar radiation on the envelope， a coupled heat and moisture transfer model which takes the climate condition of extreme heat-moisture climate area into consideration is developed. The influence of solar radiation and rain load on heat and moisture transfer is considered in boundary conditions. And the changes of moisture content are taken into account in the material parameters. And a solution method based on computer soft，COMSOL and MATLAB，is proposed to solve the coupled heat and moisture transfer model simply. The COMSOL recommended setting is as follow： the relative tolerance is 0.001； the absolute tolerance is 0.0001； the mesh must be set to very fine. The numerical results of this model agree well with HAMSTAD benchmarks.

Keywords：extreme heat-moisture climate area； coupled heat and moisture transfer； solar radiation； rain load

極端熱濕地區(qū)常年處在高溫高濕狀態(tài)，這種氣候條件會給圍護(hù)結(jié)構(gòu)隔熱隔濕帶來極大負(fù)擔(dān)。以南海西沙群島為例，其日照豐富，太陽輻射強(qiáng)烈，年平均氣溫26～27 ℃，6月份平均氣溫為28.9 ℃，一月份平均氣溫為22.9 ℃，總降雨量約為1 500 mm。南海全年氣溫高、溫差小，雨量豐沛、時(shí)空分布不均，強(qiáng)風(fēng)、臺風(fēng)、暴雨、大霧等災(zāi)害性天氣頻發(fā)，是北半球天氣氣候變化最敏感的地區(qū)之一，屬于典型的極端熱濕氣候區(qū)[1]。這樣的氣候條件將造成其傳熱、傳濕的特性發(fā)生改變，導(dǎo)致建筑空調(diào)耗電問題異常突出。因此，傳熱和傳濕過程必須被同時(shí)考慮，這樣才能更準(zhǔn)確地預(yù)測圍護(hù)結(jié)構(gòu)的熱濕遷移以及室內(nèi)舒適度。

用來描述圍護(hù)結(jié)構(gòu)熱濕耦合傳遞的數(shù)學(xué)模型多種多樣，Hou等[2]、Kong等[3]使用了水蒸汽壓力為濕驅(qū)動勢的熱濕耦合傳遞模型，F(xiàn)ang等[4]、Tomá等[5]、Min等[6]REF_Ref489776463＼r＼h＼*MERGEFORMAT、kerget等[7]使用相對濕度為濕驅(qū)動勢的熱濕耦合傳遞模型，但這些模型的基礎(chǔ)理論都相通，具體使用哪種模型需要根據(jù)實(shí)際需求來定。此外，劉向偉等[8]、張景欣等[9]、王瑩瑩[10]針對熱濕耦合模型進(jìn)行了多方面的研究，研究對象包括空調(diào)房間、計(jì)算方法、室內(nèi)環(huán)境質(zhì)量。然而，傳統(tǒng)的熱濕耦合模型針對極端熱濕地區(qū)的應(yīng)用較為稀少，尤其針對南海地區(qū)的研究更是近乎為零。筆者以溫度和毛細(xì)壓力為驅(qū)動勢來預(yù)測圍護(hù)結(jié)構(gòu)的熱量和水分傳遞，并根據(jù)極端熱濕地區(qū)的高溫高輻射多雨的氣候條件對模型進(jìn)行了優(yōu)化。該模型充分考慮了極端熱濕地區(qū)計(jì)算參數(shù)變物性的需求，使其能更加準(zhǔn)確地描述極端熱濕地區(qū)的特點(diǎn)。

1 熱濕耦合模型

濕在圍護(hù)結(jié)構(gòu)中傳遞的過程呈現(xiàn)出兩種不同相態(tài)：液相和氣相，在傳統(tǒng)計(jì)算中，液態(tài)水的傳遞一般被忽略。然而，由于極端熱濕地區(qū)的氣候特征，室外相對濕度常年在80%以上，液態(tài)水的傳輸在模型計(jì)算中不可忽略或簡化。此外，蒸汽傳輸可以分為兩部分：擴(kuò)散和對流。蒸汽擴(kuò)散由蒸汽壓梯度驅(qū)動，其系數(shù)為蒸汽滲透系數(shù)，單位為kg/（m·Pa·s）；蒸汽對流由空氣流動引起，由于極端熱濕地區(qū)室外常年處于高溫狀態(tài)，室內(nèi)很容易保持負(fù)壓狀態(tài)，空氣壓差的方向是由室內(nèi)向室外，這將減輕熱負(fù)荷。蒸汽對流過程計(jì)算較為復(fù)雜，因此，在熱工設(shè)計(jì)計(jì)算上可以忽略[11]，這樣一來，計(jì)算結(jié)果將偏安全。液態(tài)水的傳輸則經(jīng)常用毛細(xì)壓力作為驅(qū)動勢來計(jì)算，其系數(shù)為液態(tài)水傳導(dǎo)系數(shù)，單位為kg/（m·Pa·s）。

1.1 傳濕模型

極端熱濕地區(qū)圍護(hù)結(jié)構(gòu)傳濕過程是一個(gè)包含水蒸氣和液態(tài)水的傳質(zhì)過程，通過建筑圍護(hù)結(jié)構(gòu)中的蒸汽和液體遷移可以表達(dá)為[12]

2 極端熱濕地區(qū)模型的特殊性

2.1 邊界條件

任何數(shù)學(xué)模型的仿真模型都離不開邊界條件的正確設(shè)定。極端熱濕地區(qū)常年高溫多雨，相比其他地區(qū)，太陽輻射和雨水在邊界條件中的計(jì)算不可缺少。

2.2 太陽輻射與雨水的計(jì)算

圍護(hù)結(jié)構(gòu)外表面的對流換熱系數(shù)和輻射換熱系數(shù)均與外界的氣候因素（風(fēng)速、天空溫度、環(huán)境溫度）以及外表面的溫度有關(guān)，他們的值可以分別列式計(jì)算。但為了計(jì)算工作的簡便，在熱工設(shè)計(jì)計(jì)算中一般采用總換熱系數(shù)hα來表示，用于統(tǒng)一計(jì)算對流和輻射總換熱量。

中國熱工規(guī)范[16]規(guī)定，圍護(hù)結(jié)構(gòu)內(nèi)表面的總換熱系數(shù)在冬季是8.72 W/（m2·K），夏季為8.75 W/（m2·K），外表面的總換熱系數(shù)冬季為23.3 W/（m2·K），夏季為18.6 W/（m2·K）。該數(shù)值是按建筑物周圍風(fēng)速為3 m/s時(shí)實(shí)測所得。

根據(jù)圍護(hù)結(jié)構(gòu)外表面的熱平衡方程式，可以得出由圍護(hù)結(jié)構(gòu)外表面向室內(nèi)側(cè)的傳熱量q0。

2.3 材料的物性參數(shù)

極端熱濕地區(qū)建筑圍護(hù)材料常年接觸潮濕的空氣，因此，其材料的物性參數(shù)將與平常狀態(tài)下大不相同。圍護(hù)結(jié)構(gòu)材料的導(dǎo)熱系數(shù)、蒸汽滲透系數(shù)和液態(tài)水滲透系數(shù)是隨材料含濕量變化的參數(shù)，在計(jì)算中需要注意。

3 模型的數(shù)值計(jì)算與驗(yàn)證

采用COMSOL Multiphysics來求解模型。COMSOL Multiphysics會用有限元法自動對控制方程及相應(yīng)的邊界條件進(jìn)行離散，并用數(shù)值求解器對離散后的代數(shù)方程組求解，可以大大簡化求解的過程，方便研究人員將注意力集中到物理模型的研究。此外，由于熱濕耦合傳遞是一個(gè)瞬態(tài)過程，模型的相關(guān)參數(shù)隨著時(shí)間不斷變化，為了反映這一點(diǎn)，用MATLAB編寫了系數(shù)計(jì)算程序。

3.1 HAMSTAD基準(zhǔn)案例

為了提高建筑圍護(hù)結(jié)構(gòu)的性能，改進(jìn)建筑圍護(hù)結(jié)構(gòu)熱濕耦合傳遞計(jì)算模型，歐盟在2000年發(fā)起了一項(xiàng)研究項(xiàng)目HAMSTAD（Heat， Air and Moisture Standards Development）[18]，其目的在于確定多孔材料的液體水轉(zhuǎn)移特性和發(fā)展數(shù)值評估方法。該研究吸引了許多著名學(xué)校和重點(diǎn)實(shí)驗(yàn)室參與進(jìn)來，其結(jié)果可靠，并且提供了詳細(xì)的設(shè)置說明，方便模擬驗(yàn)證。

為了準(zhǔn)確地反映極端熱濕地區(qū)的特性，采用HAMSTAD-benchmark 4模型進(jìn)行模擬驗(yàn)證。HAMSTAD-benchmark 4模型是一個(gè)比較特殊的模型，該模型分析了暴露在非穩(wěn)態(tài)室外氣候和相對穩(wěn)定的室內(nèi)氣候條件下外墻內(nèi)的溫度和含濕量分布，在研究熱濕耦合傳遞的基礎(chǔ)上，進(jìn)一步探索了雨水和太陽輻射對圍護(hù)結(jié)構(gòu)的影響，同時(shí)，整個(gè)模型是氣密性的，蒸汽將以滲透的形式傳輸。模擬時(shí)長為120 h。此外，其物性參數(shù)隨含濕量的變化如表1所示，滿足研究極端熱濕地區(qū)的需要。

3.2 COMSOL軟件設(shè)置的研究

為了保證模擬的準(zhǔn)確性，先對COMSOL軟件模擬的設(shè)置對計(jì)算結(jié)果影響進(jìn)行了分析，首先進(jìn)行容差設(shè)定分析，COMSOL提供了相對容差（Rtol）和絕對容差（Atol）對誤差進(jìn)行限定。因?yàn)榻^對容差幾乎總是必須被修改，因此，沒法給出一個(gè)最好的通用缺省值，使用一個(gè)很小的絕對容差可能是一個(gè)很好的選擇，這會使相對容差成為控制誤差的決定因素，筆者在絕對容差設(shè)定時(shí)都采用比預(yù)期的變化小一個(gè)數(shù)量級。相對容差軟件默認(rèn)的是0.01，將相對容差分別設(shè)定成10-9、10-3和10-2，對結(jié)果進(jìn)行最大誤差分析。

相對容差10-3和10-9之間含濕量的最大相對誤差為0.7%，10-2和10-9之間含濕量的最大相對誤差為12.8%。從結(jié)果上看，相對容差設(shè)定成10-9或10-3較為合適，按默認(rèn)取值會造成傳濕計(jì)算的較大誤差，然而，相對容差設(shè)定成10-9的計(jì)算時(shí)間是其他兩個(gè)值的10倍左右，為了方便，相對容差推薦設(shè)定成比默認(rèn)值小10倍，即0.001。

此外，對網(wǎng)格劃分對結(jié)果的影響進(jìn)行對比分析，網(wǎng)格按極細(xì)化、較細(xì)化和常規(guī)進(jìn)行劃分，最大相對誤差都在0.1%以內(nèi)。結(jié)果上看，COMSOL軟件默認(rèn)的網(wǎng)格劃分都較細(xì)致，但由于計(jì)算時(shí)間差距不大，推薦根據(jù)計(jì)算的目標(biāo)來采用極細(xì)化或較細(xì)化的網(wǎng)格來減少網(wǎng)格劃分帶來的誤差。

此外，值得注意的是，當(dāng)材料處于濕飽和狀態(tài)，毛細(xì)壓力超過飽和毛細(xì)壓力時(shí)，雨水將不再進(jìn)入圍護(hù)結(jié)構(gòu)而是直接落下。因此，在計(jì)算的時(shí)候，需要將其分別計(jì)算。筆者通過在全局變量添加插值函數(shù)satu（w）來模擬這一過程。

因此，淋雨的濕通量為grain×satu（w）。

3.3 結(jié)果對比

圖1～圖4中TEST曲線為該模型的計(jì)算結(jié)果，其他曲線是由HAMSTAD參與學(xué)?；蚪M織提供，分別為：Chalmers University of Technology， Coordinator of WP（CTH）、TNO Building and Construction Research （TNO）、University of Leuven （KUL）、Technical University of Dresden （TUD）、Technion Israel Institute of Technology （Technion）、Eindhoven University of Technology （TUE）、National Research Council of Canada （NRC）、Fraunhofer Institute of Building Physics （IBP）。

從結(jié)果來看，無論內(nèi)表面還是外表面，該模型所代表的曲線TEST均可以準(zhǔn)確地反映出在雨水和太陽輻射影響下多層圍護(hù)結(jié)構(gòu)的熱濕耦合傳遞過程，模擬結(jié)果準(zhǔn)確地反映了太陽輻射影響下圍護(hù)結(jié)構(gòu)表面的波動，以及在雨水作用下圍護(hù)結(jié)構(gòu)的傳濕。模擬結(jié)果與基準(zhǔn)值基本一致，該模型模擬結(jié)果良好。

4 結(jié)論

以毛細(xì)壓力與溫度為驅(qū)動勢，建立了一種適應(yīng)于極端熱濕地區(qū)圍護(hù)結(jié)構(gòu)熱濕耦合傳遞模型，在該模型的基礎(chǔ)上，提出了基于計(jì)算機(jī)軟件COMSOL和MATLAB的熱濕耦合模型的求解方法，并與HAMSTAD基準(zhǔn)進(jìn)行了驗(yàn)證，驗(yàn)證結(jié)果良好。得到以下結(jié)論：

1）COMSOL模擬熱濕耦合傳遞時(shí)，相對容差推薦設(shè)定為0.001，絕對容差設(shè)定為0.000 1，網(wǎng)格按極細(xì)化劃分。

2）在極端熱濕地區(qū)，為了準(zhǔn)確反映高溫和多雨對圍護(hù)結(jié)構(gòu)的影響，提出了適宜于極端熱濕地區(qū)的熱濕耦合模型邊界條件，該方法可以準(zhǔn)確地模擬出太陽輻射和雨水對圍護(hù)結(jié)構(gòu)傳熱傳濕的影響。

3）在計(jì)算雨水負(fù)荷時(shí)，需要注意圍護(hù)結(jié)構(gòu)的吸附能力，當(dāng)圍護(hù)結(jié)構(gòu)表面達(dá)到飽和時(shí)，多余的雨水應(yīng)該不參與計(jì)算。

參考文獻(xiàn)：

[1] 劉毅. 提升南海氣象觀測能力刻不容緩——專訪中國工程院院士丁一匯[N]. 人民日報(bào)， 2015-06-21（2）.

LIU Y. Enhance the ability of the South China Sea meteorological observations without delay [N]. People's Daily， 2015-06-21（2）. （in Chinese）

[2] HOU X D， LIU F S， WANG S J. Coupled heat and moisture transfer in hollow concrete block wall filled with compressed straw bricks [J]. Energy and Buildings， 2017， 135（15）： 74-84.

[3] KONG F H， ZHANG Q L. Effect of heat and mass coupled transfer combined with freezing process on building exterior envelope [J]. Energy and Buildings， 2013， 62： 486-495.

[4] FANG X D， ANDREAS K A， PAUL P F. Methodologies for shortening test period of coupled heat-moisture transfer in building envelopes [J]. Applied Thermal Engineering， 2009， 29（4）： 787-792.

[5] TOM K， JAROSLAV K， MICHAL ， et al. Numerical analysis of coupled heat and moisture transport in masonry [J]. Computers & Mathematics with Applications， 2017， 74（1）： 229-248.

[6] MIN H G， ZHANG W P， ROBERT C. Coupled heat and moisture transport in damaged concrete under an atmospheric environment [J]. Construction and Building Materials， 2017， 143： 607-620.

[7] KERGET L， TADEU A， RAVNIK J. BEM numerical simulation of coupled heat， air and moisture flow through a multilayered porous solid [J]. Engineering Analysis with Boundary Elements， 2017， 74： 24-33.

[8] 劉向偉，陳友明，陳國杰， 等. 圍護(hù)結(jié)構(gòu)熱濕耦合傳遞模型及簡便求解方法[J]. 土木建筑與環(huán)境工程，2016， 38（4）： 7-12.

LIU X W， CHEN Y M， CHEN G J， et al.Coupled heat and moisture transfer model and simple solution method for building envelopes [J]. Journal of Civil， Architectural & Environmental Engineering， 2016， 38（4）： 7-12. （in Chinese）

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

ZHANG J X， GUO X G， CHEN Y M， et al.A time domain recursive algorithm for solving the model of coupled heat and moisture transfer in building wall [J]. Journal of Civil， Architectural & Environmental Engineering， 2015， 37（6）： 147-152. （in Chinese）

[10] 王瑩瑩. 圍護(hù)結(jié)構(gòu)濕遷移對室內(nèi)熱環(huán)境及空調(diào)負(fù)荷影響關(guān)系研究[D].西安：西安建筑科技大學(xué)， 2013.

WANG Y Y. Research on the effect of the palisade structure moisture transfer on the indoor thermal environment and air-conditioning load [D]. Xi'an： Xi'an University of Architecture and Technology， 2013. （in Chinese）

[11] 彥啟森. 建筑熱過程[M].北京：中國建筑工業(yè)出版社， 1986.

YAN Q S.Building heat process [M]. Beijing： China Architecture & Building Press， 1986. （in Chinese）

[12] A.B.雷柯夫. 建筑熱物理理論基礎(chǔ)[M]. 北京： 科學(xué)出版社， 1965.

LUIKOV A B. Theoretical basis of building thermophysics [M]. Beijing： Science Press， 1965. （in Chinese）

[13] LI Q R， JIWU R， FAZIO P. Development of HAM tool for building envelope analysis [J]. Building and Environment， 2009， 44（5）： 1065-1073.

[14] Smithsonian Meteorological Tables [Z]. Washington： Smithsonian Institution，1984.

[15] 周西華，梁茵，王小毛. 飽和水蒸汽分壓力經(jīng)驗(yàn)公式的比較[J]. 遼寧工程技術(shù)大學(xué)學(xué)報(bào)，2007（3）：331-333.

ZHOU X H， LIANG Y， WANG X M. Comparison of saturation vapor pressure formulas [J]. Journal of Liaoning Technical University， 2007（3）： 331-333.（in Cinese）

[16] 中華人民共和國建設(shè)部. 民用建筑熱工設(shè)計(jì)規(guī)范： GB 50176—93 [S]. 1993.

Ministry of Construction of the People's Republic of China. Thermal design code for civil building： GB 50176-93 [S]. 1993.（in Chinese）

[17] KNZEL H M. Simultaneous heat and moisture transport in building components： one and two-dimensional calculation using simple parameters [D]. Institute of Building Physics， Germany，1995.

[18] HAGENTOFT C E. HAMSTAD-Final report： methodology of HAM-modeling，Report R-02：8 [R]. Gothenburg， Department of Building Physics， Chalmers University of Technology， 2002： 8-50.

（編輯 胡英奎）