Evaluation of the Melting Performance in a Conical Latent Heat Thermal Unit Having Variable Length Fins

: A conical shell-tube design with non-uniform ﬁns was addressed for phase change latent heat thermal energy storage (LHTES). The shell was ﬁlled with nano-enhanced phase change material (NePCM). The cone aspect ratio of the shell and the ﬁns aspect ratio were adopted as the geometrical design parameters. The type and volume fraction of the nanoparticles were other design parameters. The investigated nanoparticles were alumina, graphite oxide, silver, and copper. The ﬁnite element method was employed to solve the natural convection ﬂow and phase change thermal energy equations in the LHTES unit. The Taguchi optimization method was utilized to maximize the melting rate in the unit. Two cases of ascending and descending conical shells were investigated. The outcomes showed that the shell-aspect ratio and ﬁn aspect ratio were the most important design parameters, followed by the type and concentration of nanoparticles. Both ascending and descending designs could lead to the same melting rate at their optimum design. The optimum design of LHTES could improve the melting rate by up to 18.5%. The optimum design for ascending (descending) design was a plain tube (a cone aspect ratio of 1.17) ﬁlled by 4.5% alumina-Bio-PCM (1.5% copper-Bio-PCM).


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.

Mathematical Model
A conical latent heat thermal energy storage (LHTES) system, filled with a nanoenhanced 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 R o = 50 mm. In the conical LHSS, the lower radius of the cross-section, i.e., R dc , is defined base on the aspect ratio, i.e., AR c = R dc / R o . Also, the length of the lowest fin, i.e., w df , is AR f × w uf in which w uf 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 nanoadditives 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 2 O 3 , GO, Ag, and Cu. Furthermore, the axisymmetric structure of the system allows the problem to be solved in 2D space.   Table 1. Thermophysical properties of the coconut oil and the nano-additives [32,33].

Properties
Bio-Based PCM Water as HTF Solid Fin Liquid Solid 103e3 Non-defined Non-defined Non-defined Table 2. Nanoparticles dispersed in the phase change material [34][35][36].
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 ( → V) and the pressure field (p) are the dependent variables. The density (ρ), volume expansion coefficient (β), and dynamic viscosity (µ) are the thermophysical properties. The gravity direction ( → g ) is toward the bottom, and its component in the horizontal direction is zero. The subscripts LNePCM shows the liquid NePCM and A m (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 T melting and T window 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 (C p ) are the thermophysical properties. VF na 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: Density: PPCM denotes the properties of the pure phase change material. Dynamic viscosity: Coefficient of thermal expansion Thermal conductivity: Effective heat capacity: 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: At the outlet tube: Also, the initial temperature of HTF, NePCM and fins is T in = 293 K 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: A LNePCM is the area that NePCM is in liquid phase, and A SNePCM is the non-melted NePCM area. To evaluate the uniformity of the temperature distribution, σ (t) can be defined as follows: T average indicates the average temperature of the PCM domain.

Solution Method
To find the fields of MVF, temperature, and velocity, we employed the numerical analysis by user defined finite element codes based on the Petrov-Galerkin approach. The governing equations were transformed into a weak form, and linear shape function (for heat equation) and quadratic Lagrange (for momentum equation) were applied. Then, the equations were integrated over mesh elements to obtain the residual equations. The residual equations were solved in the fully coupled form using the PARallel DIrect SOlver (PARDISO) [37,38] with a relative tolerance of 1 × 10 −4 . A Newtonian damping factor of 0.9 was applied to enhance the convergence of the solver. The selection of the time-step is crucial in phase change problems with natural convection effects. Here, a backward differentiation formula (BDF) was used to control the computations' accuracy and convergence automatically. The details of the utilized finite element method and solution method can be found in [39,40].
To reach the high accuracy results, the mesh density was optimized through sensitivity analysis of the numerical results to the element numbers. Here a size parameter, NE, was introduced to control the size of mesh at various regions. The NE parameter was changed in the range of 0.75-1.75. A larger NE shows more elements and thus a fanner mesh. The details of the mesh elements in each region have been reported in Table 3. The mesh type in all domains is quad, and the mesh in the fins is structured, leading to rectangular mesh in the solid-fins domain. The reset of domains was filled with un-structured quad elements. The results of optimization can be found in Figure 3. As seen, the mesh with NE = 1.25 can be selected as the optimized grid.

Validation of the Model with Literature Data
The implemented numerical model can be verified by employing a numerical analysis performed by Costa [41] and an experimental work conducted by Kumar et al. [42]. Costa [41] studied the influence of the thermal conductivity of the separators on the free convection in a partially divided cavity. Figure 4 shows the temperature distribution of the current model and that simulated by Costa [41]. The working fluid filling the cavity is air with a Prandtl number of 0.71. The left and right walls were subject to cold and hot isothermal temperatures. The domain temperatures were scaled for the temperature difference between the hot and cold walls. Thus, the isotherms are in the range of 0-1.
The natural convection Ryleigh number was 10 6 . The partition walls were ten times more thermally conductive than the fluid inside the enclosure (air), and the size of the partition was 0.6 size of the enclosure, and the distance between the partition and side wall was 0.3 enclosure size. In another validation, the liquid fraction field observed by Kumar et al. [42] which was also simulated by the current model are compared in Figure 5. The authors investigated the melting of lead in an enclosure with a height of 50 mm and an aspect ratio of 0.75 (height/width). The enclosure walls were well insulated except the right wall, which was subject to a constant heat flux (16.3 kW/m 2 ). The thermal neutron radiography method was used to map the melting interface. As seen, the melting starts from the right wall and advances toward the left. Moreover, due to the natural convection effects, the melting at the top is faster than at the bottom. In this comparison, the right side-wall of the cavity occupied by lead was exposed to a constant heat flux. The comparative analyses presented in Figures 4 and 5 confirm the reliability of the utilized code.

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 VF na (0 ≤ VF na ≤ 0.045), the ratio of the radius of the lower fin to the radius of the base fin (0.1 ≤ AR f ≤ 1.0 for ascending conical and 1.0 ≤ AR f ≤ 1.6 for descending conical design), the ratio of the lower radius of conic to the radius of the regular cylinder (0.4 ≤ AR c ≤ 1.0, 1.0 ≤ AR c ≤ 1.5), and type of nanoparticles. The type of nanoparticles could be Al 2 O 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 Tables 4 and 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. Tables 6 and 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 Tables 6 and 7. Then, the values of MVF and Total stored energy (ES) after sevenhours 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.  6 and 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 VF na = 0.045, AR f = 0.4, AR c = 1.0, and Al 2 O 3 nanoparticles for ascending design. The optimum design levels for descending design could be found in Figure 7, as follows: VF na = 0.015, AR f = 1.4, AR c = 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 Figures 6 and 7 shows the impact of a design parameter on the MVF. As seen, the most significant parameter is AR c 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 AR f , 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 (AR c ) and fin ratio (AR f ). 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.   Figures 10 and 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.  Table 9. The coordinate of temperature monitor points in optimum ascending and descending designs.
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.

Conclusions
The melting heat transfer and thermal energy storage in a conical-shell LHTES unit filled with nano-enhanced phase change materials were addressed. The melting heat transfer was enhanced by applying a non-uniform array of fins. Two possible shell-designs, ascending and descending designs, were investigated. Four types of nanoparticles were examined. For each design, the volume fraction and type of nanoparticles, the shell aspect ratio, and the fin aspect ratio were optimized using the Taguchi method for maximum melting rate. The major finding of the current research can be summarized as follows: • The most important parameter influencing the melting rate is the aspect ratio of the shell (AR c ). The optimum AR c for the ascending design was 1.0 (a cylindrical shell), and for the descending design, it was 1.17.

•
The fins arrangement was the second most important design parameter, influencing the melting rate. The best fin aspect ratios were the cases with a smooth variation of fin sizes. Thus, the middle size values of AR f = 0.4 and 1.4 for ascending and descending designs were led to optimum melting rate, respectively. The melting maps showed that the fins significantly contribute to heat transfer at the initial stages of melting, where the heat transfer is conduction dominant. However, at the final stages of melting, the natural convection circulations play an important role in melting heat transfer. The fins' arrangement is crucial since long fins could suppress the natural convection flows.

•
The type of nanoparticles was the third most important design parameter, influencing the melting rate. The optimum type of nanoparticles depends on the designed geometry. Alumina nanoparticles were the best for an ascending design, while copper nanoparticles were a better choice for a descending design. Indeed, the presence of nanoparticles influences the heat capacity and thermal conductivity of the NePCM. Thus, the type of these nano-additives should be selected carefully.

•
The presence of nanoparticles could improve the melting rate. However, the best volume fraction of nanoparticles for a descending design was only 1.5% of nanoparticles, while the maximum possible volume fraction could be 4.5%. This shows that the increase in the volume fraction of the nanoparticles is not always advantageous. The nanoparticles increase the dynamic viscosity of molten NePCM, and thus, they can weaken the natural convection flows.

Conflicts of Interest:
The authors declare no conflict of interest.