A New Numerical Approach of MHD Flow with Heat and Mass Transfer for the UCM Fluid over a Stretching Surface in the Presence of Thermal Radiation

This paper numerically investigates the magnetohydrodynamic boundary layer flow with heat and mass transfer of an incompressible upper-convected Maxwell fluid over a stretching sheet in the presence of viscous dissipation and thermal radiation as well as chemical reaction. The governing partial differential equations are transformed into a system of ordinary differential equations by using suitable similarity transformations. The resultant highly nonlinear ordinary differential equations are then solved using spectral relaxationmethod.The results are obtained for velocity, temperature, concentration, skin friction, andNusselt number. The effects of various material parameters on the flow with heat and mass transfer and the dimensionless variables are illustrated graphically and briefly discussed.


Introduction
In the past few decades, the studies of boundary layer flows of Newtonian and non-Newtonian fluids of stretching surfaces have received great attention by virtue of their numerous applications in the fields of metallurgy, chemical engineering, and biological systems.These applications include geothermal reservoirs, wire and fiber coating, food stuff processing, reactor fluidization, transpiration cooling, enhanced oil recovery, packed bed catalytic reactors, and cooling of nuclear reactors.The prime aim in extrusion is to keep the surface quality of the extricate.Coating processes demand a smooth glossy surface to meet the requirements for the best appearance and optimum properties.Sakiadis [1,2] did pioneering work on boundary layer flow on a continuously moving surface.After that many investigators discussed various aspects of the stretching flow problem (see, e.g., Chiam [3], Crane [4], Liao and Pop [5], Khan and Sanjayanand [6], Abel and Mahesha [7], and Fang et al. [8], among others).
A number of industrial fluids such as molten plastics, artificial fibers, blood, polymetric liquids, and food stuff exhibit non-Newtonian fluid behaviour.In many industrial processes, cooling continuous strips or filaments is done by drawing them through a quiescent fluid.It must be noted that during these processes, the strips are sometimes stretched or shrunk.Therefore the properties of the final product depend to a great extent on the rate of cooling.The rate of cooling can be controlled and the desired characteristics of the final product can be obtained by drawing the strips in electrically conducting fluids subjected to uniform magnetic fields.
In recent years, MHD flows of viscoelastic fluids above stretching sheets have also been studied by various researchers (Liu [9], Cortell [10], among others).This is the simplest subclass of viscoelastic fluid known as the second grade fluid.However, a non-Newtonian second grade fluid does not give meaningful industrial results for highly elastic fluids such as polymer melts, which occur at high Deborah number (Hayat et al. [11]).For theoretical results to become of any industrial use, more realistic viscoelastic fluid models such as upper-convected Maxwell model should be used in the analysis.
Alizadeh-Pahlavan et al. [12] investigated using a twoauxiliary-parameter homotopy analysis method for the problem of laminar MHD flow of an upper-convected Maxwell fluid above a porous isothermal stretching sheet.Aliakbar et al. [13] analyzed the influence of Maxwell fluids above stretching sheets.Abel et al. [14] performed an analysis to investigate the influence of MHD and thermal radiation on the two-dimensional steady flow of an incompressible, upper-convected Maxwell fluid.Motsa et al. [15] investigated the MHD boundary layer flow of an incompressible upperconvected Maxwell fluid over a porous stretching surface.Hayat et al. [16] investigated MHD flow of a more realistic viscoelastic fluid model above a porous stretching sheet.Sadeghy et al. [17] also studied MHD flows of upperconvected Maxwell model fluids above porous stretching sheet using homotopy analysis method.Raftari and Yildrim [18] approximated an analytical solution of the magnetohydrodynamic boundary layer flow of an upper-convected Maxwell fluid over a porous stretching sheet.More recently, Alinejad and Samarbakhsh [19] numerically investigated the flow and heat transfer characteristics of the incompressible viscous flow over a nonlinear stretching sheet with viscous dissipation.
Hayat et al. [20] proposed the modified decomposition method and the Padé approximants for the MHD flow over a nonlinear sheet.In recent years, scientists and engineers have endeavored to develop more accurate and fast converging numerical and/or analytical techniques.Methods such as homotopy and their hybrid techniques have been extensively applied to solve nonlinear equations (Liao [21], He [22], Domairry and Nadim [23], Shateyi and Motsa [24], Sheikholeslami et al. [25], and Ganji [26], among others).Recently, successive linearization method has been reported and successfully utilized in solving boundary value problems (see, e.g., Motsa and Shateyi [27], Motsa et al. [28], and Shateyi and Motsa [29]).
The aim of this study is to investigate the effects of thermal radiation and viscous dissipation on steady MHD flow with heat and mass transfer of an upper-convected Maxwell fluid past a stretching sheet in the presence of a chemical reaction.Secondly, we aim to use a recently developed iterative method known as spectral relaxation method (SRM), and details of this method can be found in [30][31][32].

Mathematical Formulation
We consider the steady and incompressible MHD boundary layer flow with heat and mass transfer of an electrically conducting fluid obeying UCM model over a stretching sheet in the presence of thermal radiation.The flow is generating by the stretching of the sheet by applying two equal and opposite forces along the -axis, keeping the origin fixed and considering the flow to be confined to the region  > 0. We assume that the continuous stretching sheet has a linear velocity,  = , with  as the stretching rate and  being the distance from the slit.We impose a uniform magnetic field of strength  0 along the -axis, and the induced magnetic field is negligible.This assumption is valid on a laboratory scale under the assumption of small magnetic Reynolds number, and the external electric field is zero.We also assume that the boundary layer approximations are applicable to all momentum, energy, and mass equations.Although this theory is incomplete for viscoelastic fluids, but it is more plausible for Maxwell fluids as compared to other viscoelastic fluid models (Renardy [33]).Following Sadeghy et al. [17] among others, in a two-dimensional flow, the equation of continuity, the equation of motion, and the diffusion equations can be written as where  and V are velocity components in the and -directions, respectively,  is the fluid density,  is the electrical conductivity,  0 is the uniform magnetic field,  is the temperature,   is the specific heat at constant pressure,  is the thermal conductivity,  is the concentration of the species diffusion,  is the diffusion coefficient of the diffusion species, ] is the kinematic viscosity,  is the relaxation time, and   denotes the reaction rate constant of the th-order homogeneous and irreversible reaction.The appropriate boundary conditions are By using the Rosseland diffusion approximation (Hossain et al. 1999, among other researchers), the radiative heat flux,   , is given by where  * and   are the Stefan-Boltzmann constant and the Rosseland mean absorption coefficient, respectively.We assume that the temperature differences within the flow are sufficiently small such that Using ( 6) and (7) in the last term of (3) we obtain

Similarity Transformation.
To make the problem amenable we introduce the following nondimensional quantities: where () is the dimensionless stream function and  is the similarity variable,  is the dimensionless temperature, and  is the mass concentration.The continuity equation is automatically satisfied through the variables.Then introducing the relations ( 9) into ( 2)-(3), we obtain the following nonlinear system of ordinary differential equations: Here  2 =  2 0 / and  =  are magnetic and elastic parameters, respectively, Pr =   / is the Prandtl number, Ec =  2 /    is the Eckert number,  = 4 *  3 ∞ /  is the thermal radiation parameter,  =    −1 0  −2   / 2 is the chemical reaction parameter, and Sc = /] is the Schmidt number.The boundary conditions are

Method of Solution
The Successive Relaxation Method (SRM) begins by letting so that   =   and   =   .Consequently, (10) through (12) become Proceeding in a manner similar to the Gauss-Seidel method, ( 14) and ( 15) through ( 16) are replaced by the following recursive formulae: The Chebyshev spectral collocation method will be used to solve (18) through (21).First, we replace the semi-infinite interval [0, ∞) with the closed interval [0, ], where  is sufficiently large.It is convenient to use the change of variable to map the interval [0, ] on the -axis onto the interval [−1, 1] on the -axis.On [−1, 1] we form a computational grid where are the Chebyshev collocation points.The derivative ℎ  (  ) at each collocation point is evaluated using formula where  is the Chebyshev differentiation matrix.Successive application of Chebyshev differentiation reveals the more general formula Chebyshev differentiation transforms ( 18) through (21) to discrete form: where Initial approximations needed to drive this iterative scheme must satisfy boundary conditions (13).Suitable choices are

Results and Discussion
In this section we give the SRM results for the main parameters that have significant effects on the fluid flow velocity and temperature.We remark that all the SRM results presented in this work were obtained using  = 50 collocation points, and also convergence was achieved after as few as five iterations.Also the infinity value ( ∞ ) was taken to be 50.It is also important to note that the magnetic field is taken quite strong by assigning large values of  to ensure the occurrence of steady flow near the sheet.Unless otherwise stated, the default values for the parameters are taken as  = 1,  = 0.1,  = 1, Pr = 0.71, Ec = 0.1, Sc = 0.2,  = 0.2, and  = 2.In order to validate the numerical method, it was compared with the MATLAB routine V4 which is an adaptive Lobatto quadrature iterative scheme.Table 1 presents a comparison between SRM approximate results and the V4 results for selected default values of the magnetic parameter .It can be seen from this table that there is an excellent agreement between the results from the two methods.Analyzing Table 1 shows that an increase in the magnetic field strength leads to an increase in the skin-friction but a decrease to the Nusselt number.This is physically expected as application of a transverse magnetic field produces a drag force which then reduces the flow velocity but generates heat within the fluid.Table 2 gives a comparison of the SRM results to those obtained by the V4 for different values of the elasticity parameter .We again observe that the results from the two methods agree very well giving confidence to the current proposed method.Increasing the fluid elasticity parameter leads to the increase in the skin friction coefficient but a decrease in the heat transfer coefficient.Figures 1 and 2 show the influence of the elasticity parameter  on the -velocity and V-velocity profiles, respectively.From both these figures, we observe that an increase in the elasticity parameter results in the decrease in both velocity components and a decrease in the thickness of the momentum boundary layer.
Figures 3 and 4 depict the effects of the magnetic field parameter  on the -velocity and V-velocity profiles, respectively.From these figures we clearly observe that increasing magnetism significantly reduces the thickness of the boundary layer, thereby reducing the velocity components.Physically, the application of the transverse magnetic field presents a damping effect on the flow velocity by producing a drag force that opposes the fluid motion.
Figure 5 depicts the effect of increasing the elasticity parameter  on the temperature distribution.A decrease in the streamwise velocity component, , can result in a decrease in the amount of heat transferred on the surface sheet.Similarly, a decrease in the transverse velocity component, V, means that the amount of fresh fluid which is extended from the lower-temperature region outside the boundary layer and directed towards the sheet is reduced, thereby reducing the rate of heat transfer.These two effects on the velocity components in the same direction reinforce each other.Thus, an increase in the elastic number increases the temperature distribution in the fluid as depicted in Figure 5.
Figure 6 represents the dimensionless temperature for different values of the magnetic field parameter .From this figure we clearly see that the temperature profiles increase with the increase of the magnetic field parameter.Thus the applied magnetic field tends to heat the fluid and thus reduces the heat transfer from the wall which in turn enlarges the thermal boundary layer thickness.
In Figure 7 we display the effect of the Eckert number on the temperature profiles.An increase in the values of the Eckert number is seen to increase the temperature of the fluid at any point above the sheet.Increasing the Eckert number allows energy to be stored in the fluid region as a consequence of dissipation due to viscosity and elastic deformation.
The influence of thermal radiation on the temperature profiles is shown in Figure 8.It can be seen that the thermal boundary layer thickness increases as  increases.This induces the decrease in the absolute value of the temperature gradient at the surface.Thus, the heat transfer rate at the surface decreases with increasing , thereby causing the temperature profiles to increase.
Figure 9 depicts the effect of increasing the Prandtl number on the fluid temperature distribution.An increase in  the Prandtl number as expected is seen to reduce the fluid temperature above the sheet.This is because as the Prandtl number increases, the thermal boundary layer becomes thinner.Thus the rate of thermal diffusion drops, resulting in the fluid temperature dropping as well.
Figure 10 shows how the elastic parameter  affects the concentration profiles.As the elasticity parameter reduces the flow velocity, it means that less fluid is taken away at any given point resulting in the concentration profiles increasing.The same phenomenon happens when the values of the magnetic parameter increase as can be clearly seen in Figure 11.The reduction of flow velocity as the result of increasing the strength of the magnetic field causes the fluid concentration to increase as less fluid is taken downstream at any given point.
Figure 12 presents the profiles of the concentration for selected default values of the chemical reaction parameter  ≥ 0. It is observed from this figure that an increase in the values of chemical reaction parameter leads to a decrease in the concentration profiles.The concentration boundary layer becomes thin as the reaction parameter increases.
Figure 13 shows graphically the effect of increasing the reaction-order parameter .The effect of  is seen as to increase the fluid concentration.

Conclusion
The present work analyzed the MHD flow with heat and mass transfer within a boundary layer of an upper-convected Maxwell fluid above a stretching sheet in the presence of viscous dissipation, thermal radiation, and chemical reaction.Numerical results are presented in tabular/graphical form to elucidate the details of flow with heat and mass transfer characteristics and their dependence on the various physical parameters.The accuracy of the SRM is validated against the MATLAB in-built V4 routine for solving boundary value problems.
(1) We observe that the flow velocity is decreased when the magnetic parameter increases.Also an increase in the elasticity parameter results in velocity decrements.However, both the temperature and concentration profiles are enhanced by increasing the values of the magnetic parameter as well as the elasticity parameter.
(2) The dimensionless temperature  increases with increase in the thermal radiation but decreases with increasing Prandtl number.
(3) The effect of the chemical reaction  is to decrease the fluid concentration while the concentration profiles are increased as the order of reaction  is increased.
(4) It is shown that the skin-friction increases with an increase of the magnetic parameter and elasticity parameter, but the Nusselt number and Sherwood number (not shown for brief) are found to be decreased as these parameters are increased.
(5) Finally the Nusselt number increases with increasing values of the Prandtl number but decreases as thermal radiation increases.

Figure 1 :Figure 2 :Figure 3 :Figure 4 :Figure 5 :
Figure 1: Graph of the SRM solutions of the -velocity for different values of .

Figure 6 :Figure 7 : 1 Figure 8 :
Figure 6: Variation of the temperature distribution for different values of .

2 Figure 9 :Figure 10 :
Figure 9: The effect of the Prandtl number on the temperature distribution.

Figure 11 :Figure 12 :
Figure 11: Graph of how the magnetic parameter affects the concentration distribution.

Figure 13 :
Figure 13: Graph of the SRM solutions of the concentration for different orders of the homogeneous reaction.

Table 1 :
Comparison of the SRM results of −  (0), −  (0) with those obtained by V4 for different values of the magnetic parameter.

Table 2 :
Comparison of the SRM results of −  (0), −  (0) with those obtained by V4 for different values of the elastic parameter.