Approximation of First Grade MHD Squeezing Fluid Flow with Slip Boundary Condition Using DTM and OHAM

A steady two-dimensional axisymmetric flow of an incompressible viscous fluid under the influence of a uniform transverse magnetic field with slip boundary condition is studied. An ordinary nonlinear differential equation is formed by transforming the Navier-Stokes equations using the transformation ψ(r, z) = rG(z). Differential transform and optimal homotopy analysis methods have been used to obtain the solutions by varying pertinent flow parameters. By using residuals in each case, the validity of solutions is established. Excellent results are obtained by using the proposed schemes. The influence of different parameters on the flow is shown through graphs.


Introduction
Squeezing flows induce when normal stresses or vertical velocities are applied externally by means of a mobile boundary [1].Stefan [2] carried out the fundamental research in this field.Engineers studied the analysis of Newtonian fluid squeezed between two infinite planar plates [3,4].There are many applications of squeezing flow in food industry, especially in chemical engineering [4].Polymer processing, compression, and injection modeling are the practical examples of squeezing flow.Thin Newtonian liquid films squeezing between two plates were studied by Grimm [5].Squeezing flow under the influence of magnetic field is widely applied to bearing with liquid-metal lubrication [6][7][8][9].Islam et al. [10] studied squeezing fluid flow between the two infinite parallel plates in a porous medium channel.
In case of many polymeric liquids when the weight of molecule is high, then they show slip at the boundary.The no-slip boundary condition is not applicable in this case.In many cases such as thin film problems, rarefied fluid problems, fluids containing concentrated suspensions, and flow on multiple interfaces, the no-slip boundary condition fails to work.Navier [11], for the first time, proposed the general boundary condition which demonstrates the fluid slip at the surface.The difference of fluid velocity and velocity of the boundary is proportional to the shear stress at that boundary.The proportionality constant is named the slip parameter having length as its dimension.
The slip condition is of great importance especially when fluids with elastic character are under consideration [12].Newtonian fluid was considered by Ebaid [13] to study the effects of magnetic field and wall slip conditions on the peristaltic transport in an asymmetric channel.It has great importance in medical sciences, particularly in polishing artificial heart valves and internal cavities in many manufactured parts achieved by embedding such fluids with abrasives [14].The influence of slip on the peristaltic motion of third-order fluid in asymmetric channel is studied by Hayat et al. [15].The effects of slip condition on the rotating flow of a third grade fluid in a nonporous medium are investigated by Hayat and Abelman [16].The extension of the work in [16] to a porous medium and obtaining the numerical solutions for the steady magnetohydrodynamics flow of a third grade fluid in a rotating frame is presented by Abelman et al. [17].
Various perturbation techniques can be applied to solve nonlinear differential equations [18].The limitations of these methods are based on the assumption of small parameter and there is no proper way to select this parameter.The idea of homotopy was combined with perturbation in the last decade.Liao and He did the fundamental work on it.Homotopy analysis method (HAM) was proposed by Liao [19][20][21].Homotopy perturbation method was introduced by He [22][23][24].Marinca et al. [25,26] introduced OHAM for approximate solution of nonlinear problems of thin film flow of a fourth-grade fluid down a vertical cylinder and for the study of the behavior of nonlinear mechanical vibration of electrical machines.It is observed that HPM and HAM are the special cases of OHAM [27,28].Differential transform method (DTM) was initially introduced by Zhou in 1986 [29].C ¸atal [30][31][32][33] studied solution of free vibration equations of beam on elastic soil, buckling analysis of partially embedded pile in elastic soil, analysis of free vibration of beam on elastic soil, and response of forced Euler-Bernoulli beams by using DTM.Ayaz [34] applied this method to differential algebraic equations (DAEs) of index-1.Liu and Song [35] investigated that DTM is effective in case of index-2 DAEs but not in index-3.Ayaz [36] studied the applications of two-dimensional DTM in case of partial differential equations.

Mathematical Problems in Engineering
The objective of this research paper is to use OHAM and DTM for studying MHD squeezing flow with slip boundary condition between two infinite plates approaching each other slowly.In Section 2, the basic equations are derived and the model is reduced to a nonlinear boundary value problem.In Sections 3 and 4, the basic ideas of DTM and its use for our problem are presented while the idea of OHAM and its application to our problem are in Sections 5 and 6, respectively.

Problem Formulation
Let us consider, in the presence of a magnetic field, a squeezing flow of an incompressible Newtonian fluid with constant density  and viscosity , squeezed between two large planar parallel plates separated by a small distance 2ℎ approaching each other with a low constant velocity k (Figure 1).Assume that the flow is quasisteady [1,3], and the Navier-Stokes equations governing such flow when inertial terms are retained are where ∇ denotes the material time derivative, T is the Cauchy stress tensor given by T = − +  with  = ∇k + (k)  .B is the total magnetic field given by B = B 0 +.B 0 and  represent the imposed and induced magnetic fields, respectively.The modified Ohm's law and Maxwell's equations (see [37] and the references therein), in the absence of displacement currents, are ( Here J is the electric current density,  represents the electrical conductivity,  the electric field, and   the magnetic permeability.If ,   , and  are constant,  is negligible as compared to B 0 , B is perpendicular to k so that the Reynold number is small with no electric field in the fluid flow region and then the magneto hydrodynamic force involved can be written as Assuming that the plates are nonconducting and the magnetic field is applied along the -axis.The gap distance 2ℎ between the plates changes slowly with time  for small values of the velocity k so that it can be taken as constant.An axisymmetric flow in cylindrical coordinates (, , ) with axis perpendicular to plates and  = ±ℎ at the plates.For axial symmetry, V is represented by V = (V  , 0, V  ).In view of negligible body forces with no tangential velocity, Navier-Stokes equations [1,6,10] in cylindrical coordinates are where Introducing the stream function (, ), we have, Eliminating  from (4), we have where Using the transformation (, ) =  2 (), (7) can be written as subject to the slip boundary conditions The nondimensional parameters are  * = /V/2,  * = /ℎ,  = ℎ//V, and  = B 0 ℎ√/.Omitting the * , ( 9) and ( 10) become with  = /ℎ and ,  are the Reynolds and Hartmann numbers, respectively.

Differential Transform Method
One-dimensional differential transform of a function () is defined as follows [29][30][31][32][33][34][35]: () is also called the -function.The inverse transform of () is is defined as Combining ( 12) and ( 13), we can write Usually we express () in a finite series, that is, which implies that is negligibly small.Some fundamental theorems on one-dimensional differential transform are as follows.

Analysis of Differential Transform Method
For fourth-order boundary value problem [10], with boundary conditions The differential transform is given by with transformed boundary conditions Using ( 20) and ( 21), () for  = 4, 5, . . . is evaluated to obtain the solution up to ( +1 ) as follows:

Application of DTM to Our Problem
Here, we consider an incompressible Newtonian fluid, squeezed between two large planar, parallel plates which are separated by a small distance 2ℎ and moving towards each other with velocity k.For small values of k, the gap distance 2ℎ between the plates changes slowly with time , so that it can be taken as constant, and the flow is steady as with slip boundary conditions The differential transform of ( 23) is given by Mathematical Problems in Engineering with transformed boundary conditions where  and  are unknowns to be determined later [38].Using ( 25) and ( 26), the values of G(),  = 1, 2, 3, . . ., 15 are G (2) = 0, for  = 0 (1) 7, As ,  have been nondimensionalized, so let us fix them one to find the values of  and  from the following system: which leads us to the following values: Using (22), the approximate solution of the BVP is as follows: This is a series solution and we note that the terms approach to zero if we go forward which is the beauty of DTM and other series solutions.
We choose the auxiliary function ℎ() such that where   are the convergence controlling constants to be determined.Expanding (, ) in Taylor's series about  to get the approximation substituting ( 36) into (33), and equating the coefficients of like powers of , we get the following.
Zeroth-order problem is First-order problem is Second-order problem is Generally, with boundary conditions where   ( 0 (),  1 (), . . .,  −1 ()) is the coefficient of   in the expansion of ((, )) about  as The convergence of the series (36) G will be the exact solution if  = 0, but in case of nonlinearity, it does not happen generally.To find the optimal values of the constants   , there are different methods to apply.We apply the method of least square as follows: Minimizing this function, we have where  and  are chosen in the domain of the problem for locating the suitable    ( = 1, 2, . . ., ).The approximate solution is well determined for these values of   .

Application of OHAM to Our Problem
Using  = 1,  = 1, and  = 1 in (11), we have the following.Zeroth-order problem is First-order problem is Second-order problem is (49)  Third-order problem is Fourth-order problem is Fifth-order problem is By considering the fifth-order solution, we have The residual of the problem is In order to find  1 ,  2 , we apply the method of least square as follows: Solving (56) for  1 ,  2 , we get Using these values of  1 ,  2 , the approximate solution is given by: As  ∈ [0, 1], the components of () vanish when going onward.In Table 1, the solutions obtained by OHAM and DTM along with their residual are presented.Table 2 is constructed to discuss residuals for different values of .

Numerical Results and Discussion
Here we analyzed, under a transverse magnetic field, the squeezing fluid flow between two parallel infinite plates.DTM and OHAM are applied to obtain the solution of reduced nonlinear boundary value problem.The convergence region can easily be adjusted using OHAM as compared to DTM in this case.Differential transform method which is the straight forward application of Taylor's series suffers in general and particularly in initial value problems.Beside this fact, DTM leads easily to closed form solution.OHAM does    not require discretization like other numerical and approximate methods.This technique requires less computational effort and the convergence of this method to the exact is fast enough.For computational purposes, we use Mathematica software.
In order to get an insight into the problem, it is important to discuss the influence of different flow parameters on the velocity profile () of the fluid.(2) The influence of Hartmann number  on the velocity of the fluid is shown in Figure 3.It is concluded that the velocity increases by increasing  in case of  ̸ = 0.However, from Figure 4, the velocity decreases due to an increase in  for  = 0.This observation is in agreement with the literature (see [39,Figure 5]).
(3) Figure 5 shows the influence of the slip parameter  on the fluid velocity.By increasing , it is observed that the velocity of the fluid increases.

Figure 1 :
Figure 1: A steady squeezing axisymmetric fluid flow between two parallel plates.

Figure 2 :
Figure 2: Plot shows the Graph of () for different values of .

( 1 )
Figure 2 represents the influence of Reynold's number  on the velocity of the fluid.It is observed from Figure 2 that by increasing  the velocity of the fluid decreases.

Table 1 :
DTM and OHAM solutions along with the residuals.

Table 2 :
Residuals for different values of  using DTM.