The laminar flow and heat transfer of an incompressible fluid over a stretching surface have significant applications in several technologies of process industry. A number of manufacturing industries process and include both metal and polymer sheets during the production of sheeting materials. This kind of motion of a fluid close to the rigid or stretching surface can be investigated from a standard boundary layer equation. The quality of the final product depends to a large extent on the rate of heat transfer at the stretching surface. It is, therefore, of interest to know the fluid behavior over a stretching surface which determines the rate of cooling. Hiemenz [1] first analyzed an exact solution of the stagnation point flow of a viscous fluid over a rigid surface. The axisymmetric three-dimensional cases for the rigid surface are investigated by Homann [2] . The steady flow of an incompressible viscous fluid over an elastic sheet was studied by Crane [3] where the flow being caused solely by stretching of sheet in its own plane. The temperature distribution in the flow over a stretching surface subject to uniform heat flux was studied by Dutta et al. [4] . The steady two-dimensional orthogonal stagnation point flow of an incompressible viscous fluid towards a stretching surface which is held at constant temperature, was investigated by Mohapatra and Gupta [5] . Later, Mohapatra and Gupta [6] extended his work to study the momentum and heat transfer characteristic of stagnation-point flow of a viscoelastic fluid over a stretching surface. Hassanien [7] analyzed the similarity solution for the flow and heat transfer of a viscoelastic fluid over a stretching sheet. The temperature distribution for the stagnation point flow of a micro polar fluid past a stretching sheet was investigated by Nazar et al. [8] .

In the region of adverse pressure gradient the flow in a boundary layer separates. Due to separation, it leads to increase in drag on the body immersed in the flow and adversely affects heat transfer from the surface of the body. Flow separation can be prevented by applying suction to remove the directed fluid particle from the boundary layer. On the other hand, the wall shear stress and hence friction drag is reduced by blowing. Wall suction and blowing have many applications such as turbine blade cooling, transition delay and prevention separation. The effect of suction on the Hiemenz problem was studied by Schlichting and Bussmann [9] . An approximate solution to the problem of uniform suction is given by Ariel [10] . The effect of uniform suction on the Homann problem where the flat plate oscillates in its own plane is considered by Weidman and Mahalingam [11] . The effect of uniform suction or injection on the two- or three-dimensional stagnation-point flow of a conducting fluid was given by Attia [12] [13] in the presence of an externally applied uniform magnetic field. The effect of uniform suction or injection on two-dimensional stagnation point flow towards a stretching surface with heat generation was given by Attia and Seddeek [14] .

Electrically conducting fluid like a liquid metal, electrolyte or plasma under the influence of magnetic field induces electric current. This type of fluid has many applications in industries, for example, to drive flow, induce stirring, levitation or control heat transfer and turbulence. An exact similarity solution to the MHD boundary layer equations for the steady two-dimensional flow of an electrically conducting incompressible fluid due to the stretching of a plane elastic surface in the presence of a uniform transverse magnetic field was given by Pavlov [15] . Chakrabarti and Gupta [16] studied the MHD flow of a Newtonian fluid over a stretching sheet. In the presence of uniform transverse magnetic field Andersson [17] investigated the MHD flow of a viscoelastic fluid past a stretching surface. Liao [18] developed Homotopy analysis method (HAM) which is one of the most successful and efficient methods to solve a non-linear differential equation. This method is based on a fundamental concept in topology, i.e., Homotopy (Hilton [19] ) which is widely used in numerical techniques (Grigolyuk and Shalashilin [20] ).

In this present work, HAM is employed to find the new analytic solution of the effect of uniform transverse magnetic field as well as suction or blowing on a steady two-dimensional orthogonal stagnation point flow of an electrically conducting fluid towards a permeable stretching surface. For a very small magnetic Reynolds number the induced magnetic field can be neglected. We assume the wall and free stream temperature to be constant.

Consider the two dimensional stagnation point flow of a viscous incompressible electrically conducting fluid impinging orthogonally to a permeable stretching surface coinciding with plane, the being in the region. Two equal and opposite forces are applied to the sheet along the x-axis such that the surface is stretched keeping the origin fixed as shown in Figure 1. The flow is permeated by a uniform magnetic field of strength applied transverse to the plate. The plate is assumed to be electrically non-conducting and is subject to uniform suction or blowing.

Let and are the velocity components along and directions, respectively. The velocity in the free flow is given by

where is proportional to the free stream velocity far away from the surface. Under these assumptions, along with the Boussinesq approximation and boundary layer approximations, the governing equations for the steady two-dimensional stagnation point flows by introducing Lorentz force are

where represents the temperature and, , , , and are the density, the viscosity, the kinematic viscosity, these electric conductivity of the fluid, thermal diffusivity and the specific heat of the fluid at constant pressure, respectively, and is horizontal component of the inviscid potential flow velocity over the body surface. At the stretching surface, the velocity in the direction is equal to the moving surface velocity, where is a positive constant. So boundary conditions are

where is the constant velocity at the plate with for suction and for blowing and and are constants with.

To examine the flow regime, we introduce the following similarity transformation

Here prime denotes differentiation with respect to. Using above equation, the continuity Equation (2) is satisfied automatically and Equations (3) and (4) becomes

where is the Hartman Number. Setting

In Equation (9) and equating the coefficients of and we obtain the following equations for and:

where, is the Eckert number which is a measure of viscous dissipation.

The corresponding boundary conditions (5) and (6) can be written as

Here is the suction parameter (for suction and for injection). For boundary layer flow, the wall skin friction is given by

Using as a characterized velocity, the skin friction coefficient can be defined as

Substituting (7) and (16) in (17) we get

where is the local Reynolds number.

The velocity and the temperature, can be represented in terms of a set of base function

as

and

where, and are the coefficients. Looking into the boundary condition and the base function we choose the initial guess, and as

and

The auxiliary linear operators, and can be taken as

where is an embedding parameter. The auxiliary linear operators, and satisfies

where, are arbitrary constants.

Corresponding to the governing differential Equations (8), (11) and (12), we define the non-linear operators as

If, , denotes the non-zero auxiliary parameters and, , denotes the auxiliary function, we construct the zero-th order deformation equations

subject to the boundary conditions,

where

When and we have

Note that when the embedding parameter increases from 0 to 1, , and varies from the initial approximations, and to the exact solution, and respectively. By Taylor’s theorem we express, and in the power series of the embedding parameter.

where

If the auxiliary parameters, and are properly chosen such that the series converges at then we have

Differentiating the zero-th-order Equations (34), (35), (36) -times with respect to the embedding parameter, then setting, and finally dividing by, we have -th-order deformation equation

subject to the boundary conditions

where

and

Let, , denote the special solution of Equations (52), (53) and (54). We have

where, , denotes the inverse operator of respectively. Thus the solution of the -th- order deformation equations as

Here the integral constants, are determined by the boundary conditions (55), (56) and (57) as

Therefore, it is easy to solve linear non-homogeneous Equations (52), (53) and (54) subject to the boundary conditions (55), (56) and (57) by using MATHEMATICA one after other in the order

We observed that the Equations (52), (53) and (54) consist of the auxiliary parameters, and. These auxiliary parameters are important in controlling the convergence of series solutions the functions, and. For this analysis, the admissible values of, and curves are plotted for 12th- order of approximations. Figure 2 shows the admissible values of of for 12th-order of approximation with for several values of. It is observed that curves have parallel lines segment that correspond to the regions, and for, respectively. To analyze the admissible range of and of and are displayed in Figure 3 for 12th-order of approximation with, and. It is observed that the converging ranges are and. From these figures we have examined that by choosing the values of convergence control parameter from these range we will get the convergent result up to higher decimal place.

The nonlinear boundary value problem given by Equations (8), (11) and (12) with boundary conditions (13), (14), and (15) cannot be solved in closed form. So these equations are solved by Homotopy Analysis Method (HAM) to get an analytical solution. Extensive calculation has been performed by using symbolic computation software MATHEMATICA to obtain the velocity and temperature profiles for various values of physical parameter such as Suction parameter A, Hartmann Number H, Prandtl number Pr, Eckert number E, (which characterizes viscous dissipation in the flow) and the parameter (ratio of external velocity and stretching velocity). Figure 4 shows the variation of horizontal velocity component with distance from the surface for several values of Hartmann Number H in the presence of suction parameter and i.e., stretching velocity more than the straining velocity. It can be viewed that horizontal velocity decreases with increase in Hartmann Number H. It is interesting that the flow has inverted boundary layer structure because of stretching velocity exceeds the velocity of external stream (i.e.). Whereas, in the presence of suction parameter and, the velocity increases with increase in H, as shown in Figure 5.

Table 1 shows the value of dimensionless wall shear stress i.e., skin-friction coefficient (see Equation (18)) for several values of H and in the presence of suction. It gives that for fixed value of suction parameter, the magnitude of the wall shear stress increases with increase in H and. In order to have correctness of our analytical result, we have examined the comparison result in the value of for different values of A and H in the present study and with the corresponding values by Mahapatra et al. [21] (shown in Table 1).

It is observed that during the motion of electrically conduction fluid, a certain amount of energy is stored up in the fluid as strain energy and some energy is lost due to viscous and ohmic dissipation. So in the fluid flow, energy is balanced by the internal energy, the conduction of heat, the convection of heat with the flow, the generation of heat with the flow, the generation of heat through viscous and ohmic dissipation, the strain (or deformation) energy in the fluid and energy due to wall suction effect. At a distance (Dimensionless length) from stagnation-point, the temperature distribution consists of two functions, first part of temperature distribution denoted by which depends on the Hartmann Number and suction parameter and second part denoted by which depends on viscous dissipation and Joule Heating.

In the presence of suction parameter, the variation of temperature distribution with is displayed in Figure 6 for several values of Hartman Number H with,. It is observed that for fixed value of Pr, , the temperature at point increases with increase in the Hartman Number when. But Figure 7 shows that when, decreases with increase in H. In the presence of wall

suction, viscous and ohmic dissipation and stored deformation energy in the flow, the variation of the temperature distribution given by with is shown in Figure 8, for several values of H with and. It is observed that a novel result, that in the presence of suction parameter, increases with the increase in Hartmann Number H where as in the absence of suction parameter Mahapatra et al. [22] shows that in the vicinity of the stretching surface increases with increase in H but after a certain distance from the stretching plate, the trend reverse. Figure 9 depicts that the variation of with for several values of H with and when.

The dimensionless rate of heat transfer at the surface given from Equation (10) is derived as

where. The dimensionless heat flux at the surface is computed from Equation (70) and displayed in Figure 10 and Figure 11 corresponding two distinct locations and on the stretching surface. It is observed from Figure 10 that the heat flows from the surface to fluid at the small value (for several values of H), because. On the side, for, it is seen in Figure 11 that the heat flow from the fluid to the stretching surface for several value of H because is negative. It is more interesting that heat flux is more in the presence of suction. This novel result can be interpreted physically that for small value of, the term becomes small and both ohmic dissipation and strain energy in the flow are not playing important roles to generate heat. So for small value of, heat flows from the surface to fluid since. But for the large value of, the term has significant role due combine effect of ohmic dissipation and strain energy inside the boundary layer which causes to generate sufficient heat inside the boundary layer. Thus, the temperature of the fluid near the surface exceeds the surface temperature. So heat flows from fluid to surface even when. In this case it is also observed that heat flux increase with suction for fixed value of H.

Since, for small value of, heat flows from the surface to the fluid and for the large value of, heat flows from the fluid to surface, so it can be observed in Figure 12 that at a certain value of on the stretching surface, heat flux vanishes. It is observed from this figure that value of decreases with increase in H or E. This interesting result admits of a physical interpretation. Near stagnation point, heat flows from the surface to the fluid but due to increase in the Hartmann number and Ekart Number, ohmic dissipation and strain energy will increase and that will cause to generate the heat inside the boundary layer. So that the value of will decrease with increase in H or E.

The present study describes the effect of viscous dissipation and joule heating on magneto hydrodynamic stagnation-point flow towards permeable stretching surface. To analyze the new analytic solution of this problem, we apply homotopy analysis method. By this powerful and newly developed technique, the convergence series solutions are obtained. To validate the analytic solutions of velocity distribution and temperature distribution using HAM method, we have computed the convergence regions. The HAM solutions have an excellent agreement. A novel result of this problem is that the temperature increases with the increase in Hartmann Number H at a certain distance from the stretching surface in the presence of suction parameter. Heat flows from the surface to the fluid at near stagnation point on the surface and on the hand far from the stagnation-point, heat flows from the fluid to surface due to combining effect of ohmic dissipation and strain energy inside the boundary layer.