0 引言
目前,能源危机和环境问题已经成为了各国亟需解决的问题。近年来全球二氧化碳排放量逐年增加,对生态系统造成了很大的威胁。为了减少二氧化碳的排放,我国将碳达峰、碳中和写入政府工作报告,力争在2030年前二氧化碳排放不再增长。
由于相变材料（phase change material, PCM）在相变过程中会吸收或者放出大量的热量,近年来其在热能储存方面的应用发展迅速[1,2]。相变微胶囊（microencapsulated phase change material, MPCM）使用高分子材料作为外壳,PCM作为芯材。在相变过程中MPCM不会发生形态上的变化。相变微胶囊悬浮液（microencapsulated phase change material slurry, MPCS）通常以水作为基液,将MPCM加入基液并通过物理搅拌、超声震动和添加表面活性剂的方式,使其均匀地分散在基液中形成较为稳定的两相流。MPCS中由于MPCM相变过程中吸收和放出热量以及其与周围流体发生的微对流效应可以强化MPCS的换热效果,其在强化换热和能量储存方面有着广阔的应用前景[3,4,5]。
刘臣臻[10]通过实验研究发现MPCS的对流换热系数随着MPCM质量分数的增大而增大。ZHANG等[11]将MPCS应用于换热器中以研究其对换热性能的影响,实验结果表明,由于PCM在熔化时的吸热效果,当MPCM质量分数为10%时,MPCS的传热效果比水提高了32.8%。ZHANG等[12]分别在层流、过渡流和湍流的条件下进行实验,发现在层流的情况下MPCS的换热系数小于水,在过渡流下MPCS换热接近水,在湍流下MPCS局部努塞尔数（Nux）比水高10%。QIU等[13]通过实验和数值模拟的方法,研究MPCM质量分数、粒径以及过冷度对MPCS在圆管内换热的影响,结果表明MPCM质量分数对MPCS换热效果的影响很大,而MPCM过冷度和粒径对MPCS换热效果影响相对较小。LI等[14]将石墨嵌入MPCM壳材,制备MPCS进行传热实验,在输送同样热量的情况下,使用MPCS比使用水作为传热介质时质量流量和泵功率均显著降低。WANG等[15]发现MPCS在层流转换为湍流时压降发生了明显的转变。在层流下MPCS的局部换热系数随着MPCM浓度的增大而增大,在湍流下MPCS的换热情况优于层流。YUAN等[16]的研究表明10%质量分数的MPCM浆体是最适合的载能流体,但过冷问题仍有待解决。CHEN等[17]的实验表明,在波状微通道中流动的MPCS的换热强度是水的1.7倍。
PU等[18]使用Eulerian-Eulerian模型对MPCS在Y形管内的换热特性进行了模拟分析,MPCM体积分数为12%时,其换热性能相较于纯水提高了38.9%。LIN等[19]建立多尺度模型讨论了不同粒径和浓度下MPCM浆体的传热性能。结果表明,当胶囊尺寸大于100 μm时,较小的胶囊尺寸有利于提高MPCS的换热系数,平均换热系数随浓度的增加先增大后减小,并在MPCM体积分数浓度为20%时达到最大值。YU等[20]使用数值模拟的方法将MPCS应用于提升光伏/热模块整体性能,发现在不同流速下MPCS的传热效率均高于水。LANGURI等[21]采用ANSYS Fluent软件对MPCS在盘管换热器内进行了数值模拟研究,发现接近管壁处的MPCM有着更大的速度。
1 材料制备与表征
1.1 实验材料
实验采用的MPCM以十八烷为芯材,密胺树脂为壳材,由上海儒熵新能源科技有限公司提供。十二烷基硫酸钠（sodium dodecyl sulfate, SDS）（分析纯）和ZnO（分析纯）由上海麦克林生化科技有限公司提供。
1.2 MPCS的制备
ZnO@MPCS制备流程如图1所示。以制备100 g样品为例,将1 g SDS加入93 g的离子水中加热搅拌30 min使SDS完全溶解;称取1 g ZnO加入溶液并搅拌30 min;将5 g的MPCM加入烧杯中以500 r/min的速率搅拌30 min,搅拌后的溶液放入超声波清洗机中超声分散30 min,即可得到100 g的1%ZnO@5%MPCS。
Fig. 1 Preparation flow chart of ZnO@MPCS图1 ZnO@MPCS制备流程图 |
1.3 MPCM和MPCS的表征
Fig. 2 SEM images of MPCM图2 MPCM的SEM照片 |
Fig. 3 Particle size distribution of MPCM图3 MPCM的粒径分布 |
MPCS的密度通过质量流量计（德国,KEM,TCM-3100）测得,使用导热仪（南京大展机电技术有限公司,DZDR-S）测定了MPCS的导热系数。使用差示扫描量热仪（differential scanning calorimetry, DSC）（德国,Netzsch,200F3 Maia）测量MCPM和MPCS的热物性,所有测量均在氮气下以2℃/min的加热和冷却速率进行。MPCS物性参数如表1所示。MPCM和MPCS的DSC测试结果如图4所示,MPCM表现出了较高的焓值,其熔化热为189.2 J/g,结晶热为187.3 J/g。从图4 MPCS以及ZnO@MPCS的DSC曲线可以看出,ZnO的加入并没有降低MPCS的相变潜热。
Table 1 Physical parameters of MPCS and water表1 MPCS和水物性参数 |
| 传热介质 | 密度 / (kg/m3) | 导热系数 / [W/(m∙℃)] | 熔化潜热 / (kJ/kg) | 凝固潜热 / (kJ/kg) |
|---|---|---|---|---|
| 水 | 998.2 | 0.598 | / | / |
| 5%MPCS | 991.3 | 0.513 | 8.146 | 7.444 |
| 1%ZnO@5%MPCS | 996.8 | 0.605 | 8.579 | 8.774 |
Fig. 4 DSC curves of MPCM and MPCS图4 MPCM和MPCS的DSC曲线 |
使用流变仪（英国,Malverninstruments,Kinexus Ultra）测量MPCS在不同温度下的动力黏度,结果如图5所示。可以看出MPCS和ZnO@MPCS的黏度在剪切速率大于30 S-1后趋于稳定。ZnO颗粒的加入增大了MPCS的黏度,随着温度的升高,MPCS以及ZnO@MPCS的黏度逐渐减小。
Fig. 5 Dynamic viscosity for MPCS and ZnO@MPCS at 20℃ and 35℃ versus shear rate图5 MPCS和ZnO@MPCS在20℃和35℃时黏度与剪切速率的关系 |
2 对流换热实验
2.1 实验装置和步骤
流动换热试验台布置如图6所示,实验系统包括储液水箱、质量流量计、试验段、压差变送器（杭州联测自动化技术有限公司,SIN-2051G）、调压器（上海人民电器厂,tdgc2）、加热圈、离心泵（威乐水泵系统有限公司,PUN-403EH）、板式换热器（宁波北仑博克冷却设备有限公司,K030-10F-1）、低温冷却剂循环泵（巩义市科瑞仪器有限公司,DLSB-5/10）和数据采集系统（英国,ZEEFAX,IMP-400）;试验段为光滑铜管,长度1 500 mm,外径8 mm,内径6 mm;实验装置的其他管道均为不锈钢管。
Fig. 6 Schematic diagram and physical drawing of flow heat exchange test bench图6 流动换热试验台原理图及实物图 |
耐高温绝缘加热圈均匀缠绕在铜管上,产生均匀的热流,在加热圈两端连接调压器通过调节电压来改变加热圈的加热功率。加热圈的加热功率由功率计测量。试验段的管道覆盖有40 mm厚的硅酸铝纤维保温层,以减少实验系统对环境的热量损失。在距离铜管入口分别为0.05 m、0.25 m、0.45 m、0.65 m、0.85 m、1.05 m、1.25 m和1.45 m的管壁上安装了8个K型热电偶用来测量不同位置铜管外壁温度。在铜管进口、铜管出口和板式换热器出口安装3个K型铠装热电偶,用于测量进出铜管的流体温度以及通过板式换热器冷却后的流体温度。热电偶测量数据由IMP-400数据采集仪采集,从计算机端输出。在铜管两端连接差压变送器,测量流体通过实验铜管后的压降。质量流量计安装在铜管进口,用于测量试验段的质量流量和密度。
在实验过程中,MPCS在循环水泵的作用下通过流动发展段进入质量流量计,然后流入铜管进行加热,被加热后的流体与板式换热器冷却器中的冷却水进行热交换。使得流出板式换热器的微胶囊相变材料浆料中的相变材料完全固化,最后实验流体流回储液罐,准备进行新的循环。在实验的过程中,将搅拌器调至合适的速率,以保持悬浮液的稳定性,同时减少实验中泡沫的产生。
2.2 数据处理与误差分析
圆管进口处雷诺数Re为：
$Re=\frac{{{u}_{\text{m}}}D}{v}$ (1)
实验段圆管内壁面温度通过下式计算：
${{T}_{\text{i}}}\text{=}{{T}_{\text{wo}}}-{{C}_{\text{1}}}\frac{{{q}_{\text{w}}}}{{{k}_{\text{wo}}}}\left( 1+\frac{{{C}_{\text{1}}}}{2}\frac{a{{q}_{\text{w}}}}{k_{\text{wo}}^{\text{2}}} \right)$ (2)
上式中,${{C}_{1}}\text{=}\frac{{{r}_{\operatorname{w}}}}{r_{\text{wo}}^{2}-r_{\text{w}}^{2}}\left( r_{\text{wo}}^{\text{2}}\ln \frac{{{r}_{\text{wo}}}}{{{r}_{\text{w}}}}-\frac{r_{\text{wo}}^{2}-r_{\text{w}}^{2}}{2} \right)$。其中：kwo为铜管的导热系数,其值为398 W/(m·K);a为常数,其值为 0.014 4。
实验中流体与圆管管壁间局部换热系数hx为：
${{h}_{x}}\text{=}\frac{{{q}_{\text{w}}}}{{{T}_{\text{i}}}-{{T}_{\text{b}}}}=\frac{UI/\left( \text{ }\!\!\pi\!\!\text{ }DL \right)}{{{T}_{\text{i}}}-{{T}_{\text{b,}x}}}$ (3)
层流理论Nux计算式：
$N{{u}_{x}}\text{=5}\text{.364}{{\left[ 1\text{+}{{\left( \frac{220{{x}^{\text{*}}}}{\text{ }\!\!\pi\!\!\text{ }} \right)}^{-10/9}} \right]}^{3/10}}-1.0$ (4)
式中：${{x}^{*}}=\left( x/D \right)\left( R{{e}_{\text{b}}}\cdot P{{r}_{\text{b}}} \right)$
Nux实验值计算式：
$N{{u}_{x}}=\frac{{{h}_{x}}D}{{{k}_{\text{b}}}}$ (5)
在湍流时,理论Nux通过下式计算：
$N{{u}_{x}}=0.027R{{e}^{0.8}}P{{r}^{\frac{1}{3}}}{{\left( \frac{\mu }{{{\mu }_{\text{w}}}} \right)}^{0.14}}$ (6)
MPCS的比热容计算式[25]：
${{C}_{p\text{,s}}}={{M}_{\text{wa}}}\cdot {{C}_{p\text{,wa}}}+{{M}_{\text{mp}}}\cdot {{C}_{p\text{,mp}}}+{{M}_{\text{mo}}}\cdot {{C}_{p\text{,mo}}}$ (7)
流体局部平均温度Tx计算式为：
${{T}_{x}}={{T}_{\text{b,i}}}+\frac{{{Q}_{\text{w}}}P}{\dot{m}{{C}_{p,\text{s}}}}x$ (8)
流体的吸热量Qw通过下式计算：
${{Q}_{\text{w}}}\text{=}\frac{\dot{m}{{C}_{p,\text{s}}}\left( {{T}_{\text{b,o}}}-{{T}_{\text{b,i}}} \right)}{A}$ (9)
实验中所有的热电偶都使用恒温水浴进行了校准,误差控制在 ±0.2℃,质量流量计误差在 ±0.1%,功率计最大误差为 ±2%。平均Nux不确定度为6.18%。
（1）在实验过程中,实验段的热损失主要集中在周围环境和轴向管壁上。轴向管壁的热损失导致管道进出口温度较低,使得进口处hx大于理论值。
（2）虽然实验中加热圈均匀地缠绕在管壁上,但由于加热圈的加热特性,热流并非绝对均匀,这会对管壁温度产生一定的影响。
Fig. 7 hx of water in laminar flow |
Fig. 8 hx of water in turbulent flow |
3 结果及分析
3.1 水与MPCS压差对比分析
通过压差变送器测量不同流速下水、MPCS和ZnO@MPCS的压降,结果如图9示,ZnO的加入增大了测试管道中浆体的压降。但ZnO@MPCS平均压降仅比MPCS高5.44%,ZnO的加入对水泵输送MPCS功率的影响很小。
Fig. 9 Comparison of pressure drop between MPCS and water图9 MPCS与水压降对比 |
3.2 层流下水与MPCS换热特性
为了研究层流条件下金属氧化物对MPCS传热特性的影响,对水、5%MPCS和5%MPCS@1%ZnO在入口质量流量为0.45 kg/min、加热功率为600 W的条件下进行强制对流传热实验。图10显示了试验段不同位置的内壁温度和局部流体温度沿管长的分布。由于PCM发生相变,其相变过程需要吸收大量的热量,故MPCS和ZnO@MPCS的壁温在试验段距入口处的0.25 ~ 1.05 m之间缓慢上升。在距离入口处1.05 m之后,大部分的相变材料完成相变,铜管壁温升高幅度大于0.25 ~ 1.05 m处。
Fig. 10 Inner wall temperature and local average temperature of water, MPCS and ZnO@MPCS图10 水、MPCS和ZnO@MPCS的内壁面温度和流体局部平均温度 |
图11a显示了hx与试验段轴向管长之间的关系。由于MPCM颗粒和ZnO颗粒的加入增大了流体的黏度,故MPCS和ZnO@MPCS在入口处hx小于水。在距入口处0.25 m后,相变材料开始逐渐熔化。相变材料在熔化的过程中需要吸收大量的热量,MPCS和ZnO@MPCS的黏度随着流体温度的升高而减小,流体的换热能力逐渐增强,因此MPCS和ZnO@MPCS的hx逐渐增大。在距离测试段入口1.05 ~ 1.45 m处,MPCS和ZnO@MPCS的hx先减小后增大,这是由于1.05 m后大部分PCM完成了相变,流体需要吸收的热量相较于相变过程所需的热量变少,导致hx减小;在1.3 m后,随着流动的发展,流体的温度升高、黏度减小且局部雷诺数增大,这些原因使MPCS在测试段末尾的对流换热系数增大。
Fig. 11 hx (a) and Nux (b) of water, MPCS and ZnO@MPCS图11 水、MPCS和ZnO@MPCS的hx（a）和Nux（b） |
金属氧化物的加入提高了MPCS的导热系数,使浆料与壁面间换热增强。总体来说,层流条件下5%MPCS的平均hx小于水,但加入1%ZnO的5%MPCS平均的hx比水高6.5%。ZnO@MPCS的局部换热系数在整个换热过程中均大于MPCS,这是由于ZnO的加入提升了MPCS的导热系数,而且金属颗粒的运动促进了流体间以及流体和管壁之间的热量传递。
图11b显示了Nux和轴向管长度之间的关系。可以看出,在入口处,MPCS和ZnO@MPCS的Nux小于纯水,但当PCM逐渐开始熔化后,MPCS和ZnO@MPCS的Nux大于水,并且随着流动的发展,MPCS和ZnO@MPCS的Nux均明显增大,表明MPCM和ZnO的加入增强了流体与管壁之间的对流换热强度。
3.3 湍流下水与MPCS换热特性
为进一步研究金属氧化物对MPCS传热特性的影响,在湍流下对水、5%MPCS和1%ZnO@5%MPCS进行强制对流换热实验。图12显示了水、MPCS和ZnO@MPCS的局部平均温度和铜管内壁温度。MPCM的加入使流体温度与壁面间的换热增强,流体局部平均温度接近圆管内壁面温度。加入ZnO的MPCS流体与圆管内壁温差小于未加入ZnO的MPCS。湍流下流体间的扰动更加剧烈,流体与圆管壁面间的换热增强,使流体与壁面间的温差小于层流。
Fig. 12 Inner wall temperature and local average temperature of water, MPCS and ZnO@MPCS图12 水、MPCS和ZnO@MPCS的内壁面温度和流体局部平均温度 |
图13a为湍流下hx沿轴向管长的分布情况。随着入口壁温的升高,壁温达到相变材料的熔点,靠近管壁的相变材料开始熔化,PCM在熔化过程中需要吸收更多的热量,流体在入口处的hx呈现增大的趋势。在距离圆管入口0.5 m处,有更多的PCM开始熔化,此时圆管与流体之间的hx接近最大值。随着相变材料熔化量的减少,hx逐渐减小。在测试段的末尾,流体的hx随雷诺数的增大和浆液黏度的降低而增大。湍流下MPCS的扰动更加剧烈,MPCS的平均hx较水高8.9%,加入1% ZnO的MPCS的hx较水提高了15.7%。MPCS以及ZnO@MPCS的hx均大于层流,这是由于流速的增加增大了浆体的雷诺数,使得颗粒、流体和壁面之间的换热增强。
Fig. 13 hx (a) and Nux (b) of water, MPCS and ZnO@MPCS图13 水、MPCS和ZnO@MPCS的hx（a）和Nux（b） |
MPCS和ZnO@MPCS在湍流下沿管道长度的局部Nux分布如图13b所示。5%MPCS以及1%ZnO@5%MPCS在湍流条件下均实现了传热增强。在入口处由于ZnO@MPCS的黏度较大,致使其Nux略低于纯水,随着流体温度的升高,ZnO@MPCS的黏度逐渐减小,在入口段之后,MPCS和ZnO@MPCS的Nux均高于水。
4 结论
使用ZnO作为高导热颗粒,制备具有较高导热率的MPCS并对MPCS的物理性质进行测试。设计搭建实验台,分别在层流和湍流下进行了强制对流换热实验,通过对比水,MPCS和ZnO@MPCS的局部换热系数hx和局部努塞尔数Nux,得出以下结论：
（1）制备出的1%ZnO@5%MPCS导热系数较5%MPCS提高了17.9%,1%ZnO@5%MPCS的黏度略大于5%MPCS。ZnO@MPCS具有良好的导热和储热性能。
（2）层流条件下MPCS的平均hx低于水,1%ZnO@5%MPCS平均hx比水高6.5%。湍流时, 1%ZnO@5%MPCS在相同质量流量和功率下的hx比水提高了9.6%。
（3）ZnO对MPCS局部换热的增强是由于ZnO的加入提升了MPCS的导热系数,并且加强了流体间的微对流现象。随着流动的发展,流体温度的升高,流体的黏度降低以及局部雷诺数的增大也是MPCS和ZnO@MPCS的hx增强的原因。ZnO的加入增大了MPCS的黏度,但从换热效率的角度来看,ZnO@MPCS的换热效率更高。在输送或交换同等热量的情况下,所需的流体更少,可以减少水泵的功耗。
符号表： |
| A | 铜管表面积,m2 |
|---|---|
| Cp | 比热容,kJ/(kg∙℃) |
| D | 圆管直径,m |
| hx | 局部换热系数,W/(m2∙℃) |
| I | 电流,A |
| k | 导热系数,W/(m∙℃) |
| L | 长度,m |
| M | 质量分数,% |
| m | 质量流量,kg/s |
| Nux | 局部努塞尔数 |
| P | 铜管周长,m |
| Pr | 普朗特数 |
| Q | 加热功率,W |
| Qw | 流体吸热量,J |
| qw | 热流密度,W/m2 |
| Re | 雷诺数 |
| r | 圆管半径,m |
| T | 温度,℃ |
| U | 电压,V |
| u | 流速,m/s |
| v | 流体的运动黏度,m2/s |
| x | 距离圆管入口处的距离 |
下角标： |
| b | 局部 |
|---|---|
| i | 入口 |
| m | 平均 |
| mo | 金属氧化物 |
| mp | 相变微胶囊 |
| s | 浆料 |
| w | 内壁面 |
| wa | 水 |
| wo | 外壁面 |