Effect of Chemical Reaction On Mass Transfer Over A Stretching Surface Embedded in A Porous Medium

Download as pdf or txt
Download as pdf or txt
You are on page 1of 9

International Journal of Computational Engineering Research||Vol, 04||Issue, 2||

||Issn 2250-3005 || ||February||2014|| Page 20


Effect of Chemical Reaction on Mass Transfer over a Stretching
Surface Embedded In a Porous Medium

1,
E. M. A. Elbashbeshy ,
2,
A. M. Sedki
1
,Mathematics Department, Faculty of science, Ain Shams University, Abbassia, Cairo, Egypt,
2,
Mathematics
Department, Faculty of science Jazan University, Jazan, KSA








I. INTRODUCTION
The boundary layer flows due to stretching surface in porous medium are relevant to many
engineering problems such as paper production, preparing plastic and metal sheets etc.The dynamics of the
boundary layer flow over a stretching surface originated from the pioneering work of Sakiadis [1, 2] who
initiated the study of boundary layer flow over a continuous solid surface moving with constant speed. Crane
[3] extended it to analyze the steady two dimensional boundary layer flow caused by the stretching of elastic
flat surface which moves in its plane with velocity varying linearly with distance from a fixed point. Many
authors presented some mathematical results, and good amount of references can be found in the papers by Ali
[4] and [5], Elbashbeshy [6], Ishak et al. [7] and Elbashbeshy and Bazid [8] .The studies carried out in these
papers in the case steady state flow. The unsteady state problem over a stretching surface, which is stretched
with a velocity that depends on time is considered by Anderson et al. [9], Elbashbeshy and Bazid [10] and
Ishak et al. [11]. The effects of chemically reactive solute distribution on fluid flow due to a stretching surface
also bear equal importance in engineering researches. The chemical reaction effects were studied by many
researchers on several physical aspects. The diffusion of a chemically reactive species in a laminar boundary
layer flow over a flat plate was demonstrated by Chambre and Young [12]. The effect of transfer of
chemically reactive species in the laminar flow over a stretching sheet explained by Andersson et al. [13].
Takhar et al. [14] analyzed the flow and mass transfer on a stretching sheet with a magnetic field and
chemically reactive species with n-th order reaction.

The mass transfer in boundary layer flow due to a stretching surface in porous medium also has
important applications in many industrial problems. The effect of mass transfer in laminar flow over a
stretching surface was investigated by Radwan and Elbashbeshy [15]. Radwan and Elbashbeshy [15] analyzed
the flow and mass transfer on a stretching surface with a magnetic field. Akyildiz et al. [16] reported a solution
for diffusion of chemically reactive species in a flow of a non- Newtonian fluid over a stretching sheet
immersed in a porous medium. El-Aziz [17] explained unsteady flow due to a stretching sheet with mass and
heat transfer. Recently, Krishnendu [18] studied the boundary layer flow with first order chemical reaction
over a porous flat plate. Krishnendu [19] studied the mass transfer on a continuous flat plate moving in a
parallel or reversely to a free stream in the presence of a chemical reaction. Ferdows et al. [20] investigated the
effects of order of chemical reaction on mass transfer over a linearly stretching surface. Based on the above-
mentioned investigations and applications, this paper is concerned with two-dimensional steady,
incompressible, laminar boundary layer flow of a fluid over a linearly stretching surface.
In this paper we investigate numerically the effects of chemical reaction on the steady laminar two-
dimensional boundary layer flow and mass transfer over a stretching surface embedded in porous medium.
The method of solutions based on the well-known similarity analysis together with shooting method.
ABSTRACT
An analysis is made to investigate the mass transfer over a stretching surface embedded in a porous
medium in the presence of first order chemical reaction. Using similarity transformation, the
governing partial differential equations are transformed into a set of ordinary differential equations
which solved by shooting method. Comparisons with previously published work on special cases of
the problem are performed and the results are found to be in excellent agreement. It is observed that
the local mass transfer and concentration profile are very sensitive to change in the values of
reaction rate parameter, permeability parameter and Schmidt number.

KEYWORDS: Mass transfer Stretching surface Porous medium- Chemical reaction.


Effect Of Chemical Reaction On Mass Transfer Over
||Issn 2250-3005 || ||February||2014|| Page 21
II. Formulation of the problem
Consider two dimensional steady, laminar boundary layer flow of a fluid over a linearly stretching surface
(i.e. stretched with a velocity proportional to x) embedded in porous medium with velocity
w
u and
concentration
w
C moving axially through a stationary fluid. The fluid is viscous incompressible with constant
physical properties. We assume that the fluid far away from the surface is at rest and at concentration C. .
The x axis runs along the continuous surface in the direction of motion and y axis is perpendicular to it.
The continuity, momentum and reactive concentration equations for governing the flow and concentration
distribution in the boundary layer region along the stretching surface may be written as

) 1 ( 0 =
c
c
+
c
c
y
v
x
u

) 2 (
2
2
u
k y
u
y
u
v
x
u
u
v
v
c
c
=
c
c
+
c
c

) 3 ( ) (
2
2


c
c
=
c
c
+
c
c
C C R
y
C
D
y
C
v
x
C
u

The boundary conditions are given by

= =
+ = = = =
C C u y
Ax C C v x u y
n
, 0 :
) 4 (
, 0 , : 0

where u and v are the velocity components in the x and y directions, respectively, k is the permeability of
porous medium. C is the concentration of the fluid, is the kinematic viscosity, R is a constant of first-order
chemical reaction rate, n is a power-law exponent, which signifies the change of amount of solute in the x-
direction, D is the effective diffusion coefficients and and A are constants. The equation of
continuity is satisfied if we choose a dimensionless stream function ) , ( y x such that
. ,
x
v
y
u
c
c
=
c
c
=


In order to solve Equations (1)-(4), we introduce the following similarity transformation and dimensionless
concentration.

= = =
C C
C C
n x
y x
f
n
y
w
) ( ,
) , (
) ( , q u
v

q
v
q
whereq is the similarity variable, ) , ( y x is the dimensionless stream function, ) (q u is the dimensionless
concentration. Consequently, equations (2)-(3) and the boundary conditions (4) can be written in the following
form
0 , 0 :
) 7 (
1 , 1 , 0 : 0
) 6 ( 0 ) (
) 5 ( 0
/
/
/ / //
/ 2 / // ///
= =
= = = =
= +
= +
u q
u q
u u u u
f
f f
ScR nf f Sc
Kf f ff f
C

Effect Of Chemical Reaction On Mass Transfer Over
||Issn 2250-3005 || ||February||2014|| Page 22
Where prime denotes differentiation with respect toq .
nk
K
v
= is a permeability parameter,
n
R
R
C
= is a
reaction rate parameter of the solute and
D
Sc
v
= is Schmidt parameter.
The mass transfer coefficient in terms of the Sherwood number Sh is given by

) 0 (
) 0 (
) (
/
/
u
u
v
=
=

c
c
=

e
w
R
Sh
n
x
C C
y
C
x
Sh

where
v
x u
R
w
e
= is the local Reynolds number.
III. Numerical method
The above Eqs. (5) and (6) along with boundary conditions (7) are solved by converting them to an initial
value problem. We set
q p z z f = = =
/ / /
, , u
fp z Kz p + =
2 /
(8)
] [
/
fq nz R Sc q + = u u (9)
with the boundary conditions
1 ) 0 ( , 1 ) 0 ( , 0 ) 0 (
/
= = = u f f (10)
In order to integrate (8) and (9) as an initial value problem we require a value for ) 0 ( p i.e. ) 0 (
//
f
and ) 0 (
/
u but no such values are given in the boundary. The most important factor of shooting method is to
choose the appropriate finite values of

q . In order to determine

q for the boundary value problem stated by


Eqs. (8) (9), we start with some initial guess value for some particular set of physical parameters to obtain
) 0 (
//
f and ) 0 (
/
u . The solution procedure is repeated with another large value of

q until two successive


values of ) 0 (
//
f and ) 0 (
/
u differ only by the specified significant digit. The last value of

q is finally chosen
to be the most appropriate value of the limit

q for that particular set of parameters. The value of

q may
change for another set of physical parameters. Once the finite value of

q is determined then the integration is


carried out. We compare the calculated values for
/
f and u at 10 = q (say) with the given boundary
conditions 0 ) 10 (
/
= f and 0 ) 10 ( = u and adjust the estimated values, ) 0 (
//
f and ) 0 (
/
u , to give a better
approximation for the solution. We take the series of values for ) 0 (
//
f and ) 0 (
/
u , and apply the fourth order
classical RungeKutta method with step-size h= 0.01. The above procedure is repeated until we get the results
up to the desired degree of accuracy,
5
10

.



Effect Of Chemical Reaction On Mass Transfer Over
||Issn 2250-3005 || ||February||2014|| Page 23
IV. RESULTS AND DISCUSSION
The set of non-linear ordinary differential equations (5) and (6) satisfying the boundary conditions (7)
have been solved numerically using the Mathematica method for several values of the involved parameters,
namely permeability parameter, power-law exponent, reaction rate parameter of the solute and Schmidt
parameter. In case n=0 and 0 = K our results are very similar to those considered by Uddin et al. [21], Takhar
et al. [14] and Andersson et al. [13]. To validate the numerical method used in this study, the case n=0 and
0 = K , was considered and the results for concentration gradient are compared with the numerical solution
which reported in Uddin et al. [21], Takhar et al. [14] and Andersson et al. [13]. The quantitative comparison is
shown in Table 1, and found to be in excellent agreement.

Table1 Comparison of the values of ) 0 (
/
u with that of Takhar et al. [14], Andersson et al. [13] and Uddin et
al. [21] for n=0 and 0 = K .

Sc
C
R
Present Study

Takhar et al.
[14]
Andersson et al. [13]

Uddin et al.[21]
0.1 0.1 0.149083 0.15042 0.14900 0.15057
1 0.1 0.668754 0.67044 0.66900 0.66873
1 1 1.176401 1.17761 1.17700 1.17679
10 1 3.871327 3.87469 3.88000 3.87347
10 10 10..241185 10..24283 10..25000 10..24535

0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
-0.2
0
0.2
0.4
0.6
0.8
1
1.2
n
(
d

f
)
K = 0 , 1 , 4 , 10
Sc = 1 , Rc = 1 , n = 1
q

Fig. 1: Velocity profiles f () for various values of K

/
f
Effect Of Chemical Reaction On Mass Transfer Over
||Issn 2250-3005 || ||February||2014|| Page 24
0 1 2 3 4 5
-0.2
0
0.2
0.4
0.6
0.8
1
1.2
n
(
p
h
i
)
K = 10 , 4 ,1 , 0
Sc = 1 , Rc = 1 , n = 1



Fig.2: Concentration profiles ) (q u for various values of K

The numerical solution of velocity has presented for various values of the permeability parameter K . The
reactant solute equation is solved numerically and the results are shown graphically. The velocity profiles for
various values of the permeability parameter K have been plotted in Fig. 1. From the figure it is noted that with
increase of K , the velocity for any fixed value of decreases.

0 0.5 1 1.5 2 2.5 3 3.5 4
-0.2
0
0.2
0.4
0.6
0.8
1
1.2
n
(
p
h
i
)
Sc = 1 , 2 , 5 , 11
K = 1 , Rc = 1 , n = 1



Fig. 3. Concentration profiles ) (q u for various values of Sc
q

q


Effect Of Chemical Reaction On Mass Transfer Over
||Issn 2250-3005 || ||February||2014|| Page 25
0 1 2 3 4 5 6
-0.2
0
0.2
0.4
0.6
0.8
1
1.2
n
(
p
h
i
)
Rc = 0 , 1 , 4 , 12
K= 1 , Sc = 1 , n = 1

Fig4. Concentration profiles ) (q u for various values of
C
R

Thus it is clear that the permeability parameter K opposes motion. Consequently, the momentum
boundary layer thickness reduces with the increase K in and this fact is also seen from wall shear stress
behaviour. Exhibits concentration profiles for various values of K . The value of contaminate solute at
particular value of increases with the increase of the permeability parameter K and also the concentration
boundary layer thickness increases. This implies that the permeability parameter acts to enhance the distribution
of the reaction solute on the stretching surface. This result may be useful, in the situation where the
enhancement of solute transfer from the surface is the prime important.

0 0.5 1 1.5 2 2.5 3
-0.2
0
0.2
0.4
0.6
0.8
1
1.2
n
(
p
h
i
)
n = 0 , 2 , 6 , 10
K =1 , Sc = 1 , Rc = 1

Fig. 5. Concentration profiles ) (q u for various values of n(0). with K=Sc=
C
R =1
q
q


q


Effect Of Chemical Reaction On Mass Transfer Over
||Issn 2250-3005 || ||February||2014|| Page 26
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
-0.2
0
0.2
0.4
0.6
0.8
1
1.2
1.4
n
(
p
h
i
)
n = - 6 , - 5 , - 4 , - 3 , - 1
K = 1 , Rc = 1 , Sc = 1


Fig. 6. Concentration profiles ) (q u for various values of n(<0). with K=Sc=
C
R =1

Now, we concentrate on variation the solute curves for different values of Schmidt number Sc. The
curves are drawn in the Fig. 3. The Schmidt number has major effects on the distribution of solute. The
concentration boundary layer thickness as well as the concentration at a fixed point decreases quickly with
increasing values of Sc. This is due to the fact that the rate of solute transfer from the surface increases when the
Schmidt number increases.Figure 4 is the graphical representation of concentration profiles for various values of
reaction rate parameter
C
R . It has been found that the reaction rate parameter affect the solute profiles in similar
way as that of the Schmidt number i.e., the increase of
C
R reduces both the solute boundary layer thickness and
value of the solute at fixed . So, in case of the distribution of reactive solute, the reaction rate parameter is a
decelerating agent.

Finally, Figs. 5 and 6 exhibit the concentration profiles in the boundary layer flow region for different
values of power-law exponent n. It is noticed from Fig. 5 that for the increasing values of n with n>0, the curve
representing the distribution of solute for specific value of decreases. While, in Fig. 6 the concentration profile
increases with increase in the magnitude of n with n<0 and for large negative values of n, the overshoot of solute
is observed near the surface. Thus, the effect of increase of n when the surface concentration is
n
w
Ax C C + =

is completely opposite to the effect of increase n when the surface concentration
is
n
w
x A C C / + =

where A is constant positive. Note that, the wall concentration is constant when n=0.

Table 2
Values of skin-friction coefficient ) 0 (
//
f and - ) 0 (
/
u for various values of permeability parameter K .

K
) 0 (
//
f ) 0 (
/
u -
0 -1.00000 1.45942
1 -1.414214 1.414214
4 -2.236068 1..3433
10 -3.316625 1..2781


q


Effect Of Chemical Reaction On Mass Transfer Over
||Issn 2250-3005 || ||February||2014|| Page 27
Table 3:Values of - ) 0 (
/
u for various values of Schmidt number Sc

Sc
) 0 (
/
u -
1 1.414214
2 2.063933
5 3..366092
11 5.085476

Table 4:Values of - ) 0 (
/
u for various values of reaction rate parameter
C
R

C
R
) 0 (
/
u -
0 0.892503
1 1.414214
4 2.2692285
12 3.641534

Table 5:Values of - ) 0 (
/
u for various values of n
n
) 0 (
/
u -
0 1.154956
2 1.649948
6 2.434269
10 3.061963

Tables2-5summarize the calculated values of local skin friction coefficient ) 0 (
//
f and rate of mass transfer -
) 0 (
/
u for 10 0 , 11 1 s s s s n Sc and 12 0 s s
C
R . These tables indicate that - ) 0 (
//
f increasing as
permeability parameter K increases. Since K does not occur explicitly in the diffusion equation , it effect on
mass transfer - ) 0 (
/
u - is small. A lso, these tables indicate that - ) 0 (
/
u increasing as n,
C
R and Sc increases.

V. CONCLUSIONS
In this investigation, an analysis is made to find the behaviour of the distribution of reactive solute
undergo a first order reaction in steady boundary layer flow incompressible fluid over a stretching surface
embedded in porous medium with variable surface concentration. Using the similarity transformation a set of
ordinary differential equations is obtained from the governing equations. The momentum equation and the
equation of reactive solute are solved numerically. The results show that the permeability parameter tends to
reduce the rate of flow from the wall and is broadening the solute layer. The Schmidt number and the reaction
rate parameter reduce the solute boundary layer thickness. Most, importantly, the effects of initial variable
solute distribution over a stretching surface is interesting i.e. for the increase in magnitude of n, the
concentration decreases when n>0 whereas increases when n<0.

VI. REFERENCES
[1] Sakiadis B.C., Boundary layer behavior on continuous solid surfaces: I. Boundary layer equations for two dimensional and
axisymmetric flow, AIChE J, 7(1), 26-28, 1961.
[2] Sakiadis B.C., Boundary layer behavior on continuous solid surfaces: II. Boundary layer equations on a continuous flat surface,
AIChE J., 7(1), 221-225, 1961.
[3] Crane L.J., Flow past a stretching plane, Z. Amgew Math. Phys., 21, 645-647, 1970.
[4] Ali M.E., Heat transfer characteristics of a continuous stretching surface, Warme-Und Stoffuber tragung, 29, 227-234 ,1994 .
[5] Ali M.E., On thermal boundary layer on a power law stretched surface with suction or injection, Int. J. Heat Mass Flow, 16, 280-
290, 1995.
[6] Elbashbeshy E.M.A., Heat transfer over a stretching surface with variable heat flux, J. Phys.D:Appl. Phys., 31, 1951-1955, 1998.
[7] -Ishak A., Nazar R. and Pop I., unsteady mixed convection boundary layer flow due to a stretching vertical surface , Arabian
J.Sce. Engng. , 31,165-182, 2006.
[8] Elbashbeshy E.M.A and Bazid M.A.A.,. Heat transfer over a continuously moving plate embedded in non-Darcian porous
medium, Int. J. Heat and Mass Transfer, 43, 3087-3092, 2000.
[9] Andersson H.T., Aarseth J.B. and Dandapat B.S., Heat transfer in a liquid film on an unsteady stretching surface, Int .J .Heat
Transfer, 43, 69-74 ,2000 .
[10] Elbashbeshy E.M.A and Bazid M.A.A., Heat transfer over an unsteady stretching surface, Heat Mass Transfer, 41, 1-4, 2004.
Effect Of Chemical Reaction On Mass Transfer Over
||Issn 2250-3005 || ||February||2014|| Page 28
[11] Ishak A., Nazar R. and Pop I, Heat transfer over an unsteady stretching surface with prescribed heat flux, Can. J. of Phys.86, 853-
855 ,2008 .
[12] Chambre, P.L. and Young J.D., On diffusion of a chemically reactive species in a laminar boundary layer flow, Physics of
Fluids, 1, 48-54, 1958.
[13] Andersson, H.I., Hansen O.R. and Holmedal B., Diffusion of a chemically reactive species from a stretching sheet. International
Journal of Heat and Mass Transfer, 37, 659-664, 1994.
[14] Takhar, H.S., Chamkha A.J., Nath G., Flow and mass transfer on a stretching sheet with a magnetic field and chemically reactive
species. International Journal of Engineering Science, 38, 1303-1314, 2000.
[15] Radwan A.E. and Elbashbeshy E.M.A. Mass transfer over a stretching surface with variable concentration in a transverse
magnetic field IL NUOVO CIMENTO, 105B (6), 1990.
[16] Akyildiz, F.T., H. Bellout and K. Vajravelu, Diffusion of chemically reactive species in a porous medium over a stretching
sheet. Journal of Mathematical Analysis and Application 320, 322-339, (2006).
[17] El-Aziz, M.A., Unsteady fluid and heat flow induced by a stretching sheet with mass transfer and chemical reaction, Chemical
Engineering Communications 197, 1261-1272, (2010).
[18] Krishnendu Bhattacharyya, Boundary layer flow with diffusion and first order chemical reaction over a porous flat plate subject
to suction/injection and with variable wall concentration, Chemical Engineering Bulletin, 15, 6-11, 2011.
[19] Krishnendu Bhattacharyya, mass transfer on a continuous flat plate moving in a parallel or reversely to a free stream in the
presence of a chemical reaction, International Journal of Heat and Mass Transfer, 55, 3483-3487, 2012.
[20] Ferdows M. and Qasem M. Al-Mdallal, Effects of Order of Chemical Reaction on a Boundary Layer Flow with Heat and Mass
Transfer Over a Linearly Stretching Sheet, American Journal of Fluid Dynamics, 2(6), 89-94, 2012.
[21] Uddin M.S., Bhattacharyya K., Layek G.C. and Pk W.A. , Chemically Reactive Solute Distribution in a Steady MHD Boundary
Layer Flow over a Stretching Surface, Journal of Applied Fluid Mechanics, Vol. 4, No. 4, pp. 53-58, 2011.

You might also like