【论文精选】基于场-网络耦合模型建筑内区会议室CO₂分析

官网注册后

http://www.gasheat.cn

免费下载论文

http://www.gasheat.cn/Periodical/index.html

作 者： 王雷，王砚玲，王芳，王海燕

第一作者单位： 哈尔滨工业大学建筑学院

摘自《煤气与热力》2023年9月刊

参考文献示例

王雷，王砚玲，王芳，等 .    基于场 - 网络耦合模型建筑内区会议室 CO ₂ 分析 ［J］.    煤气与热力， 2023，43（9） ：A27-A33.

1    概述

目前，国内外现行室内空气质量标准将 CO ₂ 作为重点控制污染物，其含量已成为评价室内空气质量的重要指标 ^{［ 1 ］} 。 Shendell 等人通过研究指出教室空气中 CO ₂ 等含量过高对学生完成作业效果、学习效率有明显的抑制作用，室内 CO ₂ 体积分数超过标准限值会导致学生注意力下降 ^{［ 2 ］} 。现有 CO ₂ 研究中主要针对某种类型建筑中特定房间的 CO ₂ 体积分数变化，没有考虑建筑中其他参数变化对 CO ₂ 体积分数变化的影响。

场模型是分析室内空气流动对污染物扩散作用的常用方法，对复杂建筑整体进行场模型计算分析，存在网格划分复杂、计算时间长的限制，场模型并不适用。较合理的研究方法是采用多区域网络模型 CONTAM ^{［ 3 ］} ，并对特定区域进行 CFD 建模并分析。为此， Wang 等人 ^{［ 4 ］} 用 CFD0-CONTAM 耦合模型（简称耦合模型）对建筑中空气流动进行了分析，其中 CFD0 是一种 CFD 程序，最初由 Srebric 等人 ^{［ 5 ］} 开发，由 Wang ^{［ 6 ］} 改进。 Barbosa 等 人 ^{［ 7 ］} 应用耦合模型在实验室和医院进行了污染物扩散分析。张旭涛等人 ^{［ 8 ］} 应用耦合模型对高层建筑火灾中电梯竖井烟气流动特性进行了研究。文献［ 7-8 ］都体现出耦合模型具有一定可靠度和优越性。

本研究中考虑复杂建筑中空气流动对某内区会议室 CO ₂ 体积分数变化的影响情况，该会议室无与室外相连的外窗。此会议室会存在短时间内大量人员聚集导致室内人员因高 CO ₂ 体积分数而效率降低的情况，通风效果较差，人员长期处于高 CO ₂ 体积分数的场所，健康会遭受累积性的损害，故有必要对此问题进行研究。因此，本文采用 CFD0-CONTAM 耦合模型（场 - 网络耦合模型）对哈尔滨某综合建筑进行空气流动分析，并以建筑 3 层某内区会议室为例，对空气流动和 CO ₂ 体积分数分布进行分析。

2    建筑模型

本研究中的综合办公楼位于哈尔滨市，共 5 层，总建筑面积为 10 600 m ² ，层高为 3 m ，包含一个 3 层高的中庭。运用 CONTAM 对整体建筑进行建模，建筑内各房间门窗等开关状态均考虑，建筑第 3 层平面见图 1 ，呈“ T ”形。本文研究关注的对象是位于建筑第 3 层的会议室，会议室建模见图 2 ，该会议室紧邻中庭，无与室外相通的外窗，会议室不使用时门窗均保持关闭状态，并考虑建筑物其他房间门窗对研究对象的影响。会议室体积为 90 m ³ ，长、宽分别为 10 m 、 3 m ，高 3 m ，会议室门窗各 2 个。门的中心点坐标分别为（ 1.25 m ， 0 m ， 1 m ）和（ 8.75 m ， 0 m ， 1 m ）；窗的中心点坐标分别为（ 1.25 m ， 3 m ， 1.5 m ）和（ 8.75 m ， 3 m ， 1.5 m ）。

图 1    建筑第 3 层平面（软件截图）

图 2    会议室建模

门窗采用了 CONTAM 内的泄漏面积模型和孔口模型，楼梯等选取相应模型。由于我国关于墙体泄漏特性的数据较少，而且建筑墙体具有较强的密封性，故本研究忽略墙体泄漏。参照相关文献和实测研究，空气流通路径泄漏特征值见表 1 。门窗开启时，有效泄漏面积等于门窗的实际面积。整个建筑物中的窗户尺寸均相同。

表 1    空气流通路径泄漏特征值

内区会议室不与室外相通，易存在开会时室内 CO ₂ 体积分数超过标准限值（ 1 000 × 10 ^-6 ）的情况，耦合模型将其定义为 CFD 区域，采用湍流零方程模型进行计算，网格数为 9 × 10 ⁴ ，网格已经进行无关性验证，无需进一步细化。内区空气流通路径采用的模型设置同建筑物其他房间。

本研究选择 CONTAM 为 CFD0 提供压力边界条件的耦合方法，因为该耦合方法提供了稳定收敛的耦合解 ^{［ 7 ］} 。选择均质多组分流动模型进行污染物扩散预测 ^{［ 7 ］} 。在连接 CFD0 和 CONTAM 中气流路径形成的边界上进行信息传递，信息传递内容包括流量、压差和污染物体积分数， Wang 等人 ^{［ 9 ］} 详细描述了耦合过程。 CFD0-CONTAM 耦合收敛因子设置为 0.001 。

3   CFD0-CONTAM 耦合下会议室 CO ₂ 分析

3.1   会议室空气流动及 CO ₂ 体积分数模拟分析

8 ： 00 — 9 ： 00 在会议室举行时间为 1 h 的会议，共 20 人参会，会议期间门窗保持关闭状态，在会议结束后，人员离开，将会议室门窗开启 30 min ，至此本次研究结束。对于每个参会者，按实际会议室的位置进行建模，将人理想化成为一个正六面体，即散发 CO ₂ 的污染源，六面体边长均为 0.1 m 。图 2 中参会者（正六面体）的中心均在 z=1.05 m 平面上，构成一个长方形。长方形的长边上，相邻参会者中心距 0.6 m ，长方形的短边上，相邻参会者中心距 0.55 m 。 CO ₂ 的产生速率为 0.005 2 L/s ^{［ 10 ］} ，由六面体的 6 个面均匀向外释放。模拟时考虑了热压和风压的共同作用。在过渡季，室外温度为 3 ℃，风向为 SE ，风速为 3 m/s ，大气压为 101 325 Pa ，会议室室内温度为 18 ℃。初始室外环境 CO ₂   体积分数、室内 CO ₂   体积分数取 400 × 10 ^-6 。室内外温度在模拟期间保持不变。

过渡季会议室空气速度矢量图见图 3 ，图中红点为坐标原点。在过渡季，在 x=8.55 m 截面（该截面体现会议室的门和窗），室内空气流动的方向从窗户流向门；在 y=2.05 m 截面（该截面体现 CO ₂ 污染源），在参会者附近，空气流速小；在 z=1.05 m 截面（该截面也体现 CO ₂ 污染源），可以看到 2 个旋涡， x=0~5 m ，空气旋涡呈逆时针， x=5~10 m ，空气旋涡呈顺时针，空气流动呈对称分布，是由于 2 个门窗物理参数和压强相同导致的。

图 3    过渡季会议室空气速度矢量图（软件截图）

3.2   实测与模拟对比

为了验证 CFD0-CONTAM 耦合模拟的有效性，对会议室在上述模拟工况进行了现场实测，实测期间 20 个参会者，会议期间人员静坐，有小声交谈情况。选用的实测仪器为室内温湿度及 CO ₂ 体积分数监测仪。

对于 CO ₂ 的实测点布置，结合本地区实际因素和实验条件，保证数据准确且具有代表性的条件下，选取 3 个实测点。实测点 A 位于房间水平面中心桌子上，距地面 1.05 m ，实测点 B 位于会议室墙角处，距地面 1.05 m 。实测点 A 和实测点 B 用于采集室内 CO ₂ 体积分数的变化情况。另外一个仪器布置在外窗窗沿上，监测室外 CO ₂ 体积分数和温度，根据此仪器监测数值定义模拟分析中室外环境 CO ₂ 体积分数，此数值在监测期间变化很小，室外环境 CO ₂ 体积分数取 400 × 10 ^-6 。

会议室 CO ₂ 体积分数的变化见图 4 。图 4 中所显示的模拟值为 CFD0-CONTAM 耦合下室内平均 CO ₂ 体积分数。从图 4 可以看出，在进行会议期间，会议室内 CO ₂ 体积分数实测值逐渐升高，在 8 ： 40 时，室内 CO ₂ 体积分数实测值可以达到 1 400 × 10 ^-6 ~1 500 × 10 ^-6 ，这时参会者的注意力容易下降，到 9 ： 00 时达到峰值 1 817 × 10 ^-6 ，超过标准限值 817 × 10 ^-6 ，这时参会者的工作效率将显著下降并伴有困倦状态。在会议结束后将门窗打开， 30 min 后室内 CO ₂ 体积分数实测值降到接近环境体积分数。

图 4    会议室 CO ₂ 体积分数变化

由图 4 发现，实测点 A 、实测点 B 的 CO ₂ 体积分数与模拟值的变化趋势相同，由此说明， CFD0-CONTAM 耦合模拟的结果具有较高的准确性。实测点 A 的 CO ₂ 体积分数高于实测点 B ，是由于参会者坐的位置较集中导致。模拟值低于实测值，考虑可能是人员 CO ₂ 产生速率设置有偏差，且模拟值为空间内平均值等原因。

3.3   不同季节下门窗开关状态对会议室 CO ₂ 体积分数分布影响

探究不同季节和门窗开关状态条件下会议室 CO ₂ 体积分数分布情况。会议室内无机械通风，在开会期间开启门窗更符合实际。过渡季气象参数同前，夏季室外温度为 35 ℃，会议室温度为 25 ℃；冬季室外温度为 -30 ℃，室内温度同过渡季一样，为 18 ℃。室内外温度在模拟期间保持不变。所有季节风速都设为 3 m/s ，风向为 SE 。各工况基本信息见表 2 。

表 2    模拟工况条件

工况 1 若干典型时刻会议室 CO ₂ 体积分数分布见图 5 ，在工况 1 下，过渡季会议室门窗全开启时， 20 名参会者在会议室进行 1 h 会议。从 z=1.05 m 截面 CO ₂ 体积分数分布看，在 8 ： 10 时刻，参会者周围 CO ₂ 体积分数较高，室内局部 CO ₂ 体积分数达到近 3 500 × 10 ^-6 ，此时参会者会明显感觉疲惫困倦，严重影响工作效率。基于过渡季房间空气流动方向从窗户到门，故靠近门的位置 CO ₂ 体积分数高于靠窗侧。在 9 ： 00 会议结束时，参会人员周围 CO ₂ 较均匀分布。在会议结束 10 min 后，室内 CO ₂ 体积分数普遍下降至 1 000 × 10-6 以下，会议结束 30 min 后，室内 CO ₂ 体积分数已降至环境 CO ₂ 体积分数。从 y=2.05 m 截面 CO ₂ 体积分数分布看，受室内空气流动不均匀影响，在 8 ： 10 时，人员附近 CO ₂ 体积分数较高，且室内下部区域 CO ₂ 体积分数高于上部区域。随着会议进行，室内下部空间 CO ₂ 体积分数逐渐变高，可达到 2 000 × 10 ^-6 。在 9 ： 30 时，室内 CO ₂ 已扩散至走廊和其他区域，室内 CO ₂ 体积分数降至环境体积分数。

图 5    工况 1 若干典型时刻会议室 CO ₂ 体积分数分布（软件截图）

工况 2 、 3 、 4 若干典型时刻会议室 CO ₂ 体积分数分布见图 6 。工况 2 下会议室门窗均处于关闭状态，室内空气流速较低。会议期间，在 z=1.05 m 截面上，工况 2 同工况 1 下 CO ₂ 的扩散趋势相似，在参会者周围 CO ₂ 体积分数较高。但与工况 1 下 CO ₂ 扩散范围不同，会议室门窗关闭状态下 CO ₂ 体积分数扩散范围更广。

分析过渡季（工况 1 ）、夏季（工况 3 ）、冬季（工况 4 ） 3 个不同气候条件下会议室 CO ₂ 体积分数分布。 3 个工况下会议室门窗均处于开启状态。通过空气流动分析可知，在冬季和过渡季时空气流动方向是从窗户到门，夏季的空气流动方向则是从门向窗。夏季时， z=1.05 m 截面 CO ₂ 体积分数显著高于过渡季，且夏季 CO ₂ 体积分数衰减时间长于过渡季，故在夏季时应该重点关注室内的通风问题，建议采取机械通风措施。冬季时，由于热压作用较强，在 z=1.05 m 截面， CO ₂ 体积分数均在 1 200 × 10 ^-6 以下，会议结束后， CO ₂ 体积分数衰减较快，在会议结束 10 min 后可降至环境体积分数。

图 6    工况 2 、 3 、 4 若干典型时刻会议室 CO ₂ 体积分数分布（软件截图）

3.4   门窗开启比例对会议室 CO ₂ 体积分数分布影响

从上述分析可知，室外气象条件是影响室内 CO ₂ 体积分数的重要因素，处于夏季的工况 3 较工况 1 和工况 4 ，室内 CO ₂ 体积分数较高，且衰减时间较长。建筑物第 3 层工作人员较多，门窗开启频繁，因此下文针对过渡季和夏季若干工况进行研究，探究建筑第 3 层门窗开启比例不同对会议室 CO ₂ 体积分数分布的影响。结合哈尔滨地区气象条件及使用习惯，仅考虑建筑第 3 层门窗开启比例为 50% 或 80% 。

会议室门窗空气流通路径参数不变，将建筑第 3 层所有房间的门窗数量进行统计，共 110 个外窗（与室外空气连通）， 4 个内窗（不与室外空气连通）和 96 个门。按照仅开窗、仅开门、既开门又开窗 3 种方式对建筑第 3 层所有房间的空气流通路径进行整体调节，门窗开启的部位由建筑第 3 层“ T ”形的 3 个边按比例进行选取。会议室内仍进行 1 h 会议，选取的工况见表 2 。

工况 5~8 若干典型时刻会议室 CO ₂ 体积分数分布见图 7 。工况 5 和工况 6 在会议期间 CO ₂ 体积分数变化趋势与工况 1 相似。从 z=1.05 m 截面的模拟结果看，工况 5 在 9 ： 10 时， CO ₂ 体积分数衰减较工况 1 更慢，在会议结束 30 min 后也未恢复至环境体积分数。工况 6 在 9 ： 10 时， CO ₂ 体积分数比工况 5 稍降低，但区别不明显。

图 7    工况 5~8 若干典型时刻会议室 CO ₂ 体积分数分布（软件截图）

在建筑第 3 层所有房间门均关闭的情况下，随着建筑第 3 层窗开启比例增大，从 z=1.05 m 截面看，会议室 CO ₂ 体积分数逐渐变小且衰减速度较快，在 9 ： 30 时，工况 8 条件下会议室 CO ₂ 体积分数可恢复至环境体积分数。

工况 9~12 若干典型时刻会议室 CO ₂ 体积分数分布见图 8 。由图 7 和图 8 可知，建筑第 3 层门窗开启比例均为 50% 和均为 80% 时，从 z=1.05 m 截面看， 9 ： 00 时室内高 CO ₂ 体积分数的范围比工况 1 更小。根据工况 5~10 的规律，建筑第 3 层除会议室外的其余房间门窗开启会影响会议室 CO ₂ 体积分数的变化，其余房间门窗同时开启对会议室 CO ₂ 体积分数的影响最大。因此，再对工况 3 进行分析。

在 z=1.05 m 截面，工况 3 与工况 11 和工况 12 相比， CO ₂ 体积分数较高。根据模拟数据，工况 3 实测点 A CO ₂ 体积分数约为 1 200 × 10 ^-6 ，实测点 B CO ₂ 体积分数约为 800 × 10 ^-6 。工况 11 和工况 12 在 z=1.05 m 截面上 CO ₂ 体积分数低于 1 200 × 10 ^-6 的范围更大。从工况 9 、 10 来看，当建筑第 3 层门窗开启比例为 80% 时，室内高 CO ₂ 体积分数的范围比建筑第 3 层门窗开启比例为 50% 更小。但对于夏季，空气流动方向与过渡季不同，是从室外流向室内，因此工况 11 、 12 相比，室内 CO ₂ 体积分数变化幅度不明显。

图 8 工况 9~12 若干典型时刻会议室 CO ₂ 体积分数分布（软件截图）

选取会议室内 D 点（ 2.4 m ， 1 m ， 1.1 m ）作为分析点，该点处于参会者周围， CO ₂ 体积分数变化较明显。不同工况典型时刻下 D 点 CO ₂ 体积分数变化见图 9 。可以发现在会议结束 10 min 内，仅开门工况（工况 5 、 6 ）变化速率小于仅开窗工况（工况 7 、 8 ），说明了仅将建筑第 3 层某些门打开对会议室内 CO ₂ 体积分数影响效果小于某些开窗工况。在过渡季开门开窗（工况 9 、 10 ）下，局部室内 CO ₂ 体积分数低于其余工况，说明本建筑过渡季同层其余房间开门开窗可以保证内区会议室内 CO ₂ 体积分数处于较低范围。但夏季开门开窗时（工况 11 、 12 ）， CO ₂ 体积分数仍然较高，仅依靠自然通风不能使 CO ₂ 体积分数降到较低水平，必要时应采取机械通风措施，才能够保证室内人员的良好工作状态。

图 9    不同工况典型时刻下 D 点 CO ₂ 体积分数变化

4    结论

①模拟结果与实测结果具有一致性， CFD0-CONTAM 耦合模型可靠。如提高人员 CO ₂ 产生速率、空气泄漏特征的准确性，那么耦合模型准确性将进一步提高。

②夏季建筑外窗关闭条件下，建筑物整体空气流通较弱，内区会议室容易出现 CO ₂ 体积分数高的情况。因此，在必要时，应该对建筑内区的会议室采取机械通风的措施。

③不同季节下，建筑物热压、风压作用特征不同，对会议室 CO ₂ 体积分数水平影响较大。建筑物其他房间门窗关闭情况下，冬季和过渡季相比，会议室 CO ₂ 体积分数低且衰减快。

④仅将建筑第 3 层窗开启对会议室 CO ₂ 体积分数的影响大于仅将建筑第 3 层门开启。故建筑外围护结构上空气流通面积比内围护结构上空气流通面积对会议室空气流通影响大，进而对会议室 CO ₂ 体积分数影响较大。

参考文献：

［   1  ］ ZHANG X J ， LI D D ， XIE J C ， et al. Environmental Perceptions ， Mental Performance ， and Physiological Responses of People with Respiratory Allergies Exposed to Reduced Indoor Air Quality ［ J ］ . Indoor Air ， 2021 （ 5 ） ： 1458-1472.

［   2  ］ SHENDELL D G ， PRILL R ， FISK W J ， et al. Associations Between Classroom CO ₂ Concentrations and Student Attendance in Washington and Idaho ［ J ］ . Indoor Air ， 2010 （ 5 ） ： 333-341.

［   3  ］ WANG L Z ， STUARTD W ， CHEN Q Y. Using CFD Capabilities of CONTAM 3.0 for Simulating Airflow and Contaminant Transport in and Around Buildings ［ J ］ . HVAC & R Research ， 2010 （ 6 ） ： 749-763.

［   4  ］ WANG L Z ， CHEN Q Y. Validation of a Coupled Multizone-CFD Program for Building Airflow and Contaminant Transport Simulations ［ J ］ . HVAC & R Research ， 2007 （ 2 ） ： 267-281.

［   5  ］ SREBRIC J ， CHEN Q Y ， GLICKSMANL R. Validation of a Zero-Equation Turbulence Model for Complex Indoor Airflow Simulation ［ J ］ . ASHRAE Transactions ，  1999 （ 2 ） ： 414-427.

［   6  ］ WANG L Z. Coupling of Multizone and CFD Programs for Building Airflow and Contaminant Transport Simulations （ Doctoral Dissertation ） ［ D ］ . West Lafayette ： Purdue University ， 2007 ： 26-29.

［   7  ］ BARBOSA B ， BRUM N. Validation and Assessment of the CFD-0 Module of CONTAM Software for Airborne Contaminant Transport Simulation in Laboratory and Hospital Applications ［ J ］ . Building & Environment ， 2018 ， 142 ： 139-152.

［   8  ］张旭涛，王松岭，赵加宁 .  高层建筑火灾中电梯竖井烟气流动特性的模拟研究［ J ］ .  建筑科学， 2014 （ 4 ） ： 92-97.

［   9  ］ WANG L Z ， CHEN Q Y. Theoretical and Numerical Studies of Coupling Multizone and CFD Models for Building Air Distribution Simulations ［ J ］ . Indoor Air ， 2010 （ 5 ） ： 348-361.

［ 10 ］ PERSILY A. Ventilation ， Carbon Dioxide and ASHRAE Standard 62-1989 ［ J ］ . ASHRAE Journal ， 1993 （ 7 ） ： 42-44.

（本文责任编辑：李欣雨）

同类论文

1、 新风系统对室内CO₂质量分数分布影响的模拟

2、 自然通风下教室CO₂体积分数的实测与模拟

维普免费下载《煤气与热力》论文（现刊和过刊均可）

日前，《煤气与热力》杂志社有限公司在维普网站 http://cqvip.com/ 开通论文免费下载服务，论文刊出后两个月后，可在维普网站查询，并直接免费下载。在维普网站 免费下载《煤气与热力》论文 步骤如下：

1. 在维普网站注册会员。

2. 搜索出《煤气与热力》论文，点击进入。

3. 论文免费下载界面截图见上图。点击“免费下载”，可直接下载该论文。

声明：本文著作权（版权）归《煤气与热力》杂志社所有，严禁任何微信号及媒体未经授权许可随意转载。 PS: 当然欢迎大家转发到朋友圈！

更多论文请登录煤气与热力杂志官方网站，免费注册会员阅读电子期刊。阅读步骤：登录http://www.gasheat.cn/→页面右上角注册会员→注册成功后点击《煤气与热力》→期刊索引→点击某期期刊封面即可阅读当期文章。