The non-Newtonian fluid flow is important for their numerous engineering applications. The various types of non-Newtonian fluids are nano fluid, Casson fluid, viscoelastic fluid, couple stress fluid, micro polar fluid, power-law flow, etc. These include pseudo plastic, dilatant, blood, foodstuff, slurries, cosmetics & toiletries and paints. Heat transfer in boundary layer over a stretching sheet has important applications in extrusion of plastic sheets, polymer, spinning of fibers, cooling of elastic sheets, etc. The quality of final product depends on the rate of heat transfer, and as a result the cooling procedure has to be controlled effectively. The MHD boundary layer flow of heat and mass transfer problems about an stretching sheet have become in view of its significant applications in industrial manufacturing processes such as plasma studies, petroleum industries, magneto-hydrody- namics power generator, cooling of nuclear reactors, boundary layer control in aerodynamics, glass fiber production and paper production. In recent years, due to the growing applications of non-Newtonian fluids, various researchers have done different works in this field. Astarita and Marrucci [1] and Bhome [2] have studied the behavior of non-Newtonian fluids in case of steady and unsteady flow situations. At first, Schowalter [3] [4] formulated the boundary layer flow of a Non-Newtonian fluid and also established the existence conditions of a similarity solution. A similarity solution has been done by Acrivos et al. [5] for a power-law fluid flowing along a vertical flat plate. Heat transfer of a non-Newtonian power-law fluid over a stretching surface in the presence of radiation and slip condition at the surface was investigated by Mostafa [6] . Khan and Pop [7] have studied the fundamental work on the boundary layer flow of nano fluid over a stretching sheet. Rana and Bhargava [8] performed the numerical study of heat transfer characteristic in the mixed convection flow of a nano fluid along a vertical plate in presence of heat source/sink. Rama and Goyal [9] have studied non-Newtonian nano fluid flow over a permeable sheet with heat generation and velocity slip in presence of magnetic field. Noghrehabadi et al. [10] studied the partial slip eeffect on the flow and heat transfer of nano fluids over a stretching sheet.

In this regard the aim of this work is to investigate the steady boundary layer flow and combined heat and mass transfer of a non-Newtonian fluid over an inclined stretching sheet in presence of uniform magnetic field. The numerical results are carried out for several values of the combined effects of magnetic parameter M, stretching parameter λ, Prandtl number Pr, Eckert number Ec, Schmidt number Sc, Soret number S₀, slip parameter A and Casson parameter n on velocity, temperature and concentration profiles and also the skin-friction coefficient, local Nusselt number and local Sherwood number are discussed by considering cooling plate.

Let us consider a two dimensional steady laminar viscous incompressible MHD flow of an electrically conducting fluid along an inclined stretching sheet with an acute angle (α), X-direction is taken along the leading edge of the inclined stretching sheet and Y is normal to it and extends parallel to X-axis. A magnetic field of strength B₀ is introduced to the normal to the direction to the flow. The uniform plate temperature T_w (>T_∞), where T_∞ is the temperature of the fluid far away from the plate. Let u and v be the velocity components along the x and y axis respectively in the boundary layer region. Under the above assumptions and usual boundary layer approximation, the dimensional governing equations of continuity, momentum, concentration and energy under the influence of externally imposed magnetic field [11] are:

Equation of continuity:

Momentum equation:

Energy equation:

Concentration equation:

Using free stream velocity

The Equation (2) can be written as

where u and v are the velocity components along x and y directions, T, T_w and are the fluid temperature, the stretching sheet temperature and the free stream temperature respectively while C, and are the corresponding concentrations, κ is the variable thermal conductivity of the fluid, specific heat with constant pressure, α is the angle of inclination, μ is the coefficient of viscosity, ν is the kinematic viscosity, σ is the electrical conductivity, ρ is the fluid density, β is the thermal expansion coefficient, β^* is the concentration expansion coefficient, is the magnetic field intensity, U is the free steam velocity, g is the acceleration due to gravity, D_m is the coefficient of mass diffusivity, T_m is the mean fluid temperature, K_T is the thermal diffusion ratio respectively and n is known as non-dimensional power-law index (Casson parameter) and it is an important index to sub-classified fluids like pseudo plastic fluids (n < 1) and dilatant fluids (n > 1). Also for n = 1, the fluid is simply the Newtonian fluid. The deviation of n from unity indicates that the degree of deviation from Newtonian behavior [12] .

The above equations are subject to the following boundary conditions:

To convert the governing equations into a set of similarity equations, we introduce the following similarity transformation:

The steam function as defined by

From the above transformations, the non-dimensional, nonlinear and coupled ordinary differential equations are obtained

The transform boundary conditions:

where, and are the dimensionless velocity, temperature and concentration respectively, is the similarity variable, A is the slip parameter, the prime denotes differentiation with respect to. Also

are the magnetic parameter, stretching ratio, Grashof number, modified Grashof number, Prandtl number, Sch- midt number, Eckert number, slip parameter and Soret number respectively. The important physical quantities of this problem are skin friction coefficient, the local Nusselt number and the local Sherwood number Sh which are proportional to rate of velocity, rate of temperature and rate of mass transfer respectively.

The governing partial differential equations are transformed into ordinary differential equations by using similarity transformation and stretching variable. The flow is considered under the influence of a stretching velocity and a uniform magnetic field. The governing thermal boundary layer Equation (3), concentration boundary layer Equation (4) and momentum boundary layer Equation (5) with the boundary conditions Equation (6) are transformed into a system of ordinary differential equations by the suitable local similarity transformations which are then solved numerically by using Runge-Kutta fourth-fifth order method along with shooting iteration technique. The higher order nonlinear differential equations are converted into simultaneous linear differential equations of first order and they are further transformed into initial value problem applying the shooting technique. First of all, the coupled ordinary differential Equations (7)-(9) are third order in f and second order in and φ which have been reduced to a system of seven simultaneous ordinary differential equations for seven unknowns. For the purpose of numerically solve this system of equations using Runge-Kutta method, the solution requires seven initial conditions but two initial conditions in f, one initial condition in each of and φ are known. However, the values of, and are known at. These end conditions are utilized to produce unknown initial conditions at by using shooting technique. The most important step of this scheme is to choose the appropriate finite value of.

Thus to estimate the value of, we start with some initial guess value and solve the boundary value problem consisting of Equations (7)-(9) to obtain. The solution process is repeated with another larger value of until two successive values of differ only after desired significant digit. The last value is taken as the finite value of the limit for the particular set of physical parameters for determining velocity, temperature and concentration, respectively, are in the boundary layer. After getting all the initial conditions we solve this system of simultaneous equations using fourth order Runge-Kutta integration scheme. The effects of the flow parameters on the velocity, temperature, species concentration, shearing stress, rate of heat transfer, and rate of concentration are computed, discussed and have been graphically represented in figures and tables for various value of different parameters. Now defining new variables by the equations. The higher order differential Equations (7)-(9) and boundary conditions Equation (10) may be transformed to seven equivalent first order differential equations and boundary conditions respectively are given below:

And the boundary conditions are

Numerical calculation for distribution of the velocity, temperature and concentration profiles across the boundary layer for different values of the parameters are carried out in presence of cooling of the plate (Gr > 0). For the purpose of our simulation we have chosen λ = 0.85, M = 0.3, n = 0.9, n = 1.0, n = 1.9, A = 0.5, Gr = 2.0, Gm = 2.0, Sc = 0.22, Pr = 1.0, Ec = 1.0, S₀ = 0.4, and α = 30˚ while the parameters are varied over range as shown in the figures. Figures 1-26 depict the effect of different entering parameters for (a) Pseudo plastic fluid (n < 1) (b) Dilatants fluid (n > 1) (c) Newtonian fluid for (n = 1). From Figures 1-3 it is observed that the velocity profile decreases for both Newtonian and non-Newtonian fluids with increase in magnetic parameter, M because the presence of M in an electrically conducting fluid introduces a force like Lorentz force which acts against the flow if the magnetic field is applied in the normal direction as in the present problem. As a result velocity profile is decreased. Figures 4-6 displayed the effect of stretching parameter on velocity profile. From these figures it is shown that the velocity is increased for both Newtonian and non-Newtonian fluids. The effects of Casson parameter are shown in Figure 7 and Figure 8. It is observed that the velocity is decreased for Dilatant fluids but increased in case of Pseudo plastic fluids. Therefore, it is observed that the momentum boundary layer thickness for Pseudo plastic fluid (n < 1) is greater than that of Dilatants fluid (n > 1). It happens because of plasticity of Casson fluid. The plasticity of the fluid increases while Casson parameter decreases, which causes the increment in velocity boundary layer thickness. From Figures 9-11 it is observed that the temperature profile increases for the increasing values of magnetic parameter both the case of Newtonian and non-Newtonian fluids which means that the applied magnetic field in an electrical conducting fluid is manifested as thermal energy. This heats conducting Newtonian and non-Newtonian fluids and hence increases the temperature. Besides the fluid velocity is decreased for magnetic parameter which leads to a decrease in heat convection and mass diffusion, and thus fluid temperature and species concentration in fluid increases. From Figures 12-14 it is observed that the temperature decreases up to η < 1 and increase 1 ≤ η < 5 and then decreases for the increasing values of Prandtl number in case of Newtonian and non-Newtonian fluids. Since Prandtl number is a property of fluids that is a fluid with small Prandtl number, in general, has high thermal diffusivity. Due to this property the fluid temperature is higher for a lower Prandtl number fluid and vice-versa. The increase in fluid temperature in turn increases the buoyancy force and so the fluid velocity increases. The increase in fluid velocity leads to better species diffusion; therefore, the species concentration decreases in the fluid, which is seen in Figure 26. Again, Figures 15-17 depict the temperature profile for increasing values of Eckert number. From these figures it is observed that the temperature profile increase for the effect of Eckert number because heat energy is stored within the fluid due to the frictional heating as a result enhanced the temperature. From the Figure 18-20 it is seen that the concentration profile is increased for increasing values of Soret number in case of both the Newtonian and non-Newtonian fluids. A similar kind of interplay as above for Pr, occurs for Schmidt number Sc and is seen in Figures 21-23. Schmidt number is a property of fluids that is a smaller value of Sc means better diffusion of species in fluid. Hence, for a small Schmidt number the species concentration in the fluid is higher and vice-versa for both Newtonian and non-Newtonian fluids.

Similar result has been found for the increasing effect of magnetic parameter and Eckert number in case of Dilatants fluids which are shown in Figure 24 and Figure 25. Tables 1-3 exhibit the behavior of, , and, for various values of magnetic parameter, stretching parameter, Prandtl number, Eckert number, Soret number, Casson parameter and Schmidt number. Here measures skin friction coefficient, that is, drag observed at the surface, whereas and are measures of heat flux (rate of heat transfer) and mass flux (rate of mass transfer) at the surface, respectively. Also, when skin friction is positive, it means that the fluid exerts drag on the sheet; when negative, it means that the sheet exerts drag on the fluid. From Table 1, it is observed that is decreased for increasing values of M and n and increased for increasing values of. Therefore, with the increase in the magnetic parameter M, the skin friction decreases. As discussed above, the fluid velocity decreases due to an increase in magnetic parameter. For M = 0.30, 0.35, 0.40 the sheet exerts drag on the fluid. Further, the decrease in fluid velocity with M increases fluid temperature; therefore the rate of heat transfer decrease with an increase in M which is shown in Table 2. From Table 2, it is observed that, is increased with the increasing values of M and Pr whereas reverse trend arises for Ec. Since the rate of heat transfer is negative which means that the heat is transferred from fluid to sheet. From Table 3, we see that, is decreased with the increase of M and increased for Sc, Ec and Pr.

The followings are conclusions made from the above analysis:

・ The magnitude of velocity decreases with increasing magnetic parameter of Lorentz force for both Newtonian and non-Newtonian fluids but increases for the increasing values of stretching parameter.

・ The thermal boundary layer increases for increasing values of magnetic parameter in both the cases of Newtonian and non-Newtonian fluids, meaning that the applied magnetic field in an electrical conducting fluid is manifested as thermal energy. Besides the fluid velocity decreases for magnetic parameter which leads to a decrease in heat convection and mass diffusion, and thus fluid temperature and species concentration in fluid increase. Similar results arises for the effect of Prandtl number Pr. Prandtl number is a property of fluids with small Prandtl number, and in general, it has high thermal diffusivity. Due to this property the fluid temperature is higher for a lower Prandtl number fluid and vice-versa. The increase in fluid temperature in turn increases the buoyancy force, so the fluid velocity increases. The increase in fluid velocity leads to better species diffusion; therefore, the species concentration decreases in the fluid, which is seen in Figure 26.

・ A similar kind of interplay as above for Pr occurs for Schmidt number Sc and is seen in Figures 21-23. Schmidt number is a property of fluids that is a smaller value of Sc, meaning better diffusion of species in fluid. Hence, for a small Schmidt number, the species concentration in the fluid is higher and vice versa for both Newtonian and non-Newtonian fluids. Similar result has been found for the increasing effect of magnetic parameter and Eckert number in case of Dilatants fluids which are shown in Figure 24 and Figure 25.