Effect of dual-rotation on MHD natural convection of NEPCM in a hexagonal-shaped cavity based on time-fractional ISPH method

The time-fractional derivative based on the Grunwald–Letnikove derivative of the 2D-ISPH method is applying to emulate the dual rotation on MHD natural convection in a hexagonal-shaped cavity suspended by nano-encapsulated phase change material (NEPCM). The dual rotation is performed between the inner fin and outer hexagonal-shaped cavity. The impacts of a fractional time derivative \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}α \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0.92 \le \alpha \le 1} \right)$$\end{document}0.92≤α≤1, Hartmann number Ha \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0 \le Ha \le 80} \right)$$\end{document}0≤Ha≤80, fin length \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0.2 \le L_{Fin} \le 1} \right)$$\end{document}0.2≤LFin≤1, Darcy parameter Da \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {10^{ - 2} \le Da \le 10^{ - 4} } \right)$$\end{document}10-2≤Da≤10-4, Rayleigh number Ra \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {10^{3} \le Ra \le 10^{6} } \right)$$\end{document}103≤Ra≤106, fusion temperature \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta_{f}$$\end{document}θf \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0.05 \le \theta_{f} \le 0.8} \right)$$\end{document}0.05≤θf≤0.8, and solid volume fraction \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document}φ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0 \le \varphi \le 0.06} \right)$$\end{document}0≤φ≤0.06 on the velocity field, isotherms, and mean Nusselt number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{Nu}$$\end{document}Nu¯ are discussed. The outcomes signaled that a dual rotation of the inner fin and outer domain is affected by a time-fractional derivative. The inserted cool fin is functioning efficiently in the cooling process and adjusting the phase change zone within a hexagonal-shaped cavity. An increment in fin length augments the cooling process and changes the location of a phase change zone. A fusion temperature \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta_{f}$$\end{document}θf adjusts the strength and position of a phase change zone. The highest values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{Nu}$$\end{document}Nu¯ are obtained when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha = 1$$\end{document}α=1. An expansion in Hartmann number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Ha $$\end{document}Ha reduces the values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{Nu}$$\end{document}Nu¯. Adding more concentration of nanoparticles is improving the values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{Nu}$$\end{document}Nu¯.

The nano-encapsulated phase change material (NEPCM). The NEPCMs consist of the core and the shell. The cores are made of a Phase Change Material (PCM). It solidifies or melts at a special temperature called s fusion temperature. The shell consists of a polymer. The NEPCM suspensions are a new form of hybrid nanofluids. It has several heat transfer applications 1,2 . The convection of heat and fluid flow of NEPCM have been reported in wide theoretical studies such as enclosures [3][4][5][6] , divergent heatsink 7,8 , energy storage 9 , with fins 10 , and others 11,12 . Shafee et al. 13 used Galerkin approach of the finite element method to handle the phase change process of NEPCM in a heat storage. Selimefendigil et al. 14 analyzed numerically the natural convection of CuOwater nanofluid in a square cavity with a conductive partition and a phase change material (PCM) under the effect of a uniform inclined magnetic field. The mixed convection in a phase change material-filled a square cavity under the effect of a rotating cylinder was numerically investigated by Selimefendigil and Öztop 15 . In further studies, Selimefendigil et al. 16,17 presented different numerical studies on the phase change dynamics of a 3D cylinder containing hybrid nanofluid and phase change material (PCM) by using the finite element solver.
Internal fins. As mentioned in the previous part, the use of phase change material (PCM) has acquired highly increasing concern in various engineering applications, through thermal administration systems, solar energy storage, and conservation of the energy in buildings 18,19 . But the PCMs have low thermal conductivity obstruct the heat transfer during solidification or melting processes, which break down the efficiency of energy storage. So, different methods have been suggested to enhance the heat transfer of the PCM. One of them is adding the internal fins. Different geometric parametric of the internal fins like their thickness, length, location, and inclination angle have been investigated. Ren and Chan 20 found that the longer fins were more active than the shorter fins to increase the melting rate of PCM. Sciacovelli et al. 21 examined differently shaped fins to enhance the heat inside the cylindrical cavity. They reported that the tree-shaped fin was increased the system energy efficiency by 24%. Ji et al. 22 studied the effects of fins inclination angles on the melting rate of PCM. Li and Yu 23 investigated the influence of the internal fins on the melting process, the design of dual fins, and tree-shaped fins. They illustrated how changing the fins structures better than increasing the number of fins in the offers a higher rate of heat transfer and a better energy storage capacity. More different configurations of fixed or flexible fins, and their design and geometric factors can be found in these studies [24][25][26][27][28][29] . www.nature.com/scientificreports/ Rotating cavity. The study of heat and fluid flow inside rotating cavities is significant from both theoretical as well as application sides of view. These kinds of studies are more complex due to the rotation, and different body forces have driven the flow. As examples of its applications, in astrophysical and geophysical flows, for semiconductors in manufacturing of single wafer crystal, to storage the thermal energy in rotating systems, to cooling the tools microelectronic, and so on. Earlier studies were presented to investigate the effects of the fluid flow on the flow and heat transfer rate inside the rotating enclosures [30][31][32][33][34][35][36][37][38] . Mandal and Sonawane 39 studied the flows inside a differentially heated rotating square cavity in two different formulations. They found that the force of inertia appears effect by the increasing speed of the enclosure rotation. Other studies of natural convection and surface radiation in a rotating square cavity with low rotation velocity were presented by Mikhailenko et al. 40-43 . The fractional derivatives. When the order of the derivatives and integrals are non-integers in the classical partial differential equations, then we have a new kind of differential equations called fractional differential equations (FPDEs). The importance of fractional differential equations comes from their wide applications in engineering and science. These equations can be used to simulate the problems in fluid dynamics, electrochemistry, electrodynamics, nanotechnology, astronomy sciences, and chemical physics. We can summarize the famous types of fractional derivatives as: From these types, the researchers choose the type, which is more compatible with the experimental results, when they want to solve the system's equations of the physical problems.
Many good studies were presented and published to illustrate the applications of fractional calculus to transport processes [47][48][49][50][51][52] . In general, the mesh-free nature of the ISPH method helps in handling the high deformation and fluid-structure interaction problems. So, the ISPH method is adopted for the current problem of the rotating paddle wheel inside a novel geometry of a cross-shaped cavity. The objective of this investigation is to employ the time-fractional derivative in the solving steps of the ISPH method. The dual rotation and inclined magnetic impacts on the natural convection of NEPCM embedded in a hexagonal-shaped cavity are conducted. The dual rotation between an inner fin and outer hexagonal-shaped cavity during natural convection flow can be applied in generating thermal energy from the rotating systems and cooling process of the electronic devices. The results indicated that the fractional time derivative changes the dual rotation between the inner fin and outer domain. The inserted cool fin is functioning efficiently in the cooling process and adjusting the phase change zone within a hexagonal-shaped cavity. A fusion temperature alters the strength and position of a phase change zone. The highest values of Nu are found at higher value of the time-fractional derivative (α = 1) . Increasing Hartmann number reduces the values of Nu , whilst increasing solid volume fraction enhances the values of Nu.

Mathematical analysis
The basic illustration of the present physical problem has been shown in Fig. 1. The inner fin is cooled by a temperature T c and the two rectangles in the flat walls of a cavity is heated by a temperature T h . The embedded fin is rotating clockwise, and the outer hexagonal-shaped domain is rotating anticlockwise. It is presumed the latent heat is almost 211 kJ kg −1 whilst the fusion temperature of the PCM cores is about 32 °C. The dual rotation is carrying a uniform circular velocity around the center of a cavity. Table 1 presents the physical attributes of a porous matrix and a mixture fluid. The local thermal equilibrium model is assumed amongst the mixture nanofluid and a porous medium. The dimensional governing equations 53,54 are: where (u, v) are the dimensional velocity vector, ε is a porosity, ρ is a density, g is gravity, µ is a dynamic viscosity, β is a thermal expansion, and K is a permeability. Furthermore, Eq. (3) indicates the local thermal equilibrium condition amongst the mixture fluid and the solid matrix. According to Ghalambaz et al. 53,55,56 , the mixture density is: where ρ p is a density of NEPCM particles. ρ p is calculated from a density of a core ρ c and a shell ρ s as:   53 . where χ ≈ 0.447 is a core-shell weight ratio for NEPCMs. The Specific heat capacity is calculated as: The heat capacity is: where (C p ) c,l and (C p ) s are a heat capacity of a core and a shell. The sine profile of the latent heat of phase change is: where The thermal conductivity, thermal expansion, and dynamic viscosity are: N 1 is a dynamic viscosity number and N 2 is a thermal conductivity number. The applied dimensionless quantities 54 are: The dimensionless regulating equations in Lagrangian type 53,54 are: The dimensionless parameters are Raleigh number Ra = The velocities of an inner fin: where a dimensionless angular velocity ω is kept at 2.5.
The boundary conditions: with where L ctot is a total length of the cold walls. n is a normal vector. The references 5,53,57,58 are summarizing the NEPCM and mixture fluid properties.

ISPH method
The solver steps based on the time-fractional derivative are: Step 1:

Pressure Poisson equation (PPE):
Corrected velocities: The thermal equation is: Rectangle -shapes in flat walls of a cavity: Flat walls of a cavity: Other cavity walls: The shifting technique 59 is:

Verification tests
This section checks the efficiency of the ISPH method in simulating the natural convection resultant from an inner circular cylinder. Figure 2 shows the comparison of the isotherms and streamlines between the results of Kim et al. 60 and the ISPH method at Rayleigh number Ra = 10 3 , 10 4 , 10 5 and 10 6 . It is remarked that the present results of the ISPH method agree well the reference of Kim et al. 60 . The results are presented at the steady-state and convergence criteria for the current verification is taken as 10 −6 . Further studies on the verification tests of the natural/mixed convection flow using the ISPH method can be found in references [61][62][63][64][65] .
Further validation of the ISPH method with the numerical results of for natural convection in a cavity at Rayleigh number Ra = 10 5 is shown in Fig. 3. The validation results are providing adequate confidence of the ISPH method.

Results and discussion
The current research aims to apply a fractional time derivative in the ISPH method for emulating the dual rotation between an inner cross fin and outer hexagonal-shaped domain. The hexagonal-shaped cavity is containing two rectangle heat sources on the flat walls and is suspended by NEPCM. During the executed simulations, the parameters are fixed at dimensionless angular velocity ω = 2.5 , a porosity parameter ε = 0.6 , Stefan number Ste = 0.2 , and a magnetic inclination angle γ = 45 • . The influences of a fractional time derivative (0.92 ≤ α ≤ 1) , Hartmann number (0 ≤ Ha ≤ 80) , the fin length (0.2 ≤ L Fin ≤ 1) , Darcy parameter 10 −2 ≤ Da ≤ 10 −4 , Rayleigh number 10 3 ≤ Ra ≤ 10 6 , a fusion temperature 0.05 ≤ θ f ≤ 0.8 , and solid volume fraction (0 ≤ ϕ ≤ 0.06) on the velocity field, isotherms, and mean Nusselt number Nu are discussed. The dimensionless angular velocity is lowering at ω = 2.5 to consider the natural convection mode only during the simulations. Figure 4 shows the velocity field, isotherms, and heat capacity under the variations of a fractional time derivative α. It is remarked that the dual rotation between an inner fin and outer hexagonal-shaped domain is affected by the variations on a fractional time derivative α . Consequently, the velocity field, and isotherms are varied according to the variations on α . Thus, the zone of a phase change material (PCM) is influenced by the location of the outer domain and inner shape. The current investigation reported that the PCM is changed as the fractional time derivative α is varied. Here, the factor α is playing a significant role in controlling the rotation speed between inner/outer shapes, nanofluid movements, and heat transfer in a hexagonal-shaped cavity. Figure 5 signifies the impacts of Hartmann number Ha on the velocity field, isotherms, and heat capacity. Physically, an increase in Hartmann number reduces the convection flow and suppresses the fluid flow due to the Lorentz force. In this model, due to the presence of a cold fin in the cavity's center and the two rectangle heaters in the flat walls, the contributions of the Hartmann number are less. As a result, an increment in Ha is giving a minor impact on the velocity field, isotherms, and heat capacity within a cavity. Figure 6 indicates the effects of the fin length L Fin on the velocity field, isotherms, and heat capacity. The inner fin acts as a prominent character in the cooling process inside a hexagonal-shaped cavity, the variations on the fin length L Fin are changing the nanofluid movements in a cavity. As the cool fin represents a blockage within a hexagonal-shaped cavity, an increase in L Fin declines the velocity field. Increasing the fin length L Fin augments www.nature.com/scientificreports/ the cooling area and accordingly the temperature distributions are reduced. The PCM is affected clearly by the variations in the fin length L Fin . It is observed that the expansion in the fin length controls the location of a phase change zone. Figure 7 introduces the influences of Darcy parameter Da on the velocity field, isotherms, and heat capacity. Physically, the Darcy parameter signifies the major element of the porous resistance for the fluid flow. Decreasing Da from 10 −2 to 10 −4 , leads to a decrease in the velocity's maximum by 8.69%. Further, a decrease in Da leads to a minor reduction in the temperature distributions. Thus, due to the minor change in the temperature distributions below the variations on Da , the phase change zone is affected slightly by variations on Da. Figure 8 shows the velocity field, isotherms, and heat capacity under the variations of Rayleigh number Ra . Physically, the Rayleigh number augments the buoyancy forces which powers the fluid flow and heat transfer within a cavity. The strength of the velocity field is increasing strongly as Ra increases. Further, a growth in Ra strengthens the temperature distributions in a hexagonal-shaped cavity. Thus, the heat capacity is affected clearly by an increase in Ra . The physical meaning of a high Ra is powering the buoyancy-driven flow. Figure 9 presents the fusion temperature θ f impacts on the heat capacity. It is remarked that an increment of θ f reduces a phase change zone. Further, increasing in θ f closes the phase change zone near the rectangle heaters www.nature.com/scientificreports/ of a hexagonal-shaped cavity. The physical explanation returns to the connection between a heat capacity Cr and a fusion temperature θ f . Figure 10 introduces the impacts of a solid volume fraction ϕ on the contours of the velocity field, and isotherms. Physically, adding more concentrations of the nanoparticles boosts the viscosity of the mixture fluid, and accordingly, the velocity field is decreasing. There are minor changes in the isotherms according to adding more concentrations of the nanoparticles. The fewer contributions of ϕ on the heat transfer return to the presence of a cool fin inside a hexagonal-shaped cavity. Figure 11 indicates the profiles of Nu under the effects of a fractional derivative parameter α , a fusion temperature θ f , Hartmann number Ha , and solid volume fraction ϕ . Initially, the profiles of Nu are fluctuating under the effects of pertinent parameters due to the dual rotation between an inner fin and outer hexagonal-shaped domain at the transition state. It is seen that the profiles of Nu are affected by the variations on the pertinent parameters. It is seen that a fractional derivative α = 1 gives the highest values of Nu . Second, the value of θ f = 0.2 provides the highest values of Nu and the tendency of Nu is fluctuating below the variation of θ f . Third, the values of

Conclusion
The originality of the study is emulating the natural convection of NEPCM embedded in a hexagonal-shaped cavity under the impacts of a magnetic field and dual rotation between an inner fin and outer hexagonal-shaped domain. The ISPH method is developed by including the time-fractional derivative based on the Grunwald-Letnikove derivative in the solving steps for conducting the current physical problem. The executed simulations indicated that the variations on the fractional time derivative are changing the dual rotation amongst the inner fin and outer domain. As a result, the nanofluid movements, heat transfer, and phase change material within a hexagonal-shaped cavity are affected by the variations of a fractional time derivative. The inserted cool fin is

Future work
In future work, the fractional-time derivative will adopt the most recent formulation in the fractional calculus. More respective studies in the fractional-space derivative will be researched.