第 9 卷第 1 期
2010 年 3 月
热 科 学 与 技 术
Journal of Thermal Science and Technology
Vol .9 No .1
Mar .2010
文章编号 :1671‐8097(2010)01‐0031‐05 DOI :10 .3969/
j
.issn .1671‐8097 .2010 .01 .006
同 轴 径 向 热 管 换 热 器 壳 程 流 场 数 值 模 拟
毛 建 丰 , 涂 福 炳
倡
, 许 欣 , 周 孑 民 , 张 岭
( 中南大学 能源科学与工程学院 ,湖南 长沙 410083 )
摘要 :
基于多孔介质模型和分布阻力方法 ,引入 Al‐sanea 和 Taborek 两种阻力关系式模拟同轴径向热管换
热器壳程的流场 。结果表明 :换热器壳程静压沿烟气流动方向呈线性分布 ;随入口烟气速度的增加 ,换热器阻
力损失增大 、压降增大 ;且随入口烟气速度的增加 ,压降增加的速率增大 。
关键词 :同轴径向热管 ;换热器 ;流场 ;数值模拟
中图分类号 :TK172 .4 ;TK11
+
4 文献标识码 :A
收稿日期 :2009‐10‐25 ; 修回日期 :2009‐12‐09 .
基金项目 :湖南省科技创新资助项目(2007CK3070) .
作者简介 :毛建丰 (1985‐) ,男 ,湖南双峰人 ,硕士生 ,主要研究方向为热工设备的数值仿真与优化 ;
涂福炳(1966‐) ,男 ,湖南株洲人 ,博士 ,副教授 ,主要从事热工过程 、炉窑节能及计算机仿真优化研究 .
E‐mail :tfbone@ 163 .com
0 引 言
换热器在工业烟气余热回收中得到了广泛的
应用 ,增强对流换热 ,减少流动阻力是提高换热器
综合性能的有效途径
[1]
。同轴径向热管换热器由
于其具有高导热 、等温性好 、结构紧凑等诸多优点
而广受青睐 。由于热管换热涉及相变过程 ,数值方
法直接模拟相变过程有很大的难度 ,研究者在数
值求解热管换热器的问题时提出了许多的方法来
避免这个问题 。Patankar 和 Spalding 提出分布阻
力和多孔介质孔隙率概念对管壳式换热器进行数
值模拟研究
[2]
;王定标等建立了紧凑换热器的多
孔模型 ,并提出了建立此模型的基本前提和方
法
[3]
。这些模型和方法很好地解决了壳管式热管
换热器数值模拟问题 。为给同轴径向热管换热器
的优化设计提供依据 ,本文采用文献[2‐3]的模型
和方法对其壳程流场进行数值模拟 ,求得同轴径
向热管换热器壳程压降的关系 ,并通过与实测值
的比较验证数值模拟结果的正确性 。
1 数值模拟
1 .1 数学物理模型
本文研究对象为国内某厂气态悬浮炉烟气余
热回收装置上选用的长宽高分别为 2 460 mm
×
2 380 mm
×
1 870 mm 的同轴径向热管换热器 。热
管分布见图 1 ,其中热管外径为 51 mm ,热管横向
间距为 108 mm ,纵向间距为 120 mm 。换热器由
322 根热管组成 ,管束叉排 ,每根热管外表面沿轴
线方向分布有 400 个左右的等间距 、等厚度环翅
片 。翅片高度为 26 .5 mm ,间距为 6 mm ,厚度为 1
mm 。烟气从换热器上部流入 ,下部流出 ,在换热
器内部通过热管与冷却水进行热交换 。
基于 Darcy‐Bricnkman‐forchheimer 的多 孔
介质模型假设 ,直角坐标系中以体积多孔度表示
的控制方程组可以表示为
[4]
抄
抄 x
i
(
βρ
u
i
φ
) =
抄
抄 x
i
(
β
Γ
抄
φ
抄 x
i
) + S
φ
(1)
式中 :
β
为体积多孔度 ,
ρ
为流体密度 ,u
i
为速度 ,
φ
为一般变量 ,
Γ
为有效黏度 ,S
φ
为源项
-
β
R
i
-
抄(
βρ
)
抄 x
i
,R
i
为坐标轴方向分布阻力 。
当流体在管束中流动时 ,会受到管束施加的
形状阻力和壁面摩擦阻力 ,总称为流体所受的分
布阻力 。热管外径 、表面粗糙度以及翅片的高度 、
厚度 、分布方式等对流体受到的阻力有很大的影
响 。本文在固定的热管形状尺寸和排列方式的前