错排环布圆管翅片式换热板芯的计算研究.pdfVIP

中国工程热物理学会 传热传质学 2004年学术会议 编号：043017 错排环布圆管翅片式换热板芯 的计算研究 樊菊芳 于永兴 王良壁 (兰州交通大学机电工程学院， 兰州 730070) 联系电话：0931．4956556E—mail：lbwang@mail．1习m．cn 摘要 本文运用数值模拟的方法对错排环布圆管翅片式换热板芯进行了研究t对几种结构的速度场、 温度场、局部Nu数分布、横向平均Nu数分布在300Re4800的范围内进行比较分析·为该结构与 直排圆管翅片式换热板芯作比较提供了理论依据，对此结构的进一步研究奠定基础· 关键词换热板芯；数值模拟； 同位网格；有限差分 1．引言 目前对换热器板芯的研究主要是针对直排式进行的，对环形布置换热板芯的研究未 见报道．而这种结构在实际工程中已被采用，如电力机车变压器油散热器。为了对这种 结构进行更深入的研究，本文对错排环布圆管翅片式换热板芯进行了模拟分析。 2．换热器的数学物理模型 2．1物理模型及控制方程 本文计算对象为：翅片间距，p=2．5mm、管数mn。=60、管排数^‘。-=3，4，5 的错排环布圆管翅片式换热板芯．其物理模型平面图如图1所示。对于Nro,=4的物 理模型其三维立体图如图2所示。为了计算的方便，从图1中截取0=6。的一部 分作为计算用的物理模型，如图3(a)所 示。考虑到结构的对称性及流动的特点， 从图3(a)中截取一半作为本文数值模拟的 计算区域，如图3(b)所示。 为了便于分析，我们假定：(1)流体 “ 是常物性的：(2)流动和换热处于稳态层 流非周期性的初始段；(3)体积力和耗散 项可忽略。 以下以张量的形式给出流动与换热的控制方 程： 连续性方程：O-：-(p“，J=o 图l错排环布圆管翅片式换热 板芯的物理模型 运动方程：i(p％№)=1_m—o)一÷(k=1，，3)． 运动方程：善(川音嘻)_毒舻№2，，．咚 钱啦 呶、 能量方程：虽(pcur／=三rkOT一，． 图4为使用无限插值法生成的网格示意图。通过复 杂的数学变换将直角坐标中的控制方程变换成适体坐标 F的控锘4方程，表达式为『1】： 图2：错排环布圆管翅片式 意(p∽+盖如”+兹如聊-o 鼍c一㈨+畚c。㈨+毒c。聊，=毒陪一期+击眵一刳+善眵一期+峨 2．2边界条件 入口边界条件为： u／2Uf,in，T=％5 alr —o 出口边界条件为： 塑：堡：塑：o： 出 融 融 固体壁面边界条件为： (a) 。·、 “=v=w=0，L=cD扎 对称面边界条件为： (b) 图3：计算区域示意图 誓=o，筹-o，詈吨=。册 铆 咖 3．数值求解 将适体坐标系中的控制方程，在控 制容积上积分，导出计算域面上三维通 用的微分方程离散形式聊。计算采用基 于同位网格系统的SIMPLE算法[21，并 采用Rhje．Chow[3]的方法计算界面上逆 变速度。经过网格考核，加密或稀疏网 格对数值计算结果影响很小。因此所得 的数值解是网格独立的解。 图4计算网格示意图 4．结果及分析 4．1翅片局部换热特性分析 图5为7≥2．5mm、^■=4、 N一--60时．不同Re数下翅片局部 换热系数的分布。由图可知： 在圆管前驻点位置，翅片局部 换热系数