Numerical Simulation of the Impact of the Heat Source Position on Melting of a Nano-Enhanced Phase Change Material

A 2D-symmetric numerical study of a new design of Nano-Enhanced Phase change material (NEPCM)-filled enclosure is presented in this paper. The enclosure is equipped with an inner tube allowing the circulation of the heat transfer fluid (HTF); n-Octadecane is chosen as phase change material (PCM). Comsol-Multiphysics commercial code was used to solve the governing equations. This study has been performed to examine the heat distribution and melting rate under the influence of the inner-tube position and the concentration of the nanoparticles dispersed in the PCM. The inner tube was located at three different vertical positions and the nanoparticle concentration was varied from 0 to 0.06. The results revealed that both heat transfer/melting rates are improved when the inner tube is located at the bottom region of the enclosure and by increasing the concentration of the nanoparticles. The addition of the nanoparticles enhances the heat transfer due to the considerable increase in conductivity. On the other hand, by placing the tube in the bottom area of the enclosure, the liquid PCM gets a wider space, allowing the intensification of the natural convection.


Introduction
The recent industrial revolution required a huge amount of energy, where petrol and gas have been the available energy sources for many years. However, this kind of energy source will surely vanish after few years, regardless of the environmental considerations. Humanity is obligated to find new energy sources, such as solar energy [1], that are renewable and environmentally friendly, but there is an issue with this kind of energy, which is the ability to store heat for a long period to be used when needed. In recent years, phase change materials (PCMs) have become the best choice to store thermal energy because of their high thermal capacity; however, the problem with these materials is their poor thermal conductivity.
The weakness of PCM thermal conductivity is the main obstacle during the process of charging/discharging energy in thermal energy storage systems (TES) [2]. Investigators have suggested and applied several techniques such as adding high conductive fins [3,4] and nanoparticles into the PCM [5,6], where nanoparticles can be used to enhance heat distribution/dissipation due to their high thermal conductivities. Nanofluids can be used in heat exchanger systems [7], where the heat is dissipated via the liquid. The addition of nanoparticles is considered to be an innovative and non-destructive solution to enhance the heat transfer rate. However, a high volume fraction of nanoparticles leads to higher viscosities that cause a significant pressure drop, decreased pumping power, and viscous irreversibility. The same thing can happen with PCM-thermal charge/discharge applications when free convective motions (buoyancy force) occur within the molten region. results showed that the melting rate increases by increasing the dispersion of nanoparticles, where the heat distribution into the PCM has been enhanced, making the melting rate accelerate.
The nanoparticle effect has also been analyzed during the solidification process. Some of the studies related to this subject are described in the following paragraphs.
Sheikholeslami et al. [17] numerically investigated the heat transfer during the solidification of NEPCM using different CuO nanoparticle shapes (Spherical, Platelet, Cylinder, and Brick), which were dispersed in water (used as PCM), while changing the length of the fins. The authors concluded that the solidification time may diminish by using platelet nanoparticles. Furthermore, increasing the length of the fins improved the conduction mode; consequently, the solidification rate was improved. Du et al. [18] carried out a 3D numerical study to investigate the heat transfer within paraffin (PCM) mixed with copper nanoparticles inside a cylindrical tank equipped with a spiral coil heat exchanger. The results showed that 19.6% of the solidification time was saved by using nanoparticles compared to the case of pure PCM. In addition, nanoparticle dispersion reduced the temperature non-uniformity in the tank. This indicates that the convective heat transfer has decreased because it is well known that convection depends mainly on the movement of the fluid, and the intensity of this movement may diminish with a high dynamic viscosity. It can be concluded that the amount of nanocomposite increases the viscosity, and therefore the PCM fluid cannot move freely.
Hosseinzadeh et al. [19] numerically compared two heat transfer enhancement techniques: adding hybrid-nanoparticles and adding optimized fins into a pure PCM during the discharging process. The comparison proved that the solidification time (discharging time) was less after adding optimized fins compared to dispersing nanoparticles into the PCM. The fins increased the surface of heat exchange and reduced the convective flow of the movement of the liquid PCM, causing enhancement of the performance of the solidification process. Actually, the convective flow is an undesirable factor during the solidification process [8].
Some researchers have studied the coupled effects of using nanoparticles to enhance heat transfer and applying a magnetic field to control the flow. Sheikholeslami and Mahian [20] numerically investigated the influence of magnetic force during the solidification process of an NEPCM. The results showed that, by increasing the Hartmann number from 0 to 10, the solidification rate was improved by 23.5% on average. On the other hand, by adding 4% nanoparticles to the PCM, the solidification rate was augmented by 14%. The results indicated that the magnetic force reduces the velocity of the liquid movement, which is caused by the free convection phenomena, to accelerate the solidification during the heat discharging process. Sheikholeslami [21] conducted a numerical study on the influence of magnetic force on the solidification/discharging process of CuO nanoparticles dispersed in water inside a cylindrical tank with an inner tube. The investigation proved that the solidification rate has a direct relation to the magnetic field and the volume fraction of the nanofluid. The magnetic force can act as a beneficial factor by reducing the molten motion and making the solidification process faster. Sheikholeslami [22] numerically investigated the effects of nanoparticle diameter and volume fraction on the PCM solidification process. The study proved that the solidification rate is enhanced by adding nanoparticles, and the average temperature increases by increasing the concentration and the diameter of the nanoparticles in the PCM. Sheikholeslami [23] performed a numerical study on the behavior of an NEPCM during the discharging process (solidification) inside a cylindrical tank. The work presented the impacts of the inner tube and the volume fractions of CuO nanoparticles. It has been noticed that the shape of the inner tube has a considerable effect on the rate of solidification; in addition, the peak rate of solidification was obtained with dp = 40 nm (nanoparticle diameter). Therefore, the shape of the heat source has an important effect that must be investigated, and it is recommended that the position/shape of the heat source (HTF tube) is investigated, as well as the effect of the shape/dimension/concentration of the dispersed nanoparticles. Several researchers, such as Milad Sadeghzadeh et al. [24,25], have proposed techniques to control the dynamic viscosity of nanofluids by using an innovative method based on collecting input data such as the nanofluid concentration, the flow rate, and the inlet temperature. They found that the best prediction of thermal performance was obtained by an MLP (Multilayer Perception Artificial Neural Network).
From the above-described literature, it should be mentioned that, in the solidification and melting of Nano-Enhanced Phase Change Materials (NEPCMs), the addition of nanoparticles has a positive effect on the performances of these processes. In addition, several factors are thought to optimize these processes, such as the enclosure geometry; the thermal condition (the position, intensity, etc.); and the nanoparticles diameter, shape, and thermophysical properties.
The current study was carried out to investigate a new design of NEPCM-filled enclosure (see Figure 1). The influence of the inner-tube position and the benefit of adding Al 2 O 3 nanoparticles as dispersion to the PCM were investigated. The inner tube had three different vertical positions and the nanoparticle volume fraction was ranged as 0 ≤ φ ≤ 0.06. The results are presented and discussed in terms of the liquid-solid fraction, thermal field, flow structure, and Nusselt number. The novelty of this work is based on coupling the influence of varying the position of the heat source (HTF tube) with the addition of the nanoparticles to the PCM capsule.
a cylindrical tank. The work presented the impacts of the inner tube and the volume fractions of CuO nanoparticles. It has been noticed that the shape of the inner tube has a considerable effect on the rate of solidification; in addition, the peak rate of solidification was obtained with dp = 40 nm (nanoparticle diameter). Therefore, the shape of the heat source has an important effect that must be investigated, and it is recommended that the position/shape of the heat source (HTF tube) is investigated, as well as the effect of the shape/dimension/concentration of the dispersed nanoparticles.
Several researchers, such as Milad Sadeghzadeh et al. [24,25], have proposed techniques to control the dynamic viscosity of nanofluids by using an innovative method based on collecting input data such as the nanofluid concentration, the flow rate, and the inlet temperature. They found that the best prediction of thermal performance was obtained by an MLP (Multilayer Perception Artificial Neural Network).
From the above-described literature, it should be mentioned that, in the solidification and melting of Nano-Enhanced Phase Change Materials (NEPCMs), the addition of nanoparticles has a positive effect on the performances of these processes. In addition, several factors are thought to optimize these processes, such as the enclosure geometry; the thermal condition (the position, intensity, etc.); and the nanoparticles diameter, shape, and thermophysical properties.
The current study was carried out to investigate a new design of NEPCM-filled enclosure (see Figure 1). The influence of the inner-tube position and the benefit of adding Al2O3 nanoparticles as dispersion to the PCM were investigated. The inner tube had three different vertical positions and the nanoparticle volume fraction was ranged as 0 0.06    . The results are presented and discussed in terms of the liquid-solid fraction, thermal field, flow structure, and Nusselt number. The novelty of this work is based on coupling the influence of varying the position of the heat source (HTF tube) with the addition of the nanoparticles to the PCM capsule.

Studied Configuration
The studied configurations and the fixed boundary conditions are presented in Figure 1. A cylindrical tube is located inside a rectangular cavity (30/25 mm) filled with NEPCM, and the HTF passes through the tube to ensure the melting of the NEPCM. The present investigation concerns the influence of adding Al 2 O 3 nanoparticles at various concentrations on the melting process. Three positions of the HTF tube are considered, as presented in Figure 1c,d. The configuration is considered as 2D-axisymmetric to reduce the calculation time. The proprieties of the PCM (n-Octadecane), the nanoparticles (Al 2 O 3 ), and the inner-tube wall are presented in Table 1. Table 1. Thermophysical proprieties of n-Octadecane (PCM),inner-tube wall, and Al 2 O 3 nanoparticles [15].

Mathematical Model
The heat transfer during the phase change phenomenon is investigated in this study. The mathematical formulation developed to solve this kind of problem is based on Voller's model [26]. The enthalpy method based on the Darcy law is combined with the mentioned model by considering the entire domain as a porous medium. The porosity is ranged as: 0 ≤ ε ≤ 1, where, in the liquidus zone, ε = 1, in the solidus zone, ε = 0, and in the intermediate "mushy" zone, the porosity (liquid fraction) is between 0 < ε < 1.
By assuming the liquid NEPCM as an incompressible and Newtonian fluid, the governing equations of the present model, including the continuity, modified momentum, and energy in 2D coordinates, are as follow: Continuity equation where ρ is density, µ is dynamic viscosity, and P is pressure. S a is the source term of the porosity function [27]. Based on the Carman-Kozeny model [28] and Darcy's law, S a and ∇(P) are expressed as: C is a mushy zone morphology constant and is fixed at 10 6 kg/m 3 ·s [29]. η is the liquid fraction expression of the liquid/solid zone and is used to determine the computed zone, where, in the liquid zone, η = 1, in the solid zone, η = 0, and in the mushy zone, 0 < η < 1. Where where T l and T s are the liquid and solid temperature of the PCM, respectively. By using: (T s = T m − ∆T) and (T l = T m + ∆T), where T m is the melting point of the PCM and ∆T is the melting temperature range (∆T = 2K), the liquid fraction expression becomes: where Equation (5) is used as a conditional equation in Comsol-Multiphysics. Based on the Darcy law, the momentum equation for a fluid crossing a porous medium is expressed as: where κ is the permeability, which depends on the liquid fraction, η, and S b is the Boussinesq approximation of thebuoyancy force: where ρ l is the density of the liquid NEPCM, β is the thermal expansion coefficient, and → g is the gravity acceleration vector.
Energy equation where H is the specific enthalpy: and h is the sensible enthalpy defined as: where k, h re f , T re f , c p , and ∆H are, the thermal conductivity, reference sensible enthalpy, reference temperature, specific heat capacity, and melting latent heat, respectively. The effective properties of the NEPCM dispersion are evaluated as follow [15]: Thermal conductivities of the Liquid NEPCM (k nl ) and the Solid NEPCM (k ns ): where k np , k l and k s are, respectively, the thermal conductivity of the nanomaterial, the liquid PCM, and the solid PCM, and ρ np is the density of the nanomaterial. With: By considering that the properties of the NEPCM within the liquid/solid zone are constant and temperature-independent, the function (f ) is evaluated at T = 320 K and T 0 = 273 K [15], where: β 0 = 8.4407(100φ) −1.07304 , K b = 1.381 × 10 −23 is the Boltzmann constant, and d p is the nanoparticles diameter.
Densities of the Liquid ENPCM (ρ nl ) and the Solid ENPCM (ρ ns ): where ρ l and ρ s are, respectively, the density of the nanomaterial: the liquid PCM and the solid PCM. Specific heats of the Liquid ENPCM c pnl and the Solid ENPCM (c pns ): where c pnp , c pl , and c ps are, respectively, the specific heat of the nanomaterial, the liquid PCM, and the solid PCM. Dynamic Viscosity of the Liquid ENPCM (µ nl ): where µ l is the dynamic viscosity of the liquid PCM. Latent heat of fusion of the ENPCM (∆H nl ): where ∆H is the latent heat of fusion of the pure PCM. Thermal expansion coefficient of the Liquid ENPCM (β nl ): where β np is the thermal expansion coefficient of nanoparticle. The density variation of the NEPCM: The specific heat of the NEPCM: where D is the smoothed Delta Dirac function: D = e (−(T−T m ) 2 /∆T 2 ) / ∆T √ π [8,29]. The thermal conductivity of the NEPCM: The computation grids for the considered cases are displayed in Figure 2. It is observed that the grid of the simulation area is well-refined in order to define the mushy zone transition at each time step [30]. perature, and average velocity at t = 3000 s are chosen as the sensitive variables. For all the considered geometrical configurations (inner tube at the top, central, and bottom positions) the chosen meshes are: 26,770, 27,107, and 26,794 elements, respectively. These amounts of elements were chosen to ensure results accuracy and computational time economy.

Mesh Independency Study
The computation grids for the considered cases are displayed in Figure 2. It is observed that the grid of the simulation area is well-refined in order to define the mushy zone transition at each time step [30].

Solution Procedure
The fixed-grid approach has proved its superiority in studying natural convection effects during phase change problems. Actually, it does not require any explicit treatment of the moving boundary; furthermore, it gives a more realistic and consistent solution to the studied problems [31][32][33].
The resolution procedure is performed using the constant Newtonian iteration technique (Six iterations for each time step, and the damping factor is set at 0.9; the time step is selected automatically to reach the defined error tolerance, and "0.001" is used as a tolerance value). The Backward Euler method is used as a time marching method (the backward order is between 1 and 2).

Solution Procedure
The fixed-grid approach has proved its superiority in studying natural convection effects during phase change problems. Actually, it does not require any explicit treatment of the moving boundary; furthermore, it gives a more realistic and consistent solution to the studied problems [31][32][33].
The resolution procedure is performed using the constant Newtonian iteration technique (Six iterations for each time step, and the damping factor is set at 0.9; the time step is selected automatically to reach the defined error tolerance, and "0.001" is used as a tolerance value). The Backward Euler method is used as a time marching method (the backward order is between 1 and 2).

Model Validation
A validation test was managed by comparing the present model and previously published results. Figure 3 presents the comparison of the liquid fraction obtained by the present model and the liquid fraction of the work of Arasu and Mujumdar [12]; the comparison is showing an acceptable agreement with the findings of reference [12]. A validation test was managed by comparing the present model and previously published results. Figure 3 presents the comparison of the liquid fraction obtained by the present model and the liquid fraction of the work of Arasu and Mujumdar [12]; the comparison is showing an acceptable agreement with the findings of reference [12].

Liquid Fraction, Thermal Field, and the Streamlines at Deferent Time Step
Plots of the liquid fraction, thermal field, and streamlines at different time steps (t = 0, 1, 10, 20, 30, 40, and 50 min) with different nanoparticles concentrations (φ = 0.00, 0.03, and 0.06) and for the three considered geometrical configurations (top, center, and bottom positions of the tube) are shown in Figure 4. The results show that the liquid fraction of all the nanoparticle volume fractions and the different positions of the inner tube behave generally in the same way. The liquid fraction begins to grow from the inner tube towards the adiabatic walls.
In addition, the results in Figure 4 show that the thermal distribution follows a behavior similar to that of the liquid fraction, i.e., the temperature increases from the inner tube to other parts of the container. The thermal field corresponds well to the growth of the liquid-solid fraction and to the position of the mushy zone. The thermal field and the liquid fraction take a regular shape (parallel to the shape of the inner tube) at the early stages of the melting process and change to an irregular shape at the later stages of the process. The thermal field and the mushy zone show up in a wavy form in advanced stages of the process due to the convective effects. Actually, the heat transfer is dominated by the conduction mode in earlier phases.
Streamlines in Figure 4 show the flow structure inside the liquid PCM. It is noticed that the molten PCM moves in rotating cells due to the thermal gradation effect that causes the hot parts of the liquid PCM to rise and the colder parts to fall due to the buoyancy force.

Liquid Fraction, Thermal Field, and the Streamlines at Deferent Time Step
Plots of the liquid fraction, thermal field, and streamlines at different time steps (t = 0, 1, 10, 20, 30, 40, and 50 min) with different nanoparticles concentrations (φ = 0.00, 0.03, and 0.06) and for the three considered geometrical configurations (top, center, and bottom positions of the tube) are shown in Figure 4. The results show that the liquid fraction of all the nanoparticle volume fractions and the different positions of the inner tube behave generally in the same way. The liquid fraction begins to grow from the inner tube towards the adiabatic walls.
In addition, the results in Figure 4 show that the thermal distribution follows a behavior similar to that of the liquid fraction, i.e., the temperature increases from the inner tube to other parts of the container. The thermal field corresponds well to the growth of the liquid-solid fraction and to the position of the mushy zone. The thermal field and the liquid fraction take a regular shape (parallel to the shape of the inner tube) at the early stages of the melting process and change to an irregular shape at the later stages of the process. The thermal field and the mushy zone show up in a wavy form in advanced stages of the process due to the convective effects. Actually, the heat transfer is dominated by the conduction mode in earlier phases.
Streamlines in Figure 4 show the flow structure inside the liquid PCM. It is noticed that the molten PCM moves in rotating cells due to the thermal gradation effect that causes the hot parts of the liquid PCM to rise and the colder parts to fall due to the buoyancy force.
The progress of the melting process is enhanced compared to the pure PCM case (φ = 0.00) by the addition of nanoparticles. Indeed, the melting process accelerates with the increase of the nanoparticle's concentration for all the considered positions of the inner tube.
At t = 1 min, the shapes of the thermal field and the liquid zone are formed regularly according to the shape of the inner tube, for all the geometrical cases and all the nanoparticles volume fractions. In addition, the streamlines follow trajectories that are parallel to the inner tube. When the time is up, the streamlines take a circular shape, and the flow intensity of the molten PCM becomes more significant; thus, the shape of the molten zone in the advanced stage has a wavy shape.
At 40 min and 50 min, for the top-position tube, the melt reaches an area under the inner tube characterized by a regular shape due to the dominance of the conductive heat transfer heat mode. By increasing the nanoparticle concentration, the heat transfer is enhanced, especially when it is controlled by the conduction mode. In other words, during  The progress of the melting process is enhanced compared to the pure PCM case ( φ = 0.00) by the addition of nanoparticles. Indeed, the melting process accelerates with the increase of the nanoparticle's concentration for all the considered positions of the inner tube. At t = 1 min, the shapes of the thermal field and the liquid zone are formed regularly according to the shape of the inner tube, for all the geometrical cases and all the nanoparticles volume fractions. In addition, the streamlines follow trajectories that are parallel to the inner tube. When the time is up, the streamlines take a circular shape, and the flow intensity of the molten PCM becomes more significant; thus, the shape of the molten zone in the advanced stage has a wavy shape.
At 40 min and 50 min, for the top-position tube, the melt reaches an area under the inner tube characterized by a regular shape due to the dominance of the conductive heat  The evolution of Nusselt takes the same behavior for all the considered configurations, where the Nusselt Nu increases sharply after few seconds from the running time to reach a peak value, and then becomes quasi-constant until the end of the process. Nu reaches the highest values with the cases of nanoparticles concentration equal to 0.00 (pure PCM) with all tube position cases and it takes the lowest values with the case of Nano concentration is φ = 0.06. In addition, the Nusselt values are increased by moving the inner tube downwards. As well known the Nusselt number reflects the intensity of the convection compression to conduction, from this figure it is clear that the convective heat transfer dominates the conduction since the early stage of the melting process, besides the convection effect increases as the inner tube is moved to the bottom section of the tank.

Average Nusselt Number
the highest values with the cases of nanoparticles concentration equal to 0.00 (pure PCM with all tube position cases and it takes the lowest values with the case of Nano concen tration is φ = 0.06. In addition, the Nusselt values are increased by moving the inner tub downwards. As well known the Nusselt number reflects the intensity of the convection compression to conduction, from this figure it is clear that the convective heat transfer dom inates the conduction since the early stage of the melting process, besides the convection effect increases as the inner tube is moved to the bottom section of the tank. .  Figure 6 depicts the temporal variation of the molten fraction, including the effect of the nanoparticle concentration and the inner-tube position. It is clear that the molten fraction increases with time in all cases. The blue, red, and black curves refer to the loca tion of the inner tube, i.e., the bottom, central, and top positions, respectively. The melting time decreases as the inner tube is moved to the bottom of the cavity, and the melting tim decreases by increasing the nanoparticle concentration due to the enhancement of the hea  Figure 6 depicts the temporal variation of the molten fraction, including the effects of the nanoparticle concentration and the inner-tube position. It is clear that the molten fraction increases with time in all cases. The blue, red, and black curves refer to the location of the inner tube, i.e., the bottom, central, and top positions, respectively. The melting time decreases as the inner tube is moved to the bottom of the cavity, and the melting time decreases by increasing the nanoparticle concentration due to the enhancement of the heat transfer. At the beginning (t < 500 s) of the melting process, the effect of adding nanoparticles is not very effective on the melt fraction. When the time is up, the liquid fraction is considerably increased, especially for higher nanoparticle concentrations. This is an interesting result, because the addition of low-cost nanoparticles, such as Al 2 O 3 , will lead to the amelioration of the charging process performances. In fact, for all the considered inner-tube positions, the melt fraction is approximatively 10% higher for φ = 0.06 (NEPCM) compared to the pure PCM.   Figure 7 shows the temporal variation of the thermal energy stored by the NEPCM. The stored energy increases with time in all cases, as presented in Figure 7. It can also be noted that the energy storage rate follows the same pattern of evolution of the melting rate. The rate of charge of stored energy increases with the displacement of the inner tube  Figure 7 shows the temporal variation of the thermal energy stored by the NEPCM. The stored energy increases with time in all cases, as presented in Figure 7. It can also be noted that the energy storage rate follows the same pattern of evolution of the melting rate. The rate of charge of stored energy increases with the displacement of the inner tube towards the bottom of the cavity, and with the increase of the concentration of nanoparticles. This improvement is due to the increased heat distribution in the NEPCM.  Figure 7 shows the temporal variation of the thermal energy stored by the NEPCM. The stored energy increases with time in all cases, as presented in Figure 7. It can also be noted that the energy storage rate follows the same pattern of evolution of the melting rate. The rate of charge of stored energy increases with the displacement of the inner tube towards the bottom of the cavity, and with the increase of the concentration of nanoparticles. This improvement is due to the increased heat distribution in the NEPCM.

Concluding Remarks
This numerical study aimed to investigate the effect of adding Al2O3 nanoparticles to a PCM-filled cavity equipped with an inner HTF tube. The influence of the position of the internal tube on thermal behavior has also been studied. The main findings of the study can be summarized as follows: • The melting rate increases as the position of the inner tube changes towards the bottom section of the cavity; this is because the upper region is larger and free convection is the dominating heat transfer mode in this region. In addition, conduction is the

Concluding Remarks
This numerical study aimed to investigate the effect of adding Al 2 O 3 nanoparticles to a PCM-filled cavity equipped with an inner HTF tube. The influence of the position of the internal tube on thermal behavior has also been studied. The main findings of the study can be summarized as follows:

•
The melting rate increases as the position of the inner tube changes towards the bottom section of the cavity; this is because the upper region is larger and free convection is the dominating heat transfer mode in this region. In addition, conduction is the dominating mode in the lower region, especially at the beginning of the melting process. Therefore, the amount of stored energy increases as the influence of convection increases with the displacement of the inner tube downwards.

•
The melting and charging rates are improved by increasing the volume fraction of the nanoparticles and, due to the enhancement, the thermal conductivity that ameliorates the heat transfer rate.

•
The temporal evolutions of the melting process and the amount of stored energy are quasi-similar at the start of the process because the heat input was managed by conduction. However, at later stages of the process, the heat transfer becomes dominated by convection.