The steady magnetohydrodynamic (MHD) flow and mass transfer of an incompressible, viscous, and electrically conducting fluid over a permeable flat surface stretched with nonlinear (quadratic) velocity and appropriate wall transpiration is investigated. It is shown that the problem permits an analytical solution for the complete set of equations with magnetic field influences when a fictitious presence of a chemical reaction is considered. Velocity and concentration fields are presented through graphs and discussed. The results for both skin friction coefficient and mass transfer gradient agree well with numerical results published in the literature.

1. Introduction

Momentum and mass transport generated by a linearly or nonlinearly stretching surface has extensively been treated over the last few decades due to its relevance to a wide range of industrial applications. In particular, processes involving the mass transfer effect have been recognized as important principally in chemical processing equipment [1]. The heat and mass transfer problem with or without magnetic field, suction or injection, has been analyzed by a great deal of researchers since it has many engineering applications such as production of polymer sheets, crystalline materials and glass sheets, paper and textile industries, and many others. Sakiadis [2] initiated the study of the boundary layer flow over a continuous solid surface moving with constant speed. Crane [3] was the first one who studied the stretching problem taking into account the fluid flow over a linearly stretched surface. The continuous elastic sheet moves in its own plane with a velocity varying linearly with the distance from a fixed point. Later, P. S. Gupta and A. S. Gupta [4] extended Crane’s work [3] to include suction or blowing, and Vleggaar [5] analyzed the stretching problem with constant surface temperature. After these pioneering works, analyses to the flow and heat/mass transfer over a stretching surface have intensely been made in recent years [1, 611]. Moreover, the author of the present study has also analyzed thermal radiation, suction, and viscous dissipation effects on the flow and heat transfer of a power-law fluid over a porous surface [12]. Almost all the aforementioned researches were performed taking into account either a constant value for the velocity wall or a linearly stretching sheet problem (i.e., ).

Nevertheless, it needs to be underlined that the quadratic stretching of the sheet may not necessarily be linear. In view of this, a number of studies have appeared in the literature to analyze various different nonlinearly stretching rates; a quadratic stretching function was considered in the early investigation of Takhar et al. [13]. Later, Kumaran and Ramanaiah [14] showed that a quadratic term of stretching could be superposed on the linearly stretching Crane’s [3] result by considering a suitable distribution of transpiration along the plate stretched in a nonlinear (quadratic) fashion. The above investigation was recently generalized by Weidman and Magyari [15] to continuous surfaces stretching with arbitrary polynomial velocities, which they designed as superstretching, and then the arbitrary stretching was also extended to micropolar fluids [16]. Cortell [17] proceeded to analyze momentum and heat transport on a continuous planar surface moving in a parallel stream. Mention is also made of the recent published paper by Palani and Kim [18] in which the effects of a first-order homogeneous chemical reaction in an unsteady free convective flow past a vertical plate were numerically investigated. Muhaimin and Kandasamy [19] reported interesting results for mass transfer in free, forced, and mixed convection flows along a porous wedge with magnetic field, chemical reaction, and suction/blowing influences. Further, Abdel-Rahman [20] studied thermal diffusion and MHD for Soret and Dufour’s effects on Hiemenz flow and mass transfer with variable viscosity.

Exponentially stretching sheet was studied by Elbashbeshy [21]. Further, the effects of suction/blowing and viscous dissipation on the steady nonlinear viscous flow and heat transfer over a horizontal shrinking permeable sheet were recently studied by Rohni et al. [22].

The present author has also studied [23, 24] flow and heat transfer onto a nonlinearly stretching sheet for three different types of thermal boundary conditions on the sheet: constant surface temperature (CST case), prescribed surface temperature (PST case), and prescribed heat flux (PHF case). These investigations were based on the assumption that the effects of thermal radiation and viscous dissipation were enclosed in the energy equation. Flow and heat transfer characteristics in a viscous fluid over a nonlinearly stretching sheet are recently analyzed analytically by Akyildiz and Siginer [25].

Another study which belongs to these classes of problems was recently done by Bataller [26]. He stated a true similarity problem and presented a numerical analysis in connection with the laminar boundary layer flow induced in a quiescent fluid by an impermeable stretching flat surface with velocity . Realize that, in general, , being an arbitrary positive constant (not necessarily an integer), as was very recently pointed out by Prasad and Vajravelu [27]. Among others, the quadratic stretching sheet case was devised by Khan and Sanjayanand [28]. The MHD flow of a viscous fluid over a sheet stretching with nonlinear (quadratic) velocities without suction/blowing effects at the surface was recently considered by Raptis and Perdikis [29] and Kelson [30]. These types of fluid flow problems have also been extended by Ahmad and Asghar [31] to a sheet stretching with arbitrary velocities in the case of fluid flow of a second grade fluid.

Hence, the purpose of the present paper is to extend the most recent studies [14, 29, 30] to the case of permeable wall with magnetic field influences together with mass transfer. As in Kumaran and Ramanaiah [14] was made, we also adopt both quadratic polynomial velocity and linear mass flux at the surface. However, in the previous study [14], magnetic field influences were not included; hence, the present paper is devoted to analyze MHD viscous incompressible flow and mass transfer with quadratic polynomial velocity at the sheet.

The standard approach here is to find similarity solution for the complete system of governing equations. As a result, in some few cases, the system of partial differential equations permits a similarity transformation to convert it into a system of ordinary differential equations, and then similarity solutions can be found. Although velocity and concentration fields are interdependent (through the relevant parameter ), we in the following consider the steady MHD viscous flow and mass transfer due to nonlinearly (quadratic) stretching sheet with a fictitious destructive chemical reaction effects when the full set of governing equations (i.e., continuity, momentum, and mass transfers) permits the similarity transformation for various functions of exponential type. On this assumption, a complete set of governing equations can be satisfied, and hence an exact solution for the stated problem can be found. As we will see, different from early works regarding the outlined problem (see, e.g., [29, 30]), a remarkable prescribed wall concentration is assumed which corresponds to the boundary-layer concentration solution at , together with an appropriate -dependent wall transpiration . Although these boundary conditions may be difficult to realize in practice, we are carrying out these analyses with a principal emphasis on the similarity solutions for the boundary layer MHD flow and mass transfer, which constitute the basis for real engineering problems. The present author [32] has treated the case of MHD flow and mass transfer for second grade fluid over a linearly stretching sheet with chemically reactive species when the concentration of the reactant is maintained at a prescribed constant wall value , different from the present boundary assumption in which .

The problem under consideration is governed by three parameters, that is, the magnetic parameter , the suction/blowing parameter , and the Schmidt number Sc. The reaction rate parameter ( represents a destructive chemical reaction; indicates generative chemical reaction; and is the case when a nonreactive species problem takes place) should be obtained as . It is hoped that the obtained results will not only provide useful information for applications but also serve as a completion to the related previous studies. On the other hand, in order to assess the accuracy of the presented results, several comparisons are made between the present results and previously published results, for special cases.

2. The Flow

Let us suppose steady, constant properties, laminar and two-dimensional flow of an incompressible, electrically conducting, and Boussinesq viscous fluid subject to a transverse uniform magnetic field which is applied in the positive -direction past a flat, horizontal, and porous sheet coinciding with the plane , the flow being confined to . The motion of the fluid is generated due to nonlinearly (quadratic) stretching the sheet with the application of two equal and opposite forces, which are applied along the -axis so that the wall is stretched keeping the origin fixed. The velocity of the fluid far away from the plate is equal to zero, whereas fluid suction/injection is imposed at the plate surface. The magnetic Reynolds number is considered to be small so that the induced magnetic field is negligible in comparison to the applied magnetic field. The governing equations are where and are the coordinates along and normal to the sheet, respectively, and are the velocity components in and directions, is the kinematic viscosity, is the fluid density, is the strength of a uniform magnetic field, which is applied along the -axis, and is the electric conductivity. The boundary conditions governing the flow subproblem are Here , are positive nonzero constants, and , are the stretching and transpiration sheet velocities, respectively.

The similarity variable is defined as and we also introduce the dimensionless variables and related to velocity components as Taking into account (2.4) and (2.5), one further easily can write and also Then, on substitution into (2.2), one can arrive at and we directly obtain from (2.8) the three following ordinary differential equations: where a prime denotes differentiation with respect to the independent similarity variable , and is the magnetic parameter. It is necessary to note that Kelson [30] revisited [29] in order to include (2.11). On the other hand, realize that the mass conservation equation (2.1) is automatically satisfied (see (2.6)). Equations (2.9)–(2.11) are to be solved subject to the transformed boundary conditions: It is clear from (2.11) and (2.12) that

An analytical solution of (2.9), (2.10), (2.12), and (2.13) for the cases , = 0, and was recently derived by Kelson [30]; however, the more relevant cases of and have not been treated in [30]. It is clearly necessary to note that Kelson’s [30] solution for (i.e., (9) in [30]) does not satisfy (2.11) at , and then the stated problem would still remain unsolved, and hence all the above justifies our appendix. The novelty here, compared with other works cited above, lies in the attainment of exact solutions describing MHD flow and mass transfer generated by a flat sheet undergoing nonlinear (quadratic) stretching, which we outline in the following.

Fortunately, (2.9) with the boundary conditions , , and has an exact solution in the form (see also Cortell [9]) where Thus, from (2.15) and (2.16), we get a simple exact analytical solution of (2.9) with the above associate boundary conditions and we will utilize this solution for the function .

On the other hand, (2.11) with the boundary conditions , and has also an exact solution in the form where will be given by (2.16).

Furthermore, on substitution into (2.10) of the functions and , which are given by (2.15) and (2.17), respectively, one immediately concludes that (2.10) becomes which holds when is given, precisely, by (2.16). From the above, we can see that (2.9)–(2.11), which govern our stated problem, are satisfied by (2.15), (2.16), and (2.17).

Then, a complete analytical solution for the permeable case when the sheet is stretched in a nonlinear (quadratic) fashion is found as follow: where is obtained from (2.16).

The associate boundary conditions read for , for . It is clearly necessary to note that the latter can also be obtained from (2.12)–(2.14) setting .

Based on the second equation of (2.19), the transverse velocity component of the edge of the boundary layer (i.e., the entrainment velocity) is obtained as This velocity component depends on , and flow parameters, and it is related to the amount of fluid which is entrained by the movement of the surface (dragged by the sheet).

In addition, one immediately obtains from the second equation of (2.5) and (2.19)

One also sees from (2.23) that the suction/injection imposed at the plate surface (i.e., the transpiration velocity) turns out in a -dependent blowing when is either zero or negative.

Figures 13 extract graphical information about the effects of and onto velocity profiles when  s−1 m−1s−1, and  m. Figure 1 shows the effect of the magnetic parameter onto velocity profiles when . It can be seen that the velocity decreases when the magnetic parameter increases, and then the boundary layer thickness shrinks. In Figure 2, we have plotted the velocity profiles showing the effect of the parameter at nonzero magnetic field, say . One can see from this graph that as the parameter decreases from positive to negative, the magnitude of the velocity largely increases in the boundary flow region.

On the other hand, one can analyse, from Figure 3, the influence of the parameter and sign of on velocity profiles graphically. The velocity profiles boost as the parameter decreases from positive to negative, and also the magnetic field effects acquire a major significance when is negative (i.e., variable blowing). It is also found that for strong suction () the velocity decreases in the boundary flow region. An explanation for this behavior is that the fluid is brought closer to the sheet in the case of suction, and the opposite behavior is observed, however, in the case of blowing ().

It should be noted that (2.8) reduces to the simpler (2.9) for the case when = 0 (i.e., linear stretching). For this special case and in order to asses the accuracy of the present method, we compare values of the skin friction coefficient between the present results and those of Andersson et al. [33] and Vajravelu et al. [34] for a Newtonian fluid in Table 1.

3. The Mass Transfer

The concentration field is governed by the boundary layer mass diffusion equation (see [35, 36]) which in the boundary layer approximation reduces to where is the concentration of the species of the fluid, is the diffusion coefficient of the diffusing species in the fluid, and denotes the reaction rate constant of the first-order homogeneous and irreversible reaction. The diffusing species can be either destroyed or generated in the homogeneous reaction. The concentration of the reactant will be maintained at a prescribed wall distribution at the surface and is assumed to vanish far away from the surface. As we will outline in the following, the effects onto concentration profiles of the magnetic parameter are identical to those in (3.1) due to a first-order destructive chemical reaction, which has been accounted by enclosing the last term in the right-hand side of (3.1). Thus, the appropriate boundary conditions for the concentration (3.1) are given by

From (3.1), it is seen that a peculiar analytical solution will be given in the form with a wall concentration distribution where is a constant and is the Schmidt number.

By inserting (2.19) and (3.3) into (3.1), we get where is the reaction rate parameter. Besides, replacing by into (3.5), we then find that (2.18) and (3.5) are the same. In other words, the effects of a magnetic field onto concentration distributions can now be analyzed assuming , and then there are no magnetic field effects () when a nonreactive species problem () takes place.

The nondimensional concentration will be Summarily, the solutions of the corresponding momentum transfer () and mass transfer subproblems have been given in the forms where would satisfy the following two equations:

It is consequently found that the velocity and concentration fields are interdependent so that from (3.8) we get and therefore the parameters , , and are related.

Also, from (3.9), one deduces that the concentration solutions exist only in the case of a fictitious destructive chemical reaction (i.e., reactant consuming , chemical reactions).

On the other hand, the wall transpiration and the prescribed wall concentration will be given in the form

Furthermore, we note that from the first equation of (3.7) the dimensionless velocity can be expressed as and from the third equation of (3.7) we have . It is then obvious that the velocity profile is identical to the concentration profile for Schmidt number unity (i.e., ). On the other hand, will result that the skin friction coefficient appears to be always negative and that means the stretched surface exerts a drag force on the fluid. Physically, this is reasonable since the development of the momentum boundary layer is generated solely by the stretching surface. We also note that the Schmidt number Sc gives no influences on the velocity boundary layer, which is clear from (3.7).

Bearing in mind (2.16), (2.19), and (3.9), we get

In order to more fully characterize the behavior of the concentration profiles for a viscous fluid over a nonlinear (quadratic) permeable stretching sheet, with respect to the magnetic parameter , we have depicted in Figure 4 the variation of for three values of at Sc = 0.24 (representing diffusing chemical species of H2 in air) and (impermeable surface). It can be seen that the concentration decreases when the magnetic parameter increases. Further, Figure 5 graphically confirms the later for all .

Furthermore, as expected, at nonzero magnetic field, say , Figure 6 shows that an increase in Schmidt number Sc yields a decrease in concentration profiles, whereas for fixed Sc, as the parameter decreases from positive to negative, the magnitude of the concentration increases in the boundary flow region. All of these tendencies are consistent with the results of Vajravelu et al. [34].

Finally, we also graphically present the effects of an increase in the magnetic field intensity (through parameter) onto concentration profiles for both suction and injection situations and for three selected Schmidt number Sc values, say 0.62, 0.78, and 2.62, which correspond to diffusing chemical species in air of practical interest like H2O, NH3, and propyl benzene, respectively [37]. It is found that, irrespective of Sc, an increase in magnetic field intensity causes a decrease in the concentration boundary layer thickness. Moreover, as the Schmidt number Sc increases due to a diminution in the chemical species molecular diffusivity, the thickness of the concentration boundary layer shrinks (see Figures 7, 8 and 9).

In addition, to show the excellent agreement between the present results and those earlier ones from the open literature, use is made of the relevant work by Vajravelu et al. [34], who studied steady MHD non-Newtonian power-law fluid flow over an impermeable linearly stretching sheet with species concentration and chemical reaction. For comparison, one should take the power-law index , when the fluid is simply a Newtonian fluid. In [34], the continuous stretching sheet is assumed to have species diffusion where is the distance from the slit, is a characteristic length, is a constant, is the species concentration parameter, is the concentration at the plate, and indicates the same but far away from the plate surface. Use is made in the present work of , cf. (3.4). Hence, on comparison, it should be (in [34]) and (in the present paper). Further, (3.11) with will provide the present results.

From Table 2, one can see the excellent agreement between Vajravelu et al.’s [34] results, which were obtained numerically by a finite difference scheme known as the Keller-Box method and the results of the analytical solution.

4. Conclusion

The exact solution corresponding to the motion and mass transfer of a viscous fluid over a permeable nonlinearly (quadratic) stretching sheet is established. The effects of an applied transverse magnetic field are also considered. A suitable -dependent distribution of the transpiration velocity should be applied together with a variable prescribed wall concentration . Both velocity and concentration profiles are interdependent through the relevant parameter . The current analytical study can be considered to be more general than the existing studies due to the fact that magnetic field effects can now be analyzed. The velocity and concentration profiles strongly depend upon the parameters , Sc, and and tend to zero at infinity asymptotically. The main results of the paper can be summarized as follows.(1) As the parameter decreases from positive to negative, the magnitude of the velocity largely increases in the boundary flow region. (2) The velocity and concentration decrease as increases.(3) As the parameter decreases from positive to negative, the magnitude of the concentration increases in the boundary flow region. (4) An augment in the Schmidt number Sc shrinks the concentration profiles.(5) The magnetic field effects on velocity and concentration profiles acquire a major significance when is negative (blowing).(6) For Sc = 1, the dimensionless velocity and concentration profiles are identical.(7) For the case of , the dimensionless concentration profile is smaller in comparison to the dimensionless velocity profile , and this trend is reversed for the case of .(8) The Schmidt number Sc exerts no effects on the momentum boundary layer.


Our aim in this appendix is to interpret the obtained results of Kelson [30] and to present a formulation to the stated problem (i.e., similarity solutions for MHD flow due to a permeable and nonlinearly (quadratic) stretching sheet). In [30], an exact solution of (2.9), (2.10), (2.12), and (2.13) for the case , and is presented as and one immediately sees that (A.2) does not satisfy (2.11); one can therefore point out that the stated problem would remain unsolved. Different from (A.2), our purpose for the function g is given by (2.17) with , and then . Further, we have from (2.15) which coincides with (8) of [30]. One can then see how (A.3) and (A.4) constitute a solution of (2.9)–(2.13).That is, by replacing (A.2) in [30] with (A.3), one can study the stated problem for the cases , , and . On the other hand, (A.2) and (A.3) are quite different as can be viewed graphically (see Figure 10).

Moreover, in the hypothesis that (A.3) remains unknown for us, one can expand the function in a Taylor series about , as it is outlined in the following. It is clear from (2.11) evaluated at that and incorporating the boundary condition , we get Differentiating (2.11) with respect to results in This yields Differentiating (A.7) with respect to will result in and hence (A.9) evaluated at yields Applying the Taylor series, we get Incorporating the boundary condition , all the above yield A comparison of approximate solution for sufficiently small values of , (A.12) with the exact solution (A.3), is given in Table 3, showing excellent agreement. Obviously, if the solution procedure continues, the accuracy will improve.

Bearing in mind (3.7), the velocity components , can be transcribed in the forms and also from (2.22) and (2.23), we get