RETRACTED ARTICLE: Hydrothermal analysis on MHD squeezing nanofluid flow in parallel plates by analytical method

In this paper, the heat and mass transfer of MHD nanofluid squeezing flow between two parallel plates are investigated. In squeezing flows, a material is compressed between two parallel plates and then squeezed out radially. The significance of this study is the hydrothermal investigation of MHD nanofluid during squeezing flow. The affecting parameters on the flow and heat transfer are Brownian motion, Thermophoresis parameter, Squeezing parameter and the magnetic field. By applying the proper similarity parameters, the governing equations of the problem are converted to nondimensional forms and are solved analytically using the Homotopy Perturbation Method (HPM) and the Collocation Method (CM). Moreover, the analytical solution is compared with numerical Finite Element Method (FEM) and a good agreement is obtained. The results indicated that increasing the Brownian motion parameter causes an increase in the temperature profile, while an inverse treatment is observed for the concentration profile. Also, it was found that enhancing the thermophoresis parameter results in decreasing the temperature profile and augmenting the concentration profile. Effects of active parameters have been considered for the flow, heat and mass transfer. The results indicated that temperature boundary layer thickness will increases by augmentation of Brownian motion parameter and Thermophoresis parameter, while it decreases by raising the other active parameters.


Background
Investigation of heat and mass transfer of viscous flow between two parallel plates is one of the most important and well-known academic research topics because of the wide range of its scientific and engineering approaches, such as polymer processing, systems of lubricating, cooling towers, and food processing.
Utilizing nano-scale particles in the base fluid seems to be a creative technique to improve the heat transfer rate. As known, nano-particle-containing fluids are called nanofluids. Choi 1995 was the first one ever to name the nano-particle-containing fluids as nanofluids. Sheikholeslami and Ganji (2013a) studied analytically the heat transfer of a nanofluid flow compressed between parallel plates using Homotopy Perturbation Method (HPM). They indicated that the Nusselt number is directly related to nanoparticle volume fraction along with Squeeze number and Eckert number for two separated plates, while there is an inverse relationship between the Nusselt number and Squeeze number when two plates are squeezed. Sheikholeslami and Bhatti (2017a) studied the heat transfer enhancement of nanofluid flow by using EHD. Results indicated that the effect of Coulomb force is more considerable in lower values of Reynolds number. The shape effects of nanoparticles on the natural convection nanofluid flow in a porous semi-annulus were studied by Sheikholeslami and Bhatti (2017b). Results illustrated that the maximum rate of

R E T R A C T E D A R T I C L E
heat transfer is obtained at platelet shape. Bhatti and Rashidi (2016) investigated the influences of thermos-diffusion and thermal radiation on nanofluid flow over a stretching sheet. The authors showed that the temperature profile is an increasing function of thermal radiation and thermophoresis parameters.  reviewed the micropolar nanofluid flow over a porous stretching sheet. It was found that raising the radiation parameter enhances the boundary layer thickness.  analyzed the stagnation-point flow of hybrid nanofluid over a stretching sheet. They concluded that using hybrid nanofluid instead of conventional nanofluid results in higher Nusselt numbers. Dogonchi et al. (2017) investigated the influence of thermal radiation on MHD nanofluid flow in a porous channel. Results illustrated that there is a direct relationship between the Nusselt number and nanofluid volume fraction. Recently, many authors analyzed the effect of nanoparticles in theirs studies Bhatti et al. 2017b;Saedi Ardahaie et al. 2018;.
The results of the time-dependent chemical reaction on the viscous fluid flow over an unsteady stretching sheet were considered by Abd-El Aziz (2010). Furthermore, for a certain viscous fluid between parallel disks, the magneto hydrodynamic squeezed flow was investigated by Domairry and Aziz (2009). Also, for the fluid flow between parallel plates, the Homotopy Analysis Method (HAM (Domairry and Ziabakhsh 2009a;Domairry and Ziabakhsh 2009b;Ziabakhsh et al. 2009)) was utilized by Mustafa et al. (2012).
The majority of the problems in the field of engineering especially heat transfer equations contain nonlinear equations. In this case, some of these nonlinear equations can be solved by numerical approaches, while some others are solvable using various analytical methods such as Perturbation method (PM) (Bhatti and Lu 2017), Collocation method (CM) (Rahimi et al. 2017;Atouei et al. 2015), Homotopy perturbation method (HPM), and Variational iteration method (VIM). Hence, the elimination of small parameter has been the critical issue for the scientists nowadays which has led to introduction of various ways to solve these special problems. One of these ways is the implementation of the semi-exact method called "HPM" which does not need any small parameters. The Homotopy perturbation method was suggested and modified by He (2004). This method results in a rapid convergence of the solution series in most cases. Including both the efficiency and the accuracy in solving a large number of nonlinear problems, the HPM proved itself capable in dealing with such problems. Ijaz et al. (2018) applied the HPM to investigate the effect of liquid-solid particles interaction in a wavy channel. Dogonchi et al. (2015) analyzed the sedimentation of non-spherical particles in Newtonain media using DTM-Pade approximation. It was found that enhancing the sphericity of particles results in augmenting the velocity profile. Mosayebidorcheh et al. (2016) studied the analysis of turbulent MHD Couette nanofluid flow and heat transfer using hybrid DTM-FDM. Sheikholeslami et al. (2011) investigated the rotation of MHD viscous flow along with the heat transfer between stretching and porous surfaces using HPM. Results showed that an increase in the rotation parameter along with the increase in the blowing velocity parameter and Prandtl number would cause an increase in the Nusselt number. The profiles of the variables such as the velocity, temperature, and the concentration of the nanofluids affected by the magnetic field are investigated by Uddin et al. (2014). They comprehended that the presence of magnetic field would cause a decrease in the velocity field and an increase in temperature and concentration profiles. Also, it was found that the convective heating parameter leads to augment the velocity, temperature, and concentration profiles. Uddin et al. (2014) also studied non-Newtonian nanofluid slip flow over a permeable stretching sheet. Their result indicated that the skin friction factor plays a key role in the characteristics of nanofluid flow. In addition, the chemical reaction of nanofluid in free convection in the presence of magnetic field was investigated by Uddin et al. (2015). During the study of Jing et al. (2015), they discovered that the presence of nanoparticles within the fluid can extremely increase the effective thermal conductivity of the fluid, and as a result, the heat transfer characteristics will be improved. Sheikholeslami and Ganji (2013b) examined the nanofluid flow squeezed between parallel plates utilizing Homotopy perturbation method (HPM). They reported that the Nusselt number has direct relationship with nanoparticle volume fraction, the Squeeze number and the Eckert number in the case of separated plates, while its relationship with the Squeeze number in the case of squeezed plates is vice versa. Paying attention to the nanoparticle migration, the mixed convection of alumina-water nanofluid inside a concentric annulus was investigated by Malvandi and Ganji (2016). Sulochana et al. (2016) examined the effect of transpiration on the magnetohydrodynamic stagnation-point flow of a Carreau nanofluid toward a stretching/shrinking sheet in the presence of thermophoresis and Brownian motion, numerically. They discovered that by raising the thermophoresis parameter, both the heat and mass transfer rates will be increased, whereas the Weissenberg number enlarges the momentum boundary layer thickness along with the heat and mass transfer rate. Sheikholeslami et al. (2016) studied the effect of Lorentz forces on forced-convection nanofluid flow over a stretched surface. Their results indicated that the skin friction coefficient increases by amplifying the magnetic field, while it decreases by enhancing the velocity ratio parameter. Sudarsana Reddy and Chamkha (2016) analyzed the influence of size, shape, and type of nanoparticles along with the type and temperature of the base fluid on the natural convection MHD nanofluid flow. Their results revealed that decreasing the size of the nanoparticles leads to a significant natural convection heat transfer rate. Moreover, types of nanoparticles and the base fluid also impressed the natural convection heat transfer. Mishra and Bhatti (2017) investigated the MHD stagnation-point flow over a shrinking sheet, numerically. The authors compared the accuracy of their solution with previous studies and found that a good agreement was obtained. Newly, the study of MHD flow in different geometries has attracted many attentions Hatami et al. 2014;Bhatti et al. 2018).
The main goal of the present study is to investigate the effect of Brownian motion and thermophoresis phenomenon on the squeezing nanofluid flow and heat transfer between two parallel flat plates in the presence of variable magnetic field. Both the flow and heat transfer characteristics have been examined under the effects of Squeeze number, suction parameter, Hartmann number, Brownian motion parameter, Thermophoretic parameter, and Lewis number.

Problem description and governing equations
This study is concerned with incompressible twodimensional flow of squeezing nanofluid between two parallel and movable plates at distance of h(t) = H(1 − at) 1/2 from each other. The schematic model of the problem is depicted in Fig. 1. As shown in Fig. 1, B(t) = B 0 (1 − at) −1/2 is the variable magnetic field that applied perpendicular to the plates. To simplify the problem, only the flow patterns on the left part of the channel have been mentioned. It should be noted that the flow patterns of squeezing flow are axisymmetric.
The "×" marks show that magnetic field is perpendicular to the illustrated plane. "+" and "−" indicate the positive and negative charges, respectively. T w and C w represent the temperature and concentration of nanoparticles at the bottom disk, while the temperature and concentration of nanoparticles at the upper disk are denoted by T H and C H . The upper disk at Z = h(t) can move toward or away from the motionless bottom disk with the velocity of dh . dt .
For a > 0 and a < 0, two plates are squeezed and separated, respectively. The viscous dissipation effect along with the generated heat remained intact due to the friction caused by shear forces in the flow. It should be noted that when the fluid is largely viscous or is flowing at a high speed, the dissipation effect is quite important. Knowing that the nanofluid is a two-component mixture, the following assumption has been considered: Incompressible; no-chemical reaction; with negligibleness of viscous dissipation and radiative heat transfer; nano-solid-particles and the base fluid are in thermal equilibrium and without any slip between them. The equations which govern the flow, heat, and mass transfer in viscous fluid are as follows (Hashmi et al. 2012;Turkyilmazoglu 2017): where is the velocity vector; T, P, ρ, μ, C P, K are the temperature, pressure, density, viscosity, heat capacitance, and thermal conductivity of nanofluid, respectively. Also, the operation of ∇ ! can be defined as: The boundary conditions are as follows: where u and v represent velocity components in the rand z-directions, respectively, ρ is the density, μ is the dynamic viscosity, p is the pressure, T is the temperature, C is the nanoparticles concentration, α is the thermal diffusivity, D B is the Brownian motion coefficient, T m is the mean fluid temperature, and k is the thermal conductivity. The total diffusion mass flux for nanoparticles is the last term in the energy equation which is given as a sum of the Brownian motion and thermophoresis terms. In addition, τ is the dimensionless parameter which can calculate the ratio of effective heat capacity of the nanoparticles to heat capacity of the fluid. The parameters of similarity solution are as follows: By removing the pressure gradient from Eqs. (2) and (3), then rewriting Eqs. (4) and (5), the final nonlinear equations can be obtained as follows (Turkyilmazoglu 2016): Boundary conditions are described as follows: where S is the Squeeze number, A is the suction/blowing parameter, M is the Hartmann number, Nb is the Brownian motion parameter, Nt is the Thermophoretic parameter, Pr is the Prandtl number, and Le is the Lewis number and are defined as follows: The continuity equation is identically satisfied. It should be noted that A > 0 indicates the suction of fluid from the lower disk, while A < 0 represents the injection flow.

Collocation method (CM)
Suppose we have a differential operator D acting on a function u to produce a function p (Hatami et al. 2013).
Function u can be considered as a functionũ, which is a linear combination of basic functions chosen from a linearly independent as follows: Now, we can substituteũ from Eq. (12) into the Eq. (11), generally p (x) is not the result of the operations. Therefore an error or residual will exist: The basic principle of the Collocation method is to lead an error or the residual to zero in some average sense over the domain as follows: So that the number of weight functions W i and the number of unknown constants c i (Eq.(13)) are exactly equal. The result is a set of n algebraic equations for the unknown constants c i . In collocation method, the weighting functions are obtained from the family of Dirac δ functions in the domain. That is, W i (x) = δ(x − x i ). The Dirac function is defined as follows:

Application of CM
To obtain an approximate solution for Eq. (8) in the domain 0 < η < 1, we consider the basic function to polynomial in η. The trial solution contains three undetermined coefficients and satisfies the condition for all values of c as follows: where Eq. (16) satisfies the boundary conditions of Eq. (9). The residual function (R (c1, c2, c3, η)) can be obtained by substituting Eq. (16) into Eq. (13). The residual is equal to zero only by exact solution of the problem. Here, the problem is solved by the approximate solution so that the residual stays close to zero throughout the domain 0 < η < 1. Three points are needed to find the three unknown parameters, so three specific points with approximately equal distance should be chosen in the domain. These points are: Eventually, by substitutions values of Eq. (17) into residual function R (c1, c2, c3, η), a set of three long equations with three unknown coefficients are obtained. After solving these unknown parameters (c1, c2, c3), the temperature distribution equation, Eq. (16), will be determined.

Results and discussion
In the present study, the effect of Brownian motion and Thermophoresis phenomenon on the heat and mass transfer of MHD nanofluid flow between parallel plates is investigated, and Collocation Method (CM), Homotopy Perturbation Method (HPM) along with the finite element Method (FEM) are applied to solve this problem using Maple 16 and FlexPDE 5 softwares. The influence of certain active parameters such as Squeeze number, suction parameter, Hartmann number, Prandtl number, Brownian motion parameter, Thermophoretic parameter and Lewis number on the flow and heat transfer characteristics are examined. The presented code is validated by comparing the obtained results with the results of finite element method (FEM) (Fig. 2). The comparison well showed that by implementing this code, a highly accurate solution is obtained to solve the problem.

L E
The effect of suction parameter, Squeeze number, and Hartmann number on velocity profile is shown in Fig. 3. Increasing the suction parameter would cause an increase in velocity profile due to the increase of the turbulence in the flow. It can be seen that the velocity values drop by enhancing the Squeeze number because the plate remained close to each other and limits the velocity. Also, it can be found that enhancing the Hartmann number in the flow results in augmenting the velocity profile.
Figures 4 and 5 represent the influences of suction parameter and Squeeze number on the temperature and concentration profiles, respectively. Figures depicted that increasing the suction parameter would cause a decrease in thermal boundary layer thickness and concentration profiles. Effect of Brownian motion parameter on temperature and concentration profiles is shown in Fig. 6, while the effect of Thermophoretic parameter on the mentioned profiles is examined in Fig. 7. It can be observed that increasing the Brownian motion parameter results in increasing the temperature profile, while the influence of Thermophoretic parameter on temperature profile is vice versa compared to Brownian motion parameter, whereas increasing both the Brownian motion parameter and Thermophoretic parameter individually would cause a decrease in concentration profiles which is depicted in Fig. 8. As shown in these figures, the nanoparticle temperature and concentration profiles are decreasing functions of the Hartmann number. Figure 9 shows the effect of Lewis number on temperature and concentration profiles. It is observed that an increase in temperature profile near the bottom plate and also a decrease in temperature profile near the top plate are the results of enhancing the Lewis number. Finally, it should be mentioned that higher values of nanoparticle concentration are obtained by enhancing the Lewis number.

Conclusion
The present study examines the effect of Brownian motion and Thermophoresis phenomenon on the heat and mass transfer of MHD nanofluid flow between parallel plates. To examine this problem, a number of methods such as the Collocation Method (CM), the Homotopy Perturbation Method (HPM), and the Finite Element Method (FEM) were applied. The results indicated that the outcomes of Collocation Method have the best agreement with the numerical solutions. The crucial effect of Brownian motion and thermophoresis parameter has been included in the model of nanofluid. Effects of active parameters have been considered for the flow, heat, and mass transfer. The results indicated that temperature boundary layer thickness will increase by augmentation of Brownian motion parameter and Thermophoresis parameter, while it decreases by raising the other active parameters. Also, it can be concluded that the thickness of concentration boundary layer declines by enhancing the Brownian motion parameter, while an inverse trend is observed by augmenting the Thermophoresis parameter.