^{*}

In this paper, the problem of steady laminar two-dimensional boundary layer MHD flow and heat transfer of an incompressible viscous fluid with the presence of buoyancy force and viscous dissipation over an vertical nonlinear stretching sheet with partial slip is investigated numerically. Numerical solutions of the resulting nonlinear boundary value problem in the case when the sheet stretches with a velocity varying nonlinearly with the distance is carried out. The effects of for various values of suction parameter, magnetic parameter, Prandtl number, Eckert number, buoyancy parameter, nonlinear stretching parameter and slip parameter on flow and heat transfer characteristics is investigated.

The study of flow over a stretching sheet has generated much interest in recent years in view of its numerous industrial applications such as the aerodynamic extrusion of plastic sheets, the boundary layer along a liquid film, condensation process of metallic plate in a cooling bath and glass, and also in polymer industries. Since the pioneering work of Sakiadis [

The study of two-dimensional boundary layer flow, heat and mass transfer over a porous stretching surface is very important as it finds many practical applications in different areas. To be more specific, it may be pointed out that many metallurgical processes involve the cooling of continuous strips or filaments by drawing them through a quiescent fluid and that in the process of drawing these strips, are sometimes stretched. Viscous dissipation changes the temperature distributions by playing a role like an energy source,which leads to affect heat transfer rates. The merit of the effect of viscous dissipation depends on whether the sheet is being cooled or heated.

Apart from the viscous dissipation, the Joules dissipation also acts as a volumetric heat source. Heat transfer analysis over porous surface is of much practical interest due to its abundant applications. To be more specific, heat-treated materials travelling between a feed roll and wind-uproll or materials manufactured by extrusion, glass-fiber and paper production, cooling of metallic sheets or electronic chips, crystal growing are a few practical applications of flow over a stretching sheet. In all these cases, the final product of desired characteristics depends on the rate of cooling and also the rate of of stretching. In view of all these aspects, the present work deals with the effect of viscous and Joules dissipation on MHD flow, heat and mass transfer over a porous sheet, with partial slip. Researches in these fields have been conducted by many investigators. For example, analytical results were carried out by Vajravelu and Hadjinicolaou [

Very recently, the viscous and joules dissipation and internal heat generation was taken into account in the energy equation. Sajid et al. [

A numerical reinvestigation of MHD boundary layer flow over a heated stretching sheet with variable viscosity has been analyzed by Pantokratoras [

Ishak et al. [

Hossain and Takhar [

The problem of non-linear stretching sheet for different cases of fluid flow has also been analyzed by different researchers. Vajravelu [

The study of magnetohydrodynamics of a conducting fluid finds applications in a variety of astroph-ysical and geophysical problems. The effects of magnetic field on the natural convection heat transfer have been discussed by Romig [

Hence the present study investigates the effect of viscous and Joules dissipation on MHD flow over a porous nonlinear vertical stretching sheet with viscous and joules dissipation with partial slip.

Two-dimensional, nonlinear, steady, MHD laminar boundary layer flow with heat transfer of a viscous, incompressible and electrically conducting fluid over a porous vertical stretching sheet embedded in the presence of transverse magnetic field including viscous and Joules dissipation is considered for investigation. An uniform transverse magnetic field of strength B is applied parallel to y-axis. Consider a stretching sheet that emerges out of a slit at x = 0, y = 0 and subsequently being stretched, as in a polymer extrusion process. Let us assume that the speed at a point on the plate is proportional to the power of its distance from the slit and the boundary layer approximations are applicable. In writing the following equations, it is assumed that the induced magnetic field, the external electric field and the electric field due to the polarization of charges are negligible.

Consider a steady, two-dimensional free convection flow adjacent to a nonlinear stretching vertical sheet immersed in an incompressible electrically conducting viscous fluid of temperature. The stretching velocity and the surface temperature are where a and b are constants with and.

The sheet is assumed to vary nonlinearly with the distance x from the leading edge, i.e. and. Under these conditions, the governing boundary layer equations of momentum, energy with buoyancy, viscous and Joules dissipation, with partial slip are

and are subjected to the following boundary conditions

, , atwhere

where u and v are the velocity components along the x and y axes, respectively. Further, μ, ρ, α, β, T, and g are the dynamic viscosity, fluid density, thermal diffusivity, thermal expansion coefficient, fluid temperature in the boundary layer, and acceleration due to gravity, respectively.

A common feature of all these analyses is the assumption that the flow field obeys the conventional no-slip condition at the sheet that is the velocity component parallel with the sheet becomes equal to the sheet velocity at the sheet. In certain situations, however, the assumption of no-slip does no longer apply and should be replaced by a partial slip boundary condition which relates the fluid velocity u to the shear rate at the boundary. Here L is the slip length, and y denotes the coordinate perpendicular to the surface. This slip-flow condition was first introduced by C-L. M. H. Navier more than a century ago and has more recently been used in studies of fluid flow past permeable walls, slotted plates, rough and coated surfaces, and gas and liquid flow in micro devices. The no-slip boundary condition is known as the central tenets of the Navier-Stokes theory. But there are situations wherein such condition is not appropriate. Especially, no slip condition is inadequate for most non-Newtonian fluids. For example polymer melts often exhibit macroscopic wall slip and that in general is governed by a non-linear and monotone relation between the slip velocity and traction. The fluids exhibiting boundary slip find applications in technology such as in the polishing of artificial heart valves and internal cavities. Navier suggested a slip boundary condition in terms of linear shear stress.

The momentum, and energy Equations (2), (3), and (4) can be transformed into the corresponding nonlinear ordinary differential equations by the following similarity transformation:

where, b is dimensional constant and s is the index of power law variation of temperature.

The transformed nonlinear ordinary differential equations are

Boundary conditions (4) becomes

where, , , and.

The nonlinear boundary value problem represented by Equations (7) to (9) is solved numerically using Fourthorder Runge Kutta shooting technique.

The system of non-linear ordinary differential Equations (7) and (8) together with the boundary conditions Equation (9) are similar and are solved numerically by using the fourth order of Runge Kutta integration scheme accompanied with the Shooting scheme. Making an initial guess for the values of and to initiate the shooting process is very crucial in this process. The success of the procedure depends very much on how good this guess is. Numerical solutions are obtained for several values of the physical parameters i.e. magnetic parameter M, stretching parameter m, Prandtl number Pr, slip parameter Buoyancy parameter, Eckert number (Ec) and, suction/injection parameter.

We have chosen a step size of to satisfy the convergence criterion of 10^{–6} in all cases. The maximum value of was found to each iteration loop by The maximum value of to each group of parameter is determined when the value of the unknown boundary conditions at is not changed to successful loop with error less than 10^{–6}.

In order to gain physical insight, the velocity, and temperature profiles have been discussed by assigning numerical values to the parameter, encountered in the problem i.e. numerical calculations were carried out for different values of suction parameter, magnetic parameter M, power law stretching parameter m, Prandtl number Pr, Eckert number Ec, buoyancy parameter λ, slip parameter γ, and their effects on flow and heat transfer characteristics are analysed graphically.

The influences of the magnetic parameter M on the longitudinal velocity profile is depicted in

The influence of suction parameter (< 0), over the non-dimensional longitudinal velocity profiles are shown in

The influence of injection parameter (> 0), over the dimensionless longitudinal velocity profile is shown in

Figures 4 and 5, describe respectively the behaviors of the longitudinal velocity profile and temperature profile for different values of power law stretching parameter m and it is noticed that increase in m results in decrease of longitudinal velocity profile which is more pronounced for small values of m, where as temperature profile increases with the increase of stretching parameter m. It is observed that the variation of the sheet temperature has a substantial effect on the thermal boundary layer. This effect is more pronounced when sheet temperature varies in the direction of highest stretching rate.

An increase in Prandtl number Pr is associated with a decrease in the temperature distribution which is displayed in

Dimensionless velocity profile is presented in

solutions. In fact, the amount of slip increases monotonically with γ from the no-slip solution for and towards full slip as γ tends to infinity. The latter limiting case implies that the frictional resistance between the viscous fluid and the surface is eliminated, and the stretching of the sheet does no longer impose any motion of the fluid.

In

The influence of Eckert number Ec over dimensionless temperature profiles are shown in

The velocity and temperature profiles presented in Figures 1-9, show that the far field boundary conditions are satisfied asymptotically, which support the validity of the numerical results presented.

u, v: velocity component along the x, y directions

g: acceleration due to gravity

T: temperature of the fluid inside the thermal boundary layer

T_{w}: Wall temperaturetemperature

T_{∞}: Temperature of the fluid at infinity

c_{p}: specific heat at constant pressure

v_{w}: velocity of suction

Pr: Prandtl number

: Local Grashof Number

: Local Reynolds Number

λ: buoyancy parameter

f: dimensionless stream function

m: nonlinear stretching parameter

f_{w}: dimensionless suction/injection parameter

υ: kinematic viscosity

g: Slip Parameter

s : Index of power law variation of temperature.

r: density of the fluid

b: coefficient of volume expansion

: Thermal diffusivity

k: thermal conductivity

η: similarity variable

θ: dimensionless temperature