MHD Natural Convection with Convective Surface Boundary Condition over a Flat Plate

and Applied Analysis 3 the thermal conductivity, and D is the molecular diffusivity. The respective boundary conditions are V = 0, u = N


Introduction
A review of the literature shows that to the best of our knowledge not much research has been reported on MHD flow over a flat plate with convective surface boundary conditions by applying the one parameter continuous group method.For this problem we apply similarity transformations on the partial differential equations.The transformed nonlinear coupled ordinary differential equations are solved numerically by the Keller Box method for different values of controlling parameters.
Analysis of natural phenomena usually leads to partial differential equations and nonlinear ordinary differential equations.Nonlinear differential equations appear in physics, applied mathematics, and engineering sciences.In most cases for these problems exact solutions cannot be obtained.One of the most widely used applications of nonlinear differential equations is boundary-layer problems.Fluid flow and heat transfer are a relevant problem in many industrial processes such as metal and polymer extrusion processes, glass-fiber and paper production, manufacture and drawing of plastics and rubber sheets, and crystal growing.Magnetohydrodynamics (MHD) is the flow of an electrically conducting fluid in the presence of a magnetic field.This effect is of importance in various areas of technology and engineering such as MHD flow meters, MHD power generation, and MHD pumps [1][2][3][4].The study of the interaction of conducting fluids with electromagnetic phenomena is important and such problems have received much attention from many researchers.Mukhopadhyay et al. [5], Andersson [6], Rashidi et al. [7], and Parsa et al. [8] investigated the effect of magnetic field over a stretching surface in various states.Numerical results for MHD free convection flow over a wedge in the presence of a magnetic field were presented by Watanabe and Pop [9].Kumari and Nath [10] studied unsteady MHD viscous flow and heat transfer of Newtonian fluids induced by an impulsively stretched plane surface in two lateral directions by employing the homotopy analysis method.Rashidi et al. [1] solved the governing equations of suction and injection effects on the free convection boundary-layer flow over a vertical cylinder.In addition, a complete investigation of MHD studies and their technological applications was undertaken by Moreau [11].Several interesting computational studies of reactive MHD boundary-layer flows with heat and mass transfer have appeared in recent years [12][13][14][15].Effects of anisotropic scattering on steady nonsimilar free convective radiative hydromagnetic boundary-layer flow over a diffuse reflecting surface and solution of a separate equation for the magnetic field distribution were presented by Chen [16].Ishak [17] studied steady laminar boundary-layer flow and heat transfer over a stationary permeable flat plate immersed in a uniform free stream with convective boundary condition.The problem of a vertical plate with convective boundary conditions was considered by Makinde [18].Rashidi et al. [19] presented the first and second law analyses of an electrically conducting fluid past a rotating disk in the presence of a uniform vertical magnetic field by using the homotopy analysis method (HAM) and then applied artificial neural networks (ANN) and the particle swarm optimization (PSO) algorithm to minimize the entropy generation.
The main objective in this paper is to investigate similarity solutions and scaling transformations of MHD heat and mass transfer flow of a steady viscous incompressible fluid over a flat plate with convective surface boundary conditions by using the one parameter continuous group method.A convective boundary condition instead of the commonly used constant surface temperature or constant heat flux boundary conditions is applied.The governing boundarylayer equations are transformed to a two-point boundary value problem in similarity variables, and the problem is solved numerically by the Keller Box method.The effects of governing parameters on fluid velocity, temperature, and particle concentration are investigated and shown graphically.

Mathematical Formulation of the Problem
The problem of two-dimensional steady MHD heat and mass transfer laminar flow of a viscous incompressible and electrically conducting fluid past a flat plate is considered.The  axis is taken along the plate and the  axis is normal to the plate.The gravitation acceleration vector is parallel to plate.A magnetic field of uniform strength  0 is applied perpendicular to the direction of the plate.The viscosity, thermal conductivity, and concentration diffusivity of fluid are assumed to vary linearly.The top surface of the plate is kept at uniform temperature   which is assumed to be greater than the full stream temperature  ∞ .The species concentration   at the surface is uniform and the full stream concentration is  ∞ .The bottom surface of the plate is heated by convection from a hot fluid of temperature   which provides a heat transfer coefficient ℎ  .The induced magnetic field due to the motion of the electrically conducting fluid is negligible.This assumption is valid for small magnetic Reynolds numbers.It is also assumed that the external electric field is zero and the electric field due to polarization of charges is negligible.It is also assumed that the pressure gradient and viscous and electrical dissipation are neglected.
The physical configuration and schematic of the problem are shown in Figure 1.It is known that this is a type of Falkner-Skan flow.Furthermore the following assumptions are considered: (i) fluid has constant kinematic viscosity and thermal diffusivity and the Boussinesq approximation may be adopted for steady laminar flow, (ii) the particle diffusivity is constant, (iii) the concentration of particles is sufficiently dilute that particle coagulation in the boundary layer is negligible, and (iv) the magnetic Reynolds number is small so that the induced magnetic field is negligible in comparison with the applied magnetic field.Under these assumptions the governing Prandtl boundary-layer equations in dimensional form are as follows (see Kays et al. [20] and White [21]): where  and V are the velocities in the  and  directions, respectively,  is the temperature within the boundary layer,  ∞ is the temperature far away from the plate,  is the species concentration,  is the acceleration due to gravity,   is the volumetric coefficient of thermal expansion,   is the volumetric coefficient of concentration expansion,  is the thermal conductivity, and  is the molecular diffusivity.
The respective boundary conditions are where   is the velocity over the plate that should be in the form   =   .This condition will be imposed later.(), (), and  () are variable viscosity, thermal conductivity, and molecular diffusivity, respectively; the dimensions of  1 are (velocity) −1 and the dimension of  1 is length.It is assumed that the temperature dependent viscosity and thermal conductivity vary linearly and are given by (see Seddeek and Salem [22]) where  ∞ and  ∞ are the constant undisturbed viscosity and undisturbed thermal conductivity,  1 is a constant with  1 > 0, and  is a constant which depends on the fluid.It is assumed that the concentration diffusivity varies linearly and is given by (see Seddeek and Salem [22]) where   is the constant concentration diffusivity.
The following dimensionless variables are introduced: where Re is the Reynolds number,  is the characteristic length,  is the dimensionless temperature variable, and  is the dimensionless concentration variable.Introducing the stream function  such that  = / and V = −/, continuity equation ( 1) is satisfied identically and ( 2)-( 4) now yield The boundary conditions are In the above equations the parameters are defined as where Re is the Reynolds number, Sc is the Schmidt number,  is the magnetic parameter, Pr is the Prandtl number of the fluid,  is the thermal conductivity parameter,  is the viscosity parameter,   is the concentration diffusivity parameter,  is the velocity slip parameter, and  is the thermal slip parameter.

Application of Group Transformations
Determining similarity solutions of ( 9)-( 10) is equivalent to determining invariant solutions of these equations under a particular continuous one parameter group (Hamad et al. [23] and Kandasamy et al. [24]).Thus we search for a transformation group from the elementary set of oneparameter scaling transformations as one of the techniques that are defined by the following group which is called G1: Here ( ̸ = 0) is a parameter of the group and the 's are arbitrary real numbers whose connection will be determined by our analysis.The transformations listed in (12) may be treated as point transformations which transform the coordinates (, , , , ,   ,   ,   ) to ( * ,  * ,  * ,  * ,  * ,  *  ,  *  ,  *  ) .
The system ( 9)-( 10) remains invariant under the group transformation G1.Hence we have the following relationships among the parameters, namely, From boundary conditions (10), these will be invariant if Solving ( 14) and ( 15), we obtain With these relations the boundary conditions remain invariant.
The set of transformations G1 in (12) then reduces to Using a Taylor series expansion in powers of , retaining terms up to first order, and neglecting higher powers of  results in The characteristic equations are Solving the above characteristic equations gives Substituting ( 20) into ( 9)-( 10) yields Here Gr = (  −  ∞ )  0 / 2 ∞ and Gc = (  −  ∞ )  0 /  2 ∞ are Grashof numbers based on temperature and on concentration, respectively.
The corresponding boundary conditions are To obtain a similarity solution for the energy equation, the quantity  must be independent of  and for this to occur the heat transfer coefficient ℎ  must be directly proportional to () −1/2 .

Parameters of Physical Interest.
We are interested in the friction factor   , Nusselt number Nu, and Sherwood number Sh, respectively.Physically,   indicates wall shear stress and Nu indicates the rate of heat transfer whilst Sh indicates the rate of mass transfer.These quantities may be conveniently determined from By substituting (10) and ( 19) into (23), we obtain From (24) it can be shown that the skin friction factor   , the Nusselt number Nu, and the Sherwood number Sh are proportional to the numerical values   (0), −  (0), and −  (0), respectively.

The Keller Box Method
Equation ( 21) subject to boundary conditions ( 22) is solved numerically using a very efficient finite difference scheme known as the Keller Box method.The details of this method are described in Cebeci and Bradshaw [25] and Na [26].For more information refer to Keller [27,28].

Results and Discussion
Applying scaling group transformations to analyze the governing equations and the boundary conditions, the two independent variables are reduced by one.Consequently the governing equations reduce to a system of nonlinear ordinary differential equations with the appropriate boundary conditions.The transformed momentum, energy, and concentration equation ( 21) subject to the boundary conditions (22) were solved numerically by using the Keller Box method.We obtained velocity, temperature, and concentration profile graphs for different values of governing parameters.
Figures 2, 3, and 4 show the effects of the viscosity parameter and the magnetic parameter on the velocity, temperature, and concentration distributions, respectively.The velocity distribution decreases with increasing  and  whereas they have no significant effect on the temperature and concentration distributions.This behavior can be predicted from (21) and also the physical definition of parameters  and , since the viscosity and magnetic parameters only appear in the momentum equation.Figures 5, 6, and 7 show the effect of the thermal and mass Grashof numbers on the velocity, temperature, and concentration distributions, respectively.Physically, since the thermal Grashof number (Gr) is the ratio of buoyancy to viscous forces in the boundary layer, increasing its value suggests an increase in the buoyancy forces relative to the viscous forces and this is clearly reflected in the progressive increase in the velocity of the flow.Increase in the mass transfer Grashof number (Gc) yields a similar effect on the velocity of the flow.Moreover, the reverse trend is seen for the temperature and concentration distributions.Figures 8, 9, and 10 illustrate the influence of the thermal conductivity parameter  and Prandtl number Pr on the velocity, temperature, and concentration distributions, respectively.It is observed that the velocity and temperature distributions increase with increasing thermal conductivity parameter and decrease with increasing Prandtl number.This is in agreement physically since the thermal boundarylayer thickness decreases with increasing Pr.The thermal conductivity parameter  and the Prandtl number Pr have no significant effect on the concentration distribution and this can be predicted from (21).The effects of the concentration diffusivity parameter   and the Schmidt number Sc on the velocity, temperature, and concentration distributions are shown in Figures 11-13.The velocity and concentration distributions increase with increasing concentration diffusivity parameter whereas they decrease with increasing Schmidt number.Since Schmidt number is the ratio of viscosity to diffusivity, this behavior can be predicted.From Figure 12, the concentration diffusivity parameter and the Schmidt number have no significant effect on the temperature distribution.In Figures 14-16 effects of the velocity slip parameter () and the thermal slip parameter () are depicted.In Figure 14 it is observed that velocity distribution increases with increasing velocity slip parameter and decreases with increasing thermal slip parameter.From Figure 15 we observe that as the velocity slip parameter and the thermal slip parameter increase the temperature distribution decreases.
Figure 16 shows that the concentration distribution decreases with increasing velocity slip parameter and increases with increasing thermal slip parameter.In some of the velocity profiles an overshoot of the velocity profile is observed.This depends on the boundary conditions.In other words, since there is not a no-slip condition on the plate, a larger velocity rather than free stream velocity can exist inside the boundary layer.With attention to boundary conditions (5) an overshoot of the velocity is likely.

Conclusions
A numerical study based on the Keller Box method for MHD heat and mass transfer flow of a steady viscous Table 1: Numerical results of   (0),   (0), and   (0) for different values of the parameters Gr, Gc, , and Pr when  = 1.0,  = 1.0,   = 1.0,Sc = 1.0;  = 1.0, and  = 1.0.(1) The thickness of the velocity boundary layer decreases with an increase in viscosity parameter , magnetic field , Schmidt number Sc, and thermal slip parameter .

Gr
(2) The thickness of the velocity boundary layer increases with an increase in thermal Grashof number Gr, mass transfer Grashof number Gc, thermal conductivity parameter , concentration diffusivity parameter   , and velocity slip parameter .
(3) The thickness of the thermal boundary layer decreases with an increase in thermal Grashof number Gr, mass transfer Grashof number Gc, Prandtl number Pr, velocity slip parameter , and thermal slip parameter .
(4) The thickness of the thermal boundary layer increases with an increase in thermal conductivity parameter .
(5) The thickness of the concentration boundary layer decreases with an increase in thermal Grashof number Gr, mass transfer Grashof number Gc, Schmidt number Sc, and velocity slip parameter .
(6) The thickness of the concentration boundary layer increases with an increase in concentration diffusivity parameter   and thermal slip parameter .

Figure 1 :
Figure 1: Physical configuration and the coordinate system.

Figure 2 :
Figure 2: Effects of the viscosity parameter  and the magnetic parameter  on the dimensionless velocity.

Figure 3 :
Figure 3: Effects of the viscosity parameter  and the magnetic parameter  on the dimensionless temperature.

Figure 4 :
Figure 4: Effects of the viscosity parameter  and the magnetic parameter  on the dimensionless concentration.

1 Gr = 1 ;Figure 7 :Figure 8 :
Figure 7: Effects of the thermal Grashof number Gr and the mass Grashof number Gc on the dimensionless concentration.

1 ;Figure 9 :Figure 10 :
Figure 9: Effects of the thermal conductivity parameter  and the Prandtl number Pr on the dimensionless temperature.

Figure 11 : 5 ; 1 D c = 1 ;Figure 12 :
Figure 11: Effects of the concentration diffusivity parameter   and the Schmidt number Sc on the dimensionless velocity.

Figure 13 :Figure 14 :
Figure 13: Effects of the concentration diffusivity parameter   and the Schmidt number Sc on the dimensionless concentration.

Figure 15 :Figure 16 :
Figure 15: Effects of the velocity slip parameter  and the thermal slip parameter  on the dimensionless temperature.
a flat plate has been performed.We have investigated the effects of various governing parameters, namely, the viscosity parameter , the magnetic field , thermal Grashof number Gr, mass transfer Grashof number Gc, thermal conductivity parameter , Prandtl number Pr, concentration diffusivity parameter   , Schmidt number Sc, velocity slip parameter , and thermal slip parameter  on flow field and heat transfer characteristics.The following conclusions can be made.