送风参数对老化箱内部温湿度场的影响研究
Study on influence of air supply parameters on temperature and humidity field inside aging chamber
通讯作者:
收稿日期: 2023-10-16 修回日期: 2023-12-14
基金资助: |
|
Received: 2023-10-16 Revised: 2023-12-14
作者简介 About authors
孙 敏(1998—),女,湖南娄底人,硕士生,从事热力学仿真研究,E-mail:
关键词:
Keywords:
本文引用格式
孙敏, 卢丰源, 赵宇轩, 王青春, 陈忠加.
SUN Min, LU Fengyuan, ZHAO Yuxuan, WANG Qingchun, CHEN Zhongjia.
随着机械设备的应用场景不断拓展,机械设备电子元件的老化问题逐渐凸显。通过老化箱筛选出性能稳定的电子元件是目前最常见的一种老化筛选方法。开展老化筛选的重要前提是保证老化箱内部流场的均匀性。CFD(computational fluid dynamics,计算流体力学)仿真模拟是一种高效的数值模拟方法[1],可有效探究老化箱内部温湿度场的分布情况。
近年来,国内外学者针对烘房、发电机房等设施的送风参数优化开展了大量研究。王振文等[2]对热泵烘房的送风工艺参数进行了优化,结果表明,排风速度为6 m/s时热泵烘房内的气流分布较均匀。曹如玥[3]对蜂箱调温过程进行了数值模拟,结果显示:送风速度对蜂箱流体域内湍流强度的影响不显著;通风孔径对湍流强度的影响最为显著,其次是通风位置。林海等[4]综合考虑了多项工艺参数,通过仿真优化得到:当初始温度为25 ℃时,送风温度为207 ℃、送风速度为5 m/s的参数组合下汽车车身烘房的能耗最低。魏亚兴[5]选取了10个对矿井通风巷热交换产生影响的因素,以温度作为评价指标开展了正交试验,结果显示,进风温度对巷道终点温度的影响较大,且两者呈正相关。杜林昕等[6]为探究进风量、进风角度和进风口位置对烤房内部气流分布情况的影响,设计了相应的正交试验,结果表明,上述3个送风参数改变均会对烤房内部的气体速度产生显著影响,影响的主次顺序为进风口位置、进风角度和进风量。Li等[7]选择送风口高度、送风口宽度、送风速度和热源强度作为影响因素,通过开展正交试验来优化发电机房的通风性能,经极差分析和方差分析可知,送风口高度对发电机房内部的流场和通风性能有显著影响。
综上所述,送风温度、送风口位置、送风口高度对烘箱等热循环设备内部流场的影响显著。基于此,本文将探究送风温度、送风速度和送风角度这3个参数对老化箱内部温湿度场的影响。通过设计9组正交试验来系统性地探讨不同送风温度、送风速度和送风角度下老化箱内部温度场、湿度场的分布情况,旨在为老化箱工艺参数的优化设计提供理论依据。
1 老化箱内部流场计算理论
1.1 湍流模型
一般情况下,对于老化箱内部的流体流动,通常认为雷诺数Re
式中:V为流道内流体的平均速度,m/s;ρ为流体密度,kg/m3;L为特征长度,m;μ为动力黏度,Pa·s;Dr为流道特征尺寸,m;
老化箱的送风口为矩形,对于非圆形管道内的流动,
1.2 基本控制方程
在流体力学领域,流体流动和传热传质均遵循三大基本守恒定律,即质量守恒、动量守恒和能量守恒定律。
1)质量守恒方程。
流体流动在时间和空间尺度上的求解方法基于质量守恒方程。质量守恒是指流入和流出某个微元体的流体质量相等,可表示为:
式中:t为时间,s;u、v、w为x、y、z方向上流体的速度分量,m/s;SM 为质量源项,kg/(m3·s)。由于老化箱内会发生水的相变,因此引入水蒸气质量源项,表示水由气相转化为液相所产生的质量变化。
通过前文分析可知,老化箱内部的流体为不可压缩流体,即流体密度为常数,因此质量守恒方程可简化为:
2)动量守恒方程。
动量守恒是指一个微元体的始末动量变化量等于该运动过程中微元体所受力之和。Naiver-Stokes方程(N-S方程)是描述黏性不可压缩流体动量守恒的运动方程。本文老化箱内部的不可压缩流体在x、y、z方向上的动量方程[9]可简化为:
式中:p为微元体上的压力,N;
3)能量守恒方程。
能量守恒是指流动过程中流体的总能量保持不变。能量守恒方程可表示为:
式中:cp 为比热容,kJ/(kg·℃);T为温度,℃;k为流体传热系数,W/(m2·℃);ST 为黏性耗散项[10],W。
1.3 多相流模型
2 老化箱送风参数的正交试验方案
表1 老化箱送风参数的正交试验因素水平表
Table 1
水平 | 因素 | ||
---|---|---|---|
A(送风温度)/ ℃ | B(送风速度)/(m/s) | C(送风角度)/(º) | |
1 | 50 | 8 | -10 |
2 | 70 | 10 | 0 |
3 | 90 | 12 | 10 |
表2 老化箱送风参数的正交试验方案
Table 2
序号 | 因素 | 方案 | ||
---|---|---|---|---|
A/℃ | B/(m/s) | C/(º) | ||
1 | 50 | 8 | -10 | A1B1C1 |
2 | 50 | 10 | 0 | A1B2C2 |
3 | 50 | 12 | 10 | A1B3C3 |
4 | 70 | 8 | 0 | A2B1C2 |
5 | 70 | 10 | 10 | A2B2C3 |
6 | 70 | 12 | -10 | A2B3C1 |
7 | 90 | 8 | 10 | A3B1C3 |
8 | 90 | 10 | -10 | A3B2C1 |
9 | 90 | 12 | 0 | A3B3C2 |
3 老化箱内部流场仿真模型构建
3.1 老化箱三维模型
图1
图2
3.2 网格划分
本文的老化箱三维模型较为规则,以及采用结构化网格的计算方式耗时短且计算效率高,故采用结构化网格对老化箱三维模型进行划分。考虑到网格数量过多会增加计算时长,在保证仿真结果精度不受影响的前提下,选择合适尺寸的网格单元对老化箱三维模型进行划分,最终的网格数量为160万个,如图3所示。
图3
3.3 边界条件设置
此外,虽然在设计老化箱箱体壁面时加入了保温材料,但是箱内的热量仍会通过壁面传递到外界,故将壁面的热边界条件设为固定温度壁面。
4 老化箱内部流场仿真结果分析
4.1 正交试验结果分析
图4
图4
老化箱内距底面600 mm处水平面的流线图
Fig.4
Streamline diagram of horizontal plane at a distance of 600 mm from bottom inside aging chamber
图5
图5
老化箱内距底面600 mm处水平面的温度云图
Fig.5
Temperature cloud map of horizontal plane at a distance of 600 mm from bottom inside aging chamber
由图4可以看出,在9组送风方案下,老化箱内部的气流运动趋势基本一致,即气流从送风口进入老化箱后向下移动,在接触到壁面后往上散开,总体上气流运动大致呈上下来回移动。通过对比发现,在A3B1C3、A3B2C1、A3B3C2这3组方案下,老化箱内距底面600 mm处水平面上的气流分布相对较均匀。
对不同送风方案下老化箱内部流场的仿真结果进行处理,得到对应的老化箱性能评价指标,结果如表3所示。
表3 不同送风方案下老化箱的性能评价指标
Table 3
序号 | 方案 | 能量利用系数 | 温度不均匀系数 | 相对湿度不均匀系数 |
---|---|---|---|---|
1 | A1B1C1 | 0.909 6 | 0.077 1 | 0.125 2 |
2 | A1B2C2 | 1.058 6 | 0.078 8 | 0.119 7 |
3 | A1B3C3 | 1.045 0 | 0.080 2 | 0.127 5 |
4 | A2B1C2 | 1.028 1 | 0.054 3 | 0.104 8 |
5 | A2B2C3 | 1.029 1 | 0.054 1 | 0.106 4 |
6 | A2B3C1 | 0.938 8 | 0.050 9 | 0.106 8 |
7 | A3B1C3 | 1.061 9 | 0.040 2 | 0.092 0 |
8 | A3B2C1 | 0.983 4 | 0.038 8 | 0.098 7 |
9 | A3B3C2 | 1.051 3 | 0.039 5 | 0.099 0 |
4.2 正交试验结果极差分析
表4 各因素对老化箱能量利用系数影响的极差分析结果
Table 4
分析指标 | 因素A | 因素B | 因素C |
---|---|---|---|
K1 | 3.013 2 | 2.999 6 | 2.831 8 |
K2 | 2.996 0 | 3.071 1 | 3.138 0 |
K3 | 3.096 6 | 3.035 1 | 3.136 0 |
k1 | 1.004 4 | 0.999 9 | 0.943 9 |
k2 | 0.998 7 | 1.023 7 | 1.046 0 |
k3 | 1.032 2 | 1.011 7 | 1.045 3 |
R | 0.033 5 | 0.023 8 | 0.102 1 |
优水平 | A3 | B2 | C2 |
图6
图6
老化箱能量利用系数的影响因素位级图
Fig.6
Level diagram of influencing factors of energy utilization coefficient of aging chamber
表5 各因素对老化箱温度不均匀系数影响的极差分析结果
Table 5
分析指标 | 因素A | 因素B | 因素C |
---|---|---|---|
K1 | 0.236 1 | 0.171 6 | 0.166 8 |
K2 | 0.159 3 | 0.171 7 | 0.172 6 |
K3 | 0.118 5 | 0.170 6 | 0.174 5 |
k1 | 0.078 7 | 0.057 2 | 0.055 6 |
k2 | 0.053 1 | 0.057 2 | 0.057 5 |
k3 | 0.039 5 | 0.056 9 | 0.058 2 |
R | 0.039 2 | 0.000 3 | 0.002 6 |
优水平 | A3 | B3 | C1 |
图7
图7
老化箱温度不均匀系数的影响因素位级图
Fig.7
Level diagram of influencing factors of temperature non-uniformity coefficient of aging chamber
表6 各因素对老化箱相对湿度不均匀系数影响的极差分析结果
Table 6
分析指标 | 因素A | 因素B | 因素C |
---|---|---|---|
K1 | 0.372 4 | 0.322 0 | 0.330 7 |
K2 | 0.318 0 | 0.324 8 | 0.323 5 |
K3 | 0.289 7 | 0.333 3 | 0.325 9 |
k1 | 0.124 1 | 0.107 3 | 0.110 2 |
k2 | 0.106 0 | 0.108 3 | 0.107 8 |
k3 | 0.096 6 | 0.111 1 | 0.108 6 |
R | 0.027 5 | 0.003 8 | 0.002 4 |
优水平 | A3 | B1 | C2 |
图8
图8
老化箱相对湿度不均匀系数的影响因素位级图
Fig.8
Level diagram of influencing factors of relative humidity non-uniformity coefficient of aging chamber
综合以上分析可得:送风角度为影响老化箱能量利用系数的主要因素;送风温度对温度不均匀系数的影响最为显著,且对相对湿度不均匀系数的影响也十分显著,均呈线性相关性,即送风温度的提高能有效提高箱内温度和相对湿度分布的均匀性;送风速度对3个评价指标的影响均不明显。
4.3 正交试验结果方差分析
为进一步分析送风温度、送风速度和送风角度这3个因素对老化箱能量利用系数、温度不均匀系数和相对湿度不均匀系数的影响程度,利用方差分析法对正交试验结果进行处理。定义显著度P,当P<0.01时,表明该因素对评价指标的影响极为显著,用“**”表示;当P<0.05时,表明该因素对评价指标的影响显著,用“*”表示;当P>0.05时,表明该因素对评价指标无显著影响,用“—”表示。同时,引入因素贡献率来量化各因素对不同评价指标的影响程度,其定义为某一因素的离差平方和与总离差平方和之比。F值为正态分布表中的边界值。
表7 各因素对老化箱能量利用系数影响的方差分析结果
Table 7
方差来源 | 离差平方和 | 自由度 | 均方 | F值 | P值 | 显著性 |
---|---|---|---|---|---|---|
总和 | 0.024 508 103 | 8 | 0.012 254 051 | |||
因素A | 0.001 930 196 | 2 | 0.000 965 098 | 1.881 572 | 0.347 033 | - |
因素B | 0.000 852 056 | 2 | 0.000 426 028 | 0.830 591 | 0.546 272 | - |
因素C | 0.020 700 009 | 2 | 0.010 350 004 | 20.178 550 | 0.047 218 | * |
误差 | 0.001 025 842 | 2 | 0.000 512 921 |
图9
图9
各因素对老化箱能量利用系数的贡献率
Fig.9
Contribution rate of each factor to energy utilization coefficient of aging chamber
表8 各因素对老化箱温度不均匀系数影响的方差分析结果
Table 8
方差来源 | 离差平方和 | 自由度 | 均方 | F值 | P值 | 显著性 |
---|---|---|---|---|---|---|
总和 | 0.002 390 041 | 8 | 0.001 195 019 | |||
因素A | 0.002 376 960 | 2 | 0.001 188 480 | 1 128.303 797 | 0.000 886 | ** |
因素B | 0.000 000 247 | 2 | 0.000 000 123 | 0.117 089 | 0.895 184 | - |
因素C | 0.000 010 727 | 2 | 0.000 005 363 | 5.091 772 | 0.164 156 | - |
误差 | 0.000 002 107 | 2 | 0.000 001 053 |
图10
图10
各因素对老化箱温度不均匀系数的贡献率
Fig.10
Contribution rate of each factor to temperature non-uniformity coefficient of aging chamber
表9 各因素对老化箱相对湿度不均匀系数影响的方差分析结果
Table 9
方差来源 | 离差平方和 | 自由度 | 均方 | F值 | P值 | 显著性 |
---|---|---|---|---|---|---|
总和 | 0.001 243 421 | 8 | 0.000 621 709 | |||
因素A | 0.001 177 727 | 2 | 0.000 588 863 | 35.002 774 | 0.027 776 | * |
因素B | 0.000 023 087 | 2 | 0.000 011 543 | 0.686 150 | 0.593 067 | - |
因素C | 0.000 008 960 | 2 | 0.000 004 480 | 0.266 297 | 0.789 704 | - |
误差 | 0.000 033 647 | 2 | 0.000 016 823 |
图11
图11
各因素对老化箱相对湿度不均匀系数的贡献率
Fig.11
Contribution rate of each factor to relative humidity non-uniformity coefficient of aging chamber
5 实验验证
为验证上述仿真结果的可靠性与有效性,对老化箱内部的温度和湿度进行测量,并对比实验结果与仿真结果。
5.1 实验方案
图12
图12
老化箱实物图及箱内测点布置示意
Fig.12
Physical drawing of aging chamber and layout of measuring points in chamber
1)为保证测量分布全面,分别在老化箱工作室内进行3个水平面的布点测量。第1个水平面设置在距最底层样品架上方10 mm处,第2个水平面设置在工作室内的几何中心处,第3个水平面设置在距离顶面1/10工作室高度的位置处。
2)根据GB/T 10586—2006中的要求,共布置13个温度测点(I,II,…,XIII)和4个相对湿度测点(甲、乙、丙、丁),测点分布在3个平面上,中心测点位于工作室几何中心处,其余测点位于距离工作室壁面1/10边长处。
温度和湿度的测量步骤如下:实验开始时送风口持续输送温度为80 ℃和90 ℃、相对湿度为80%和90%的气流,使老化箱内部的温度和相对湿度达到规定值并稳定30 min;安装在相应测点处的温湿度传感器每间隔1 min测量一组温湿度数据并存储在数据采集器中,共测量30组数据。
5.2 实验结果
当送风温度分别为80 ℃和90 ℃时,老化箱内各测点温度的实测值与仿真值的对比如表10所示。由表10可知,老化箱内温度实测值与仿真值一致,两者之间的最大误差为2.3 ℃,误差在合理范围内,由此验证了仿真结果的可靠性与有效性。当送风温度为80 ℃时,上层4个测点的最大温度差为1.7 ℃,中层5个测点的最大温度差为1.5 ℃,下层4个测点的最大温度差为0.9 ℃;当送风温度为90 ℃时,上层4个测点的最大温度差为1.4 ℃,中层5个测点的最大温度差为1.1 ℃,下层4个测点的最大温度差为0.4 ℃。通过计算可得,当送风温度分别为80 ℃和90 ℃时,老化箱内的平均温度分别为78.0 ℃和87.7 ℃,箱内温度与实际送风温度有一定差距,这是因为老化箱存在热耗散现象,热量通过壁面或未密封完全的空隙向外界传递热量。
表10 老化箱内部温度的实测值与仿真值对比 (℃)
Table 10
测点 | 送风温度为80 ℃ | 送风温度为90 ℃ | ||
---|---|---|---|---|
实测值 | 仿真值 | 实测值 | 仿真值 | |
I | 78.6 | 76.7 | 86.9 | 86.5 |
II | 77.3 | 77.8 | 87.0 | 87.0 |
III | 76.9 | 76.1 | 88.3 | 87.1 |
IV | 78.5 | 78.6 | 87.9 | 86.7 |
V | 77.2 | 78.1 | 88.1 | 85.8 |
VI | 78.7 | 77.8 | 87.1 | 86.7 |
VII | 77.6 | 78.2 | 88.2 | 86.8 |
VIII | 77.8 | 78.0 | 87.7 | 87.2 |
IX | 78.1 | 77.5 | 87.2 | 85.8 |
X | 78.8 | 76.7 | 88.1 | 86.5 |
XI | 78.4 | 78.9 | 87.8 | 87.0 |
XII | 77.9 | 78.7 | 88.2 | 88.1 |
XIII | 78.5 | 77.8 | 87.9 | 86.7 |
表11 老化箱内部相对湿度的实测值与仿真值对比 (%)
Table 11
测点 | 相对湿度为80% | 相对湿度为90% | ||
---|---|---|---|---|
实测值 | 仿真值 | 实测值 | 仿真值 | |
甲 | 80.3 | 81.6 | 89.9 | 89.4 |
乙 | 80.5 | 80.5 | 90.0 | 90.0 |
丙 | 80.0 | 81.0 | 90.0 | 90.0 |
丁 | 80.2 | 81.2 | 90.1 | 90.4 |
6 结 论
本文选择送风温度、送风速度、送风角度这3个因素设计了三因素三水平正交试验,研究了各个因素对老化箱内部温度均匀性、相对湿度均匀性和速度均匀性的影响,相较于单因素研究更具综合性和实际意义。
通过对正交试验结果进行极差分析和方差分析可知,对于能量利用系数,各因素的影响次序为送风角度、送风温度、送风速度,即送风角度的影响最为显著。对于温度均匀性,各因素的影响次序为送风温度、送风角度、送风速度,即送风温度的影响最为显著,其对温度不均匀系数的贡献率高达99.45%。对于相对湿度不均匀系数,各因素的影响次序为送风温度、送风速度、送风角度,即送风温度的影响最为显著。综上,针对3个不同的评价指标,最优送风方案分别为:送风温度为90 ℃、送风速度为10 m/s、送风角度为0°,送风温度为90 ℃、送风速度为12 m/s、送风角度为-10°,以及送风温度为90 ℃、送风速度为8 m/s、送风角度为0°。
实验验证结果表明,老化箱内温度、湿度的仿真值与实测值之间的差异较小,其中温度的最大误差为2.3 ℃,相对湿度的最大误差为1.3%,误差均在合理范围内,由此验证了仿真结果的可靠性与有效性。
参考文献
机械式自动垂直钻具执行机构内部流场规律研究
[J].DOI:10.3785/j.issn.1006-754X.2022.00.044 [本文引用: 1]
Research on internal flow field law of mechanical automatic vertical drilling tool actuator
[J].DOI:10.3785/j.issn.1006-754X.2022.00.044 [本文引用: 1]
热泵烘房结构及参数优化仿真设计
[J].DOI:10.6041/j.issn.1000-1298.2020.S1.055 [本文引用: 1]
Optimal simulation design of structure and parameter in heat pump drying room
[J].DOI:10.6041/j.issn.1000-1298.2020.S1.055 [本文引用: 1]
基于CFD的蜂箱调温过程数值模拟及实验研究
[D].
Numerical simulation and experimental research of beehive temperature adjustment process based on CFD
[D].
基于CFD的汽车车身烘房节能研究
[J].DOI:10.3969/j.issn.1000-4998.2010.10.025 [本文引用: 1]
Research on energy saving of oven for automobile body based on CFD
[J].DOI:10.3969/j.issn.1000-4998.2010.10.025 [本文引用: 1]
矿井通风井巷热交换规律数值试验研究
[D].
Numerical test research on the law of mine ventilation roadway heat exchange
[D].
基于CFD的密集式烤房气流均匀性研究
[J].DOI:10.3969/j.issn.1673-3142.2021.12.027 [本文引用: 1]
Research on air flow uniformity of intensive barn based on CFD
[J].DOI:10.3969/j.issn.1673-3142.2021.12.027 [本文引用: 1]
Analysis and optimization of air distribution and ventilation performance in a generator hall using an innovative attached air supply mode
[J].
可压缩流体的瓦斯抽采管路流阻计算方法研究
[J].煤,DOI:10.3969/j.issn.1005-2798.2016.04.003 [本文引用: 1]
Study on calculation method of gas drainage pipeline flow resistance of compressible fluid
[J].DOI:10.3969/j.issn.1005-2798.2016.04.003 [本文引用: 1]
基于强迫对流的热风干燥烘房送风速度及温度优选
[J].DOI:10.11975/j.issn.1002-6819.2022.z.005 [本文引用: 1]
Air
velocity and temperature optimization of hot air drying room based on forced convection[J].DOI:10.11975/j.issn.1002-6819.2022.z.005 [本文引用: 1]
化学分析实验室换气次数数值模拟分析
[J].DOI:10.3969/j.issn.1003-0344.2021.07.018 [本文引用: 1]
Research on numerical simulation of air exchange rate in chemical analysis laboratory
[J].DOI:10.3969/j.issn.1003-0344.2021.07.018 [本文引用: 1]
高炉风口回旋区燃烧特性的数值模拟研究
[D].DOI:10.1016/j.fuel.2021.122490 [本文引用: 1]
Numerical study of combustion characteristics in the raceway of an ironmaking blast furnace
[D].DOI:10.1016/j.fuel.2021.122490 [本文引用: 1]
基于机械虹膜机构的新型摆动泵及其流场分析
[J].DOI:10.3969/j.issn.1001-4551.2023.10.017 [本文引用: 1]
Flow characteristics of new swing pump based on mechanical iris mechanism
[J].DOI:10.3969/j.issn.1001-4551.2023.10.017 [本文引用: 1]
高低温环境试验箱循环风机优化设计
[J].DOI:10.3969/j.issn.1004-7204.2022.06.041 [本文引用: 1]
Optimization of blower in enviornmental test chamber
[J].DOI:10.3969/j.issn.1004-7204.2022.06.041 [本文引用: 1]
双轴型热空气交换试验箱的研制
[J].DOI:10.3969/j.issn.1671-5381.2022.4.hccllhyyy202204035 [本文引用: 1]
Development of biaxial hot air exchange test chamber
[J].DOI:10.3969/j.issn.1671-5381.2022.4.hccllhyyy202204035 [本文引用: 1]
基于温度和速度均匀性的侧送风烘房设计及仿真
[J].DOI:10.11975/j.issn.1002-6819.2021.19.003 [本文引用: 1]
Design and simulation of side air supply drying room based on temperature and velocity homogeneity
[J].DOI:10.11975/j.issn.1002-6819.2021.19.003 [本文引用: 1]
侧送风气流组织设计计算
[J].DOI:10.3969/j.issn.1673-0232.2013.08.238 [本文引用: 1]
The design and calculation of side air supply airflow organization
[J].DOI:10.3969/j.issn.1673-0232.2013.08.238 [本文引用: 1]
Computational fluid dynamic (CFD) investigation of thermal uniformity in a thermal cycling based calibration chamber for MEMS
[J].
热风干燥烘房气流速度场均匀性优化
[J].DOI:10.3969/j.issn.1005-2895.2019.01.016 [本文引用: 1]
Optimization of airflow velocity field uniformity in hot air drying kiln
[J].DOI:10.3969/j.issn.1005-2895.2019.01.016 [本文引用: 1]
单侧供热汽车烘干房气流组织数值模拟研究
[J].DOI:10.3969/j.issn.1000-2324.2021.01.015 [本文引用: 1]
Study on numerical simulation for airflow distribution in automotive baking room with a single side heating
[J].DOI:10.3969/j.issn.1000-2324.2021.01.015 [本文引用: 1]
/
〈 |
|
〉 |
