Magneto Hydrodynamics Stagnation Point Flow of a Nano Fluid over an Exponentially Stretching Sheet with an Effect of Chemical Reaction , Heat Source and Suction / Injunction

A numerical investigation is carried out on the effects of heat source suction and viscous dissipation on Magneto hydrodynamics boundary layer flow of a viscous, steady and incompressible fluid. The flow is assumed to be over on exponentially stretching sheet. The governing system of partial differential equations has been transformed into ordinary differential equation using similarity transformation. Keller box method is simulated on the dimensionless system of differential equations. The skin friction coefficient and the heat and mass transfer rates are very significant parameters that are computed, analysed discussed in detail.


Introduction
Stagnation point is of considerable importance in most of industrial applications.The stagnation point is the one where fluid becomes stationary.In certain situations, flow is stagnated by a solid wall while in others; there is a line interior to a homogeneous fluid domain or the interface between two immiscible fluids [1]- [3].A good

Mathematical Formulation
Consider a steady, two-dimensional boundary layer stagnation-point flow of a Nano fluid over an exponentially stretching sheet, the stretching and free stream velocities are assumed to be of the forms ( ) e x l w u x a = , and ( ) e x l u x b ∞ = , respectively.Where a > 0 and b are constants, x is the coordinate measure along the stretching surface and l is length of the sheet.A non-uniform transverse magnetic field of strength ( ) = is imposed parallel to the y-axis, where Bo is a uniform magnetic field strength.It is assumed that the induced magnetic field due to the motion of an electrically conducting field is negligible.Further it is also assumed that the external electric field is zero and the electrical field due to polarization of charges is negligible [22].
Figure 1 shows that the temperature T and the Nano-particles fraction C take forms Tω(x) and Cω(x), respectively whereas the ambient values of temperature T ∞ and Nano-particle fraction C ∞ are obtained when y tends to infinity.
The governing boundary layer equations of the conservation Law of mass, momentum, energy and concentration in the flow as follows: Here u and v are the velocity components in the x and y directions respectively, µ is the viscosity, ρ is the density of the base field and σ is the electrical conductivity.is the radiation flux.The Rosseland approximation is defined as [23] [24]; where σ * is the Stefan-Boltzmann constant and k * is the mean absorption coefficient.It is assumed that the temperature difference between the free stream T ∞ and local temperature T is small enough expanding T 4 in a Taylor series about T ∞ and neglecting higher order terms results; After substituting Equations ( 5) and (6) in Equation ( 3), it will be reduced to ( ) ( ) the subjected boundary conditions are , , , at 0 e , , , as The prescribed temperature and concentration on the surface of the sheet is assumed to be ( ) where T 0 , C 0 are the reference temperature and concentration respectively, now, the non-linear partial differential equations for the purpose of a stream function ( ) where the continuity Equation ( 1) is satisfied identically.A similarity transformation is defined as [20] follows as such Equation (10), Equations ( 2), ( 4) and ( 7) reduce to the following system of nonlinear ordinary differential equations.
( ) where ( ) Here, prime denote the differentiation with respect to , B η is the velocity ratio parameter, v is the kinematic viscosity of the fluid, P r is the Prandtl number, Le is the Lewis number, M is the magnetic parameter, A is the heat source parameter, S is the suction parameter.where Nb is the Brownian motion parameter and Nt is the thermophoresis parameter, Ec Eckert number, γ is chemical reaction parameter; The corresponding boundary conditions Equation ( 8) are transformed into The parameters of practical interest in the formulated problem are velocity, heat and mass transfer respectively, which are presented in terms of Skin friction C f , Nusselt number Nu and Sherwood numbers Sh.Using the transformed variables (10), the non-dimensional expressions for the Skin friction coefficient

Numerical Procedure
The equations ( 11)-( 14) subject to the boundary conditions ( 15) are solved numerically using an implicit finite-difference scheme known as Keller box method.The method has the following four basic steps.1) Reduce equations ( 11)-( 14) to first order equations; 2) Write the difference equations using central differences; 3) Linearize the resulting algebraic equations by Newton's method and write them in matrix-vector; 4) Use the Block-tridiagonal elimination technique to solve the linear system.The system of ordinary differential equations ( 11)-( 13) has been solved numerically using Keller-box method.From the numerical computation, the main physical quantities of interest namely the local skin friction coefficient, the local Nusselt number and the local Sherwood number are obtained and the results are presented in Table 1, Table 2 and Table 3.
From Table 2, it is observed that with the increase in heat source parameter, there is a decrease in rate of heat transfer and increase in mass transfer.
From Table 3, it is observed that with increase in suction parameter, there is a increase in rate of heat transfer, mass transfer and skin friction coefficient.
Figure 2 shows the effects of the suction parameter "s" on the flow field velocity f′(η) for three different values of the velocity ratio parameter B, B = 0.1, 1 and 2.
When B = 0.1 the velocity profile and the boundary layer thickness decrease with an increases in "s".When s = 2.0, higher value of the Lorentz force further reduces the velocity and consequently the thickness of boundary         layer reduces.When B = 1, there is no influence of source on f′(η) and f′(η) attains a constant value of "1" for any values of η indicating that there is no boundary layer of fluid, as shown by a thick (-) line in Figure 2. Figure 3 and Figure 4 demonstrate the effects of suction parameter on temperature and concentration.It is clear that the fluid velocity decreases significantly with increasing values of suction parameter while it is found to enhance with blowing.The presence of suction would result in the reduction of the thickness of the boundary layer.So the Nano fluid temperature and concentration decreases with an increase in the suction parameter.
Figure 5 shows effects of the heat source parameter on temperature.Increasing the heat source parameter from −0.3 to 0.5 increases the temperature profile.Because the presence of source of heat enhances thermal energy, as a result of this, temperature profile increases.Figure 6 shows the heat source parameter significantly decreases the concentration profile.
Figure 7 and Figure 8 shows the effects of viscous dissipation parameter "Ec" on Temperature and concentration profile.Eckert number is the ratio of the kinetic energy of the flow to the boundary Layer enthalpy differ-   ence.The effect of viscous dissipation on flow field is to increase the energy, yielding a greater fluid temperature and as consequence greater buoyancy force.The increase in the buoyancy force due to an increase in the dissipation parameter enhances the temperature.While concentration decreases with an increase in the viscous dissipation parameter.
The influence of chemical reaction parameter "γ" on concentration profile is shown in Figure 9. Concentration decreases with an increase in the chemical reaction parameter indicating that the nanoparticle volume fraction decreases with the increase of chemical reaction parameter, while effect of chemical reaction parameter γ is   not significant on the temperature profile.

Conclusions
A numerical study corresponding to the flow and heat transfer in a steady flow region of Nano fluid over an exponential stretching surface and effects of heat source, suction parameter and Eckert number are examined and discussed in detail.The main observation of the present study as follows: 1) It is interesting to note that an increase in the suction parameter reduces the velocity profile and increase the temperature and concentration; 2) An increase in the heat source parameter increases the temperature profile while it reduces the concentration profile; 3) An increase in the Eckert number increases the temperature profile while it reduces the concentration profile.

Figure 1 .
Figure 1.Physical flow model and coordinate system.

(
capacitance of the nanoparticle.D B is the Brownian diffusion coefficient, D T is the thermophoresis diffusion coefficient and r q is the local wall Shear stress, ρ is the fluid density.In heat transfer at a boundary within a fluid, the Nusselt number (Nu) is the ratio of convection to conductive heat transfer across the boundary The Sherwood number (Sh) is also called the mass transfer Nusselt number.It represents the ratio of convective to diffusive mass transport Reynolds number based on the stretching velocity.

Figure 2 .
Figure 2. Velocity profile against η for different values of S.

Figure 3 .
Figure 3. Temperature profile against η for different values of S.

Figure 4 .
Figure 4. Concentration profile against η for different values of S.

Figure 5 .
Figure 5. Temperature profile against η for different values of A.

Figure 6 .
Figure 6.Concentration profile against η for different values of A.

Figure 7 .
Figure 7. Temperature profile against η for different values of Ec.

Figure 8 .
Figure 8. Concentration profile against η for different values of Ec.

Figure 9 .
Figure 9. Concentration profile against η for different values of γ.

Table 1 .
Comparison of the reduced Nusselt number

Table 2 .
Values of the reduced Nusselt number

Table 3 .
Values of the reduced Nusselt number