Non-equilibrium Model for Nanofluid Free Convection Inside a Porous Cavity Considering Lorentz Forces

In current article, transportation of CuO nanoparticles through a porous enclosure is demonstrated. The enclosure has complex shaped hot wall. Porous media has been simulated via two temperature equations. Magnetic force impact on nanofluid treatment was considered. Control volume based finite element method has been described to solve current article in vorticity stream function form. Single phase model was chosen for nanofluid. Nanofluid characteristics are predicted via KKL model. Roles of solid-nanofluid interface heat transfer parameter (Nhs), porosity, Hartmann and Rayleigh numbers have been illustrated. Outputs illustrated that conduction mode reduces with augment of Ra. Increasing magnetic forces make nanofluid motion to decrease. Temperature gradient of nanofluid decreases with augment of Nhs. Reducing porosity leads to enhance in Nusselt number.


Non-equilibrium Model for Nanofluid
Solar power collectors and drying technologies are two common uses for porous enclosure. In some application, researchers should consider two-temperature model. Alsabery et al. 1 demonstrated nanoparticletransportation in a tilted porous cavity. They indicated that convective flow is significantly influenced by the permeable layer augmentation. Lu et al. 2 investigated about nanofluid radiative three dimensional flow containing gyrotactic microorganism with anisotropic slip. They considered the effect of activation energy. Zaimi et al. 3 illustrated nanofluid boundary layer movement on a porous plate. Sheikholeslami and Shehzad 4 showed the two temperature model for nanoparticle migration inside a permeable medium. They revealed that Nu increases with decrease of porosity. Khan et al. 5 reported the nanofluid mixed convection over an oscillating vertical plate. They utilized Laplace transform method to solve the governing equations.
Sheikholeslami and Shehzad 6 investigated the role of radiation on nanoparticle treatment. They found that Nu decreases with reduce of radiation parameter. Haq et al. 7 used carbon nanotubes to improve convective heat transfer over plate with slip flow. Carbon Nanotubes has been dispersed in to engine oil by Haq et al. 8 is examined of magnetic forces. Tripathi et al. 9 illustrated viscous dissipation and Hall effects on nanofluid rotating flow. Selimefendigil and Oztop 10 depicted impact of inclination on hydrothermal behavior. They found that titled angle can be used as control parameter.
Promvonge et al. 11 applied new way in a duct to improve the thermal characteristics. Sheikholeslami 12 described the impact of electric filed on nanofluid free convection. He proved that Nusselt number enhances by adding electric field. Aman et al. 13 illustrated the nanofluid thermal improvement in migration of CNTs nanoparticles. Sheikholeslami and Seyednezhad 14 illustrated nanofluid Electrohydrodynamic flow in a permeable enclosure. Different applications of Fe3O4-water nanofluid were categorized by Sheikholeslami and Rokni 15 . Akbar et al. 16 showed the role of Hartmann flow on nanoparticles migration in a duct. Najib et al. 17 demonstrated the impact of chemical factor on flow style. They found that the Nu augments with augment of curvature. Various articles were been available about nanoparticle migration through porous media [18][19][20][21] .
Current publication is about nanoparticle migration in a porous enclosure with two temperature model via CVFEM considering magnetic force. Results illustrate the roles of significant parameters on contours.

CVFEM and Explanation
Formulation. According to existence of magnetic force and two temperature model for porous medium the basic formulas are:         (ρC p ) nf , (ρβ) nf , ρ nf , σ nf and k nf , μ nf can define as 22 : p nf p p p f    where φ is nanofluid volume fraction. Required characteristics and parameters are illustrated in Tables 1 and 2 22 .
Considering following definitions:       Nu loc and Nu ave are: Nu dr 1 2 (20) ave l oc 0 2 CVFEM. The innovative in which triangular element is used and upwind method is applied for advection term is CVFEM (Fig. 1(b)). Gauss-Seidel is the name of the method which is used for final step as mentioned in ref. 23 .

Mesh Independent Test and Validation
Obviously, the final outputs should not alter by changing mesh size. So, this test should be done for various cases as illustrated Table 3. Also, we should be sure about accuracy of written code by applying this code for previous published problem. Table 4 and Fig. 2 illustrates nice accuracy [24][25][26] .

Results and Discussion
In current article, migration of nanoparticle in a porous medium which is described by new porous model is presented. Numerical method is applied to display the roles of the porosity (ε = 0.3 to 0.9), Rayleigh number (Ra = 100,500 and 10 3 ), Hartmann number (Ha = 0 to 20) and Nhs (Nhs = 10 to 1000).
Influences of Nhs, Ra, ε and Ha on isotherms for solid (θ s ), streamlines (Ψ) and nanofluid (θ nf ) were depicted in Figs 3, 4, 5, 6, 7 and 8. When buoyancy force is weak, conduction mode is significant and the impacts of other variables are negligible. In this case, (θ nf ) contours are similar to (θ s ) contours. As Ra increase, thermal plume can be seen close to the vertical centerline and (θ nf ) contours become more complex but (θ s ) contours have no changes. Velocity of nanofluid decrease with augment of magnetic forces and the (θ nf ) contours are stratified with enhance of Ha. (Ψ max ) enhances with augment of Nhs due to stronger convective flow. As ε enhances, the pores volume through the enclosure augments. Therefore, convective flow becomes stronger. Besides, the impact of ε on isotherms is as same as Nhs. Figures 9 and 10  Coefficient of determination is 0.98 for this correlation. Also, in this formula we have: Ha* = 0.1Ha, Nhs* = 0.001Nhs, Ra* = 0.001Ra. Nu ave is a reducing function of Nhs because temperature gradient decreases with augment of this factor. Decreasing porosity makes Nu to enhance which is similar of Ha impact. Nu ave decreases with reduce of Rayleigh number.

Conclusions
In current research, migration of CuO nanoparticles is simulated via Non-equilibrium. Innovative method is applied to show the impacts of porosity, buoyancy, magnetic forces and the Nhs. Results show that Nu ave reduces with augment of Ha, ε, Nhs. Convective flow enhances with increase of Nhs and ε while it reduces with enhance of magnetic force. When Ra = 1000, ε = 0.3, increasing Nhs leads to 17.06 percent decrement of Nu ave in absence of magnetic field. Impact of Nhs is negligible in existence of Lorenz force. Also Nu ave decreases about 53% with augment of Hartmann number.