纵向涡发生器强化传热的研究进展.docxVIP

纵向涡发生器强化传热的研究进展 紧凑换热器广泛应用于空调、空调、乐器冷却等领域。由于气体的传热性能很差，如何加强气体侧的换热是学术界和工程教学的中心主题。垂直流器是传热表面的一种特殊结构。如果水流流经，旋转方向对应于方向的二次流，即垂直流。传统观点，垂直漩涡可以降低水流的热边界层，使其混合并增加，从而加强传热。由于空气动力学设备的结构简单易用，因此很容易在换热大师的翅膀上安装。因此，在过去的十年里，对垂直旋转换向的研究有所增加。请参考文献中关于纵向轴向运动的实验和数值研究方法的讨论。 现有文献的研究有以下不足:首先,没有同时考虑纵向涡发生器的厚度及因加工的原因而在通道壁上产生的小孔;其次,计算通道的进、出口边界条件的处理方法也不同,对于所使用的边界条件的合理性也缺乏论述;最后,关于纵向涡强化对流换热的根本机制究竟是什么,没有统一的观点.一般认为是由于纵向涡破坏了边界层的发展,使流体旋转、增加混合所造成,没有从速度场与温度场协同观点来分析. 本文在考虑纵向涡发生器的厚度及其下面的小孔的情况下,采用延长计算区域的方法,研究Re、涡发生器的冲角、纵向涡发生器的形状对通道平均Nu(Num)和平均阻力系数的影响,并对计算结果应用场协同原理来分析,以进一步揭示纵向涡强化换热的内在机理. 1 数值控制方程 本文所计算的Re限在3 000以下,使流动处于层流状态,同时假定流动过程是稳态的.所用的控制方程包括三维稳态常物性的连续性方程、动量方程、能量方程(无内热源).这些控制方程在文献中已有充分介绍,此处从略.所计算的矩形通道的边界条件说明如下. 1.1 带小翼小翼的通道结构及边界条件 对于矩形通道,选一对矩形小翼作为纵向涡发生器来进行研究.通道长度方向为x,高度方向为y,宽度方向为z.由于对称性取通道的一半作为计算对象(见图1).通道宽B为160 mm,高H为40 mm,长L为400 mm.因为矩形小翼常常是从壁面上冲出来的,所以其厚度就是翅片的厚度.目前,实际应用的翅片换热器的翅片间距和翅片厚度之比大约为10左右.因此,当通道高度为40 mm,矩形小翼的厚度b就取为4 mm.矩形小翼高h为20 mm,即通道高度的一半,长l为40 mm.带小翼的通道结构见图1a,矩形小翼在通道中的位置见图1b.计算区域从通道入口向上游延长100 mm(为通道高度的2.5倍),从通道出口向下游延长600 mm(为通道高度的15倍).为了进行比较,对没有纵向涡发生器的相同通道也要进行计算.边界条件汇总如下. 进口采用纵向速度、温度均匀,横向速度为0的边界条件.出口采用充分发展边界条件.上、下壁面采用无滑移、恒壁温边界条件.前、后对称面上为对称边界条件. 因为这里模拟的通道是多片翅片组成的多个通道中的一个,每片翅片上都冲有一个小孔,所以上、下两壁面上的小孔处取为周期性边界条件.矩形小翼的壁面取为无滑移边界条件,温度与壁面温度相同.计算区域延长段的所有侧面均取为对称边界条件. 1.2 单元稳定性计算方法 本计算采用商业软件FLUENT进行计算.在矩形通道中,矩形小翼与主流成一定夹角,使得计算区域的几何形状变得复杂,所以将矩形小翼附近区域分成12个小块,分别定义各边上的节点数以控制网格的疏密,使翅片壁面附近和矩形小翼附近网格较密.矩形小翼附近的网格分布如图2所示.控制方程离散采用控制容积积分法,对流项离散采用二阶迎风格式,扩散项采用中心差分格式,压力和速度的耦合采用SIMPLEC算法.收敛性判据采用迭代过程中单元不平衡质量余量的绝对值之和与迭代之初前5次的最大余量的绝对值之比小于10-7.满足这个条件后,动量方程、能量方程的相对余量都已小于10-8,且温度场不再发生变化.所计算的介质为空气,Pr为0.71. 2 平均温度计算方法 在数据处理中,无量纲量的特征速度为通道入口平均速度,特征尺寸为通道高度H的2倍.计算通道平均对流换热系数的温差采用壁温与空气进、出口平均温度之差.协同角的定义为速度与温度梯度之间的夹角,这里反映全场协同性的协同角取其体积平均值.f为平均摩擦阻力系数,p为断面平均静压. 2.1 不同冲角的分布 作为比较的基础,首先计算了没有纵向涡发生器的光通道在Re分别为800、1 600、2 400和3 000时的4个算例.根据文献,在涡发生器的冲角β小于45°时,通道内主要产生纵向涡,而β大于45°后,横向涡的成分越来越多,β等于90°时,则全部为横向涡.在纵向涡发生器的尺寸、位置不变的情况下,分别计算了带矩形小翼的通道在矩形小翼的冲角为15°、30°、45°、60°和90°,Re为800、1 600、2 400和3 000时的各种情况下相应的物理量,并把计算结果与光通道(带下标“0”)的相应值进行对比,结果如图3和4所示. 由图3可见,Re增加,Nu也