Domestic Hot Water Storage Tank Utilizing Phase Change Materials (PCMs): Numerical Approach

Abstract

Thermal energy storage (TES) is an essential part of a solar thermal/hot water system. It was shown that TES significantly enhances the efficiency and cost effectiveness of solar thermal systems by fulfilling the gap/mismatch between the solar radiation supply during the day and peak demand/load when sun is not available. In the present paper, a three-dimensional numerical model of a water-based thermal storage tank to provide domestic hot water demand is conducted. Phase change material (PCM) was used in the tank as a thermal storage medium and was connected to a photovoltaic thermal collector. The present paper shows the effectiveness of utilizing PCMs in a commercial 30-gallon domestic hot water tank used in buildings. The storage efficiency and the outlet water temperature were predicted to evaluate the storage system performance for different charging flow rates and different numbers of families demands. The results revealed that increases in the hot water supply coming from the solar collector caused increases in the outlet water temperature during the discharge period for one family demand. In such a case, it was observed that the storage efficiency was relatively low. Due to low demand (only one family), the PCMs were not completely crystallized at the end of the discharge period. The results showed that the increases in the family’s demand improve the thermal storage efficiency due to the increases in the portion of the energy that is recovered during the nighttime.

1. Introduction

Phase change materials (PCMs) are an important topic in research and industry. Due to the intermittent supply characteristics of solar energy, energy storage represents a potential solution by storing the energy during the daytime to use it during the nighttime. PCMs can store large amounts of heat at constant temperature while the material changes phase or state. When PCMs are coupled with solar technology, efficiency of traditional heating, ventilation, and air conditioning (HVAC) systems may potentially increase. The majority of applications for PCMs are for space heating/cooling and providing domestic hot water for buildings. Also, PCMs have high energy density and latent heat, and they are cost-effective. There are several applications in which PCMs were implemented as shown in

Table 1. Phase change material (PCM) applications.

Application	References
Thermal storage of solar energy	[5,6,7]
Heating and sanitary hot water	[8,9]
Cooling	[10,11,12,13,14,15,16]
Thermal comfort in vehicles	[17]
Solar power plants	[18,19,20,21]
Cooling of engines	[22]
Thermal protection of electronic devices	[23]
Spacecraft thermal system	[24]

.

Phase change material (PCM) is an environmentally friendly material used to improve building energy consumption and indoor thermal comfort [1,2].

An experimental study investigated a heat pump utilizing a thermal energy storage (TES) tank [3]. In their research, it was found that a PCM storage tank has 14.5% better performance. Moreover, the PCM storage tank improved indoor temperature stability within comfort 20.65% longer in time compared to the conventional water tank. Also, a thermal storage system was installed in a one family house [4], where sodium acetate trihydrate (SAT) was used as a PCM.

The performance of thermal storage systems was analyzed [25]. A south-oriented wall was used as thermal storage with phase change materials embedded in the wall. The thermal storage system can store solar radiation up to 6–8 h after solar irradiation; this has effects on the stability of the daily temperature swings (up to 10 °C).

Zalba et al. [26] conducted an experimental study to store outdoor cold during nighttime and to release it indoors during daytime using PCMs. The results revealed that the system was successfully tested for 1000 cycles without any degradation in the performance of the system.

Kenneth [27] developed a solar system that utilizes a PCM in domestic houses in the United Kingdom (UK). The system consisted of solar flat plate collectors, which put the energy into a storage tank and PCM-filled panels. Calcium chloride was used as the PCM with a melting point of 29 °C. The results revealed that the use of PCMs reduced energy consumption by 18–32%.

The literature shows a limited number of investigations on PCM thermal storage tanks coupled with a solar thermal collector to provide domestic hot water or building heating. In the present paper, a numerical model was developed for a 30-gallon hot water tank utilizing n-eicosane PCM as the thermal store medium. The solar energy was used to charge the storage tank during the daytime (9 h) and then the thermal energy was recovered to provide hot domestic water during the discharge period (15 h).

2. Numerical Model Description

In this study, a three-dimensional numerical model was created using finite element techniques (COMSOL Multiphysics) for a 30-gallon domestic hot water thermal storage tank connected to a photovoltaic thermal collector. The tank height and diameter were 1.15 m and 0.46 m, respectively, as shown in Figure 1. PCMs are present in small cylindrical containers inside the tank as shown in Figure 2a. The tank contains a spiral heat exchanger to exchange the heat from the circulating water to the thermal storage medium as shown in Figure 2b. The diameter of the heat exchanger pipe was 0.02 m, and the length of the heat exchanger was 1.1 m with a spiral diameter of 0.28 m. The length and diameter of the PCM containers were 0.2 m and 0.05 m, respectively. The tank consisted of 32 PCM cylindrical containers per row with a total of five rows. The properties of the PCM used in the model are presented in

Table 2. PCM properties.

n-Eicosane	Solid Phase	Liquid Phase
Density (kg/m3)	856	778
Specific heat (kJ/(kg∙K))	2.136	2.1336
Thermal conductivity (W/m∙K)	0.35	0.15
Melting point	36.4 °C
Latent heat (kJ/kg)	247.3

.

2.1. Governing Equations

The assumptions on which the governing equations were based prior to the formulation of the model were as follows:

(1)

The fluid flow is Newtonian and incompressible;

(2)

No heat is generated inside tank solid domains;

(3)

The variation of thermo-physical properties of the PCM can be neglected.

When these assumptions were taken into consideration, a system of governing equations to describe the heat transfer and fluid flow were solved and coupled. In the present study, an attempt was made to solve the Navier–Stokes equations for the fluid flow through the internal spiral heat exchanger and a free convection fluid flow inside the tank. In addition, the energy equations for all domains including the tank solid domain, the PCM domain, and the fluid domain were solved. It is important to note that the Navier–Stokes equations and the energy equations were coupled. The coupling was established for each time step using the velocity field obtained from Navier–Stokes equations as an input to evaluate the convective heat transfer term in the energy equation.

The Navier–Stokes equation and the energy equation were solved numerically using COMSOL Multiphysics [3] as follows:

Momentum equation along x-direction:

$$\rho \left(u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}+w\frac{\partial u}{\partial z}\right)=-\frac{\partial p}{\partial x}+\mu \left(\frac{{\partial }^{2}u}{\partial x^2}+\frac{{\partial }^{2}u}{\partial y^2}+\frac{{\partial }^{2}u}{\partial z^2}\right);$$

(1)

Momentum equation along y-direction:

$$\rho \left(u\frac{\partial v}{\partial x}+v\frac{\partial v}{\partial y}+w\frac{\partial v}{\partial z}\right)=-\frac{\partial p}{\partial y}+\mu \left(\frac{{\partial }^{2}v}{\partial x^2}+\frac{{\partial }^{2}v}{\partial y^2}+\frac{{\partial }^{2}v}{\partial z^2}\right);$$

(2)

Momentum equation along z-direction:

$$\rho \left(u\frac{\partial w}{\partial x}+v\frac{\partial w}{\partial y}+w\frac{\partial w}{\partial z}\right)=-\frac{\partial p}{\partial z}+\mu \left(\frac{{\partial }^{2}w}{\partial x^2}+\frac{{\partial }^{2}w}{\partial y^2}+\frac{{\partial }^{2}w}{\partial z^2}\right)+\rho g.$$

(3)

The continuity equation for this simulation can be expressed as

$$\left(\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}+\frac{\partial w}{\partial z}\right)=0.$$

(4)

The energy equation was as follows:

$${\left(\rho c_p\right)}_{}\left(u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}+w\frac{\partial T}{\partial z}\right)=\mathrm{K} \left(\frac{{\partial }^{2}T}{\partial x^2}+\frac{{\partial }^{2}T}{\partial y^2}+\frac{{\partial }^{2}T}{\partial z^2}\right),$$

(5)

where $\rho$ represents the water density, c_p represents the fluid specific heat, p represents the pressure, u, v, and w represent the coordinates of a velocity field vector, T represents the temperature, $\mu$ represents the dynamic viscosity of the fluid, and K represents the thermal conductivity.

For the phase change material,

$$c_p=\frac{1}{\rho }\left(\theta {\rho }_{phase1}{c}_p{}_{phase1}+\left(1-\theta \right){\rho }_{phase2}{c}_p{}_{phase2}\right)+L·\frac{\partial {\alpha }_m}{\partial T},$$

(6)

$${\alpha }_m=\frac{1}{2}\frac{\left(1-\theta \right){\rho }_{phase2}-\theta {\rho }_{phase1}}{\theta {\rho }_{phase1}+\left(1-\theta \right){\rho }_{phase2}},$$

(7)

where $\theta$ represents the PCM solid fraction, L represents the latent heat of the phase change material, phases 1 and 2 represent the solid and liquid phases, respectively, and $\frac{\partial {\alpha }_m}{\partial T}$ represents the melting fraction per degree of temperature.

2.2. Boundary Conditions

Boundary conditions can be divided into two types: heat transfer boundary conditions and fluid flow boundary conditions. The thermal boundary conditions involve the inlet water temperature to the heat exchange, which was 50 °C during the charging period over the daytime (9 h) and 20 °C during the discharge period over the nighttime (15 h). Moreover, the thermal boundary conditions included the insulated wall of the tank outer surface (∇.T = 0).

The fluid flow boundary conditions included the inlet velocity (U) at the inlet portion (assuming a flat profile), with no pressure constraints at the outlet portion and walls (e.g., no slip condition) at the remainder of the surface (see Figure 2a). It is important to note that, in the present study, the flow rates during the charging period were set to 2 L/m, 3 L/m, and 4 L/m. In addition, the flow rate of the demand of hot water during the discharge period was set to match the typical hot water demand for one, two, three, and four families, as shown in Figure 3 [28]. A typical family consists of two adults and two children [28]. As shown in Figure 3, the demand was low over the first 5 h of the day and increased dramatically over the rest of the day.

3. Results and Discussion

3.1. Effect of Charging Flow Rates

In this section, the flow rates over the charging period (9 h) were set to 2 L/m, 3 L/m, and 4 L/m. The demand of domestic hot water flow rate over the discharge period (during the remaining 15 h) was fixed to match the typical hot water demand for one family, as illustrated in Figure 3. Figure 4 shows the storage tank inlet and outlet temperatures versus time.

As shown in Figure 4, the inlet temperature was 50 °C during the charging period (9 h), which came from the photovoltaic thermal (BIPVT/T) collector. Then, the inlet water temperature decreased to 20 °C during the discharge period. The amount of thermal energy was calculated as follows:

$$q_{stored}=m·c_p\left(T_{in}-T_{out}\right)\Delta t;$$

(8)

$$q_{discharged}=m·c_p\left(T_{out}-T_{in}\right)\Delta t.$$

(9)

It was noted that, over the charging period, the amount of heat stored ($q_{stored}$) increased upon increasing the hot water supply flow rate at given demand.

As a result, the average outlet water temperatures during the discharge period increased with the increase in the hot water flow rate during the 9-h charging period, as shown in Figure 4.

Table 3. Tank storage efficiency.

Hot Water Supply Flow Rate	Heat Input (kJ)	Heat Extracted (kJ)	Storage Efficiency (%)
2 L/m	22,307.5	7869	35%
3 L/m	21,417	8253	38.5%
4 L/m	21,738	8511	39%

shows the amount of heat being stored and extracted during charging and discharging periods, respectively. It was observed that the maximum storage efficiency was 39%. The thermal storage efficiency is defined as a ratio between the energy extracted and the energy injected into the tank. This means that only 35–39% of stored energy was successfully extracted during the discharge period. Figure 5 and Figure 6 explain the reason for such low storage efficiency. Figure 5 shows the melting fraction of the PCM cylinders; when it is equal to one, it means that the PCMs are completely melted. As shown in this figure, PCMs were not melted completely after the 9-h charging period. This means that there was a capability of the storage tank to absorb more heat during the charging period.

Moreover, Figure 6 shows that the PCMs were not solidified/frozen completely after the discharge period (24 h). It was observed that about 50% of PCMs were not completely solidified/frozen. This means that the tank still had some available heat to be used. To reach better efficiency, the hot water demand should be increased in order to absorb the remaining stored heat, allowing the PCM to be completely solidified/frozen.

In order to better understand the melting phenomenon, Figure 7 shows the melting fraction of a single PCM cylinder. As shown in this figure, around half of the cylinder was melted completely after the charging period. Moreover, after completing one cycle (e.g., 24 h), the PCM cylinder still had a liquid phase up to 0.7, as shown in Figure 7b.

Figure 8a shows the mid-plane temperature contours at the end of charging period. As shown in this figure, the temperature of the circulating water coming from the solar collector was about 50 °C, and the average temperature of the tank was about 34 °C. Also, Figure 8b shows the temperature contours after 24 h, where the average temperature was about 30 °C.

3.2. Effect of Number of Families

In order to enhance the storage efficiency, domestic hot water demand was increased to match the typical demand of two, three, and four families. Figure 9a shows the outlet temperatures from the tank versus time for the demand of one, two, three, and four families and 2 L/m of hot water supply during the charging period. It was noted that the outlet temperature decreased with the increase in the number of families. The increase in hot water demand during the discharge periods (number of families) meant increases in the thermal energy to be extracted/recovered from the thermal storage tank, as shown in Figure 9b. However, for the demand of four families, it was observed that the domestic hot water temperature over the discharge period decreased dramatically until it reached below 30 °C. In such a case, a supplement back-up system is needed to maintain the outlet domestic water above 30 °C.

In order to quantify the thermal storage efficiency,

Table 4. Thermal storage efficiencies for multiple families.

Hot Water Supply Flow Rate	Number of Families	Heat Input (kJ)	Heat Extracted During Night (kJ)	Storage Efficiency (%)
2 L/m	1	22,307.5	7869	35%
2 L/m	2	22,307.5	13,415.96	60%
2 L/m	3	22,307.5	16,643.48	74%
2 L/m	4	22,307.5	18,288.89	82%

shows the amount of heat that was injected into the thermal storage tank during the charging period (9 h) and the amount of heat that was absorbed/recovered during the discharge periods. It was observed that, for a given hot water supply, increasing the number of families increased the efficiency from 35% for one family to 82% for four families. The four families used the whole store of energy and, because of the high demand rate, a big portion of the stored energy was recovered. As a result, the storage efficiency was increased up to 82%.

4. Conclusions

This paper presented a CFD numerical code of a domestic hot water tank utilizing a phase change material as a storage medium. The following conclusions were made:

• The increases in the hot water supply during the charging periods increased the storage efficiency from 35% to 39%.
• At given hot water supply, increasing the number of families increased the efficiency from 35% for one family to 82% for four families.
• At given hot water supply, the heat extracted over the nighttime increased from 7869 kJ to 18,288.89 kJ upon increasing the demand from one family to four families.

Further developments and future work on the present topic involve the introduction of nanofluid to enhance the thermal storage efficiency of the tank, the use of different PCMs and melting temperatures, and the calculation of the electricity consumption to determine the energy efficiency. Moreover, the effect of inlet water profile (cycling temperature) will be investigated as an extension of the present work.