# MHD Slip Flow of Newtonian Fluid past a Stretching Sheet with Thermal Convective Boundary Condition, Radiation, and Chemical Reaction

Authors: Reda G. Abdel-Rahman
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

An analysis is carried out to study the problem of heat and mass transfer flow over a moving permeable flat stretching sheet in the presence of convective boundary condition, slip, radiation, heat generation/absorption, and first-order chemical reaction. The viscosity of fluid is assumed to vary linearly with temperature. Also the diffusivity is assumed to vary linearly with concentration. The governing partial differential equations have been reduced to the coupled nonlinear ordinary differential equations by using Lie group point of transformations. The system of transformed nonlinear ordinary differential equations is solved numerically using shooting techniques with fourth-order Runge-Kutta integration scheme. Comparison between the existing literature and the present study was carried out and found to be in excellent agreement. The effects of the various interesting parameters on the flow, heat, and mass transfer are analyzed and discussed through graphs in detail. The values of the local Nusselt number, the local skin friction, and the local Sherwood number for different physical parameters are also tabulated.


---

## Body

## 1. Introduction

Investigations of laminar boundary flow of an electrically conducting fluid over a moving continuous stretching surface are important in many manufacturing processes, such as materials manufactured by polymer extrusion, continuous stretching of plastic films, artificial fibers, hot rolling, wire drawing, glass fiber, metal extrusion and metal spinning, cooling of metallic sheets, or electronic chips. The first and foremost work regarding the boundary-layer behavior in moving surfaces in quiescent fluid was considered by Sakiadis [1]. Heat and mass transfers on stretched surface with suction and injection was introduced by Fox et al. [2]. P. S. Gupta and A. S. Gupta [3] studied the same problem for linearly stretching sheet. Heat transfer past a moving continuous plate with variable temperature was studied by Soundalgekar and Murty [4] and Grubka and Bobba [5]. Subsequently, many researchers [6–12] worked on the problem of moving or stretching plates under different situations.MHD viscous flow over a stretching sheet in the presence of slip velocity was studied by many authors, such as Martin and Boyd [13], Fang and Lee [14], Ariel et al. [15], Andersson [16], Wang [17, 18], Mukhopadhyay and Anderson [19], Fang [20], and Hayat et al. [21]. Fang et al. [22] studied MHD viscous flow over a permeable shrinking sheet. They observed that the velocity at the wall increased with slip parameter. Mahmoud and Waheed [23] included the effects of slip and heat generation/absorption on MHD mixed convection flow of a micropolar fluid over a heated stretching surface. Recently, Hamad et al. [24] investigated heat and mass transfer over a moving porous plate with hydrodynamic slip and thermal convective boundary conditions. They found that the wall slip velocity and wall shear stress are a decreasing function of the slip parameter.Thermal radiation effects on an electrically conducting fluid arise in many practical applications such as electrical power generation, solar power technology, nuclear reactors, and nuclear waste disposal (see Mahmoud [25] and Chamkha [26]). Radiation effects on boundary layer flow with and without applying a magnetic field have been investigated by many authors [27–33] under certain conditions.Heat and mass transfer problems with a chemical reaction have received a considerable amount of attention due to their importance in many chemical engineering processes, for example, food processing, manufacturing of ceramics, and polymer production. On the other hand, the study of heat generation or absorption in moving fluids is important in problems dealing with chemical reactions and those concerned with dissociating fluids. Possible heat generation effects may alter the temperature distribution and, consequently, the particle deposition rate in nuclear reactors, electronic chip, and semiconductor wafers. Many investigators [34–40] displayed the effects of heat generation and chemical reactions on heat transfer over a stretching sheet under various conditions. Recently, thermal convective surface boundary conditions were used by several authors such as Ishak [41], Makinde [42], Yao et al. [43], Hamad et al. [24], and Ferdows et al. [44].Lie group point of transformations can be used to generate similarity transformations. It reduces the number of independent variables of the partial differential equations and keeps the system and associated initial and boundary conditions invariant. Reviews for the mathematical theory and applications of Lie group point of transformations to differential equations may be found in the texts by Olver [45], Bluman and Kumei [46], Hill [47], and Ibragimov and Kovalev [48]. This technique has been applied by many researchers to solve different flow phenomena over different geometries (see [49–52]).Motivated by the previous works, the aim of the present work is to study the problem of heat and mass of Newtonian fluid with thermal convective boundary condition, slip velocity, radiation, heat generation/absorption, and first-order chemical reaction taking into account the viscosity dependent on the temperature and the thermal diffusivity dependent on the concentration. This problem has not been studied before despite its applications in many engineering processes, such as nuclear waste, polymer production, dispersion of chemical pollutants through water-saturated soil, geothermal energy extractions, plasma studies, liquid metal fluids, and power generation systems.


## 2. Mathematical Formulation of the Problem

Consider laminar boundary layer flow heat and mass transfer of viscous incompressible Newtonian fluid over moving stretching sheet as shown in Figure1. A magnetic field of nonuniform B(x-) is applied parallel to the y- axis. The bottom surface of the plate is heated by convection from a hot fluid of temperature T-f, and this generates a heat transfer coefficient hf(x-), and T-∞ is the ambient fluid temperature. The mass concentration of the species at the plate is also maintained uniform C-w which is also higher than the ambient fluid concentration C-∞. Figure 1 shows the coordinates and flow model. Under these assumptions and boundary layer approximations, the governing equations are
(1)∂u-∂x-+∂v-∂y-=0,(2)u-∂u-∂x-+v-∂u-∂y-=1ρ∂∂y-(μ(T-)∂u-∂y-)-σ0ρB02(x-)u-,(3)u-∂T-∂x-+v-∂T-∂y-=κρcp∂2T-∂y-2-1ρcp∂qr∂y-+Q0(x-)ρcp(T--T-∞),(4)u-∂C-∂x-+v-∂C-∂y-=∂∂y-(D(C-)∂C-∂y-)-K(x-)(C--C∞),
subject to the boundary conditions
(5)u-=Uw(x-)+λ0(x-)μ(T-)∂u-∂y-,v-=Vw(x-),κ∂T-∂y=-hf(x-)[T-f-T-(x-)],C-=C-w,aty-=0,u-⟶0,T-⟶T-∞,C-⟶C-∞,asy-⟶∞,
where u- and  v- are the velocities in x- (along the sheet) and y- (normal to the sheet) directions, respectively, Vw(x-) is the mass transfer velocity, μ is the coefficient of dynamic viscosity, ν=(μ/ρ) is the kinematic coefficient of viscosity, D(C-) is the variable diffusivity, σ0 is the electric conductivity, ρ is the density of the fluid, cp is the specific heat at constant pressure, κ is the thermal conductivity, λ0(x-) is the velocity slip factor having dimension (velocity)-1, Uw(x-) is the velocity of moving plate, Vw(x-) is the suction or injection velocity through the surface, T- is the temperature, C- is the species concentration, K(x-) is the rate of chemical reaction, and Q0(x-) is the heat generation/absorption rate. The subscripts denote wall conditions and free stream conditions, respectively.Figure 1
Flow configuration.Using Rosseland approximation for radiation [53], we can write qr=-(4σ1/3κ*)(∂T-4/∂y-), where σ1 is the Stefan-Boltzmann constant and κ* is the mean absorption coefficient. Assuming that the temperature difference is within the flow such that T-4 may be expanded in a Taylor series and expanding T-4 about T-∞ and neglecting higher orders, we get T-4≅4T-∞3T--3T-∞4. Therefore, (3) becomes
(6)u-∂T-∂x-+v-∂T-∂y-=κρcp∂2T-∂y-2+16σ1T-∞33κ*ρcp∂2T-∂y-2+Q0(x-)(T--T-∞).We assume that the temperature dependent viscosity is varying linearly and is given by (see [54])
(7)μ(T-)=μ∞[a0+b0(T-f--T-)],
where μ∞ is the coefficient of viscosity far away from the plate and b0 and a0 are constants depend on the nature of the fluid. We choose a0=1.Also here, we assume that the concentration diffusivity varies linearly as (see [55])
(8)D(C-)=Dm[1+d(C--C-∞)],
where Dm is the constant concentration diffusivity and d is a constant.We introduce the nondimensional parameters(9)x=x-L,y=y-LRe,u=u-U0,v=v-U0Re,Uw(x)=U-w(x-)U0,Vw(x)=V-w(x-)U0,θ=T--T-∞T-f-T-∞,ϕ=C--C-∞C-w-C-∞,
with L being characteristic length, Re=LU0/ν is the Reynolds number, and U0 is a reference velocity. Also, we introduce stream function ψ defined as
(10)u=∂ψ∂y,v=-∂ψ∂x.
The continuity equation (1) is satisfied identically and (2), (4), and (6) yield
(11)G1=∂ψ∂y∂2ψ∂x∂y-∂ψ∂x∂2ψ∂y2-[1+ζ(1-θ)]∂3ψ∂y3+ζ∂θ∂y∂2ψ∂y2+M*B02(x)∂ψ∂y=0(12)G2=∂ψ∂y∂θ∂x-∂ψ∂x∂θ∂y-1Pr(1+4R3)∂2θ∂y2-γ*Q0(x)θ=0,(13)G3=∂ψ∂y∂ϕ∂x-∂ψ∂x∂ϕ∂y-1Sc(1+nϕ)∂2ϕ∂y2-nSc(∂ϕ∂y)2+α*K(x)ϕ=0.
The boundary conditions (5) now become
(14)∂ψ∂y=Uw(x)+λ*λ0(x)[1+ζ(1-θ)]∂2ψ∂y2,∂ψ∂x=V*Vw(x),∂θ∂y=-h*hf(x)[1-θ],ϕ=1,aty=0,(15)∂ψ∂y⟶0,θ⟶0,ϕ⟶0,asy⟶∞,
where ζ=b0(T-f-T-∞), n=d(C-w-C-∞),  M*=σ0L/ρU0, γ*=L/ρcpU0,  λ*=μ∞Re/L,  V*=1/Re, h*=L/κU0Re, R=4σ1T-∞3/κκ*, Sc=ν/Dm, Pr=ρνCp/κ, and α*=L/U0. The quantities Pr, ζ, n, R, and Sc are the Prandtl number, the viscosity parameter, the concentration diffusivity parameter, the radiation parameter, and the Schmidt number, respectively.


## 3. Symmetries of the Problem

In this section, we apply the Lie group of transformations which leaves the equations (11)–(13) invariant. Firstly, we consider Lie group of transformations with independent variables x and y and dependent variables ψ, Θ, and Φ for the problem
(16)x*=x*(x,y,ψ,θ,ϕ;ε),y*=y*(x,y,ψ,θ,ϕ;ε),ψ*=ψ*(x,y,ψ,θ,ϕ;ε),Θ*=Θ*(x,y,ψ,θ,ϕ;ε),Φ*=Φ*(x,y,ψ,θ,ϕ;ε),
where ε is the group parameter. The infinitesimal generator of the group (16) can be expressed in the following vector form:
(17)Γ=X∂∂x+Y∂∂y+Ψ∂∂ψ+Θ∂∂θ+Φ∂∂ϕ,
in which X, Y, Ψ, Θ, and Φ are infinitesimal functions of the group variables. Then, the corresponding one-parameter Lie group of transformations is given by
(18)x*=eεχ(x)=x+εX(x,y,ψ,θ,ϕ)+O(ε2),y*=eεχ(y)=y+εY(x,y,ψ,θ,ϕ)+O(ε2),ψ*=eεχ(ψ)=ψ+εΨ(x,y,ψ,θ,ϕ)+O(ε2),Θ*=eεχ(θ)=Θ+εΘ(x,y,ψ,θ,ϕ)+O(ε2),Φ*=eεχ(ϕ)=Φ+εΦ(x,y,ψ,θ,ϕ)+O(ε2).The action of the infinitesimal generatorΓ is extended to all derivatives appearing in equations (11)–(13) through the third prolongation
(19)Γ(3)=Γ+Ψ[x]∂∂ψx+Ψ[y]∂∂ψy+Θ[x]∂∂θx+Θ[y]∂∂θx+Φ[x]∂∂ϕx+Φ[y]∂∂ϕy+η1(xy)∂∂ψxy+Ψ[yy]∂∂ψyy+Θ[yy]∂∂θyy+Φ[yy]∂∂ϕyy+Ψ[yyy]∂∂ψyyy,
where Ψ[x],Ψ[y],Θ[x],…,Ψ[yyy] are the infinitesimal functions corresponding the derivatives. The infinitesimal generator Γ is a point symmetry of (11)-(12) if
(20)Γ(3)(Gi)∣Gi=0=0,i=1,23.
Since the coefficients of Γ do not involve derivatives, we can separate (20) with respect to derivatives and solve the result over a determined system of linear homogeneous partial differential equations known as the determining equations. After straightforward calculations for the determining equations, we find the solutions X, Y, Ψ,  Θ, and Φ of the determining equations in the form
(21)X=(c1-c2)x-c3,Y=-c2y+g(x),Ψ=c4+c2ψ,Θ=Φ=0,
where cj (j=1,2,3,4) are arbitrary constants and g(x) is the infinite parameter Lie group transformation, and the functions satisfy the following ordinary differential equations:
(22)[(c1-c2)x+c3]dB0(x)dx-c2B0(x)=0,(23)[(c1-c2)x+c3]dQ(x)dx-2c2Q(x)=0,(24)[(c1-c2)x+c3]dK(x)dx-2c2K(x)=0,
which directly give the form for the functions B(x), Q(x), and K(x),
(25)B(x)=A1(x+c3c1-c2)c2/(c1-c2),Q0(x)=A2(x+c3c1-c2)2c2/(c1-c2),K(x)=A3(x+c3c1-c2)2c2/(c1-c2),
where A1, A2,  and A3 are arbitrary constants and c1≠c2.Imposing the restrictions from boundaries and from the boundary conditions on the infinitesimals, we obtaing(x)=c4=0.What remains is to require invariance of the data which must be held on the boundary surface. This requirement means(26)Γ(2)[ψy-Uw(x)-λ*λ0(x)[1+ζ(1-θ)]ψyy]=0,whenψy(x,0)=Uw(x)+λ*λ0(x)[1+ζ(1-θ(x,0))]×ψyy(x,0),Γ(1)[ψx-V*Vw(x)]=0,whenψx(x,0)=V*Vw(x),Γ(1)[θy+h*hf(x)[1-θ]]=0,whenθy(x,0)=-h*hf(x)[1-θ].Examining the aforementioned conditions, we obtain the following differential equations:(27)[(c1-c2)x+c3]dλ0(x)dx-c2λ0(x)=0,[(c1-c2)x+c3]dUw(x)dx-(c1+c2)Uw(x)=0,[(c1-c2)x+c3]dVw(x)dx-c2Vw(x)=0,[(c1-c2)x+c3]dhf(x)dx-c2hf(x)=0,
which have the following solutions:
(28)Uw(x)=A4(x+c3c1-c2)(c1+c2)/(c1-c2),λ0(x)=A5(x+c3c1-c2)-c2/(c1-c2),Vw(x)=A6(x+c3c1-c2)c2/(c1-c2),hf(x)=A7(x+c3c1-c2)c2/(c1-c2),
where A4, A5, A6, and A7 are arbitrary constants.Having determined the infinitesimalsX, Y, Ψ, Θ, and Φ, the similarity variables are found by solving the characteristic equations
(29)dx(c1-c2)x+c3=dy-c1y=dψc1ψ=dθ0=dϕ0.
Solving (29), we have the following similarity transformations (absolute invariants):
(30)η=y(x+c3c1-c2)c2/(c1-c2),ψ=(x+c3c1-c2)c1/(c1-c2)f(η),θ=θ(η),ϕ=φ(η).Substituting the transformation (30) into (11)–(15) yields the following nonlinear system of ordinary differential equations for f(η), θ(η), and φ(η):
(31)(1+ζ(1-θ))f′′′+c1c1-c2ff′′-c1+c2c1-c2f′2-ζθ′f′′-Mf′=0,(1+4R3)θ′′+Pr(c1c1-c2fθ′+Qθ)=0,1Sc(1+nφ)φ′′+c1c1-c2fφ′+nScφ′2-K0φ=0.
The boundary conditions now become
(32)f′=1+λ[1+ζ(1-θ)]f′′,f=fw,θ′=-γ[1-θ],φ′=1,atη=0,f′=0,θ=0,φ=0,asη⟶∞,
where Q=γ*A2 is the heat generation/absorption parameter, M=A1M* is the magnetic parameter, K0=α*A3 is the chemical reaction parameter, λ=λ*A5  is the slip parameter,  γ=h*A7 is the Boit number, and fw=V*A6 is the suction/injection parameter and  A4=1.


### 3.1. Particular Cases

(i)
Forζ=R=M=Q=K0=c2=c3=0,   (3) and boundary conditions (32) reduce to (23)–(25) and boundary conditions (26) of Hamad [24].
(ii)
Forζ=n=Q=K0=M=c3=λ=0 and c2=(-1/2)c1, (31) and boundary conditions (32) reduce to (17)–(19) and boundary conditions (20)  of Ferdows  [44] when Gr=Gc=m=M=0.


### 3.2. Physical Quantities of Interest

The parameters in which we are interested for our problem are the local skin friction coefficientfx, the local Nusselt number Nux, and the local Sherwood number Shx, which are defined, respectively, as
(33)12Cfx(Rex)1/2=-(1+ζθ(0))f′′(0),Nux(Rex)-1/2=-θ′(0),Shx(Rex)-1/2=-φ′(0),
where Rex=(x-u-/ν) indicates the local Reynolds number.


## 3.1. Particular Cases

(i)
Forζ=R=M=Q=K0=c2=c3=0,   (3) and boundary conditions (32) reduce to (23)–(25) and boundary conditions (26) of Hamad [24].
(ii)
Forζ=n=Q=K0=M=c3=λ=0 and c2=(-1/2)c1, (31) and boundary conditions (32) reduce to (17)–(19) and boundary conditions (20)  of Ferdows  [44] when Gr=Gc=m=M=0.


## 3.2. Physical Quantities of Interest

The parameters in which we are interested for our problem are the local skin friction coefficientfx, the local Nusselt number Nux, and the local Sherwood number Shx, which are defined, respectively, as
(33)12Cfx(Rex)1/2=-(1+ζθ(0))f′′(0),Nux(Rex)-1/2=-θ′(0),Shx(Rex)-1/2=-φ′(0),
where Rex=(x-u-/ν) indicates the local Reynolds number.


## 4. Numerical Solution and Discussion

The transformed system of nonlinear ordinary differential equations (31) with the boundary conditions (32) has been solved numerically by using a fourth-order Runge-Kutta integration along with shooting techniques. In order to verify the accuracy of our results, we have compared the computed values of -f′′(0) with the previously published works by Hamad et al. [24], and Hayat et al. [21] with fw=0 and γ→∞ at different values of λ as shown in Table 1. The results of -f′′(0), -θ′(0) and -φ′(0) for Pr=0.72, Sc=0.22, n=0.2, γ=0.5, and fw=-0.5 for various values of λ are compared with those obtained by Hamad et al. [24] and are listed in Table 2. Moreover, Table 3 represents the comparison of -φ′(0) between our results and the results obtained by Hamad et al. [24] for Pr=0.72, λ=0.2,   γ=0.5,   and fw=0.5 for various values of n. Also, the numerical results of -f′′(0) have been compared with those obtained by Ferdows et al. [44] for M=0at different values of fw as shown in Table 4. In all comparisons, it is found that our results are in excellent agreement with the previous published results. Further, the values of the local skin friction, the local Nusselt number, and the local Sherwood number for different values of γ, K0, Q, λ, M, and ζ are listed in Table 5. It is noted from this table that the local skin friction increases with increasing ζ, heat absorption (Q<0), and M, while it decreases by increasing γ,λ, and heat generation (Q>0). The local Nusselt number decreases with the increase of λ,   M and heat generation (Q>0), while the local Nusselt number increases with the increase of ζ, heat absorption (Q<0), and γ. The local Sherwood number increases with the increase of ζ, heat absorption (Q<0), and K0, while it decreases with the increasing of γ, λ, M, and heat generation (Q>0).Table 1
Comparison of skin friction coefficient-f′′(0) for various values of slip parameter  λ  when  ζ=R=M=R=M=Q=K0=c2=c3=fw=0, and  γ→∞.
λ
Hayat et al. [21]
Hamad  et al. [24]
Present work
0
1
1
1
0.1
0.87208247
0.87208247
0.87208247
0.2
0.776377
0.77637707
0.776377
0.5
0.591195
0.59119548
0.591195
1
0.430162
0.43015970
0.430162
2
0.283981
0.28397959
0.283981
5
0.144841
0.14484019
0.14484
10
0.081249
0.08124198
0.081249
20
0.043782
0.04378834
0.043782
50
0.018634
0.01859623
0.018634
100
0.009581
0.00954997
0.009581Table 2
Comparison of values of-f′′(0),-θ′(0), and -φ′(0) for  ζ=R=M=R=M=Q=K0=c2=c3=0,  Pr=0.72, fw=-0.5, Sc=0.22, n=0.2, γ=0.5, and  ηmax=20.
λ
-
f
′
′
(
0
)
-
θ
′
(
0
)
-
φ
′
(
0
)
Hamad et al. [24]
Present work
Hamad et al. [24]
Present work
Hamad et al. [24]
Present work
0
0.78077641159
0.780776
0.1738614482
0.173859
0.1148467203
0.1148467
0.5
0.4925323862
0.492532
0.1478033629
0.147789
0.0972062425
0.0972061
1
0.3683505453
0.368350
0.1317723331
0.131737
0.0878003658
0.0878003
1.5
0.2969706178
0.296970
0.1201962831
0.120196
0.0815754873
0.0815754Table 3
Values of-φ′(0) for  ζ=R=M=R=M=Q=K0=c2=c3=0, Pr=0.72, fw=0.5, λ=0.2,  γ=0.5, and  ηmax=20.
n
Sc
=
0.22
Sc
=
0.78
Hamad et al. [24]
Present work
Hamad et al. [24]
Present work
0
0.2357864723
0.235787
0.7103150065
0.710316
0.5
0.1622076843
0.162208
0.4969259049
0.496926
1
0.1243641685
0.124365
0.3857860094
0.385787
1.5
0.1013742450
0.101374
0.3166941114
0.316694Table 4
Values of-f′′(0) for various values of  fw for ζ=n=Q=K0=c3=λ=0, c2=(-1/2)c1,   and  M=0.
f
w
Cortell [12]
Ferdows et al. [44]
Present work
-
0.75
0.453521
0.453523
0.453523
-
0.5
0.518869
0.518869
0.51887
0
0.677647
0.677648
0.677648
0.5
0.873627
0.873643
0.87363
0.75
0.984417
0.984439
0.98444Table 5
Values of(-1/2)Cfx(Rex)1/2,
-Nux(Rex)-1/2, and -Shx(Rex)-1/2 for fw=0.5,  Pr=0.72, Sc=0.6, n=0.5, and R=1.
ζ
K
0
γ
λ
M
Q
(
-
1
/
2
)
C
f
x
(
Re
x
)
1
/
2
-
Nu
x
(
Re
x
)
-
1
/
2
-
Sh
x
(
Re
x
)
-
1
/
2
0
0.01
1
0.3
0.5
-
0.2
0.8282097
0.267065009
0.2817571148
2
0.01
1
0.3
0.5
-
0.2
0.98684111
0.272061144
0.30316651
5
0.01
1
0.3
0.5
-
0.2
1.16050976
0.275202928
0.31591854
0.1
0
1
0.3
0.7
-
0.5
0.89348981
0.329592869
0.263842318
0.1
1
1
0.3
0.7
-
0.5
0.89348981
0.329592869
0.70701644
0.1
3
1
0.3
0.7
-
0.5
0.89348981
0.329592869
1.128126581
1
0.3
0.1
0.3
0.5
-
0.2
0.9886912570
0.07894766720
0.4729295509
1
0.3
0.5
0.3
0.5
-
0.2
0.938324634
0.213201753
0.46994994
1
0.3
1
0.3
0.5
-
0.2
0.913264332
0.270045367
0.468340615
0.1
0.5
0.1
0
0.7
-
0.2
1.29926459
0.079169100
0.56364235
0.1
0.5
0.1
1
0.7
-
0.2
0.53978525
0.077353309
0.52577635
0.1
0.5
0.1
2
0.7
-
0.2
0.34699071054
0.07673167148
0.5143743566
0.1
0.6
0.7
0.5
0
-
1
0.5697224247
0.08671881627
0.5920198931
0.1
0.6
0.7
0.5
1
-
1
0.80801328749
0.08644826115
0.5705347139
0.1
0.6
0.7
0.5
3
-
1
1.0116231621
0.08626331255
0.5563999142
0.7
0.5
0.8
1
0.5
-
0.8
0.5409262133
0.335939548322
0.52997592226
0.7
0.5
0.8
1
0.5
-
0.3
0.535413351
0.265641294131
0.52980646921
0.7
0.5
0.8
1
0.5
0
0.5252295045
0.15291844567
0.529531543678
0.7
0.5
0.8
1
0.5
0.02
0.52331016216
0.13341074189
0.52948110939
0.7
0.5
0.8
1
0.5
0.1
0.51234845001
0.037287564758
0.52938666933The profiles for velocity, temperature and concentration are shown in Figures2–23, respectively, with various values of the parameters. Figures 2, 3, and 4 show the dimensionless velocity profiles f′(η), the dimensionless temperature profiles and the dimensionless mass profiles in the boundary layer for various values of the viscosity parameter ζ. These figures depict that increase in the value of ζ the dimensionless velocity decreases near the surface but increases as larger distances, while the dimensionless temperature and the dimensionless mass in the boundary layer region decrease as the viscosity parameter increases.Figure 2
Velocity distribution for various values ofζ.Figure 3
Temperature distribution for various values ofζ.Figure 4
Mass distribution for various values ofζ.Figure 5
Velocity distribution for various values offw.Figure 6
Temperature distribution for various values  offw.Figure 7
Mass distribution for various values  offw.Figure 8
Velocity distribution for various values ofγ.Figure 9
Temperature distribution for various values ofγ.Figure 10
Mass distribution for various values  ofK0.Figure 11
Velocity distribution for various values ofλ.Figure 12
Temperature distribution for various values  ofλ.Figure 13
Mass distribution for various values  ofλ.Figure 14
Velocity distribution for various values ofM.Figure 15
Temperature distribution for various values ofM.Figure 16
Mass distribution for various values ofM.Figure 17
Mass distribution for various values ofn.Figure 18
Velocity distribution for various values ofPr.Figure 19
Temperature distribution for various values ofPr.Figure 20
Velocity distribution for various values ofQ.Figure 21
Temperature distribution for various values  ofQ.Figure 22
Velocity distribution for various values ofR.Figure 23
Temperature distribution for various values  ofR.Figures5–7 exhibit the effect of the suction or injection parameter fw on the dimensionless velocity profiles, the dimensionless temperature, and the dimensionless mass profiles. It is shown from Figure 5 that the suction parameter decreases the velocity indicating the usual fact that suction stabilizes the boundary layer growth, while injection increases the velocity in the boundary layer region indicating that injection helps the flow penetrate more into the fluid. In Figure 6, it is found that the temperature decreases as the suction parameter increases. This means that larger suction leads to faster cooling of the plate, while the temperature increases as the injection parameter increases; this is because of the heat transfer from fluid to surface. Figure 7 illustrates that the mass decreases as the suction parameter increases and increases as the injection parameter increases.The dimensionless velocity profiles for different values of the Boit numberγ are described in Figure 8. From this figure, we observe that an increase in γ leads to increase in f′ near the wall but it decreases with larger distances, while Figure 9 shows that the dimensionless temperature decreases as the Boit number γ increases.Effect of the chemical reaction parameterK0 on the dimensionless concentration profiles is shown in Figure 10. It is seen that the dimensionless concentration of the fluid decreases with the increase of K0.Figures11, 12, and 13 illustrate the influence of slip velocity parameter λ on the velocity, temperature, and concentration profiles, respectively. It is obvious that slip parameter decreases the velocity f′ and increases both dimensionless temperature θ and the dimensionless concentration φ.Figures14, 15, and 16 show the velocity, temperature, and concentration distributions for various values of the magnetic field M. We can notice that increasing M decreases the velocity f′ and increases the temperature θ and concentration; this result qualitatively agrees with the expectations; this is because the application of a transverse magnetic field to an electrically conducting fluid gives rise to a resistive-type force called the Lorentz force. This force has the tendency to slow down the motion of the fluid in the boundary layer and to increase its temperature and concentration.Figure17 illustrates the effect of the concentration diffusivity parameter n on the concentration profiles; one can observe that the concentration profiles increase as n increases.The effects of the Prandtl numberPr on the dimensionless velocity and temperature dimensionless velocity are illustrated in Figures 18 and 19, respectively. From Figure 18, it can be seen that as the Prandtl number Pr increases, the dimensionless velocity decreases near the surface but increases with larger distances. Figure 19 illustrates that the dimensionless temperature decreases with increasing Pr; this is because, when the Prandtl number increases the thickness of the thermal boundary layer decreases and, hence, the temperature decreases.The effects of the presence of heat source (Q>0) or heat sink (Q<0) in the boundary layer on the velocity and temperature profiles are presented in Figures 20 and 21, respectively. From Figure 20, one sees that the dimensionless velocity decreases with the heat source (Q>0) increasing, while it increases with the heat sink (Q<0) increasing. The presence of heat source in the boundary layer generates energy which causes the temperature of the fluid to increase. On the other hand, the presence of heat sink in the boundary layer absorbs energy which causes the temperature of the fluid to decrease, and this corresponds with the observation in Figure 21, in which we note that the dimensionless temperature increases with the heat source (Q>0) increasing, while it decreases with the heat sink (Q<0) increasing.For different values of radiation parameterR, the dimensionless velocity profiles are plotted in Figure 22. It is obvious that velocity decreases with the increase in R. The effects on the temperature profiles are presented in Figure 23. From this figure, we observed that, as the value of R increases, the temperature profiles increase, with an increase in the thermal boundary layer thickness.


## 5. Conclusions

In this study, the symmetries of the problem of heat and mass transfer flow over a moving permeable flat stretching sheet in the presence of thermal radiation, variable viscosity, a uniform transverse magnetic field, chemical reaction, heat generation/absorption, suction/injection, convective boundary condition, and slip velocity are obtained using Lie group analysis. Numerical solutions of the resulting system of nonlinear ordinary differential equations are obtained by using the shooting method coupled with Runge-Kutta scheme. The effects of various parameters on the dimensionless velocity, temperature, and concentration profiles have been studied graphically. The influence of some parameters on the local skin friction, local Nusselt number, and local Sherwood number is tabulated. It was found that with the increase of the convective parametern, injection parameter fw<0, and the slip velocity λ, the absolute concentration increases, while it decreases with the increase in the magnetic parameter M, viscosity parameter ζ, suction parameter fw>0, and chemical reaction parameter K0.



---

*Source: 359817-2013-09-26.xml*