歐陽(yáng)煥英

中鐵第五勘察設(shè)計(jì)院集團(tuán)有限公司 北京 102600

隨著全球經(jīng)濟(jì)的快速發(fā)展和城鎮(zhèn)化水平的提高，城市建筑能耗逐年增長(zhǎng)，加上北方地區(qū)的供暖能耗占比較大，隨之導(dǎo)致的環(huán)境污染問(wèn)題日趨嚴(yán)重，由此近些年國(guó)家大力提倡建筑節(jié)能和環(huán)保，并積極推廣綠色建筑在全國(guó)范圍內(nèi)的應(yīng)用。我國(guó)幅員遼闊，太陽(yáng)能資源較為豐富，而充分利用太陽(yáng)能有效降低建筑供暖能耗正成為主要方向之一。以太陽(yáng)房為例可有效利用太陽(yáng)能供暖，節(jié)能環(huán)保，但其數(shù)值研究和成果多基于二維或三維的簡(jiǎn)化模型，很少對(duì)整體房屋實(shí)體建模，并進(jìn)行三維非穩(wěn)態(tài)數(shù)值模擬分析。

近些年來(lái)，國(guó)內(nèi)外許多學(xué)者都對(duì)太陽(yáng)房的數(shù)值模擬進(jìn)行了研究。2013年，Basak Kundakci Koyunbaba和Koray Ulgen等人[1]通過(guò)建立太陽(yáng)房的二維模型，得出模擬結(jié)果預(yù)測(cè)了BIPV型太陽(yáng)房墻體系統(tǒng)的溫度和速度分布。沈嬌和李德英[2]等人采用Fluent二維穩(wěn)態(tài)模擬改變太陽(yáng)房通風(fēng)口的尺寸大小，分析了太陽(yáng)房?jī)?nèi)溫度場(chǎng)和速度場(chǎng)的分布，得出結(jié)論合理選擇太陽(yáng)房墻體的幾何特征參數(shù)，有利于提高太陽(yáng)房墻體的集熱性能。2018年，趙志強(qiáng)[3]等人建立了相變材料的太陽(yáng)房墻體進(jìn)行二維非穩(wěn)態(tài)傳熱狀態(tài)空間模型，模擬數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)吻合較好。

1 太陽(yáng)房模型建立與驗(yàn)證

鑒于ANSYS模擬軟件諸多的優(yōu)越性和其他硬件設(shè)施的比選，本文擬采用ANSYS16.0進(jìn)行模型建立和模擬計(jì)算。選取基于青藏線試驗(yàn)工程的太陽(yáng)能房屋的物理模型進(jìn)行整體建模，進(jìn)行三維非穩(wěn)態(tài)模擬，并基于青藏線太陽(yáng)能房屋的測(cè)試數(shù)據(jù)驗(yàn)證本文才用的數(shù)值模型和數(shù)值方法的有效性，以期得到特朗勃墻建筑研究更加合理有效的數(shù)值方法。太陽(yáng)房物理模型及結(jié)構(gòu)部件如下：

圖1 太陽(yáng)房系統(tǒng)模型示意圖

1.1 數(shù)值模型

1.1.1 物理模型

青藏線課題組前期對(duì)站后試驗(yàn)工程的太陽(yáng)能房及其性能進(jìn)行了實(shí)驗(yàn)測(cè)試和分析，太陽(yáng)房為對(duì)流環(huán)路式系統(tǒng)。該試驗(yàn)工程建筑房間正南方向布置?；诖?，本文模擬房間簡(jiǎn)化物理模型如圖1所示，目標(biāo)房間沒(méi)有與室外相連的墻體，長(zhǎng)5.1m，寬3.6m，高3.4m；空氣通道寬度為150mm，房間附有窗戶，太陽(yáng)房墻體厚度為240mm。

1.1.2 數(shù)學(xué)模型

采用的數(shù)學(xué)模型為三大方程：質(zhì)量守恒方程、動(dòng)量守恒方程和能量守恒方程。簡(jiǎn)化為浮力驅(qū)動(dòng)的自然對(duì)流傳熱問(wèn)題，自然對(duì)流采用標(biāo)準(zhǔn)的Standard κ-ε湍流模型，無(wú)其他內(nèi)熱源并且為不可壓縮流體，忽略流體中的粘性耗散，同時(shí)認(rèn)為Boussinesq假設(shè)成立。

1.2 邊界條件

1 ）太陽(yáng)房的室外溫度定義為逐時(shí)變化，通過(guò)udf簡(jiǎn)單編寫。

2 ）玻璃蓋板與窗戶：采用第三個(gè)邊界條件，設(shè)置傳熱系數(shù)K和熱流都為常量；窗戶：采用第三個(gè)邊界條件，設(shè)置對(duì)流傳熱系數(shù)K為常數(shù)，由太陽(yáng)輻射通過(guò)窗戶照射到室內(nèi)地面的熱流邊界條件在后面給出。

3 ）集熱墻左側(cè)：采用第一個(gè)邊界條件，定義熱流逐時(shí)變化，將測(cè)試數(shù)據(jù)編譯到udf中。

4 ）集熱墻右側(cè)：定義為Couple邊界條件。表示此交界面兩側(cè)能夠?qū)帷?/p>

5 ）交界面interface和樓板：太陽(yáng)房模型中定義的interface接觸面均采用絕熱邊界條件。樓板也使用絕熱邊界條件。

6 ）其他墻面設(shè)置：地面，由于從窗戶照進(jìn)室內(nèi)地面上的熱量一部分被地面吸收，另一部分反射給其他墻面，已知窗戶的透射率為0.8，地面吸收率為0.3，反射率為0.7。由窗戶的透射率和地面的反射率求得周圍墻面的得熱量，根據(jù)面積比例得出各墻面的熱流量，將墻面設(shè)置為第一個(gè)邊界條件，即定熱流邊界條件。地面定義為絕熱邊界條件。

7 ）定義其他條件：采用初始測(cè)試數(shù)據(jù)，初始溫度設(shè)置為268K，初始空氣流速為0m/s。為使計(jì)算結(jié)果盡量接近實(shí)際值，殘差曲線和能量曲線結(jié)果必須在允許的誤差范圍內(nèi)，其殘差曲線＜10-3，能量曲線＜10-6。為了保證殘差曲線在該范圍內(nèi)，經(jīng)過(guò)多次調(diào)試，模擬最終采用的時(shí)間步長(zhǎng)定為0.1s，時(shí)步為4000。

2 數(shù)值方法驗(yàn)證

為了驗(yàn)證數(shù)值模型和方法，選取青藏線測(cè)試數(shù)據(jù)中1月17號(hào)的9:35、10:35、12:35、16:00、16:35幾個(gè)時(shí)刻的測(cè)試數(shù)據(jù)進(jìn)行對(duì)比，對(duì)比參數(shù)有：上通風(fēng)口溫度和風(fēng)速、下通風(fēng)口溫度和風(fēng)速和室內(nèi)溫度。

驗(yàn)證結(jié)果如圖2~6所示。從圖中可知，模擬數(shù)據(jù)基本上符合和實(shí)驗(yàn)數(shù)據(jù)，上下通風(fēng)口的風(fēng)速和溫度的曲線趨勢(shì)一致，大部分?jǐn)?shù)據(jù)誤差在10%左右，最大誤差也在20%以內(nèi)；從圖中可以發(fā)現(xiàn)模擬結(jié)果的上下通風(fēng)口的風(fēng)速大于實(shí)驗(yàn)風(fēng)速，分析原因是：在靠近壁面的區(qū)域，湍動(dòng)能被強(qiáng)烈的耗減，而模擬中采用κ-ε為一種高雷諾數(shù)模型，對(duì)該靠近壁面的區(qū)域不再適用，所以導(dǎo)致模擬數(shù)據(jù)的速度普遍高于測(cè)試數(shù)據(jù)中的速度。由于模擬結(jié)果的山下通風(fēng)口的風(fēng)速高于實(shí)驗(yàn)結(jié)果的風(fēng)速，所以由下通風(fēng)口帶入上通風(fēng)口的熱量在相同時(shí)間內(nèi)就多，所以導(dǎo)致上通風(fēng)口溫度也偏高；從圖中可知，下通風(fēng)口的模擬溫度略基本上趨近于測(cè)試溫度，只有在太陽(yáng)輻射最高的時(shí)刻12:35和早晚時(shí)刻兩者才有些許誤差。模擬的室內(nèi)溫度稍低于測(cè)試的室內(nèi)溫度，分析原因?yàn)椋河纱皯粽者M(jìn)室內(nèi)的熱量，沒(méi)有直接施加給室內(nèi)環(huán)境，而是給地面和周圍墻壁一個(gè)很小的定熱流，導(dǎo)致由窗戶照進(jìn)室內(nèi)的熱量經(jīng)過(guò)室內(nèi)地面反射和吸收之后，熱量所剩無(wú)幾，致使室內(nèi)溫度低于測(cè)試值。綜上所述，本文模擬與青藏線工程測(cè)試數(shù)據(jù)吻合較好。

圖2 上風(fēng)口風(fēng)速對(duì)比

圖3 上風(fēng)口溫度對(duì)比

圖4 下風(fēng)口風(fēng)速對(duì)比

圖5 下風(fēng)口溫度對(duì)比

圖6 室內(nèi)溫度對(duì)比

3 結(jié)論

1 ）本文采用分段函數(shù)自定義的方法對(duì)逐時(shí)室外溫度進(jìn)行定義，將每十分鐘的室外溫度值用VS軟件進(jìn)行編輯，采用UDF將逐時(shí)邊界條件導(dǎo)入到ANSYS中，計(jì)算特定時(shí)刻的模擬結(jié)果；玻璃蓋板的左側(cè)和窗戶設(shè)置為第三類邊界條件；玻璃蓋板右側(cè)設(shè)置為第一類邊界條件；集熱墻的左側(cè)采用熱流的分段函數(shù)自定義；集熱墻的設(shè)置為采用Couple邊界條件，即考慮墻體的導(dǎo)熱；樓板采用第一類邊界條件，其他墻體采用定熱流的絕熱邊界條件。

2 ）本文將青藏線工程測(cè)試的數(shù)據(jù)導(dǎo)入特朗勃墻模型中進(jìn)行計(jì)算，并與測(cè)試結(jié)果進(jìn)行了對(duì)比，完成了對(duì)數(shù)值模型和方法的有效性驗(yàn)證，為以后太陽(yáng)房三維模擬提供了方法借鑒。