嚴樂 郭愛賓 鄧全亮 高維嘉

北京可視化智能科技股份有限公司

0 引言

計算流體力學（Computational Fluid Dynamics，CFD）軟件可有效地模擬計算流體的特性。由于其含有豐富的計算模型，能適應(yīng)較多場合，被越來越多地應(yīng)用到室外風場模擬、暖通空調(diào)系統(tǒng)的深化設(shè)計中，是建筑設(shè)計階段的有力工具[1-5]。然而在利用CFD軟件進行深化設(shè)計時，邊界條件的確定至關(guān)重要，它直接影響了計算結(jié)果的合理性和準確性。

本文以北京某公司總部大樓的多個相連通的大廳為算例，首先，根據(jù)建筑的物理參數(shù)、夏季基礎(chǔ)空調(diào)設(shè)計參數(shù)、夏季空調(diào)系統(tǒng)設(shè)計方案等，選取所要研究的物理模型及數(shù)學計算模型，借助CFD-FLUENT軟件對空調(diào)系統(tǒng)設(shè)計方案溫度分布、速度分布及熱舒適性進行三維數(shù)值模擬研究。研究探討了不同邊界條件對氣流組織設(shè)計方案的結(jié)果影響，分析優(yōu)選出較為完善方案并將其應(yīng)用到實際工程中，然后，對所研究的高大空間會議廳進行現(xiàn)場實際測試，并將完善后的氣流組織模擬所得結(jié)果與實際測試結(jié)果作對比驗證，以期為實際暖通空調(diào)的深化設(shè)計提供邊界條件的選擇參考方法。

1 物理及數(shù)值模型

1.1 物理模型

本文以北京某公司總部大樓的多個相連通的大廳為算例（圖1），該連通大廳由3部分組成，即南廳、東北廳和西北廳，南廳由1層起，東北廳地下1層起，西北廳2層起，南廳和東北廳在1層相通，東北廳和西北廳在2層通過2組連廊相通。計算對象為此三部分大廳連通形成的內(nèi)部空調(diào)空間，該空調(diào)空間南北向最大尺寸為即90 m，東西向最大尺寸為58 m，高度方向最大尺寸為25 m。三個大廳頂部分別設(shè)有1組天窗。

圖1 多連通大廳物理模型圖

表1 多連通大廳6組送風口參數(shù)

從空調(diào)設(shè)計角度，根據(jù)前期冷負荷計算，該空間共設(shè)置6個回風口、6組送風口，這6組送風口的參數(shù)見表1：第一組設(shè)置在南廳南半部分（hall-downward），其頂部在每隔2m布置有55個20cm口徑的圓形射流風口；第二組設(shè)置在南廳的北側(cè)墻面，為側(cè)送風口（2F-sidewall1），共15個；第三組設(shè)置在南廳和東北廳連通的1樓頂部（1F-downward），每隔2m布置有22個20cm口徑的圓形射流風口；第四組設(shè)置在東北廳和西北廳的2組連廊底部（CB-downward），共14個；第五組設(shè)置在2組連廊側(cè)面（CB-sidewall），共28個；第六組設(shè)置在西北廳的西側(cè)墻面（2F-sidewall2），共21個。

1.2 計算模型

本算例中的數(shù)學模型模擬計算采用k-ε兩方程紊流模型。在夏季，通過該大廳頂部的三組天窗的輻射傳熱和導熱形成的冷負荷為主要組成部分，CFD-Fluent軟件提供了四種不同的熱輻射模型（表2），本算例選取DO模型，以得到更準確的計算結(jié)果。

表2 四種熱輻射模型特點

表3 物理參數(shù)

表4 天窗熱物性參數(shù)

本算例計算夏至日室內(nèi)的空氣情況，在采用DO模型進行計算時，需要算例所處地點時間的太陽輻射照度、角度、經(jīng)緯度及時區(qū)等，本算例位于北京，夏至日時外界參數(shù)如表3所示。

2 不同邊界條件的分析

2.1 邊界條件1

送風18℃，連廊向下風口45°角對吹，回風自由出流，該廳四周房間為空調(diào)區(qū)域，故模擬計算時為絕熱條件，南廳南墻和北廳北墻選擇建筑設(shè)計節(jié)能規(guī)范中的墻壁厚度和傳熱系數(shù)，計算結(jié)果如圖2～6所示，可以看出，南廳距地面1.5 m高度速度場較為合理，但東北廳相距地面1.5 m高處的速度分布不均勻，并且局部過大，會有不適感。從溫度場來看，東北廳連廊處溫度較高。

圖2 南廳1.5 m高處速度場

圖3 東北廳B1層-4.5 m高處速度場

圖4 東北廳連廊處剖面速度矢量圖

圖5 1.5 m高處溫度場

圖6 東北廳連廊處剖面溫度場

2.2 邊界條件2

與邊界條件1相比，邊界條件2中，所有結(jié)構(gòu)參數(shù)參考節(jié)能規(guī)范，南、北墻為玻璃幕墻，加載人員發(fā)熱量，使得邊界條件的選取更接近實際，并且為了解決case1中溫度速度不均的問題，調(diào)整連橋處送風方向為45°相背向外，圖7～11為計算結(jié)果。

圖7 南側(cè)玻璃幕墻加載前后南廳地面溫度對比

圖8 連廊送風角度調(diào)整前后大廳剖面度度場對比

圖9 西北廳剖面速度場

圖10 大廳剖面溫度場

圖11 三個廳相對高度1.5m溫度場

從圖7中可以看出，當南墻改為玻璃幕墻后，南廳近玻璃幕墻且為太陽投射地面處溫度明顯有所升高，廳內(nèi)人員活動區(qū)域溫度分布較為均勻但稍高；圖8對比了連廊送風角度調(diào)整前后大廳剖面度度場，可看出，調(diào)整送風角度向外后人員活動空間速度分布更為均勻合理，并且從圖9可看出西北廳剖面速度場也較為合理；圖10為東北廳、西北廳兩廳東西向大廳剖面溫度場，人員活動空間溫度在26～29℃，整體剖面溫度分層現(xiàn)象明顯；圖11位三個廳相對高度1.5 m處的溫度場，為20～29℃，平均溫度26℃，在可接受的范圍之內(nèi)。

3 實測驗證

2016年7月，也是算例中總部大樓運行第一年，本文筆者組織測試人員對該大連通大廳進行了溫度的實測（邊界條件與case2相同，對應(yīng)圖11，約每隔10 m見方，共選擇50個溫度實測點，實測與計算結(jié)果對比如圖12），結(jié)果表明，溫度的實測結(jié)果和計算結(jié)果吻合度較高，說明了邊界條件選擇準確時，模擬計算結(jié)果的準確性更高。

圖12 實測與計算的溫度對比圖

4 結(jié)論

1）CFD中熱輻射模型的DO模型對小尺度到大尺度輻射計算都適用，且可計算非－灰度輻射和散射效應(yīng)，現(xiàn)今的計算機能力可以適應(yīng)稍微大的計算量。在計算時，需要的參數(shù)較多，但參數(shù)選擇準確時，計算結(jié)果也更為準確；

2）通過對本工程夏季空調(diào)工況下實測值和數(shù)值模擬數(shù)值的對比分析，可以驗證CFD數(shù)值模擬技術(shù)的可靠性，并證明CFD Fluent軟件能夠為高大空間類建筑空調(diào)系統(tǒng)優(yōu)化方案設(shè)計、預(yù)測氣流組織分布、評價熱舒適指標等。

3）一個良好的氣流組織和送風的角度密不可分，本文通過對不同送風角度的氣流組織型式進行數(shù)值模擬并定性分析，選定比較合理的氣流組織方案，為實際工程設(shè)計提供有利參考；

4）通過兩種邊界條件的計算對比，可看出在利用CFD Fluent軟件進行深化設(shè)計時，邊界條件的選擇一定要越接近實際越好，越可以準確輔助設(shè)計人員進行空調(diào)的優(yōu)化設(shè)計。本算例提供的一系列參數(shù)可以為實際暖通空調(diào)的深化設(shè)計提供邊界條件的選擇參考?？梢灶A(yù)測，在實際運行中，基于CFD Fluent軟件和準確的邊界條件的加載，能夠更好地預(yù)測空調(diào)運行狀態(tài)對空間舒適度的改變，以指導系統(tǒng)更合理、更節(jié)能運行。