正六边形蜂窝蓄热体换热动态特性数值模拟研讨.pdfVIP

中国热物理学会 传热传质学 学术会议论文 编号：083016 正六边形蜂窝蓄热体换热动态特性 数值模拟研究 何阳 苑中显 北京工业大学 环境与能源工程学院 教育部传热强化与过程节能重点实验室, 北京 100124 Tel: E-mail: yy7014@ 摘要：采用 Fluent 软件对正六边形通道的蜂窝蓄热体传热进行了数值模拟，考察了蓄热体的换热和阻 力特性，研究了蓄热体达到稳定工作状态的时间与烟气入口速度、蜂窝体厚度的关系，以及蓄热体截 面的温度和速度分布等。通过数值模拟得出：蓄热体达到稳定工作状态的时间与蓄热体的厚度基本呈 线性增加关系。在蓄热期通道内的局部传热系数沿蓄热体长度逐渐减小。气体与固体表面之间的换热 在边长中点附近最强，在六个角点附近换热最弱。 关键词：蓄热体，数值模拟，换热特性，温度分布 1 引言 在工业炉上普遍应用的高温空气燃烧技术主要是通过高效蓄热系统来实现节能和 降低NOx排放量的双重目的[1]。在这种系统中,蓄热体是十分关键也是最能体现技术含 量及工业制造水平的部件，它直接影响蓄热室的大小及热效率和经济效益的高低。评价 它的指标主要有传热性能，压力损失特性等。 对于蜂窝蓄热体的研究，Willmott在大量实验研究的基础上,给出了蓄热体传热计算 的准则关联式[2],但是它们仅适用于严格限定的参数范围内。Klein在薄壁条件下，忽略 蜂窝体在垂直气流方向的温度变化，假设蓄热体导热只在平行气流方向发生的前提下， 提出薄壁模型并求得解析解[3]。Shen和Worek 同时考虑了平行于气体流动方向的蓄热体 导热和垂直于气体流动方向的导热(即二维导热)，用有限差分的方法得到了周期准稳态 情况下的数值解[4]。Ren和Wang在探讨蓄热体长度方向导热的同时，研究了热物性随温 度变化的影响[5]。传统的蜂窝陶瓷为三角形和四边形，实践发现，锐角和直角与钝角相 比更容易引起应力集中。正六边形格孔结构与传统的三角形和四边形相比,应力分布均 匀,结构稳定，不易发生堵塞。当受到应力作用时,应力能从其它两个方向扩散出去，因 此适合于蓄热式换热器的蜂窝陶瓷体[6]。所以对正六边形格孔结构蜂窝体内的换热、压 力损失特性等的研究十分必要。 本文以正六边形蜂窝体为对象，应用FLUENT软件对其进行三维非稳态数值模拟研 究。 2 计算模型及数学描述 2.1 物理模型 正六边形蜂窝体的结构如图 1 所示。由于格孔分布均匀，本文以正六边形半壁厚单 通道蜂窝体为例，来研究蜂窝体内的换热情况。蓄热体材质选用莫来石。由于流过蓄热 体的烟气和空气物性相近，本文采用空气作为传热介质，气体物性按线性变化。单元通 道结构尺寸如图 2 所示，孔边长为 2mm，厚度为 0.8mm,长度为 300mm，由于通道是轴 对称图形，通道内的温度、压力、速度等也是对称的，因此选取通道的一半为流场计算 区域，这样可以大大减少计算网格数量。 图 1 正六边形蜂窝体 图2 单元通道尺寸（单位 mm ） 2.2 换向条件及边界条件 蓄热室从室温条件下开始工作要经历启动过程和稳定工作状态两个阶段，本文模拟 的蓄热室初始温度为300K。在蓄热期，流过格孔的高温烟气将热量传递给蜂窝蓄热体， 蓄热体温度升高；在放热期，低温的空气以相反的方向流过格孔并获得热量，使空气被 预热到较高的温度之后进入炉膛。蓄热体与气体之间的换热反复进行，最后进入准稳定 换热状态。当相邻的两个蓄热周期结束时烟气出口温度相差小于0.3K时，则认为蓄热体 已经达到稳定工作状态。本文换向时间定为30s 。 计算区域外部采用绝热边界条件，对称面为对称边界条件，气体与蜂窝体接触的面 采用耦合求解。烟气入口采用速度入口，入口温度为1200K，入口速度为4m/s ，出口为 压力出口，回流温度为300K,且蓄热期与放热期应保证质量流量相等。并由计算得出蓄 热体内气流的雷诺数小于光滑