墻體內(nèi)熱濕耦合過程的時域遞歸展開算法

摘要：利用時域遞歸展開算法對墻體內(nèi)熱濕耦合傳遞方程進行求解。以木板為例，應用該算法進行了熱濕耦合傳遞的分析計算，在時間域和空間域上分別運用遞歸展開法和控制容積法進行離散，從而得到遞歸形式的線性方程組，運用MATLAB軟件對這一過程進行求解。計算結果與有限差分算法、解析解計算結果以及實驗數(shù)據(jù)吻合良好，表明該算法能夠用于求解多孔介質(zhì)熱濕耦合傳遞模型。時間步長的改變對計算結果影響較小，可通過增加時間步長方法來減少工作量。

關鍵詞：熱濕耦合傳遞；時域分析；時間步長

中圖分類號：TU111.19 文獻標志碼：A 文章編號：1674-4764（2015）06-0147-06

Abstract：A time domain recursive algorithm was used to solve the model of coupled heat and moisture transfer equation through building wall. The calculation was carried out by using algorithm in wood.The recursive algorithm and the finite volume algorithm were separately used in the time and space domains respectively to obtain the linear equations of recursive form. The results were compared with that of finite difference algorithm and showed that the result of algorithm agreed with that of the finite difference algorithm， analytical solution and experimental data. The algorithm can be used to solve coupled heat and moisture transfer model in porous media. The change of the time step has slight effects on results and increasing time step would reduce the calculating time .

Key words：coupled heat and moisture transfer； time domain analysis； the sizes of time step

在熱濕氣候地區(qū)，墻體的熱濕遷移對圍護結構的熱工性能、建筑能耗和室內(nèi)環(huán)境有著十分重要的影響[1-2]。因此，研究多孔材料的熱濕遷移，預測材料內(nèi)部溫度和濕度的分布對改善圍護結構的熱工性能、抑制霉菌生長、提高圍護結構的耐久性具有重大意義。而建立多孔材料熱濕傳遞模型和采用高精度的算法是研究這些問題的基本前提。

許多學者采用的熱濕耦合傳遞模型為Luikov模型[3]。該模型可以就建筑材料和建筑結構在真實外界條件下長期變化的特性做出模擬，奠定了液態(tài)水和水蒸氣同時傳遞模型的理論基礎。一些學者對Luikov方程進行了改進和完善[4-5]，得到了一維瞬態(tài)熱濕耦合傳遞方程的解析解。對于熱濕耦合傳遞模型的數(shù)值解法，被廣泛采用的主要是有限元法和有限差分法[6-8]。Kerget等[9]采用基于邊界元法的數(shù)值模型對多孔介質(zhì)非線性熱濕耦合傳遞進行了模擬，此方法減少了問題的維數(shù)，但應用范圍以存在相應微分算子的基本解為前提。郭興國等[10]運用有限容積法對多層墻體熱濕耦合模型進行了求解，此類方法時間步長選擇受到一定的限制。Mendes等[11]運用有限差分方法對方程進行離散，并用MTDMA（Multitridiagonal-Matrix Algorithm）對離散后的方程進行求解。李鳳志等[12]利用時域遞歸展開算法求解了織物熱濕耦合方程，該方法的主要優(yōu)點是時間步長的選取幾乎不受限制。

筆者建立了多孔介質(zhì)一維熱濕傳遞控制方程，在時域上采用了遞歸展開算法，有效克服了有限差分法對時間步長的限制問題，在空間域上采用有限容積法對方程進行離散，并且用MATLAB編程進行求解。計算結果與有限差分方法的計算結果以及文獻[5]中的解析解、實驗數(shù)據(jù)進行比較，證明了此解法的正確性。研究結果表明，該方法能夠用于墻體熱濕耦合遷移方程的求解。

1 熱濕耦合模型的建立

1.1 控制微分方程

在計算過程中，也可預先設置展開階數(shù)的上限為max，改變max的大小，可以控制求解精度以及程序運行的時間。

3 分析與討論

將上述求解方法應用到求解建筑材料熱濕耦合過程中，模型計算采用的材料為木板，邊界條件及初始條件按照文獻[5]中設置，模型其他計算參數(shù)主要根據(jù)文獻[13]選取，如表1所示。計算結果與有限差分方法、文獻[5]中解析解進行了對比，經(jīng)過有限差分離散后的方程分別用直接解法和MTDMA算法進行求解。

圖3為木質(zhì)建筑構件的中點溫度隨時間變化的曲線，對比可以看出，時域遞歸算法預測材料溫度的變化與文獻[5]中的解析解計算的結果具有高度一致性，平均相對誤差為0.1%。有限差分算法中的直接解法及MTDMA算法計算的結果是一致的，與時域遞歸算法相比，平均相對誤差均為1.4%。在大約1 500 s后，溫度均穩(wěn)定在83.5 ℃。圖4為建筑構件中點濕度勢隨時間變化曲線，對比可知：采用時域遞歸算法和有限差分算法得到的曲線幾乎重合，平均相對誤差僅為0.074%，經(jīng)過計算，在大約48 h后，濕度勢穩(wěn)定在了24.5 °M。與文獻[5]中的結果對比，平均相對誤差為2.4%，數(shù)值較與有限差分算法對比時的誤差增加。分析原因，可能是由于解析解對模型做了一定的假設和簡化導致的[5]。

如圖5、6所示為該方法與Qin等[16]實驗結果的對比，計算條件與實驗條件相同，實驗天數(shù)為5周，材料屬性見文獻[16]。圖5給出了墻體溫度分布，計算結果與實驗結果吻合較好，相對誤差小于10%。圖6為墻體水蒸氣含量分布，從圖中可以看出計算結果與實驗結果吻合良好，相對誤差小于7%。

通過與數(shù)值解、解析解以及實驗對比，該算法誤差較小且誤差均在可以接受的范圍內(nèi)，證明了該算法的正確性，可以用來求解墻體內(nèi)熱濕耦合傳遞模型。

為了研究采用該算法時間步長的選取對計算精度的影響，將時域遞歸展開算法與有限差分的隱式格式進行了對比，結果如圖7、8所示。計算時室外溫度35 ℃，相對濕度80%，室內(nèi)溫度25 ℃，相對濕度60%，材料的厚度為0.1 m，其余計算參數(shù)見表1，總時間設定為240 h。材料的含濕量與相對濕度的函數(shù)關系參考文獻[17]。

圖7（a）為采用時域展開算法預測材料中點溫度隨時間變化曲線。從圖中可以看出，隨著運行時間的增加，溫度變化較快。在2×104 s時，溫度穩(wěn)定在30.06 ℃，所以在圖中僅顯示20 h內(nèi)溫度隨時間變化曲線。對比圖5（a）和（b）可以看出，有限差分格式在取不同的時間步長時，曲線波動較大，而圖5（a）的時域遞歸展開算法對結果幾乎沒有影響。

圖8（a）為采用時域展開算法預測的相對濕度隨時間變化曲線。從圖中可以看出隨著運行時間的增加，相對濕度逐漸增加，并且有趨于穩(wěn)定的趨勢。延長計算的總時間，相對濕度穩(wěn)定在了70%。對比圖8中的（a）和（b），可以得出：用時域遞歸展開算法計算得到的相對濕度幾乎不受時間步長的影響，圖線基本重合。而有限差分算法對時間步長的選擇具有一定的依賴性。綜合圖3～8可以得出：溫度比含濕量能夠較快進入穩(wěn)定的狀態(tài)，材料的熱擴散系數(shù)大于濕擴散系數(shù)。

4 結 論

利用時域遞歸展開算法對熱濕耦合傳遞模型進行了求解，在時域上采用了遞歸展開算法進行離散，在空間域上采用控制容積法，用MATLAB編程計算。通過與有限差分算法進行對比、分析可以得到：

1）材料中點溫度和濕度勢隨時間變化曲線的大致走勢是一致的，驗證了算法的正確性。

2）增加展開階數(shù)可提高計算結果的精度，減少展開階數(shù)可減少工作量。

3）時域遞歸展開算法的計算精度幾乎不受時間步長的影響。因此，可以適當?shù)脑黾訒r間步長來減少工作量，節(jié)約計算所需時間。

4）可以將自適應的精細算法與有限差分方法結合起來，解決計算精度要求高的物理問題。

參考文獻：

[1]孔凡紅， 王懷柱. 圍護結構熱質(zhì)耦合傳遞對室內(nèi)環(huán)境的影響[J]. 熱科學與技術， 2010， 9（3）： 212-217.

Kong F H， Wang H Z. Effect of envelope heat and mass coupled transfer on indoor environment [J]. Journal of Thermal Science and Technology， 2010， 9（3）： 212-217. （in Chinese）

[2]孔凡紅， 鄭茂余， 韓宗偉， 等. 新建建筑圍護結構熱質(zhì)傳遞對室內(nèi)溫濕度環(huán)境的影響[J]. 建筑科學， 2008，24（8）： 94-97.

Kong F H， Zheng M Y， Han Z W， et al. Influences of envelope heat and moisture transfer of new building on indoor temperature and humidity [J]. Building Science， 2008， 24（8）： 94-97. （in Chinese）

[3 ]Luikov A W. Heat and mass transfer in Capillary-Porous bodies [M]. Pergamon Press， Oxford， 1966.

[4 ]Qin M H， Belarbi R， At-Mokhtar A， et al. An analytical method to calculate the coupled heat and moisture transfer in building materials [J]. International Communications in Heat and Mass Transfer， 2006， 33（1）：39-48.

[5 ]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）： 1354-1360.

[6 ]Qin M H， Belarbi R， At-Mokhtar A. Coupled heat and moisture transfer in multi-layer building materials [J]. Construction and Building Materials， 2009， 23： 967-975.

[7 ]Qin M H， At-Mokhtar A， Belarbi R. Two-dimensional hygrothermal transfer in porous building materials [J]. Applied Thermal Engineering， 2010， 30： 2555-2562.

[8 ]郭興國. 熱濕氣候地區(qū)多層墻體熱濕耦合遷移特性研究[D]. 湖南：湖南大學， 2010： 32-34.

Guo X G. Research on coupled heat and moisture transfer characteristics of multilayer walls in hot and humid climate [D]. Changsha： Hunan University， 2010： 32-34. （in Chinese）

[9 ]Skerget L， Tadeu A. BEM numerical simulation of coupled heat and moisture flow through a porous solid [J]. Engineering Analysis with Boundary Elements， 2014， 40： 154-161.

[10]郭興國， 陳友明. 熱濕氣候地區(qū)墻體內(nèi)部冷凝分析[J]. 土木建筑與環(huán)境工程， 2011， 33（5）： 121-124.

Guo X G， Chen Y M. Aanlysis of condensation in porous wall exposed to hot humid climate [J]. Journal of Civil， Architectural and Environmental Engineering， 2011， 33（5）： 121-124. （in Chinese）

[11]Mendes N， Philippi P C. Multitridiagonal-matrix algorithm for coupled heat transfer in porous media： stability analysis and computational performance [J]. Journal of Porous Media， 2004， 7（3）： 193-211.

[12]李鳳志， 李毅， 曹業(yè)玲. 求織物熱濕耦合方程的控制體-時域遞歸展開算法[J]. 南京航空航天大學學報， 2009， 41（3）： 319-323.

Li F Z， Li Y， Cao Y L. Control volume-time domain recursive method for solving coupled equations of heat and moisture transfer in fabrics [J]. Journal of Nanjing University of Aeronautics and Astronautics， 2009， 41（3）： 319-323. （in Chinese）

[13]Chang W J， Weng C I. An analytical solution to coupled heat and moisture diffusion transfer in porous materials [J]. International Journal of Heat and Mass Transfer， 2000， 43： 3621-3632.

[14]Qin M H， Belarbi R. Development of an analytical method for simultaneous heat and moisture transfer in building materials utilising transfer function method [J]. Journal of Materials in Civil Engineering， 2005， 17（5）： 492-497.

[15]Yang H T， Gao Q. A Precise time stepping scheme to solve hyperbolic and parabolic heat transfer problems with radiative boundary condition [J]. Heat and Mass Transfer， 2003， 39： 571-577.

[16]Qin M H， Belarbi R， At-Mokhtar A， et al. Nonisothermal moisture transport in hygrocopic building materials： modeling for the determination of moisture transport coefficients [J]. Transport Porous Media， 2008， 72（2）： 255-271.

[17]Hameury S. Moisture buffering capacity of heavy timber structures directly exposed to an indoor climate： a numerical study [J]. Building and Environment， 2005， 40： 1400-1412.

（編輯 胡英奎）