壁面温度对室内空气流动过程的影响.docVIP

壁面温度对室内空气流动过程的影响 摘 要：为了分析在进风位置、进风温度以及射流方向不同的条件下，房间壁温对室内空气流动的影响，利用CFD计算、流动可视化和定量测量相结合的办法研究了不同工况下室内空气的流动形态、速度以及模型整体换热情况。结果表明：进风温度变化时，壁面温度对水平送风方式下的空气流动形态影响较大，对垂直送风方式影响较小；进风温度的降低引起水平进风方式的自由射流区缩短，右侧垂直进风方式下射流向右侧壁面偏转，同时浮升力引起的自然对流作用增强，换热温差加大，增加了围护结构漏热量。 关键词：空调通风；流动形态；CFD模拟；可视化实验 通风空调房间内的空气温度、流速以及污染物浓度等对建筑物能耗、室内空气品质和人体健康至关重要。自20世纪40年代以来，研究人员运用不同方法来预测机械通风房间内的气流组织结构[1-6]。Chen Qinyan[3,4]利用CFD模拟研究了不同的通风方式、风口结构对房间内空气流动的影响；Fitzner[5]通过原型实验得出室内气流主要受风口结构、位置，空气流速以及房间热负荷的影响，而排风口位置对其影响较小；Zhang[6]实验研究了水平送风条件下，不同进口风速、温度以及房间热负荷对室内空气速度和湍流特性的影响。然而，在进风口位置、进风射流方向以及进风温度不同的条件下，房间壁面温度对室内空气流动形态的影响规律未作具体研究。本文通过CFD模拟和可视化模型实验相结合的办法研究了不同的风口位置、射流方向以及进风温度下，房间壁面温度对室内空气流动过程的影响。 1 实验模型及数值模拟方法 图1 实验模型结构示意图 设定如图1所示的坐标系统，模型的在x、y、z方向的长度分别为1.2 m、0.6 m（x=0 m）4 mm导热铝板组成的复合加热板，通过调节输入电压可以控制实验所需的壁温，而模型的其他外表面通过粘贴厚度为100 mm的保温材料来模拟绝热边界条件。本文讨论A、B两种不同位置的进风方式，风口尺寸都为240 mm×24 mm，且以z=0.4 m的平面前后对称。其中进风口A设在x=1.2 m的平面，水平中心线距模型顶面0.1 m；进风口B设置在y=0.6 m的顶面，右沿距模型x=1.2 m的平面0.1 m；排风口设在y=0.6 m的顶面，左沿距x=0 m的平面0.1 m。 通常情况下房间内的平均流速小于0.15 m/s，只有射流区域的流速达到相对容易进行测量分辨的0.4 m/s以上。为此，房间内部的风速、温度测量点主要布置在射流方向上。对于A型进风方式，7个测量点布置在y=0.3 m、z=0.4 m的直线上，其中1点、7点x坐标分别为1.088 m、0.112 m，其余5点x坐标在1点和7点之间均匀布置；对于B型进风方式，5个测量点均匀布置在x=1.088 m、z=0.4 m的直线上。 表1给出本实验的研究工况，加热表面的温度为40 ℃、风量为0.0078 m3/s。 表1实验研究工况 工况序号 Ⅰ Ⅱ Ⅲ Ⅳ 风口位置 A A B B 进风温度/℃ 22 14.5 22 14.5 本文采用FLUENT6.0商业软件进行计算，入流湍动能kin为来流平均动能的1％，湍流耗散εin可按照入口湍流雷诺数Ret以及湍流特征长度之间的关系来选取[7]。利用SIMPLEC算法进行压力——速度解耦，PRESTO！ QUICK离散格式处理控制方程；采用RNG k-ε湍流模型，用壁面函数法考虑边界对流动的影响。在划分网格时，尽可能采用正六面体网格单元，固体壁面相邻的第一个网格宽度限定为为1mm，最宽一级的网格宽度≤8cm，在进风口及壁面附近区域加密网格，将相邻两网格间距的尺度比严格控制地控制在0.8～1.25。 2、结果及分析 2．1 CFD计算及可视化实验结果 工况Ⅰ、工况Ⅱ下的z方向中心截面（z=0.4 m，下同）的CFD模拟结果、流动可视化照片如图2所示。可以看出，壁面温度对模型内空气的流动形态有明显的影响。进风温度较高时，自由射流区较长，射流进入到1/2的宽度（x=0.6 m的平面）后，开始向下弯曲，射流能够冲刷到左侧壁面，并且在0.3 mx0.6 m、0 my0.3 m范围内形成一个逆时针涡；而进风温度较低时，射流进入1/4的宽度后，便开始向下弯曲，射流冲刷到底面，不能形成明显的涡。产生这种差别的原因在于当进风温度较高时，空气受到方向向下的重力作用较小；反之，进风温度较低时，空气受到的重力作用较大。 图2-a 工况Ⅰ下的流动形态的计算结果 图2-b 工况Ⅰ下的流动形态的可视化结果 图2-b 工况Ⅱ下流动形态的计算结果 图2-b 工况Ⅱ下的流动形态的可视化结果 Ⅲ、工况Ⅳ下的CFD模拟结果以及流动可视化照片。对于垂直送风方式，壁面温度对模型内空气的流动形态影响较小。只是当进风温度的较低时