Application of Daftardar Jafari Method to First Grade MHD Squeezing Fluid Flow in a Porous Medium with Slip Boundary Condition

and Applied Analysis 3 Table 2: Absolute residuals for different values of Re keeping Me and γ fixed. z Absolute residuals for different values of Reynolds number (Re = 1.5) (Re = 2) (Re = 3.5) (Re = 4.0) (Re = 5.0) 0. 0. 0. 0. 0. 0. 0.1 8.72 × 10 1.63 × 10 7.04 × 10 1.44 × 10 4.40 × 10 0.2 4.53 × 10 9.03 × 10 3.92 × 10 8.03 × 10 2.45 × 10 0.3 1.98 × 10 3.82 × 10 1.70 × 10 3.52 × 10 1.10 × 10 0.4 3.67 × 10 5.52 × 10 2.62 × 10 5.47 × 10 1.75 × 10 0.5 4.50 × 10 4.24 × 10 2.27 × 10 4.83 × 10 1.60 × 10 0.6 4.05 × 10 2.04 × 10 1.40 × 10 2.00 × 10 1.01 × 10 0.7 3.00 × 10 6.20 × 10 6.44 × 10 1.43 × 10 5.03 × 10 0.8 1.60 × 10 7.30 × 10 2.52 × 10 6.00 × 10 2.12 × 10 0.9 7.31 × 10 4.40 × 10 8.62 × 10 2.05 × 10 8.00 × 10 1. 3.00 × 10 4.00 × 10 3.00 × 10 7.00 × 10 3.00 × 100 Table 3: Absolute residuals for different values of Me keeping Re and γ fixed. z Absolute residuals for different values of Hartmann number (Me = 1.5) (Me = 2) (Me = 3) (Me = 4.0) (Me = 5.0) 0. 0. 0. 0. 0. 0. 0.1 1.30 × 10 1.40 × 10 2.00 × 10 1.43 × 10 1.20 × 10 0.2 6.40 × 10 6.55 × 10 9.32 × 10 7.30 × 10 6.14 × 10 0.3 3.00 × 10 2.30 × 10 3.50 × 10 3.00 × 10 2.40 × 10 0.4 6.00 × 10 3.00 × 10 4.44 × 10 4.00 × 10 3.13 × 10 0.5 8.00 × 10 2.30 × 10 3.20 × 10 3.00 × 10 2.32 × 10 0.6 8.00 × 10 1.45 × 10 1.60 × 10 1.40 × 10 1.20 × 10 0.7 5.51 × 10 7.60 × 10 6.20 × 10 5.31 × 10 5.00 × 10 0.8 3.13 × 10 3.10 × 10 2.12 × 10 2.00 × 10 1.55 × 10 0.9 1.50 × 10 9.00 × 10 7.00 × 10 4.00 × 10 4.40 × 10 1. 6.00 × 10 1.10 × 10 2.20 × 10 1.10 × 10 1.10 × 100 Introducing stream function ψ(r, z) [12] and the generalized pressure p = ρu + p, (4) and (5) are reduced to ∂p ∂r − ρ r2 ∂ψ ∂r E 2 ψ = μ r ∂ ∂z E 2 ψ − A r ∂ψ ∂z ∂p ∂r − ρ r2 ∂ψ ∂z E 2 ψ = − μ r ∂ ∂r E 2 ψ + A r ∂ψ ∂r , (7) where E 2 = ∂ 2 ∂r2 − 1 r ∂ ∂r + ∂ 2 ∂z2 . (8) Eliminating p from (7) and using the transformation ψ(r, z) = r 2 y(z), we have y (iV) (z) − ( 1 k + ηB 2


Introduction
Movable boundary stressed by vertical velocity causes squeezing flow [1].The most practical examples of squeezing flows are polymer processing, compression, and injection modeling.The preliminary results in this field are carried out by Stefan [2].The Newtonian squeezing fluid flow between two parallel plates is studied in [3].The study of thin Newtonian liquid films squeezed between two plates is investigated by Grimm (Figure 1) [4].Squeezing flow, under the influence of magnetic field, has many applications in the field of chemical engineering and is experimented by many researchers [5][6][7][8].
Molecules near the boundary show slip at the boundary when the weight of molecule is high.In case of thin film problems, rarefied fluid problems, fluids containing concentrated suspensions, and flow on multiple interfaces, no-slip boundary condition does not work.The general boundary condition, which shows the fluid slip at the surface, is introduced by Navier [9].The difference of the velocities of fluid and of the boundary, according to Navier, is proportional to the shear stress at that boundary.This proportionality constant is called slip parameter.
The slip condition has its significance when fluids with elastic character are considered [10].In medical sciences, particularly in polishing artificial heart valves and internal cavities, the slip condition is highly considered [11].The effect of slip condition on MHD squeezing fluid flow is discussed in [12].For rotating flow of a third grade fluid in a nonporous medium, the effects of slip condition are investigated by Hayat and Abelman [13].The same model with the extension of porous medium is argued by Abelman et al. [14].Newtonian fluid was considered by Ebaid [15] to study the effects of magnetic field and wall slip conditions on the peristaltic transport in an asymmetric channel.
There are many applications of magnetohydrodynamics (MHD) of an electrically conducting fluid in geophysics, astrophysics, engineering, and other industrial areas.More to the point, the flows of electrically conducting fluid through porous medium have acquired exceptional rank and have been the spotlight of interest of many researchers.The precise applications are initiated in the study of ground water flow, irrigation problems, crude petroleum recovery, heat-storage beds, thermal and insulating engineering, chromatograph, and chemical catalytic reactors.The idea of DJM is introduced by Daftardar-Gejji and Jafari [16,17].This scheme (DJM) is computer-friendly and effective to deal with highly nonlinear boundary value problems [18].The convergence of DJM can be found in [19] where fractional boundary value problems with Dirichlet boundary conditions are considered.To compare different analytic solutions for squeezing fluid flow between two parallel plates, DJM is successfully used [20].Using this scheme, excellent results are obtained in [21] for the solution of fifth-and sixth-order boundary value problems.
The objective of this research paper is to use DJM for studying, with slip boundary condition, first grade MHD squeezing fluid flow in a porous medium bounded by two parallel plates approaching each other slowly.The basic idea of the proposed scheme is discussed and the modeled boundary value problem is solved through it.The results along with residuals, for different values of parameters, are shown in Tables 1, 2, 3, 4, and 5.The effects on velocity profile of the fluid, by the influence of different parameters, are shown in the graphs.

Mathematical Model
The Navier-Stokes equations for such fluid flow between two large parallel plates separated by 2ℎ distance are where  denotes the constant density,  is viscosity, k is the velocity of plates, ∇ denotes the material time derivative,  denotes the pressure,  is permeability,  is the total magnetic field given by  =  0 + , and  denotes the electrical conductivity while  is given by In the total magnetic field ,  0 and  are imposed and induced magnetic fields respectively.Consider  as negligible as compared to  0 , , , and  constant and B is perpendicular to k with no electric field in the fluid flow region; the involved magnetohydrodynamic force can be written as Here  is the electric current density.Magnetic field is applied along the -axis.In cylindrical coordinates, the components of the velocity k are given by k -component: where Introducing stream function (, ) [12] and the generalized pressure p =  + , (4) and ( 5) are reduced to where Eliminating p from (7) and using the transformation (, ) =  2 (), we have with boundary conditions Using the nondimensional parameters and omitting the * , (9) and ( 10) can be written as with  = /ℎ.Re and Me are Reynolds and Hartmann numbers, respectively.

Analysis of Daftardar Jafari Method
Consider fourth-order nonlinear boundary value problem with boundary conditions Using the operator  =  4 / 4 , ( 16) can be written as Initial conditions are used when the inverse operator  −1 is applied on (18) to get Abstract and Applied Analysis 5 In series form (19) can be written as are defined in (14).The components of () can be found as For the solution of our problem, we rewrite ( 12), using the inverse operator  −1 which is a fourfold integral with initial conditions (i.e., at  = 0), as In series form, (22) reduces to The last equation gives the components of ().For the nonlinear term () = () ⋅   (), DJM introduces     which are already defined in (14).The recursive process for the components of () is given by Following the above recursive process, the respective components are obtained as follows: )) = 1, )) = 0. (27) The above system of nonlinear equations is solved by using Mathematica software to get  = 1.45539,  = 0.754966 (29)

Numerical Results and Discussion
The velocity profile of MHD squeezing fluid flow between two parallel plates in porous medium with slip boundary is deduced by using Daftardar Jafari method.All the calculations, tables, and figures are carried out by using Mathematica    scheme.This technique requires less computational effort and the convergence of this method to the exact solution is fast enough.The comparison of DJM with other techniques can be found in [12,21].For getting the insight of the problem, we discuss the velocity profile () of the fluid as follows.
(1) In Figure 2  (3) Figures 5, 6, and 7 are special cases when the slip parameter  vanishes.From Figure 5 we conclude that the velocity increases with equal increase in Me and Re.In the same way the increase in velocity is observed in Figures 6 and 7 while increasing Me and Re, keeping Reynolds and Hartmann numbers fixed, respectively.

Figure 1 :
Figure 1: Squeezing flow of Newtonian fluid between two large parallel plates.

Me = 10 ,Figure 2 :
Figure 2: Velocity profile for different values of Reynold and Hartmann numbers with Re ≤ Me.

7 . 0 .Figure 3 :Re = 6 ,Figure 4 :Re
Figure 3: Velocity profile for different values of slip parameter and Reynold and Hartmann numbers with Re ≥ Me.

Figure 5 :
Figure 5: Velocity profile for different values of Reynold and Hartmann numbers with Re = Me and  = 0.

Figure 6 :
Figure 6: Velocity profile for different values of Hartmann number with Re = 1 and  = 0.

Figure 7 :
Figure 7: Velocity profile for different values of Reynold number with Me = 1 and  = 0.

( 2 )
the influence of Hartmann and Reynolds numbers with constant slip parameter is shown.By decreasing Hartmann number with Me ≥ Re, we observe from Figure 2 that the velocity of the fluid decreases.The influence of Reynolds and Hartmann numbers for different values of slip parameter is shown in Figure 3.It can be seen in Figure 3 that by increasing Reynolds and Hartmann numbers with Re ≥ Me and decreasing the slip parameter , the velocity decreases.The same effect on the velocity profile is shown in Figure 4 by decreasing Reynolds and Hartmann numbers with Re < Me.

Table 1 :
Absolute residuals for different values of  keeping Me and Re fixed.

Table 2 :
Absolute residuals for different values of Re keeping Me and  fixed.

Table 3 :
Absolute residuals for different values of Me keeping Re and  fixed.

Table 4 :
Absolute residuals for different values of Me and Re with Re > Me keeping  fixed.

Table 5 :
Absolute residuals for different values of Me, Re, and  with Re < Me.