CuO–Water Nanofluid Magnetohydrodynamic Natural Convection inside a Sinusoidal Annulus in Presence of Melting Heat Transfer

Impact of nanofluid natural convection due tomagnetic field in existence of melting heat transfer is simulated using CVFEM in this research. KKLmodel is taken into account to obtain properties of CuO–H2Onanofluid. Roles of melting parameter (δ), CuO–H2O volume fraction (φ), Hartmann number (Ha), and Rayleigh (Ra) number are depicted in outputs. Results depict that temperature gradient improves with rise of Rayleigh number and melting parameter. Nusselt number detracts with rise of Ha. At the end, a comparison as a limiting case of the considered problem with the existing studies is made and found in good agreement.


Introduction
Melting process has various uses such as application like in thermocouple, heat exchangers, and heat engines.Chamkha and Ismael [1] examined the Cu-H2O nanofluid mixed convection in an enclosure with lid wall.Also, Ellahi et al. [2] reported the similar study for nanofluid over permeable wedge in mixed convection.Harikrishnan et al. [3] investigated melting behavior of composite PCM for heating application.Nanofluid melting in a pipe with external fins has been studied by Mat et al. [4].Manikandan and Rajan [5] utilized hybrid nanofluid for heat transfer augmentation.Koca et al. [6] utilized Ag-H2O nanofluid for natural circulation loops.Shyam and Tiwari [7] investigated the coiled heat exchanger by means of nanofluid.Das [8] reported the radiative flow in existence of melting heat.
Sheremet et al. [9] illustrated the free convective flow of ferrofluid in a rotating enclosure.Wavy duct in existence of Brownian forces has been examined by Shehzad et al. [10].Sheikholeslami et al. [11] illustrated different uses of nanotechnology in their article paper.Abbas et al. [12] demonstrated the nanofluid flow through a horizontal Riga plate.Chamkha and Rashad [13] reported the nanoparticle migration on porous cone.Ellahi et al. [14] analyzed particle shape effects on Marangoni convection boundary layer flow of a nanofluid.They considered the Lorentz forces impact in governing equations.Malvandi et al. [15] reported nanofluid flow inside a channel in existence of Lorentz forces.Garoosi et al. [16] investigated performance of heat exchanger via nanofluid.They indicated that optimum volume fraction of nanoparticle exists to obtain maximum Nusselt number.Mesoscopic method has been utilized by Sheikholeslami and Ellahi [17] for a three-dimensional problem.
This research intends to present the impact of melting heat transfer on free convection of ferrofluid in the presence of Lorentz forces.CVFEM is selected to find the outputs.Roles of melting parameter, CuO-water volume fraction, and Hartmann and Rayleigh numbers are presented.the outer one is melting surface ( =   ).Other walls are adiabatic.Horizontal magnetic field has been applied.The enclosure is field with nanofluid.

Governing Equation and Simulation
3.1.Governing Formulation.2D steady convective nanofluid flow is considered in existence of constant magnetic field.The fluid flow is laminar and incompressible.The flow is steady and Newtonian.The viscous dissipation is negligible in this study.The effects of Brownian force and thermophoretic force are not taken into condition.The flow also follows the Boussinesq assumption.The prevailing equations under these constraints can be written as [20] V  +   = 0, (  ) nf ,  nf , () nf , and  nf are calculated as [17,18] nf ,  nf are calculated via KKL model [19]: Properties and needed parameters are provided in Tables 1  and 2 [19].Vorticity and stream function should be used to eliminate pressure source terms: Dimensionless quantities are introduced as follows: The final formulae are Boundary conditions are   = 0.0 on other walls,  = 0 on outer wall,  = 1.0 on inner wall, Ψ = 0.0 on all walls except melting surface (7) and in melting surface, we have where dimensionless and constants parameters are illustrated as It should be mentioned that  is related to Stefan numbers.
Local and average Nusselt over the cold wall can be calculated as follows:  3.2.Numerical Procedure.CVFEM uses both benefits of two common CFD methods.This method uses triangular element (see Figure 1(b)).Upwind approach is utilized for advection term.Gauss-Seidel approach is applied to solve the algebraic equations.Further notes can be found in [21].

Grid Independency and Validation
Outputs should not rely on mesh size.Therefore, several grids should be tested.For example, as shown in Table 3, a grid size of 71 × 211 can be selected.The correctness of CVFEM code for nanofluid natural convective heat transfer is demonstrated in Figure 2 ( [18]).This CVFEM code has good accuracy for magnetohydrodynamic flow as depicted in Table 4 [20].

Results and Discussion
Nanofluid flow in a half sinusoidal annulus due to magnetic field in presence of melting surface is examined. nf ,  nf of CuO-water nanofluid are estimated by means of KKL model.Graphs and tables are depicted for different amounts of CuO-H 2 O volume fraction ( = 0 to 0.04), melting parameter ( = 0 to 0.2), Rayleigh number (Ra = 500 to 5000), and Hartmann number (Ha = 0 to 40).Impact of adding CuO nanoparticles in water on velocity and temperature contours is depicted in Figure 3. Temperature gradient decreases with augment of .|Ψ max | augments with adding nanoparticles because of increment in the solid movements.In presence of melting heat transfer and magnetic field, effect of adding nanoparticles on isotherms becomes negligible.
Figures 4 and 5 depict the impact of Rayleigh and Hartmann numbers in absence of melting heat transfer.There  is only one eddy in streamline.In low Rayleigh number, conduction mechanism is dominant.As Ra increases, the distortion of isotherms is enhanced close to the hot wall.Adding magnetic field makes isotherms become parallel.In presence of melting heat transfer the primary eddy diminishes and increasing Lorentz forces generates three layers for streamline.Increasing melting parameter augments the bottom eddy.
where Ha * = 0.1Ha and Ra * = 0.001Ra.As melting parameter augments, temperature gradient is enhanced and in turn Nusselt number is enhanced.Increasing buoyancy forces leads the thermal boundary layer thickness to reduce.So Nusselt number increases with enhancement of Ra.As Hartmann number augments, isotherms become parallel.Therefore Nusselt number has opposite relationship with Ha.

Conclusions
Nanofluid free convection due to Lorentz force in existence of melting surface is reported.Combination of FEM and FVM is utilized to solve the PDEs.KKL model is considered for nanofluid.Roles of Hartmann number, CuO-water volume fraction, Rayleigh number, and melting parameter are presented.Outputs depict that temperature gradient improves with augment of melting parameter and Rayleigh number.Adding magnetic field makes the temperature gradient reduce due to domination of conduction mechanism in high Hartmann number.

Figure 1 Figure 1 :
Figure 1 depicts the geometry, boundary condition, and sample element.The inner wall is hot wall ( =  ℎ ) and

Figure 6
Figure 6 illustrates the impact of , Ra, and Ha on Nu ave .The formula for Nu ave corresponding to active parameters is Nu ave = 0.799 + 0.305 + 0.19Ra * − 0.16Ha *

Figure 6 :
Figure 6: Impacts of , Ha, and Ra on Nu ave .

Table 4 :
Nu ave for various Gr and Ha at Pr = 0.733.