八边形翼攻角对翅片管空气侧流动与换热特性的影响_赵辉.pdf

第 39 卷第 1 期2023 年 2月上海电力大学学报JournalofShanghaiUniversityofElectricPowerVol39，No 1Feb2023DOI:10 3969/j issn 2096 8299 2023 01 008收稿日期:2022-03-08通信作者简介:阴继翔(1964)，女，博士，副教授。主要研究方向为对流传热的强化及新型换热设备的开发。E-mail:459619332 qq com。基金项目:国家自然科学基金青年基金(51706150);山西省科技重大专项(MD2016-02)。八边形翼攻角对翅片管空气侧流动与换热特性的影响赵辉，阴继翔，李志斌，李娜英(太原理工大学 电气与动力工程学院，山西 太原030024)摘要:在传统翅片管上加装八边形翼涡发生器以强化传热。利用 FLUENT 软件对八边形翼翅片管模型进行了数值模拟，研究了不同雷诺数(e)下八边形翼攻角(10，20，30，40)对翅片管空气侧流动与换热特性的影响。结果表明，八边形翼攻角对管后尾迹区换热影响显著;与平直翅片管相比，加装八边形翼的翅片管空气侧在不同攻角、e 下的平均努塞尔数(Nu)提高了131%431%，相应的阻力系数(f)增加了 6 9%61 1%;通过对综合传热性能评价因子(PEC)的分析发现，在 e 为862 3 735 时，八边形翼攻角为 20的翅片管综合传热性能最佳。关键词:八边形翼;翅片管换热器;涡发生器;数值模拟;综合传热性能中图分类号:TK124文献标志码:A文章编号:2096 8299(2023)01 0047 07Air-side Flow and Heat Transfer Characteristics of Fin-tube withOctagonal Winglet at Different Attack AnglesZHAO Hui，YIN Jixiang，LI Zhibin，LI Naying(School of Electrical and Power Engineering，Taiyuan University of Technology，Taiyuan Shanxi030024，China)Abstract:Octagonal winglet vortex generators are installed on the traditional finned tubes to en-hance heat transfer The numerical simulation of the finned tube model with octagonal winglet iscarried out by using FLUENT software，and the effect of the attack angles(10，20，30，40)ofthe octagonal winglet on the air-side flow and heat transfer of the finned tubes is studied at differ-ent eynolds numbers(e)The results show that the attack angle of the octagonal winglet has asignificant effect on the heat transfer in the wake area behind the tube Compared with the flatfinned tube，the average Nusselt number(Nu)of the air side of the finned tube with octagonalwinglet is increased by 13 1%43 1%at different attack angles and eynolds numbers，and thecorresponding resistance coefficient(f)is increased by 6 9%61 1%By analysizing the PECfactor，the finned tube with the attack angle of 20 for the octagonal winglet has the best compre-hensive heat transfer performance in the range of eynolds number 862 to 3 735Key words:octagonal winglet;finned tube heat exchanger;vortex generator;numerical simu-lation;comprehensive heat transfer performance上海电力大学学报2022 年翅片管换热器广泛应用于能源动力、石油化工和制冷工程等诸多领域。通常管外流体为空气，由于空气的热物性导致翅片管换热器热阻主要集中在空气侧，因此如何改善空气侧换热是科技人员长期关注的焦点问题。涡发生器技术作为强化换热的一种有效手段，其原理是主流流体流经涡发生器在其后部产生二次流，并在主流流体推动作用下形成三维螺旋运动，增强流体扰动的同时减薄甚至破坏流体边界层，提高流体换热系数，从而达到强化换热的目的。近年来，研究者引入涡发生器强化换热器换热，对涡发生器的形状、摆放位置和几何尺寸进行了大量的数值模拟和实验研究。何雅玲等人 1 采用数值模拟的方法分析了矩形翼对翅片管式换热器流动和传热的影响，以空气侧平均换热系数、压降为参数对矩形翼的攻角、数目、安装位置进行了优化;苏石川等人 2 采用三维数值模拟方法对加装三角翼的翅片管换热器的流动和传热特性进行了研究，结果发现涡发生器的存在改善了流道中速度场和温度场的协同性;王成刚等人 3 研究了涡发生器攻角对翅片传热性能的影响，发现涡发生器攻角为45的翅片平均努塞尔数(Nu)最大，且随着攻角增加，阻力系数随之增加;胡万玲等人 4 研究了三角形、矩形和梯形涡发生器对换热器传热及阻力的影响，发现三角形涡发生器能够使换热器获得最佳综合换热性能;田丽亭等人 5 应用三维数值模拟方法对加装三角翼涡发生器的波纹翅片管换热器的流动换热特性进行了研究，发现三角翼的存在不仅改善了尾迹区的换热，还大大加强了三角翼下游管排壁面的换热;郭磊等人 6 利用 FLUENT 软件对加装梯形翼的翅片进行了数值模拟分析，并与未加装涡发生器的翅片进行对比，结果显示梯形翼的存在有效延长了空气与基管的接触时间，同时空气流过梯形翼后形成了旋转轴与主流方向一致的纵向涡旋，强化了翅片的传热;党伟等人 7 用数值方法分析了曲面矩形涡发生器高度对换热器传热和流动的影响，并指出当涡发生器高度为 1 7 mm 时产生的二次流强度最大;闵春华等人 8 在矩形翼的基础上提出了一种新型的八边形翼，通过实验比较了加装两种翼后矩形通道内的流动和传热性能，发现八边形翼强化换热效果比矩形翼好，且阻力没有明显增加，由此表明八边形翼是一种性能优良的涡发生器。综上所述，加装涡发生器可以增强翅片管的换热性能。目前对于涡发生器的研究主要集中在矩形、三角形和梯形翼上，而对加装八边形翼的翅片管空气侧流动和换热特性研究的较少。因此，本文按常规涡发生器“渐缩式”安装方式将其加装在翅片管基管两侧，并改变八边形翼与主流方向的夹角(即攻角)，应用 FLUENT 软件模拟研究4 种攻角(10、20、30、40)的八边形翼翅片管在 6 种入口风速(1 6 m/s)下的流动和换热特性，并与平直翅片管(未加装涡发生器)进行对比，旨在寻求综合传热性能最优的翅片管结构。1几何模型及网格划分八边形翼翅片管模型计算域如图 1 所示。考虑到翅片结构以及通道内流体流动与换热的对称特性，取上、下翅片厚度的一半以及其通道内的流体作为研究区域;空气沿正 x 方向流入，为了保证入口流速均匀，防止出口受回流影响，将换热器进、出口分别延长翅片间距的 10 倍和 35 倍(由于页面限制，延长段未按实际尺寸画出)。综合上述所有区域为计算域，其范围是:30 x168，0y12 6，0z3 2。换热基管取 3 排并以叉排方式排列。八边形翼采取冲压方式从翅片表面冲出，因此会在翅片上留下与八边形翼面积相同的孔洞。图 1 中，S1为纵向管间距，S2为横向管间距，D 为基管外径，L 为翅片纵向长度，为翅片厚度，Tp为翅片间距。未加装八边形翼的基本翅片管参数见表 1。图 2 为攻角示意图，攻角 为八边形翼竖直平面与主流方向的夹角。攻角作为决定八边形翼位置关系的关键参数，其角度值对翅片管空气侧的流动和换热过程有明显影响。本文旨在寻找使翅片管达到最佳综合传热性能的攻角。八边形翼是在矩形翼的基础上切除其 4 个顶角得来的，4 个顶角大小相等。八边形翼长、高以及顶角的纵切边、横切边分别用 a，d，c，d 表示，其中 a=5 4 mm，b=1 8 mm，c=d=0 5 mm。其结构如图 3 所示。此外，因为翼采用从翅片表面冲压的方式形成，所以其厚度与翅片厚度相同。84赵辉，等:八边形翼攻角对翅片管空气侧流动与换热特性的影响图 1八边形翼翅片管模型计算域表 1基本翅片管参数单位:mmS1S2DLTp2125 29 4630 23图 2攻角示意图 3八边形翼涡发生器结构为了获得高质量的网格，利用 GAMBIT 软件对计算域进行分块网格划分。八边形翼和管壁区域附近空气流动剧烈，温度梯度较大，采用四面体非结构化网格进行划分，并进行适当的加密处理，以保证计算的精度;进、出口延长区采用较粗的六面体结构化网格进行划分，以降低计算成本。2数学模型2 1控制方程假设计算域中为常物性不可压缩空气的三维稳态流动，不考虑翅片、八边形翼和管壁之间的辐射换热，忽略管壁和翅片的接触热阻。采用 k-模型处理湍流流动，速度与压力的耦合采用SIMPLE算法完成，动量方程和能量方程的离散采用二阶迎风格式。相关控制方程见文献 9。2 2边界条件入口采用速度边界条件(1 6 m/s)，温度为293 K;出口为自由出流边界条件;圆管设为恒壁温(350 K)。加装八边形翼时，翅片厚度中剖面采用温度周期性变化边界条件，翅片两侧表面采用绝热边界条件;与空气接触的翅片进出口端面、翅片表面、八边形翼表面以及孔洞壁面为耦合面，其温度由翅片导热以及和空气对流换热耦合求得;孔洞94上海电力大学学报2022 年表面采用周期性边界条件;其余采用对称边界条件。无八边形翼时，翅片厚度中剖面采用绝热边界条件，其余边界条件与加装八边形翼的翅片条件相同。2 3数据处理采用 FLUENT 软件模拟计算，雷诺数 e、定性温度 tf、努塞尔数 Nu、通道内阻力系数 f 和综合性能评价因子(Performance Evaluation Critevia，PEC)10 的定义如下。e=uinDv(1)式中:uin 入口空气的速度，m/s;v 空气的运动黏度，m2/s。tf=tw tm(2)tm=(tw tin)(tw tout)lntw tintw t()out(3)tout=AoututdAAoutudA(4)式中:tf 流体平均温度，K;tw 翅片、八边形翼以及换热基管表面平均温度，K;tm 对数平均温差，K;tin、tout 流体通道进、出口截面平均温度，K;u 空气出口 x 方向的速度，m/s;A 出口截面积，m2。Nu=hD(5)式中:h 翅片、八边形翼以及换热基管表面平均对流换热系数，W/(m2K);空气的导热系数，W/(mK)。f=2pDu2mL(6)式中:p 流道进出口压降，Pa;um 流道最小截面处的平均流速，m/s;空气密度，kg/m3。PEC=NuNuPFfPF()f13(7)式中:NuPF、fPF 相同工况下平直翅片管表面平均努塞尔数及通道内阻力系数。上述式(1)、式(5)、式(6)中，v 取