低气压下翅片管换热器空气侧传质特性实验研究
2021-10-12苑云潇柳建华张嘉文
苑云潇 张 良 柳建华 张嘉文 刘 旗
(1 上海理工大学制冷技术研究所 上海 200093;2 空调国际上海有限公司 上海 200001)
翅片管换热器作为空气热湿处理系统的重要设备,得到了广泛应用。进入二十一世纪,人类在高海拔地区活动范围的拓展以及对城市高空区域的利用使涉及翅片管换热器的空气热湿处理设备开始在各种高海拔地区普及应用,极大地拓宽了翅片管换热器工作环境的气压范围。
伴随着环境气压的改变,湿空气物性和空气流场结构均会发生显著变化:随着环境气压降低,湿空气密度线性下降,湿空气承受水蒸气的能力增加,湿空气的汽化潜热变大。而空气物性的变化及空气流场结构的变化会进一步影响翅片管换热器空气侧的传质性能,因此研究低环境气压下翅片管换热器空气侧的传质特性对于翅片管换热器在高海拔地区应用范围的拓展具有重要意义。
目前人们对于换热器空气侧传质特性的研究主要集中在常压和湿工况下影响参数的分析以及传质特性预测模型的建立[1-3]方面,其中对影响参数的研究主要包括翅片间距[4-6]、换热器管排数[7-9]、空气侧流速[4-5,10-11]及入口空气相对湿度[12-14]等。而关于低气压环境下换热器空气侧传质特性的研究较少:Liu Qi等[15]对翅片管换热器在低环境气压下的传热传质特性进行了实验研究,研究环境压力为40~101.3 kPa,讨论了环境压力的变化对换热器空气侧传热因子、传质因子及刘易斯因子的影响,结果表明环境压力降低时,传热传质性能均下降。Jia Runze等[16]对板翅式换热器在低环境气压和湿工况环境下空气侧传热和传质特性进行了性能实验和模拟,研究气压范围为61~101 kPa,实验结果显示,与常压相比低气压环境下板翅式换热器的潜热和质量传递过程均明显减弱,模拟结果显示低压环境下换热器性能的下降,主要与传质过程中浓度边界层的剧烈变化有关。两者所研究换热器参数单一,实验工况少,对低气压下影响参数的研究不足。因此,有必要研究相关参数对低气压下翅片管换热器空气侧传质特性的影响,为低气压环境下换热器的设计以及传质特性预测模型的建立提供一定的数据基础。
本文实验研究了两个具有不同翅片间距的翅片管换热器在不同入口空气工况下空气侧传质特性,分析环境气压、迎面风速、翅片间距及入口空气相对湿度对换热器空气侧传质特性的影响。
1 实验系统
图1所示为低气压实验系统,主要包括空气循环系统和水循环系统。空气循环系统主体为一个不锈钢板构造的封闭式循环风道,包括横截面为外径377 mm的圆形循环段以及横截面为500 mm×500 mm的矩形实验段,主要功能为:内部设置风机为空气循环提供动力并调节空气侧流速,采用真空泵和气压调节阀调节空气循环系统内的环境气压,设置电加热器和加湿器调节实验用换热器空气侧入口空气干湿球温度和相对湿度,设置传感器对空气相关参数进行测量采集。水循环系统主要包括冷水机组、恒温水箱及水泵,为实验用换热器提供一定流量的恒温冷媒水。
1加湿器;2电加热器;3真空泵;4进气阀;5循环风机;6实验换热器;7整流孔板;8标准喷嘴;9恒温水箱;10冷水机组。图1 实验系统Fig.1 Experimental system
空气循环系统内环境气压的调节由真空泵和气压调节阀协同控制,利用真空泵将环境气压快速抽至设定工况气压值附近,然后调节气压调节阀,实现系统内气压在设定工况值的精确控制。空气经过实验换热器后被冷却除湿,设置电加热器和加湿器补偿该部分的温、湿度变化,同时通过电加热器和加湿器控制调节实验换热器入口截面处空气的温、湿度在工况设计值,相对湿度的测量精度为±0.8%。为保证换热器入口处气流均匀性,在换热器前设计孔板对空气进行整流。空气侧流速采用喷嘴进行测量,测量段包括整流孔板、取压孔和标准喷嘴,安装和测试过程符合ASHRAE 41.2标准[17]。
实验研究的换热器几何参数如表1所示,实验工况如表2所示。
表1 翅片管换热器几何参数Tab.1 Geometry parameters of tested heat exchanger
表2 实验工况Tab.2 Experimental working conditions
2 实验数据处理
实验工况涉及翅片管换热器的析湿工况,因此采用基于焓差的Threlkeld[18]方法进行数据处理,并采用Colburn传质因子jm对翅片管换热器空气侧传质特性进行描述。重要参数计算过程如下。
1)实验中换热器两侧流体总传热量采用两侧流体换热量的平均值表示:
Qa=Ga(ia,in-ia,out)
(1)
Qw=Gwcp,w(Tw,out-Tw,in)
(2)
(3)
2)基于焓差的两侧流体总传热量:
Qavg=Uo,wA0Δilm
(4)
(5)
在翅片管换热器的换热过程中,热阻主要来源于空气侧和水侧的对流换热热阻,因此Uo,w可简化为式(6)进行计算:
(6)
对于式(6)中的管内水侧对流换热系数hw,采用Gnielinski[19]公式计算:
(7)
fw=(1.58lnRew-3.28)-2
(8)
(9)
对于式(6)中的3个系数b、bw,m、bw,t为饱和湿空气焓值曲线在对应温度时的斜率。可由下式进行计算:
(10)
(11)
(12)
(13)
3)空气侧显热换热量:
Qs=haA0η0(Ta,m-Ta,sur)
=Gacp,a(Ta,in-Ta,out)
(14)
翅片总效率η0:
(15)
对于全湿工况下的翅片效率η,采用S.Y.Liang等[20]提出的翅片效率计算模型进行计算,对于部分湿工况下的翅片效率,采用Ma Xiaokui等[21]提出的翅片效率计算模型计算,可得到全湿工况和部分湿工况下的翅片效率。
在式(14)中只有传热系数和翅片效率两个未知数,结合求出的翅片效率,即可得到传热系数。实际求解过程中,对式(14)以及翅片效率计算公式进行迭代计算,从而获得空气侧传热系数ha,以便进行空气侧传质系数的计算。
4)全湿工况下,Threlkeld提出的空气侧传质系数求解方程为:
(16)
采用W.Pirompugd等[22]提出的全湿工况下的传质系数改进计算式:
(17)
(18)
(19)
翅片表面平均温度对应的饱和湿空气焓值im,sur:
(20)
5)部分湿工况下,空气侧传质系数计算式为:
(21)
对式(17)、式(21)进行计算可得空气侧传质系数hm。
无量纲数Colburn传质因子jm:
(22)
计算得到传质因子jm表征空气侧传质特性。
根据Moffat[23]法对实验参数进行误差分析,结果如表3所示。
表3 实验参数精度表Tab.3 Accuracy of experimental parameters
3 实验结果分析
3.1 实验系统精度验证
利用实验台在常压下的实验结果对实验台精度进行验证,以保证实验数据的可靠性。目前对翅片管换热器空气侧对流传质的预测模型研究较少,其中可用于平翅片管换热器的模型主要有W.Pirompugd等[1]换热器空气侧传质因子关联式:
(23)
(24)
将常压下的实验数据代入关联式中进行对比,结果如图2所示。传质因子关联式预测值与实验值的最大相对偏差为21.7%,平均相对偏差为7.4%,有97.9%的预测值与实验值的偏差在±20%以内。结果表明,本实验装置在常压下传质因子的实验结果与预测模型较为吻合,实验装置精度较高,实验数据可靠。
图2 实验台精度验证Fig.2 Accuracy verification of test bed
3.2 环境气压的影响
实验中,在改变环境压力时,通过加湿器向湿空气中喷入水蒸气保持湿空气的相对湿度不变。湿空气的含湿量随着环境压力的变化如图3所示。由图3可知,湿空气的含湿量随着环境压力的下降而增大。
图3 湿空气的含湿量随着环境压力的变化Fig.3 The variation of moisture content of wet air with ambient pressure
图4所示为换热器管排数为两排,换热器入口空气相对湿度为50%时,传质因子随换热器迎面风速、环境气压及翅片间距的变化。由图4可知,相同风速条件下,传质因子随着环境压力的降低而增大。研究工况下,环境压力从40 kPa升至100 kPa时,空气侧传质因子下降18.2%~23.6%,且在风速较小时,传质因子受气压变化的影响更大。分析其原因,在保持入口空气干球温度和入口空气相对湿度不变的情况下,环境气压降低,湿空气的含湿量不断升高,从而导致传质速率增加,传质因子增大。
图4 不同环境气压下传质因子随换热器迎面风速的变化Fig.4 Variation of jm with air velocity under different ambient pressure
3.3 迎面风速的影响
由图4可知,在实验工况下,随着换热器迎面风速的增大,传质因子逐渐减小,变化趋势与Wang Chichuan等[24]和N.H.Kim等[9]在常压下的实验结果趋势相符。研究工况下,迎面风速从0.5 m/s增至4 m/s,传质因子减小81.3%~83.2%。且在较低环境气压中,传质因子随迎面风速的变化更显著。原因在于,迎面风速的增大虽然可以增大空气侧对流传质系数,但空气侧速度场、温度场及浓度场的协同性下降,导致空气侧对流传质系数增长趋势不显著,而由于传质因子表征的是空气单位质量流量时的传质能力,因此空气侧传质因子反而下降。
3.4 翅片间距的影响
由图4可知,环境压力为40、60、80、100 kPa时,对于相同的迎面风速,随着翅片间距的增大,换热器空气侧传质因子均减小,且该减小趋势随着迎面风速的增加逐渐减小,与W.Pirompugd等[22]在常压下的实验结果趋势相符。在迎面风速为0.5 m/s时,翅片间距从2.5 mm增至3 mm,换热器空气侧传质因子降低4.7%~15.9%。该变化趋势出现的主要原因在于,迎面风速减小时,更小的翅片间距有利于换热器空气侧整个流动区域保持稳定层流状态,可以抑制换热管后部涡流区域的影响范围。此外,当翅片间距增大时,冷凝水在翅片间管外壁形成的冷凝水膜会越来越厚,从而减弱了空气侧的传质特性。
3.5 入口相对湿度的影响
图5所示为换热器管排数为两排,换热器翅片间距为2.5 mm时,传质因子随换热器迎面风速、环境气压及换热器入口空气相对湿度的变化。
图5 不同环境气压下入口空气相对湿度对空气侧传质因子的影响Fig.5 Effect of inlet air relative humidity on jm under different ambient pressure
由图5可知,在环境压力为40、60、80、100 kPa时,对于相同的迎面风速,随着入口空气相对湿度的增大,换热器空气侧传质因子均减小,该趋势与马小魁等[12]在常压下的研究结果相符,且该减小趋势随迎面风速的减小更加显著。分析其原因,随着入口空气相对湿度的增大,换热器外表面的冷凝水也越来越多,导致换热器表面冷凝水膜的厚度增大,减弱了换热器空气侧的传质特性。
图5中,迎面风速为0.5 m/s时,当环境气压为100 kPa,相对湿度由90%减至50%,传质因子增加45.4%;当环境气压为80 kPa,相对湿度由90%减至50%,传质因子增加46.7%;当环境气压为60 kPa,相对湿度由90%减至50%,传质因子增加53.3%;当环境气压为40 kPa,相对湿度由90%减至50%,传质因子增加62.8%。主要原因是在研究的低环境气压工况时,环境气压降低使传质因子增大,而相对湿度减小也使传质因子增大,两者对传质因子的影响叠加,使得随着环境气压的降低,相对湿度的变化对空气侧传质因子的影响更加显著。
4 结论
本文实验研究了低气压下翅片管换热器空气侧传质特性,分析了环境气压、换热器迎面风速、翅片间距以及入口空气相对湿度对换热器空气侧传质特性的影响,在所研究工况下,得到如下结论:
1)换热器空气侧传质因子随着环境压力的升高而下降。环境压力从40 kPa升至100 kPa时,空气侧传质因子下降18.2%~23.6%。
2)换热器空气侧传质因子随着换热器迎面风速的增大而减小。换热器迎面风速从0.5 m/s增至4 m/s,传质因子减小81.3%~83.2%。
3)换热器空气侧传质因子随着翅片间距的增大而减小。在迎面风速为0.5 m/s时,翅片间距从2.5 mm增至3 mm,换热器空气侧传质因子降低4.7%~15.9%。
4)换热器空气侧传质因子随着相对湿度的减小而增大。在迎面风速为0.5 m/s时,当相对湿度由90%减至50%,传质因子增加45.4%~62.8%。
5)在较低迎面风速和较低环境气压下,翅片间距和换热器入口空气相对湿度对换热器空气侧传质特性的影响更加显著。
本文受上海市部分地方院校能力建设专项计划(16060502600)资助。(The project was supported by the Capacity Building Plan for some Non-military Universities and Colleges of Shanghai Scientific Committee (No.16060502600).)
符号说明
D——换热管管径,mm
Fp——翅片间距,mm
P——管间距,mm
N——管排数
T——温度,K
u——速度,m/s
umax——最窄面风速,m/s
p——压力,Pa
jm——传质因子
Q——换热量,W
Qavg——换热器两侧流体换热量平均值,W
Qs——空气侧显热换热量,W
G——质量流率,kg/(m2·s)
i——焓值,kJ/kg
is,in——入口水温对应的饱和空气焓值,kJ/kg
is,out——出口水温对应的饱和空气焓值,kJ/kg
ia,sur——翅片表面温度对应饱和空气焓值,kJ/kg
ia,m——空气平均焓值,kJ/kg
is,m——水侧平均温度对应的饱和空气焓值,kJ/kg
ifg——饱和水蒸气焓值,kJ/kg
cp——定压比热容,J/(kg·K)
Uo,w——基于焓差的总传热系数,kg/(m2·s)
A0——换热器空气侧总换热面积,m2
At——换热管外表面积,m2
Af——翅片表面面积,m2
b——换热管壁温与管内平均水温间饱和空气焓值曲线斜率,J/(kg·K)
bw,m——翅片表面液膜平均温度对应的饱和空气焓值曲线斜率,J/(kg·K)
bw,t——换热管表面液膜平均温度对应饱和空气焓值曲线斜率,J/(kg·K)
hw——水侧传热系数,W/(m2·K)
ha——空气侧显热对流传热系数,W/(m2·K)
hm——空气侧传质系数,kg/(m2·s)
Lef——刘易斯因子
f——摩擦系数
Re——雷诺数
Pr——普朗特数
Sca——施密特数
w——绝对含湿量,g/kg
wm,sur——翅片表面饱和空气绝对含湿量,g/kg
n2——雷诺数指数
Δilm——对数平均焓差,kJ/kg
λ——导热系数,W/(m·K)
ρ——密度,kg/m3
μ——动力黏度,N·s/m2
η——翅片效率
η0——整体翅片效率
ε——系数
δ——翅片厚度,mm
下标
a——空气
w——冷却水
out——出口
in——入口
m——平均
o——外侧
i——内侧
t——横向
l——纵向
sur——翅片表面
s——饱和