Next Article in Journal
Analysis of Transient Interactions between a PWR Nuclear Power Plant and a Faulted Electricity Grid
Next Article in Special Issue
Conical Two-Phase Swirl Flow Atomizers—Numerical and Experimental Study
Previous Article in Journal
Smart Organizations as a Source of Competitiveness and Sustainable Development in the Age of Industry 4.0: Integration of Micro and Macro Perspective
Previous Article in Special Issue
Numerical Study on the Cavitation Flow and Its Effect on the Structural Integrity of Multi-Stage Orifice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation of a Fast-Charging Porous Thermal Energy Storage System Saturated with a Nano-Enhanced Phase Change Material

1
Metamaterials for Mechanical, Biomechanical and Multiphysical Applications Research Group, Ton Duc Thang University, Ho Chi Minh City 758307, Vietnam
2
Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City 758307, Vietnam
3
Young Researchers and Elite Club, Yasooj Branch, Islamic Azad University, Yasooj 7591493686, Iran
4
Faculty of Mechanical and Energy Engineering, Shahid Beheshti University, Tehran 1983969411, Iran
5
Department of Mechanical Engineering, Shahid Chamran University of Ahvaz, Ahvaz 61355, Iran
6
Department of Mechanical Engineering, College of Engineering at Wadi Addwaser, Prince Sattam Bin Abdulaziz University, Wadi Addwaser 11991, Saudi Arabia
7
Department of Mechanical Engineering, Faculty of Engineering, University of Khartoum, Khartoum 11111, Sudan
8
Mechanical & Automotive, School of Engineering, RMIT University, Bundoora, VIC 3083, Australia
9
Department of Mechanical and Aeronautical Engineering, Clarkson University, Potsdam, NY 13699, USA
*
Author to whom correspondence should be addressed.
Energies 2021, 14(6), 1575; https://doi.org/10.3390/en14061575
Submission received: 27 January 2021 / Revised: 24 February 2021 / Accepted: 6 March 2021 / Published: 12 March 2021
(This article belongs to the Special Issue Computational Heat Transfer and Fluid Mechanics)

Abstract

:
The melting of a coconut oil–CuO phase change material (PCM) embedded in an engineered nonuniform copper foam was theoretically analyzed to reduce the charging time of a thermal energy storage unit. A nonuniform metal foam could improve the effective thermal conductivity of a porous medium at regions with dominant conduction heat transfer by increasing local porosity. Moreover, the increase in porosity contributes to flow circulation in the natural convection-dominant regimes and adds a positive impact to the heat transfer rate, but it reduces the conduction heat transfer and overall heat transfer. The Taguchi optimization method was used to minimize the charging time of a shell-and-tube thermal energy storage (TES) unit by optimizing the porosity gradient, volume fractions of nanoparticles, average porosity, and porous pore sizes. The results showed that porosity is the most significant factor and lower porosity has a faster charging rate. A nonuniform porosity reduces the charging time of TES. The size of porous pores induces a negligible impact on the charging time. Lastly, the increase in volume fractions of nanoparticles reduces the charging time, but it has a minimal impact on the TES unit’s charging power.

1. Introduction

The design of thermal energy storage (TES) systems has improved significantly in recent years. Phase change materials (PCMs) have been implemented in various practical applications such as buildings materials [1], air heating and cooling in buildings [2], cooling of electronic components [3], recovering low-temperature industrial waste heat [4], and automotive applications [5].
Most of the advances in TES system design have been around improving the response time of storage units during the charging and discharging process and the synthesis of new composite PCM materials. For instance, using expanded graphite additives [6], tree-like fins [7], stepped fins [8], and heat pipes [9] are some of the recent techniques. Moreover, thermal conductive structures such as three-dimensional (3D) porous diamond foams [10] and continuous diamond–carbon nanotube foams [11] have also shown promising performance.
Many researchers attempted using a frame of thermal conductive metals and porous metal foams, saturating them with phase change material. The thermal conductive structure of the metal foam carries out the heat to and from PCM inside the pores. Thus, the composite PCM–metal foam channels the heat from the heat sources and improves the thermal conductivity. Open metal foams allow some degree of natural convection heat transfer in molten areas of the TES unit. The natural convection heat transfer is an important phenomenon that transfers thermal energy by advection in a molten region. Considering open metal foams and natural convection effects, Talebizadehsardari [12] utilized a composite of metal foam and PCMs and designed an air heater for domestic application. The results indicated that the geometrical shape of air passages induces dominant effects on the PCM unit’s thermal behavior. Sardari et al. [13] used an aluminum foam–PCM composite TES unit to absorb a wall-mounted radiator’s heat during off-peak loads. Later, the TES unit releases the absorbed heat when the primary heating system turns off. The authors showed that using the metal foam shortens the charging time by 95% compared to a simple PCM. Interestingly, the increase in metal foam porosity (97%) produces a positive influence and reduces the TES unit’s charging/discharging time.
Zhao et al. [14] compared the advantages of using fins with those of using metal foams in reducing the melting time of a PCM in a shell-and-tube shape TES unit. They found an optimum design for the fins. They then compared the melting time of using the optimum fin with a case made using the same amount of metal foam instead of the fins. They reported that properly designed fins could be as advantageous as metal foam.
Using nanoparticles is another approach to synthesize nano-enhanced phase change materials (NePCMs). Zhang et al. [15] dispersed copper oxide nanoparticles in RT28 PCM to improve the charging heat transfer rate in a wavy-channel TES unit. The presence of nanoparticles improved the heat transfer rate and reduced the charging time. Bondareva et al. [16] examined the melting of an alumina–paraffin NePCM inside a copper radiator. This investigation showed that the nanoparticle presence increases the viscosity of molten PCM, limiting the mobility of liquid PCM.
On the other hand, nano-additives enhance the thermal conductivity of the PCM. Thus, using NePCM was most advantageous at an early melting heat transfer stage, where the conduction regime was dominant. It can be concluded that the advantage of using nanoparticles mainly depends on the internal structure of the TES enclosure and dominant mechanisms of heat transfer.
All of the above studies used a uniform metal foam to improve the thermal conductivity of composite PCM. Mahdi et al. [17] employed a cascade (multiple segments) of PCM–metal foam composites to enhance a TES unit’s discharging time. They also employed 5% nanoparticles to further enhance the thermal conductivity of the system. Using the multiple segments and nanoparticles, they reduced the discharging time by 94%.
The literature studies showed that the natural convection effects and the thermal conductivity of composite PCMs could significantly influence a TES unit’s phase change behavior. A dense porous structure (low porosity) increases the composite thermal conductivity but suppresses convection flows. This is while a high porosity metal foam allows convection flows. Thus, a nonuniform porous metal foam could be advantageous from both points of view. The present study investigates the melting process of a PCM embedded in a nonuniform metal foam. The influence of using various concentrations of nanoparticles, porous densities, pores sizes, and porous gradients on the full charging time of the TES unit are addressed. The Taguchi optimization approach was employed to systematically find an optimum design of nonuniform porosity. In particular, the study addressed the following research questions: (1) Does a nonuniform porosity reduce the melting time for a fixed amount of metal foam? (2) What is the impact of nanoparticles and pore size on the melting time of TES?

2. Governing Equations and Boundary Conditions

A schematic view of the physical and computational domains with dimensions of the latent heat thermal energy storage unit is depicted in Figure 1a,b. The TES unit is a cylinder with a longitudinal inner tube containing heat transfer fluid (HTF). The cylinder is filled with an inhomogeneous porous medium in which porosity and permeability vary along the z-direction. The porosity in the direction of the z-axis could be increasing, constant, or decreasing. Biobased coconut oil containing the CuO nanoparticles fills the pores. Table 1 lists the properties of coconut oil, CuO, copper metal foam, and water.
The modeling was done under the following assumptions: (1) the melted NePCM behaves as an incompressible Newtonian fluid, and the flow in the pores is laminar; (2) nanoparticles do not settle, and they are homogeneously dispersed in the PCM; (3) expansion of the NePCM is neglected; (4) the properties of the HTF in tube are considered to be constant and Newtonian. Here, a linear porosity profile is assumed as
ε ( z ) = a z + b ,
where Equation (1) shows the linear variation of porosity which is proposed in the present study. As seen, the porosity varies linearly in the vertical direction. The values of porosity at the lower and upper walls are
ε ( z = 0 ) = ε 0 = b ,
ε ( z = L ) = ε L .
According to the values of porosity on the lower and upper walls, and the average porosity of the domain, i.e., ε a v g , the b parameter is calculated as follows:
b = ε a v g a L 2 .
The equations for the transient melting process include mass conservation, momentum balance, and energy conservation equations [20,21,22].
Mass conservation:
U = 0 ;
Momentum equations:
1 ε ( ρ l n p c m U ) t + 1 ε ρ l n p c m ( U ) U 1 ε + p I = ( μ l n p c m 1 ε U ) ρ l n p c m β l n p c m g ( θ m e l t i n g θ ) + ( M ( θ ) μ l n p c m σ ) U
where U , p, and θ are the velocity vector, pressure, and temperature, respectively. Here, ε, β, μ, and g are the porosity, volume expansion coefficient, viscosity, and gravity acceleration, respectively. The subscript lnpcm denotes the molten PCM, and θmelting is the phase change temperature. The linear Boussinesq model is employed to apply the buoyancy force, and the Forchheimer term was neglected since the natural convection velocities are very small. The source term M(θ) is a function of temperature (θ) and is zero in liquid PCM, but it rises to large values in the solid PCM. The large values of M(θ) induce significant resistance forces to the fluid motion and force the velocities to zero in a solid PCM. The source term M(θ) is defined as
M ( θ ) = B v c ξ 2 ( θ ) 2 ξ ( θ ) + 1 ξ 3 ( θ ) + 10 3 ξ ( θ ) = { 0 θ < θ m e l t i n g θ w i n d o w / 2 0.5 θ m e l t i n g θ θ w i n d o w θ m e l t i n g 0.5 θ w i n d o w < θ < θ m e l t i n g + 0.5 θ w i n d o w 1 θ > θ m e l t i n g + 0.5 θ w i n d o w
where Bvc is a large value, which intensifies the magnitude of the temperature-dependent term. The porous permeability, σ, as a function of porosity, is calculated as follows [23,24]:
σ = 0.73 × 10 3 ( 1 ε ) 0.224 R m 2 ( R z R m 1 ) 1.11 R z R m 1 = 1.3924 ( 1 ε ) 3 π 1 1 e 25 ( ε 1 ) R m = 0.254 × 10 7 χ 1   ( P P I )
where Rz and Rm are the pore characteristics, as introduced in [23,24], and χ is the pore density in pore size per inch (PPI).
Energy conservation [25]:
[ ξ ( θ ) [ ( ρ c p ) e f f _ l n p c m ( ρ c p ) e f f _ s n p c m ] + ( ρ c p ) e f f _ s n p c m ] θ t + ( ρ c p ) l n p c m U . θ = . ( κ e f f , n e p c m θ ) + ( V F n a 1 ) ρ p p c m l p p c m ε ξ ( θ ) t
where
( ρ c p ) e f f , l n p c m ( s n p c m ) = ( 1 ε ) ( ρ c p ) s m + ε ( ρ c p ) l n e p c m ( s n e p c m ) .
The subscript snepcm denotes the characteristics of the NePCM in solid state, and sm refers to the solid matrix of metal foam. The term ξ ( θ ) denotes the liquid volume fraction and can be changed in each element. The variation of the liquid fraction represents the variation of the stored/released latent heat. The thermal conductivity of the porous medium saturated with the NePCM is considered as the weighted average of the two media and is computed as
κ e f f , n e p c m = ξ ( θ ) κ e f f , l n p c m + ( 1 ξ ( θ ) ) ξ e f f , s n p c m ,
in which,
κ e f f , l n p c m ( s n p c m ) = Ξ 1 Ξ 1 ,
Ξ 1 = [ κ l n p c m ( s n p c m ) + ( ( π π ε 3 ) 0.5 1 ε 3 ) ( κ s m κ l n p c m ( s n p c m ) ) ] × [ κ l n p c m ( s n p c m ) + ( 1 ε 3 ) ( κ m f κ l n p c m ( s n p c m ) ) ]
Ξ 2 = κ l n p c m ( s n p c m ) + [ 4 3 ( 1 ε 3 π ) 0.5 ( 1 ε ) + ( π π ε 3 ) 0.5 ( 1 ε ) ] ( κ m f κ l n p c m ( s n p c m ) ) ,
The governing equations for flow and heat transfer of HTF passing the tube are laminar convection heat transfer.
Mass conservation:
U = 0 ;
Momentum Balance:
ρ H T F U t + ρ H T F ( U ) U + p I = ( μ H T F U ) ;
Energy conservation:
( ρ c p ) H T F θ t + U . ( ( ρ c p ) H T F θ ) = . ( κ H T F θ ) ,
where subscript HTF shows the heat transfer fluid.
The thermophysical properties of the nano-enhanced phase change material can be defined as weighted functions of the pure PCM and nanoparticle properties. These properties are tabulated in Table 2. The initial and boundary conditions are listed in Table 3.
Sensible energy, latent energy, and total energy stored in the unit are
E S ( t ) = V ( ( 1 V F n a ) ρ p p c m ε ( z ) l p p c m L a t e n t   e n e r g y + ( ρ C p ) e f f _ n e p c m ( θ θ i n i t i a l ) S e n s i b l e   e n e r g y T o t a l   e n e r g y   s t o r e d )   d V ,
where V is the total volume of the NePCM domain. The melt volume fraction is
M V F ( t ) = S ln e p c m S ln e p c m + S s n e p c m .
The charging power of the TES unit is defined as
C P = E S | M V F = 1 t | M V F = 1 .
Moreover, the uniformity of the temperature field can be measured using
T U 2 ( t ) = ( V ( θ θ m e a n ) 2 d V ) / ( V d V ) ,
where θ m e a n is the mean temperature of the NePCM domain.

3. Numerical Approach and Grid Dependency

The model equations described above were solved by employing the finite element approach. A mesh independence test was performed to determine a balance between the precision and performance of the numerical simulation. Structural meshes with rectangular elements were employed for both computational domains. Initially, a mesh with 75 × 75 elements was considered. Then, five levels of finer meshes were used to test the solution accuracy and grid optimization. The simulation results for different grids are depicted in Figure 2. Figure 2a,b show, respectively, the mesh size’s influence on the melting volume fraction (MVF) and local temperature at a specified point. It is seen that the temperature field is more sensitive to the mesh size. On the basis of Figure 2b, a mesh size with 125 × 125 elements was selected to satisfy the balance between the solution accuracy and simulation efficiency.
To guarantee the reliability of the implemented model and the employed numerical solution, we studied the liquid fraction field of biobased coconut oil contained in a uniform porous rectangular enclosure by conducting the two-dimensional (2D) numerical simulation. The predicted liquid fraction field was compared with the experimental observations of Al-Jethelah et al. [18]. In this validation test, a net heat flux was imposed on the cubic cavity’s left side wall, and the insulated boundary conditions were set at the other surfaces. The porosity of copper metal foam used in the container was 0.92, and the corresponding permeability was 3.3142 × 10−7 m2. The liquid fraction field obtained from the present model is in good agreement with the experimental field observation of Al-Jethelah et al. [18], as shown in Figure 3.
To verify the computational model for application to an inhomogeneous porous domain, a square cavity occupied with an inhomogeneous metal foam utilized in [26] was studied. Here, the insulation boundary conditions were assigned to the cavity’s upper and lower surfaces. However, the left and right vertical side boundaries were kept at high and low temperatures. The isotherms and streamlines of the present study are in excellent agreement with those of Xiong et al. [26], as illustrated in Figure 4. In this analysis, porosity was a linear function of the y-coordinate such that εavg=0.7. In addition, the nondimensional parameters were as follows: Ra = 106 and Da = 10−1.

4. Results and Discussion

Here, the Taguchi optimization method was adopted to systematically investigate the impact of design parameters on the melting rate of PCM in the TES unit. The volume fraction of nanoparticles (VFna), average porosity of porous medium ( ε a v g ), intensity of porosity gradient (a), and pore size per inch (PPI) were selected as design parameters. These parameters are known as control factors in the context of the Taguchi method. The aim was to minimize the full melting time (charging time). Thus, the charging time was selected as the target parameter. Since the time should be minimized, the target function was selected as “the lower, the better” for the Taguchi method. Following the Taguchi method, the control factors were divided into several levels. Here, five levels were selected for each factor. The details of the levels and factors are summarized in Table 4. The porosity gradient was selected in the range of −0.6 to +0.6, which covered a negative and positive gradient distribution. The zero gradient case represents uniform porosity.
The total combination of four factors and five levels resulted in 54 combination variables. A full melting process should be computed for each combination, which is computationally impractical. The Taguchi method uses an orthogonal table to probe the possible solution space efficiently. Here, we adopted the L25 design. The L25 design selects only 25 unique combinations of possible designs out of 54 possibilities. These selected combinations are summarized in Table 5. The simulations were carried out for all 25 cases listed in Table 5, and the corresponding values for the time of full melting ( t | M V F = 1 ), stored energy (ES), and TES unit power (CP) at full charge are reported in Table 5. The Taguchi method was then used to compute the signal to noise (S/N) ratio values. The S/N ratio indicates the robustness of a factor level to possible noises. Thus, a factor with the highest S/N value is a promising candidate to produce an optimum design (lowest charging time). Using the values of t | M V F = 1 from Table 5, the S/N of the Taguchi method was computed, and the results are summarized in Table 5.
Following the standard Taguchi method, Table 6 shows the S/N and rank values of each level and control factor. Figure 5 shows a graphical representation of S/N values in Table 6. The highest value of S/N for each control parameter shows the promising level for the minimum charging time. Table 6 also shows the rank of each factor with respect to the variation of melting time. As seen, the average porosity ( ε a v g ) was the most effective parameter influencing the melting time, followed by the volume fraction of nanoparticles. The porosity gradient (a) and pore size (PPI) were ranked third and fourth. Here, δ indicates the maximum difference between the computed S/N values of each factor. A higher δ value denotes a more influential factor.
Figure 5 shows that Levels 4, 1, 4, and 4 corresponded to the highest S/N ratios for nanoparticle volume fraction (VFna), average porosity ( ε a v g ), porosity gradient (a), and pore size (PPI), respectively. The approach for computing the S/N ratio was “smaller is better”; thus, a case with a lower melting time produced a larger S/N ratio value. Using Table 4, these levels can be read as VFna = 0.04, ε a v g = 0.8, a = 0.3, and PPI = 25. The Taguchi method estimated a melting time of 8784 s for this design parameter combination. A simulation was performed for this specific combination to confirm the reduction in the TES unit’s charging time. The simulation outcomes showed a full melting (charging) time of 9097 s, which is slightly higher than the estimated time of 8784 s estimated using the Taguchi method. However, a comparison between the computed melting time and the 25 cases of Table 5 confirms that this melting time was the smallest. Thus, the proposed design according to the Taguchi method was selected as the optimum design. The details of the optimum design are shown in Table 7.
Table 5 also shows that case 5 resulted in the highest melting time (15,900 s). Considering the investigated range of control parameters, the optimum melting time of 9097s was 74% better than the highest melting time. It should be noted that case 5 contained no nanoparticles and had the maximum average porosity of 0.9. Moreover, case 21 with a full melting time of 9150 s provided a melting time that was close to the optimum case. The only difference between these two cases was the porosity gradient.
The results of Table 5 were also used to develop a linear relationship for design parameters and the charging time of the TES. That is,
Time + 8892 243.0   V F n a + 1047.0   ε a v g 9.0   a + 45.0   P P I .
Here, 16 more cases were adopted to further explore the impact of the variation of the design parameters on the optimum design. The details of selected cases are presented in Table 8. Figure 6a shows the time evaluations of MVF during the charging process for various values of nanoparticle volume fraction. It is seen that the increase in the concentration of nanoparticles slightly increased the MVF. The differences were only visible in the final stages of charging.
Figure 6b depicts the predicted temperature of point A in the enclosure. The temperature in the solid region increased sharply until it reached the fusion temperature. Then, the temperature slightly increased until the entire PCM changed to liquid at this point. The temperature then increased again. The increase in nanoparticle concentration slightly increased the temperature of point A in the liquid region due to the improved thermal conductivity of NePCM and better heat transfer between the hot liquid in the tube and the PCM inside the enclosure. The melting interfaces during the melting process were very close; thus, they were not plotted.
Figure 7 shows the phase change interfaces ( = 0.5 ) for various average porosity ( ε a v g ) values and different melting times. It is seen that the melting interface advanced in the enclosure toward the right wall with time.
Interestingly, the melting started at the bottom since the HTF fluid at the entry was hot, and the convection heat transfer coefficients at the undeveloped region of the HTF tube was also larger. Thus, the heat transfer between HTF liquid and NePCM at the bottom was more intense. The advancement of the melting interface accelerated as the porosity decreased (mass of metal foam increased). This was due to the increase in thermal conductivity of metal foam, which enhanced the heat transfer rate. It should be noted that the inlet temperature of HTF fluid was constant; hence, a higher composite thermal conductivity at the PCM side increased the heat transfer rate. Thus, the amount of stored heat in an enclosure increased for a higher composite thermal conductivity. Therefore, as seen in Figure 7 at t = 7800 s, there was a clear difference between the melting interfaces. Figure 8a illustrates the time history of MVF for various average porosities. The average porosity induced notable changes in the values of MVF. In agreement with the melting interfaces of Figure 7, this figure shows that the increase in average porosity raised the MVF.
Figure 8b depicts the temperature uniformity inside the PCM domain. The enclosure temperature distribution experienced a slightly higher temperature gradient during the melting process when average porosity was maximum. The rise in temperature was due to the lower composite thermal conductivity and large temperature gradients. At the early stages of heat transfer, the NePCM was in a solid state, and there was no practical molten region or natural convection circulations. As the MVF increased and the molten region grew, natural convection flows occurred. During the phase change, the temperature uniformity stayed almost constant. When the molten region grew significantly, the natural convection circulation started, and the temperature gradients also increased.
Figure 9 is plotted to show the impact of pore sizes on MVF and temperature evaluations of point A. It is seen that the variation of pore size induced a negligible impact on the MVF. There was a slight temperature variation for point A around 6000 s. Half of the enclosure was in a molten state, and there was a dominant natural convection heat transfer flow. In a microscopic view, the overall heat, which diffused to the composite PCM (foam and PCM), first channeled through the pore walls, and then the PCM inside the pores could absorb it. A decrease in the pore size could provide a larger contact surface and a better melting rate. However, since the amount of foam was low (80% or higher porosity), the heat transfer rate was limited by the foam’s capacity to conduct heat into the solid PCM region. Since the amount of foam was limited (constant porosity), the condition heat transfer through the foam was also limited. Thus, the change in pore size could impact the heat transfer, but it did not play a dominant role. Moreover, a smaller pore size led to a smaller permeability, which reduced the natural convection effects.
The variation of pore sizes could affect the permeability of the medium, according to Equation (4); thus, the impact of pore size on the flow and heat transfer can be boosted in convective dominant heat transfer regimes. Since the influence of the pore size (PPI) on the MVF was minimal, the melting interfaces were not plotted for the sake of brevity.
Figure 10 display the dependency of the porosity gradient on the evaluation of melting interfaces at various time steps. Interestingly, the increase in a shifted the melting interface toward the left wall and resulted in an increased melting rate. According to Equation (1), the increase in a meant more metal foam at the bottom and less metal foam at the top. Thus, the PCM’s thermal conductivity at the bottom of the enclosure was higher than that at the top. At the bottom and in the initial stages of phase change, the dominant mechanism of heat transfer was conduction. Thus, the increase in thermal conductivity enhanced the heat transfer. As the phase change continued, the molten region grew, and convection heat transfer occurred in the enclosure’s top regions. In natural convection flows, the increase in porosity of metal foams increased the medium’s permeability. Consequently, the liquid PCM could circulate more easily in the top area, contributing to convection heat transfer. The relative difference in full charging of TES for a = 0.3 (optimum case) and a = 0.6 (case 12, Table 8) was only 0.5%, and Figure 10 shows the close competition of these two cases at 7800 s. Case a = 0.3, which was the adopted optimum design, led to a better melting at the final stage of thermal energy storage (above 90% MVF) compared to case 12 in Table 8.
Figure 11 illustrates the time history of MVF and temperature uniformity for various values of porosity gradients. Figure 11a shows that an increase in a improved the MVF during the middle stages of phase change. This finding is in agreement with the advancement of the melting interfaces, as observed in Figure 10. However, in the final stages of charging, the case a = 0.3 provided better performance than the case of a = 0.6. The temperature uniformity for the initial stages of the melting process was almost independent of the a parameter. However, in the middle stages, where the molten region grew, the increase in a parameter promoted the temperature nonuniformity. In the middle stages, the increase in a promoted temperature nonuniformities. A high value of a led to a stronger heat transfer flow and natural convection flows.
Nine points in the NePCM–metal foam region were selected to investigate the temperature variations during the melting process. These points were located at a1(2Ri, 0.1L/3), a2(2Ri, L/2), a3(2Ri, 2.9L/3), b1(5Ri, 0.1L/3), b2(5Ri, L/2), b3(5Ri, 2.9L/3), c1(9Ri, 0.1L/3), c2(9Ri, L/2), and c3(9Ri, 2.9L/3). Figure 12 shows the variation of temperatures at these nine points in the PCM domain during the melting process. All points started with an initial temperature of 293 K followed by a linear increase in temperature. The linear segment showed pure conduction in a super-cold solid PCM before it reached a melting temperature. Then, there was a flattened segment, in which the temperature slightly increased. This was the phase change stage, where the PCM absorbed the energy in the form of latent heat. Then, there was a semi-linear increase in the temperature, which corresponded to the heat transfer in the molten region with no phase change. The time variations of temperatures for all cases showed that the impact of the porosity gradient (a) on temperature profiles was minimal in the solid state and during melting phase change. However, after the phase change, the molten NePCM circulated in the enclosure, enhancing the effect of the porosity gradient.
Here, a1, b1, and c1 were the points placed at the bottom of the enclosure. Point a1 was next to the tube wall, point b1 was in the middle, and c1 was next to the insulated shell wall. Since a1 was the closest point next to the tube, it experienced a sharp temperature rise at the beginning of the charging process. In contrast, c1 showed a smooth temperature variation since it was far from the tube. Point c3 was placed near the top-right corner of the enclosure, which was the last place to be melted. Thus, during the melting process, this point remained at the fusion temperature.
Figure 13 compares the full charging power of TES units for various values of design factors. Figure 13a shows that the variation of nanoparticle concentration slightly changed the unit power. This is because the presence of nanoparticles improved the heat transfer rate; they also reduced the heat capacity of the storage unit. Figure 13b depicts an increasing trend of power (CP) as the porosity gradient (a) increased. Both cases of a = 0.3 and a = 0.6 showed almost similar charging powers. The increase in average porosity (Figure 13c) reduced the charging power. As the average porosity increased, the heat capacity of the TES also increased. However, a reduction in the composite PCM’s effective thermal conductivity weakened the heat transfer rate and consequently increased the charging time. The variation of pore sizes had a negligible impact on the TES power since PPI variation did not induce notable changes in either the unit’s heat capacity or the melting time.

5. Conclusions

The charging process of a TES unit filled with a nonuniform metal foam was studied. The metal foam was saturated with coconut oil–CuO NEPCM, and the porosity was linearly changed from bottom to top. Hot water flowed in the heat transfer tube. Finite element method FEM was applied to solve the governing partial differential equations for momentum and phase change energy. The Taguchi optimization method was employed to find an optimum design for the TES unit. The results led to the following conclusions:
  • The average porosity was the dominant factor influencing the charging time and power. Thus, in the design of TES units improved by metal foams, the first step should be selecting the average porosity of the metal foam. The uniformity of the porous medium contributed to heat transfer, and a positive value of porosity gradient (a > 0) could reduce the charging time of the TES unit. A positive porosity gradient corresponded to low porosity at the bottom and high porosity at the top of the PCM enclosure.
  • The volume fraction of nanoparticles was the second most important parameter affecting the TES unit’s phase change behavior. The increase in nanoparticle volume fraction smoothly decreased the charging time.
  • The porous medium’s average pore size had a negligible influence on the charging time. This design factor was the least critical parameter among design variables.

Author Contributions

Conceptualization, M.G., S.A.M.M., and P.T.; methodology, M.G., S.A.M.M., H.S., F.S., O.Y., K.I., G.A., and P.T.; software, S.A.M.M., H.S., and F.S.; validation, S.A.M.M. and M.G.; formal analysis, S.A.M.M., H.S., and F.S.; investigation, M.G., S.A.M.M., H.S., F.S., O.Y., K.I., G.A., and P.T.; resources, M.G. and O.Y.; writing—original draft preparation, M.G., S.A.M.M., H.S., F.S., O.Y., K.I., G.A., and P.T.; writing—review and editing, M.G., S.A.M.M., H.S., F.S., O.Y., K.I., G.A., and P.T.; visualization, S.A.M.M., H.S., and F.S.; supervision, M.G. Authorship must be limited to those who have contributed substantially to the work reported. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pomianowski, M.; Heiselberg, P.; Zhang, Y. Review of thermal energy storage technologies based on PCM application in buildings. Energy Build. 2013, 67, 56–69. [Google Scholar] [CrossRef]
  2. Iten, M.; Liu, S.; Shukla, A. A review on the air-PCM-TES application for free cooling and heating in the buildings. Renew. Sustain. Energy Rev. 2016, 61, 175–186. [Google Scholar] [CrossRef]
  3. Sahoo, S.K.; Das, M.K.; Rath, P. Application of TCE-PCM based heat sinks for cooling of electronic components: A review. Renew. Sustain. Energy Rev. 2016, 59, 550–582. [Google Scholar] [CrossRef]
  4. Du, K.; Calautit, J.; Eames, P.; Wu, Y. A state-of-the-art review of the application of phase change materials (PCM) in Mobilized-Thermal Energy Storage (M-TES) for recovering low-temperature industrial waste heat (IWH) for distributed heat supply. Renew. Energy 2021, 168, 1040–1057. [Google Scholar] [CrossRef]
  5. Jaguemont, J.; Omar, N.; Bossche, P.V.D.; Mierlo, J. Phase-change materials (PCM) for automotive applications: A review. Appl. Therm. Eng. 2018, 132, 308–320. [Google Scholar] [CrossRef]
  6. Bahraseman, H.G.; Languri, E.M.; East, J. Fast charging of thermal energy storage systems enabled by phase change materials mixed with expanded graphite. Int. J. Heat Mass Transf. 2017, 109, 1052–1058. [Google Scholar] [CrossRef]
  7. Liu, X.; Huang, Y.; Zhang, X.; Zhang, C.; Zhou, B. Investigation on charging enhancement of a latent thermal energy storage device with uneven tree-like fins. Appl. Therm. Eng. 2020, 179, 115749. [Google Scholar] [CrossRef]
  8. Nakhchi, M.; Esfahani, J. Improving the melting performance of PCM thermal energy storage with novel stepped fins. J. Energy Storage 2020, 30, 101424. [Google Scholar] [CrossRef]
  9. Tiari, S.; Qiu, S.; Mahdavi, M. Discharging process of a finned heat pipe–assisted thermal energy storage system with high temperature phase change material. Energy Convers. Manag. 2016, 118, 426–437. [Google Scholar] [CrossRef]
  10. Zhang, L.; Zhou, K.; Wei, Q.; Ma, L.; Ye, W.; Li, H.; Zhou, B.; Yu, Z.; Lin, C.-T.; Luo, J.; et al. Thermal conductivity enhancement of phase change materials with 3D porous diamond foam for thermal energy storage. Appl. Energy 2019, 233–234, 208–219. [Google Scholar] [CrossRef]
  11. Hu, N.; Li, H.; Wei, Q.; Zhou, K.; Zhu, W.; Zhang, L.; Li, S.; Ye, W.; Jiao, Z.; Luo, J.; et al. Continuous diamond-carbon nanotube foams as rapid heat conduction channels in composite phase change materials based on the stable hierarchical structure. Compos. Part B Eng. 2020, 200, 108293. [Google Scholar] [CrossRef]
  12. Talebizadehsardari, P.; Mohammed, H.I.; Mahdi, J.M.; Gillott, M.; Walker, G.S.; Grant, D.; Giddings, D. Effect of airflow channel arrangement on the discharge of a composite metal foam-phase change material heat exchanger. Int. J. Energy Res. 2021, 45, 2593–2609. [Google Scholar] [CrossRef]
  13. Sardari, P.T.; Babaei-Mahani, R.; Giddings, D.; Yasseri, S.; A Moghimi, M.; Bahai, H. Energy recovery from domestic radiators using a compact composite metal Foam/PCM latent heat storage. J. Clean. Prod. 2020, 257, 120504. [Google Scholar] [CrossRef]
  14. Zhao, C.; Opolot, M.; Liu, M.; Bruno, F.; Mancin, S.; Hooman, K. Numerical study of melting performance enhancement for PCM in an annular enclosure with internal-external fins and metal foams. Int. J. Heat Mass Transf. 2020, 150, 119348. [Google Scholar] [CrossRef]
  15. Zhang, X.; Sheikholeslami, M.; Jafaryar, M.; Sheremet, M.A.; Shafee, A.; Babazadeh, H. Simulation for melting of paraffin for saving energy with utilize of nanoparticles. J. Mol. Liq. 2020, 313, 113574. [Google Scholar] [CrossRef]
  16. Bondareva, N.S.; Gibanov, N.S.; A Sheremet, M. Melting of nano-enhanced PCM inside finned radiator. J. Phys. Conf. Ser. 2018, 1105, 012023. [Google Scholar] [CrossRef]
  17. Mahdi, J.M.; Mohammed, H.I.; Hashim, E.T.; Talebizadehsardari, P.; Nsofor, E.C. Solidification enhancement with multiple PCMs, cascaded metal foam and nanoparticles in the shell-and-tube energy storage system. Appl. Energy 2020, 257, 113993. [Google Scholar] [CrossRef]
  18. Al-Jethelah, M.; Ebadi, S.; Venkateshwar, K.; Tasnim, S.; Mahmud, S.; Dutta, A. Charging nanoparticle enhanced bio-based PCM in open cell metallic foams: An experimental investigation. Appl. Therm. Eng. 2019, 148, 1029–1042. [Google Scholar] [CrossRef]
  19. Choi, S.-K.; Kim, S.-O.; Lee, T.-H. Computation of the Natural Convection of Nanofluid in a Square Cavity with Homogeneous and Nonhomogeneous Models. Numer. Heat Trans. Part A Appl. 2013, 65, 287–301. [Google Scholar] [CrossRef]
  20. Nield, D.A.; Bejan, A. Convection in Porous Media; Springer: New York, NY, USA, 2013; Volume 3. [Google Scholar]
  21. Shahsavar, A.; Ali, H.M.; Mahani, R.B.; Talebizadehsardari, P. Numerical study of melting and solidification in a wavy double-pipe latent heat thermal energy storage system. J. Therm. Anal. Calorimet 2020, 141, 1785–1799. [Google Scholar] [CrossRef]
  22. Yang, X.; Yu, J.; Xiao, T.; Hu, Z.; He, Y.-L. Design and operating evaluation of a finned shell-and-tube thermal energy storage unit filled with metal foam. Appl. Energy 2020, 261, 114385. [Google Scholar] [CrossRef]
  23. Mahdi, J.M.; Nsofor, E.C. Melting enhancement in triplex-tube latent heat energy storage system using nanoparticles-metal foam combination. Appl. Energy 2017, 191, 22–34. [Google Scholar] [CrossRef]
  24. Sardari, P.T.; Mohammed, H.I.; Giddings, D.; Walker, G.S.; Gillott, M.; Grant, D. Numerical study of a multiple-segment metal foam-PCM latent heat storage unit: Effect of porosity, pore density and location of heat source. Energy 2019, 189, 116108. [Google Scholar] [CrossRef]
  25. Xu, Y.; Li, M.-J.; Zheng, Z.-J.; Xue, X.-D. Melting performance enhancement of phase change material by a limited amount of metal foam: Configurational optimization and economic assessment. Appl. Energy 2018, 212, 868–880. [Google Scholar] [CrossRef]
  26. Xiong, Q.; Poor, H.Z.; Izadi, M.; Assareh, E. Natural heat exchange in inhomogeneous porous medium using linear and quadratic porosity distribution. Int. J. Therm. Sci. 2021, 161, 106731. [Google Scholar] [CrossRef]
Figure 1. The physcial model of the storage unit: (a) cylindrical thermal energy storage (TES) system filled with phase change material (PCM) and metal foam, along with a heat transfer fluid tube in the center; (b) the computational domain.
Figure 1. The physcial model of the storage unit: (a) cylindrical thermal energy storage (TES) system filled with phase change material (PCM) and metal foam, along with a heat transfer fluid tube in the center; (b) the computational domain.
Energies 14 01575 g001
Figure 2. The variation of (a) melting volume fraction (MVF), and (b) the temperature at selected points with coordinates (5Ri, 0.5L) for various cases of grid size for ε a v g = 0.84 , V F n a = 0.04 , a = 0.6 , and P P I = 30 .
Figure 2. The variation of (a) melting volume fraction (MVF), and (b) the temperature at selected points with coordinates (5Ri, 0.5L) for various cases of grid size for ε a v g = 0.84 , V F n a = 0.04 , a = 0.6 , and P P I = 30 .
Energies 14 01575 g002
Figure 3. A comparison of melting fields simulated in the present simulations and the experimental study conducted by Al-Jethelah et al. [18].
Figure 3. A comparison of melting fields simulated in the present simulations and the experimental study conducted by Al-Jethelah et al. [18].
Energies 14 01575 g003
Figure 4. Comparisons of the predicted streamlines and isotherms with the simulations of Xiong et al. [26].
Figure 4. Comparisons of the predicted streamlines and isotherms with the simulations of Xiong et al. [26].
Energies 14 01575 g004
Figure 5. Mean values of the signal-to-noise (S/N) ratios for all the levels of the controlling parameters. Optimum case: A = 5, B = 1, C = 4, and D = 4.
Figure 5. Mean values of the signal-to-noise (S/N) ratios for all the levels of the controlling parameters. Optimum case: A = 5, B = 1, C = 4, and D = 4.
Energies 14 01575 g005
Figure 6. The time variation of (a) MVF, and (b) temperature distribution for various nano-additive volume fractions for Ro = 0.1 m, Ri = 0.01 m, L = 0.3 m, ε a v g = 0.800 , a = 0.3 , and P P I = 25 .
Figure 6. The time variation of (a) MVF, and (b) temperature distribution for various nano-additive volume fractions for Ro = 0.1 m, Ri = 0.01 m, L = 0.3 m, ε a v g = 0.800 , a = 0.3 , and P P I = 25 .
Energies 14 01575 g006
Figure 7. Melting interface at four steps of melting progress for various average porosities for V F n a = 0.04 , a = 0.3 , and P P I = 25 .
Figure 7. Melting interface at four steps of melting progress for various average porosities for V F n a = 0.04 , a = 0.3 , and P P I = 25 .
Energies 14 01575 g007
Figure 8. The variation of (a) MVF, and (b) temperature uniformity for the various average porosities as a function of time for V F n a = 0.04 , a = 0.3 , and P P I = 25 .
Figure 8. The variation of (a) MVF, and (b) temperature uniformity for the various average porosities as a function of time for V F n a = 0.04 , a = 0.3 , and P P I = 25 .
Energies 14 01575 g008
Figure 9. The variation of (a) MVF, and (b) temperature distribution at the various cases of pore per inch of the metal matrix as a function of time for V F n a = 0.04 , a = 0.3 , and ε a v g = 0.800 .
Figure 9. The variation of (a) MVF, and (b) temperature distribution at the various cases of pore per inch of the metal matrix as a function of time for V F n a = 0.04 , a = 0.3 , and ε a v g = 0.800 .
Energies 14 01575 g009
Figure 10. Melting interface at four steps of melting for various porosity gradients for V F n a = 0.04 , P P I = 25 , and ε a v g = 0.800 .
Figure 10. Melting interface at four steps of melting for various porosity gradients for V F n a = 0.04 , P P I = 25 , and ε a v g = 0.800 .
Energies 14 01575 g010
Figure 11. The variation of (a) MVF, and the (b) temperature distribution at various gradients of porosity as a function of time for V F n a = 0.04 , P P I = 25 , and ε a v g = 0.800 .
Figure 11. The variation of (a) MVF, and the (b) temperature distribution at various gradients of porosity as a function of time for V F n a = 0.04 , P P I = 25 , and ε a v g = 0.800 .
Energies 14 01575 g011
Figure 12. Variation of temperature at selected points during the charging process for optimum condition.
Figure 12. Variation of temperature at selected points during the charging process for optimum condition.
Energies 14 01575 g012
Figure 13. Variation of charging power during the charging process for optimum designs as a reference and the variation of (a) the volume fraction of the nano-additives, (b) the gradient of porosity for optimum conditions, (c) the average porosity, and (d) the number of pores per inch of the metal matrix.
Figure 13. Variation of charging power during the charging process for optimum designs as a reference and the variation of (a) the volume fraction of the nano-additives, (b) the gradient of porosity for optimum conditions, (c) the average porosity, and (d) the number of pores per inch of the metal matrix.
Energies 14 01575 g013
Table 1. Thermophysical properties of the coconut oil and the nanoadditives [18,19].
Table 1. Thermophysical properties of the coconut oil and the nanoadditives [18,19].
PropertiesCoconut Oil (Measured)Heat Transfer FluidCuO
Nanoadditives
Metal Foam
Solid Liquid
θ m e l t i n g (K)297Not applicableNot applicableNot applicableNot applicable
θ w i n d o w (K)2Not applicableNot applicableNot applicableNot applicable
lppcm (J·kg–1)-Not applicable1.03 × 104Not applicableNot applicable
ρ (kg·m–3)920993.7391465008900
μ l p p c m (N·s·m−2)Not applicable7.05 × 10–43.26 × 10–2Not applicableNot applicable
Cp (J·kg·K–1)375041782010540386
κ (W·m–1·K–1)2.28 × 10–16.23 × 10–11.66 × 10–118380
Table 2. Properties of the nano-enhanced phase change material.
Table 2. Properties of the nano-enhanced phase change material.
Density ρ n e p c m = ρ p p c m V F n a ρ p p c m + V F n a ρ n a
ρ p p c m ( θ ) = ξ ( θ ) ρ l p p c m ξ ( θ ) ρ s p p c m + ρ s p p c m
(18)
Dynamic viscosity μ l n p c m = ( μ l p p c m 0.4 μ l p p c m 0.4 V F n a ) 2.5 (19)
Coefficient of thermal expansion ρ l n p c m β n p c m = ρ l p p c m β l p p c m V F n a ρ l p p c m β l p p c m + V F n a ρ n a β n a (20)
Thermal conductivity κ l n p c m ( s n p c m ) = κ l p p c m ( s p p c m ) Π 1 Π 2
Π 1 = ( κ n a + 2 κ l p p c m ( s p p c m ) ) 2 V F n a ( κ l p p c m ( s p p c m ) κ n a )
Π 2 = ( κ n a + 2 κ l p p c m ( s p p c m ) ) + V F n a ( κ l p p c m ( s p p c m ) κ n a )
(21)
Heat capacity ( ρ c p ) l n e p c m ( s n e p c m ) = ( ρ c p ) l p p c m ( s p p c m ) + V F n a ( ρ c p ) , n a ( ρ c p ) l p p c m ( s p p c m ) (22)
Table 3. Initial and boundary conditions.
Table 3. Initial and boundary conditions.
Initial condition U r , n e p c m = U z , n e p c m = 0 ,   θ n e p c m = 293 K (23)
Interface of the tube wall and NePCM domain ( κ θ r ) H T F = ( κ θ r ) e f f , n e p c m ,   θ H T F = θ n e p c m (24)
Entrance of the tube θ H T F = 293 K ,   U r , H T F = 0 ,   U z , H T F = 0.01 m / s (25)
Outlet of the tube U r , H T F = 0 ,   ( θ z ) H T F = ( U z z ) H T F = 0 (26)
The right side surface of the NePCM domain U r , n e p c m = U z , n e p c m = 0 ,   ( θ r ) n e p c m = 0 (27)
Top and bottom surface of the NePCM domain U r , n e p c m = U z , n e p c m = 0 ,   ( θ z ) n e p c m = 0 (28)
Table 4. The range and levels of control factors.
Table 4. The range and levels of control factors.
FactorsDescriptionLevel 1Level 2Level 3Level 4Level 5
A V F n a (Volume fraction)0.000.010.020.030.04
B ε a v g (Average of porosity)0.8000.8250.8500.8750.900
Ca (Gradient of porosity)−0.6−0.30.00.30.6
DPPI (Pore per inch of the metal matrix)1015202530
Table 5. Taguchi L25 orthogonal table corresponding to range and levels of control parameters. PPI, pore size per inch; ES, stored energy; CP, unit power.
Table 5. Taguchi L25 orthogonal table corresponding to range and levels of control parameters. PPI, pore size per inch; ES, stored energy; CP, unit power.
Experiment NumberControl ParametersFull Melting CharacteristicsS/N Ratio
VFna ε a v g aPPI t | M V F = 1 ES
(kJ)
CP
(kW)
10.000.800−0.61010,200933.7450.09154−80.1720
20.000.825−0.31510,500965.1320.09192−80.4238
30.000.8500.02011,1001011.9590.09117−80.9065
40.000.8750.32512,1501071.5430.08819−81.6915
50.000.9000.63015,9001192.0820.07497−84.0279
60.010.800−0.3209900943.9270.09535−79.9127
70.010.8250.02510,350984.3720.09511−80.2988
80.010.8500.33010,9501025.5800.09366−80.7883
90.010.8750.61012,7501109.6300.08703−82.1102
100.010.900−0.61514,8501072.5190.07222−83.4345
110.020.8000.0309600950.4170.09900−79.6454
120.020.8250.31010,050987.7570.09828−80.0433
130.020.8500.61511,1001050.6370.09465−80.9065
140.020.875−0.62012,150975.1150.08026−81.6915
150.020.900−0.32512,9001023.7900.07936−82.2118
160.030.8000.3159300953.1800.10249−79.3697
170.030.8250.62010,0501007.8660.10029−80.0433
180.030.850−0.62511,250952.1280.08463−81.0231
190.030.875−0.33011,850997.4530.08417−81.4744
200.030.9000.01012,9001054.9670.08178−82.2118
210.040.8000.6259150963.7740.10533−79.2284
220.040.825−0.63010,500928.2510.08840−80.4238
230.040.850−0.31010,950967.6760.08837−80.7883
240.040.8750.01511,7001018.2230.08703−81.3637
250.040.9000.32013,2001085.4840.08223−82.4115
Table 6. The S/N and rank values of the control factors.
Table 6. The S/N and rank values of the control factors.
LevelsVFna ε a v g aPPI
1−81.44−79.67−81.35−81.07
2−81.31−80.25−80.96−81.10
3−80.90−80.88−80.89−80.99
4−80.82−81.67−80.86−80.89
5−80.84−82.86−81.26−81.27
δ0.623.190.490.38
Rank2134
Table 7. The optimum values of the controlling parameters.
Table 7. The optimum values of the controlling parameters.
FactorsTime for MVF = 1Estimated Time Fir MVF = 1ES
(kJ)
CP
(kW)
VFnaεavgaPPI90978784s931.5080.10240
0.040.80.325
Table 8. Table of 16 cases for further analysis around the optimum design.
Table 8. Table of 16 cases for further analysis around the optimum design.
Experiment No.CaseParameter StudyControl Parameters t | M V F = 1 ES
(kJ)
CP
(kW)
VFna ε a v g aPPI
11VFna0.000.8000.3259473978.6960.10195
220.010.8000.3259400966.3400.10226
330.020.8000.3259327965.6800.10219
440.030.8000.3259254953.1600.10249
55εavg0.040.8250.3259861975.1620.09850
660.040.8500.325106691012.8750.09378
770.040.8750.325116871039.6200.08886
880.040.9000.325131021085.4540.08223
99a0.040.800−0.6259796909.1590.09183
10100.040.800−0.3259530919.0700.09574
11110.040.8000.0259317937.6230.09922
12120.040.8000.6259150963.7740.10533
1313PPI0.040.8000.3109180952.5080.10242
14140.040.8000.3159180952.5190.10242
15150.040.8000.3209180952.5990.10243
16160.040.8000.3309180952.5770.10243
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ghalambaz, M.; Mehryan, S.A.M.; Shirivand, H.; Shalbafi, F.; Younis, O.; Inthavong, K.; Ahmadi, G.; Talebizadehsardari, P. Simulation of a Fast-Charging Porous Thermal Energy Storage System Saturated with a Nano-Enhanced Phase Change Material. Energies 2021, 14, 1575. https://doi.org/10.3390/en14061575

AMA Style

Ghalambaz M, Mehryan SAM, Shirivand H, Shalbafi F, Younis O, Inthavong K, Ahmadi G, Talebizadehsardari P. Simulation of a Fast-Charging Porous Thermal Energy Storage System Saturated with a Nano-Enhanced Phase Change Material. Energies. 2021; 14(6):1575. https://doi.org/10.3390/en14061575

Chicago/Turabian Style

Ghalambaz, Mohammad, S.A.M. Mehryan, Hassan Shirivand, Farshid Shalbafi, Obai Younis, Kiao Inthavong, Goodarz Ahmadi, and Pouyan Talebizadehsardari. 2021. "Simulation of a Fast-Charging Porous Thermal Energy Storage System Saturated with a Nano-Enhanced Phase Change Material" Energies 14, no. 6: 1575. https://doi.org/10.3390/en14061575

APA Style

Ghalambaz, M., Mehryan, S. A. M., Shirivand, H., Shalbafi, F., Younis, O., Inthavong, K., Ahmadi, G., & Talebizadehsardari, P. (2021). Simulation of a Fast-Charging Porous Thermal Energy Storage System Saturated with a Nano-Enhanced Phase Change Material. Energies, 14(6), 1575. https://doi.org/10.3390/en14061575

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop