标准办公室热环境数值模拟

标准办公室热环境的数值模拟

摘要：办公室是目前工作人员停留时间较长的场所，文章采用cfd技术，对标准办公室在不同气流组织与空调工况下的热环境进行数值模拟,以便为提高办公室的热舒适性设计提供参考依据。

关键词：办公室；cfd；数值模拟

中图分类号： c931.4 文献标识码：a 文章编号：

numerical simulation of thermal environment on standard office

li xiaofei, cheng haifeng

anhui architecture design institute of anhui university of architecture and industry

abstract: office as a place for missionary staying, in order to supply some references in improving thermal comfort of office place, simulated the standard office under different air flow and air conditioning circumstances by using cfd technology.

key word: office; cfd; numerical simulation

0 前言

随着办公建筑的大量兴建和人们对室内空气品质要求的提高，办公室室内空气的热舒适性变的更加重要。作为工作人员长期停留的公共场所，从暖通空调专业来说，应保持健康舒适的室内热环境。随着计算机技术的发展，空调领域中利用cfd技术来模拟气流组织，

热舒适性以及优化系统设计的应用越来越广泛。本文采用cfd方法，对标准办公室不同气流组织和冬夏季空调工况进行了数值模拟，为以后办公建筑的空调设计提供参考依据。

1 数学模型

1.1 室内零方程模型

除了遵循三大守恒定律以外，本文的数学模型采用室内零方程喘流模型（indoor zero equation）。室内零方程模型是在室内空气自然对流和混合对流的直接数值模拟dns (directly numerical simulation)结果的基础上提出的湍流模型[1]，该模型针对房间内非等温流rayleigh 数范围（2.6～3.0×1010），认为涡粘系数正比于流体密度、当地速度和距壁面最近之距离，比例系数由直接数值模拟的结果拟合而得

(1)

其中：为室内空气密度；为当地平均速度；l为长度标尺(取离壁面最近的距离)，0.03784为经验值。该模型少求解2个微分方程，而仅求解关于质量、动量和能量守恒的5个微分方程，故计算最省时间。

1.2 boussinesq密度假设

浮升力模型采用非线性boussinesq假设[2]：

(2)

2 几何模型及边界条件

2.1 模拟对象及网格划分

模拟对象为某标准办公室，办公室尺寸长、宽、高分别为5m、4m、3m，模拟条件为室内四个工作人员，坐姿模式，高度为1.1m；四台电脑，三盏灯具，送风口尺寸为0.8m×0.2m，回风口尺寸为0.9m×0.3m，南外墙上有外窗。具体几何模型见图1所示，模型网格划分采用六面体结构化网格，在各个热源以及送、回风口处局部加密，见图2所示。

图1 模拟几何模型

图2 网格划分

2.2 边界条件的简化

为方便计算，作如下简化和假设[3]：

(1)室内空气流速低，可视为不可压缩流体且符合boussinesq 密度假设；

(2)四面外墙仅南外有温差传热，设定为第一类热边界条件；南外窗同设为第一类边界条件，给定温度值，其余维护结构设为绝热；

(3)人体、电脑和灯具采用第二类热边界条件，给定热通量，人员在办公环境下属于极轻活动，全热散热量取134w[4]，灯具取45w/台，电脑150w/台；

(4)不考虑门空气渗透的影响；

(5)考虑重力因素的影响，沿z向设定重力加速度的大小和方向[5]。

3 模拟结果与分析

3.1 夏季工况模拟结果

1.假定室外气温为36℃，气流组织为侧送侧回，送风口风速

2.0m/s，送风温度为15℃（夏季工况1），取过房间中心线y=2m和人体脚踝以及站姿z=0.1m（脚踝面）、z=1.7m（站姿面）剖面来分析，模拟结果见图3所示，根据图3分析可知：

（1）由图3可以看出：整个流场属于典型的侧送风贴附式射流，空气由送风口送出，起始端风速较大，沿射流方向下沉并且温度和速度缓慢衰减，在送风达到南外墙时折回经回风口排出，流场上部为紊流区，下部为单向流动区。热源附近区域温度明显要高且对室内流场有显著的扰动作用，房间中心过道由于扰动较小，温度和速度均低于房间其他区域；室内存在温度梯度但分层现象不明显，总体来说，离送风口和回风口区域温度较低，外墙和外窗对温度场和速度场有明显的“排斥”作用。

（2）在假定的送风条件下，温度场在房间内变化较大，在整个人员工作区域（1.7m以下位置）温度处于25℃以下，脚踝处温度在21℃左右，地板到站姿面之间的温度梯度在3℃以下，人体的热感觉差异不明显；同时，速度场在房间的区域变化同样显著，过道、房间上部区域、热源、送回风口处温度值较高，脚踝面除过道处处在0.3m/s以下，站姿模式下除过道处都处在0.4m/s以下。综合温度和速度分布，人体的热感觉较为舒适，符合现行规范舒适性空调活动区风速≤0.3m/s的规定[6]。

(a) y=2m室内速度场分布

(b) y=2m室内温度场分布

(c) z=1.7m室内速度场分布

(d) z=1.7m室内温度场分布

(e) z=0.1m室内速度场分布

(f)z=0.1m室内温度场分布

图3 夏季工况1模拟结果图

2.调换送回风口的位置，设定送风温度18℃，送风风速1.5m/s （夏季工况2），模拟结果见图4所示，根据图4分析可知：（1）由图4得到：竖直方向上温度存在分层，下部区域温度低上部温度高；温度较高区域分布在南外墙和窗、热源、天花板处，温度低的区域处在房间的下部、过道处；热源和维护结构的传热对室内流场扰动较大；房间的热量被带到上部区域经回风口排出，有类似置换通风的作用；房间下部区域为速度场的单向流动区，送回风口、热源附近、过道的速度值较大。

（2）温度场竖直方向上存在热力分层，脚踝面在热源附近温度较高24℃左右，地板面到脚踝面的温度梯度很小，并且脚踝面0.1m 到人体坐姿高度1.1m的温差小于3k，满足国际标准的推荐值；其余区域温度在20.5℃以下；工作区域（1.7m以下位置）18℃-22℃之间。从速度场显示可以看出，热源对速度场的影响比上侧送的效果更为显著；脚踝面的速度值不高于0.5m/s，站姿面（1.7m高度）则普遍在0.2以下，同样满足现行规范的要求，由此可以看出，送风口下置能改善室内的热舒适性环境。

（3）对比上侧送风口的工况可以看出，送风口下置，室内热量