空调房间室内热的环境的三维数值模拟及实验的研究.docVIP

PAGE 483空调房间室内热环境的三维数值模拟及实验研究山东建筑工程学院 马秀力?，赵菊，肖勇全摘 要：本文通过建立相应的数学物理模型，利用PHOENICS对空调房间的室内热环境进行了数值模拟，得出了空气温度、流速的分布图，并对模拟房间进行了实验测试，对模拟结果与实测结果进行了分析、比较。分析结果表明模拟值与实测值的吻合度很好，这表明所建立的模型是正确的，模拟结果是可信的。关键词：计算流体力学 数值模拟 PHOENICS 实验测试 数值分析0 引言计算流体动力学简称CFD(Computational Fluid Dynamics)在暖通空调领域的应用是运用流体动力学的基本原理，通过建立数学物理模型，根据提供的合理的边界条件和参数，对空调区域内气流的速度场、温度场、压力场等进行模拟计算。而室内空气的速度场、温度场又是空调房间室内气流组织设计及空调房间室内舒适环境评价的基础。建筑室内的气流分布和温度分布是体现舒适和卫生的空气环境的主要指标之一。如果能够在规划和设计阶段即可正确而详细的预测和评价室内气流组织及温度分布，不仅可以实现现代暖通空调系统的优化设计和运行管理，提高室内热舒适性和室内空气品质，而且对系统乃至整个建筑物的节能也具有重要的指导意义。在国外，CFD数值模拟技术已被广泛应用在工程领域，成为设计上不可缺少的技术手段；而在国内CFD技术的研究和应用也已逐步受到重视。本文就是利用CFD模拟计算软件PHOENICS，对立柜式空调房间的温度场、速度场进行预测和分析，并通过实验测试加以比较验证。1 房间简介 模拟计算的空调房间为本学院的一个多媒体会议室，房间大小为11m×7.8m×3.6m，位于学院主体楼一层的北侧，除北墙为外墙外，其余都为内墙，南面两个门朝向走廊，北面有三个窗户，模拟时门、窗都是关闭的。房间内的布置如图1所示，房间东、西两侧各对称安放一台“美的”世纪星立柜式空调，此外讲台上还有多媒体设备，因为一直没有开，不考虑其散热量，因此没有在图上画出。物体大小、位置都按照实际物体尺寸和在房间内的实际位置布置。图中人所在的位置为实验测试时测点的位置，在实际测试时整个房间是空的，并没有如图所示的人存在，所以在模拟过程中没有考虑内热源。同时由于窗户都在背阴处，没有太阳辐射，而且是关闭的，所以在模拟时与北外墙相同，并没有考虑额外的传热。 图1 模拟房间俯视图2 流体力学模型 空气在房间里的流动是一个很复杂的物理过程，涉及到气体流动、传热传质等多方面的内容，因房间内气体流动是湍流过程，所以湍流模型是空调房间气流模拟的基础。本文选用的是应用于紊流自然对流时的模型。在模拟过程中作以下假设：室内空气为低速不可压缩气体；气流流动为紊态流动；忽略能量方程中由于粘性作用引起的能量耗散。室内无内热源。数值模拟中的基本控制方程包括连续性方程、动量方程、能量方程和k方程和方程，可表示为[1]：连续性方程： （1）动量方程： （2）能量方程： （3）k方程： （4） EMBED Equation.3 方程： （5）式中：为流速矢量在三个坐标上的分量，分别是，，2，3分别代表三维直角坐标系的三个方向。是一数学符号，是体积膨胀系数，是紊流粘性系数，，为层流粘性系数，。流体的速度矢量在三个坐标上的分量分别为，，2，3分别代表三维直角坐标系的三个方向。在这些控制方程中包括6个经验常数，其取值如表1所示。计算时采用稳态，方程（1）～（5）中左侧第一项为0。 表1 紊流模型经验常数[2，3]常数值0.091.441.920.9～1.01.01.33. 边界条件PHOENICS中的辐射换热模型选用的是基于辐射换热基本公式：斯蒂芬-波尔兹定律的IMMERSOL模型[3]，模型方程描述如下：表示固体间辐射温度时，方程为： （7）式中：为综合辐射率，代表单位体积内获得或损失辐射能量的净值，W/m2；为第一、二项流体表面平均辐射力，W/m2；为介质的吸收率；为介质的散射系数；为用于计算的两个固体辐射表面之间的距离，m。该模型的几何尺寸、物性参数、边界条件由PHOENICS软件Menu对话框输入，模拟时考虑墙体内表面辐射的影响，以边界条件的形式输入。在计算过程中，引用Boussinesq假设，即认为流体密度的变化仅对浮升力项产生影响。进行室内气流模拟时，最主要的边界条件是壁面的温度，它是通过人机对话从ATTRIBUTES对话框中直接输入的。壁面温度的取值来自于笔者在该房间的测试结果，该值为每个墙