045不同门窗开度下自然通风对室内空气品质及能耗的影响

门窗不同开度下自然通风

对室内空气品质及能耗的影响

清华大学建筑技术科学系石卉高鹏杨旭东

摘要一直以来，自然通风对于能耗的影响是大家普遍关注和研究的焦点。近年来，由于住宅装修污染问题的加重，人们对于室内空气品质问题的关注程度和要求也在日益地提高。自然通风对室内空气品质及住宅能耗的综合影响研究是很有意义的。本文将研究不同门窗开度情况下住宅的自然通风对室内空气品质和能耗影响。

首先以建筑热模拟模型（DeST）、多区域网络通风模型及多区域网络污染物浓度模拟模型为主要分析工具，搭建了一个分析自然通风、室内空气品质和能耗的综合平台；并利用计算流体力学软件研究并总结了平开类型开口不同开度下的流量系数变化规律；进而将上述成果应用于北京某住宅的实例中，通过实例的模拟，分析不同门窗开度组合下的自然通风状况以及其对室内空气品质及能耗的影响。

关键词开度、流量系数、自然通风、室内空气品质、能耗

1背景

人们关注自然通风对于能耗的影响是由来已久的事情, 而近年来，由于住宅装修污染问题的加重，人们对于室内空气品质问题的关注度和要求也在日益地提高。这种对能耗和空气品质的双重需求，最终决定着实际生活中的如何开关门窗的问题。而很多时候，这两种需求对自然通风的期望并不一致，比如在炎热的夏季和寒冷的冬季，自然通风在降低室内能耗方面是应当尽量避免，而在保证室内空气品质方面是应当大力提倡的。所以，不同门窗开启状态下自然通风对室内空气品质及住宅能耗的综合影响研究是一个十分有意义和重要的课题。

要回答自然通风对二者的综合影响这一问题，准确地确定不同门窗开启状态下的自然通风量显得尤为重要。模拟计算是人们研究自然通风的有利工具，目前常见的自然通风软件有基于N-S方程的计算流体力学相关软件，比如Fluent，PHOENICS等，其特点是能详细的刻画房间中流场，有较高的精度，而缺点是模拟时间长，代价高，不太适合于长期动态模拟问题；另外还有基于多区域网络模型的CONTAMW，COMIS等。这些软件在对多房间的自然通风进行长期动态模拟方面见长，但没有考虑与热的耦合，输入参数较多，且对最后模拟结果的准确性影响较大[1]。所以总的来说，在自然通风量的软件方面，还需要人们根据需求进一步完善和提高。

本文首先利用现已成熟的建筑热模拟软件DeST和多区域网络模型（通风模型，污染物传播模型）

搭建综合分析室内自然通风、能耗、空气品质的模拟平台；针对多区域网络模型输入参数对模拟计算结果准确性影响较大的情况，利用数值计算流体力学相关软件，对建筑中普遍存在的平开类型开口阻力特性规律进行模拟分析。最后以北京一户住宅为例，模拟不同门窗开度组合的情况下的室内空气流通情况，分析其对空气品质和能耗的影响，并从中得到一定结论。

2 模拟平台简介

此综合模拟平台由两部分组成：DeST-vent 和IAQ 模型。

2.1 DeST-vent

DeST 是由清华大学建筑技术科学系开发的建筑热模拟软件，目前在中国已经得到广泛应用。其最近开发的DeST-vent 通风版本中。采用onion 耦合的方法[2,3]综合考虑了通风与热之间的耦合作用：通风模型将计算结果输入热模拟模型，后者再将计算结果反馈回前者，如此循环，直到前后两次计算结果之差满足精度要求，再进入下一时刻的计算。

DeST-vent 通风计算采用多区域网络模型，将每个房间视为一个节点，同一节点内空气混合均匀，空气温度、密度、压力和污染物浓度等参数均匀一致；门、窗、缝隙等视为支路。所有的节点和支路组成流体网络。这样对于有N+1个节点（其中1个节点代表室外），B 条支路的网络，每个节点流量平衡，每个支路均满足伯努利方程，这样可以得到N+B 个方程来求解N+B 个未知数。具体方程见下：

视为常温常压下空气，根据理想气体状态方程：

i i i i i i RT V

P V m ==ρ

(1) 式中： i P 为节点参考压力； i V 为节点体积；

i T 为节点空气温度，采用绝对温标；

R 为理想气体常数[J/kg K] 对于常温常压下的干空气，

气体常数R air ＝287 J/kg K 。

由于建筑通风都在中低压下进行的，空气流动可以考虑为不可压缩流体的连续流动，所以对于通风网络中的任意一个节点i ，根据连续性方程都有：

i j

i

j F f

=∑- (2)

其中i j f -为与该节点相连的各支路质量流量，i F 为该节点的净流量，忽略动能的变化。对自然通风网络中任意一条支路e ：

当e 为内部支路时，根据恒定总流能量方程有：

gz f f s gz P gz P P i j j i j i e j j i i e ij )()()(ρρρρ-+=---=∆-- (3)

当e 为外部开口，i 为室外时：

图 1 通风支路示意图

2

2

)()()(wind

wind

e

p

i j j i j i e j j e w i i e ij v C gz f f s gz P P gz P P ρρρρρ+-+=--+-=∆-- (4) 其中wind v 为室外风速；e

w P 为室外风场引起的支路e 压力增加值；e

p C 为支路e 建筑表面风压系数；

z 为开口e 中心距离房间地面的高度；e

S 为支路e 的阻力系数。

应用上式，以矩阵和向量的形式表达，考虑风压作用以及支路两侧节点位能差，建筑自然通风网络模型方程如下：

⎪⎪

⎩

⎪⎪⎨⎧-+=∆∆==⎰wind T

P gdZ F F S P P

P A AF ρ0 (5) 式中参数定义如下：A － 网络关联矩阵；F － 支路流体质量流量，kg/s ； P － 网络节点静压，Pa ；P ∆－ 支路两端压力差，Pa ；wind P －室外来流对支路产生的压力增加值，Pa ；S － 支路阻力系数，22

kg s

Pa ⋅；Z － 支路垂直投影高度，m ；ρ－ 空气密度，随温度不同而不同，kg/m 3；g － 自

由落体加速度，2

s m 。

整个系统有N 个未知数，N 个方程，求解此方程组，即可获得各个支路的流量，同时可以得到房

间节点、通风网络各个节点的静压。

2.2 IAQ 模型

以DeST-vent 中计算得到的空气流量为输入参数，可以计算得到室内的污染物浓度。污染物输运方程、污染物散发方程及其边界条件如下：

)()()()(t R t S t C Q t C Q dt dC V i i j

i j i j j i j i

i

-+-=∑∑-- (6)

⎪⎪

⎭⎫ ⎝⎛∂∂∂∂

=∂∂j m m j

m x C D x t C (7)

s

m KC C = (材料-空气界面的边界条件) (8)

()room s m m

m

C C h x C

D -=∂∂- (空气边界条件) (9)

C C m = at t=0 (材料中污染物初始浓度) (10)

t ——时间；C i ——节点i 中污染物浓度； V i ——节点i 的体积；Q j-i ——从节点j 到i 的体积流量；

S i ——节点i 中污染源的散发效率；R i ——节点 i 中汇的吸附效率。C m ——材料中VOC 浓度； C 0 —