The 3-D numerical simulation of heat transfer process for multi-start spiral pipe
-
摘要: 随着节能减排的生态发展理念深入人心,如何最大限度地提高能源的利用率,成为时下研究的热点. 文中以多头螺旋管为研究对象, 并以当量直径的光管作为基准管, 运用有限元分析软件ANSYS 对模型的传热过程进行数值模拟,得出了换热过程中温度场与流场的分布状况,而后结合传热学中的边界层理论和协同场理论对结果进行综合分析.结果表明:在等热流的情况下,多头螺旋管的换热性能明显优于光管,相对光管,多头管内流体的出口温升相对提高了2.6 倍,出口速度相对下降50 %;其特殊的内外凸筋结构,是强化传热的主要原因,同时,凸筋造成流体的有涡流动,也使得多头螺旋管管壁不易结垢.Abstract: The utilization of energy has become the current hotspot with the development concept of energy saving. This article studies the 3-d numerical simulation of heat transfer process for multi-start spiral pipe by using the light pipe with equivalent diameter as the benchmark tube. The heat transfer process of models is simulated through the ANSYS, obtaining the distribution of the temperature and flow field in the heat transfer process. The results are analyzed comprehensively by the boundary layer theory and collaborative field theory in heat transfer. The results show: the heat transfer performance of multi-start spiral pipe is significantly better than that of the light pipe with the same heat flux. In comparison with the light pipe, the exit temperature of the fluid for the multi-start spiral pipe is increased to 2.6 times and the exit speed is reduced to 50 %. The main reason of heat transfer enhancement lies in its internal and external rib structure. The walls of multi-start spiral pipe are not liable for scaling for fluid vortex flow due to its rib structure.
-
0 引言
多头螺旋管作为近几年来兴起的新型高效换热元件,已然成为传热行业研究的热点,相比对称型良好的传统换热管材,如光管、翅片管等,螺旋管内的流体在进行管内外对流换热的过程中,涉及一个更为复杂的旋转紊流的三维流动[1].考虑到多头螺旋管复杂、特殊的结构及换热过程流体复杂的流动状态,文中选用数值模拟的方法对其换热过程进行研究,通过模拟计算,得到螺旋管内流体紊流流动过程中的温度场与流场分布,为其强化换热的研究提供理论依据.
1 换热过程的有限元模型
实体模型是建立有限元模型的基础,考虑到多头螺旋管非轴对称的结构特性,文中选用了三维实体建模软件Pro/E 来建立多头螺旋管的模型,然后以特定格式导入到ANSYS 中进行数值模拟分析.表 1示出了实验管材的基本结构参数.
表 1 多头螺旋管结构参数根据表 1 给出的管材结构参数,建立多头螺旋管的实体模型,选用当量直径为内径20 mm, 管长为200 mm 的光管,由于文中研究的是管内紊流流体的流场及温度场变化,所以在建立模型时只需建立流体的实体模型,这使得模型得到了科学的简化.图 1及图 2 分别为螺旋管换热过程的几何模型及有限元模型[2].
2 边界条件及载荷施加
固体壁面处确定为无滑移条件,为了简化计算模型,假定换热过程,液体的物性参数为常值,管外的高温空气以恒定的壁面热流率形式施加[3-5],具体的载荷及流体物性参数如表 2 所示.
表 2 载荷及流体物性参数3 模拟结果及分析
3.1 温度场结果分析
图 3及图 4 分别示出了光管及六头螺旋管换热过程中的温度分布,从图 3、图 4 不难看出,低温的水在流动过程中,温度逐渐升高,越靠近管壁,温度越高,但由于流体存在轴向的速度,各个部分的温度并不相同,沿着流动方向温度分布呈现峰状,低温流体区域逐步递减,并且相对光管,螺旋管的峰形更尖锐,这也从侧面说明螺旋管内流体的换热更剧烈,效率更高.临近出口端,流体的温度分布基本趋向均匀.
图 5和图 6 将温度分布以数据的形式反映在数学曲线上,其中横坐标代表从入口端到出口端的管的轴向长度,纵坐标代表对应管长截面上的温度值.设温升的相对提高为ε,Tg为光管中流体从入口端到出口端的温升,Tl为螺旋管内流体入口端与出口端的温差[6].则根据曲线中的数据不难得到:
相对光管,螺旋管内流体的温升相对提高了近2.6 倍.这从客观数据上证明了螺旋管的强化传热作用.
3.2 流场结果分析
图 7及图 8 表征了光管和螺旋管中流体运动过程的速度分布,从图 7、图 8 可以看出,沿管材的径向呈现速度由近壁处向中心逐渐增大,且近壁处速度升幅剧烈,中间部位速度变化较为平缓,速度几近峰值.相对光管,螺旋管内紊流的流体,紊流的无序特征变得明显,特别是在近壁处,流体已出现明显的流体回旋运动,即所谓的有涡流动[7],正是由于这种三维有涡流动的存在,使得螺旋管内壁不易结垢,从而起到自除垢的作用.
图 9和图 10 以数值曲线的方式反映了出口流体沿径向上合速度的变化,不难发现无论光管还是螺旋管,其管内流体的合速度分布均沿管的轴向成中心对称,相比之下,螺旋管内流体近壁处速度变化更剧烈.根据曲线上的速度值可得出光管出口流体的平均速度约为0.6 m/s, 螺旋管出口流体的平均速度约为0.3 m/s.这主要是因为螺旋管旋转的凸筋结构增加了流体流动的阻力,但是在流体流动过程中,螺旋管正是由于其具有这种特殊的螺旋凸筋结构,破坏了流体传热的边界层,才使得其强化传热强度得到了提高[7-14],同时,根据过增元院士的场协同理论可知,速度矢量与热流矢量的夹角越小,速度场与热流场的协同性越好,则其强化传热性能越好.相对光管,螺旋管内流体靠近管壁处其速度表现为轴向、径向和周向上三维流动,这种有涡流动大大改善了速度场和热流场的协同效应,从而起到明显的强化传热作用[15-19].
4 结论
(1)在等热流的情况下,多头螺旋管的换热性能优于光管,以文中研究的六头螺旋管为例,其流体出口温升较之基准管相对提高近2.6 倍.
(2)流体流动过程中,流体流速的变化与到管壁的距离成反比,近壁处速度变化剧烈.相比光管内流体,多头螺旋管内流体流动阻力更大,其出口速度相对降低了50 %.
(3)多头螺旋管特有凸筋结构的存在,一方面破坏了流体传热的边界层,并且改善了速度场与热流场的协同性,使得强化传热效率得到提高,另一方面,有涡流动对管壁的冲刷作用,使得多头螺旋管管壁不易结垢.
-
表 1 多头螺旋管结构参数
表 2 载荷及流体物性参数
-
[1] 林宗虎,汪军,李瑞阳. 强化传热技术[M]. 北京:化学工业出版社,2006. [2] 王咏梅,康显丽,张丽萍. Pro/ENGINEER Wildfire 4.0 中文版基础教程[M]. 北京:清华大学出版社,2008. [3] 彭洁. 螺旋槽管换热过程的三维数值模拟[D]. 秦皇岛:燕山大学,2006. [4] 朱冬生,郭新超,刘庆亮. 扭曲管内传热及流动特性数值模拟[J]. 流体机械,2012,40(2):63-67. http://www.cnki.com.cn/Article/CJFDTOTAL-LTJX201202013.htm [5] 许洋,党沙沙,胡仁喜. ANSYS11.0/FLOTRAN 流场分析实例指导教程[M]. 北京:机械工业出版社,2009. [6] 明宇. 具有螺旋形插件换热管的流场与温度场分析[D]. 沈阳:东北大学,2008. [7] 翟庆良. 湍流新理论及其应用[M]. 北京:冶金工业出版社,2009. [8] 过增元,黄素逸. 场协同原理与强化传热新技术[M]. 北京:中国电力出版社,2004. [9] Bergles A E.Heat transfer enhancement-encouragement and accom- modation of high heat fluxes[J]. Journals of Heat Transfer,1995,119: 8-19. http://cn.bing.com/academic/profile?id=1987971803&encoded=0&v=paper_preview&mkt=zh-cn
[10] Lisenhard J H V. Review of heat transfer augmentation in turbulent flows[J]. Applied Mechanics Reviews,1998,51(2):B19.
[11] 贾瑞宣. 螺旋槽管强化传热研究[D]. 北京:华北电力大学,2003. [12] Islamg Y. Effect of rounding of protructing edge on convection heat transfer in a converging-diverging channel[J]. International Commu- nication in Heat and Mass Transfer, 2008,35(1):643-647. http://cn.bing.com/academic/profile?id=2020383064&encoded=0&v=paper_preview&mkt=zh-cn
[13] 邹华生,钟理,伍钦. 流体力学与传热[M]. 广州:华南理工大学出版社,2004. [14] 周瑞,何海澜,白明爽. 波纹管内湍流传热数值模拟[J]. 石油化工设备,2011, 40(6):49-51. http://www.cnki.com.cn/Article/CJFDTOTAL-SYSB201106012.htm [15] 过增元. 对流换热的物理机制及其控制[J]. 科学通报,2001, 45(19):2118-2122. http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB200019019.htm [16] Wang S, Guo Z Y, Li Z X. Heat transfer enhancement by using metallic filament insert in channel flow[J]. IJ Heat Mass Transfer , 2001, 44:1373-1378. doi: 10.1016/S0017-9310(00)00173-3
[17] ZHONG Jian-hua, JIANG Li-ming, FENG Kai. The 3-D numerical simulation of heat transfer process of finned copper tube [J]. Advanced Materials Research, 2012, 393-395:412-415. https://www.researchgate.net/publication/282585117_The_3-D_Numerical_Simulation_of_Heat_Transfer_Process_of_Finned_Copper_Tube
[18] 车晓梅,陈伟庆,赵江,等. 水口类型对大方坯连铸结晶器内流场和温度场的影响[J]. 有色金属科学与工程,2012, (2)1:32-36. http://ysjskx.paperopen.com/oa/DArticle.aspx?type=view&id=201201008 [19] 靳遵龙,李赵,董其伍,等. 螺旋槽管管内强化传热的三维数值模拟[J]. 郑州大学学报:工学版,2011,32(6):5-8. http://www.cnki.com.cn/Article/CJFDTOTAL-ZZGY201106003.htm