# 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

Authors: S. Shateyi; G. T. Marewo
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

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 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.


---

## Body

## 1. 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 two-auxiliary-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 upper-convected 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 upper-convected 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–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 thex-axis, keeping the origin fixed and considering the flow to be confined to the region y>0. We assume that the continuous stretching sheet has a linear velocity, u=bx, with b as the stretching rate and x being the distance from the slit. We impose a uniform magnetic field of strength B0 along the y-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
(1)∂u∂x+∂v∂y=0,(2)u∂u∂x+v∂u∂y+λ1[u2∂2u∂x2+v2∂2u∂y2+2uv∂2u∂x∂y]=ν∂2u∂y2-σB02ρu,(3)u∂T∂x+v∂T∂y=kρcp∂2T∂y2+μρcp(∂u∂y)2-1ρcp∂qr∂y,(4)u∂C∂x+v∂C∂y=D∂2C∂y2-Kn(C-C∞)n,
where u and v are velocity components in the x- and y-directions, respectively, ρ is the fluid density, σ is the electrical conductivity, B0 is the uniform magnetic field, T is the temperature, cp is the specific heat at constant pressure, k is the thermal conductivity, C is the concentration of the species diffusion, D is the diffusion coefficient of the diffusion species, ν is the kinematic viscosity, λ is the relaxation time, and Kn denotes the reaction rate constant of the nth-order homogeneous and irreversible reaction. The appropriate boundary conditions are
(5)u(x,0)=Uw(x)=bx,v(x,0)=0,T(x,0)=T∞+Ts(xl)2,C(x,0)=C∞+Cs(xl)2,u(x,∞)=0,T(x,∞)=T∞,C(x,∞)=C∞.
By using the Rosseland diffusion approximation (Hossain et al. 1999, among other researchers), the radiative heat flux, qr, is given by
(6)qr=-4σ*3Ks∂T4∂y,
where σ* and Ks 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
(7)T4≈4T∞3T-3T∞4.
Using (6) and (7) in the last term of (3) we obtain
(8)∂qr∂y=-16σ*T∞33Ks∂2T∂y2.


### 2.1. Similarity Transformation

To make the problem amenable we introduce the following nondimensional quantities:(9)η=bνy,ψ=bνxf(η),θ(η)=T-T∞Tw-T∞,ϕ(η)=C-C∞Cw-C∞,
where f(η) 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:
(10)f′′′+ff′′-f′2-M2f′+β(2ff′f′′-f2f′′′)=0,(11)(1+43R)θ′′+Pr(fθ′-2f′θ)+PrEcf′′2=0,(12)ϕ′′+Sc[fϕ′-2f′ϕ-γϕn]=0.
Here M2=σB02/ρb and β=λb are magnetic and elastic parameters, respectively, Pr=μcp/k is the Prandtl number, Ec=bl2/cpTs is the Eckert number, R=4σ*T∞3/kKs is the thermal radiation parameter, γ=KnC0n-1bn-2xn/l2 is the chemical reaction parameter, and Sc=D/ν is the Schmidt number. The boundary conditions are
(13)f(0)=0,f′(0)=1,θ(0)=1,ϕ(0)=1,f′⟶0,θ⟶0,ϕ⟶0,asη⟶∞.


## 2.1. Similarity Transformation

To make the problem amenable we introduce the following nondimensional quantities:(9)η=bνy,ψ=bνxf(η),θ(η)=T-T∞Tw-T∞,ϕ(η)=C-C∞Cw-C∞,
where f(η) 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:
(10)f′′′+ff′′-f′2-M2f′+β(2ff′f′′-f2f′′′)=0,(11)(1+43R)θ′′+Pr(fθ′-2f′θ)+PrEcf′′2=0,(12)ϕ′′+Sc[fϕ′-2f′ϕ-γϕn]=0.
Here M2=σB02/ρb and β=λb are magnetic and elastic parameters, respectively, Pr=μcp/k is the Prandtl number, Ec=bl2/cpTs is the Eckert number, R=4σ*T∞3/kKs is the thermal radiation parameter, γ=KnC0n-1bn-2xn/l2 is the chemical reaction parameter, and Sc=D/ν is the Schmidt number. The boundary conditions are
(13)f(0)=0,f′(0)=1,θ(0)=1,ϕ(0)=1,f′⟶0,θ⟶0,ϕ⟶0,asη⟶∞.


## 3. Method of Solution

The Successive Relaxation Method (SRM) begins by letting(14)p=f′
so that p′=f′′ and p′′=f′′′. Consequently, (10) through (12) become
(15)p′′+fp′-p2-M2p+β(2fpp′-f2p′′)=0,(16)(1+43R)θ′′+Pr(fθ′-2pθ)+PrEcp′2=0,(17)ϕ′′+Sc[fϕ′-2pϕ-γϕn]=0.
Proceeding in a manner similar to the Gauss-Seidel method, (14) and (15) through (16) are replaced by the following recursive formulae:
(18)fr+1′=pr,fr+1(0)=0,(19)pr+1′′+fr+1pr+1′-pr2-M2pr+1p′′r+1+β(2fr+1prpr′-fr+12pr′′)=0,pr+1(0)=1,pr+1(∞)=0,(20)1Pr(1+4R3)θr+1′′+fr+1θr+1′hh-2pr+1θr+1+Ecpr+1′2=0,θr+1(0)=1,θr+1(∞)=0,(21)ϕ′′+Sc[fϕr+1′-2pr+1ϕr+1-γϕrn]=0,ϕr+1(0)=1,ϕr+1(∞)=0.
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,L], where L is sufficiently large. It is convenient to use the change of variable
(22)η(ξ)=1+ξ2L
to map the interval [0,L] on the η-axis onto the interval [-1,1] on the ξ-axis. On [-1,1] we form a computational grid -1=ξN<ξN-1<⋯<ξ0=1, where
(23)ξj=cos(πpN),j=0,1,…,N,
are the Chebyshev collocation points. The derivative h′(ξj) at each collocation point is evaluated using formula
(24)h′(ξj)=∑k=0N‍Djkh(ξk)⟹h′=Dh,
where D is the Chebyshev differentiation matrix. Successive application of Chebyshev differentiation reveals the more general formula
(25)dqdξqh=Dqh.
Chebyshev differentiation transforms (18) through (21) to discrete form:
(26)A1fr+1=B1,fr+1(ξN)=0,A2pr+1=B2,pr+1(ξN)=1,pr+1(ξ0)=0,A3Θr+1=B3,θr+1(ξN)=1,θr+1(ξ0)=0,A4Φr+1=B4,ϕr+1(ξN)=1,ϕr+1(ξ0)=0,
where
(27)A1=D,B1=pr,A2=D2+diag[fr+1]D-M2I,B2=pr2-β(2fr+1prpr′-fr+12pr′′),A3=(1+43R)D2+Pr(diag[fr+1]D-2diag[pr+1]),B3=-PrEcpr+1′2,A4=D2+Sc(diag[fr+1]D-2diag[pr+1]),B4=ScγΦrn.
Initial approximations needed to drive this iterative scheme must satisfy boundary conditions (13). Suitable choices are
(28)f0(η)=1-e-η,p0(η)=e-η,θ0(η)=e-η,ϕ0(η)=e-η.


## 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 usingN=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 M to ensure the occurrence of steady flow near the sheet. Unless otherwise stated, the default values for the parameters are taken as M=1, β=0.1, R=1, Pr=0.71, Ec=0.1, Sc=0.2, γ=0.2, and n=2. In order to validate the numerical method, it was compared with the MATLAB routine bvp4c which is an adaptive Lobatto quadrature iterative scheme. Table 1 presents a comparison between SRM approximate results and the bvp4c results for selected default values of the magnetic parameter M. 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 1
Comparison of the SRM results of-f′′(0), -θ′(0) with those obtained by bvp4c for different values of the magnetic parameter.
-
f
′
′
(
0
)
-
θ
′
(
0
)
M
b
v
p
4
c
SRM
b
v
p
4
c
SRM
0.0
1.42811186
1.42811186
0.57932373
0.57932373
1.0
1.87763199
1.87763199
0.40978550
0.40978550
2.0
2.71897650
2.71897650
0.27809236
0.27809236Table2 gives a comparison of the SRM results to those obtained by the bvp4c 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 u-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.Table 2
Comparison of the SRM results of-f′′(0), -θ′(0) with those obtained by bvp4c for different values of the elastic parameter.
-
f
′
′
(
0
)
-
θ
′
(
0
)
β
b
v
p
4
c
SRM
b
v
p
4
c
SRM
0.0
1.39382112
1.39382112
3.22830067
3.22830067
0.2
1.46105462
1.46105462
3.20679682
3.20679682
0.4
1.52344386
1.52344386
3.18653107
3.18653107Figure 1
Graph of the SRM solutions of theu-velocity for different values of β.Figure 2
Effect of theβ on the v-velocity profile.Figures3 and 4 depict the effects of the magnetic field parameter M on the u-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 3
Influence of the magnetic parameterM on the u-velocity profile.Figure 4
The variation of thev-velocity component for different values of M.Figure5 depicts the effect of increasing the elasticity parameter β on the temperature distribution. A decrease in the streamwise velocity component, u, 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 5
Graph of the SRM solutions of the temperature distribution for different values ofβ.Figure6 represents the dimensionless temperature for different values of the magnetic field parameter M. 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.Figure 6
Variation of the temperature distribution for different values ofM.In Figure7 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.Figure 7
Graph of the SRM solutions of the temperature profile for different values of Ec.The influence of thermal radiation on the temperature profiles is shown in Figure8. It can be seen that the thermal boundary layer thickness increases as R 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 R, thereby causing the temperature profiles to increase.Figure 8
Influence of thermal radiation of the fluid temperature.Figure9 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 9
The effect of the Prandtl number on the temperature distribution.Figure10 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 10
Graph of the SRM solutions of the concentration profile for different values ofβ.Figure 11
Graph of how the magnetic parameter affects the concentration distribution.Figure12 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 12
Influence of the chemical reaction on the concentration distribution.Figure13 shows graphically the effect of increasing the reaction-order parameter n. The effect of n is seen as to increase the fluid concentration.Figure 13
Graph of the SRM solutions of the concentration for different orders of the homogeneous reaction.


## 5. 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-builtbvp4c 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 n 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.



---

*Source: 670205-2013-09-17.xml*