呂　潔， 張　瓊， 周　勃， 王　慧

(沈陽工業(yè)大學 建筑與土木工程學院， 沈陽 110870)

?

自然通風作用下計算機房夏季熱環(huán)境數(shù)值模擬*

呂潔， 張瓊， 周勃， 王慧

(沈陽工業(yè)大學 建筑與土木工程學院， 沈陽 110870)

為了研究自然通風作用下計算機房夏季熱環(huán)境能否滿足學生舒適度的需要，通過對計算機房實地測量，利用測量數(shù)據(jù)建立數(shù)學模型，并以部分測量參數(shù)作為模型的邊界條件，利用CFD軟件對計算機房進行數(shù)值模擬.將模擬結果與實測數(shù)據(jù)進行對比，表明計算機房內溫度場和速度場的數(shù)值模擬結果與實測結果趨于一致.通過分析模擬結果，找到計算機房通風不良的位置，并對該處的設計缺陷提出相應建議，對計算機房的熱舒適性做出分析，且對機房空調設計給出具體建議.結果表明，自然通風作用下計算機房夏季熱環(huán)境不能滿足學生舒適度需要，且計算機房應當安裝空調.

自然通風； 計算機房； 降低能耗； 熱環(huán)境； 舒適度； 計算流體力學； 數(shù)值模擬； 實地測量

能源是人類生存和發(fā)展的重要基礎資源，開發(fā)和利用可再生能源對環(huán)境保護、減少溫室氣體排放、優(yōu)化我國能源組織結構等方面具有十分重要作用[1].事實證明，建筑能耗與環(huán)境密切相關，找到合理的方式減少對空調的依賴是實現(xiàn)低能耗建筑的有效途徑，盡可能使用自然通風，有利于降低建筑能耗[2].計算機房放置大量的電腦設備，當多臺機器穩(wěn)定運行時會產生一定的熱量，室內通風不良將在一定程度上影響學生上課質量.因此，合理地利用自然通風取代空調設備的使用，不僅能夠降低能耗，還能提高學生在計算機房上課的舒適性[3].

計算流體力學CFD方法常被用來求解氣流及溫度分布[4].CFD數(shù)值模擬方法成本低，可對建筑自然通風情況進行預測，為建筑設計的完善提供依據(jù)[5].本文基于自然通風基本原理，通過數(shù)值模擬探討該地區(qū)夏季炎熱條件下自然通風是否滿足舒適度需要，為計算機房的建筑設計和布置提供實際指導.

1　CFD基本原理

CFD是一種由計算機模擬流體流動、傳熱及相關傳遞現(xiàn)象的系統(tǒng)分析方法和工具.CFD的基本思想是把原來在時間域和空間域上連續(xù)的物理量場用一系列離散點上變化值的集合來代替，通過一定的原則和方式建立起反映這些離散點上場變量值的集合及場變量之間關系的代數(shù)方程組，求解代數(shù)方程組獲得場變量的近似解.

CFD計算主要包括前處理、求解和后處理3個部分.前處理階段需要進行的工作包括定義所求問題的幾何計算域；將計算域劃分為多個互不重疊的子區(qū)域，形成由單元組成的網格；對所要研究的物理或化學現(xiàn)象進行抽象，選擇相應的控制方程；定義流體的屬性參數(shù)；為計算域邊界處的單元指定邊界條件；對于瞬態(tài)問題，指定初始條件.目前各CFD軟件采用的求解方法為采用簡單函數(shù)近似表示未知的流量變動；將近似函數(shù)代入流動控制方程和所得到的數(shù)學式進行離散化；求解代數(shù)方程.后處理階段需要進行的工作包括，由于計算機技術的不斷進步，CFD軟件提供的數(shù)據(jù)可視化技術和工具越來越多，如計算域和網格顯示、壓力云圖、溫度云圖、速度云圖、速度矢量圖等[6].

2　自然通風作用下CFD算法驗證及比較

2.1試驗測試

2.1.1目的與方法

本次試驗的目的在于精準測量計算機房室內熱環(huán)境，為接下來數(shù)值模擬提供邊界參數(shù)和試驗驗證參數(shù)，并為以后采取具體措施改善計算機房的熱環(huán)境提供參考.

一天內不同時刻機房內具有不同的熱環(huán)境，為了獲得具有代表性的數(shù)值，選擇在下午兩點進行數(shù)據(jù)測量.

2.1.2測試及儀器

本次測量中，考慮到不能影響學生正常上課，選擇在周末進行.測量過程中有兩位同學做畢業(yè)設計，將其作為測量對象考慮在內.本次測量的參數(shù)包括室內溫度、風速、相對濕度以及電腦和人的表面溫度.

測試儀器主要采用JT-IAQ室內熱環(huán)境舒適度測試儀和紅外線測溫儀.其中，JT-IAQ主要用來監(jiān)測、記錄、計算和顯示室內環(huán)境數(shù)據(jù)，并自動計算WBGT、PMV、PPD等重要熱環(huán)境參數(shù)，為暖通空調及建筑熱環(huán)境等研究和檢測提供準確可靠的數(shù)據(jù)，可同步測試萬向微風速、黑球溫度、濕球溫度、空氣溫度和相對濕度等，多通道同時采集、自動顯示和記錄數(shù)據(jù)，具有測試精度高，穩(wěn)定性好等特點.本次測量設置時間間隔為60 s.測試儀器的詳細參數(shù)如下：空氣溫度，測量范圍和精度為5～60 ℃，±0.5 ℃；空氣濕度，測量范圍和精度為0～100%RH，±1.5%RH；風速，0.05～5 m/s，±(0.03 m/s+2%讀數(shù)).

為了便于模擬，方便觀察，測點的布置必須滿足要求.由于在熱環(huán)境模擬中，溫度場占有重要地位，因此，測點的布置主要考慮溫度傳感器的位置.測點分布平面圖如圖1所示.

圖1　測點布置平面圖

2.2物理模型

本文選取典型的自然通風計算機房，室內未設置空調.計算機房長10.8 m(對應y軸)，寬14.6 m(對應x軸)，高3.7 m(對應z軸).房間內有3排桌子，78臺電腦，2名學生，前后共16個窗子.由于每排兩臺電腦之間距離很小，于是將兩臺電腦簡化為一個長方體；測量時房門一直處于關閉狀態(tài)，因此，將房門簡化為墻.機房的布置和物體的幾何尺寸如圖2和表1所示.

圖2　房間的幾何模型

表1　各個模型的幾何尺寸

2.3數(shù)學模型

為了使問題簡化，本文做出如下假設：

1) 機房內空氣低速，定常，不可壓縮并且符合Boussinesq假設[7]；

2) 機房內空氣為穩(wěn)態(tài)湍流；

3) 忽略固體壁面和物體之間熱輻射；

4) 空氣為輻射透明介質；

5) 不考慮漏風影響，認為室內密封良好.

機房內氣流流動屬于湍流流動，采用FLUENT軟件進行數(shù)值模擬，采用k-ε兩方程模型，控制方程有連續(xù)性方程、動量方程、能量方程、湍流脈動動能方程(k方程)和湍流能量耗散率方程(ε方程).

2.4邊界條件和模擬計算

2.4.1邊界條件

人和電腦為室內熱源，其中，人體體表溫度為32 ℃，實測電腦表面溫度為40 ℃.南窗設為VELOCITY-INLET，實測入口風速為0.74m/s，溫度為29.2 ℃，北窗設為OUTFLOW，其余邊界設為WALL.

2.4.2模擬計算

CFD模擬是從微觀角度，針對某一區(qū)域或房間，利用質量、能量及動量守恒等基本方程對流場模型進行求解，分析其空氣流動狀況[8].計算機房內空氣流動屬于湍流流動，采用標準的兩方程模型(標準k-ε模型).本文模擬自然通風，考慮到浮升力的作用，因此，在設置流體材料時選擇Boussinesq假設進行計算，選中gravity選項，在z方向輸入-9.81m[9].速度和壓力耦合采用SIMPLE算法[10].在離散版面中，壓力采用BodyForceWeighted，動量、能量和其余兩選項都選擇一階迎風格式.方程的收斂條件設為流動方程相對誤差在10-3，能量方程相對誤差在10-6.

2.4.3計算機房模型模擬結果

下午2點z=0.9m處的空氣流速和溫度分布如圖3、4所示.

圖3　z=0.9 m的空氣流速分布圖

圖4　z=0.9 m的空氣溫度分布圖

下午2點y方向空氣流速分布如圖5～8所示.

圖5　y=1.2 m的空氣流速分布圖

圖6　y=4.1 m的空氣流速分布圖

圖7　y=7.4 m的空氣流速分布圖

圖8　y=9.4 m的空氣流速分布圖

下午2點y方向空氣溫度分布如圖9～12所示.

圖9　y=1.2 m的空氣溫度分布圖

圖10　y=4.1 m的空氣溫度分布圖

圖11　y=7.4 m的空氣溫度分布圖

圖12　y=9.4 m的空氣溫度分布圖

3　試驗驗證及分析

3.1模擬值與實測值對比分析

為了對模擬結果做進一步的驗證，將現(xiàn)場測量結果與模擬結果進行比較.選擇y=4.1 m，z=0.9 m位置，按照測量點x=2 m，x=7 m，x=13 m處測得的空氣流速和溫度與模擬數(shù)值進行比較，結果如表2所示.

表2　實測值與模擬值比較

由表2可以看出，大部分模擬值與測量值吻合較好，相對誤差在10%以內，說明建立的物理模型是正確的.相對誤差較大的原因為，就空氣溫度而言，在x=7 m處，從圖9可以看出，人體散熱使附近空氣溫度升高，使得模擬值比實際測量值要高；就空氣流速而言，由于采用自然通風，機房內的氣流主要受到機房外氣流影響，并且室外的空氣方向和流速不定，對機房內的空氣流速影響較大.

3.2模擬結果分析

3.2.1速度場模擬結果分析

從空氣速度分布圖可以看出：

1) 在水平方向上，如圖3所示，從3個方面進行分析：

① 空氣流速.機房的空氣進出口兩處空氣流速相對較大，最高可達0.22 m/s，通風效果較好；在中間區(qū)域能明顯看出空氣流速急劇下降趨近于零，由此可以看出，該區(qū)域的通風效果是最不利的.

② 窗間墻.在氣流入口側，中間的寬窗間墻阻擋作用較大，對其正對座位有較大影響，空氣流速接近零.其他窗間墻較窄，在臨窗處有一定阻擋作用，但座位處空氣流速在0.11～0.16 m/s之間，影響不大.

③ 過道.兩邊的過道雖然都沒有阻擋，但過道寬窄差別導致流速有差別，較寬的過道氣流通暢，較窄的則較差.

針對上述因素提出如下建議：首先，該類型房間適宜采用較窄的窗間墻；其次，過道采用相同寬度，同時在對應窗戶中間增加小過道；最后，應當避免電腦桌成排擺放，并且電腦桌之間應留出一定間隙.

2) 在豎直方向上，隨著y值的增大，由南向北方向依次如圖5～8所示，其中，圖5～8分別在水平方向距離x軸1.2、4.1、7.4和9.4 m，圖6、7分別位于中間兩個過道位置.從上述圖中可以看到明顯的條紋分布層，可知空氣流速在豎直方向上分布相當不均勻，可從以下兩個方面進行分析：

① 最高風速.圖5～8最高風速出現(xiàn)位置依次為進風口、計算機上方空間、出風口.在經過這些位置時，最高風速從0.32 m/s降為0.096 m/s，這是由于空氣由入口進到房間后受到電腦的阻擋，氣流被限制通過，導致流速降低.

② 最低風速.從圖中可知，最低風速出現(xiàn)的z=0.4 m到z=0.8 m之間，介于電腦桌面和地面中間位置.以電腦為分界線，電腦以上空間空氣速度在0.05 m/s以上，以下位置空氣流速明顯降低，大部分約在0.006 m/s，少部分空氣流速為零.這是由于機房內電腦擺放過密，以至于氣流極少能到達電腦桌面以下的位置，另外電腦桌上有擋板，不利于氣流通過.

針對上述因素提出如下建議：首先，將窗臺降低，增大進風量；其次，將電腦桌上擋板去掉，電腦桌之間改變排列方式，減少對下部空氣阻擋，將氣流影響范圍增大.

3.2.2溫度場模擬結果分析

從空氣溫度分布圖可以看出，在水平方向上，由于機房內計算機大部分處于工作狀態(tài)，散發(fā)出一定的熱量，使得機房內溫度比室外溫度高.空氣從室外由窗口進入機房，由于室外溫度比室內低，經過一定的熱交換使得進風側空氣溫度降低，最低溫度為30.2 ℃.中間較寬的窗間墻對空氣雖然有一定阻擋作用，但該處空氣溫度仍有一定降低，為30.3 ℃，與最低溫度差距0.1 ℃.隨著水平方向的深入，溫度降低程度減少，并逐漸穩(wěn)定在30.4 ℃.

在豎直方向上，隨著y值的增大，在分別距離x軸1.2、4.1、7.4和9.4 m處取計算機房的空氣溫度分布圖，如圖9～12所示，其中，圖10、11分別在兩過道位置.從圖中能看到明顯的條紋分布層，由此可見，室內溫度分布很不均勻.從最高溫度和最低溫度分布范圍分析可知：

1) 最高溫度.除圖9最高溫度出現(xiàn)在寬窗間墻阻擋的位置外，圖10～12最高溫度主要出現(xiàn)在電腦和電腦桌擺放的區(qū)域，且隨著y值的增大，最高溫度由30.29 ℃增加到30.43 ℃.由此可知，寬窗間墻比窄窗間墻阻擋區(qū)域的溫度高出0.05 ℃，可見寬窗間墻不利于通風.

由于電腦散熱造成其所在區(qū)域溫度比其余位置高出0.1～0.2 ℃，隨著水平方向深入，高溫范圍逐漸增大，室外空氣已經不能有效帶走余熱，降溫效果不明顯.

2) 最低溫度.低溫區(qū)域主要分布在氣流入口處，隨著y值增大，低溫區(qū)域集中于遠離電腦的位置，且最低溫度從30.24 ℃增加到30.31 ℃，可見室外空氣通過熱交換吸收室內部分熱量，但影響不大.

3.3舒適度分析

根據(jù)美國ASHRAE55-2013標準，建議地面上方0.1～1.8 m之間垂直溫差不超過3 ℃.從圖10可以看出，在0.1 m處溫度在29.5～30.1 ℃之間，在1.8 m處溫度在29.8～30.4 ℃之間，垂直溫差最大為0.9 ℃，符合標準要求.

PMV-PPD是評價熱環(huán)境的重要指標，ISO7730對PPD的推薦值為0.1，對PMV的推薦值為-0.5～+0.5.由儀器JT-IAQ室內熱環(huán)境舒適度測試儀測量可得，計算機房內PMV的測量值在1.3～1.7之間，PPD的測量值在40.3～61.8之間.經過與ISO7730推薦值對比可知，PPD、PMV均大于推薦值，由此可知，計算機房的熱環(huán)境不能滿足學生舒適度要求.

4　結　論

本文利用FLUENT軟件對計算機房進行數(shù)值模擬，得到了室內空氣的速度場和溫度場，并與實際測量結果進行比較，分析得出模擬溫度場與實際測量結果擬合較好，模擬速度場和實際測量結果相對誤差較大，并對誤差大的原因做出了如下分析：

1) 計算機房存在建筑設計缺陷，建議將空氣入口處中間較寬的窗間墻變窄，同時降低窗臺；

2) 電腦桌的擺放存在問題，應該加大桌子間距，使氣流順利通過；

3) 自然通風的降溫作用有限，應減少機器臺數(shù)或設置空調，否則不能滿足計算機房學生舒適度要求.

另外，建議計算機房安裝空調，且應采用工位送風或下送上回形式，同時應設置機械排風系統(tǒng)，排風口應同時布置在電腦區(qū)域和頂部.由于本文重點研究自然通風對計算機房熱環(huán)境影響，因此，關于計算機房空調設計未作過多贅述和模擬分析.

[1]孫秀剛，石俊龍，程大磊.吉林可再生能源建筑應用示范項目測評分析 [J].沈陽工業(yè)大學學報，2015，37(4)：469-474.

(SUN Xiu-gang，SHI Jun-long，CHENG Da-lei.Eva-luation and analysis for building application demonstration projects of renewable energy in Jilin province [J].Journal of Shenyang University of Technology，2015，37(4)：469-474.)

[2]Taleb H M.Natural ventilation as energy efficient solution for achieving low-energy houses in Dubai [J].Energy and Buildings，2015，99：284-291.

[3]楊天文.自然通風建筑室內氣流組織優(yōu)化研究 [J].四川建筑科學研究，2012，38(4)：312-316.

(YANG Tian-wen.Optimization study on building indoor airflow organization on the condition of natural ventilation [J].Sichuan Building Science，2012，38(4)：312-316.)

[4]江向陽.廣東地區(qū)某建筑夏季室內熱環(huán)境的CFD仿真評價 [J].廣州建筑，2012，40(1)：38-41.

(JIANG Xiang-yang.CFD simulation of summer indoor thermal environment for a building in Guangdong province [J].Guangzhou Architecture，2012，40(1)：38-41.)

[5]李靜.典型布局下的建筑室內自然通風 [J].建筑技術，2012，43(12)：1124-1126.

(LI Jing.Indoor natural ventilation in typical layout [J].Architecture Technology，2012，43(12)：1124-1126.)

[6]孫幫成，李明高.ANSYS FLUENT 14.0仿真分析與優(yōu)化設計 [M].北京：機械工業(yè)出版社，2013.

(SUN Bang-cheng，LI Ming-gao.ANSYS FLUENT 14.0 the simulation analysis and optimization design [M].Beijing：Mechanical Industry Press，2013.)

[7]陶文銓.數(shù)值傳熱學 [M].西安：西安交通大學出版社，2001.

(TAO Wen-quan.Numerical heat transfer [M].Xi’an：Xi ’an Jiaotong University Press，2001.)

[8]龔波.自然通風的設計策略及模擬分析 [J].建筑熱能通風空調，2004，23(5)：30-34.

(GONG Bo.The design strategy and simulation of natural ventilation [J].Building Energy and Environment，2004，23(5)：30-34.)

[9]陳雨，許志浩，馬國川.關于自然通風CFD 算法幾點探討 [J].制冷與空調，2011，25(1)：78-81.

(CHEN Yu，XU Zhi-hao，MA Guo-chuan.The study about CFD algorithms of natural ventilation [J].Refri-geration and Air Conditioning，2011，25(1)：78-81.)

[10]董志超，婁君.某辦公室熱環(huán)境CFD模擬研究 [J].制冷與空調，2011，25(1)：102-106.

(DONG Zhi-chao，LOU Jun.Study on thermal environment of office with CFD simulation [J].Refrigeration and Air Conditioning，2011，25(1)：102-106.)

(責任編輯：鐘媛英文審校：尹淑英)

Numerical simulation of summer thermal environment in computer room under action of natural ventilation

Lü Jie, ZHANG Qiong, ZHOU Bo, WANG Hui

(School of Architectural and Civil Engineering, Shenyang University of Technology, Shenyang 110870, China)

In order to study whether the summer thermal environment in computer room under the action of natural ventilation can meet the need in the comfort degree of students, a mathematical model was established with the measured data in the computer room. The partial measurement parameters were used as the boundary conditions of the model, and the numerical simulation for the computer room was carried out with CFD software. Through comparing the simulation results with the measured data, it is noted that the simulated results of temperature field and velocity field in the computer room tend to be consistent with the measured results. Through analyzing the simulated results, the location with poor ventilation in the computer room could be found, and the corresponding suggestions for the design flaws were given. In addition, the thermal comfort of computer room was analyzed, and the specific suggestions for air conditioning design of computer room were given. The results show that the summer thermal environment in computer room under the action of natural ventilation can not meet the need in the comfort degree of students, and an air conditioning should be installed in the computer room.

natural ventilation; computer room; energy consumption reducing; thermal environment; comfort degree; computational fluid dynamics(CFD); numerical simulation; field measurement

2015-08-27.

中國博士后科學基金資助項目(2014M560220).

呂潔(1967-)，女，遼寧朝陽人，副教授，主要從事暖通空調應用等方面的研究.

10.7688/j.issn.1000-1646.2016.01.20

TU 111.1

A

1000-1646(2016)01-0115-06

*本文已于2015-12-07 16∶16在中國知網優(yōu)先數(shù)字出版. 網絡出版地址： http：∥www.cnki.net/kcms/detail/21.1189.T.20151207.1616.006.html