辦公建筑室內(nèi)環(huán)境數(shù)值模擬*

周 勃， 宋勝波， 董立國

(1. 沈陽工業(yè)大學 建筑與土木工程學院， 沈陽 110870； 2. 特靈空調(diào)系統(tǒng)(中國)有限公司， 沈陽 110013)

隨著我國經(jīng)濟的不斷發(fā)展，人們對室內(nèi)空氣品質(zhì)的要求不斷提高.大量分析表明，人的一生中有超過80%的時間是在室內(nèi)度過的[1]，室內(nèi)空氣品質(zhì)的優(yōu)劣與人類健康密切相關，因此，選擇合適的送風方式對辦公類建筑室內(nèi)熱舒適性和空氣品質(zhì)的控制特別重要.

近些年來，計算流體力學(CFD)技術在室內(nèi)環(huán)境模擬方面得到了重大發(fā)展[2-6].利用計算流體力學技術進行數(shù)值模擬具有投資費用少和所用時間短等特點.關于甲醛污染物分布情況前人做了許多工作[7-9].本文以沈陽市某辦公類建筑中一會議室為研究對象，利用計算流體力學技術對會議室內(nèi)環(huán)境進行數(shù)值模擬，在不同送風方式下對室內(nèi)溫度場、污染物(甲醛)濃度場和速度場分布情況進行了數(shù)值分析.

1 物理模型和數(shù)學模型

1.1 物理模型

會議室在不同送風方式下物理模型如圖1、2所示，其空間結構為：長度為X=7.5 m，寬度為Y=6.8 m，高度為Z=2.7 m，北墻為外墻，東、南、西墻為內(nèi)墻.窗戶在北墻上，窗戶尺寸為高×寬=1.7 m×0.9 m.本文著重分析同側上送風下回風和頂送風下回風兩種不同送風方式對室內(nèi)環(huán)境的影響情況.模擬工況中，人員、桌子、送風口和回風口尺寸等相關參數(shù)如表1所示.

圖1 側送風房間的幾何模型Fig.1 Geometrical model for side air supply room

圖2 頂送風房間的幾何模型Fig.2 Geometric model for top air supply room

1.2 數(shù)學模型

本文使用CFD軟件進行數(shù)值模擬，計算過程中不考慮空氣與室內(nèi)各面摩擦情況，會議室內(nèi)空氣為穩(wěn)態(tài)湍流，不可壓縮且符合Boussinesq假設，采用標準的k-ε模型.會議室內(nèi)通風的數(shù)學模型控制方程包括連續(xù)性方程、動量方程和能量方程[10]，其表達式如下所示.

表1 室內(nèi)人員、桌子、送風口、回風口的相關參數(shù)Tab.1 Relevant parameters for indoor personnel，desk， air supply outlet and return air inlet m

連續(xù)性方程：

(1)

動量方程：

(2)

能量方程：

(3)

由于本文模擬了污染物(甲醛)從桌子表面向會議室內(nèi)擴散的過程，應用到組分傳輸模型，因此，控制方程還應包括組分質(zhì)量守恒方程，其表達式為

(4)

2 模擬結果及分析

影響室內(nèi)空氣品質(zhì)的因素有很多，包括送回風口的位置、數(shù)量、形式以及送風口的風速和送風溫度等.本文主要研究影響會議室內(nèi)熱舒適性和空氣品質(zhì)的主要因素，即送風口位置的變化對會議室內(nèi)溫度場、污染物(甲醛)濃度場和速度場的影響.本文截取Z=1.2 m(距離地面1.2 m為會議室內(nèi)人員正坐時呼吸區(qū)域的高度)和X=3.75 m截面來分別顯示會議室水平面和立面上溫度、污染物(甲醛)濃度和風速的分布情況，通過對比分析選擇出相對較好的空調(diào)系統(tǒng)送風方式.

2.1 不同送風方式對室內(nèi)溫度場的影響

圖3～6分別為側送風方式和頂送風方式下Z=1.2 m、X=3.75 m截面處的溫度場分布情況.由圖3可知，側送風情況下Z=1.2 m截面上溫度場分布相較頂送風方式下差異性較大，局部區(qū)域溫度低于23 ℃，會議室內(nèi)人員會有明顯冷感；靠近北墻一側由于受室外溫度影響，溫度分布趨于均勻.由圖4可知，頂送風情況下室內(nèi)溫度較均勻，會議室內(nèi)活動區(qū)域溫度保持在23.7～25 ℃的范圍內(nèi)，室內(nèi)舒適性較好.雖然存在局部區(qū)域溫度低于23 ℃，但范圍較小，對整體的影響不會太大.由圖5、6可知，采用側送風方式，在X=3.75 m截面上會議室內(nèi)西墻附近溫度較低，會議室后排就坐的人員會有明顯的冷感，不滿足ISO7730中規(guī)定的當人在固定位置從事輕工作時，有效溫度為23～26 ℃的要求.而頂送風方式下，由于送風口直接將空氣送至教室中間，因此，送風口下方部分區(qū)域溫度會稍低一些.由于頂送風情況下會議室內(nèi)溫度大部分處于23.9～25.5 ℃范圍之內(nèi)，因此滿足室內(nèi)舒適性要求.

圖3 Z=1.2 m側送風溫度分布Fig.3 Temperature distribution under sideair supply mode with Z=1.2 m

圖4 Z=1.2 m頂送風溫度分布Fig.4 Temperature distribution under topair supply mode with Z=1.2 m

2.2 不同送風方式對室內(nèi)污染物(甲醛)濃度場的影響

室內(nèi)甲醛污染帶來的危害已經(jīng)越來越被人們所熟知[11].圖7、8分別為側送風和頂送風方式下在Z=1.2 m截面處的甲醛污染物濃度分布情況.由圖7可知，人員在會議室開會時的呼吸區(qū)域和其他局部區(qū)域內(nèi)甲醛污染物濃度均未超過0.1 mg/m3，室內(nèi)甲醛污染物濃度滿足《室內(nèi)空氣質(zhì)量標準》(GB/T18883)規(guī)定的室內(nèi)空氣中甲醛含量不得超過0.1 mg/m3限值的情況.由于側送風方式空氣在送達西墻后出現(xiàn)下沉，室內(nèi)空氣循環(huán)擾動從而帶動局部甲醛污染物運動，因此，可以明顯看出會議室內(nèi)西墻附近區(qū)域內(nèi)甲醛污染物濃度值低于東墻側污染物濃度值.由圖8可知，頂送風方式下Z=1.2 m截面上甲醛污染物濃度分布情況相對優(yōu)于側送風方式，室內(nèi)污染物濃度最高1.055×10-7mg/m3，且相對高濃度區(qū)域較側送風要小.兩種送風方式下靠近北墻一側甲醛污染物濃度明顯高于靠近南墻一側，這是因為北墻一側與室外空氣接觸，有熱量透過北墻傳入室內(nèi)，從而加速了室內(nèi)靠近北墻一側桌子內(nèi)污染物(甲醛)的散發(fā).圖9、10分別為側送風和頂送風方式下在X=3.75 m截面處的甲醛污染物濃度分布情況.從圖中可以看出，兩種送風方式基本上均能夠滿足室內(nèi)空氣中對甲醛污染物含量的要求.由于甲醛污染物是從桌子表面散發(fā)出來，因此從圖中可知桌子表面甲醛濃度較大，隨著與桌面距離的不斷增大，甲醛濃度逐漸減少.

圖5 X=3.75 m側送風溫度分布Fig.5 Temperature distribution under sideair supply mode with X=3.75 m

圖6 X=3.75 m頂送風溫度分布Fig.6 Temperature distribution under topair supply mode with X=3.75 m

圖7 Z=1.2 m側送風甲醛濃度場Fig.7 Formaldehyde concentration under sideair supply mode with Z=1.2 m

圖8 Z=1.2 m頂送風甲醛濃度場Fig.8 Formaldehyde concentration under topair supply mode with Z=1.2 m

圖9 X=3.75 m側送風甲醛濃度場Fig.9 Formaldehyde concentration under sideair supply mode with X=3.75 m

本文對室內(nèi)甲醛污染物濃度的模擬結果與前人對室內(nèi)空氣中甲醛污染物濃度的監(jiān)測和調(diào)查結果進行了比較[12-13]，結果基本吻合.

2.3 不同送風方式對室內(nèi)速度場的影響

會議室工作區(qū)域內(nèi)送風速度分布情況的好壞是衡量熱舒適性的一個重要指標.如果工作區(qū)域內(nèi)氣流運動速度過快，不僅會影響工作區(qū)域內(nèi)溫度的分布，而且會使人有明顯吹風感，但是會提高甲醛的稀釋速度.圖11～14分別為兩種送風方式下Z=1.2 m、X=3.75 m截面處的風速分布情況.

圖10 X=3.75 m頂送風甲醛濃度場Fig.10 Formaldehyde concentration under topair supply mode with X=3.75 m

圖11 Z=1.2 m側送風速度場Fig.11 Velocity field under side air supplymode with Z=1.2 m

圖12 Z=1.2 m頂送風速度場Fig.12 Velocity field under top air supplymode with Z=1.2 m

在側送風方式下，氣流從側送風口射出沿射流方向送至人員工作區(qū)域，由圖11可知，由于氣流速度隨輸送距離的增加逐漸減小，再加上受重力影響，在氣流到達西側墻體時送風速度普遍小于0.089 84 m/s，Z=1.2 m截面上最大風速為0.251 5 m/s，工作區(qū)域內(nèi)風速均小于0.3 m/s，符合空調(diào)室內(nèi)有關速度的設計要求.為了與側送風方式下室內(nèi)的速度分布情況進行對比，本文截取了頂送風方式下Z=1.2 m截面模擬圖.頂送風情況下氣流從送風口直接送至會議室內(nèi)人員工作區(qū)，工作區(qū)域內(nèi)氣流的速度處于0.026 74～0.147 1 m/s之間，遠遠低于0.25 m/s，符合ASHRAE國際標準.由圖13、14對比分析可以得出，X=3.75 m截面上頂送風方式下送風速度更為均勻且均小于0.3 m/s.

圖13 X=3.75 m側送風速度場Fig.13 Velocity field under side air supplymode with X=3.75 m

圖14 X=3.75 m頂送風速度場Fig.14 Velocity field under top air supplymode with X=3.75 m

兩種不同送風方式下室內(nèi)溫度場、污染物(甲醛)濃度場及速度場的數(shù)值對比如表2所示.

表2 不同送風方式下不同參數(shù)值的對比Tab.2 Comparison in different parameter values under different air supply modes

3 結 論

本文利用計算流體力學技術模擬沈陽市某辦公類建筑中一會議室內(nèi)不同送風方式下的溫度場、污染物(甲醛)濃度場和速度場，并對模擬結果進行比較分析，得出以下結論：

1) 頂送風情況下所得室內(nèi)溫度場分布絕大部分處于23.7～25 ℃之間，滿足夏季室內(nèi)人員舒適度要求.但側送風方式會造成會議室內(nèi)西墻附近溫度偏低，局部區(qū)域內(nèi)溫度低于23 ℃，不利于會議室內(nèi)工作人員的思考和注意力的集中.

2) 會議室采取頂送風和側送風兩種送風方式都未出現(xiàn)會議室內(nèi)局部區(qū)域甲醛污染物濃度超出限值0.1 mg/m3的情況，在頂送風方式下Z=1.2 m截面上甲醛污染物濃度最高為1.055×10-7mg/m3，相比較而言頂送風方式下室內(nèi)甲醛濃度分布較為均勻，頂送風更有利于維持室內(nèi)良好的空氣品質(zhì).

3) 側送風與頂送風兩種送風方式下室內(nèi)風速均小于0.3 m/s，但頂送風方式下會議室內(nèi)空氣流速更趨于均勻，室內(nèi)風速基本處于0.013～0.093 m/s之間，會議室內(nèi)人員感覺更舒適.

綜上所述，會議室內(nèi)采用頂送風方式可更好地為室內(nèi)人員創(chuàng)造舒適良好的環(huán)境.

[1] 張寅平.室內(nèi)空氣質(zhì)量控制：新世紀的挑戰(zhàn)和暖通空調(diào)人的責任 [J].暖通空調(diào)，2013，43(12)：1-7.

(ZHANG Yin-ping.Indoor air quality control：the challenges and responsibilities of HVAC researchers in the new century [J].Journal of HV&AC，2013，43(12)：1-7.)

[2] 劉婷婷.某數(shù)據(jù)機房氣流組織模擬及運行優(yōu)化 [J].暖通空調(diào)，2015，45(3)：71-75.

(LIU Ting-ting.Air distribution simulation and operation optimization for a data room [J].Journal of HV&AC，2015，45(3)：71-75.)

[3] 王燁，張文霞，胡文婷.室內(nèi)環(huán)境參數(shù)對室外氣象參數(shù)瞬時變化的動態(tài)響應研究 [J].重慶大學學報，2015，38(3)：8-14.

(WANG Ye，ZHANG Wen-xia，HU Wen-ting.Dynamic response study of the indoor environmental parameters for variable outdoor meteorological conditions [J].Journal of Chongqing University，2015，38(3)：8-14.)

[4] 呂潔，張瓊，周勃，等.自然通風作用下計算機房夏季熱環(huán)境數(shù)值模擬 [J].沈陽工業(yè)大學學報，2016，38(1)：115-120.

(Lü Jie，ZHANG Qiong，ZHOU Bo，et al.Numerical simulation of summer thermal environment in compu-ter room under action of natural ventilation [J].Journal of Shenyang University of Technology，2016，38(1)：115-120.)

[5] 秦怡，李麗，周孝清，等.送風速度對層式通風房間空氣品質(zhì)影響的數(shù)值模擬 [J].制冷與空調(diào)，2016，30(3)：350-354.

(QIN Yi，LI Li，ZHOU Xiao-qing，et al.Numerical simulation on the influence of different supply air velocity on air quality under stratum ventilation [J].Refrigeration and Air-Conditioning，2016，30(3)：350-354.)

[6] 呂潔，吳亞平，周勃，等.低屋面橫向通風牛舍傾斜擋風板流場數(shù)值模擬 [J].沈陽工業(yè)大學學報，2015，37(6)：700-704.

(Lü Jie，WU Ya-ping，ZHOU Bo，et al.Numerical simulation for flow field of sloping wind deflector in low profile cross ventilated cow barn [J].Journal of Shenyang University of Technology，2015，37(6)：700-704.)

[7] 岳高偉，陸夢華，賈慧娜.室內(nèi)污染物擴散的通風優(yōu)化數(shù)值模擬 [J].流體機械，2014，42(4)：81-85.

(YUE Gao-wei，LU Meng-hua，JIA Hui-na.Numerical simulation of indoor pollutant diffusion on ventilation optimization [J].Fluid Machinery，2014，42(4)：81-85.)

[8] 駱偉，劉榮華，邱漢峰.自然通風作用下室內(nèi)甲醛濃度分布數(shù)值模擬分析 [J].制冷與空調(diào)，2013，27(5)：504-508.

(LUO Wei，LIU Rong-hua，QIU Han-feng.Numerical simulation and analysis of indoor formaldehyde concentration distribution under the action of nature ventilation [J].Refrigeration and Air-Conditioning，2013，27(5)：504-508.)

[9] 張金萍，于水靜，宋夢堃.家具和服裝市場室內(nèi)甲醛和PM2.5污染水平的測試研究 [J].建筑科學，2016，32(6)：21-26.

(ZHANG Jin-ping，YU Shui-jing，SONG Meng-kun.Study on the pollution level of formaldehyde and PM2.5 in indoor environment of furniture and clothing markets [J].Building Science，2016，32(6)：21-26.)

[10]周俊杰.FLUENT工程技術與實例分析 [M].北京：中國水利水電出版社，2010.

(ZHOU Jun-jie.FLUENT engineering technology and case analysis [M].Beijing：China Water & Power Press，2010.)

[11]王罡.室內(nèi)裝修空氣污染現(xiàn)狀與防治措施 [J].山西建筑，2014，40(26)：196-197.

(WANG Gang.The status and prevention measures of indoor decoration air pollution [J].Shanxi Architecture，2014，40(26)：196-197.)

[12]劉延秋，唐功臣，陶勇，等.某部營房室內(nèi)空氣中甲醛污染狀況調(diào)查 [J].解放軍預防醫(yī)學雜志，2008，26(1)：35-36.

(LIU Yan-qiu，TANG Gong-chen，TAO Yong，et al.Investigation on formaldehyde pollution in indoor air of a certain barracks [J].Journal of Preventive Medicine of Chinese People’s Liberation Army，2008，26(1)：35-36.)

[13]崔凱杰，古金霞，侯瑞，等.室內(nèi)空氣甲醛污染狀況及其影響因素分析 [J].南開大學學報(自然科學版)，2013，46(2)：28-31.

(CUI Kai-jie，GU Jin-xia，HOU Rui，et al.Indoor air pollution of formaldehyde and its affecting factors analysis [J].Acta Scientiarum Naturalium University Nankaiensis，2013，46(2)：28-31.)