0 引言

换热设备广泛应用于化工、石油、能源、冶金、建筑等各个领域.据统计，在现代化工企业中40%~60% 的设备属于换热设备.随着生产规模的扩大和能源问题的日益突出，制造体积小、耗能低、流动阻力小、高效紧凑的换热器将具有重要意义.而强化传热是实现这一目标的根本途径^[1].

翅片管作为强化传热的有效元件在换热器中的使用越来越广泛，特别是在两侧换热系数相差较大（例如一侧为液体，一侧为气体），其强化传热效果更突出.翅片管的散热性能一方面依赖于其翅片高度、翅片厚度和翅片间距^[2]；另一方面依赖于翅片管周围流体的流动状态.利用有限元分析软件ANSYS中的流体传热分析模块（FLOTRAN CFD）对翅片管换热过程进行二维数值模拟，研究其热交换过程和传热效果，通过对正交实验结果的直观分析，优化其结构参数，进而指导其加工.

1 实验研究方法

目前，翅片管多用于空调的冷凝器、蒸发器的换热元件，本文以常用于冷凝器中的低翅片管为研究对象，根据其实际的传热方式，选定气-气的两组分对流传热模型，即模拟翅片管的实际工况，管内为高温R22制冷剂气体，管外为冷空气，以逆向对流的方式进行对流换热.

选取对翅片管换热影响较大的3个结构参数作为研究的出发点，分别为翅片间距（P），翅片厚度（T）及翅片高度（H），每个因素取3种水平，忽略交互作用，采用正交表L₉（3⁴）来安排实验^[3]；翅片管其他主要结构参数如下：内径取φ16 mm，底壁厚（基管厚度）为0.8 mm.由于翅片管具有轴对称的结构，故此处选用沿管轴向横截面的上半部分作为研究对象，如图 1所示 （图 1中箭头表示为流动方向）：

其中A₁为管内流通的R22 制冷剂气体，A₂为翅片管，A₃为管外流通的冷空气，在ANSYS中建立上述模型并进行数值模拟分析；通过对正交实验结果的直观分析，得出各结构因素对翅片管换热、压降等性能影响的主次关系，进而确立其综合性能评价指标，从而完成对低翅翅片管结构的优化^[4].

2 边界条件及载荷施加

管外空气外边界L₁确定为无滑移条件，即此处V_x=0、V_y=0，管轴线L₂为对称线，故此处V_y=0，为了简化计算模型，假定换热过程中的流体的物性参数为常值，具体的载荷及物性参数如表 1、表 2所示^[5].

3 数值模拟过程 3.1 流动状态模拟

在ANSYS中对翅片管的传热过程进行数值模拟时，先是模拟其流动状态，求解出流动方程，进而可以得到动能、速度、压力等因素的分布场，在流体分析的基础上，加载热载荷，进行温度分析，从而得出温度分布情况^[6].对翅片管管内外流体流动状态进行数值分析时，由其雷诺系数（管内：Re =2.423×10⁵ ； 管外：Re =7042）可判断管内外流体均为湍流流动，由其马赫数可知流体为不可压缩流体^[7-8]，故求解翅片管管内外流体流动状态时，将流体设置为不可压缩紊流状态，进行稳态求解.

求解流动状态时，需对计算过程的收敛性及稳定性进行监测，如图 2所示，为求解计算的收敛监测曲线，由此可知，整个求解过程是收敛的，求解结果可靠^[9].

3.2 流动模拟结果

图 3及图 4分别示出了光管与某型号翅片管换热过程中管内外流体的压力分布，从图 3可以看出光管管内外流体的压力分布比较有规律，流体在流动过程中，压力逐渐减小；在保持进口速度不变的情况下，光管的压降比翅片管的压降要小很多，即翅片管上的翅片在增加散热面积的同时也增大了流体的流动阻力，从而使位于管外有翅片一侧的流体的压力成倍增大.由图 4可看 出，翅片管管外入口端流体压力变化剧烈，特别是在靠近第一片翅片的地方，流体压力变化异常明显，总体趋势是由管壁向边界从大到小线性变化，靠近边界处压力趋于稳定.

翅片管管内外流体的速度分布情况类似于其压力分布情况，管内流体流动较规律，处于层流状态，靠近管壁的地方，流速小，由管壁向管轴方向速度不断变大，管轴处的速度近似等于流体入口速度.由图 5可知管外流体的流速分布较复杂，入口处第一翅片周围的流体流速变化剧烈，是典型的紊流状态，翅片间流体流动缓慢，流速小，但已出现明显的回旋运动，即所谓的有涡流动^[10]，正是因为这种有涡流动的存在，增强了翅片管的对流换热效率，并且使得翅片管外侧不易结垢，起到一定的自除垢作用.在由翅片往边界的径向上，流速有以下变化特点：小→大→小.随着流动的发展，沿管轴方向，翅片周围的流体流速逐渐减小.

3.3 温度模拟结果

图 6所示为管内出口端管内外流体的温度分布情况，由图 6可知，管外流体由于刚流入，温度还很低，温度基本上为初始温度，翅片管的温度也近乎等于管外流体的温度.图 7为管内出口端流体沿管径方向的温度变化曲线，横坐标表示在管径方向上，距离管壁的距离长度，纵坐标表示对应的温度.由图 6及图 7可以看出管内流体的温度分布呈现如下特点：在径向方向上，靠近管壁处流体温度变化剧烈，越靠近管轴，温度变化越缓慢，在轴向方向上，随着流动的发展，流体的低温区域越来越大，平均温度逐渐降低.

4 实验结果分析 4.1 正交模拟实验的直观分析

利用上述方法，可以分别得到各组实验的管内流体出口平均温度及管外压降、速度等指标的数据^[11]，表 3示出了实验方案及温降的结果分析，表 3中，因素A 为翅片间距，取3种水平：1 mm、1.5 mm、2 mm，因素B 为翅片厚度，分别取0.25 mm、0.3 mm、0.35 mm3种水平，因素C为翅片高度，分别取0.8 mm、1.4 mm、2.0 mm3种水平.

分析表 3可以得出如下结论：对翅片管换热效能影响最大的因素是A（翅片间距），次之是C（翅高），影响最小的是B（翅厚）.利用表 3同样能得到各因素与流体压力、速度间的关系，据此可以将因素水平作为横坐标，以它的试验指标的平均值k_i 为纵坐标，得到因素与指标的关系图 ——趋势图^[3]，如图 8~图 10所示.

从图 8可以看出，翅片管管内流体的温降受翅片间距的影响最大，翅片间距在1~2 mm范围内，温降随翅距的增加而增大，翅距超过2 mm时，温降随翅距的增大而减小，（此实验中温降的大小代表翅片管换热效能的高低）；流体温降随翅片厚度的增大而减小，随翅片高度的增大而增大.由图 9可知，管外流体的压降受翅片高度的影响最明显，压降随翅高的增加而成倍数变大（此实验压降大小代表翅片管工作噪音大小）^[12]，翅片间距对压降的影响有以下特点：翅距在1~2 mm范围内，压降随之增大而增大，在2~2.5 mm压降随翅距的增大而减小；压降 基本不受翅片厚度的影响.管外流体出口流速受各因素的影响如图 10所示：翅片高度对出口流速的影响最明显，与压降类似，流速随翅高的增大而迅速增大；翅片间距在1~2.5 mm范围增大时，流速先是随之减小，然后增大；翅片厚度同样对流速也影响很小.

4.2 模拟结果分析

由文献^[13]可知：翅片间距增大，翅片管的换热面积减小，使平均传热系数减小，即翅片管的换热效能随翅距的增大而减弱，这与本文的结论（即翅距在1~2 mm范围时，翅片管的换热效能随翅距的增大而加强）相矛盾，查阅相关资料获知，翅片周围的流体存在边界层，当翅片间距小于或等于两倍边界层厚度时，相邻翅片间的流体近似处于静止状态，传热形式由对流传热变成了流体的导热，那么就大大降低了传热效果^[14-15].当翅距在1~2 mm变化时，翅距越大，翅片间流体流动趋于活跃，传热方式由导热变成了对流传热，传热系数增大，当翅距超过2 mm时，管外传热方式固定为对流传热，换热面积成为影响换热效能的主要因素，翅距变大，单位长度翅片管的翅片换热面积变小，换热效能降低^[16].

4.3 换热特性实验研究

为了验证模拟结果的可靠性，本文还采取相应的热交换器性能实验的实验数据作对比，实验装置示意图见图 11.

图 11 冷凝器性能测试实验装置示意图

对规格为内径φ16 mm、管壁厚0.8 mm和长度100 mm的光管及文中模拟所用型号的翅片管（管长为100 mm），在相同入口温度（冷流体283 K，热流体314 K），相同入口流量（2.17 m3/h），相同流动形式（逆流）下通过上述实验装置获得管内R22制冷剂的出口温度，见表 4.

经计算，管内流体出口温度的实验值与模拟值之间的误差，光管误差为3.6 %，翅片管的误差在10 %~15 %范围，造成误差的原因可能是管内外流体有少许发生了液化相变及生成有限元模型时，网格划分的不够精确.分析表 4的实验值可以进一步确认翅片管较光管有更好的强化传热效果，且翅片管的结构参数对其换热效果有明显的影响.

5 结论

（1）管内外流体初始条件相同的情况下，翅片管的换热性能优于光管，以实验中A3B2C3规格的翅片管为例，在100 mm流程内，其管内流体的出口温降较之光管相对提高0.96 K.

（2）管内外流体的流动：管内流体流动规律、简单，翅片管管外流体流动情况比光管的复杂的多，翅片高度对流体的流动影响显著，其次是翅片间距.

（3）翅片间距对翅片管换热性能影响明显，在1~2 mm范围内，翅距越大，换热性能越强；翅片高度越大，换热效果越好，但压降也随之成倍数增长，翅片厚度越小，换热越好，综合换热与压降两方面，得到最优化方案为A3B1C2，即翅片间距为2 mm，翅片厚度为0.25 mm，翅片高度为1.4 mm.