壁挂式空调制冷与制热运行的数值模拟.doc

壁挂式空调制冷与制热运行的数值模拟摘要：本文利用湍流K-方程计算了壁挂式空调的制冷与制热运行时的三维速度与温度场，分别比较了制冷与制热运行时三个不同坐标平面内的速度与温度分布，为分析空调的设计运行参数提供依据。关键词：空调器，CFD，流场 1 引言空调设计的最终目的是以经济技术合理的系统设计以及设备选型实现所要求的室内气候环境（温度、气流、污染物浓度等的分布）。为实现对这些环境参数的合理控制，有必要把握其分布特征，CFD计算是除模型实验以外的可详细研究室内气流分布特征的重要手段1 2 3。相比传统的模型实验与经验公式预测流体的流动与传热，CFD计算具有成本低、速度快、资料完备等优点。很多学者对空调的气流组织进行了较深入地研究，这些研究大多采用稳态的模型，其研究大多在不考虑室内与室外的热源，房间内部主要是空态的，与实际运行时的情况有一定的区别。本文的研究是比较壁挂式空调在窗辐射、室内热源与墙体热源的情况下，对室内三维流场、速度场进行计算模拟，分析室内的流场与温度场。2 几何模型为了简化实际问题，便于分析，在建立数学模型前对室内气体的流动作如下假设：1）、室内气体满足牛顿内摩擦定律，为牛顿流体；2）、室内气体温度变化不大，密度可视为常数；3）、室内气体的流动形式为稳态紊流；4）、在紊流中心区，忽略能量方程中由于粘性作用而引起的能量耗散。以传统室内挂壁式空调为研究对象，对其温度场进行数值模拟计算。建立房间三维立体模型，长、宽、高分别为5m、4m、3m。见图1所示，空调安装高度为2.5 m。X、Y、Z三方向坐标也表示于图中。设定空调的原始参数如下：几何尺寸：800mm300mm200mm出风角度：45；出风速度：送风速度2.6m/s； 出风温度：制冷时18C。制热时18C。根据已知参数，建立适当的房间有限元模型。为使结果尽可能接近真实情况，边界条件设定为：六面墙处于绝热状态，制冷运行时的壁面固定温度为30C（制热时为5C），空调本身不发热，房间内部无任何热源。图1 空调送风示意图数学模型采用k双方程模型4模型的基本微分方程包括连续方程、动量方程、紊动能k方程及紊动能耗散率方程。分别表示如下：k方程：方程： 式中： 和 分别为容积分数平均密度和分子粘性系数； 为紊流粘性系数，它可由紊动能 和紊动耗散率 求出： 式中， 0.09为经验常数。 和 分别为k和的紊流普朗特数，其中 =1.0， =1.3。 1.44和 1.92为方程常数。G为由平均速度梯度引起的紊动能产生项，定义为： 数值求解算法采用控制容积法，在每个控制容积单元中对微分方程进行积分，再将积分方程线性化，把控制方程离散为可以数值求解的代数方程，就可以得到相应各未知变量，如压力、速度、紊动能k及其耗散率等变量的代数方程组，然后再对方程组进行求解，就可以求出各未知变量及其空气的温度（或焓值和）。 3 模拟结果与分析有了以上的数学模型和物理边界条件，笔者模拟了该空调房间的流场和温度场，并选取了其中具有代表性的断面进行分析，模拟结果见以下各图（图2图7）。a)速度矢量分布 b)温度分布图2 空调制冷运行时对称立面(XY面)速度和温度分布a)速度矢量分布 b)温度分布图3 空调制冷运行时空调下方YZ面速度和温度分布a)速度矢量分布 b)温度分布图4 空调制冷运行时XZ面速度和温度分布a)速度矢量分布 b)温度分布图5 空调制热运行时XZ面速度和温度分布a)速度矢量分布 b)温度分布图6 空调制热运行时空调下方YZ面速度和温度分布a)速度矢量分布 b)温度分布图7 空调制热运行时XZ面速度和温度分布由上述空调的两种运行状况可以看出：在空调制冷运行时由于冷空气自身的下坠，靠近地面的水平面的温度较低，由于空调射流的自身卷吸作用，加速冷热空气的混合，从而使其室内温度场分布较为均匀。在空调制热运行时，由于热空气的浮力作用，热空气难以到达地面附近，在人体垂直方向上存在温度梯度，头部比脚部温度高，这是由于冷空气较热空气密度大，向下趋势明显，从而导致房间下部要比高处降温快。温度差会造成人体局部性热感不适，随着时间推移温度梯度会逐步渐小。4 结论本文利用湍流K-方程计算了壁挂式空调的制冷与制热运行时的三维速度与温度场，分别比较了制冷与制热运行时三个不同坐标平面内的速度与温度分布，空调的制冷与制热运行的两个工矿况中，制冷运行中出现了冷风下坠现象，在制热运行中有部分热空气，未经充分热交换而直接回流，能量利用不充分，温度均匀性也不好。利用CFD计算可以改善室内的气流组织，通过改变空调不同气流出口速度，不同的出口角度从而得到合理的气流组织形式，这是采用传统的射流理论无法做到的。The numerical simulation study temperature distribution and fluid flow of Wall Hanging Air-condition Huang Zhikai Huang Deshuang Zhang Jun1.Anhui Industry Polytechnic College 2. Institute of Intelligent Machines, Chinese Academy of Sciences. Abstract: This paper is an application of CFD(computer fluid dynamics) in the turbulent flow of air conditioned room in which the numerical model is introduced. Based on theoretically analysis of physical and mathematical model, the K-two-equation turbulence model are adopted to simulate the characteristics of temperature distribution and fluid flow inside the air-conditioned room. The results indicate that CFD is an effective way to predict temperature distribution and fluid flow inside the air-conditioned room and it also provide basis for air-condition system optimal design.Keywords: air-condition, CFD; fluid flow参考文献1、汤广发等，室内气流数值计算及模型试验M，长沙：湖南大学出版社，19892、Davidson, L. And Nielsen, P. V. Large Eddy Simulations of the Flow in a Three Dimensional Ventilated Room. Proc Of ROOMVENT96 on Air Distribution in Rooms, Japan, 19963、Lannder B E，Spalding D BMathematical Models of TurbulenceM.Academic Press，London and New Tork，1972：90110. 4、张智力，吴喜平，CFD基本算法及其在暖通空调领域