New Developments in Fluid Mechanics and Its Engineering Applications
View this Special IssueResearch Article  Open Access
S. Shateyi, G. T. Marewo, "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", Mathematical Problems in Engineering, vol. 2013, Article ID 670205, 8 pages, 2013. https://doi.org/10.1155/2013/670205
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
Abstract
This paper numerically investigates the magnetohydrodynamic boundary layer flow with heat and mass transfer of an incompressible upperconvected 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 relaxation method. The results are obtained for velocity, temperature, concentration, skin friction, and Nusselt 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.
1. Introduction
In the past few decades, the studies of boundary layer flows of Newtonian and nonNewtonian 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 nonNewtonian 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 nonNewtonian 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 upperconvected Maxwell model should be used in the analysis.
AlizadehPahlavan et al. [12] investigated using a twoauxiliaryparameter homotopy analysis method for the problem of laminar MHD flow of an upperconvected 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 twodimensional steady flow of an incompressible, upperconvected 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 upperconvected 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 upperconvected 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–32].
2. 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 . 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 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 twodimensional flow, the equation of continuity, the equation of motion, and the diffusion equations can be written as where and are velocity components in the  and directions, respectively, is the fluid density, is the electrical conductivity, 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 thorder 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 StefanBoltzmann 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
2.1. 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 and are magnetic and elastic parameters, respectively, is the Prandtl number, is the Eckert number, is the thermal radiation parameter, is the chemical reaction parameter, and is the Schmidt number. The boundary conditions are
3. 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 GaussSeidel 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 semiinfinite interval with the closed interval , where is sufficiently large. It is convenient to use the change of variable to map the interval on the axis onto the interval on the axis. On 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
4. 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 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 , , , , , , , and . In order to validate the numerical method, it was compared with the MATLAB routine which is an adaptive Lobatto quadrature iterative scheme. Table 1 presents a comparison between SRM approximate results and the 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 skinfriction 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 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 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 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, , means that the amount of fresh fluid which is extended from the lowertemperature 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 . 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 reactionorder parameter . The effect of is seen as to increase the fluid concentration.
5. Conclusion
The present work analyzed the MHD flow with heat and mass transfer within a boundary layer of an upperconvected 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 inbuilt 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 skinfriction 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.
Acknowledgments
The authors wish to acknowledge financial support from the University of Venda and NRF.
References
 B. C. Sakiadis, “Boundarylayer behaviour on continuous solid surfaces: I. Boundarylayer equations for twodimensional and axisymmetric flow,” AIChE Journal, vol. 7, no. 1, pp. 26–28, 1961. View at: Publisher Site  Google Scholar
 B. C. Sakiadis, “Boundarylayer behaviour on continuous solid surfaces: II. Boundarylayer equations for twodimensional and axisymmetric flow,” AIChE Journal, vol. 7, no. 2, pp. 221–225, 1961. View at: Publisher Site  Google Scholar
 T. C. Chiam, “Hydromagnetic flow over a surface stretching with a powerlaw velocity,” International Journal of Engineering Science, vol. 33, no. 3, pp. 429–435, 1995. View at: Google Scholar
 L. J. Crane, “Flow past a stretching plate,” Zeitschrift für Angewandte Mathematik und Physik, vol. 21, no. 4, pp. 645–647, 1970. View at: Publisher Site  Google Scholar
 S. Liao and I. Pop, “Explicit analytic solution for similarity boundary layer equations,” International Journal of Heat and Mass Transfer, vol. 47, no. 1, pp. 75–85, 2004. View at: Publisher Site  Google Scholar
 S. K. Khan and E. Sanjayanand, “Viscoelastic boundary layer flow and heat transfer over an exponential stretching sheet,” International Journal of Heat and Mass Transfer, vol. 48, no. 8, pp. 1534–1542, 2005. View at: Publisher Site  Google Scholar
 M. S. Abel and N. Mahesha, “Heat transfer in MHD viscoelastic fluid flow over a stretching sheet with variable thermal conductivity, nonuniform heat source and radiation,” Applied Mathematical Modelling, vol. 32, no. 10, pp. 1965–1983, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 T. G. Fang, C. F. F. Lee, and J. Zhang, “The boundary layers of an unsteady incompressible stagnationpoint flow with mass transfer,” International Journal of NonLinear Mechanics, vol. 46, no. 7, pp. 942–948, 2011. View at: Publisher Site  Google Scholar
 I. C. Liu, “Flow and heat transfer of an electrically conducting fluid of second grade in a porous medium over a stretching sheet subject to a transverse magnetic field,” International Journal of NonLinear Mechanics, vol. 40, no. 4, pp. 465–474, 2005. View at: Publisher Site  Google Scholar
 R. Cortell, “Flow and heat transfer of an electrically conducting fluid of second grade over a stretching sheet subject to suction and to a transverse magnetic field,” International Journal of Heat and Mass Transfer, vol. 49, no. 1112, pp. 1851–1856, 2006. View at: Publisher Site  Google Scholar
 T. Hayat, Z. Abbas, and M. Sajid, “Series solution for the upperconvected Maxwell fluid over a porous stretching plate,” Physics Letters A, vol. 358, no. 56, pp. 396–403, 2006. View at: Publisher Site  Google Scholar
 A. AlizadehPahlavan, V. Aliakbar, F. VakiliFarahani, and K. Sadeghy, “MHD flows of UCM fluids above porous stretching sheets using twoauxiliaryparameter homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 2, pp. 473–488, 2009. View at: Publisher Site  Google Scholar
 V. Aliakbar, A. AlizadehPahlavan, and K. Sadeghy, “The influence of thermal radiation on MHD flow of Maxwellian fluids above stretching sheets,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 3, pp. 779–794, 2009. View at: Publisher Site  Google Scholar
 M. S. Abel, J. V. Tawade, and J. N. Shinde, “The effects of MHD flow and heat transfer for the UCM fluid over a stretching surface in presence of thermal radiation,” Advances in Mathematical Physics, vol. 2012, Article ID 702681, 21 pages, 2012. View at: Publisher Site  Google Scholar
 S. S. Motsa, T. Hayat, and O. M. Aldossary, “MHD flow of upperconvected Maxwell fluid over porous stretching sheet using successive Taylor series linearization method,” Applied Mathematics and Mechanics, vol. 33, no. 8, pp. 975–990, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 T. Hayat, Z. Abbas, M. Sajid, and S. Asghar, “The influence of thermal radiation on MHD flow of a second grade fluid,” International Journal of Heat and Mass Transfer, vol. 50, no. 56, pp. 931–941, 2007. View at: Publisher Site  Google Scholar
 K. Sadeghy, H. Hajibeygi, and S. Taghavi, “Stagnationpoint flow of upperconvected Maxwell fluids,” International Journal of NonLinear Mechanics, vol. 41, no. 10, pp. 1242–1247, 2006. View at: Publisher Site  Google Scholar
 B. Raftari and A. Yildirim, “The application of homotopy perturbation method for MHD flows of UCM fluids above porous stretching sheets,” Computers and Mathematics with Applications, vol. 59, no. 10, pp. 3328–3337, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 J. Alinejad and S. Samarbakhsh, “Viscous flow over nonlinearly stretching sheet with effects of viscous dissipation,” Journal of Applied Mathematics, vol. 2012, Article ID 587834, 10 pages, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 T. Hayat, Q. Hussain, and T. Javed, “The modified decomposition method and Padé approximants for the MHD flow over a nonlinear stretching sheet,” Nonlinear Analysis: Real World Applications, vol. 10, no. 2, pp. 966–973, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 S. J. Liao, Beyond Perturbation: Introduction to Homotopy Analysis Method, CRC Series: Modern Mechanics and Mathematics, Chapman & Hall/CRC Press, Boca Raton, Fla, USA, 2003. View at: MathSciNet
 J. H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 34, pp. 257–262, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 G. Domairry and N. Nadim, “Assessment of homotopy analysis method and homotopy perturbation method in nonlinear heat transfer equation,” International Communications in Heat and Mass Transfer, vol. 35, no. 1, pp. 93–102, 2008. View at: Publisher Site  Google Scholar
 S. Shateyi and S. S. Motsa, “A new approach for the solution of threedimensional magnetohydrodynamic rotating flow over a shrinking sheet,” Mathematical Problems in Engineering, vol. 2010, Article ID 586340, 15 pages, 2010. View at: Publisher Site  Google Scholar
 M. Sheikholeslami, H. R. Ashorynejada, D. D. Ganji, and A. Yildirim, “Homotopy perturbation method for threedimensional problem of condensation film on inclined rotating disk,” Scientia Iranica, vol. 19, no. 3, pp. 437–442, 2012. View at: Publisher Site  Google Scholar
 D. D. Ganji, “A semianalytical technique for nonlinear settling particle equation of motion,” Journal of HydroEnvironment Research, vol. 6, no. 4, pp. 323–327, 2012. View at: Publisher Site  Google Scholar
 S. S. Motsa and S. Shateyi, “Successive linearisation solution of free convection nondarcy flow with heat and mass transfer,” in Advanced Topics in Mass Transfer, M. ElAmin, Ed., pp. 425–438, InTech, Rijeka, Croatia, 2011. View at: Google Scholar
 S. S. Motsa, Y. Khan, and S. Shateyi, “Application of piecewise successive linearization method for the solutions of the Chen chaotic system,” Journal of Applied Mathematics, vol. 2012, Article ID 258948, 12 pages, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 S. Shateyi and S. S. Motsa, “Variable viscosity on magnetohydrodynamic fluid flow and heat transfer over an unsteady stretching surface with Hall effect,” Boundary Value Problems, vol. 2010, Article ID 257568, 20 pages, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 S. S. Motsa, P. G. Dlamini, and M. Khumalo, “Solving hyperchaotic systems using the spectral relaxation method,” Abstract and Applied Analysis, vol. 2012, Article ID 203461, 18 pages, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 S. S. Motsa, P. Dlamini, and M. Khumalo, “A new multistage spectral relaxation method for solving chaotic initial value systems,” Nonlinear Dynamics, vol. 72, no. 12, pp. 265–283, 2013. View at: Publisher Site  Google Scholar
 S. S. Motsa and Z. G. Makukula, “On spectral relaxation method approach for steady von Kármán flow of a ReinerRivlin fluid with Joule heating, viscous dissipation and suction/injection,” Central European Journal of Physics, vol. 11, no. 3, pp. 363–374, 2013. View at: Publisher Site  Google Scholar
 M. Renardy, “High Weissenberg number boundary layers for the upper converted Maxwell fluid,” Journal of NonNewtonian Fluid Mechanics, vol. 68, no. 1, pp. 125–132, 1997. View at: Google Scholar
Copyright
Copyright © 2013 S. Shateyi and G. T. Marewo. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.