1. 引言
空调在我们的生活中变成了不可或缺的一部分,空调实验室就是专门用来检测新研发的空调的性能。实验室能为需要检测的制冷设备提供国家标准规定的实验工况,温度场分布的均匀性将会影响到制冷设备测试的稳定性和精确性。对于环境实验室的实时流场的检测比较复杂,检测周期较长,同时由于环境室内空气流速比较小,室内轻微的扰动都有可能影响到最后的测试结果,不利于实验室内的气流组织结构的分析与改进 [1]。同时气流组织的分布情况又会影响到被测机工作周围环境的温湿度精度、风速以及局部温差等的主要参数的变化。CFD (Computational Fluid Dynamics)即计算流体力学,余晓明等 [2] 用CFD模拟软件针对焓差实验室内气流组织的分布情况进行数值模拟研究,在此基础上对测试环境室内的气流组织结构进行优化设计;宋金蔚等 [3] 运用CFD模拟计算、定性理论分析和实验测试研究三种方法,针对某高速列车内气流组织的分布及送风均匀性进行了深入研究。CFD方法可方便快速地对各种工况进行模拟 [4],被应用于暖通空调等领域进行模拟分析 [5] [6]。
本文以某企业已建成的空调实验室为研究对象,布置合适的测点测量速度和温度值,建立模型利用CFD软件对该实验室不同送风速度下的温度场和速度场进行模拟。
2. 模型的建立
2.1. 物理模型
实验室总体尺寸为5300 mm (宽) × 6700 mm (长) × 4550 mm (高),如图1所示,同时设置双开平门,尺寸大小为:(1000 + 1000) mm (宽) × 2200 mm (高),门上带防雾中空观察窗。另外房间设有两个风阀,风阀1尺寸2500 mm (宽) × 1250 mm (高)、风阀2尺寸2000 mm (宽) × 2500 mm (高),送风腔尺寸为800 mm (宽) × 6700 mm (长) × 1050 mm高 × 2。
Figure 1. Structural diagram of side air outlet form (air duct supply)
图1. 侧出风形式结构示意图(风道送风)
2.2. 数学模型
大量试验表明:空调室内的气流绝大部分是湍流。基于此本文采用湍流模型对室内的气流组织进行模拟,标准k-ε方程模拟精度高,相对符合实际情况。标准k-ε模型一般只用于完全湍流的流动过程的模拟,通常情况下忽略分子黏性的影响。其湍动能与耗散率方程如下:
(1)
(2)
式中:G——湍动能产生;
——可压湍流膨胀影响耗散率的大小;
——湍流粘性系数;
其中默认值常数:
,
,
,
。
2.3. 边界条件
1) 进口边界速度入口边界条件;
2) 出口边界压力出口边界条件;
3) 壁面边界条件壁面上对速度取无滑移条件,绝热壁面。
本文所用到的边界条件为velocity-inlet、pressure-outlet、wall、porous-jump。循环风机出口(即实验室入口)和被测机进、出口均选择Magnitude,Normal to Boundary定义入口流速,输入速度矢量的大小就可以定义速度,速度值以实测为准。选择pressure-outlet作为实验室回风口处的边界条件,出口压力设为一个大气压。实验室内四周壁面以及室内隔板等均选择wall边界条件,因不直接与室外进行热交换,视为绝热面,wall并不需要指定任何参数,只需默认即可。
3. 模型验证
实验室区域距离地面0.8 m、1.6 m定义A、B两平面,这两个平面是放置被测机的主要空间,也是在测试中对温度以及速度的均匀性要求较高的区域,因此本文对实验室区域内所选取的这两个平面进行分析。由于本文后面的优化主要针对制冷工况 [7] (T = 35℃)下的侧出风被测机实验室内部温度场及速度场的分布,所以模型验证也只针对制冷工况(T = 35℃)的送风模型。将模拟得到速度分布云图的数据与实验室相同测点的数据进行对比,来验证模拟方法是否可以用于对实验室温度场均匀性优化。
图2是制冷工况下(T = 35℃)现有风道送风结构下侧出风被测机的模拟结果云图,分别显示了A、B平面上的温度场、速度场分布情况。A、B内平面内被测机出风口附近温度变化比较明显,主测试区域温度梯度变化A平面比较平缓,B平面内温度变化梯度较大。由于本实验是在被测机运行状态下进行测试,被测机出风温度较高,与室内空气的换热也会影响实验室内温度场的分布。在实际的实验过程中,需要保证被测机入口侧的温度场分布的均匀性;从图2可以看出A、B平面内墙角处温度偏低,虽然A、B平面温度波动范围较大,但是被测机进风口平面处测点温度波动范围均在±0.5℃之间。对图2(c)、图2(d) A、B平面内的速度云图分析,A平面内的空气流动相对稳定,B平面速度变化梯度比A平面大,两平面内主测试区域速度均在0.52 m/s~1.56 m/s区间内变化,符合国标要求。
(a) A平面等温线图 (b) B平面等温线图(c) A平面等速线图 (d) B平面等速线图
Figure 2. Isotherms and isotherms of each plane under the air model of T = 35˚C orifice side
图2. T = 35℃孔板侧出风模型下各平面的等温线和等速线图
此次实验采用精度较高的美国TSI风速仪VELOCICALC 9535进行温度和速度测量。在A、B平面,进行测点布置,每个平面平均布置17个测点。1~8号测点布置在靠近被测机进风口处的位置,空气流动速度在0.6 m/s~1.4 m/s内波动,9~16号测点布置在靠近墙壁的位置,该区域空气流动速度大,而且有较明显的波动,但空气流动速度基本在1.8 m/s以内变化;由表1可看出,模拟计算得到的数据与实验测得的数据误差很小,绝对误差的最大值为0.11,除3号和7号测点绝对误差超过0.1,大部分绝对误差控制在0.02 m/s~0.07 m/s内,相对误差基本在10%以下。综合考虑各方面因素,模拟计算应用在实际试验过程中是可行的,模拟计算得到的结果基本可以如实反映实际试验中实验室内部温度场和速度场的分布。
由于本文测量数据较多,只选取B平面部分测点测得的实验数据与模拟结果中相应测点的数据进行对比分析。被测机进风口处选择1、3、5、6和8号测点,靠近墙壁处的位置选取9、11、13、15和17号测点。从图3中我们可以看出B平面相同位置的测点温度存在温差,温差最大值为1.4 K,位于被测机出风口的1号测点,在被测机进风口处两平面温差最小。而且,B平面靠近墙壁的位置温度梯度均大于靠近被测机位置,温度变化更加明显。B平面实际测量值与模拟计算值差值控制在0.05 K~0.25 K,由于实际试验测量存在误差,因此认为应用的物理模型和计算方法可以应用于实验室温度场与速度场的优化研究,并且最后得到的计算结果与实际情况相似。此次对比分析对模拟计算的可行性进行了验证,可以使用该物理模型和计算方法对实验室内部的温度场和速度场进行验证分析,模拟结果基本与试验情况相似。
Table 1. Comparison of measured and simulated values of wind speed in each plane
表1. 各平面风速实测值与模拟值对比
4. 数值计算结果及分析
对于实验室内空气流场的分布而言,风速的大小是影响室内流场均匀度分布的一个重要参数。通过改变空柜送风速度,使实验室测试区域风速产生变化,进而改变实验室温度场均匀性。原送风口的风速为3 m/s,在A平面内的测量区域内风速最低到0.4 m/s,清华大学的朱颖心教授提出测试区域的风速在0.5 m/s以上时,机械送风的效果才与实际送风的效果基本相同 [8],但测试区域的风速最高不宜超过2 m/s。在第一个变化模型中将风口面积缩小12%,将风速由3 m/s提高到3.5 m/s,在第二个变化模型中将风口面积缩小20%,将风速由3 m/s提高到4 m/s。
(a) A平面速度(b) B平面温度
Figure 3. Comparison of measured and simulated values of velocity and temperature in plane a and plane B
图3. A平面及B平面速度和温度实测值和模拟值得比较
图4是制冷工况下(T = 35℃)送风速度为3.5 m/s时侧出风被测机的模拟计算结果云图;图5是制冷工况下(T = 35℃)送风速度为4 m/s时,侧出风被测机的模拟计算结果云图,分别对不同风速时A、B两个平面上的温度场、速度场进行分析。结果如下图所示:
(a) A平面等温线图 (b) B平面等温线图(c) A平面等速线图(d) B平面等速线图
Figure 4. Isotherms and isotherms of each plane when the air supply speed is 3.5 m/s
图4. 送风速度为3.5 m/s时各平面等温线和等速线图
(a) A平面等温线图 (b) B平面等温线图(c) A平面等速线图 (d) B平面等速线图
Figure 5. Isotherms and isotherms of each plane when the air supply speed is 4 m/s
图5. 送风速度为4 m/s时各平面等温线和等速线图
通过图4(a)、图4(b)与图5(a)、图5(b)对不同送风速度时实验室温度场变化进行对比分析。当送风速度提高到3.5 m/s时,A、B平面测试区域的温度场更加均匀,根据此变化可以认为适当的提高风速,可以对实验室温度场产生影响,使得温度场更加均匀。但当送风速度达到4 m/s时,A、B平面的温度波动相对于3.5 m/s时又变大,均匀性降低。该现象证明,当送风速度提高时,温度场的变化并不是一个均匀变化的过程,而是在风速提高到一定程度后,温度变化的情况由不太均匀到均匀,然后再到变化剧烈的过程。
通过图4(c)、图4(d)与图5(c)、图5(d)对不同送风速度时实验室速度场变化进行对比分析。当送风速度为3.5m/s时,测试区域风速相对于之前有所提高,主要测试区域风速在0.5 m/s~2 m/s之间波动,气流速度完全符合实验室要求。但当送风速度提高到4 m/s时,在B平面靠近墙壁处部分区域风速远远超过了2.5 m/s,已经不符合对被测机进行测试的要求。
5. 结论
本文以现有的风道送风焓差实验室为基础,在实验室内部两个不同高度的平面上布置测点,在制冷工况下进行测试,运用CFD软件对实验室稳定后的工况进行模拟,将模拟结果与之前实际测得的实验室温度和速度数据进行对比,验证该数值模拟计算方法的可行性。对于实验室内气流组织而言,风速的大小是影响室内流场均匀度分布的一个重要参数,因此提出改变风速来优化温度场和速度场的均匀性,通过前面的实验和模拟可得出以下结论:
1) 根据现有实验室的结构,建立相对应的三维模型,并运用数值模拟软件进行仿真计算,将得到的模拟云图结果与实验测试数据进行对比。通过对比结果的分析,验证了模拟方法的可行性,
2) 通过改变送风速度来改善实验室内部温度场均匀性时发现当送风速度为3.5 m/s时,实验室内部流场均匀性与元送风速度3 m/s相比得到明显提升,但当继续增大送风速度时,温度场变化逐渐剧烈,速度也会超出实验室要求的最大风速。因此,小幅度地提高风速会使实验室内部流场更加均匀,但当风速达到一定值时,就会使实验室温度场均匀性下降。