Numerical Solution of MHD Flow of Micropolar Fluid with Heat and Mass Transfer towards a Stagnation Point on a Vertical Plate

The paper investigates the numerical solution of problem of magnetohydrodynamic (MHD) micropolar fluid flow with heat and mass transfer towards a stagnation point on a vertical plate. In this study, we consider both strong concentrations (n = 0) and weak concentrations (n = 1/2). The governing equations have been transformed into nonlinear ordinary differential equations by applying the similarity transformation and have been solved numerically by using the finite difference method (FDM) and analytically by using (DTM). The effects of various governing parameters, namely, material parameter, radiation parameter, magnetic parameter, Prandtl number, Schmidt number, chemical reaction parameter and Soret number on the velocity, microrotation, temperature and concentration have been computed and discussed in detail through some figures and tables. In order to verify the accuracy of the present results, we have compared these results with the analytical solutions by using the differential transform method (DTM) and the multi-step differential transform method (MDTM). It is observed that this approximate numerical solution is in good agreement with the analytical solution.


Introduction
The hydromagnetic stagnation point flows past heated or cooled bodies have attracted many researchers due to their applications in most of the engineering and natural processes.Examples include plasma studies, blood flow problems, MHD generators, and the cooling of an infinite metallic plate in a cooling bath.The classical two dimensional stagnation point flow on a flat plate first studied by Hiemenz [1] was extended to axisymmetric case by Homann [2].The problem of two dynamic impinging stagnation flows of two fluids of different densities was investigated by Wang [3] and it was shown that the flow field depends heavily on the ratios of viscosity and density.The influence of an external magnetic field on Hiemenz flow was studied by Ariel [4].Two dimensional MHD steady stagnation point flow towards a stretching surface was analyzed by Mahapatra and Gupta [5].Ishak et al. [6] theoretically studied the similarity solutions for the steady magnetohydrodynamic flow towards a stagnation point on a vertical surface immersed in an incompressible micropolar fluid.Lokendra Kumar [7] has analyzed the effect of MHD flow of micropolar fluid towards a stagnation point on a vertical stretching sheet.Ishak et al. [8] investigated the steady two-dimensional stagnation point flow of an incompressible micropolar fluid towards a vertical stretching sheet.
The theory of micropolar fluids which is originally formulated by Eringen [9] can be used to explain the flow of crystals, animal blood, paints, polymers, etc.The theory introduces new material parameters, an additional in dependent vector field the microrotation-and new constitutive equations which must be solved simultaneously with the usual equations for Newtonian flow.Ramachandran et al. [10] studied laminar mixed convection in two-dimensional stag nation flows around surfaces.He considered both cases of an arbitrary wall temperature and arbitrary surface heat flux variations and found that they are versed flow developed in the buoyancy opposing flow region, and dual solutions are found to certain range of the buoyancy parameter.Hassanien et al. [11] extended Ramachandran's work to micropolar fluid.They considered both assisting and opposing flows, but the existence of dual solutions was not reported [12].Devi et al. extended the problem posed by Ramachandran et al. [10] to the unsteady case, and they found that dual solution exists for a certain range of the buoyancy parameter when the flow is opposing.Similar problem for steady and unsteady cases, for a vertical surface immersed in a micropolar fluid was investigated by Lok et al. [13] [14].Existence of dual solutions was reported in [13] only for the opposing flow regime.Mahapatra et al. [15] examined heat transfer in stagnation-point flow on stretching sheet with viscous dissipation effect.Attia [16] studied the hydromagnetic stagnation point flow on porous stretching sheet with suction and injection.El-dabe et al. [17] studied the peristaltic motion of incompressible micropolar fluid through a porous medium in a two-dimensional channel under the effects of heat absorption and chemical reaction in the presence of magnetic field.Effects of thermal radiation on magnetohydrodynamic (MHD) flow of a micropolar fluid towards a stagnation point on a vertical Plate were studied by Olanrewaju et al. [18].Pop et al. [19] examined the flow over a stretching sheet near a stagnation point taking radiation effect.Olanrewaju and Adesanya [20] studied the effects of radiation and viscous dissipation on stagnation flow of a micropolar fluid towards a vertical permeable surface.Manjoolatha et al. [21] studied the radiation and mass transfer effects on MHD flow of a micropolar fluid towards a stagnation point on a vertical stretching sheet.
The main objective of the present work is to study the numerical solutions of problem of magnetohydrodynamic (MHD) micropolar fluid flow with heat and mass transfer towards a stagnation point on a vertical plate.In this study, we consider both strong concentrations (n = 0) and weak concentrations (n = 1/2).The governing equations have been transformed into nonlinear ordinary differential equations by applying the similarity transformation and have been solved numerically by using the finite difference method (FDM) and analytically by using (DTM).The effects of various parameters of the problem on fluid velocity, microrotation, temperature and concentration have been computed and discussed in detail through some figures and tables.Comparisons with analytical solution by using (DTM) of present work are performed and show that the present results apply with results of analytical solution, as well as the comparisons with previous published results and are found to be in good agreement.

Mathematical Formulation
We consider the steady, two-dimensional flow of an incompressible electrically conducting micropolar fluid near the stagnation point on a vertical heated stretching sheet placed in the (y = 0) of a cartesian coordinate system with the x-axis along the sheet as shown in Figure 1.The fluid occupies the half plane (y > 0).It is assumed that the velocity of the flow external to the boundary layer U, the temperature of the sheet w T and concentration of sheet w C are proportional to the distance x from the stagnation point, i.e. , and where a, b and c are positive constants, w T T ∞ > with T ∞ being the uniform temperature of the fluid and w C C ∞ > with C ∞ being the uniform concentration of the fluid.A uni- form magnetic field of strength 0 B is assumed to be applied in the positive y-direction normal to the plate.The magnetic Reynolds number of the flow is taken to be small enough so that the induced magnetic field is negligible.
Under the usual boundary layer approximation, the governing equations are given as follows: ( ) The boundary conditions are as follows: ( ) 0, 0, , , , at 0 , 0, , as where u and v are the velocity components along the x and y direction, g is the acceleration due to gravity, T is the fluid temperature, w T is the plate temperature, T ∞ is the ambient temperature, C is the fluid concentration in the boundary layer, w C is the plate concentration, C ∞ is the ambient fluid concentration, ν is the kine- B is the magnetic field of constant strength, σ is the electric conductivity of the fluid, r q is the radiative heat flux, p c is the specific heat at constant pressure of the fluid, m D is the mass diffusion coef- ficient, m T is the mean fluid temperature, T k is the thermal diffusivity ratio, s c is concentration susceptibility, and r k is the chemical reaction coefficient.Furthermore, , , , , , and k j N µ ρ γ α are respectively the dynamic viscosity, vortex viscosity (or the microrotationviscosity), fluid density, microinertia density, microrotation vector (or angular velocity), spingradient viscosity and thermal diffusivity.
We follow the work of many authors by assuming that 1 2 2 is the material parameter.This assumption is involved to allow the field of equations predicts the correct behavior in the limiting case when the microstructure effects become negligible and the total spin N reduces to the angular velocity (see Ishak et al. [6]).By using the Rosseland approximation (Brewster [22]), the radiative heat flux r q is given by given by where * σ is the Stefan-Boltzmann constant and * k is the mean absorption coefficient.It should be noted that by using the Rosseland approximation, the present analysis is limited to optically thick fluids.If temperature differences within the flow are sufficiently small, then Equation ( 7) can be linearized by expanding 4   T into the Taylor series about T ∞ , which after neglecting higher order terms takes the form In view of Equations ( 7) and ( 8), Equation ( 4) reduces to where k is the constant thermal conductivity.

Method of Solutions
To solve system of Equations ( 1)-( 5), we will consider the following similarity transformations: where η is the independent similarity variable, ( ) f η the dimensionless stream function, ( ) h η the dimen- sionless micorotation, ( ) θ η the dimensionless temperature and ( ) η ∅ the dimensionless concentration.Further, ϕ is the stream function which is defined in the usual way where prime denotes differentiation with respect to η .By using Equation (10) in the Equations ( 1)-( 5) and ( 9), we get the following ordinary differential equations in dimensionless form: ( ) where is the thermal buoyancy parameter, ( ) ± sign has the same meaning as in Equation ( 2).For ‫ג‬ > 0, buoyancy forces act in the direction of the mainstream and fluid is accelerated in the manner of a favorable pressure gradient (assisting flow).When ‫ג‬ < 0, buoyancy forces oppose the motion, retarding the fluid in the boundary layer, acting as an adverse pressure gradient (opposing flow), is the solutal buoyancy parameter, ( ) is the local Reynolds number, ( ) is the Eckert number,  ( ) Also, the subjected boundary conditions Equation ( 6) will now take the form: The skin-friction coefficient f C , local Nusselt number ux N and Sherwood number hx S are defined as: In which , , w w w q j τ are given by ( ) Substituting by similarity transformations and Equation (18) into Equation ( 17), we get

The Differential Transform Method
The differential transformation of an analytical function ( ) f t for one variable is defined as [23].
( ) ( ) where f(t) is the original function and F(k) is the transformed function.
The differential inverse transformation of F(k) is defined as: Combining Equations ( 20) and ( 21), we obtain From Equations ( 20)-( 22), it can be seen that the differential transformation method is derived from Taylor's series expansion, but the method does not calculate the derivatives representatively.However, the relative derivatives are calculated by an iterative way which is described by the transformed equations of the original function.For implementation purposes, the function f(t) is expressed by a finite series and Equation ( 21) can be written as By Equation ( 20), the following theorems can be deduced:

Basic Concepts of the Multi-Step Differential Transform Method (MDTM)
When the DTM is used for solving differential equations with the boundary condition at infinity or problems that have highly non-linear behavior, the obtained results were found to be incorrect (when the boundary-layer variable go to infinity, the obtained series solutions are divergent).Besides that, power series are not useful for large values of the independent variable.To overcome this shortcoming, the MDTM that has been developed for the analytical solution of the differential equations is presented in this section.For this purpose, the following non-linear initial-value problem is considered: , , , , 0 subject to the initial conditions ( ) ( ) . Let [0, T] be the interval over which we want to find the solution of the initial-value problem (24).In actual applications of the DTM, the approximate solution of the initial value problem (24) can be expressed by the following finite series: The multi-step approach introduces a new idea for constructing the approximate solution.Assume that the in- of equal step size h = (T/M) by using the nodes t m = mh.The main ideas of the MDTM are as follows.First, we apply the DTM to Equation (24) over the interval [0, t 1 ], we will obtain the following approximate solution: using the initial conditions ( ) ( ) , where 0 t in Equation ( 24) is replaced by 1 m t − .The process is repeated and generates a sequence of approximate solutions ( ) , for the solution f(t): where . In fact, the MDTM assumes the following solution: ( ) The new algorithm, MDTM, is simple for computational performance for all values of h.It is easily observed that if the step size h = T, then the MDTM reduces to the classical DTM.As we will see in the next section, the main advantage of the new algorithm is that the obtained series solution converges for wide time regions and can approximate non-chaotic or chaotic solutions.

Analytical Solution Be the MDTM
By applying the MDTM to Equations ( 12)-( 15), the following recursive relations in each sub-domain ( ) where The differential transformed boundary conditions in Equation ( 16) to: where 1 2 1 1 1 , , , , f f h t c are constants.These constants are computed from the boundary condition.Moreover, substituting Equation (33) into Equations ( 29)-(32) and by using the recursive method, we can calculate other values of Hence, substituting all ( ) ( ) ( ) , ( ), and 23), we obtain series solutions.
The velocity, microrotation, temperature and concentration distributes achieved with the aid of MATHAMATICA application software.

Numerical Solution
The governing ordinary differential Equations ( 10)-( 12) are being highly nonlinear and difficult to solve analytically.To solve these coupled equations, we use a finite difference based numerical algorithm.Following Ashraf et al. [24], we reduce the order of Equation ( 12) by one with the help of the substitution: Equations ( 12)-( 15) in view of Equation ( 34) can be written as: The boundary conditions: Now, we can write Equations ( 35)-(38) in the following finite difference method form: ( ) ( ) ( ) The velocity, microrotation, temperature and concentration distributes at all interior nodal points computed by successive applications of the above finite difference equations and these are achieved with the aid of MATLAB application software.

Results and Discussion
In this section, we present our results in tabular and graphical form is shown in Figures 2-16 to illustrate the influence of physical parameters such as, material parameter K, magnetic parameter M, thermal buoyancy ‫,ג‬ solutal buoyancy parameter δ , radiation parameter a R , Prandtl number r P , Dufour number u D , Schmidt- number c S , chemical reaction r K and Soret number o S , on the velocity, microrotation, temperature and con- centration when n = 0 and n = 1/2 for both the cases of assisting and opposing flows.
Figure 2 and Figure 3 illustrate the effect of the material parameter K on the normal velocity profiles ( ) It is seen that the normal velocity profile decreases with an increasing in the material parameter K for both assisting and opposing flows.It is clear from the figure that the normal velocity profile is higher for assisting flow when compared with the normal velocity profile for opposing flow.Figures 4-7 represent the velocity profiles for the flow parameters.It is observed that the velocity profile increases with increasing the magnetic parameter M, while it decreases with increasing the material parameter K for both assisting and opposing flows.It is clear from the figures that the velocity is higher for assisting flow when compared with the velocity for opposing flow.In order to verify the accuracy of the numerical solution of present work by using (FDM) we have compared   and Olanrewaju et al. [18].It is observed that this numerical solution in excellent agreement with the results previous obtained values as shown in Table 1 and Table 2.
Finally, Tables 3-5 show comparison between numerical solution by using (FDM) and the analytical solu-   for some values of K r .
tion using (DTM).It is observed that this approximate numerical solution is in excellent agreement with the results of the analytical solution of by using (DTM).[18] and our results for K = 0, M = 0, ‫ג‬ = 1 and R a = 0.   [18] and our results for K = 0, M = 0, ‫ג‬ = 1 and R a = 0.

Conclusions
In this work, we have obtained the numerical solution of problem of magnetohydrodynamic (MHD) micropolar fluid flow with heat and mass transfer towards a stagnation point on a vertical plate.In this study, we consider both strong concentrations (n = 0) and weak concentrations (n = 1/2).The resulting partial differential equations which describe the problem are transformed into ordinary differential equations by using a similarity transformation and then solved numerically by using the finite difference method (FDM) and analytically by using DTM.A representative set of numerical results for velocity, microrotation, temperature and concentration profiles is presented graphically and discussed.The figures and tables clearly show that the results by using (FDM) are in excellent agreement with previously published works and with the results of the analytical solution by using DTM.The important results for this study are summarized as follows: • The normal velocity profile decreases with the increase of the material parameter K for both assisting and opposing flows.• The velocity distribution increases with the increase of M, while it decreases with the increase of K for both assisting and opposing flows.• The microrotation distribution increases with increase of K for assisting flow, while it increases with an increase of K in the vicinity of the surface but the reverse happens away from surface for opposing flows, when n = 1/2.• The microortation decreases with an increase K for both assisting and opposing flows when n = 0.
• The temperature distribution increases with the increase of R a and D u , while it decreases with the increases of

Figure 1 .
Figure 1.Physical model and coordinate system.
the Chemical reaction parameter, ( )

Figure 3 .
Figure 3. Normal velocity profiles ( ) f η for some values of K when n = 0.

Figure 6 .
Figure 6.Velocity profiles ( ) f η for some values of K when n = 0.5.

Figure 7 .
Figure 7. Velocity profiles ( ) f η for some values of K when n = 0.

Figure 8 .
Figure 8. Angular velocity profiles ( ) h η for some values of K when n = 0.5.

Figure 9 .
Figure 9. Angular velocity profiles ( ) h η for some values of K when n = 0.5.

Figure 10 .
Figure 10.Angular velocity profiles ( ) h η for some values of K when n = 0.

Figure 11 .
Figure 11.Temperature profiles ( ) θ η for some values of P r .

Figure 12 .
Figure 12.Temperature profiles ( ) θ η for some values of R a .

Table 3 .
Comparison between the previous obtained values of ( ) ( )

Table 4 .
Comparison between the previous obtained values of ( ) ( )

Table 5 .
Comparison between the previous obtained values of