1. Introduction
Offering a high energy density and nearly isothermal operation, a latent heat thermal energy storage (LHTES) system is often a preferable choice over other thermal energy storage systems. LHTES systems use a phase change material (PCM) as the storage medium. Thermal energy is stored when the solid PCM undergoes a phase change to liquid, and then when the phase change process is reversed, i.e., the molten PCM solidifies, the stored energy is released. The performance of an LHTES system might be degraded due to the low thermal conductivity of the PCMs; therefore, implementation of heat transfer enhancement techniques may become necessary. The heat transfer enhancement techniques can be divided into two categories: (1) improving the thermal physical properties of the PCM by dispersion of highly conductive micro- and nano-sized particles and (2) using highly conductive inserts such as metal fins, heat pipes, and porous structures [
1,
2,
3].
In recent years, the dispersion of highly conductive nanoparticles into the PCM has been utilized to improve the unit performance during the charging and discharging processes [
4]. Various nanoparticles such as pure metals [
5,
6], metal oxides [
7,
8], and carbon-based nanoparticles [
9] have been considered. The results of these studies have denoted that the volume fraction of nanoparticles plays an important role in controlling the melting and solidification processes. The impact of graphene nanoparticles on the melting process of various PCMs in a square cavity was studied by Kant et al. [
10]. It was found that using nanoparticles improves the conductive heat transfer while weakening the natural convection heat transfer. The overall results, however, indicated that there is a reduction in the melting time.
Lin and Al-Kayiem [
11] have experimentally evaluated the effects of incorporating Cu nanoparticles on the thermal characteristics of a paraffin wax during the melting and solidification processes. They have reported that the dispersion of copper nanoparticles with a concentration of 2.0% by weight leads to a 46.3% increase in the thermal conductivity of the composite PCM.
A numerical study performed by Mahdi and Nsofor [
12] indicated that the discharging time is reduced by 8–20% of alumina nanoparticles with volume concentrations of 8–13%, respectively, are used. The melting of paraffin in the presence of copper nanoparticles in an elliptical capsule was numerically studied by Boukani et al. [
13]. Three volumetric nanoparticle concentrations of 0%, 1%, and 3% and three container aspect ratios of 2, 1, and 0.5 were employed. The results revealed that, at a fixed aspect ratio, the presence of nanoparticles improves the performance of the system.
As mentioned earlier, fins improve the distribution of heat through high thermal conductive extended surfaces, and thus, they reduce a thermal resistance in an LHTES. Gharebaghi and Sezai [
14] studied the influence of aluminum fin arrays on the melting process of paraffin wax in a rectangular container. The results showed that the melting time and the stored energy were impacted by the spacing and the height of the fin arrays. The performance of a fin and tube LHTES system was experimentally investigated by Rahimi et al. [
15]. The heat transfer fluid (HTF) flowed inside a spiral tube, and plate fins were used to increase the heat transfer rate. The results showed that the presence of fins notably reduces the melting time, and the closer the fins, the shorter the melting time. Yang et al. [
16] simulated the phase change heat transfer in an annularly finned LHTES. The authors discussed the fin quantity, height, and thickness for optimal charging of a shell-tube LHTES unit. The authors considered the contribution of natural convection on the overall melting. It was found that the presence of fins diminishes the charging time by 65%.
In addition to plate and circular fins, the use of longitudinal fins has also been studied by researchers [
17,
18]. In a research study conducted by Kazemi et al. [
19], the effect of the longitudinal fin angle was evaluated. The results showed that a lower melting time was obtained when the angle was increased from 60° to 120° for the case with three fins. However, when two fins were used, increasing the angle from 45° to 150° deteriorated the performance of the system. The effect of fin angle was also investigated by Deng et. al. [
20]. In an experimental study carried out by Rathod and Banerjee [
21], the augmentation in heat transfer was studied in a shell and tube heat exchanger. Three longitudinal fins were installed at the outer surface of the inner tube to increase the heat transfer surface. The results showed that the decreased percentage in melting time relies on the heat transfer fluid (HTF) inlet temperature.
The available studies in the literature have reported that combining heat transfer enhancement techniques can lead to a better heat transfer performance of the storage unit [
22,
23]. Studies can be found using a combination of fins and heat pipes [
24], nanoparticles and fins [
25,
26], nanoparticles and metal foams [
27], metal foams, and heat pipes [
28]. Ren et al. [
29] studied the application of triangular fins in a rectangular LHTES with and without the presence of nanoparticles in the PCM. It was reported that using triangular fins was more efficient than solely using nanoparticles. In addition, the arrangement with long and narrow fins that were appropriately spaced provided a higher energy storage efficiency. Li et al. [
30] studied the melting of graphene nanoplatelets in 1-tetradecanol in a rectangular enclosure. They investigated the impact of various mass fraction of nanoparticles and enclosure aspect ratio on the melting behavior of nano-enhanced phase change material (NePCM). They found that a slender enclosure could be better for application of NePCMs.
Darzi et al. [
31] studied the influence of copper nanoparticles in three horizontal double pipes LHTES systems. The three configurations included a circular shell with an elliptical inner pipe, a cylindrical inner pipe, and a finned cylindrical inner pipe. The effect of the number of fins, as well as the orientation of the ellipse, were evaluated for different concentrations of nanoparticles. Mahdi and Nsofor [
22] studied the combined effects of nanoparticles and longitudinal fins in a shell and tube LHTES unit. Fins were added both on the outer surface of the inner pipe and the inner wall of the shell. It was found that better solidification time was achieved for the case with only fins when compared to the cases with only nanoparticles or with a combination of both.
The literature review of the available studies shows that natural convection heat transfer plays a role in the system’s performance, especially during the melting process. It should be noted that the significance of these effects is influenced by various parameters such as PCM type, container configuration, and container orientation. Shell and tube latent heat thermal energy storage systems are one of the common types of LHTES systems and are considered for the current study. The heat transfer fluid flows inside the inner tube, while the PCM is located in the shell. As the container is located vertically, the melting process will be governed by natural convection heat transfer. Therefore, due to the upward motion of the molten PCM, more melting will occur at the top of the container as compared to the lower section. To provide a more uniform and accelerated melting process, a conical shell is considered. To further accelerate the melting process, combinations of fins and dispersed nanoparticles are used as the heat transfer enhancement technique.
As mentioned earlier, fins contribute to the enhanced heat transfer process by providing more surface area. The level of improvement varies based on the design of the shell and the fins geometry. In some cases, the presence of the fins may dampen the natural convection flows, which may suppress some of the potential improvements. On the other hand, using dispersed nanoparticles increases the effective thermal conductivity of the NePCM, but it also increases the viscosity of NePCM, which consequently limits the liquid motion and natural convection. By considering these factors, finding an ideal configuration of an LHTES system and enhancement techniques is a crucial task. Thus, the present study aims to address the melting heat transfer in a conical shape LHTES and its optimum design for the shell’s aspect ratio, fins’ aspect ratio, type of nanoparticles, and volume fractions of nanoparticles for the first time.
2. Mathematical Model
A conical latent heat thermal energy storage (LHTES) system, filled with a nano-enhanced phase change material (NePCM), is illustrated in
Figure 1. The general shape of the LHTES’ shell is a conical shape that could be a descending or ascending conical enclosure. The fins can also be uniform, ascending, or descending. Based on the shape of the LHTES unit shell here, two types of the general design were adopted. A view of the ascending conical and the descending conical LHTES units are depicted in
Figure 1a,c, respectively. The cylindrical design is a special case of each of these two designs and is illustrated in
Figure 1b. Assuming the LHTES volume is constant, the radius of the cross-section of the system along the
z-direction can be descending, constant, or increasing. Furthermore, the height of the connected copper fins can be changed according to the shape of the system. The height of the LHTES is
L = 200 mm, and the radius of the normal system with a uniform circular cross-section (depicted in
Figure 1b) is
Ro = 50 mm. In the conical LHSS, the lower radius of the cross-section, i.e.,
Rdc, is defined base on the aspect ratio, i.e.,
ARc =
Rdc/
Ro. Also, the length of the lowest fin, i.e.,
wdf, is
ARf ×
wuf in which
wuf is the length of identical fins. It is worth noting that the volume of all fins is considered a constraint.
Figure 2 shows the structural specification of the unit. Water, as the heat transfer fluid has a high temperature, enters the tube and the energy is transferred to the composite NePCM through the tube wall. The NePCM contains the highly thermal conductive nano-additives added to the coconut oil as the base PCM. The thermophysical properties of the base PCM, HTF, and the fins are listed in
Table 1. As tabulated in
Table 2, the nano additives dispersed in the PCM can separately be Al
2O
3, GO, Ag, and Cu. Furthermore, the axisymmetric structure of the system allows the problem to be solved in 2D space.
For the NePCM domain, the enthalpy-porosity technique is used to model the melting flow. The controlling equations, i.e., the equations describing mass, momentum, and energy conservations, are:
where the velocity vector (
) and the pressure field (
p) are the dependent variables. The density (ρ), volume expansion coefficient (
β), and dynamic viscosity (
μ) are the thermophysical properties. The gravity direction (
) is toward the bottom, and its component in the horizontal direction is zero. The subscripts
LNePCM shows the liquid NePCM and
Am (O(10
6)) is a large number to procurance the value of the source term. The liquid fraction (
s) is a temperature dependent-function:
controls the velocity components in the momentum domain. The symbols
and
are the melting and window temperatures.
where temperature (
T) is the dependent variable. Subscripts
SNePCM and
PPCM refer to the properties of the solid NePCM, and the pure PCM (without nanoparticles). Thermal conductivity coefficient (
k), latent heat of phase transition (
L), and heat capacity (
Cp) are the thermophysical properties.
VFna represents the volume fraction of nanoparticles.
For the solid fin domain, the energy balance can be defined as the following:
where the subscript
s denotes the solid fin.
Assuming the fluid flow passing the tube is laminar and incompressible, the governing equations are:
The subscript HTF points to heat transfer fluid (water) in the tube. The thermophysical characteristics of the nano-enhance phase change material are reached based on the relations listed below:
PPCM denotes the properties of the pure phase change material.
Coefficient of thermal expansion
The relevant initial and boundary conditions for the above equations are as follows.
At the interface of the tube and NePCM domain:
At the entrance of the tube:
Also, the initial temperature of HTF, NePCM and fins is
where the subscripts
r and
z are the coordinates in vertical (axial) and horizontal (radial) directions, respectively. At the vertical outer wall of the NePCM domain:
At the horizontal walls of the NePCM domain:
The pertinent initial condition of the NePCM zone:
In the LHTES unit, the total stored energy is defined as the following:
The liquid fraction of NePCM during melting flow is expressed as:
is the area that NePCM is in liquid phase, and
is the non-melted NePCM area. To evaluate the uniformity of the temperature distribution, σ (
t) can be defined as follows:
indicates the average temperature of the PCM domain.
4. Results and Discussion
The present study aims to maximize the melting rate at seven hours of charging time. Here, the Taguchi optimization method is utilized to optimize four design parameters for the descending and ascending conical designs. The design parameters are the volume fraction of nanoparticles
(
), the ratio of the radius of the lower fin to the radius of the base fin (
for ascending conical and
for descending conical design), the ratio of the lower radius of conic to the radius of the regular cylinder (
,
), and type of nanoparticles. The type of nanoparticles could be Al
2O
3, GO (graphene oxide), Ag, and Cu. The amount of molten PCM is selected as the target design. Thus, “the higher, the better” design strategy was set for the Taguchi method. Each of the design parameters is divided into four levels. The details of the design parameters and their levels are summarized in
Table 4 and
Table 5. These tables show the design for ascending and descending conical shells, respectively.
The combination of four design parameters and four levels leads to 4
4 possible design combinations for each shell geometry. Simulating all possible combinations requires a tremendous computational budget. Thus, Taguchi uses an orthogonal table to reduce the required simulations. Here, the standard orthogonal L16 table has been selected for each design geometry.
Table 6 and
Table 7 show the L16 design cases for each of the ascending and descending conical shells, respectively.
The melting phase change simulations were executed for each of the design cases of L16
Table 6 and
Table 7. Then, the values of MVF and Total stored energy (
ES) after seven-hours of charging (melting) were reported in the Tables. The Taguchi method was used to compute the values of the S/N ratio based on the values of MVF. Since the target is “the higher, the better,” a larger value of S/N ratio shows a better-designed case. The Taguchi method uses these S/N ratios to find the best level of each design parameter and reach a final optimum design.
Figure 6 and
Figure 7 show the values of the S/N ratio for each level of each design parameter for two geometrical designs of ascending and descending conical shells. From
Figure 6, it can be found that the maximum S/N values correspond to
VFna = 0.045,
ARf = 0.4,
ARc = 1.0, and Al
2O
3 nanoparticles for ascending design. The optimum design levels for descending design could be found in
Figure 7, as follows:
VFna = 0.015,
ARf = 1.4,
ARc = 1.17, and Cu nanoparticles. The optimum design parameters and the estimated values of MVF by the Taguchi method are summarized in
Table 8. Numerical simulations were also executed for these optimum cases, and the computed MVF is reported in
Table 8. The simulated cases show an MVF = 0.92 after seven hours of thermal charging for both ascending and descending designs. These values are greater than equal to the simulated case of L16 tables. Thus, they have been selected as optimum designs.
It should be noted that the range of the variation of S/N in
Figure 6 and
Figure 7 shows the impact of a design parameter on the MVF. As seen, the most significant parameter is
ARc for both cases. This parameter indicates the ratio of the lower radius of conic to the radius of the regular cylinder. The next important design parameter is
ARf, which represents the ratio of the radius of the lower fin to the radius of the base fin. The other less important design variables are the type of nanoparticles and the volume fraction of nanoparticles. Thus, special attention should be made to selecting shell ratio (
ARc) and fin ratio (
ARf). Then, the LHTES unit’s charging rate can be fine-tuned by selecting the type and volume fraction of nanoparticles.
For ascending conical shell geometry, the minimum
MVF corresponds to case 16
Table 6 with
MVF = 0.75. Thus, the optimum case with
MVF = 0.92 increases the melting rate by 18.5% compared to this case. The lowest
MVF = 0.76 corresponds to case 13 of
Table 7 for a descending design. Since the optimum case is
MVF = 0.92, a 17.3% increase of
MVF can be achieved by an optimum design.
Figure 8 shows the melting maps of NEPCM for the optimum case of ascending conical shell. The results are reported at seven different time steps during the melting process.
Figure 9 depicts the isotherm maps for the same snapshots of
Figure 6. As seen, the melting commences around the tube wall and fins. Then, it advances in the enclosure toward the adiabatic shell. At t = 3750 s, small circulation cells between fins can be seen in molten regions. The circulation extends to a large portion of the enclosure by the advancement of melting, and a single central circulation occurs at t = 15,000 s and beyond. The solid zones quickly disappear at the top of the enclosure due to the circulation flows and long fins at the top. The only un-melted region at 26,250 s is at the bottom of the enclosure. This region slowly melts down as time passes.
Figure 9 shows that the HTF in the tube is mostly at a uniform hot temperature. Only at the initial charging stages, a smooth low-temperature boundary layer next to the tube wall can be seen. At 7500 s and larger, the heat transfer rate at the PCM sides drops notability, and the temperature of HTF in the tube is uniformly hot. It is interesting that the temperature of the molten PCM reaches the HTF temperature at the top regions of the enclosure, where there is a small amount of solid PCM at the bottom.
Figure 10 and
Figure 11 show melting fraction maps and the temperature in the optimum design of an ascending conical shell. As seen in this design, the fins are long at the bottom and short at the top. The shell is also wide at the bottom and narrow at the top. In the descending case, the same as the ascending one, the melting commences from the tube wall and around the fins and advances into the enclosure toward the shell. However, in this design, there is a narrow space between the tube and con-shell at the top. This narrow space breaks the main circulation flow into two distinct circulation flows. The PCM melts uniformly next to the cone shell. Thus, a smaller piece of solid PCM can be seen at the bottom at t = 22,500 s compared to the ascending design. Moreover, the hot molten PCM is in the middle of the enclosure and advances toward the corners.
Figure 12 shows the temperature variation of eight monitor points in each of the ascending and descending designs. The monitor points are depicted in
Table 9. The points start from the bottom and move upward. As seen, all of the temperature points are initially at a supper cold temperature of 293 K. Then, the temperature sharply increases to about 295 K. This sharp rise of temperature corresponds to the solid region’s conduction heat transfer. Then the phase change occurs, and the rise of the temperature limits to the phase change temperature bond. In this situation, the heat transfer does not increase the temperature notably, and the heat stores in the PCM in the form of latent heat. Once the PCM at the monitor point melts down, the temperature rises again. The rise of temperature at the end of the melting process is due to the increase of temperature in the molten PCM. The monitor points at the top of the enclosure tend to quickly reach about 313 K, which is the temperature of HTF. The monitor points show a short phase change history for the case of ascending design. The melting maps also showed that the bottom of the enclosure remains at solid state during most of the charging time.
Figure 12b shows the monitor points have almost similar behavior since the melting front advances in the enclosure uniformly toward the enclosure shell.
Figure 13a–c show the melted portion (MVF), total stored energy (ES), and temperature non-uniformity (σ) during the melting process, respectively. The results are plotted for two ascending and descending designs. As seen, the descending designs lead to a higher MVF during the charging process.
Only at the final stages of melting, the ascending design is better than the descending design. Thus, for applications in which the full charging is not very important and about 85% of charging is adequate, then the descending design can be selected. However, when the full melting is the main target of design, the ascending design should be selected.
Figure 13b shows that the total stored energy for the case of ascending design is also higher around the final stage of charging (just before full charging). The temperature non-uniformity for the ascending design is higher than the case of the descending design at the end of the charging process. This is because the top of the ascending design at the end of the charging process is at a hot temperature and subject to natural convection while its bottom is not melted yet. Thus, the LHTES unit is subject to high-temperature differences.
The results of
Figure 6 show that the presence of nanoparticles could induce a minor impact on MVF, and the increase of the concentration of nanoparticles does not improve MVF monastically. These findings are in agreement with the results of the experimental investigation of Li et al. [
30]. The experimental results show that the presence of nanoparticles increases both thermal conductivity and dynamic viscosity. Thus, using NePCMs could benefit conduction-dominant regimes where the heat transfer is under the direct influence of thermal conductivity while the dynamic viscosity is not important. However, as the melting advances and natural convection flows occur the dynamic viscosity influences the intensity of liquid circulation. Thus, the growth of viscosity due to the presence of nanoparticles could slow down the natural convection effects. It should also be noted that the natural convection effects are under the influence of geometrical design (enclosure shape and fins), and therefore, the combination of these effects should be taken into account in the heat transfer design of LHTES systems.
Currently, various types of NePCMs are being synthesized by various researchers. Some of these works have been reviewed by Leong et al. [
43] and Nižetić [
44]. The analysis of synthesized NePCM samples shows that the presence of nanoparticles changes the host PCM’s thermophysical properties. It should be noted that the heat transfer performance of LHTES systems is a nonlinear function of the thermophysical properties and their geometrical design. Thus, numerical simulations/experimental studies are demanded to investigate the advantage/drawback of using NePCMs. The present study took into account the combination of geometrical factors and the presence of nanoparticles to analyze the heat transfer behavior of LHTES systems for various concentrations and types of nanoparticles.