作者：张毅         

   本文基于单蓄热介质套管式蓄热器蓄热过程中换热流体沿流向的温度分布特点，提出了多相变材料套管式蓄热器的概念，在尽可能保证高品位热能的条件下，提高相变蓄热过程传热速率。通过Fluent软件研究了蓄热过程流体温度变化、蓄热材料温度分布等特点，为太阳能热利用蓄热系统设计提供基础数据。

1  物理问题及数学模型

1.1物理模型

  多相变材料蓄热器结构如图1(a)所示，蓄热器为圆柱形，蓄热材料被封装在换热流体管道与壳体之间，换热流体在管道内自下向上流动，通过换热管与蓄热材料进行热交换。沿换热体的流动方向，相变材料蓄热腔分为前、后两部分，前端部分的蓄热材料为切片石蜡[12]，其蓄热腔体积占总体积的四分之三，后端部分的蓄热材料为石蜡C16[13]，其蓄热腔体积占总体积的四分之一。两种蓄热材料的物性参数见表1[13] 、表2[12]所示。为便于分析，对物理模型做出如下假设：

  1)各相变材料各向同性；

  2)换热管热阻忽略不计；

  3)除密度外，各相变材料各物性参数为常数，且 固相与液相参数彼此独立；

  4)各相变材料密度变化仅考虑动量方程中与体积有关的项，其余各项中的密度为常数，且满足Boussinesq假设。

  根据上述假设，选取蓄热系统中的某一蓄热单元作为研究对象，如图1 (b)所示。蓄热单元长度为600mm，半径为16mm，换热管半径为4mm。

1.2数学模型

  本文采用Solidification/Melting模型[q建立相变材料的融化模型，其计算区域为图1(b)中的ABCD区域，考虑相变材料溶化后自然对流的影响，基于上述假设，其二维轴对称连续性方程、动量方程及能量方程如下：

1. 连续性方程

  2  网格化分与参数设置

  本文采用Gambit软件建立物理模型并对其进行 网格划分，网格采用四边形单位，其网格数为9600。能量方程、动量方程的对流项离散采用二阶迎风差分格式，压力梯度项采用PRESTO格式，压力和速度的耦合采用SIMPLE算法。相变材料的密度、比热容、导热系数以及粘滞系数采用分段多项式表示（温度的函数）：采用连续性方程和能量方程的残差控制方程的迭代次数，计算中设置连续性方程为10-4，能量方程的残差为10-8，时间步长选取0.2s。

3  结果分析

  图2给出了多相变材料蓄热器与单相变材料蓄热器液相率的比较，由图2可知，与普通套管式相比，采用两种蓄热材料后融化时间缩短约46min，熔化时间缩短了近15%，并且，从图2中可以看出，两种蓄热材料的液相率变化接近一致，可以说二者近似同步熔化，石蜡C16熔化稍稍快些。选取150min与300min两时间点，分析换热流体断面平均温度沿流向的变化，如图3所示，可以看出，在150min时，对于单相变材料蓄热容器，第一段温降最大，以后依次逐渐减小，越接近出口处，温降越来越小，沿换热流体流向，换热流体温度的降低是非线性的，即：越靠近出口端，温度变化越小；而使用多种蓄热材料后，在x=480mm处断面平均温度陡降，直接改变了曲线原来的走势。并且，x=420mm之前的各断面平均温度并没有发生很大变化，说明分段后前四分之三的换热效果并没有受到很大影响，并且，在相变材料改变处，温度梯度大幅增加，在最后四分之一段，温差增大了2倍多二当蓄热过程进行到300min时，使用多种相变材料的蓄热期内蓄热材料熔化基本完全，此时显热蓄热成为主导，可以看到断面平均温度大大升高，说明了显热蓄热效果远不及潜热蓄热效果。

  图4给出了不同截面处断面平均温度随时间的变化曲线，可以看出，熔化刚开始的一小段时间内，温度  变化较快，这是因为相变材料与内导热管直接接触的部分只存在金属导热热阻，热量很快传递到相变材料中，使得这部分相变材料熔化速度很快。然后，传热方式变为导热与对流换热，并且由于相变材料导热率很低，致使金属导热热阻可以忽略不计。对比图4(d)、图(e)、图(f)中可以看出，在选用另一种蓄热材料以后，后四分之一段断面平均温度明显有所降低，并且两条曲线存在交点。将两种情况下换热流体的出口温度进行比较，如图4(f)所示。可以看出，换热流体入口与出口温差明显增大，换热流体出口温度降低5℃左右，温差增大近38%。在流速相同的情况下，温差的增大说明新型套管式蓄热容器的换热效率明显增强，蓄热过程得到了很好的强化。同时，图4(f)中两条曲线存在交点，大约在18500s左右，此时，多相变材料蓄热器内相变材料已近乎全部熔化，出口温升迅速加大。对比图4(d)、图4(e)、图4(f)发现两曲线交点逐渐向右偏移，这种现象产生的主要原因是：单一相变材料出口四分之一段温差很小（如图3所示），而使用两种相变材料后出口四分之一段温差明显增强（如图3所示）．而且熔化加快，断面温升提前，所以出现了交点向右偏移的现象。对比图4(a)、图4(b)、图4(c)发现在熔化开始后大部分时间内，断面平均温度基本相同，这是因为前四分之三段相变材料相同，但不同点是一段 时间后出现分离，分离原因为使用多相变材料后熔化 时间缩短，相变材料熔化到一定程度后率先开始显热蓄热。

  选取30min、60min、100min和200min四个时间段，绘出液相率分布图，如图5所示。图5(a)反映出在蓄热过程中，入口段蓄热材料在0～60min熔化快于其他处蓄热材料，这是由于入口段换热流体温度最高所致，与图5(b)中前四分之三段熔化规律相同，伴随着蓄热材料中自然对流作用的加强，这种优势逐渐减弱，对于图5(b)中的后四分之一段，并没有出现这种规律，是因为后四分之一段较短，温差不大所致。同时图5(a)、图5(b)两图都反映出蓄热过程中左下角与右下角蓄热材料最难熔化，产生这种现象的原因为：换热表面蓄热材料受热熔化后，温度较高的液态蓄热材料在浮升力的作用下向上流动，所以上部蓄热材料最先完全熔化，例如图5(a)中200min液相率图所示；同时温度相对较低的液态蓄热材料向下流动，致使右下角与左下角蓄热材料熔化最慢，如图5(b)中200min液相率图所示。本文中自然对流对相变材料熔化影响与文献[14]中实验规律一致。

4  结论

本文提出了针对换热流体温度降选择不同熔点的蓄热材料这一概念，首先研究了普通光管套管蓄热器换热流体的温度降，针对换热流体在后半段温度降低很慢这一现象，将后四分之一段的蓄热材料换成熔点较低的其他蓄热材料。计算了两种情况下，相变蓄热材料的融化时间、换热流体沿轴向的温度分布并进行了对比分析。模拟结果表明：当换热流体入口与口温差较大时，合理采用不同熔点的蓄热材料能够大大加快蓄热速率。沿换热流体温度降，合理选择熔点不同的蓄热材料，对于相变蓄热有极大的帮助。

5摘  要：本文基于单蓄热介质套管式蓄热器蓄热过程中换热流体沿流向的温度分布特点提出了多相变材料套管式蓄热器，并采用数值模拟的方法研究了该蓄热器的蓄热过程特性，分析了蓄热过程中换热流体沿流向的温度变化特性，蓄热材料的熔化及温度分布特性。与单蓄热介质（切片石腊）相比，多相变材料（切片石腊和石蜡C16）蓄热器换热流体出口温度降低了38%，且达到相同液相率时其蓄热时间缩短了15%。结果表明，在蓄热时，相变材料的熔化温度随换热流体温度降低而相应降低将有利于提高蓄热速率、增大单位时间内的蓄热能力：