能源研究与信息  2023, Vol. 39 Issue (4): 242-249   PDF    
柱状颗粒填充床的储能性能分析
兰健, 郭雪岩     
上海理工大学 能源与动力工程学院,上海 200093
摘要:填充床储能是一种很有发展前景的热能储存技术,它具有可降低存储成本和提高太阳能热系统开发效率等优点。研究人员多采用球形的储能单元,而圆柱体在储能填充床换热中有其独特的优势,因此基于圆柱形和拉西环形两种柱状颗粒,建立了一种潜热储能填充床的三维模型,采用数值模拟的方法分别研究两种柱状颗粒组成的填充床的储能性能,分析了储能填充床的直径比对其性能的影响。研究表明,填充床直径比越大,其储能性能越好。同时研究了圆柱形储能单元高度和拉西环形储能单元孔径对储能性能的影响。结果表明,在研究范围内,由高度为3 mm的圆柱形储能单元和孔半径1.50 mm的储能单元分别组成的填充床储能速率最高。
关键词储能     热能     相变材料     储能单元    
Analysis of energy storage performance of a packed bed with cylindrical particles
LAN Jian, GUO Xueyan     
School of Energy and Power Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China
Abstract: Packed bed is a promising heat storage way with the advantages of reduced storage cost and improved development efficiency of solar thermal system. The spheric energy storage unit was frequently adopted by researchers, while the cylindrical one also had some unique advantages in energy storage. Based on two column particles such as cylinder and Raschig ring, a three-dimensional packed bed model for latent heat storage was established. Their energy storage performance in the packed bed was studied by numerical simulation. Effect of the diameter ratio of bed to particle on its performance was studied. Results show that the larger the diameter ratio of bed to particle, the better the energy storage performance. Effects of cylinder height and Raschig ring pore size in the energy storage unit on energy storage performance were studied. Results show that within the investigated range, the energy storage unit in the packed bed with the cylinder height of 3 mm or pore size of 1.50 mm had the highest energy storage rate.
Key words: energy storage     thermal energy     phase change materials     energy storage unit    

目前,全球正在大规模发展可再生能源。未来20年,全球范围内的可再生能源装机总量会每年以2.3%的速度增长。而我国更是在2006年就已经通过并且实施了《中华人民共和国可再生能源法》1,以推动可再生能源的发展。可再生能源也有无法避免的劣势——间歇性和波动性,这使得其面临较大的挑战。可再生能源装机容量在发达国家的高速增长对电价市场产生了较大的影响。在一些国家,由于太阳能和风力资源比较充足,进而出现了“负电价”2-3。与此同时,由于可再生能源装机容量的增长,基荷电厂需要频繁地进行功率调节,这无疑会损耗其运行寿命4

为保证可再生能源高速、健康发展,需要一种令人满意的技术来储存能量,并将其转化为所需形式的能量。填充床储能是一种很有发展前景的热能储存技术,它具有可降低存储成本和提高太阳能热系统开发效率等优点。储能技术按照储能方式主要分为三种:显热储能、相变储能和热化学储能。顾名思义,显热储能即只通过材料的显热来储存能量,因此材料为非相变材料;而潜热储能则同时利用储能材料的显热和潜热,故其材料为相变材料;热化学储能通过储能材料的反应热来实现能量储存。三种储能方式中,热化学储能的储能密度远大于其他两种,但其系统组成和维护较复杂,使得其在实际应用中有许多的限制,无法发挥其储能密度大的优势5-7。而其他两种储能方式中,显热储能的储能密度小于潜热,而相比于热化学储能,潜热储能系统的组成和维护较简单,其在太阳能发电、工业余热利用等相关领域得到广泛应用8-11

储能填充床由储能单元与换热流体组成。在已有换热流体中12-17,熔盐具有使用温度范围广、成本相对较低等优势。在储能单元方面,对于相变储能填充床,有学者18-19使用价格相对低廉且易于获取的石蜡作为相变材料,对储能填充床进行了相关研究,得出了较好的结果。Peng等20利用NaNO2制成的相变胶囊作为储能单元进行了相关研究,其储能性能明显优于石蜡,但由于未与非相变材料形成直接进行对比,无法体现出可量化的优势。同时,有研究者21利用陶瓷和石蜡制成了一种多颗粒相变复合材料,研究了不同混合比下储能填充床的性能,在两者之间得到一个较优混合比。Li等22以混合相变材料作为储能单元,分别进行了数值模拟和实验研究。除此之外,Li等23还制作出了两种不同粒径的相变胶囊,将大、小粒径的相变胶囊分层堆叠进行相关研究。大、小粒径相变胶囊的分层堆叠,对换热流体的流动起到了优化作用,同时也可以将大、小粒径的相变胶囊与换热流体的换热进行对比,得出粒径对其换热的影响。

储能填充床实现储能功能主要依靠储能单元与换热流体间的换热以及冷、热流体之间的换热,两部分共同作用使得储能量达到最大。熔盐在中高温填充床储能领域被认为是现阶段最佳换热流体之一。因此,提高储能填充床的性能,需要从储能单元与换热流体间的换热入手。首先,需要考虑的是储能单元的结构。在实际研究过程中,考虑到填充床内储能单元堆积过程中的简易性,研究人员多采用球形储能单元,而圆柱形储能单元则很少被使用,实际上,圆柱形储能单元由于结构上的特殊性,在堆积时有明显的方向性,并且结构上几何元素多,堆积的效果与球形储能单元明显不同。本文基于柱状颗粒建立了一种潜热储能填充床的三维模型,对圆柱与拉西环两种结构储能单元进行相关研究。

1 物理及数学模型 1.1 物理模型

本文填充床使用管径D=30 mm的圆管,圆管长度为150 mm,其中在进、出口处分别空出30 mm和20 mm,以消除进、出口效应,剩余100 mm为填充段长度,使用Partopour24的方法随机填充储能单元。填充床几何模型如图1所示。在填充床的物理模型生成过程中,为了避免储能单元之间以及储能单元与壁面之间的接触,将储能单元体积减小至其原始体积的99%。本文使用的网格类型为多面体网格,该网格具有较好的适应性。除此之外,在相同的结构下,多面体网格的数量相较于传统四面体网格的更少,在不影响计算结果准确性的基础上,对于计算资源要求更低,网格局部如图2所示。本研究主要采用结构为圆柱形、拉西环形。本研究中储能单元的相变模型为凝固与熔化模型。该模型中假设相变材料由其他固体非相变材料包裹,其厚度忽略不计。本研究中固体非相变材料为陶瓷,陶瓷在换热过程中不会发生相变,故相变材料发生熔化后不会与换热流体混合。下文中除特殊说明外,圆柱的底面直径与高均为6 mm,拉西环在此基础上孔半径为1.50 mm。

图 1 填充床几何模型 Fig.1 Geometric model of packed bed

图 2 网格局部 Fig.2 Local grid enlargement
1.2 数学模型

在本文采用的数值研究中,基于求解器ANSYS Fluent软件,采用kω SST模型作为湍流模型,使用SIMPLE算法二阶离散格式对压力速度进行耦合。采用速度入口作为入口边界条件,出口为压力出口,忽略填充床外壁面的散热,故设置其为绝热壁面。本研究使用的流体熔盐25密度为1 865 kg·m−3,比热容为1 510 J·kg−1·K−1,导热系数为0.571 W·m−1·K−1,黏度为3.22×10−3 Pa·s;相变材料为三元熔盐22,其密度为2 310 kg·m−3,比热容为1 540 J·kg−1·K−1,导热系数为1.69 W·m−1·K−1,潜热为273 kJ·kg−1。其中,换热流体为熔盐,相变材料为三元熔盐,其由Li2CO3、K2CO3、Na2CO3分别以32%、35%、33%的质量分数制成,其相变温度为668.1 K22

凝固与熔化模型中,其焓值H为显焓$h$与潜焓$\Delta H$之和,即

$\qquad H = h + \Delta H $ (1)
$\qquad h = {h_{\rm{r}}} + \displaystyle \smallint \nolimits_{{T_{\rm{r}}}}^T {c_p}{\rm{d}}T $ (2)
$\qquad \Delta H = \eta R $ (3)

式中:${h_{\rm{r}}}$为参考焓值;${T_{\rm{r}}}$为参考温度;${c_p}$为比热容;$R$为潜热;$\eta $为液相分数;T为材料温度 。

$\eta $表达式为

$\qquad \left\{ \begin{array}{l} \eta = 0\;\left( {T < {T_{\rm{g}}}} \right)\\ \eta = 1\;\left( {T > {T_{\rm{h}}}} \right)\\ \eta = \dfrac{{T - {T_{\rm{g}}}}}{{{T_{\rm{h}}} - {T_{\rm{g}}}}}\;\left( {{T_{\rm{g}}} < T < {T_{\rm{h}}}} \right) \end{array} \right. $ (4)

式中:${T_{\rm{g}}}$为凝固温度,${T_{\rm{h}}}$为熔化温度。

1.3 数值模型验证

由于填充床本身无法验证,本文利用Whitaker26提出的球体对流传热关联式来验证模拟结果。由于该关联式只针对球体周围的对流换热,而对圆柱与拉西环结构无法直接进行对比,所以采用单个球体的结构来验证数值计算方法的合理性。式(5)为关联式,式(6)为平均努塞尔数计算式,分别为

$ N{u_{\rm{m}}} = 2 + \left( {0.4R{e^{1/2}} + 0.06R{e^{2/3}}} \right)P{r^{0.4}}{\left( {\frac{{{u_\infty }}}{{{u _0}}}} \right)^{1/4}} $ (5)
$\qquad Nu = \frac{{kl}}{\lambda } $ (6)

式中:$k$为对流换热系数;$l$为等效球直径;$\lambda $为三元熔盐导热系数;Re为雷诺数,其值基于球直径5 mm来计算;Pr为普朗特数;${u_\infty }$为主流速度;${u_0}$为入口速度;$N{u_{\rm{m}}}$为关联式计算出的努塞尔数;$Nu$为计算式得出的努塞尔数。

$\qquad k=\dfrac{Q}{A\Delta t} $ (7)

式中:$Q$为换热量;$A$为表面积;$ \text{Δ}{t} $为温差。

图3为利用本文的数值方法基于单个球体得到的计算值与关联式计算值的差别。图中可以看出,在不同雷诺数下,两个方法的努塞尔数差别不大,基本吻合。

图 3 单颗粒努塞尔数与雷诺数的关系验证 Fig.3 Relationship between Nusselt number and Reynolds number verified with single particle

本研究中,储能填充床运行温度为573 ~ 773 K,储能截止温度为772 K,放能截止温度为574 K,均以出口位置储能单元温度为准。储能填充床所吸收/释放能量为填充床进、出热通量绝对值的差值。

2 结果与讨论 2.1 储能单元在储放能过程温度变化规律

在储放能过程中,入口速度均为0.001 m·s−1,填充床运行温度为573 ~ 773 K。图4为不同轴向位置拉西环温度随时间t的变化。在储能过程中,拉西环储能单元由左往右逐渐熔化。储能结束时间为334.38 s。每个储能单元的储能过程分为三个阶段——固体显热吸热、相变潜热储热、液体显热吸热。在固体显热吸热阶段,储能单元吸收热量后温度升高;在相变潜热储热阶段,相变储能单元因相变潜热吸收大量热量,温度基本不变;在液体显热吸热阶段,相变材料熔化过程结束并且持续吸收热量,温度继续升高。储能过程结束后,储能单元温度保持在773 K。在放能过程中,储能单元由右往左逐渐凝固。该过程与储能过程类似。从图中可以看出,z=40 mm和z=120 mm两个轴向位置的储能单元在储能、放能过程中相变潜热阶段的温度变化相对于其他位置均不明显,原因是这两个位置在相应的过程中均处在换热流体的入口附近,换热温差相对较大,从而换热速率较快。在Li等23的数值计算结果中,其潜热吸热(放热)阶段时长随着与进、出口距离的减小而减少,但并未减少到无法明显看出的程度,因为其距离进、出口最近的位置x/H0为0.25(x为填充段到未填充段的距离,H0为填充段总长度),本文中x/H0为0.1,距离进、出口更近;其次,其填充床以及储能单元尺寸远大于本文中的尺寸,潜热阶段所需能量更多,故该阶段时长更长,因此其潜热阶段更为明显。

图 4 不同轴向位置拉西环温度随时间的变化 Fig.4 Changes of Raschig ring temperature with time at different axial positions
2.2 不同直径比下储能性能

颗粒与填充床的直径比决定了颗粒分布的复杂程度,进而影响填充床内流场的混乱程度。分别研究直径比D/d=3.0、5.0、7.5(d为颗粒底面直径)的拉西环填充床(拉西环高度为6 mm,d分别为10、6、4 mm,孔半径为1.50 mm)。三种填充床堆积模型如图5所示。

图 5 不同直径比下填充床堆积模型 Fig.5 Packing models with different diameter ratios of bed to particle

图6为不同直径比填充床z=80 mm处储能单元温度随时间的变化和填充床的储能时间和储能量。三种填充床空隙率分别为0.59、0.60、0.77,储能速率分别为167.18、240.76、267.66 W。三种填充床的储能量、储能时间与空隙率均成反比。综合考虑储能时间与储能量,填充床的储能速率随着空隙率的增大而增大。直径比为7.5的填充床储能单元数量多,且孔半径与拉西环底面直径比大,故其换热面积大。在前75 s左右,由于换热速率高,前段储能单元吸热量多,在填充床中段的换热温差较小,则前期该处储能单元升温慢。75 s后,三种填充床前段储能单元均换热充分,中段温差基本一致,对于直径比为7.5的填充床,其储能单元换热面积大,温度快速上升。

图 6 不同直径比下填充床的储能过程及储能性能 Fig.6 Energy storage process and performance of packed bed with different diameter ratios of bed to particle
2.3 不同储能单元结构下储能性能 2.3.1 不同圆柱高度下储能性能

圆柱储能单元的高度会直接影响100 mm的填充段中储能单元的个数以及分布,进而影响填充床的流场。因此,对由高度h分别为3、6、9 mm的圆柱储能单元(底面直径均为6 mm)组成的填充床进行研究。不同圆柱高度下填充床堆积模型如图7所示。

图 7 不同圆柱高度下填充床堆积模型 Fig.7 Packing models under different cylinder heights

图8为不同圆柱高度填充床z=80 mm处储能单元温度随时间的变化与填充床的储能时间和储能量。随着圆柱高度的增大,储能量基本呈降低的趋势,三种结构填充床填充段的空隙率分别为0.479、0.484、0.495,储能速率分别为228.73、204.77、204.75 W,这意味着高度为3 mm的填充储能单元所占体积最大,而单位体积下储能单元所能储存的能量比换热流体的大,因此高度为3 mm的填充床储能量最多。如果将储能量与储能时间的比值作为储能速率,则储能速率随着孔隙率的增大先减小后基本不变。由圆柱高度为3 mm的颗粒组成的填充床储能单元个数多,换热面积大。前期储能单元换热速率较高,前段储能单元吸热量大,则中段储能单元换热温差小,因此前期该处储能单元升温慢。之后,三种填充床前段储能单元均换热充分,中段温差基本一致,对于圆柱颗粒高度为3 mm的填充床,其内部换热面积大,温度上升较快。

图 8 不同圆柱高度下填充床的储能过程及储能性能 Fig.8 Energy storage process and performance of packed bed with different cylinder heights sizes
2.3.2 不同拉西环孔半径下的储能性能

通过改变拉西环孔半径可在保持储能单元数基本不变的情况下改变填充床的空隙率,进而影响填充床内部的换热。因此,对孔半径r分别为0.75、1.50、2.25 mm的拉西环(拉西环底面直径和高度均为6 mm)组成的填充床进行相关研究。三种孔径下填充床堆积模型如图9所示。

图 9 不同孔半径下填充床堆积模型 Fig.9 Packing models under different pore sizes

图10为不同孔半径填充床z=80 mm处储能单元温度随时间的变化和填充床的储能时间与储能量,其中孔半径为0表示圆柱储能单元。四种填充床空隙率分别为0.48、0.52、0.60、0.77,储能速率分别为204.77、211.12、240.76 、219.50 W,对比孔半径为0、0.75 mm的储能单元组成的填充床,发现孔半径为0.75 mm的换热面积大,换热速率高,则在储能截止时间时,末段颗粒与出口之间的换热流体温度高于孔半径为0的填充床,因此0.75 mm的填充床储能量较大。同样的,孔半径为1.50、0.75 mm的填充床也有这种情况,但孔半径为2.25 mm的填充床中储能单元的体积占比过小,这部分能量小于储能单元之间能量差,所以孔半径为1.50 mm的填充床储能量较大。就储能速率而言,由孔半径为1.50 mm的拉西环组成的填充床储能速率最高。

图 10 不同孔径下填充床的储能过程及储能性能 Fig.10 Energy storage process and performance of packed bed with different pore sizes
3 结 论

采用数值模拟方法分别对圆柱和拉西环结构储能单元组成的填充床进行相关研究,采用储能时间、储能量来表示填充床性能,分析了填充床直径比与储能单元结构参数对填充床性能的影响,得出如下结论:

(1) 随着填充床直径比的增大,储能速率逐渐加快。具体而言,直径比为7.5的填充床的储能单元个数约为直径比为3的填充床的个数的7.4倍。相较于直径比为3的填充床,直径比为7.5的填充床储能时间减少了45.0%,但同时储能量仅减少了11.9%,综合两点考虑,其储热速率提高了60.1%。

(2) 圆柱储能单元高度的增加对填充床储能量的影响较小,储能时间先增加后基本不发生变化,综合而言,储能速率先减小后基本不变,所以在圆柱高度研究范围内,高度为3 mm的圆柱组成的填充床储能速率最高。

(3) 随着储能单元孔径的增大,储能量呈先增后减的趋势,但储能时间相差不大,故在孔径研究范围内,储能速率同样出现先增后减的规律,其中孔径为1.50 mm的拉西环组成的填充床储能速率最高。

参考文献
[1]
薛惠锋, 王海宁. 《中华人民共和国可再生能源法》的实施回顾及展望[J]. 中外能源, 2010, 15(3): 33-36.
[2]
DE LAGARDE C M, LANTZ F. How renewable production depresses electricity prices: evidence from the German market[J]. Energy Policy, 2018, 117: 263-277. DOI:10.1016/j.enpol.2018.02.048
[3]
VALITOV N. Risk premia in the German day-ahead electricity market revisited: the impact of negative prices[J]. Energy Economics, 2019, 82: 70-77. DOI:10.1016/j.eneco.2018.01.020
[4]
LOCATELLI G, BOARIN S, FIORDALISO A, et al. Load following of Small Modular Reactors (SMR) by cogeneration of hydrogen: a techno-economic analysis[J]. Energy, 2018, 148: 494-505. DOI:10.1016/j.energy.2018.01.041
[5]
PRIETO C, COOPER P, FERNÁNDEZ A I, et al. Review of technology: thermochemical energy storage for concentrated solar power plants[J]. Renewable and Sustainable Energy Reviews, 2016, 60: 909-929. DOI:10.1016/j.rser.2015.12.364
[6]
CHEN X Y, ZHANG Z, QI C G, et al. State of the art on the high-temperature thermochemical energy storage systems[J]. Energy Conversion and Management, 2018, 177: 792-815. DOI:10.1016/j.enconman.2018.10.011
[7]
吴娟, 龙新峰. 太阳能热化学储能研究进展[J]. 化工进展, 2014, 33(12): 3238-3245.
[8]
SHAO J J, DARKWA J, KOKOGIANNAKIS G. Review of phase change emulsions (PCMEs) and their applications in HVAC systems[J]. Energy and Buildings, 2015, 94: 200-217. DOI:10.1016/j.enbuild.2015.03.003
[9]
SHUKLA A, BUDDHI D, SAWHNEY R L. Solar water heaters with phase change material thermal energy storage medium: a review[J]. Renewable and Sustainable Energy Reviews, 2009, 13(8): 2119-2125. DOI:10.1016/j.rser.2009.01.024
[10]
SHARMA R K, GANESAN P, TYAGI V V, et al. Developments in organic solid-liquid phase change materials and their applications in thermal energy storage[J]. Energy Conversion and Management, 2015, 95: 193-228. DOI:10.1016/j.enconman.2015.01.084
[11]
TAO Y B, HE Y L. A review of phase change material and performance enhancement method for latent heat storage system[J]. Renewable and Sustainable Energy Reviews, 2018, 93: 245-259. DOI:10.1016/j.rser.2018.05.028
[12]
PACHECO J E, SHOWALTER S K, KOLB W J. Development of a molten-salt thermocline thermal storage system for parabolic trough plants[J]. Journal of Solar Energy Engineering, 2002, 124(2): 153-159. DOI:10.1115/1.1464123
[13]
YANG X P, YANG X X, QIN F G F, et al. Experimental investigation of a molten salt thermocline storage tank[J]. International Journal of Sustainable Energy, 2016, 35(6): 606-614. DOI:10.1080/14786451.2014.930465
[14]
HOFFMANN J F, FASQUELLE T, GOETZ V, et al. A thermocline thermal energy storage system with filler materials for concentrated solar power plants: experimental data and numerical model sensitivity to different experimental tank scales[J]. Applied Thermal Engineering, 2016, 100: 753-761. DOI:10.1016/j.applthermaleng.2016.01.110
[15]
ZANGANEH G, PEDRETTI A, ZAVATTONI S, et al. Packed-bed thermal storage for concentrated solar power–Pilot-scale demonstration and industrial-scale design[J]. Solar Energy, 2012, 86(10): 3084-3098. DOI:10.1016/j.solener.2012.07.019
[16]
BELLAN S, ALAM T E, GONZáLEZ-AGUILAR J, et al. Numerical and experimental studies on heat transfer characteristics of thermal energy storage system packed with molten salt PCM capsules[J]. Applied Thermal Engineering, 2015, 90: 970-979. DOI:10.1016/j.applthermaleng.2015.07.056
[17]
杨小平, 杨晓西, 丁静, 等. 高温填充床蓄热过程传热特性的实验研究[J]. 工程热物理学报, 2014, 35(3): 529-532.
[18]
SENTHIL R, SUNDARAM P, KUMAR M. Experimental investigation on packed bed thermal energy storage using paraffin wax for concentrated solar collector[J]. Materials Today:Proceedings, 2018, 5(2): 8916-8922. DOI:10.1016/j.matpr.2017.12.326
[19]
REGIN A F, SOLANKI S C, SAINI J S. An analysis of a packed bed latent heat thermal energy storage system using PCM capsules: numerical investigation[J]. Renewable Energy, 2009, 34(7): 1765-1773. DOI:10.1016/j.renene.2008.12.012
[20]
PENG H, DONG H H, LING X. Thermal investigation of PCM-based high temperature thermal energy storage in packed bed[J]. Energy Conversion and Management, 2014, 81: 420-427. DOI:10.1016/j.enconman.2014.02.052
[21]
RADY M. Thermal performance of packed bed thermal energy storage units using multiple granular phase change composites[J]. Applied Energy, 2009, 86(12): 2704-2720. DOI:10.1016/j.apenergy.2009.04.027
[22]
LI M J, JIN B, MA Z, et al. Experimental and numerical study on the performance of a new high-temperature packed-bed thermal energy storage system with macroencapsulation of molten salt phase change material[J]. Applied Energy, 2018, 221: 1-15. DOI:10.1016/j.apenergy.2018.03.156
[23]
LI M J, JIN B, YAN J J, et al. Numerical and Experimental study on the performance of a new two-layered high-temperature packed-bed thermal energy storage system with changed-diameter macro-encapsulation capsule[J]. Applied Thermal Engineering, 2018, 142: 830-845. DOI:10.1016/j.applthermaleng.2018.07.026
[24]
PARTOPOUR B, DIXON A G. Computationally efficient incorporation of microkinetics into resolved-particle CFD simulations of fixed-bed reactors[J]. Computers & Chemical Engineering, 2016, 88: 126-134.
[25]
YANG X P, CAI Z D. An analysis of a packed bed thermal energy storage system using sensible heat and phase change materials[J]. International Journal of Heat and Mass Transfer, 2019, 144: 118651. DOI:10.1016/j.ijheatmasstransfer.2019.118651
[26]
WHITAKER S. Forced convection heat transfer correlations for flow in pipes, past flat plates, single cylinders, single spheres, and for flow in packed beds and tube bundles[J]. AIChE Journal, 1972, 18(2): 361-371. DOI:10.1002/aic.690180219