| 多头螺旋管换热过程的三维数值模拟 |
多头螺旋管作为近几年来兴起的新型高效换热元件,已然成为传热行业研究的热点,相比对称型良好的传统换热管材,如光管、翅片管等,螺旋管内的流体在进行管内外对流换热的过程中,涉及一个更为复杂的旋转紊流的三维流动[1].考虑到多头螺旋管复杂、特殊的结构及换热过程流体复杂的流动状态,文中选用数值模拟的方法对其换热过程进行研究,通过模拟计算,得到螺旋管内流体紊流流动过程中的温度场与流场分布,为其强化换热的研究提供理论依据.
1 换热过程的有限元模型实体模型是建立有限元模型的基础,考虑到多头螺旋管非轴对称的结构特性,文中选用了三维实体建模软件Pro/E 来建立多头螺旋管的模型,然后以特定格式导入到ANSYS 中进行数值模拟分析.表 1示出了实验管材的基本结构参数.
| 表1 多头螺旋管结构参数 |
![]() |
| 点击放大 |
根据表 1 给出的管材结构参数,建立多头螺旋管的实体模型,选用当量直径为内径20 mm, 管长为200 mm 的光管,由于文中研究的是管内紊流流体的流场及温度场变化,所以在建立模型时只需建立流体的实体模型,这使得模型得到了科学的简化.图 1及图 2 分别为螺旋管换热过程的几何模型及有限元模型[2].
![]() |
| 图 1 几何模型 |
![]() |
| 图 2 有限元模型 |
2 边界条件及载荷施加
固体壁面处确定为无滑移条件,为了简化计算模型,假定换热过程,液体的物性参数为常值,管外的高温空气以恒定的壁面热流率形式施加[3-5],具体的载荷及流体物性参数如表 2 所示.
| 表2 载荷及流体物性参数 |
![]() |
| 点击放大 |
3 模拟结果及分析 3.1 温度场结果分析
图 3及图 4 分别示出了光管及六头螺旋管换热过程中的温度分布,从图 3、图 4 不难看出,低温的水在流动过程中,温度逐渐升高,越靠近管壁,温度越高,但由于流体存在轴向的速度,各个部分的温度并不相同,沿着流动方向温度分布呈现峰状,低温流体区域逐步递减,并且相对光管,螺旋管的峰形更尖锐,这也从侧面说明螺旋管内流体的换热更剧烈,效率更高.临近出口端,流体的温度分布基本趋向均匀.
![]() |
| 图 3 光管的温度分布剖面图 |
![]() |
| 图 4 螺旋管的温度分布剖面图 |
图 5和图 6 将温度分布以数据的形式反映在数学曲线上,其中横坐标代表从入口端到出口端的管的轴向长度,纵坐标代表对应管长截面上的温度值.设温升的相对提高为ε,Tg为光管中流体从入口端到出口端的温升,Tl为螺旋管内流体入口端与出口端的温差[6].则根据曲线中的数据不难得到:
![]() |
| 图 5 光管内流体的温度变化 |
![]() |
| 图 6 螺旋管内流体的温度变化 |
|
相对光管,螺旋管内流体的温升相对提高了近2.6 倍.这从客观数据上证明了螺旋管的强化传热作用.
3.2 流场结果分析图 7及图 8 表征了光管和螺旋管中流体运动过程的速度分布,从图 7、图 8 可以看出,沿管材的径向呈现速度由近壁处向中心逐渐增大,且近壁处速度升幅剧烈,中间部位速度变化较为平缓,速度几近峰值.相对光管,螺旋管内紊流的流体,紊流的无序特征变得明显,特别是在近壁处,流体已出现明显的流体回旋运动,即所谓的有涡流动[7],正是由于这种三维有涡流动的存在,使得螺旋管内壁不易结垢,从而起到自除垢的作用.
![]() |
| 图 7 光管的速度分布剖面图 |
![]() |
| 图 8 螺旋管的速度分布剖面图 |
图 9和图 10 以数值曲线的方式反映了出口流体沿径向上合速度的变化,不难发现无论光管还是螺旋管,其管内流体的合速度分布均沿管的轴向成中心对称,相比之下,螺旋管内流体近壁处速度变化更剧烈.根据曲线上的速度值可得出光管出口流体的平均速度约为0.6 m/s, 螺旋管出口流体的平均速度约为0.3 m/s.这主要是因为螺旋管旋转的凸筋结构增加了流体流动的阻力,但是在流体流动过程中,螺旋管正是由于其具有这种特殊的螺旋凸筋结构,破坏了流体传热的边界层,才使得其强化传热强度得到了提高[7-14],同时,根据过增元院士的场协同理论可知,速度矢量与热流矢量的夹角越小,速度场与热流场的协同性越好,则其强化传热性能越好.相对光管,螺旋管内流体靠近管壁处其速度表现为轴向、径向和周向上三维流动,这种有涡流动大大改善了速度场和热流场的协同效应,从而起到明显的强化传热作用[15-19].
![]() |
| 图 9 光管内流体速度的变化 |
![]() |
| 图 10 螺旋管内流体速度的变化 |
4 结论
(1)在等热流的情况下,多头螺旋管的换热性能优于光管,以文中研究的六头螺旋管为例,其流体出口温升较之基准管相对提高近2.6 倍.
(2)流体流动过程中,流体流速的变化与到管壁的距离成反比,近壁处速度变化剧烈.相比光管内流体,多头螺旋管内流体流动阻力更大,其出口速度相对降低了50 %.
(3)多头螺旋管特有凸筋结构的存在,一方面破坏了流体传热的边界层,并且改善了速度场与热流场的协同性,使得强化传热效率得到提高,另一方面,有涡流动对管壁的冲刷作用,使得多头螺旋管管壁不易结垢.
| [1] | 林宗虎, 汪军, 李瑞阳. 强化传热技术[M]. 北京: 化学工业出版社 , 2006. |
| [2] | 王咏梅, 康显丽, 张丽萍. Pro/ENGINEER Wildfire 4.0 中文版基础教程[J]. 北京:清华大学出版社, 2008. |
| [3] | 彭洁. 螺旋槽管换热过程的三维数值模拟[D]. 秦皇岛:燕山大学,2006. |
| [4] | 朱冬生, 郭新超, 刘庆亮. 扭曲管内传热及流动特性数值模拟[J]. 流体机械, 2012, 40(2): 63–67. |
| [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. |
| [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. |
| [13] | 邹华生, 钟理, 伍钦. 流体力学与传热[M]. 广州: 华南理工大学出版社 , 2004. |
| [14] | 周瑞, 何海澜, 白明爽. 波纹管内湍流传热数值模拟[J]. 石油化工设备, 2011, 40(6): 49–51. |
| [15] | 过增元. 对流换热的物理机制及其控制[J]. 科学通报, 2001, 45(19): 2118–2122. |
| [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. |
| [18] | 车晓梅, 陈伟庆, 赵江, 等. 水口类型对大方坯连铸结晶器内流场和温度场的影响[J]. 有色金属科学与工程, 2012(2): 32–36. |
| [19] | 靳遵龙, 李赵, 董其伍, 等. 螺旋槽管管内强化传热的三维数值模拟[J]. 郑州大学学报:工学版, 2011, 32(6): 5–8. |
2012, Vol. 3













