于水,姚鈺鋒,李銀宗

(沈陽建筑大學市政與環(huán)境工程學院,沈陽 110168)

世界正在成為一個城市環(huán)境[1],包括中國在內的全球新興地區(qū)目前正在以人類歷史上前所未有的速度和規(guī)模經歷城市化[2]。在過去二十年中,中國經歷了一個快速城市化的時期,城市人口從1978年的18%上升到1998年的29%,2016年達到57%[3]。不斷出現(xiàn)的土地資源限制、環(huán)境壓力、交通擁堵和資源分配不公平等問題嚴重威脅著城市生活環(huán)境。然而,地下空間將能夠提供25%～40%的額外空間,而不會占用寶貴的地表空間資源[4]。城市地下空間(UUS)開發(fā)已成為解決城市發(fā)展過程中出現(xiàn)的交通、環(huán)境和土地使用問題的一個有吸引力的解決方案[5]。

面對地下空間的不斷發(fā)展,對地下建筑的能效[6]、溫濕度和室內氣候控制的要求將會愈發(fā)嚴格。然而,由于降水和融雪產生的水分、土壤/回填的高相對濕度(RH)以及地下水的存在,地下建筑受到土壤比地面更大的水分應力。土壤和防水層的存在限制了建筑物向外干燥的能力。因此,相較于地面上的建筑,地下建筑的圍護結構不得不面臨更加潮濕的問題。排水不充分、設計不當?shù)膰o結構或新建結構會逐漸累積大量水分[7],這可能導致機械和物理性能的變化,降低承載能力和整體熱性能[8]。為了準確捕捉建筑圍護結構對室內環(huán)境的影響,必須建立瞬態(tài)熱、空氣和水分(HAM)傳遞模型,研究通過建筑圍護結構的耦合熱、空氣與水分傳遞現(xiàn)象。

現(xiàn)階段中外對于地下建筑的研究多集中于室內熱濕環(huán)境控制、圍護結構傳熱分析、通風等方面,且相關熱濕耦合傳遞研究也多集中于地上建筑。Vajdi等[9]綜述了地下建筑的建筑能源性能、熱舒適性和室內空氣質量性能的相關研究。華永等[10]對地下全覆蓋式和半覆蓋式的辦公建筑進行熱濕環(huán)境測試,并使用ANSYS模擬分析優(yōu)化后的全覆蓋房間內溫濕度及APMW(預計適應性平均熱感覺指標)。劉妙坤等[11]通過建立三維非穩(wěn)態(tài)傳熱模型分析地下建筑熱流及溫度場變化并與3種常見方法進行了對比。Zoras[12]綜述了與土壤耦合的建筑結構的傳熱預測方法,并將其分為四類:分析/半分析方法、數(shù)值方法、手動方法和設計指南。

與地下建筑的室內邊界條件不同,確定位于土壤側的外邊界條件顯得尤為困難。對于位于較深位置的地下建筑,土壤溫度可以認為是恒定的[13]。然而,由于受到諸多條件的影響,淺層地下建筑的外部邊界條件的確定變得更具挑戰(zhàn)性:①外部氣候條件,包括太陽輻射、遮蔭、降水、積雪等;②地下水位高度的影響;③土壤的濕熱性質、成分組成、水分含量和凍融,以及④由于高濕度對濕熱特性產生的影響。Asphaug等[14]通過對十篇地下建筑熱濕傳遞研究的仔細審查后,認為使用全尺寸物理測量對地下室的濕熱模擬缺乏徹底的驗證。

為優(yōu)化相關理論數(shù)值模型以提高地下建筑的濕熱模擬準確性并填補地下建筑熱濕研究相關研究空缺,以地下建筑為研究對象,在大量梳理中外地下建筑熱濕傳遞研究的基礎上,綜述了多孔建筑材料熱濕耦合傳遞機理和現(xiàn)階段的不足,回顧研究中的熱濕傳遞數(shù)值模型和計算工具及方法并對比分析了邊界條件設置方法,總結存在的問題并作出展望,為后續(xù)相關研究發(fā)展提供參考。

1 熱濕耦合傳遞機理

1.1 水分傳遞

在多孔材料中,水分轉移包括液體和氣體的水分轉移[15]。吸濕材料中的水分傳輸是由幾個作用于相(液體和蒸汽)運動的力學驅動的。根據(jù)不同的驅動力,液體和蒸汽被分開處理。它們共同得出了水分的傳輸方程?？紤]到水汽條件下的水轉移,可以根據(jù)孔徑和分子密度設想兩種情況:一是分子擴散可用于描述半徑大于平均自由程的孔隙中的蒸汽轉移;二是當孔隙如此之小,以至于分子和孔壁之間的碰撞比分子之間的碰撞更頻繁時,稱為滲出或努森傳輸。動力學氣體理論通過基本上包含3種擴散勢(質量分數(shù)、溫度和總壓力)的方程描述多組分氣體混合物中分子的擴散[16]。由質量分數(shù)差異引起的擴散通常稱為菲克擴散。假設總壓力梯度和基于溫度梯度的熱擴散(索爾特效應)可以忽略不計。由于理想氣體中組分的質量分數(shù)與其在總壓力中的份額[17]和引入水蒸氣擴散阻力系數(shù)之間存在比例關系,因此可以使用水蒸氣分壓計算多孔介質中的菲克擴散,可表示為

(1)

在等溫條件下,該蒸汽壓力梯度很容易轉化為蒸汽濃度梯度,符合Fick的原始定律[18],該定律也是從等溫實驗中推導出來的[19]。水蒸汽的分壓被濃度梯度代替,可表示為

(2)

式(2)中:gvap為矢量蒸汽通量密度,g/(cm2/s);v為質量流量系數(shù);α為允許額外路徑長度的曲折系數(shù);a為材料的體積空氣含量,cm3/cm3;ρ為水蒸氣密度,g/cm3。

液相水分轉移[20]通過吸附或解吸過程中的毛細吸力或重力進行操作。在早期的菲利普方程中,液體通量分為3個分量,即溫度梯度、濕度梯度和重力梯度。假設由于熱梯度和重力,液體中的水分轉移可以忽略不計,達西定律可以描述液體通量密度。

(3)

式(3)中:gw為含水量的流動密度,kg/(m2/s);Dw為毛細管輸送系數(shù),m2/s;w為含水量,kg/m3。

1.2 傳熱

根據(jù)能量守恒定律,受控元件中的焓變化等于進入受控元件的能量流(由液態(tài)水傳輸和蒸汽傳輸引起的熱傳導和焓流[21]。傳熱方程可表示為

(4)

式(4)中:ρd為干燥材料的密度,kg/m3;c為潮濕材料的比熱,J/(kg·K);λ為材料的導熱系數(shù),W/(m·K);Lv為蒸發(fā)潛熱,J/kg;gv為水蒸氣擴散通量密度,kg/(m2/s)。

在多孔介質的一般傳熱方程中,通常忽略顯熱部分,只考慮潛熱,也選擇忽略液體擴散。

2 熱濕傳遞數(shù)值方程

質量守恒和能量守恒原理被用于建立熱量和水分傳遞模型。水分運動的驅動力是多孔建筑圍護結構和材料的不同濕熱模型之間的共同因素。水分含量、毛細壓力、水蒸氣分壓、相對濕度和空氣含水量是當前濕熱模型中使用的一些水分驅動力。表1[22-33]中給出了各類熱濕耦合傳遞數(shù)值方程應用于地下建筑的相關研究。

表1 地下建筑熱濕模擬研究[22-33]Table 1 Hygrothermal simulation studies of underground buildings[22-33]

2.1 Philip and De Vries模型

1957年,Philip等[34]提出描述多孔材料中溫度梯度和濕度梯度聯(lián)合作用下的濕度和傳熱方程。方程的建立是基于毛細作用的蒸汽擴散和液體運動的經典機制上發(fā)展而來的。水分和熱量傳輸方程的主要驅動力分別為體積含水量和溫度,描述溫度梯度和濕度梯度聯(lián)合作用下多孔介質中水分運動的一般微分方程如式(5)所示。

(5)

Dθ=Dθ,liq+Dθ,vap

(6)

DT=DT,liq+DT,vap

(7)

式中:θ為體積含水量,m3/m3;t為時間,s;Dθ為等溫水分擴散系數(shù),m2/s;DT為溫度水分擴散系數(shù),m2/(s·K);T為溫度,K;Dθ,liq、Dθ,vap分別為等溫水分擴散系數(shù)和水蒸氣擴散系數(shù)(與水分梯度有關),m2/s;DT,liq、DT,vap分別為溫度水分擴散系數(shù)和水蒸氣擴散系數(shù)(與溫度梯度有關),m2/(s·K);Kh為導水率(一般與重力有關),m/s;z為高度,m。

傳熱方程可表示為

(8)

式(8)中:C為體積熱容,J/(kg·K);等號右側第二項表示由水分梯度引起的蒸餾效應。

2.2 Künzel模型

1995年,Künzel[35]考慮到非物理驅動力或多層建筑構件中不連續(xù)的電位被認為不適合建模任務,以相對濕度和溫度作為勢推導出建筑構件中的熱濕耦合傳遞模型。不考慮總壓差或重力引起的蒸汽和液體對流,以及液體流動引起的焓變,它們對熱平衡只有不可忽略的影響。

(9)

(10)

2.3 Zhang模型

2004年,張華玲[36]參考Pederson[37]的研究,基于建筑多孔圍護結構任意一點的濕遷移可以描述為Jm=Jv+Jl,其中Jm、Jv和Jl分別為水分、水蒸氣和液態(tài)水擴散通量,kg/(m2/s),提出了相對于孔隙飽和水分含量的相對濕度(φ)和溫度(T)為驅動勢多孔材料熱濕傳遞過程的理論模型。模型同時考慮了水蒸氣和液態(tài)水的擴散遷移過程。濕傳遞方程為

(11)

(12)

(13)

熱傳遞方程為

(14)

(15)

式(15)中:L(T)為水蒸氣潛熱,J/kg;λm為多孔材料導熱系數(shù),W/(m·K);λapp為材料表觀導熱系數(shù),W/(m·K)。

2.4 Matsumoto模型

1978年,Matsumoto[38]給出了同時建立熱濕傳遞方程是整體分析的控制方程,用于描述混凝土和地面中熱量和水分的同時傳遞。這是在沒有固相(冰)的情況下,假定水的液相和氣相的局部平衡而導出的方程為

(16)

(17)

式中:μ為相對于自由水的水化學勢,J/kg;ψ為水分含量,m3/m3;λ′μ水化學勢梯度的總導濕系數(shù),kg/m·s·(J/kg);λ′T為溫度梯度下的導濕系數(shù),kg/m·s·K;λ′μg為水化學勢梯度下氣相水分傳導系數(shù),kg/m·s·(J/kg);λ′Tg為溫度梯度下的氣相導濕系數(shù),kg/m·s·K;ρw、ρ分別為液態(tài)水和土壤密度,kg/m3;c為含有水分的土壤的熱容,J/(kg·K);r為汽化潛熱,J/kg。

2.5 Grunewald模型

Delphin的歷史可以追溯到1987年,根據(jù)文獻[39]中的方程所開發(fā)的瞬態(tài)一維傳熱代碼。雖然Delphin 5/CHAMPS-BES中使用的物理模型與Delphin4中的基本相同,但平衡方程系統(tǒng)的公式略有變化。引入揮發(fā)性有機化合物的質量平衡方程以及鹽運移和結晶過程的平衡方程[40],但不討論這些熱濕傳遞之外的方程。根據(jù)Grunewald[41]的研究,濕分和能量守恒方程為

(18)

(19)

2.6 Milly模型

1982年,Milly[42]基于Philip等[34]提出的土壤中熱和濕耦合的現(xiàn)象方程,使用變量ψ和T對方程進行了一些擴展和調整。所提出的方程組適用于非均質土和滯回土。包括植物根系吸收水分的可能性。將質量平衡方程展開到自變量矩陣頭ψ(m)和溫度T(K)得

(20)

式(20)中:ρv、ρl分別為液態(tài)水和蒸汽密度,kg/m3;Dψv為隱含定義的傳遞系數(shù),m/s;DTv為水力蒸汽擴散系數(shù),m2/s·K;θa為體積空氣含量,m3/m3;θ為體積含水量,m3/m3;K為土壤基質的水力傳導率,m/s。

土壤中熱傳遞的最終表達式為

(21)

式(21)中:C=Cd+clρlθ+(caρa+cvρv)θa為土壤的體積熱容量,J/(kg·K);Cd為干燥土壤材料的體積熱容,J/(kg·K);cl、ca、cv分別為液態(tài)水、干燥空氣和水蒸氣的比熱;T0為任意參考溫度,℃。

在目前的研究中,已經使用了各類驅動勢來計算水分傳輸。在上述方程中,Künzel模型和Zhang模型使用相對濕度作為計算水分傳輸?shù)尿寗恿?。然?相對濕度的缺點是,當其高于100%時,結果不準確且不可靠。Dong等[43]通過Fortran和COMSOL對Liu模型和Künzel模型的充分驗證認為Künzel模型不適用于過吸濕區(qū)域。因此,在模擬土壤等具有高濕度的多孔介質時,使用相對濕度作為水分驅動勢并不是十分合適。在Philip and De Vries模型、Matsumoto模型和Milly模型中以土壤作為參考多孔介質,分別采用了含水量、水化學勢和基質水頭作為驅動勢。然而,非物理驅動力或多層建筑構件中不連續(xù)的驅動勢被認為不適用于建模任務,驅動勢應是易于測量且較為熟悉的量。Janssen[44]對比分析4種常用于多孔建筑材料的3種水分驅動勢:對數(shù)毛細壓力(Lpc)毛細壓力(pc)和相對濕度。他認為這3種驅動勢中的任何一種都沒有一致的優(yōu)勢,且在不同飽和條件下顯示出不同的性能強弱。因此,很明顯,在對地下建筑圍護結構(包括土壤區(qū)域)進行濕熱模擬時,難以使用單個水分驅動勢來描述水分傳輸,這將不可避免地導致計算錯誤。為了減少計算成本同時避免土壤域的高濕度問題,兩步法廣泛應用于地下建筑的熱濕研究中。首先,通過對地下圍護結構和周圍土壤的熱模擬,確定土壤剖面的溫度分布。其次,對圍護結構進行濕熱模擬,將地下室墻外表面的溫度隨時間變化作為相對濕度接近100%時的外邊界溫度條件。Janssen等[45]在研究地下室通過土壤的熱損失時,選擇忽略圍護結構,嘗試將封套視為熱阻,并將地下圍護結構看作對土壤域是不透水的。相比之下,Ma等[25]在研究地面以下的圍護結構時,與土壤域接觸的邊界被認為是溫度連續(xù)但圍護結構外部不透水。兩步法雖然能夠適當估計地下建筑的熱損失,但可能并不足以用于濕熱性能的研究。簡單地將建筑的外邊界條件設為100%,顯然忽略了降水、積雪以及地下水等外界因素對于土壤水分分布的影響。

3 熱濕耦合模擬工具

在建筑物理領域,濕熱模擬被廣泛應用于預測及建筑材料、構件和整個建筑的濕熱性能[46]。2005年,加拿大抵押和住房公司(CMHC)研發(fā)了45種計算機化濕熱建模工具[47]。據(jù)Delgado等[48]的研究,到2013年濕熱數(shù)值工具至少有57種。然而,由于它們之間存在一些根本性的差異,如不同的數(shù)學模型、不同的驅動勢[49]、材料屬性的數(shù)量和實現(xiàn)方式,以及邊界條件(氣候荷載)并使用不同的數(shù)值求解方式,會根據(jù)邊界條件、模擬的意圖或需要解決的問題,可能導致不同的結果[50]。因此,選擇正確的模擬工具對于地下建筑熱濕性能預測和評估十分重要。表2中介紹了3個現(xiàn)階段相關研究中較為常用的熱濕模擬軟件。

表2 常見濕熱模擬工具總結Table 2 Summary of common used hygrothermal simulation tools

3.1 節(jié)點法

節(jié)點法基于這樣一個原則,即每個建筑分區(qū)都是一個單獨的體積,定義為一組處于統(tǒng)一狀態(tài)的物理變量。節(jié)點通常表示房間、地板或大廳等,但它也可能反映更復雜的內容,如內部負荷,如占用、設備或暖通空調系統(tǒng)散發(fā)的熱量。因此,一個區(qū)域基本上相當于一個具有自身物理參數(shù)的節(jié)點,如溫度、濕度、壓力等。對于整個建筑體積的每個節(jié)點,可以使用矩陣計算方法求解熱濕平衡方程。

節(jié)點法是一種一維的計算方法,常被用于建筑物能源仿真軟件,如TRNSYS,EnergyPlus以及基于特定熱濕耦合傳遞方程發(fā)展而來的用于準確描述建筑墻壁各層內的水分擴散和傳熱的濕熱軟件,如WUFI和Delphin。由于基于區(qū)域間連續(xù)性方程,在數(shù)學上相較于有限體積/元建模更易求解且簡單的網(wǎng)格劃分使得計算時間更少。對于極端寒冷的氣候,Wang等[51]使用節(jié)點法進行濕熱模擬,以分析防潮閣樓結構。他們認為室內濕度負荷水平對不通風閣樓房間的濕熱條件有著較大的影響。Moujalled等[52]在墻和圍護結構兩個尺度上對建筑物的濕熱性能進行了數(shù)值評估,在實際的氣候條件下找到更真實的麻石灰混凝土濕熱特性。

WUFI是一個基于Künzel等[53]在Fraunhofer建筑物理研究所開發(fā)的濕熱包絡計算模型的整體模型。它以相對濕度作為的驅動力描述了多孔材料內水分源和匯、毛細作用、擴散以蒸汽帶解吸的液態(tài)水遷移過程并使用CVM來推導數(shù)學模型的數(shù)值解。Straube[34]使用WUFI 1D對3個不同高度位置的地下室墻體進行了一維濕熱模擬,使用測量的溫度變化最為外界條件,并假設土壤RH為100%。然而由于二維熱流的顯著影響和缺少混凝土干燥機制的描述,使得模擬并不成功。Pallin等[29]使用WUFI模擬了20 m深的一維土壤柱,地面平均溫度假設為9.5 ℃,在20 m深處典型土壤年溫度變化為0.2%。通過分析草的紋理,建筑物等對傳熱系數(shù)的影響分析后,認為實測和模擬之間有較好的擬合。但模擬的過程中缺少:①土壤表面降水的液體吸收;②下邊界的水分轉移;③土壤表面?zhèn)鳠嵯禂?shù)(隨多種因素變化);④積雪問題。Torres等[31]和Pazderka等[32]都使用了WUFI-2D 對墻基通風的實驗結果進行了驗證,不同的是前者的模擬結果是相對濕度而后者為水分含量,兩者的模擬結果都很好的驗證了實驗結果。Goldberg等[54]對2D的地下室墻和周圍的土壤進行了濕熱模擬。將土壤分為3個區(qū)域,以研究不同的土壤類型和測量的水分含量的影響。但WUFI未能求解水分輸運方程。

3.2 計算流體動力學(CFD)法

在節(jié)點法中空氣被認為是均勻的,而在存在低空氣循環(huán)區(qū)域的情況下,基于混合空氣近似的HAM模型可能會提供錯誤的結果。實際上,在HAM傳輸方法中,只有墻的一部分可能用于計算,而其余部分實際上可能處于非活動狀態(tài)。而計算流體力學(CFD)方法基于空氣體積的精細空間離散化和空氣中質量和能量的精確守恒方程的模型,可以高精度地計算規(guī)定體積內的溫度和蒸汽壓力標量場。ANSYS FLUENT和COMSOL Multiphyciscs等其他CFD工具提供了各種各樣的數(shù)值解方法且不局限于建筑物的能源模擬。

Comsol Multiphyciscs的多孔建筑材料熱濕傳遞模塊中,相對濕度(RH)梯度提供了水分傳遞的可能性,但Comsol的PDE模塊還允許實現(xiàn)由任何類型的數(shù)學方程定義的傳輸現(xiàn)象,從而增加其作為模擬工具的增強功能[55]。Ma等[25]使用了COMSOL PDE模塊建立二維的熱濕瞬態(tài)方程,來模擬遮地圍護結構的熱濕特性,研究表明,COMSOL可以實現(xiàn)溫度場與濕度場的耦合并有效解決非線性問題。他們認為在研究多孔材料傳熱時,在某些情況下應當考慮熱濕耦合傳熱模型。Saaly等[27]使用COMSOL Multiphysics研究了在寒冷氣候下,考慮土壤中的凍融循環(huán)通過土壤可變導熱系數(shù)當溫度下降到零度以下是,不同保溫措施下的地下圍護結構熱性能。研究說明了如何使用COMSOL利用簡化的凍結土瞬態(tài)熱傳導控制方程模擬二維和三維的地下室墻壁、地板和周圍土壤傳熱以及實現(xiàn)土壤的變化的導熱系數(shù)。韋明等[30]利用COMSOL Multiphysics系數(shù)形式偏微分方程模塊構建了以相對濕度為驅動勢的熱濕耦合方程,研究了在時變邊界條件下離壁式襯砌和貼壁式襯砌的深埋地下工程的熱濕傳遞過程。結果表明,空氣夾層的存在有利于降低地下工程的熱負荷以及混凝土墻體表面與室內空氣的傳濕量。

然而,CFD方法通常以幾千到幾百萬個網(wǎng)格對定義的整體進行劃分。對于每個控制體積,可以使用有限元或有限體積法求解守恒方程?；贑FD法[大渦模擬(LES)或直接數(shù)值模擬(DNS)]求解N-S(Navier-Stokes)方程是該類軟件的基礎之一。當需要對建筑物進行詳細的2D和3D建模的,精細的網(wǎng)格會導致數(shù)值求解需要消耗大量的時間成本。此外,在沒有流體動力學、性質和求解設置的相關知識下使用它會十分困難。Piscia等[56]使用耦合能量平衡模擬(ES)和CFD的方法研究溫室夜間氣候以減少CFD的較高計算成本。

3.3 方法對比分析

節(jié)點法在數(shù)學上比有限體積/單元建模更容易求解,而且因為不需要復雜的網(wǎng)格和計算步驟,所以使用起來也更簡單且計算時間更少。Yildiz[57]通過文獻檢索發(fā)現(xiàn)WUFI在學術研究中作為濕熱模擬工具是使用最多的程序。Kim[58]綜述了用于建筑物濕熱評估的不同濕熱工具,確定了WUFI和Delphin作為高級專業(yè)軟件可以評估建筑中熱量、水分和空氣并提供1、2、3維模型。但同時他們也指出了節(jié)點法由于數(shù)學模型的局限性無法解決建筑組件中的空氣流動和地下水的吸收。Nusser等[59]也認為、限制對控制方程的訪問是WUFI的缺點并在COMSOL的偏微分方程接口中實現(xiàn)了WUFI的物理方法。Defo等[60]對比Delphin、WUFI、hygIRC和COMSOL在不同氣候條件下的對于墻體熱濕響應后,認為WUFI相較于COMSOL傾向于在建筑結構積累更多的水分。Ozolins等[61]通過對比兩種不同的多層壁的實現(xiàn)結果來證明WUFI Plus和COMSOL數(shù)值模擬的差異,結果表明,COMSOL的計算結果與實驗結果吻合良好,且當面對較高濕分風險時,COMSOL能夠更好的估算多層建筑中的霉菌生長風險。

在WUFI的模擬過程中,降水通常是作為一個邊界條件,通過降水在外部邊界上施加濕度負荷并影響最接近邊界的網(wǎng)格然后再進一步在組件內分布。然而當模擬土壤以及地下建筑圍護結構等高相對濕度的問題區(qū)域時,它將不再有進一步的水分緩沖能力;因此,多余的水分就被忽略了。且對于主要相對濕度為100%RH的土壤和地下圍護結構,Delphin和WUFI面對密集的網(wǎng)格集中,計算速度會大大減慢。因此,復雜的計算模型需要較長的計算時間,并往往導致不穩(wěn)定的收斂。顯然,使用WUFI和Delphin模擬土壤的熱濕耦合可能不是最好的選擇。因為這些工具是為多孔建筑材料、組件和建筑圍護結構量身定制的。Delphin和WUFI都是由特定的偏微分方程發(fā)展而來,在數(shù)據(jù)輸入(氣候數(shù)據(jù)、材料屬性)上為用戶提供較高的便利,但在模型更改或組合上通常不是十分靈活。而類似于COMSOL這類商用計算軟件則可以更加靈活的模擬任何耦合偏微分方程系統(tǒng)。同時,大多數(shù)其他的HAM模擬軟件在都面臨著工具自身無法解決的難題。暖通空調組件和控制器(秒量級)和建筑響應(小時級)之間在時間常數(shù)上的差異會導致較長的模擬時間[62]。Janssen等[28]所述在降雨事件開始時,非常小的時間增量是必要的,以實現(xiàn)對降水的吸收和排水的精確模擬。在其他時候可能需要使用較大的時間步,來減少計算時間。此外,由于缺乏集總參數(shù)工具(包含內部2D和3D有限元方法(FEM)功能)導致在建筑模型中集成基于2D和3D幾何模型(建筑物的氣流或HAM響應)時出現(xiàn)問題[63]。Sally等[27]使用COMSOL對地下建筑進行了2D和3D的熱模擬。Ma等[25]研究認為,COMSOL可以有效地耦合溫度場和濕度場。然而,COMSOL這類軟件往往需要用戶提供大量的知識、資源在求解器中來實現(xiàn)以及氣候數(shù)據(jù)、邊界條件和材料屬性不是預先定義的,往往需要在其他專業(yè)軟件或是文獻中獲取。同時復雜的三維模型計算往往會導致模擬時間過長而無法實際應用。此時,采用耦合軟件法可能是一個有效的解決方法。Zhai等[64]將這種“協(xié)同模擬”方法描述為靜態(tài)或動態(tài)耦合?？娦∑降萚65]利用MATLAB和COMSOL協(xié)同模擬了地下工程圍護傳熱簡化計算方法和熱負荷的動態(tài)變化。Schellen等[66]利用COMSOL的Simulink功能聯(lián)合MATLAB進行了熱濕耦合方程的構建,研究了在地下水的影響下地基中的水分上升。Van Schijndel等[67]通過研究得出SimuLink和COMSOL的聯(lián)合仿真可以解決大量的綜合熱濕耦合(HAM)問題。上述的研究所涉及的耦合軟件法模擬中的幾何模型以及邊界條件都較為簡單。相關的文獻中此方法應用于地下建筑這類較為復雜的情況中的研究較少,還有待進一步研究。

4 邊界條件

除熱濕傳遞耦合方程的物理模型外,邊界條件對于數(shù)學模型求解結果的輸出質量也至關重要。對于地下建筑面臨著內外邊界上的各種水分問題挑戰(zhàn)。圖1[68]中顯示了半地下室可能面對的主要水分問題。

對于節(jié)點法和CFD法,內部氣候、外部環(huán)境和組件之間的邊界條件是至關重要的。控制體積(流體或固體域)時邊界處界面相互作用的適當預處理對于有效建模至關重要。表3中介紹了相關研究中對于不同邊界條件的設置。

表3 邊界條件設定總結Table 3 Summary of Boundary Condition Setting

4.1 室內邊界條件

溫度和相對濕度的動態(tài)值定義室內環(huán)境氣候,也稱為室內環(huán)境。常見的室內邊界條件設置經驗如下。

(1)假設地下空間由暖通空調系統(tǒng)進行空氣調節(jié),室內溫度和相對濕度設置為恒定值。

(2)使用溫濕度探測器測量一定時間內的室內溫度和相對濕度分布。

(3)以年為周期波動的溫濕度,且夏季最高氣溫不得高于30 ℃,冬季不得低于10 ℃,相對濕度控制在50%～80%。

然而,Shi等[22]假設最初的室內溫度和相對濕度分別為20 ℃和50%后,通過在不同換氣次數(shù)下為期兩周的初步模擬,以減少初始條件對模擬結果帶來的影響。Fedorik等[26]根據(jù)《建筑水和濕氣處理指南》(RIL 107—2012)[69]從室外條件根據(jù)建筑物的使用情況形成的,具有恒定的室內溫度和額外的水蒸氣質量濃度(3個單獨的水平)。

4.2 室外環(huán)境條件

室外環(huán)境條件直接影響地面附近的土壤區(qū)域和地表附近的地下建筑。對于處于地下較深位置的地下建筑,室外大氣環(huán)境條件的影響往往會被忽略。而淺層地下建筑的室外大氣環(huán)境條件確定為:當?shù)卦O計參考標準年、氣象站的氣象數(shù)據(jù)、測量。然而,簡單地使用大氣溫度和相對濕度作為邊界條件不能滿足土壤實際外部氣候環(huán)境的復雜性?？偨Y了從上述文獻中獲得的室外氣候使用經驗:①考慮地下室墻壁的干燥能力時,不應忽略太陽輻射,因為太陽輻射會導致水分向內轉移;②土壤表面降水的液體吸收;③下邊界處的水傳輸(無限地面);④與小時數(shù)據(jù)相比,日平均氣候數(shù)據(jù)的使用對土壤中的熱量和水分傳輸沒有顯著影響;⑤地表熱平衡包括顯熱交換H、太陽和長波輻射Rt,以及通過蒸發(fā)LE和降水HP傳遞的顯熱和潛熱,而水分平衡包括降水P和蒸發(fā)E;⑥物理域的上表面暴露于太陽和長波輻射、對流熱和水分傳遞;⑦相鄰房屋的陰影可能導致模擬和測量的溫度差異。

充分確定地下室圍護結構的地面以下的外部邊界條件非常具有挑戰(zhàn)性,因為這些條件隨地下高度而變化,并取決于若干不同因素(如土壤成分、不同含水量和濕熱特性、圍護結構的熱阻、室內溫度和外部氣候因素)。通過實際測量獲得的數(shù)據(jù)可能是準確的,但在布置儀器時,需要移除地下室周圍的土壤,對于周圍有遮擋的建筑物來說,移除土壤以獲得實際測量值變得十分困難。因此,邊界條件(或初始條件)可通過相鄰土壤和地下圍護結構的兩步法模擬確定。這可以通過以下方式實現(xiàn)。

(1)通過由平均溫度和靜水力矩水頭初始化的一維熱濕模擬計算溫度和矩陣水頭分布,從而初始化二維土壤完整的熱濕耦合模擬。

(2)計算地下土壤溫度分布,并根據(jù)輸入數(shù)據(jù)模擬樣品之間的熱濕環(huán)境,作為相應吸濕建筑材料的內部和外部邊界。

(3)僅模擬整個幾何體的溫度分布,作為后續(xù)建筑圍護結構熱和濕度模擬的熱邊界條件。

(4)僅使用熱模擬,但考慮了可變熱導率。

Janseen等[28]的研究中初步的一維模擬需要持續(xù)15年才能確保穩(wěn)定的周期解。隨后的二維模擬鏡像需要10～15年的模擬間隔才能達到穩(wěn)態(tài)周期狀態(tài)。通過長時間模擬來獲得真實的初始條件和熱濕傳遞方程中的大量參數(shù),計算效率往往是熱濕模擬中最為基本的關注點。因此,在后續(xù)研究中,熱模擬主要用于確定土壤溫度分布,作為地下圍護結構的熱外部邊界條件。Geving等[68]使用兩步方法對地下室圍護結構改造進行二維濕熱模擬。傳熱模擬包括了較大的土壤去。由于使用了飽和土壤的恒定熱特性并忽略了太陽輻射和積雪,導致冬季地下室墻體熱損失和冷凝風險的保守估算。Fedorik等[26]在研究不同年代混凝土地下室墻體設計的多種翻新策略的影響中也使用相似的計算方法。Saally等[27]雖僅對地下建筑進行了熱模擬,但比較了當土壤的熱特性(如熱導率、熱容和密度)恒定時以及當土壤的熱學特性隨凍融變化時的熱損失,研究發(fā)現(xiàn),相較于采用恒定的土壤熱特性時,地下建筑的能量損失會增加34%。這是因為當孔隙水開始凍結時,土壤的熱導率增加,經過地下室墻壁和樓板傳遞到土壤的熱量增加。Janssen等[28]將耦合熱模擬與線性熱模擬進行了比較,發(fā)現(xiàn)土壤表面溫度的大耦合振幅,存在顯著偏差;土壤導熱系數(shù)受到含水量的影響而發(fā)生變化;液體濕傳遞的顯熱對流。盡管土壤區(qū)域的簡化熱模擬是一種容易實現(xiàn)的方法,但由于水分的高度影響,不能更好地模擬土壤中的溫度分布。然而,在現(xiàn)階段的研究中還尚未解決土壤的全熱和水分模擬或簡化熱模擬對地下圍護結構的熱和熱行為的影響。相對于建筑物的熱和吸濕耦合,土壤的熱和濕耦合顯然需要研究者們更多的知識(土壤科學)和時間,因此通過簡化土壤域來關注建筑物圍護結構可能是更好的方法。

5 結論

正確模擬地下圍護結構的濕熱特性是一項具有挑戰(zhàn)性的工作。與地上建筑不同,地下建筑被土壤緊密包裹,這使得難以通過直接測量獲得地下圍護結構外部的邊界條件。然而,使用模擬獲得的邊界條件通常不能還原實際的熱濕邊界條件。往往不正確的模擬可能導致不充分的結論,這可能導致在建筑在建造和執(zhí)行節(jié)能改造時采取不當?shù)男袨?。通過對上述對數(shù)值方程、計算工具和邊界條件設置的分析,得出以下結論。

(1)地下建筑的熱濕傳遞與環(huán)境因素(室外氣候、地下水、土壤種類和室內熱濕負荷等)和自身構造相關,眾多因素影響下導致地下建筑的熱濕傳遞難以預測和計算。目前中外學者對于多孔介質熱濕耦合傳遞理論已給出了較為成熟的方程,但是上述數(shù)值模型為針對單一區(qū)域(多層建筑圍護結構或土壤)而建立。因此,未來的工作中應當開發(fā)耦合土壤域和建筑域的數(shù)值理論模型。

(2)傳統(tǒng)的基于熱濕耦合機理發(fā)展出的節(jié)點法計算工具,難以滿足地下建筑復雜的計算情況,如幾何模型中較大的土壤域等?；谟嬎懔黧w力學法的COMSOL等先進的商用模擬軟件可能是未來研究地下建筑熱濕耦合傳遞的首選工具。其對各類的通量的詳細分析、體積網(wǎng)格化可考慮到不同物理參數(shù)的強梯度以及結合了質量守恒、動量和能量方程以及流體動力學狀態(tài)方程,能夠滿足地下建筑的復雜環(huán)境。同時,COMSOL與其他軟件接口(如MATLAB)的耦合可以解決計算時間較長的問題并為后續(xù)數(shù)值模型發(fā)展提供。

(3)地下建筑熱濕耦合模擬的挑戰(zhàn)在于土壤成分和水分的不確定性、土壤表面的雨水沉淀和水分遷移以及較大土壤域帶來的巨大計算成本。兩步法的簡化雖然能夠適當估計熱損失,但可能不足以用于濕熱性能的研究。因此,通過忽略土壤域的存在,類似于用于確定地上建筑暖通空調設計和保溫設計的氣候數(shù)據(jù)的標準設計年,地下建筑的外部邊界條件可以根據(jù)土壤種類、氣候類型、地下建筑深度、建筑的服務用途、地下建筑的使用情況和地下建筑類型等各類因素給出類似規(guī)范來確定,以便于地下建筑中未來的熱濕模擬、結構設計和相關暖通空調設計。

(4)在論文篩選過程中,發(fā)現(xiàn)外國相關研究大多集中地下室,地下辦公室等人員會長期居住或逗留的地下建筑。而中國相關研究主要集中在有遮蔽的深地下工程(避難所、地下襯砌構造)或公共建筑(地下停車場)。這可能是因為中國的半地下室使用相較歐美作為生活和工作用途較少。然而,面對經濟發(fā)展和城鎮(zhèn)建筑面積不斷增加[70],建筑地下空間將會面臨更多地在生活和工作上的需求。因此,中國住宅半覆蓋式地下建筑的熱濕性能受到環(huán)境因素(土壤、氣候條件、地下水等)的影響有待進一步研究。