MHD Forced Convective Laminar Boundary Layer Flow from a Convectively Heated Moving Vertical Plate with Radiation and Transpiration Effect

A two-dimensional steady forced convective flow of a Newtonian fluid past a convectively heated permeable vertically moving plate in the presence of a variable magnetic field and radiation effect has been investigated numerically. The plate moves either in assisting or opposing direction to the free stream. The plate and free stream velocities are considered to be proportional to whilst the magnetic field and mass transfer velocity are taken to be proportional to where is the distance along the plate from the leading edge of the plate. Instead of using existing similarity transformations, we use a linear group of transformations to transform the governing equations into similarity equations with relevant boundary conditions. Numerical solutions of the similarity equations are presented to show the effects of the controlling parameters on the dimensionless velocity, temperature and concentration profiles as well as on the friction factor, rate of heat and mass transfer. It is found that the rate of heat transfer elevates with the mass transfer velocity, convective heat transfer, Prandtl number, velocity ratio and the magnetic field parameters. It is also found that the rate of mass transfer enhances with the mass transfer velocity, velocity ratio, power law index and the Schmidt number, whilst it suppresses with the magnetic field parameter. Our results are compared with the results existing in the open literature. The comparisons are satisfactory.


Introduction
Heat and mass transfer is an important flow transport process which is driven by a combination of both solute and thermal diffusion effects [1]. This type of problem is encountered in fields such as geology, oceanography and astrophysics as well as practical applications such as the design of energy efficient building components, control of the pollutant spread in groundwater and compact heat exchangers [2][3]. Makinde [4] studied similarity solution of MHD flow with heat and mass transfer over a moving vertical plate. He included magnetic field and thermal convective surface boundary condition in his study. In another paper, the same author [5] presented similarity solutions of the hydromagnetic mixed convective flow along a vertical plate in a porous medium which includes radiation and the heat generation effect. Magnetohydrodynamic flows have received the attention of many investigators because of their many engineering applications. In metallurgy, for example, some processes require the cooling of continuous strips by drawing them through an electrically conducting fluid which is subject to a magnetic field. This allows the cooling to be controlled and final product with the requisite properties obtained [6]. Many researchers have investigated MHD flows past various geometries in porous as well as clear media. For example, Pal and Mondal [7] studied the flow and heat transfer effects of an electrically conducting fluids such as liquid metals; water mixed with a little acid and other equivalent substances in the presence of a magnetic field. Herdrich et al. [8] pointed out possible applications of magnetically controlled plasmas in space technology. The similarity solutions for MHD flow with heat transfer over a wedge considering variable viscosity and thermal conductivities effects was investigated by Seddeek et al. [9]. Alam et al. [10] discussed the effects of suction and thermophoresis on steady MHD combined free-forced convective heat and mass transfer flow over an inclined radiative plate. Rahman and Salauddin [11] dealt with the effects of variable electric conductivity and viscosity on hydromagnetic heat and mass transfer flow. Hendi and Hussain [12] presented solution of MHD Falkner-Skan flow over a porous surface by homotopy analysis method. Bararnia et al. [13] obtained analytical solutions of the Falkner-Skan equations with heat transfer using the homotopy perturbation method.
The forced convective flow past a moving plate in a moving fluid has numerous application in engineering including the cooling of polymer sheets and plates and conveyor belts. Thermal transport due to a continuously moving plate through quiescent or moving fluid find applications in processes such as plastic extrusion and continuous casting. Therefore, the knowledge of the flow and thermal fields at the vicinity of the plate is required for ensuring, the quality of the manufactured products [14]. Klemp and Acrivos [15] presented a method for integrating the boundary layer equations through a region of reverse flow. They applied the method to uniform flow past a parallel flat plate. Similar problems have also been studied by Weidman [16], Ishak et al. [17] amongst others.
Heat and mass transfer due to forced convection occurs in chemical process, biochemical process and engineering. Ali [18] investigated the effect of temperature dependent viscosity on laminar mixed convective flow and heat transfer for a continuously moving vertical isothermal surface and obtained local similarity solutions. According to previous studies, the forced convection flow over a moving plate has engineering applications such as liquid films in condensation process and in aerodynamics. Similarity solutions about moving plate were investigated by many authors. Boundary layer flow along a moving plate parallel to a moving stream was studied by Ishak et al. [19]. A similar problem also was studied by Ishak [20] in 2010. Magyari [21] studied the moving plate thermometer with uniform velocity using the Merkin transformation method. Hoernel [22] reported the similarity solutions for steady MHD flow near the forward stagnation point of two-dimensional moving axisymmetric bodies.
The flow of an electrically conducting fluid with thermal radiation occurs for example in electrical power generation and solar power technology. The radiative heat mass transfer also occurs in a variety of geophysical and engineering applications such as migration of moisture through air contained in fibrous insulations, nuclear reactors, nuclear waste disposal and others [23]. The temperature as well as heat transfer may affect radiation, at high operating temperature. Radiation plays an important role in controlling heat transfer in certain processes. Its effects on boundary layer flow with or without magnetic field have been investigated by amongst others: Duwairi and Duwairi [24], Cortell [25], Bataller [26], El-Kabeir et al. [27], Ishak et al. [28], Bakier [29], Cortell [30] etc. All of the above mentioned researchers limited their investigation to either isothermal or isoflux thermal boundary conditions. It seems appropriate to us to generalize the previous investigations to a thermal convective boundary condition instead of an isothermal or isoflux boundary condition.
The present study attempts to pinpoint the combined effect of the transpiration velocity and the radiation on the magnetohydrodynamic forced convective flow of a viscous incompressible Newtonian fluid past a convectively heated moving vertical plate in the moving free stream which have by far not been elucidated in the literature. In an effort to achieve this goal, we shall use a linear group of transformation to develop the similarity transformations (not using the existing transformations) and hence the corresponding similarity representation of the governing equations, before being solved numerically by the Runge-Kutta-Fehlberg fourthfifth order numerical method. The effect of the controlling parameters the velocity ratio, mass transpiration, radiation, magnetic field, power law index, convective heat transfer on the flow, heat and the mass transfer characteristics are investigated and analyzed.

Mathematical Modeling of the Problem
Consider a 2-D MHD forced convective flow along a moving permeable radiating convectively heated vertical flat plate as shown in Fig. 1. (i, ii, iii respectively represent velocity, temperature and concentration boundary layers). It is assumed that an electrically conducting cold fluid having an electric conductivity s in the presence of a magnetic field of variable strength B( x x) is flowing along the flat plate. It is considered that the plate is stretched with a power law velocity and the magnetic field is normal to the plate. Further, the left surface of the plate is assumed to be heated by convection from a hot fluid of temperature T f and this generates a variable heat transfer coefficient h f . The temperature and concentration at the wall are T w and C w . Ambient temperature and concentration are T ? and C ? . Fluid properties are assumed to be invariant. The particle coagulation, magnetic Reynolds number, electric field due to polarization of charges and Hall effects are assumed to be negligible. Under these assumptions the governing transport equations in dimensional form are [22,31].
subject to the boundary conditions Here u u e ( x x): the potential velocity, u u w ( x x): the velocity of the plate, v w ( x x): the velocity normal to the plate with v w ( x x)w0 for injection (blowing), v w ( x x)v0 for suction and v w ( x x)~0 corresponds to an impermeable plate. As usualn~m=r: the kinematic coefficient of viscosity, r: the density of the fluid, s: the electric conductivity, p: the pressure, a~k r c p : the thermal diffusivity, k: the thermal conductivity, D m : the mass diffusivity, s 1 : the Stefan-Boltzmann constant, k 1 : the Rosseland mean absorption coefficient.

Nondimensionalization
Introducing the following dimensionless boundary layer variables with L being the characteristic length, u ? is the reference velocity and Re~u ? L n is the Reynolds number. We assume the following where B 0 , u 0 , v 0 , u ? are constants [22]. Here m (not necessarily an integer) is the power law exponent, u 0 w0 or u 0 v0 mean the plate is moving along positive or along negative direction of x xaxis and u 0~0 represent a stationary plate. Introducing the stream function y so as to reduce the number of the dependent variables as well as the number of the equations. y is defined by into Eqs. (2-5), we have, Ly L y The boundary conditions are Here M is the magnetic field parameter, R is the radiation parameter, Pris the Prandtl number, Sc is the Schmidt number, which are defined as

Scaling Symmetries
The independent and dependent variables are scaled as where A, a i (i~1,2,:::,6)are constant [32][33][34]. Values of a i are sought such that the form of the Eqs.
Solving (15), we have We then seek the functions having the same form before and after the transformation, (i.e., absolute invariants). It is seen from Eqs. (14) and (16) that This combination of the variables is hence an invariant under this group of transformation and so is an absolute invariant. We denote this functional relationshp by g~y Using the same techniques, other absolute invariants are Here g is the similarity independent variable, f (g), h(g), w(g) are similarity dependent variables and h f À Á 0 is the constant heat transfer coefficient.
Substituting Eqs. (18) and (19) into Eqs. (9)-(12), lead to the following ordinary differential equations subject to the boundary conditions where V~u 0 =u ? (constant) is the velocity ratio parameter with V w0 means that motion of the plate is in the same direction as the free stream velocity and V v0 means the plate is moving in the opposite direction of the free stream and V~0 stands for static plate. The case 0vV v1 means the speed of the plate is less than that of the free stream fluid and V w1 mean the speed is greater. V~1 is the case when the plate and the fluid move with the same velocity. Herec~h f À Á 0 L=k ffiffiffiffiffiffi Re p is the convective heat transfer parameter and fw~2v 0 ffiffiffiffiffiffi Re p =u ? (1zm) is the suction/injection parameter, primes denote ordinary differentiation with respect to g. We have successfully transformed a set of partial differential equations into a set of ordinary differential equations with relevant boundary conditions. The solutions of the transformed equations are much easier than the original equations.

Physical Quantities
The quantities we are interested in are the Skin friction factorC f x x (wall shear stress), the local Nusselt number (heat transfer rates) Nu x x and the local Sherwood number (mass transfer rates) Sh x x respectively. They can be derived from the following definitions: where t w q w , q m are the wall shear stress, the wall heat and the wall mass fluxes, and are defined as , q w~{ k L T L y y y y~0 , q m~{ D m LC L y y y y~0 By substituting from Eqs. (6), (18) and (19) into Eq. (25), we get where Re x x~ u u e x x n is the local Reynolds number.
In the absence of energy Eq. (29), we have, subject to the boundary conditions where M 2~2 M (1zm) . Note that Eq. (31) is the MHD Falkner-Skan equation which has been recently investigated by Hendi and Hussain [12]. For stationary impermeable isothermal plate (V~fw~0, c??) and for purely hydrodynamic flow (M~0), Eqs. (28) h''zPr f h'~0: ð34Þ subject to the boundary conditions The parameters band m are related through the expression b~2m=mz1: The parameter b is the Hartee pressure gradient parameter.

Numerical Solution
Similarity Eqs. (20)- (22) with boundary conditions (23) were solved numerically using a fourth-fifth order Runge-Kutta-Fehlberg numerical method. This method has been used in several convective heat and mass transfer problems. The step size Dg~0:001 was used while obtaining numerical solution with g max~1 0 and the accuacy to the sixth decimal place was found sufficient for convergence. The asymptotic boundary conditions given by Eq. (23) were replaced by using a value of 15 for the similarity variable g max as follows.

Results and Discussion
Our numerical results are compared with Aziz [35] in Table 1, Jafar et al. [38] in Table 2, Cebeci and Bradshaw [36] and Yih [40] in Table 3, Hendi and Hossain [12] and Abbasbandy and Hayat [39] in Table 4, and Barania et al. [13] and Rajagopal et al. [42] in Table 5, Hendi and Hossain [12] and White [41] in Table 6 to justify the accuracy of the method used. The comparisons are satisfactory. Figures 2, 3, 4, 5, 6, 7 exhibit the emerging parameters effect on the flow and heat and the mass transfer characteristic. Fig. 2 (a) depicts the effect of the magnetic field and velocity ratio parameters whilst Fig. 2 (b) exhibits the effect of the power law index and the transpiration velocity parameters on the dimensionless velocity. The magnetic field leads to decrease the velocity and then the boundary layer thickness shrinks when the plate moves faster than the free stream in the same direction (V~2). This is due to the fact that applied magnetic force act like a drag force which reduces the velocity. It is interesting to note that in the case of the stationary plate (V~0) the magnetic field causes to elevate not retard the velocity (Fig. 2 (a)) (as in classical MHD boundary layer flow. Physically this is because for stationary plate in the moving free stream, application of the magnetic field which is moving with the free stream has the tendency to induce the motive force which induces the motion of the fluid and lessens the boundary layer thickness. Clearly, an increase in the velocity  It is found from Fig. 2 (b) that the velocity suppresses with the rising of the nonlinear power law parameter m~0 for both permeable (fw=0) and solid (fw~0) convectively heated plates. This figure is drawn by taking Pr~Sc~1which means that the momentum, heat and mass diffuse at the same rate. This result is consistent with the know results in the literature. It is also found that the velocity enhances with the suction and suppresses with the injection for both constant free stream (m~0) and quadratically varying free stream (m~2). An illustration for this trend is that the fluid is brought nearer to the plate in the case of suction (fww0), and the reverse behavior is noticed for the case of injection (fwv0).Note that this figure is drawn for the plate which moves in the same direction of free stream with velocity faster than free stream.
The influences of the radiation, mass transfer velocity, magnetic filed and the velocity ratio parameters on the dimensionless temperature profiles are shown in Fig. 3. The temperature reduces with the injection and rises with the suction both for the radiating and non-radiating moving plate. The thermal radiation is found to increase the temperature both for permeable and solid plate (Fig. 3  (a)). This is because the inclusion of the thermal radiation increases the thermal diffusivity (Eq. 2) which then increases the temperature and hence reduces the heat transfer rate. This is what we observe from Fig. 6a. It is noticed from Fig. 3 (b) that the velocity ratio leads to reduce the dimensionless temperature for both the hydrodynamic and magnetohydrodynamic boundary layer flows. This phenomenon occurred when the plate and free stream move  in the like direction with the plate velocity faster than the free stream velocity or for stationary plate. Figure. 3 (b) reveals that the dimensionless temperature reduces with the elevating of the magnetic field parameter in case of stationary (V~0) permeable (fw~0) flat plate for uniform free stream (m~0). Physically this is because magnetic field enhances the flow velocity, thermal energy is replaced by kinetic energy and this causes reduction of the fluid temperature through the boundary layer. It is further revealed that the dimensionless temperature enhances with the elevating of the magnetic field parameter for the plate which moves in the like direction of the free stream with higher velocity than free stream (V~2,4). Physically this is due to the fact that as Mincreases the Lorentz force opposes the flow which increases the friction between the fluid layers and hence enhances the temperature.
It is noticed from Fig. 4 (a), that both the Schmidt number and the velocity ratio parameters reduce the concentration both for stationary and moving plates. The plate is moving either in the same or in the opposite direction to the free stream velocity. This is because as the Schmidt number Sc increases due to a reduction in the molecular diffusivity of chemical species for a fixed value of kinematic viscosity, the thickness of the concentration boundary layer shrinks. Figure 4 (b) shows the variations of the nondimensional concentration with the suction/injection and magnetic field parameters. Like, the temperature, concentration is found to be decreased with the magnetic field parameter for stationary (permeable/solid) plate. As usual, suction increases whilst injection decreases the concentration. Figure 5 shows the effect of the magnetic field, power law index and suction/injection parameters on the dimensionless friction factor when the plate is either stationary or moves along positive xaxis. An inspection of this Fig. reveals that friction factor increases with magnetic filed, index parameter as well as with mass transfer velocity. This trend is true for both the stationary and moving plates. A close observation of Fig. 5 a and 5 b show that friction factor for stationary plate is higher than the moving plate.
The effect of the convective heat transfer, radiation, mass transfer velocity, velocity ratio, magnetic field and Prandtl number on the dimensionless heat transfer rates is exhibited in Fig. 5. It is noticed that the rate of heat transfer increases with the convective heat transfer parameter for both the permeability and solid plate ( Fig. 6 (a)). It is further noticed that the heat transfer rate rises with the transpiration velocity. The thermal radiation leads to a reduction of the heat transfer rate for both permeable and impermeable convectively heated plates. Figure 6 (b) reveals that the rate of heat transfer is a monotonic increasing function of the  velocity ratio and Prandtl number for both the magnetohydrodynamic and purely hydrodynamic boundary layer flow. It is further seen that the heat transfer rate is an increasing function of the magnetic field parameter for both stationary and moving permeable radiating plate. Finally, Fig. 7 exhibits the effect of the Schmidt number, mass transfer velocity, magnetic field, velocity ratio and power law index parameters on the dimensionless mass transfer rates. From Fig. 7 (a), it is observed that dimensionless the mass transfer rate enhances with Schmidt number for both permeable and solid plate in case of both magnetohydrodynamic and purely hydrodynamic boundary layer flow. The rate of mass transfer is suppressed with the increasing value of the magnetic field whilst it increases for increasing of the velocity ratio and power law index parameters.

Conclusions
In this paper we studied the problem of heat and mass transfer along a moving permeable radiating vertical flat plate in the moving free stream with the convective surface boundary condition numerically. Similarity transformations have been developed by a linear group of transformation to transform the transport equations into a corresponding system of nonlinear ordinary differential equations. These equations are then solved numerically using the Runge-Kutta-Fehlberg fourth-fifth order numerical method under Maple 13. The main results can be summarized as below: N The velocity decreases with pressure gradient, mass transfer velocity and the magnetic field parameter.
N Thermal radiation leads to increase temperature and decrease heat transfer rates.
N The temperature, concentration and mass transfer rates decrease with the magnetic field.
N The rate of heat transfer elevates with mass transfer velocity, convective heat transfer, Prandtl number, velocity ratio and magnetic field parameters.
N The rate of mass transfer enhances with the mass transfer velocity, velocity ratio, power law index and the Schmidt number, whilst it falls as the magnetic field increases.