陳勇 王軍 龍恩深四川大學建筑與環(huán)境學院

建筑室內(nèi)熱舒適性隨通風條件變化的敏感性研究

陳勇 王軍 龍恩深四川大學建筑與環(huán)境學院

針對典型建筑房間，通過建立物理模型與數(shù)學模型并結(jié)合CFD模擬方法，研究室內(nèi)熱舒適性隨混合通風、置換通風和地板孔板通風三種通風方式的送風條件變化的敏感性，得到人體前后側(cè)0.2m處空氣齡、溫度場、PMV、PPD分布及變化特征。研究結(jié)果將為不同通風方式下合理選擇送風條件以改善室內(nèi)人體熱舒適性提供理論指導。

混合通風 置換通風 地板孔板通風 熱舒適

隨著生活水平的提高，人們對室內(nèi)環(huán)境的關(guān)注不僅局限于最初的冷熱感覺，而更多地側(cè)重于人體的整體熱舒適性。針對室內(nèi)熱舒適問題，Salmanzadeh提出了影響室內(nèi)人體微環(huán)境熱羽流動的數(shù)值模型[1～2]；Mundt提出了置換通風時的溫度分布預(yù)測模型[3]；Rees提出了預(yù)測室內(nèi)特性參數(shù)節(jié)點模型[4～5]。將面熱源和體熱源虛擬為點熱源來簡化問題的描述，并獲得點源羽流流量和熱分層高度解的簡潔形式[6～7]。Kurazumi分析了人體周圍在不同環(huán)境條件下的對流換熱系數(shù)[8]。本文通過建立典型房間的物理模型并利用CFD方法對混合通風、置換通風、地板孔板通風在不同送風條件下室內(nèi)的人體微環(huán)境熱舒適狀況變化敏感性進行了分析，研究結(jié)果將為合理確定送風條件以提高室內(nèi)熱舒適性提供理論指導。

1 研究對象與模型

1.1 研究對象與物理模型

本文所研究的房間尺寸為2.5m（X）×2.3m（Y）× 1.5m（Z），如圖1所示?；旌贤L及置換通風的送風口和回風口尺寸均為0.3m×0.3m?？装逅惋L口為直徑0.01m圓孔，共330個。室內(nèi)人員為高強度工作狀態(tài)，人體發(fā)熱量為170W，分布在房間中央。房間壁面設(shè)定為絕熱條件。

圖1 物理模型

1.2 數(shù)學模型

建筑室內(nèi)空氣流動過程滿足以下假設(shè)：①室內(nèi)空氣為不可壓縮流體，即認為流體體積變化不足以影響計算結(jié)果的精度，此時體積的變化可以忽略不計；②室內(nèi)空氣視為連續(xù)介質(zhì)，即服從連續(xù)性定律，在單位時間內(nèi)各斷面上流過的空氣質(zhì)量不變；③空氣的所有熱物理性參數(shù)均視為常數(shù)，但密度隨溫度的變化遵循Boussinesq假設(shè)。對于室內(nèi)三維空氣穩(wěn)態(tài)流動過程，其控制方程為：

1）連續(xù)性方程

2）動量方程

3）能量方程

在上述室內(nèi)空氣過程的控制方程組中方程個數(shù)小于未知數(shù)個數(shù)，方程不封閉。因此，還需借助湍流模型建立湍流動力粘度μt與湍流時均參數(shù)之間的關(guān)系，從而使方程封閉。為了充分考慮人體熱浮力的影響，本文選擇RNG κ-ε模型使控制方程組封閉。

4）κ方程

5）ε方程

1.3 模型的數(shù)值求解

控制方程的求解首先采用基于交錯網(wǎng)格的控制容積法進行離散，離散過程中擴散項則采取中心差分格式，壓力項選用body-force-weighted格式，動量方程與能量方程的離散格式選用一階迎風格式，壓力速度耦合采用SIMPLE算法。同時，區(qū)域離散滿足網(wǎng)格獨立性條件。

2 結(jié)果分析與討論

2.1 空氣齡特征分析

本文主要針對混合通風（MV）、置換通風(DV)、地板孔板通風（BV）三種通風方式在不同送風速度和不同送風溫度條件下進行了模擬研究。不同送風速度下，混合通風、置換通風、地板孔板通風送風溫度分別為16℃、19℃、17℃。不同送風溫度下，混合通風、置換通風、地板孔板通風送風速度分別為1.5m/s、0.3m/s、1.2m/s。選取人體前后側(cè)0.2m處離地面0.2m至2.0m處空氣齡作為比較分析對象。不同送風速度（V）和送風溫度（T）的空氣齡分別如圖2、圖3所示。

圖2 三種送風方式在不同送風速度時空氣齡分布

圖3 三種送風方式在不同送風溫度時空氣齡分布

從圖中可以看出，混合通風條件下人體前后側(cè)空氣齡相差很小。由于混合通風誘發(fā)的強制對流主導整個流場，送風速度較大時，在人體呼吸區(qū)以上空間空氣齡比呼吸區(qū)域小。此時空氣齡不受送風溫度的影響而受送風速度的影響較大；同時，當送風速度增大相同數(shù)值時，空氣齡并不是同等變化，這說明在速度為1.0m/s左右存在一個臨界風速，當速度再增大時，對于空氣新鮮程度的改善作用很小。

另一方面，當為置換通風且在不同送風條件下時，人體前后側(cè)空氣齡相差很小。在垂直方向上，從人體足部到頭部的空氣齡逐漸變大，即置換通風方式下，人體活動區(qū)域比非活動區(qū)空氣新鮮。空氣齡隨送風速度的增大而減小。此外，送風溫度對人體0.5m以下區(qū)域空氣齡影響較大，主要是由于置換通風方式下送風速度較小，隨著送風溫度的增加，空氣浮升力增大，從而到達人體下部區(qū)域的新鮮空氣減少。

再次，在地板孔板通風方式下，送風速度和送風溫度對于空氣齡的影響都很小，相差10s以內(nèi)。由于送風引起的強制對流方向和人體熱羽流的自然對流方向相同，送風速度增大必然引起空氣齡的減小。同理，空氣溫度越高，熱浮力越強，則空氣上升越快。因此，地板孔板送風方式下，送風溫度越高，空氣齡越小。

2.2 溫度場特征分析

溫度是影響人體舒適度的一個重要參數(shù)，選取人體前后側(cè)0.2m處離地面0.2m～2.0m處溫度作為比較分析對象。三種不同通風方式在不同送風溫度和送風速度條件下的溫度分布情況分別如圖4、圖5所示。

圖4 三種送風方式在不同送風速度時溫度分布

圖5 三種送風方式在不同送風溫度時溫度分布

影響人體周圍溫度分布的主要因素為送風射流引起的強制對流和人體熱羽流的自然對流。從圖中可以看出，在混合通風條件下人體前后側(cè)溫度差別較小。在人體呼吸截面以下區(qū)域，溫度場的分布主要受到人體熱羽流的影響，由于熱羽流形成的自然對流不斷向上升，所以溫度隨高度的增加而不斷升高。而在呼吸截面以上區(qū)域，混合通風的強制對流越來越強，自然對流作用被削弱，所以人體周圍溫度峰值出現(xiàn)在呼吸截面。同時，送風速度越低，人體周圍溫度越高。當風速達到1.0m/s后，增大風速對溫度影響較小，說明在1.0m/s左右存在一個臨界風速。

其次，當在置換通風方式下時，人體熱羽流形成的自然對流不斷向上流動。同時，置換通風引發(fā)的強制對流空氣也隨高度的增加而溫度逐漸升高。因此，在置換通風方式下人體周圍空氣的溫度隨高度的增加而不斷升高。送風速度越低，人體周圍的空氣溫度越高，這是空氣對人體散熱吸收更加充分的結(jié)果。

此外，地板孔板送風方式與置換通風相似，人體熱羽流形成的自然對流不斷向上流動。同時，地板孔板通風引發(fā)的強制對流空氣也隨高度的增加而溫度逐漸升高。因此，在地板孔板通風方式下人體周圍空氣的溫度隨高度的增加而不斷升高。與其它兩種送風方式不同之處在于地板孔板送風引起的強制對流能夠迅速到達人體周圍，因此送風速度對于人體周圍溫度場分布的影響非常小。同時，可以看出地板孔板通風方式下，人體周圍垂直方向上溫度差很小。若單從溫度場分布來看，應(yīng)選擇地板孔板通風，溫差越小，人體感覺越舒適。

2.3 熱舒適指標特征分析

選取人體前后側(cè)0.2m處離地面0.2m至2.0m處的PMV和PPD分布作為比較分析對象。

2.3.1 PMV值分析

從圖6～7可以看出，混合通風方式下，人體周圍PMV值受溫度場分布影響很大。在人體軀干部分前后側(cè)周圍溫度相對較高，所以PMV值越大，人體感覺越熱。同樣可看出，風速在0.5m/s變化到1.0m/s過程中，PMV值變化較大。當風速達到1.0m/s后增大風速對PMV值的影響很小。所以風速對PMV值的影響在1.0m/s左右同樣存在一個臨界風速。同時，送風溫度越高，PMV值越大。

圖6 三種送風方式在不同送風速度時PMV分布

圖7 三種送風方式在不同送風溫度時PMV分布

其次，置換通風條件下的PMV分布情況與混合通風方式下的結(jié)果相似，都是人體軀干部分前后側(cè)周圍數(shù)值相對較高。送風速度越小，PMV值越大。同時，送風溫度越高，PMV值越大，人體感覺越熱。

再次，在地板孔板通風方式下，由于人體溫度分布不均勻，人體軀干部位溫度與周圍空氣溫差較小，造成人體軀干部位周圍比頭部和足部周圍PMV值小。PMV值隨送風速度的增加并不是單調(diào)的變化，送風速度為1.0m/s時PMV值最大。同樣，送風溫度越高，人體周圍PMV數(shù)值越大，人體感覺越暖和。

2.3.2 PPD值分析

從圖8～9可以看出，混合通風方式下，人體前后側(cè)PPD值相差較小。相同通風條件下，人體軀干部分周圍PPD值較小，不滿意率較低。同時，送風速度增大相同大小數(shù)值時，PPD值并不是同等變化，說明在1.0m/s左右存在一個臨界風速。送風溫度越小，PPD值越小，即人體不滿意率越小。故混合通風方式下應(yīng)選擇高溫、低速條件下進行送風。

圖8 三種送風方式在不同送風速度時PPD分布

圖9 三種送風方式在不同送風溫度時PPD分布

其次，置換通風方式下，送風溫度和送風速度對于PPD值的影響都非常小。同時，人體周圍的絕大部分PPD值都小于10%，并且波動范圍很小，也就是說大部分送風條件下人群中只有小于10%的人感覺不滿意，能滿足絕大部分人對熱舒適性的要求。

再次，地板孔板通風方式下，人體頭部和足部周圍區(qū)域比人體熱舒適性不滿意率低。送風速度對PPD值的影響較小，三種送風速度相比，當速度為1.0m/s時，人體周圍PPD值最小，人體熱舒適性不滿意率最低。同時，送風溫度越高，人體熱舒適性不滿意率越低。

3 結(jié)論

1）三種通風方式下人體前后側(cè)空氣齡差別很小?；旌贤L送風速度較大時人體活動區(qū)域以上部分空氣比活動區(qū)域空氣齡小。同時，在混合通風方式下，送風速度在1.0m/s左右存在一個臨界風速。置換通風方式下，人體0.5m以下區(qū)域空氣齡受送風溫度的影響較大。地板孔板通風方式下，人體空氣齡隨高度、送風速度、送風溫度的影響都很小。

2）混合通風、置換通風方式下人體微環(huán)境溫度場的分布受送風速度影響較大。混合通風方式下，受送風強制對流影響，溫度在呼吸截面出現(xiàn)拐點。同時，對溫度場分布的影響，送風速度在1.0m/s左右存在一個臨界風速。置換通風和地板孔板通風方式下，送風速度對溫度分布的影響較小。同時，人體周圍垂直溫度逐漸上升，但是置換通風溫差比地板孔板通風溫差大。

3）由于溫度分布、氣流組織的影響，混合通風、置換通風方式下PMV值從足部到頭部都是先增加后減少而地板孔板送風是先減小后增加?；旌贤LPPD值在不同送風溫度、速度下從足部到頭部都是先減小后增大，而地板孔板通風方式下PPD值先增大后減小。相同送風速度下，送風溫度越高，不滿意率越低。實際工程中可以適當提高送風溫度以降低不滿意率。

[1]M Salmanzadeh.Computational modeling of effects of thermal plume adjacent to the body on the indoor airflow and particle tran -port[J].Journal of Aerosol Science,2012,53:29-39

[2]Naiping Gao.CFD Study on micro-environment around human body and personalized ventilation[J].Building and Environment, 2004,39:795-805

[3]Elisabeth Mundt.The performance of displacement ventilation systems[J].Building and Environment,1996,10:23-40

[4]Hensen J L M,Hamelinck M J H.Energy simulation of displacem -ent ventilation in offices[EB/OL].[2005-10-04]

[5]Simon J Rees,Philip Haves.A nodel for displacement ventilation and chilled ceiling system in office space[J].Indoor Air,1996,5: 181-201

[6]Yuguo Li.Buoyancy-driven natural ventilation in a thermally strat -ified one-zone building[J].Building and Environment,2000,35: 207-214

[7]Mundt Z.Displacement ventilation systems-convection flows and temperature gradient[J]Building and Environment,1995,30:129-133

[8]Yoshihito Kurazumi,Lauris Rezgals.Convective heat transfer coe -fficients of the human body under forced convection from ceiling [Z].J Ergonomics issn:2165-7556.

Va ria tion Se ns itivity of Indoor The rm a l Com fort w ith Ve ntila tion Conditions

CHEN Yong,WANG Jun,LONG En-shen
College of Architecture and Environment,Sichuan University

As to typical building space,its physical model and mathematic model were established.Moreover, combined with CFD,variation sensitivity of indoor thermal comfort with ventilation conditions was determined for mixing ventilation,displacement ventilation and floor perforated panel ventilation.Meanwhile,distribution and variation characteristics of air of age,temperature,PMV and PPD at 0.2 meter front and back of human body were obtained.The results provides theoretical guidance for reasonably choosing air supply conditions to improving indoor human thermal comfort under different ventilation patterns.

mixing ventilation,displacement ventilation,floor perforated panel ventilation,thermal comfort

1003-0344（2015）03-005-4

2014-3-28

陳勇（1989～），男，碩士研究生；四川省成都市四川大學建筑與環(huán)境學院（610065）；E-mail:986179554@qq.com

國家自然科學基金資助項目（51308361）；四川省科技計劃項目（2014GZ0133）