1. 引言
为了使空间内部流场更加均匀,国内外学者对此进行了相关的研究。肖飚 [1] 等研究发现,带扰流件的单风口流场均匀性有所提高。Sadrizadeh S等 [2] 对手术室内垂直送风气流分布进行研究,结果表明增大送风气流的速度可以使得通风效率得到提高,但也会使层流与湍流发生混合,进一步影响室内流场的稳定性。Ye W B [3] 提出了管道均匀气流分布的设计原则,使用有限体积法对设计结果进行验证。Li Tianyu等 [4] [5] 建立变截面管道的运动学模型,通过挤压管道的局部截面得到了一种三维弯管,对弯管进行数值模拟,并且利用商业上常使用的有限元程序进行验证,最终形成了变壁厚管道模型。杨志恒 [6] 对列车空调实验室内的气流组织进行模拟,结果表明,纵向送风方式下风速的变化会对实验室内气流组织的均匀性产生影响,7.9 m/s是最佳的纵向送风风速;可以采用均匀送风方式来改善气流组织分布。滕琴 [7] 介绍了两种不同的均匀送风的形式,通过列举实际应用时的例子,给出了相应的计算方法和考虑因素,分别阐述了两种设计方法的适用场所和考虑条件。肖婷等 [8] 利用CFD模拟方法对改变截面风管的送风结构的流场分布进行研究,模拟结果验证了改变截面风管有利于管内静压分布和出风均匀性。戚新秋 [9] 通过CFD模拟与实验得出了6种送风结构在相同送风条件下某一平面流场的分布特征,得出60˚斜坡送风风道效果最佳,且该结构下送风速度在2~6 m/s的范围用于实验室均匀性最好。宋金蔚等 [10] 提出了在室内增设调节板的方法,用CFD模拟进行对此结构进行验证并且与实验结果进行对比,发现改进后的室内气流组织状况良好。
以上研究对于机房空调实验室气流组织研究较少,对于工艺性空调气流组织的研究多数集中于送风参数、内部结构、模拟方法等影响因素与流场分布的关系,而对于外部结构优化对流场影响的研究较少,没有充分考虑到实验室实际建造尺寸对送风风道设计的要求。本文在上风道方形送风结构基础上提出了10˚倾斜挡板送风结构、12˚倾斜挡板送风结构、14˚倾斜挡板送风结构以及等风量连续斜板送风结构,通过CFD模拟以及实验验证得出适合于工程实际的上风道送风结构。
2. 实验室结构与不同送风风道结构
室外侧实验室的内部尺寸为:10,200 mm (长) * 7000 mm (宽) * 5400 mm (高),送风腔尺寸为6800 mm (长) * 1200 mm (宽) * 2200 mm (高) * 2。在回风口设置2个风阀,当被测机为侧出风时,实验室采用两侧上风道送风、侧回风方式,如图1所示,开启风阀2、关闭风阀1,气流流经空柜进入上侧静压风道,再由静压风道经送风口进入室内,最后经过风阀2回到空柜。
Figure 1. Structural diagram of side air outlet form
图1. 侧出风结构示意图
如图2所示,机房空调实验室5种不同送风风道结构示意图(Pi,Vi分别为各截面的平均静压,平均风速),结构A为方形送风结构,结构B、结构C与结构D分别为在方形送风结构中增设10˚、12˚以及14˚的倾斜挡板。结构E为等风量连续斜板送风风道结构(各水平板高度分别为2.20 m、1.687 m、1.137 m、0.576 m)。
3. 机房空调实验室的数值模拟
3.1. 数值模型建立
本文对流场均匀性的研究是在稳定工况下进行的,因此研究对象可视为稳态、定常流动。为了简化数学物理模型,需要做出以下假设 [1] [11]:
1) 实验室外层结构简化为平面,即不存在厚度问题;2) 由实验室外围保温结构以及观察窗采取一定的保温措施,忽略太阳进入实验室内的辐射;3) 将实验室内的尘埃等污染物视为无质量的物质,对气流的流动不产生影响;4) 将被测机视为内热源。
本文对于非等温送风室内的气流分布选用湍流模型,并且视为定常流动,采用K-e两方程模型作为求解方程模型 [12] [13] [14]。
3.2. 边界条件
本文根据GB/T 19413-2010《计算机和数据处理机房用单元式空气调节机》 [15]、GB/T 17758-2010《单元式空气调节机》 [16] 和被测机的要求,边界条件设置为下:
1) 送风温度设置为35℃,送风速度设为3 m/s;
2) 出口压力设为一个大气压;
3) 由于流过壁面的空气流速速度为零,假设与壁面平行和垂直方向的速度都为零,墙壁和室内隔板等均选择wall边界条件,因不直接与室外进行热交换,视为绝热面,wall并不需要指定任何参数。
3.3. 不同送风风道结构模拟分析
3.3.1. 不同送风结构的出风口模拟分析
对5种不同送风风道结构进行模拟,如图3所示Z = 0.65 m平面速度流场矢量图(X方向表示长度、Y方向表示高度、Z方向表示宽度。
结构A出风口风速在1.0~5.0 m/s,中心区域的速度低于地面速度。结构B在风道两侧增加了10˚的倾斜挡板,出风口风速在1.5~5.0 m/s,送风气流与风道板的角度变大。结构C每个送风口的最大风速范围在4.0~5.0 m/s,差距减小、风速更加均匀。结构D与结构C相比,风速更加均匀。结构E第一个风口的速度较大,达到5 m/s,第二个和第三个风口的速度近似相等,第四个风口速度最小。方形风道内部设置倾斜挡板增大空气流动阻力改变流动状态,使前后出风口流速更加均匀一致,有益于四个出风口的均匀出风,其中增设14˚倾斜挡板的风道结构流场均匀性最好。
3.3.2. 制冷工况下不同送风结构模拟结果
对5种不同送风风道结构在制冷工况(T = 35℃)下选取两个平面,Y = 0.8平面和Y = 1.5平面的温度场和速度场的分布云图,模拟结果如下图所示(图4~8分别为结构(A)~(E)的模拟结果图)。
1) 结构A
图4(a)、图4(b)可看出Y = 0.8平面和Y = 1.5平面墙角温度较低,最低307.5 K,被测机制冷工况下,因为被测机出风影响,出风口位置的温度场发生剧烈变化。主测试区域,Y = 0.8平面温度场变化较Y = 1.5平面平缓,Y = 0.8平面的温度变化范围为307.51~308.45 K,Y = 1.5平面的温度变化范围307.506~308.50 K,误差都在±0.5 K之间,均满足国标要求。
图4(c)、图4(d)中Y = 0.8平面测试区域的速度变化范围在0.26~1.04 m/s之间,Y = 1.5平面测试区域的速度变化范围在0.52~1.56 m/s之间,Y = 1.5平面的气流扰动较大,而Y = 0.8平面的速度场变化平缓。总体分析,平面内存在风速较小和风速不均匀两个问题。
(a) Y = 0.8平面等温线图 (b) Y = 1.5平面等温线图 (c) Y = 0.8平面等速线图 (d) Y = 1.5平面等速线图
Figure 4. Structure A air duct isotherm and isovelocity diagram for each plane when T = 35˚C
图4. T = 35℃结构A风道各平面等温线和等速线图
2) 结构B
图5(a)、图5(b)中受到被测机制冷运行的影响,Y = 1.5平面温度场的变化较Y = 0.8平面剧烈,在被测机出口的相同位置,Y = 1.5平面的温度高于Y = 0.8平面,被测机周围的温度场变化均在308 ± 0.5 K,符合国标对制冷工况温度场均匀性的要求。增加倾斜挡板,相当于改变了风道截面积,图5(c)、图5(d)中,Y = 0.8平面的速度场变化范围为0.52~1.3 m/s,Y = 1.5平面的速度场变化范围为0.78~1.56 m/s,与结构A相比结构B的速度场变化平缓。
3) 结构C
对比图6(a)、图6(b)和图5(c)、图5(d)中温度场的分布发现,结构B与结构C的温度场分布基本一致。
对图6(c)、图6(d)分析发现,气流受到挡板的影响,出风口最高风速超过2.5 m/s,随着高度的下降,速度逐渐减小。从两个平面的速度场分布发现两边墙角的速度偏低,与A结构相比低速区域少。被测机的中心区域,Y = 0.8平面内的风速在0.52~1.56 m/s,Y = 1.5平面的风速范围为0.78~1.56 m/s,与结构A相比,风速均匀性有很大的提高。
(a) Y = 0.8平面等温线图 (b) Y = 1.5平面等温线图 (c) Y = 0.8平面等速线图 (d) Y = 1.5平面等速线图
Figure 5. Structure B air duct isotherm and isovelocity diagram for each plane when T = 35˚C
图5. T = 35℃结构B风道各平面等温线和等速线图
(a) Y = 0.8平面等温线图 (b) Y = 1.5平面等温线图 (c) Y = 0.8平面等速线图 (d) Y = 1.5平面等速线图
Figure 6. Structure C air duct isotherm and isovelocity diagram for each plane when T = 35˚C
图6. T = 35℃结构C风道各平面等温线和等速线图
4) 结构D
Y = 0.8平面被测机周围1.5米范围内温度变化范围为307.529~308.2 K,Y = 1.5平面内温度变化范围为307.522~308.473 K,对比于结构C增设12˚挡板,结构D增设14˚挡板温度变化范围没有太大的影响,但是温度场均匀性的范围提高。对图7(c)、图7(d)进行分析发现,在被测机1.5 m范围内,Y = 0.8平面内的风速在0.78~1.30 m/s,Y = 1.5平面的风速范围为0.78~1.56 m/s,相对于B结构和结构C,速度的变化范围缩小,表明增设14˚挡板会使得速度场均匀性更好。
5) 结构E
由图8(a)、图8(b)温度场分布可知,Y = 0.8平面被测机周围1.5米范围内温度变化范围为307.53~308.16 K,Y = 1.5平面内温度变化范围为307.522~308.25 K,对比于结构D,等风量送风结构对温度场的影响较小,但是温度场稳定的范围增大。
由于风道中加了不同高度的挡板,风速增大,最大为2.6 m/s,在被测机1.5 m范围内,Y = 0.8平面内的风速在0.78~1.30 m/s,Y = 1.5平面的风速范围为0.78~1.56 m/s。与结构D相比,温度场变化范围以及速度场的变化范围没有太大的优化,考虑的工程实施的简便性,结构D更适用于工程实际。
(a) Y = 0.8平面等温线图 (b) Y = 1.5平面等温线图 (c) Y = 0.8平面等速线图 (d) Y = 1.5平面等速线图
Figure 7. Structure D air duct isotherm and isovelocity diagram for each plane when T = 35˚C
图7. T = 35℃结构D风道各平面等温线和等速线图
(a) Y = 0.8平面等温线图 (b) Y = 1.5平面等温线图 (c) Y = 0.8平面等速线图 (d) Y = 1.5平面等速线图
Figure 8. Structure E air duct isotherm and isovelocity diagram for each plane when T = 35˚C
图8. T = 35℃结构E风道各平面等温线和等速线图
5种不同结构送风风道对测试区域温度场均匀性影响较小,主要影响集中在速度场上,方形风道内部设置倾斜挡板,改变空气流动状态,有益于四个出风口的均匀出风,其中增设14°倾斜挡板的风道结构流场均匀性最好。
4. 实验验证
4.1. 实验内容
为了验证结构D模拟结果的可信性以及优化结构的准确性,需要对被测机周围实际的速度和温度进行测量,现场结构D送风风道如图9所示。热电偶测点在测试区域的中间区域布置,本次实验选取Y = 0.8 m、Y = 1.5 m两个平面布置测点,并且两个平面布置相同数量的测点,布置方式一致 [17],测点布置如图10。每个测点的风速测量采用TESTO425热敏风速仪。制冷工况(T = 35℃,Tw = 24℃)达到稳态后,进行数据采集与分析。
Figure 10. Layout of measuring points in the air-conditioning laboratory
图10. 实验室测点布置
4.2. 实验与模拟对比分析
图11所示为结构D (方形送风风道内增加14˚倾斜挡板)的Y = 0.8 m,Y = 1.5 m平面流场的实测及模拟对比图。由图11(a)可知,结构D温度的模拟与实测值最大误差为0.13%。图11(b)中所示,速度实测与模拟结果的最大误差点在1号测点处,差值为0.34 m/s,除去该点处误差较大外,其他处试验与模拟速度值的曲线拟合程度较高。
如图11(c)所示,结构D温度的模拟与实测值的最大误差为0.17%。图11(d)中显示,速度的实测值与模拟值均呈现明显的起伏波动,实测最大值为1.37 m/s,最小值为0.60 m/s;模拟结果的最大值为1.56 m/s,最小值为0.52 m/s。速度的实测与模拟结果的最大误差点在47号测点处,差值为0.53 m/s,误差在可接受范围之内。
(a) Y = 0.8平面温度实测与模拟对比 (b) Y = 0.8平面速度实测与模拟对比(c) Y = 1.5平面温度实测与模拟对比 (d) Y = 1.5平面速度实测与模拟对比
Figure 11. Comparison of simulated with measured for Y = 0.8 m and Y = 1.5 m planes
图11. Y = 0.8 m、Y = 1.5 m平面的模拟与实测对比
由模拟仿真以及试验验证可以发现结构D (增设14˚倾斜挡板)可以满足人工环境试验室温度场和速度场均匀性的要求,并且相对于其他几种结构速度变化范围较小,温度场分布更加均匀,相比于结构E (等风量连续斜板送风结构),结构D在工程上施工便利,经济型更好。
5. 结论
1) 方形风道内增设倾斜挡板,送风气流与风道板的角度变大,提高流场均匀性,有益于出风口均匀出风。
2) 变截面积送风风道使实验室温度场更加均匀,满足温度场 ± 0.5 K的精度要求,在主测试区域,结构D温度场的温差幅度在±0.3 K,满足高精度的测试要求。
3) 结构D风道送风,Y = 0.8平面所选测点温度实测值与模拟值最大的误差为0.13%,速度的实测值与模拟值最大误差为0.34 m/s;Y = 1.5平面温度的实测值与模拟值最大误差为0.17%,速度的实测值与模拟值的最大误差为0.53。结构D温度场和速度场的均匀性更好,在工程上易于实现,因此工程上推荐该结构的设计。