A numerical Procedure for PCM Thermal Storage Design in Solar Plants

A numerical Procedure for PCM Thermal Storage Design in Solar Plants

Cammarata G. Monaco L. Cammarata L. Petrone G.

Department of Industrial Engineering, University of Catania, Viale A. Doria 6, 95123 Catania, Italy

Corresponding Author Email: 
gpetrone@dii.unict.it
Page: 
105-110
|
DOI: 
https://doi.org/10.18280/ijht.310214
Received: 
N/A
| |
Accepted: 
N/A
| | Citation

OPEN ACCESS

Abstract: 

This paper is aimed to propose an analytical and numerical procedure for design thermal storage systems employing Phase Change Materials PCM as medium. The use of PCM for thermal energy storage is a really challenging opportunity because of their high heat capacity. It is assessed that heat capacity per unit volume of a PCM medium can be 5-14 times more than traditional sensible heat storage materials such as water, masonry or rock. In addition, PCM absorbs and release heat at a nearly constant temperature. Therefore, PCM represent a valuable choice for those energy systems needing a thermal storage section, such as solar thermo-dynamic plants. Latent heat storage is based on the heat absorption or release when the storage material undergoes a phase change from solid to liquid or vice versa. A complete understanding of the phase change phenomenon involves an analysis of the various processes that accompany it. The most important of these processes, from a macroscopic point of view, is the heat transfer process. This is complicated by the release, or absorption, of the latent heat of fusion at the "moving" solid-liquid interface. A numerical procedure to define design criteria for engineered thermal storage units is proposed in this study. It is mainly based on the following steps: i) assessment of the total mass/volume of a chosen PCM by the energy balance solution; ii) analytical computation of the transient position of the phase front by a transcendental equation solution applied to outward-directed melting/freezing of a hollow cylinder due to a temperature/thermal flux imposed at the internal radius; iii) determination of the geometrical properties of the storage unit (length, number and relative arrangement of tubes in the array); iv) FE-based simulation of a portion of the storage unit for estimating the total transient time needed by the melting/freezing process and the energy effectively stored by the designed system. Numerical simulations were carried-out adopting the enthalpy formulation for the energy equation, that allows to solve one single equation for heat transfer both in solid and liquid phase. The applied methodology was firstly validated by comparison with literature results. As test-case for the numerical tool built-up, we analysed the melting process of a Lithium Nitrate storage unit serving a 1 MW Direct Steam Generation solar plant. The present work has been developed in the framework of the FREeSun, supported by a grant from the Italian Ministry of the Economic Development.

  References

[1] Lewis, R.W., Nithiarasu, P., Seetharamu, K.N. (2004) Fundamentals of the Finite Element Method for Heat and Fluid Flow. John Wiley & Sons

[2] Henze, R.H., Humphrey, J.A.C. Enhanced heat conduction in phase-change thermal energy storage devices (1981) International Journal of Heat and Mass Transfer, 24 (3), pp. 459-474. doi: 10.1016/0017-9310(81)90054-5

[3] Ramachandran, N., Gupta, J.P., Jaluria, Y. Thermal and fluid flow effects during solidification in a rectangular enclosure (1982) International Journal of Heat and Mass Transfer, 25 (2), pp. 187-194. doi: 10.1016/0017-9310(82)90004-7

[4] Sparrow, E.M., Larson, E.D., Ramsey, J.W. Freezing on a finned tube for either conduction-controlled or natural-convection-controlled heat transfer (1981) International Journal of Heat and Mass Transfer, 24 (2), pp. 273-284. doi: 10.1016/0017-9310(81)90035-1

[5] Kemink, R.G., Sparrow, E.M. Heat transfer coefficients for melting about a vertical cylinder with or without subcooling and for open or closed containment (1981) International Journal of Heat and Mass Transfer, 24 (10), pp. 1699-1710. doi: 10.1016/0017-9310(81)90078-8

[6] Crank, J. (1984) Free and Moving Boundary Problems. Oxford University Press

[7] Bonacina, C., Comini, G., Fasano, A., Primicerio, M. Numerical solution of phase-change problems (1973) International Journal of Heat and Mass Transfer, 16 (10), pp. 1825-1832. doi: 10.1016/0017-9310(73)90202-0

[8] Alexiades, V., Solomon, A.D.(1993) Mathematical Modelling of Melting and Freezing Processes. Hemisphere publishing corporation

[9] Rao, P.R., Sastri, V.M.K.Efficient numerical method for two-dimensional phase change problems (1984) International Journal of Heat and Mass Transfer, 27 (11), pp. 2077-2084. doi: 10.1016/0017-9310(84)90194-7

[10] Voller, V.R. Implicit finite-difference solutions of the enthalpy formulation of stefan problems (1985) IMA Journal of Numerical Analysis, 5 (2), pp. 201-214. doi: 10.1093/imanum/5.2.201

[11] Kim, C.-J., Kaviany, M. A numerical method for phase-change problems (1990) International Journal of Heat and Mass Transfer, 33 (12), pp. 2721-2734. doi: 10.1016/0017-9310(90)90206-A

[12] Mampaey, F. A stable alternating direction method for simulating multi‐dimensional solidification problems (1990) International Journal for Numerical Methods in Engineering, 30 (4), pp. 711-728. doi: 10.1002/nme.1620300411

[13] Voller, V.R. Fast implicit finite-difference method for the analysis of phase change problems (1990) Numerical Heat Transfer, Part B: Fundamentals, 17 (2), pp. 155-169. doi: 10.1080/10407799008961737

[14] Lee, S.L., Tzong, R.y. An enthalpy formulation for phase change problems with a large thermal diffusivity jump across the interface (1991) International Journal of Heat and Mass Transfer, 34 (6), pp. 1491-1502. doi: 10.1016/0017-9310(91)90291-L

[15] Raw, W.Y., Lee, S.L. Application of weighting function scheme on convection-conduction phase change problems (1991) International Journal of Heat and Mass Transfer, 34 (6), pp. 1503-1513. doi: 10.1016/0017-9310(91)90292-M

[16] Rabin, Y., Korin, E. An efficient numerical solution for the multidimensional solidification (or melting) problem using a microcomputer (1993) International Journal of Heat and Mass Transfer, 36 (3), pp. 673-683. doi: 10.1016/0017-9310(93)80043-T

[17] Clavier, L., Arquis, E., Caltagirone, J.P., Gobin, D. A fixed grid method for the numerical solution of phase change problems (1994) International Journal for Numerical Methods in Engineering, 37 (24), pp. 4247-4261. doi: 10.1002/nme.1620372408

[18] Lamberg, P. Approximate analytical model for two-phase solidification problem in a finned phase-change material storage (2004) Applied Energy, 77 (2), pp. 131-152. http://www.elsevier.com/inca/publications/store/4/0/5/8/9/1/index.htt doi: 10.1016/S0306-2619(03)00106-5

[19] Ismail, K.A.R., Gonçalves, M.M.Thermal performance of a pcm storage unit (1999) Energy Conversion and Management, 40 (2), pp. 115-138. doi: 10.1016/S0196-8904(98)00042-9

[20] Petrone, G., Cammarata, G. Numerical simulation of phase change materials melting process (2013) Computational Thermal Sciences, 5 (3), pp. 227-237. http://www.dl.begellhouse.com.ezproxy3.lhl.uab.edu/download/article/6be7d6f2021c84f7/CTS0503%286%29-6360.pdf doi: 10.1615/ComputThermalScien.2013006360

[21] Wang, Y., Amiri, A., Vafai, K. An experimental investigation of the melting process in a rectangular enclosure (1999) International Journal of Heat and Mass Transfer, 42 (19), pp. 3659-3672. doi: 10.1016/S0017-9310(99)00024-1