中心被堵波纹管流动及换热特性的数值模拟研究.pdf

中国工程热物理学会 传热传质学 第十二届学术会议 编号：063195 中心被堵波纹管流动及换热特性 ‘的数值模拟· 王小佳 曾敏 张宏伟 陶文铨 西安交通大学能源与动力工程学院，动力I[程多相流国家重点实验室，陕西，西安，710049 Tel：029 jtu．edu．cn E—mail：wqtao@mail．x 摘要：本文在二维圆柱轴对称坐标系下采用微分方程法生成适体网格对中心被堵波纹管的流动及换热 特性进行了数值模拟。模拟时堵住的固体区和流体区分别生成两套网格，搭接成一套网格气固耦合整 场求解。将所开发程序应用于圆管中心被堵模型的计算中，同解析解进行对比，结果吻合良好。计算 中心被堵波纹管的结果表明：在相同流量下，中心堵管的存在起到强化换热作用，同时引起阻力增大。 其换热得到强化的机理可通过场协同原理进行解释。 关键词：中心被堵波纹管，圆柱轴对称适体坐标，气固耦合，数值模拟 1．引言 迄今为止，管壳式换热器广泛的应用于各种工业领域。光管由于低成本，易加工等 优点在管壳式换热器中的应用是最广的，但存在换热面积较小、管壁易结垢、换热能力 较差等缺点。作为直壁圆管的管壳式换热器的取代品，波纹管式换热器具有传热系数高、 换热量大、有效抑制污垢IlJ等优点，目前在热电、化工、轻工、城市的集中供暖系统等 众多领域中发挥出越来越大的作用。字波等心1通过对纵向内波纹翅片管加中心堵管的方 法强化换热，受此思想的启发，笔者将中心堵管应用到波纹管中，采用二维圆柱轴对称 坐标系下的椭圆形微分方程生成的适体坐标，对中心 被堵波纹管内的流动及换热特性进行数值模拟。相关 的实验研究将在后继工作中进行。 2．物理模型、控制方程及边界条件 2．1物理模型的描述’ 本文计算所采用的物理模型的几何尺寸如下：波 面型线用多项式(1)表示。中心堵管为直壁圆管，半径尺J。=O．2×咫，因参与换热， 其物性参数按铜取定。流动工质采用物性参数恒定的空气。 Y=I I／1000(1) 图2给出了整个计算模型的儿何示意图，其中的万表示进出口各延K的一个?仁控制 容积。 基金项日：国家自然科学幕金项H 1089 ut也t 图2整个计算模型的几何小恿图 。 2．2控制方程及边界条件． 对于动量及能量方程可归纳成统一的表达式： 一 (2) V·(py≯)=V·(F，V≯)+墨 ， 对于材方程，式(2)中r，为刁·墨为一印／苏；对于1，方程，r，为77，墨为 一劲lar-刀vlr2；对于能量方程，F．为r／／Pr，邑为0。’ 相关边界条件如下： 入1：3边界条件： 流体区u=“肿1，=0，t=t。=25．2uC； 固体区Ⅳ=0，v=0，Ot／Ox=0； (3) 出口边界条件： 流体区速度和温度均采用局部单项化； 固体区“=0，v=0。OtlOx=0： (4) 壁面边界条件： “=0，v=0，q。=25．038W／mz； (5) 轴线边界条件： u=0，v=0，Ot／Or=0；