开缝翅片管翅式换热器翅片侧换热表面换热流动性能的数值模拟.docVIP

开缝翅片管翅式换热器翅片侧换热表面换热流动性能的数值模拟 2008年3月7日 - 工程流体网 如果您对本文所提及的软件、服务等信息感兴趣，请填写表单，如有相关资料，会发送给您，谢谢关注! 赵震 李光熙 西安交通大学热工教研室 摘要 本文对由表面开缝翅片组成的管翅式换热器翅片侧换热表面在层流情况下的换热和阻力性能进行了数值模拟。计算模型网格生成采用商用软件gambit，数值计算采用商用软件FLUENT。对开缝翅片在不同流体进口流速情况下的换热和阻力特性进行了比较，同时得到了有关情况下翅片表面换热和阻力特性的规律性结果。 关键词 开缝翅片,管翅式换热器 1.引言 管翅式换热器常用于制冷、空调行业中有空气或其他液体介质存在的换热场合。这种换热器一般换热较强的介质在管内流动，壳侧热阻比较大，所以常常是传热过程强化的首选目标，通常的做法是在空气侧（壳侧）加翅片，改变翅片的几何形状和几何尺寸。因此，对于管翅式换热器中由不同的翅片组成的换热表面的换热及阻力性能进行研究，具有十分重要的意义。但是采用实验的方法，设备复杂且往往成本较高，因而采用数值计算的方法来模拟其流动和换热情况越来越受到人们的重视[1-3]。而开缝翅片又是各种强化换热翅片中很重要的一种类型[4]，所以对于表面开缝翅片组成的管翅式换热器翅片侧换热表面换热流动性能进行数值模拟有着很重要的意义。本文的目的旨在通过数值模拟给出在一种由开缝翅片组成的管翅式换热器的换热表面的流动和换热规律，对新型翅片的开发和换热器的设计提供帮助。 2.物理问题及数学描写 对于一种平翅片换热器（在本文中取两排翅片管），考虑到问题的对称性，取如图中阴影所示的部分为计算区域。在此平翅片的基础上进行不同数目的开缝，得到本文所研究的开缝翅片（图1）。 图1 物理问题示图（单位mm） 本文中计算的为三维稳态层流常物性强制对流模型，其的数学描写如下： (1) (2) 不同变量在方程中(2)的具体形式如表1所示： 表1 方程(2)中的 及 ? u方程 u v方程 v w方程 w 能量方程 T ? 表1 方程(2)中的 及 在计算时，方程对流项采用二阶迎风格式。此功能在FLUENT中通过solve-controls －solution对话框实现。 对开缝翅片进行计算时，进口流速及温度均匀。容积两侧除翅片管外壁以外均为对称边界条件，上下表面、翅片管外壁以及开缝时冲出的小面均为恒定壁温。以上边界条件均在gambit中实现。需要注意的是：在这种情况的数值模拟中，用到了周期性边界条件，即换热器中每一片翅片的几何形状都是完全一样的，这样，在开缝处就可以在垂直方向以周期性边界条件定义，开缝面定义为流体，由于定义为周期性边界条件的目标（本文中为垂直于来流方向的不同小面）要求均为成对且网格生成完全一致，故先将需要定义的目标在gambit中 link faces for mesh，生成网格后后成对定义，再读入FLUENT中。此功能在用FLUENT和gambit软件生成网格时实现。对图1中所示的开缝翅片生成如图2的计算网格。 图2 开缝翅片生成网格示意图 3.有关准则数的定义 本文中Reynolds数、Nusselt数和阻力系数f的定义如下： （3） （3）式中各变量的意义如下： um 计算容积内沿流动方向最窄截面的流速 d3 翅片管外径 va 空气的运动粘性系数 ho 管外换热系数 Q 恒壁温表面换热量 Tw 壁面温度 Tm 空气进出口平均温度 S 换热面积 ka 空气导热系数 △p 进出口压降 ρa 空气密度 L 计算容积沿流速方向的长度 △t取对数平均温差，即 ，其中 ， ? 4.计算结果及分析 本文计算中流体为空气，进口温度为303K,速度由0.5m/s逐渐增加到10m/s。翅片上下表面、翅片管外壁取313K。 1）不同流速下开缝翅片传热特性的比较 保持翅片几何形状不变，改变流体进口流速，其传热性能结果如图4所示。从结果可以看出，流