贾润泽等
摘要:对设计的桥式翅片换热器空气侧的传热和阻力性能进行了实验研究,将大量的实验数据进行了线性回归,得出了在实验雷诺数范围内传热和阻力性能关联式及特性曲线.对比可知,在相同泵功情况下桥式翅片换热器比相同尺寸的平直翅片换热器具有更高的传热性能.同时,对以上两种翅片空气侧的温度场和速度场进行了数值模拟,并利用场协同原理对模拟结果进行了分析.分析结果表明,桥式翅片换热器具有更高传热性能的根本原因在于翅片的桥式布置能有效地改善翅片温度场和速度场的协同性.
关键词:桥式翅片;换热器;实验研究;数值模拟;场协同
中图分类号:TM925.12 文献标识码:A
多年来,人们对散热器翅片强化传热开展了大量研究,并提出了很多有效的强化传热的方法,其中将翅片开缝就是其中一种[1].这是因为在翅片上开缝可以破坏边界层的生长,从而抑制传热系数沿流动方向上的降低.同时将翅片的开缝位置加工成拱桥状,形成所谓的桥式开缝翅片,其在不减少换热面积的情况下进一步增大气流的扰动,这样有可能会进一步提高开缝翅片空气侧的传热性能,但同时也会带来负面效应,比如增大空气阻力,增加功率消耗等.
目前对于开缝翅片的传热及阻力特性的研究比较多[2-5].李惠珍等[6]对2排X型双向开缝翅片进行了实验研究和模拟计算,并将实验结果与同尺寸的单向开缝翅片和平直翅片进行了比较,研究结果表明开缝翅片的传热性能远高于平直翅片,且X型的双向开缝翅片的性能要好于单向开缝翅片.Yun J Y等[7]的研究显示翅片节距、开缝翅片的缝宽、缝长、缝高、开缝数及缝的分布等都不同程度地影响其流动和传热性能.屈治国等[8]在对平直开缝翅片的传热特性进行了数值模拟的基础上,进一步应用场协同理论对计算结果进行了分析,计算结果显示在速度场和温度场协同比较差的区域开缝要比在场协同比较好的区域开缝对传热强化更有效.综合上面分析发现,对开缝翅片的已有研究主要集中在纯开缝形式的翅片上,而对于桥式开缝翅片的流动及传热特性的研究还有待进一步的开展.
本文采用实验研究和数值模拟的方法,对桥式翅片换热器和同尺寸的平直翅片换热器进行了研究,并利用场协同理论对计算结果进行分析,为桥式翅片换热器进一步的优化设计提供了依据.
1实验装置及元件参数
实验测试是在标准的风洞实验台上进行的,测试采用热平衡法.实验台由风洞、整流栅、恒温水箱、电加热器、循环水泵、水流量测量、水温控制、压差测量、阀门、换热器、温度传感器、空气流量测量、空气阻力测量、风量调节装置、风机、变频器等组成.实验台工作原理如图1所示.
实验过程中,通过控制恒温水浴出口阀门的开度来控制流经换热器的水量,翅片侧空气的流量则通过改变风机的转速(2 000~5 200 r/min ,每隔400 r/min 取一个转速)来加以控制,因此测试共有9个工况点.为提高测量的准确性,每次测量均在工况稳定40 min后进行,各工况的热平衡偏差控制在5% 以内.如图1所示,各点的温度依据空气及热水位置的不同,分别采用水银温度计和热阻式温度传感器两种方式进行测量.在空气侧,空气的流量采用毕托管配倾斜式微压计进行测量计算,同时,在换热器前后的风管上设置静压测嘴测量空气通过换热器的流通压降.在水侧,通过恒温水浴加热使换热器的入口水温稳定在60 ℃,同时在换热器进出口处设置压力表对热水通过换热器的前后的压降进行测量.为保证测量的准确性,水流量的测量采用称重法.
2实验误差分析及数据整理
实验测试前为验证实验台的准确性,用光管进行了阻力和传热特性实验,所得的实验结果与经典的光管布劳修斯公式f=0.316 4Re -0.25及公认的SiederTate 公式计算式Nu=0.027Re 0.8Pr 1/3(μ/μw)0.14的计算结果进行了对比,试验结果最大偏差均控制在±6%以内,表明本实验台具有足够的精确度.
实验时管外翅片侧为空气强制对流换热,管内侧为热水与管内壁的强制对流换热,再加上管自身的导热,因此翅片管的热阻共有3部分组成.应用热阻分离法[9],空气侧的对流换热系数表达式如下:
3实验结果
利用最小二乘法对实验数据进行拟合整理,得出桥式翅片管的传热和阻力特性实验关联式如下:
传热系数关联式为:Nu=6.173Re 0.324 4, 适用范围为1.59×102< Re<8.2×102 ,最大拟合偏差±1.94%.
气侧阻力系数关联式为:f=0.289 5Re 0.192 5,适用范围为1.59×102< Re<8.2×102 ,最大拟合偏差 ±3.51%.
平直翅片管和桥式翅片管的传热和阻力特性试验曲线如图3所示.从图中可知,随着Re增大,翅片空气侧的对流换热系数h也随之增大,风阻系数f也增大;桥式翅片管的对流换热能力高于平直翅片,但风阻系数也略高.在实验的Re范围内,桥式翅片的传热能力比平直翅片平均高52%,风阻系数平均高25%.通过比较可知,桥式翅片是一种有效强化空气的传热方式.
4数值模拟
4.1数值模拟方法
利用有限元分析软件Fluent对以上2种换热器的翅片单元进行仿真,由于桥式翅片管换热器在几何结构上沿横向和纵向具有周期性和对称性,因此可以取相邻翅片间的对称区域为计算单元,如图4所示.利用三维造型软件solidworks建立的实体模型导入Gambit进行网格划分,Gambit中实体建立遵循点—线—面—体的关系,通过此顺序将夹层空间建立实体,根据流道空间形状,将实体划分为六面体网格,有利于计算、收敛和提高精度.该计算模型为三维稳态层流不可压缩流动,因此算法采用simple算法,模型选用laminar.在计算区域的选取上,沿流动方向上进口取为1.5 L,出口取为5 L以使出口回流影响最小,入口边界设置为velocityinlet,出口边界为pressureoutlet,流场的四周边界设
为symmetry,圆管与流场和翅片的交界面设为wall.计算中假设管壁温度为固定温度,翅片的温度
受到空气流动和翅片管导热的影响,需要计算确定,因而这是一个对流与导热的耦合问题,对此问题的边界条件处理方法可参见文献[8].至此,完成在Gambit中网格划分和定义边界.
4.2数值计算结果与台架实验结果的对比
桥式翅片管通过两种方式所得的换热系数随迎面风速的变化曲线如图5所示.从趋势上看,数值计算结果与实验结果在总体上保持一致,翅片侧的对流换热系数在实验风速范围内皆随迎面风速的增加而增加.从定量上进行分析,u0=0.8 m/s是个分水岭,当u0>0.8 m/s时,实验值高于计算值,当迎面风速u0<0.8 m/s时,实验值低于计算值;最大误差出现在u0=2.2 m/s时,误差为9.8%,说明数值计算结果和实验结果符合良好.
4.3数值计算结果分析
本文同时对桥式翅片和相同尺寸的平直翅片进行了数值模拟,计算得到的翅片单元温度场和速度场如图6,图7所示.沿流动方向上对流换热系数及温度梯度和速度梯度夹角随迎面风速的变化曲线如图8,图9所示.
从图6,图7可以清楚地看到,翅片形式设计成桥式后使得边界层在翅片台阶处被破坏,有效地减薄了厚度,因而对流传热显著强化.还可以从场协同理论[10]得到进一步理解,场协同原理指出,在相同的温度和速度的边界条件下,使速度场和温度梯度的夹角变小(即提高场的协同性),可有效提高换热效果.由图6,图7可知在翅片的进口段,速度与温度梯度的方向几乎一致,因而速度与温度梯度的协同性很好,换热强烈,但在流经翅片间隙时,两矢量的夹角增加,速度与温度梯度的协同性明显变差,因而局部对流换热系数下降,如图8中的曲线所示.
根据速度和温度梯度夹角的计算公式[11]:
由上式计算得到的沿流场方向上协同角随迎面风速的变化曲线如图9所示.由图9可知,平直翅片在整个流场的平均协同角要高于桥式翅片,而且随着迎面风速的增加有增大的趋势.因此可以看出,桥式翅片相对于平直翅片具有更好的温度场和速度场的协同性,这也是图8中桥式翅片空气侧对流换热系数沿流动方向下降较慢的原因.
就阻力损失而言,如图10所示,桥式翅片使流动阻力增加,因而桥式翅片的进出口压力损失高于相应的平直翅片,但在迎面风速0.6~2.8 m/s的范围内,桥式翅片的压力损失的增加幅度小于翅片空气侧对流换热系数的增加幅度,综合评价得出在相同泵功下桥式翅片的传热性能比平直翅片高出约50.3%.
5结论
1)实验测试了桥式翅片换热器的传热性能和阻力性能,得出了该类型翅片换热器的传热和阻力性能关联式.
2)与平直翅片实验结果对比,桥式翅片换热器的换热性能更优,在风阻系数增加25%的情况下,平均传热系数高出约52%.
3)通过数值模拟,得出了相同尺寸的桥式和平直翅片计算单元区域的温度场和速度场分布情况,及流动方向上对流换热系数沿程变化曲线和不同迎面风速下的压力损失及协同角的变化曲线.并用场协同原理对仿真结果进行分析,分析结果表明翅片桥式布置能有效改善翅片温度场和速度场的协同性,因而相同泵功下大幅度增加了翅片的传热性能.
参考文献
[1]胡俊伟,丁国良. 开缝翅片压降和换热特性的数值模拟[J]. 上海交通大学学报,2004,38(10): 1639-1642.
HU Junwei, DING Guoliang. The numerical simulation of pressure drop performance and heat transfer performance of strip fin[J]. Journal of Shanghai Jiaotong University,2004, 38(10): 1639-1642.(In Chinese)
[2]金巍巍,屈治国,张超超,等. 空调蒸发器用管翅式换热器开缝翅片的数值设计[J]. 工程热物理学报,2006, 27(4): 688-690.
JING Weiwei, QU Zhiguo, ZHANG Chaochao,et al. Numerical design of slotted fin of tube fin heat exchanger used in evaporator of air conditioning[J]. Journal of Engineering Thermophysics, 2006, 27(4): 688-690. (In Chinese)
[3]李红智, 罗毓珊, 王海军, 等. 增压空冷器用开缝翅片的结构优化设计[J]. 西安交通大学学报, 2008, 42(11): 1355-1360.
LI Hongzhi,LUO Yushan,WANG Haijun,et al. Configuration optimization design for slotted fin surfaces in pressurized air cooler[J]. Journal of Xi'an Jiaotong University,2008,42(11): 1355-1360. (In Chinese)
[4]李红智, 罗毓珊, 王海军, 等. 空间波动参数在开缝翅片结构优化设计中的应用[J]. 工程热物理学报, 2009,30(3):445-448.
LI Hongzhi,LUO Yushan,WANG Haijun, et al. Space fluctuation parameters and its application in configuration optimization design for slotted fin surfaces[J]. Journal of Engineering Thermophysics, 2009,30(3):445-448. (In Chinese)
[5]熊伟, 罗毓珊, 王海军, 等. 二种开缝翅片特性的试验及数值模拟比较研究[J]. 化学工程, 2009, 37(1): 11-14.
XIONG Wei,LUO Yushan,WANG Haijun,et al. Experiment and numerical simulation for characteristics of two slotted fins[J]. Chemical Engineering, 2009, 37(1): 11-14. (In Chinese)
[6]李惠珍, 屈治国, 程永攀, 等. 开缝翅片流动和传热性能的实验研究及数值模拟[J]. 西安交通大学学报, 2005, 39(3): 229-232.
LI Huizhen, QU Zhiguo, CHENG Yongpan, et al. experimental and numerical study on heat transfer and fluid flow characteristics of slotted finandtube heat transfer surfaces[J]. Journal of Xi'an Jiaotong University, 2005, 39(3): 229-232. (In Chinese)
[7]YUN J Y, LEE K S. Influence of design parameters on the heat transfer and flow friction characteristics of the heat exchanger with slit fins[J]. International Journal of Heat and Mass Transfer, 2000, 43(14): 2529-2539.
[8]屈治国, 何雅玲, 陶文铨. 平直开缝翅片传热特性的三维数值模拟及场协同原理分析[J]. 工程热物理学报, 2003, 24(5):825-827.
QU Zhiguo, HE Yaling, TAO Wenquan. 3d numerical simulation on heat transfer performance of slit fin surfaces and analysis with field synergy principle[J]. Journal of Engineering Thermophysics, 2003, 24(5):825-827. (In Chinese)
[9]杨世铭,陶文铨. 传热学 [J]. 第3版.北京:高等教育出版, 1998: 207-212.
YANG Shiming,TAO Wenquan. Heat transfer [M]. 3rd ed.Beijing: Higher Education Press, 1998: 207-212. (In Chinese)
[10]GUO Z Y, LI D Y, WANG B X. A novel concept for convective heat transfer enhancement[J]. International Journal of Heat and Mass Transfer, 1998, 41(14): 2221-2225.
[11]过增元.对流换热的物理机制及控制:速度场与热流场的协同[J]. 科学通报,2000,45(19):2118-2122.