Front. Phys.Frontiers in PhysicsFront. Phys.2296-424XFrontiers Media S.A.95702510.3389/fphy.2022.957025PhysicsOriginal ResearchEffects of Branched Fins on Alumina and N-Octadecane Melting Performance Inside Energy Storage SystemWeera et al.Significance of Different Branched FinsWeeraWajaree1*ManeengamApichit2SaeedAbdulkafi Mohammed34AissaAbderrahmane5GuedriKamel6YounisObai7MarzoukiRiadh89AsogwaKanayo K.101Department of Mathematics, Faculty of Science, Khon Kaen University, Khon Kaen, Thailand2Department of Mechanical Engineering Technology, College of Industrial Technology, King Mongkut’s University of Technology North Bangkok, Bangkok, Thailand3Department of Mathematics, College of Science, Qassim University, Buraydah, Saudi Arabia4Department of Mathematics, College of Education, Hodeidah University, Al-Hudaydah, Yemen5Laboratoire de Physique Quantique de la Matière et Modélisation Mathématique (LPQ3M), University of Mascara, Mascara, Algeria6Mechanical Engineering Department, College of Engineering and Islamic Architecture, Umm Al-Qura University, Mecca, Saudi Arabia7Department of Mechanical Engineering, College of Engineering at Wadi Addwaser, Prince Sattam Bin Abdulaziz University, Wadi Addwaser, Saudi Arabia8Chemistry Department, College of Science, King Khalid University, Abha, Saudi Arabia9Chemistry Department, Faculty of Sciences of Sfax, University of Sfax, Sfax, Tunisia10Department of Mathematics, Nigeria Maritime University, Okerenkoko, Nigeria
Edited by:Animasaun I. L., Federal University of Technology, Nigeria
Reviewed by:Abayomi Samuel Oke, Adekunle Ajasin University, Nigeria
H. A. Kumara Swamy, Presidency University, India
N. Keerthi Reddy, Ulsan National Institute of Science and Technology, South Korea
*Correspondence: Wajaree Weera, wajawe@kku.ac.th
This article was submitted to Interdisciplinary Physics, a section of the journal Frontiers in Physics
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
The importance of Phase change material (PCM) energy storage systems is no longer new in the industry. However, the influence of using branched fins inside the energy storage system on the melting process of alumina nanoparticles and n-octadecane has not been reported in the literature. Consequently, the outcome of a study on the numerical simulation for optimizing the melting performance of a PCM in various tubes, including those with branching fins is presented in this report. Four examples were assessed in relation to a suspension of alumina nanoparticles and n-octadecane paraffin that contains heated fins. A numerical technique based on the Galerkin finite element method (GFEM) was used to solve the dimensionless governing system. The average liquid percentage over the flow zone in question was computed. The primary results indicated that altering the number of heated fins might affect the flow structures, the system’s irreversibility, and the melting process. Case four, with eight heated fins, likewise produces the greatest average liquid fraction values and completes the melting process in 850s. Additionally, when 6% nano-enhanced PCM was used instead of pure PCM, the melting process is accelerated by 28.57 percent.
melting processphase change materialGFEM/XFEMfinsnusselt numberBackground Information
Energy storage technologies have risen in importance as a result of the unreliability of wind and solar energy sources, as well as their development in the efficiency and utilization of energy systems. Other PCM applications by [1] are of high importance in the industry due to the heat transfer mechanism of different PCM salts during phase change and of liquid salts. Also, using one storage module, the inclusion of PCM in solar cookers to extend usage time, and inclusion in a paint-drying system to recover exhaust heat are other important. [2] investigated thermal energy storage using curved-slab containers filled with phase change material and exposed to convective boundary conditions, taking into ac-count heat transfer fluid flows between the containers. One of the results of improving the performance of the latent heat thermal energy storage system is that a significant quantity of storage capacity may be obtained at greater efficiency. [3] found that the temperature range is attributed to possible CO2 emissions reductions that may be achieved by substituting cooling applications in residential and traditional heating. Diversifying energy supplies through the usage of suitable energy storage solutions and renewable energy sources, according to [4], are yardsticks for coming closer to energy sustainability.
[5] said that hot water tanks, which have a considerable influence on the tank’s heat loss and internal thermal stratification, are a classic example of the most often utilized thermal energy storage. In research by [6] on the importance of electromagnetic nanomaterials coating processes with intricate chemical reactions, the nanoparticle concentration magnitude was discovered to be capable of decreasing with an increase in the chemical reaction. [7] investigated the effects of nanoparticles on the dynamics of a fluid in a 3D cubic with two rotating cylinders in the center of the enclosure exposed to Lorentz force, angular velocity under various circumstances, and mixed convection transport phenomena. It was determined that when angular velocity increased, cumulative energy and average temperature decreased. [8] explored mixed convection of phase change material (PCM) in a 3D cavity filled with a revolving cylinder when a temperature gradient occurred between the hot right vertical wall and the cold left vertical wall while the remaining walls were deemed adiabatic. Phase change materials (PCM) are one of the most effective and active topics of study in terms of long-term heat energy storage and thermal management due to their great qualities. It is worth noting that phase change materials (PCM) can be used in conjunction with solar collectors to store excess solar energy while also regulating the temperature of photovoltaic solar collectors.
Interactions between cell/tissue systems and nanoparticles are no longer a surprise. According to [9], there is a link between respiratory illness, environmental exposure, and nano-sized particles. [10] recently published the results of another investigation on the addition of multi-walled carbon nanoparticles to the dynamics of water within a vertical Cleveland Z-staggered cavity, which demonstrated that raising the Reynold number increases the inertial force. [11] revealed that an optimal mix of nanoparticles is another essential feature in increasing thermal transfer in an annular geometry comprising nanoliquids with variably heated borders. According to [12], adding nanoparticles to the base fluid reduces the strength of fluid flow and the rate of heat transmission. Due to the incorporation of nanoparticles in the host fluid of water, [13] discovered that shallow annular enclosures provide higher thermal performance with little entropy formation. According to [14], the great thermal conductivity of alumina nanoparticles permits the nanoparticles to effectively diffuse heat from the host base fluids via Brownian motion. Aluminum oxide is useful to improve the rheological and filtration properties of drilling fluid. Some hybrid nanoliquids flow owing to buoyancy across a non-uniformly heated annulus, as examined by [15], Oke [16], [17], [18], and [19]. It was discovered that the proper mix of nanoparticles is a critical criterion for improving heat transfer. [20] emphasized the significance of grapheme nanoparticles due to their monomolecular layer of carbon atoms that are bonded by their remarkable and unique structural, optical, and electrical capabilities.
According to the findings of different studies on ternary-hybrid nanofluids by [21,22], [23], and [24], the viscosity and thermal conductivity of base fluids may be modified by the inclusion of nanoparticles. [25] investigated Darcy–Forchheimer nanoliquid convective flow in an odd-shaped cavity loaded with a multi-walled carbon nanotube-iron (II, III) oxide hybrid nanofluid when the walls are adiabatic and the internal and external cavity boundaries are isothermally is at low and high temperatures. In view of this, nothing is known on computational of alumina and n-octadecane melting performance inside energy storage system when the branched fins are two vertical but opposite fins, two horizontal but opposite fins, two vertical and two horizontal fins, and eight branched fins that are distributed. Research questions to be answered in this report are
1. At 500 s, what is the effect of the number of fins, liquid fraction, and Bejan number on the PCM melting process?
2. In the absence of volume fraction, what is the variations of temperature, liquid fraction, and Bejan number during the PCM melting process as time goes on?
3. At t=500s, how does increasing nanoparticle concentration, liquid fraction, and Bejan number affects alumina and n-octadecane melting performance in-side energy storage system?
4. What is the impact of fins number on the impact of fins number on liquid fraction, average Bejan number, average Nusselt number, and temperature distribution?
Research Methodology
The present computational model is shown in Figure 1. The heat transfer fluid is alumina and n-octadecane. Four distinct configurations of fins are modeled. The details of the four cases are shown in Figure 2. All fluids begin with a solidus temperature. The fined and cylinders were hot and adiabatic, respectively. The properties of nanoparticles and PCM are outlined in Table 1.
Computational setup of LTESS (A) 3-D model (B) 2-D cross-section.
Different cases of branched fins.
Properties of PCM and alumina.
Property
Al2O3
n-octadecane
ρ[kg/m3]
3,970
770
β×105[K−1]
0.85
91
k[w/mK]
40
0.157
L[j/kg]
−
242,9
Fusion [C]
−
28×103
μ×103[ Pa.s ]
−
3.79
Cp[j/kgK]
765
2,189
Sources [7] and [22].
Problem Formulation
The GFEM technique was employed to simulate the transient flow that could be described as Newtonian and laminar. Boussinesq estimation was utilized to account for the gravitational force effect. The continuity (Eq. 1), momentum in x and y direction (Eqs 2, 3), and energy (Eq. 4) equations read [26] and [27].∇⋅V→=0(∂v∂t+V→⋅∇v)=vC(λ−1)2ε+λ3+1ρnf(−∇P+μnf∇2v)+1ρnf(ρβ)nfg(T−Tref)∂u∂t+V→⋅∇u=uC(λ−1)2ε+λ3+1ρnf(−∇P+μnf∇2u)(ρCp)nf∂(ρLλ)nf∂t+(ρCp)nf∂T∂t−knf∇2T=−(ρCp)nfV→⋅∇TThe boundary conditions associated with Eqs 1–4 for the inner hot wallU=V=0,θ=1For the outer wallU=V=0,∂θ∂X=0We considered ε=10−3,C=105. A single-phase model was used to predict attributes Refs. [16 - 17](ρCp)f−1(ρCp)nf=(1−ϕ)+ϕ(ρCp)s(ρCp)f−1,ρnf=ϕρs+ρf(1−ϕ),μnf=μf(1−ϕ)2.5(ρβ)nf=ϕ(ρβ)s+(1−ϕ)(ρβ)f,(ρL)f=(ρL)nf(1−ϕ),knf=2kf+2ϕ(ks−kf)+kpkp−ϕ(ks−kf)+2kfkf
Enthalpy
The specific heat capacity of the PCM changes significantly with its phase. The specific enthalpy of the PCM is thus determined as a function of its temperature by [28] ash=href+∫TretT(Cp)nfdTLiquid fraction is introduced asλ={1T≤TlT−TsTl−TsTs≤T≤Tl,He=h+λL0T≤TsThe formula of Sgen,total,Sgen,th,Sgen,fandSgen,m by [29]Sgen,total=Sgen,th+Sgen,f+Sgen,mSgen,total=knfT2[(∂T∂x)2+(∂T∂y)2]+μnfT{2[(∂ux∂x)2+(∂uy∂y)2]+(∂ux∂y+∂uy∂x)2}+NμσnfσfHa2V2WhereNμ=μfT0kf(αfL(∆T))2Where Sgen,thdenotes the entropy production resulting from heat transfer irreversibility (HTI), Sgen,fdenotes the entropy production resulting from fluid friction irreversibility (FFI), and Sgen,m denotes the entropy production resulting from magnetic field impacts. And consequently the mean Bejan numbers be determined asBe=Sgen,thSgen,total
GFEM Approach
The transformed coupled Eqs 1–4, which comprise both the motion and heat transport processes as well as the requisite boundary conditions, are handled using the Galerkin Finite Element Method. First, on a non-uniform structural grid, the weak forms of the governing equations are presented and discretized. After that, mathematical software is used to model the outcomes and gives a detailed description of the procedure. It is crucial to keep in mind that the governing equations and their related constraints were solved by employing the Galerkin finite element technique. The programming environment is divided into triangle-shaped sections. Triangular Lagrange finite elements of various orders are employed on all the flow variables within the computational domain. By substituting the governing equations for the approximations, the residue is obtained.|Γi+1−ΓiΓi|≤10−6
Finally, using further numerical findings from [30], the validation of the current code is produced and illustrated in Figure 3. Based on this graph, we may be confident in our findings.
Mesh grid and present model in comparison with [30].
Results and Discussion
The findings of studying the melting impacts on the motion of a suspension containing phase change material (PCM) are depicted and analyzed in this section. In this case, the heat transfer fluid is alumina and n-octadecane, and the motion area is a square tube with cross-section wings. For different heating scenarios, The isotherms lines, streamlines, Bejan number, and molten liquid of several cylinders are explored, including cylinders with two horizontal fins, cylinders with two vertical fins, cylinders with four fins, and cylinders with eight heated wings. Temporal Variations between 100 and 1000 s are investigated, and the nanoparticles volume percentage values are considered between 0 and 8%. The mean values of the molten PCM, Bejan number Beavg, and Nusselt number Nuavg versus the progression of time are graphically displayed over a broad range of the studied parameters to offer a complete examination. The condition of the completely melted (β=1) may also be utilized to stop the calculations. For different scenarios of inner heating, Figure 4 shows the isotherms lines, streamlines, Bejan number, and molten PCM percentage. It should be noticed that the temperature characteristics are focused around the fins in all instances, indicating a cold area around the outer cubic lowest. Case 4 (eight wings) had the highest temperature distributions, indicating a reduction in the above-mentioned cold area at the lowest.
Effect of the number of fins, liquid fraction, and Bejan number during the PCM melting process at 500 s.
It has also been noted that as the number of fins increases, the temperature differences decrease, so the temperature drop and thermal performance rate decrease. In addition, when the number of fins is improved, the velocity values clearly decrease. Physically, increasing the number of fins increases the complexity of the motion region, which increases motion resistance. Simultaneously, the characteristics of the Bejan number reveal that increasing the number of heated fins diminishes temperature gradients, and therefore fluid friction irreversibility takes over. In addition, the melting area appears in the top half of the area for all of the variables investigated, and the number of heated wings increases the volume of the melting zone. Figure 5 depicts the characteristics of isotherms lines, streamlines, Be number, and mol-ten PCM percentage as time passes. Case 3 is employed in these calculations, using an inner cylinder with four wings. The findings showed that at the start of the computations (short time values), isotherms lines, streamlines, and Bejan number distributions occur around the inner heated zone, suggesting a non-active area towards the outer limits. The liquid begins to convey and disperse the temperature across the area as time passes. As a result, a suitable thermal domain with a greater velocity percentage towards the lowest of the outer borders is produced at time = 600. Furthermore, when comparing the fluid friction irreversibility to the heat transmission irreversibility, the higher the time value, the greater the dominance of the fluid friction irreversibility at the bottom. When the time values grow, the mushy zone appears across the whole flow domain.
Variations of temperature, liquid fraction, and Bejan number during. the PCM melting process as time goes on when ϕ=0.
Figure 6 depicts the isotherms lines, streamlines, Bejan number, and liquid fraction distributions as a function of the volume fraction parameter. In this situation, inner heated cylinders with four wings are employed. Due to the rise in the viscosity of the mixture, a poor convective transport is seen at higher values of. The data show that as it expands, the velocity and temperature gradients decrease. At low values of, the Be number occurs around the fins rather than the bottom borders. On the contrary, increasing increases the mushy zone inside the motion region until the circumstances are entirely melted at 0.04. Figure 7 and Figure 8 show the mean molten PCM rate, Beavg, and Nuavg profiles as a function of the heated wings number, duration parameter, and concentration parameter.
Effect of nanoparticle concentration, liquid fraction, and Bejan number. during the PCM melting process at t=500s.
The impact of fins number on liquid fraction, Average Bejan number, Average Nusselt number, and temperature distribution.
The impact of nanoparticles concentration on molten liquid, Average Nusselt number (Nuavg) and Average Nusselt number Average Nusselt number (Beavg).
The findings indicated that example 4, which assumes eight heated fins, had the highest β values. However, the Beavg and Nuavg values decline when the number of heated wings increases. Furthermore, the average heat transfer rate decreases when the temperature gradients reduce. Furthermore, greater values of cause the irreversibility of heat transmission to take precedence over the irreversibility of fluid friction. Finally, increasing the nanoparticles concentration improves the mushy zone, increasing the β.
Conclusion
The current research numerically investigated the influence of inserting different numbers of branched fins inside energy storage system on the melting process of alumina and n-octadecane. Based on the number of heated fins, four instances were considered. The shaky scenario was taken into account, and fully melted situations were anticipated. The governing system was solved using the finite element technique (GFEM) and the Poisson pressure equation. The following are some of the most important findings:
1. As the heated fins increase due to the enhancement in the buoyancy-convective situation, the temperature, velocity, and Bejan number distributions increase. In scenario 4, the melted area was also regulated throughout the majority of the flow domain.
2. Isotherms lines, streamlines, and molten PCM occur around the inner heated forms at short values of time, whereas as time passes, a nice isothermal and melted motion domain was produced.
3. As it rises, the dynamic viscosity of the mixture increases, and as a result, the velocity of the mixture decreases.
4. When compared to heat transmission irreversibility, the irreversibility owing to fluid friction becomes more dominant with time.
5. Case four, with eight heated fins, likewise produces the greatest average liquid fraction values and completes the melting process in 850 s.
6. When 6% nano-enhanced PCM was used instead of pure PCM, the melting process was found to accelerate by 28.57 percent.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author Contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Funding
Authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through research group under grant number R.G.P.2/224/43. The authors would like to thank the Deanship of Scientific Research at Umm Al-Qura University for supporting this work by Grant Code: (22UQU4331317DSR21).
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
ReferencesZalbaB.Marı́nJMCabezaL. F.MehlingH.Review on thermal Energy Storage with Phase Change: Materials, Heat Transfer Analysis and Applications. Appl Therm Eng (2003) 23(3):251–83. 10.1016/s1359-4311(02)00192-8KoizumiHJinY. Performance Enhancement of a Latent Heat thermal Energy Storage System Using Curved-Slab Containers. Appl Therm Eng (2012) 37:145–53. 10.1016/j.applthermaleng.2011.11.009Pereira da CunhaJEamesP. Thermal Energy Storage for Low and Medium Temperature Applications Using Phase Change Materials – A Review. Appl Energy (2016) 177:227–38. 10.1016/j.apenergy.2016.05.097DinçerIRosenM A.Renewable Energy Systems with Thermal Energy Storage. Therm Energy Storage (2021) 521–55. 10.1002/9781119713173.ch7CabezaL. F.Advances in Thermal Energy Storage Systems: Methods and Applications. Advanc in Therm Energy Storage Systems (2021) 37–54. 10.1016/b978-0-12-819885-8.00002-4ShahidA.BhattiM. M.Anwar BégO.AnimasaunI. L.JavidK.Spectral Computation of Reactive Bi-directional Hydromagnetic Non-Newtonian Convection Flow from a Stretching Upper Parabolic Surface in Non-Darcyporous Medium. Int J Mod Phys B (2021) 35(29):2150294. 10.1142/s0217979221502945AbderrahmaneA.HatamiM.YounisO.MouradA.Effect of Double Rotating Cylinders on the MHD Mixed Convection and Entropy Generation of a 3D Cubic Enclosure Filled by Nano-PCM. Eur Phys J Spec Top. (2022). 10.1140/epjs/s11734-022-00586-7Al-KouzW.AissaA.Uma DeviS.PrakashM.KolsiL.MoriaH.Effect of a Rotating cylinder on the 3D MHD Mixed Convection in a Phase Change Material Filled Cubic Enclosure. Sustain Energy Technologies and Assessments (2022) 51:101879.10.1016/j.seta.2021.101879TortorellaS.KaragiannisT. C.The Significance of Nanoparticles in Medicine and Their Potential Application in Asthma. In: MaulikNKaragiannisT, editors. Molecular Mechanisms and Physiology of Disease. New York, NY: Springer (2014). p. 247–75. 10.1007/978-1-4939-0706-9_10RasoolG.SaeedA. M.LareA. I.AbderrahmaneA.GuedriK.VaidyaH.Darcy-Forchheimer Flow of Water Conveying Multi-Walled Carbon Nanoparticles through a Vertical cleveland Z-Staggered Cavity Subject to Entropy Generation. Micromachines (2022) 13:744. 10.3390/mi13050744ReddyN. K.SwamyH. A. K.SankarM.Buoyant Convective Flow of Different Hybrid Nanoliquids in a Non-uniformly Heated Annulus. Eur Phys J Spec Top. (2021) 230(5):1213–25. 10.1140/epjs/s11734-021-00034-ySankarM.SwamyH. A. K.DoY.AltmeyerS.Thermal Effects of Nonuniform Heating in a Nanofluid‐filled Annulus: Buoyant Transport versus Entropy Generation. Heat Trans (2022) 51(1):1062–91. 10.1002/htj.22342SwamyH. A. K.SankarM.ReddyN. K.Analysis of Entropy Generation and Energy Transport of Cu-Water Nanoliquid in a Tilted Vertical Porous Annulus. Int J Appl Comput Math (2022) 8(1):1–23. 10.1007/s40819-021-01207-ySongY.-Q.ObideyiB. D.ShahN. A.AnimasaunI. L.MahrousY. M.ChungJ. D.Significance of Haphazard Motion and Thermal Migration of Alumina and Copper Nanoparticles across the Dynamics of Water and Ethylene Glycol on a Convectively Heated Surface. Case Studi Thermal Engineering (2021) 26:101050. 10.1016/j.csite.2021.101050Keerthi ReddyN.SankarM.Buoyant Convective Transport of Nanofluids in a Non-uniformly Heated Annulus. J Phys Conf. Ser(2020) 1597(1):012055. 10.1088/1742-6596/1597/1/012055OkeA. S.Heat and Mass Transfer in 3D MHD Flow of EG-Based Ternary Hybrid Nanofluid over a Rotating Surface. Arab J Sci Eng (2022) 1-17. in-press. 10.1007/s13369-022-06838-xOkeA. S.MutukuW. N.Significance of Viscous Dissipation on MHD Eyring-Powell Flow Past a Convectively Heated Stretching Sheet. Pramana J Phys (2021) 95(4):1–7. 10.1007/s12043-021-02237-3ReddyN. K.SwamyH. A.SankarM.Buoyant Convective Flow of Different Hybrid Nanoliquids in a Non-uniformly Heated Annulus. European Physical Journal Special Topics (2021) 230(5):1213–25. 10.1140/epjs/s11734-021-00034-yMebarek-OudinaF.Keerthi ReddyN.SankarM.Heat Source Location Effects on Buoyant Convection of Nanofluids in an Annulus. In: Advances in Fluid Dynamics. Singapore: Springer (2021). p. 923–37. 10.1007/978-981-15-4308-1_70SharmaR.RajuC. S.AnimasaunI. L.SanthoshH. B.MishraM. K.Insight into the Significance of Joule Dissipation, thermal Jump and Partial Slip: Dynamics of Unsteady Ethelene Glycol Conveying Graphene Nanoparticles through Porous Medium. Nonlinear Engineering (2021) 10(1):16–27. 10.1515/nleng-2021-0002AnimasaunI. L.YookS.-J.MuhammadT.MathewA.Dynamics of Ternary-Hybrid Nanofluid Subject to Magnetic Flux Density and Heat Source or Sink on a Convectively Heated Surface. Surfaces and Interfaces (2022) 28:101654. 10.1016/j.surfin.2021.101654AnimasaunI. L.ShahN. A.WakifA.MahantheshB.SivarajR.KorikoO. K.Ratio of Momentum Diffusivity to Thermal Diffusivity: Introduction, Meta-Analysis, and Scrutinization. New York: Chapman and Hall/CRC (2022). ISBN-13: 978-1032108520, ISBN-10: 1032108525, ISBN9781003217374.SaleemS.AnimasaunI. L.YookS.-J.Al-MdallalQ. M.ShahN. A.FaisalM.Insight into the Motion of Water Conveying Three Kinds of Nanoparticles Shapes on a Horizontal Surface: Significance of Thermo-Migration and Brownian Motion. Surfaces and Interfaces (2022) 30:101854. 10.1016/j.surfin.2022.101854CaoW.I.L.A.YookS.-J.V.A.O.JiX.Simulation of the Dynamics of Colloidal Mixture of Water with Various Nanoparticles at Different Levels of Partial Slip: Ternary-Hybrid Nanofluid. International Communicatioin Heat and Mass Transfer (2022) 135:106069. 10.1016/j.icheatmasstransfer.2022.106069Al-KouzWAissaAKoulaliAJamshedWMoriaHNisarKSMHD Darcy-Forchheimer Nanofluid Flow and Entropy Optimization in an Odd-Shaped Enclosure Filled with a (MWCNT-Fe3O4/water) Using Galerkin Finite Element Analysis. Sci Rep (2021) 11(1):22635–15. 10.1038/s41598-021-02047-ySheikholeslamiM.JafaryarM.ShafeeA.LiZ.Hydrothermal and Second Law Behavior for Charging of NEPCM in a Two Dimensional thermal Storage Unit. Chinese Journal of Physics (2019) 58:244–52. 10.1016/j.cjph.2019.02.001MouradA.AissaA.SaidZ.YounisO.IqbalM.AlazzamA.Recent Advances on the Applications of Phase Change Materials for Solar Collectors, Practical Limitations, and Challenges: A Critical Review. Journal Energrgy Storage (2022) 49:104186. 10.1016/j.est.2022.104186BouzennadaT.MechighelF.FilaliA.GhachemK.KolsiL.Numerical Investigation of Heat Transfer and Melting Process in a PCM Capsule: Effects of Inner Tube Position and Stefan Number. Case Stud Thermal Engineering (2021) 27:101306. 10.1016/j.csite.2021.101306AbderrahmaneA. MHD Hybrid Nanofluid Mixed Convection Heat Transfer and Entropy Generation in a 3-D Triangular Porous Cavity with Zigzag Wall and Rotating Cylinder. Mathematics (2022) 10(5):769. 10.3390/math10050769TanF. L.HosseinizadehS. F.KhodadadiJ. M.FanL.Experimental and Computational Study of Constrained Melting of Phase Change Materials (PCM) inside a Spherical Capsule. Int J Heat and Mass Transfer (2009) 52:3464–72. 10.1016/j.ijheatmasstransfer.2009.02.043NomenclatureCp(J/kg.K)