# Unsteady Natural Convection Flow past an Infinite Cylinder with Thermal and Mass Stratification

Authors: Ashish Paul; Rudra Kanta Deka
Journal: International Journal of Engineering Mathematics
 (2017)

---

## Abstract

This paper presents an analytical solution of unsteady one-dimensional free convection flow past an infinite vertical circular cylinder in a stratified fluid medium. The dimensionless coupled linear governing partial differential equations are solved by Laplace transform technique for unit Prandtl number and Schmidt number. Effects of various physical parameters are presented with graphs. Numerical values of boundary layer thickness for different parameters are presented in table. Due to the effects of thermal and mass stratifications, the velocity, temperature, and skin friction, Nusselt number shows oscillatory behaviour at smaller times and then reaches steady state at larger times.


---

## Body

## 1. Introduction

Natural convection flows with heat and mass stratification are frequently encountered in nature. These types of problem over vertical cylinder have wide range of applications in the field of science and technology such as startup of chemical reactors and emergency cooling of nuclear fuel elements. In glass and polymer industries, hot filaments are considered as vertical cylinder and cooled as they pass through the surrounding environment. Free convective flows driven by temperature and concentration difference have been studied extensively. When both the temperature and concentration differences occur simultaneously, the free convective flow can become quite complex.Gebhart and Pera [1] analysed the steady combined buoyancy effects on vertical natural convection flows. Bottemanne [2] studied the combined effect of heat and mass transfer in the steady laminar boundary layer of a vertical cylinder placed in still air for Pr = 0.71 and Sc = 0.63. Chen and Yuh [3] presented a numerical study of steady heat and mass transfer processes near cylinder with uniform wall heat and mass fluxes and wall temperature. Their study covered a wide range of radii and Prandtl numbers. Velusamy and Garg [4] studied natural convection adjacent to a heat generating vertical cylinder. Ganesan and Rani [5] presented a numerical solution for the transient natural convection flow over a vertical cylinder under the combined buoyancy effect of heat and mass transfer. Numerical analysis of two-dimensional unsteady natural convective flow past semi-infinite vertical cylinder with heat and mass transfer under different physical situations was studied by Ganesan and Loganathan [6–8]. Rani [9] presented a numerical study on transient natural convection along vertical cylinder with variable surface temperature and mass diffusion. Devakar et al. [10] presented closed form solution for Poiseuille flow, Couette flow, and generalized Couette flows of an incompressible couple stress fluid between two concentric circular cylinders with slip boundary condition. Recently, Deka et al. [11] presented the analytical investigation of one-dimensional unsteady natural convection flow past an infinite vertical cylinder with heat and mass transfer under the effect of constant heat flux at the surface of the cylinder. They have shown that the velocity and temperature increase unboundedly with time, while the concentration approaches steady state at larger times.In recent times many researchers have shown interests in the study of stratification effects on transient natural convective flows along vertical bodies under various physical situations. Takhar et al. [12] presented a numerical study of natural convection boundary layer flow over a continuously moving vertical surface immersed in a thermally stratified medium by an implicit finite difference scheme. Again, Takhar et al. [13] investigated the natural convection flow past a vertical cylinder embedded in a thermally stratified high-porosity medium. They solved the coupled nonlinear partial differential equations by finite difference as well as perturbation technique and found that separation of flow occurs for some values of stratification parameter. Loganathan and Ganesan [14] presented a numerical study of free convective flow of a viscous incompressible fluid past a moving, semi-infinite vertical cylinder with constant temperature and mass diffusion in a thermally stratified medium by employing a finite difference scheme of Crank-Nicolson type.Shapiro and Fedorovish [15] presented analytical solution of one-dimensional laminar natural convection along an infinite vertical plate by introducing the pressure work term and the ambient thermal stratification in the thermodynamic energy equation for the case of unit Prandtl number. They have shown that thermal stratification provides a negative feedback mechanism: warm fluid rises, expands, and cools relative to the environment, whereas cool fluid subsides, compresses, and warms relative to the environments. Later on, Shapiro and Fedorovich [16] carried out study on natural convection in a stably stratified fluid along vertical plates and circular cylinders seeking solutions in the form of harmonic oscillators.In recent times, the effect of double stratifications, namely, thermal stratification and mass stratification, has been considered by different researchers. For example, Cheng [17] studied the coupled heat and mass transfer by natural convection near a vertical wavy surface in a non-Newtonian fluid saturated porous medium with thermal and mass stratification and obtained solutions by collocation method. Recently, Srinivasacharya and Reddy [18], Srinivasacharya and RamReddy [19], Rathish Kumar and Krishna Murthy [20], and Neagu [21] have investigated numerically on flow past plates and wavy surfaces taking double stratifications into account.Deka and Paul [22, 23] presented the analytical investigation of transient free convection flow past an infinite moving vertical cylinder in a stably stratified fluid including thermal stratification by employing Laplace transform technique. Very recently Deka and Paul [24] presented analytical investigation to study the effects of thermal stratification and mass stratification on natural convection heat and mass transfer over moving vertical cylinder. This motivates undertaking this study. This paper presents an analytical investigation of one-dimensional free convective flow past a stationary infinite vertical cylinder with combined effects of thermal and mass stratification. The unsteady nondimensional governing linear equations are solved by Laplace transform technique for the case of unit Prandtl number and unit Schmidt number. Solutions are presented in closed form and this is always necessary for validating numerical models. Also solutions of unsteady state for larger time are compared with the solutions of steady state.


## 2. Mathematical Analysis

Consider an unsteady, laminar, and incompressible viscous flow past an infinite vertical cylinder of radiusr
0 with constant temperature and concentration in presence of thermal and mass stratification. The z-axis of the cylinder is taken vertically upward along the axis of the cylinder and the radial coordinate r is taken normal to the cylinder as shown in Figure 1. The physical model and coordinate system of the flow problem is shown in Figure 1. Upon commencement of the transient, we consider the fluid moving up from the leading edge (z
=
0) parallel to the cylinder as a wave, in front of which the velocity, temperature, and concentration are only functions of the time and the radial distance r from the cylinder. Behind the wave there must be a dependence on the vertical coordinate, z. The basic premise in this work is that convective effects will begin at a position, z, as soon as fluid which was initially located at the leading edge rises to this position, regardless of the distance, r, away from the cylinder at which it first arrives. Since the surface temperature and concentration above the leading edge are uniform with z, the temperature of the fluid and concentration may be assumed to be independent of z. In addition, the vertical velocity, u, must be independent of z and from the continuity equation, the velocity normal to the plate is seen to be zero, except that the temperature and concentration of the ambient fluid are function of the vertical distance z only. At time t
′
>
0, the uniform temperature (T
0) and concentration (C
0) are specified at the surface of the cylinder. Viscous dissipation terms have been neglected. All derivatives in the direction parallel to the cylinder are zero, except d
T
∞
′
/
d
z and d
C
∞
′
/
d
z termed as thermal stratification and mass stratification, respectively. Here, T
∞
′ and C
∞
′ are the temperature and concentration of the undisturbed fluid. It is to be noted that initially the fluid may not be stratified, but upon commencement of the transient the fluid gets self-stratifications. Then following Boussinesq’s approximation, the one-dimensional equations for momentum, energy, and concentration are as follows:(1)
∂
u
∂
t
′
=
ν
r
∂
∂
r
r
∂
u
∂
r
+
g
β
T
′
-
T
∞
′
+
g
β
∗
C
′
-
C
∞
′
,
∂
T
′
∂
t
′
=
α
r
∂
∂
r
r
∂
T
′
∂
r
-
γ
u
,
∂
C
′
∂
t
′
=
D
r
∂
∂
r
r
∂
C
′
∂
r
-
ξ
uwith initial and boundary conditions as(2)
t
′
≤
0
:
u
=
0
,
T
′
=
T
∞
′
C
′
=
C
∞
′
∀
r
t
′
>
0
:
u
=
0
,
T
′
=
T
0
,
C
′
=
C
0
a
t
r
=
r
0
u
⟶
0
,
T
′
⟶
T
∞
′
,
C
′
⟶
T
∞
′
a
s
r
⟶
∞
.Figure 1
Sketch of the geometry.The physical variables involved above are mentioned in Nomenclature. However, the particular variableγ, the main concern in our study, is the ambient thermal stratification parameter [
≡
d
T
∞
′
(
z
)
/
d
z for Boussinesq flow of liquids or gases, ≡
d
T
∞
′
(
z
)
/
d
z
+
g
/
C
p for a perfect gas with pressure work term retained], and z is height. Thus γ is the combination of thermal stratification and compression, where d
T
∞
′
(
z
)
/
d
z stands for thermal stratification and g
/
C
p for compression. In an adiabatic environment, d
T
∞
′
z
/
d
z
+
g
/
C
p
=
0. In that case d
T
∞
′
(
z
)
/
d
z
=
-
g
/
C
p and (-
g
/
C
p) is called the adiabatic temperature gradient and is the largest rate at which the temperature can decrease with height without causing instability. The stability of the atmosphere is determined according to γ
>
0 (stable), γ
=
0 (neutral), and γ
<
0 (unstable). It should be noted, however, that the compression work term is generally quite small (since g = 9.8 m s−2, C
p = 1004 J kg−1 K−1) and that the main interest in our solutions will probably be in the effects of temperature stratification. We have retained the compression work term so that numerical convection models can be developed for examining the validity of the solution. Also, it can be predicted that the compression, being the additive one to thermal stratification, also plays an important role in the solution of temperature and vertical velocity as well. A similar variable ξ
=
d
C
∞
′
(
z
)
/
d
z is termed as mass stratification. We shall also study the simultaneous effect of mass stratification in addition to thermal stratification.Introducing the nondimensional quantities(3)
R
=
r
r
0
,
U
=
u
r
0
ν
,
t
=
t
′
ν
r
0
2
,
θ
=
T
′
-
T
∞
′
z
T
0
-
T
∞
′
z
,
ϕ
=
C
′
-
C
∞
′
z
C
0
-
C
∞
′
z
,
P
r
=
ν
α
,
S
c
=
ν
D
G
r
=
g
β
r
0
3
T
0
-
T
∞
′
z
ν
2
,
G
c
=
g
β
∗
r
0
3
C
0
-
C
∞
′
z
ν
2
,
S
=
γ
r
0
T
0
-
T
∞
′
z
,
K
=
ξ
r
0
C
0
-
C
∞
′
zthe governing equations (1) reduce to(4)
∂
U
∂
t
=
1
R
∂
∂
R
R
∂
U
∂
R
+
G
r
θ
+
G
c
ϕ
,
∂
θ
∂
t
=
1
P
r
1
R
∂
∂
R
R
∂
θ
∂
R
-
S
U
,
∂
ϕ
∂
t
=
1
S
c
1
R
∂
∂
R
R
∂
ϕ
∂
R
-
K
Uand the corresponding initial and boundary conditions in nondimensional form are(5)
t
≤
0
:
U
=
0
,
θ
=
0
,
ϕ
=
0
∀
R
t
>
0
:
U
=
0
,
θ
=
1
,
ϕ
=
1
a
t
R
=
1
U
⟶
0
,
θ
⟶
0
,
ϕ
⟶
0
a
s
R
⟶
∞
.


## 3. Solution Technique

To solve the governing nondimensional unsteady equations (4) subject to initial and boundary conditions (5), we apply Laplace transform technique for the case of unit Prandtl number and Schmidt number, as for arbitrary Prandtl number or Schmidt number, the Laplace transform technique leads to a difficult inverse transformation step (integrand of the Bromwich integral becomes a complicated multivalued function). Laplace transforms of (4) with initial conditions of (5) with P
r
=
S
c
=
1 become(6)
d
2
U
¯
d
R
2
+
1
R
d
U
¯
d
R
-
p
U
¯
+
G
r
θ
¯
+
G
c
ϕ
¯
=
0
,
(7)
d
2
θ
¯
d
R
2
+
1
R
d
θ
¯
d
R
-
p
θ
¯
-
S
U
¯
=
0
,
(8)
d
2
ϕ
¯
d
R
2
+
1
R
d
ϕ
¯
d
R
-
p
ϕ
¯
-
K
U
¯
=
0
,where p is the parameter of the Laplace transformation.Solutions of (6), (7), and (8) subject to the boundary conditions (5) are obtained as (the processes of determining U
¯, θ
¯, and ϕ
¯ are presented in Appendix)(9)
U
¯
=
G
r
+
G
c
2
i
L
K
0
R
p
-
i
L
p
K
0
p
-
i
L
-
K
0
R
p
+
i
L
p
K
0
p
+
i
L
,
θ
¯
=
G
c
K
-
S
L
2
K
0
R
p
p
K
0
p
+
S
G
r
+
G
c
2
L
2
K
0
R
p
-
i
L
p
K
0
p
-
i
L
+
K
0
R
p
+
i
L
p
K
0
p
+
i
L
,
ϕ
¯
=
G
r
S
-
K
L
2
K
0
R
p
p
K
0
p
+
K
G
r
+
G
c
2
L
2
K
0
R
p
-
i
L
p
K
0
p
-
i
L
+
K
0
R
p
+
i
L
p
K
0
p
+
i
L
.Inverse Laplace transforms of (9) give the expressions of velocity, temperature, and concentration profile [following Deka and Paul [24], Carslaw and Jaeger [25]] as follows: (10)
U
=
G
r
+
G
c
2
L
i
K
0
R
-
i
L
K
0
-
i
L
-
K
0
R
i
L
K
0
i
L
+
2
G
r
+
G
c
π
L
∫
0
∞
e
-
V
2
t
V
2
sin
⁡
L
t
+
L
cos
⁡
L
t
V
4
+
L
2
Γ
R
,
V
V
d
V
(11)
θ
=
G
c
K
-
S
L
2
1
+
2
π
∫
0
∞
e
-
V
2
t
Γ
R
,
V
d
V
V
+
S
G
r
+
G
c
2
L
2
K
0
R
-
i
L
K
0
-
i
L
+
K
0
R
i
L
K
0
i
L
+
2
S
G
r
+
G
c
π
L
2
∫
0
∞
e
-
V
2
t
V
2
cos
⁡
L
t
-
L
sin
⁡
L
t
V
4
+
L
2
Γ
R
,
V
V
d
V
,
(12)
ϕ
=
G
r
S
-
K
L
2
1
+
2
π
∫
0
∞
e
-
V
2
t
Γ
R
,
V
d
V
V
+
K
G
r
+
G
c
2
L
2
K
0
R
-
i
L
K
0
-
i
L
+
K
0
R
i
L
K
0
i
L
+
2
K
G
r
+
G
c
π
L
2
∫
0
∞
e
-
V
2
t
V
2
cos
⁡
L
t
-
L
sin
⁡
L
t
V
4
+
L
2
Γ
R
,
V
V
d
V
,where L
2
=
S
G
r
+
K
G
c. It is to be noted that i
=
-
1 and in expressions U, θ, and ϕ above, the expressions with complex quantities appear along with their conjugates, thereby resulting in real quantities. This kind of appearance is also present in our forgoing analysis. MATHEMATICA is used to deal with these situations during computations.


### 3.1. Skin Friction

The nondimensional skin frictionτ
=
-
∂
U
/
∂
R
R
=
1 is obtained from the velocity profile (10) as(13)
τ
=
G
r
+
G
c
2
L
i
-
i
L
K
1
-
i
L
K
0
-
i
L
-
i
L
K
1
i
L
K
0
i
L
+
2
G
r
+
G
c
π
L
∫
0
∞
e
-
V
2
t
V
2
sin
⁡
L
t
+
L
cos
⁡
L
t
V
4
+
L
2
Γ
1
V
V
2
d
V
.


### 3.2. Nusselt Number

The nondimensional Nusselt numberN
u
=
-
∂
θ
/
∂
R
R
=
1 is obtained from the temperature profile (11) as(14)
N
u
=
2
G
c
K
-
S
π
L
2
∫
0
∞
e
-
V
2
t
Γ
1
V
d
V
+
S
G
r
+
G
c
2
L
2
i
L
K
1
i
L
K
0
i
L
+
-
i
L
K
1
-
i
L
K
0
-
i
L
+
2
S
G
r
+
G
c
π
L
2
∫
0
∞
e
-
V
2
t
V
2
cos
⁡
L
t
-
L
sin
⁡
L
t
V
4
+
L
2
Γ
1
V
V
2
d
V
.


### 3.3. Sherwood Number

The nondimensional Sherwood numberS
h
=
-
∂
ϕ
/
∂
R
R
=
1 is obtained from the concentration profile (12) as(15)
S
h
=
2
G
r
S
-
K
π
L
2
∫
0
∞
e
-
V
2
t
Γ
1
V
d
V
+
S
G
r
+
G
c
2
L
2
i
L
K
1
i
L
K
0
i
L
+
-
i
L
K
1
-
i
L
K
0
-
i
L
+
2
S
G
r
+
G
c
π
L
2
∫
0
∞
e
-
V
2
t
V
2
cos
⁡
L
t
-
L
sin
⁡
L
t
V
4
+
L
2
Γ
1
V
V
2
d
V
,where (16)
Γ
R
,
V
=
J
0
R
V
Y
0
V
-
Y
0
R
V
J
0
V
J
0
2
V
+
Y
0
2
V
,
Γ
1
V
=
J
1
V
Y
0
V
-
Y
1
V
J
0
V
J
0
2
V
+
Y
0
2
V
.


## 3.1. Skin Friction

The nondimensional skin frictionτ
=
-
∂
U
/
∂
R
R
=
1 is obtained from the velocity profile (10) as(13)
τ
=
G
r
+
G
c
2
L
i
-
i
L
K
1
-
i
L
K
0
-
i
L
-
i
L
K
1
i
L
K
0
i
L
+
2
G
r
+
G
c
π
L
∫
0
∞
e
-
V
2
t
V
2
sin
⁡
L
t
+
L
cos
⁡
L
t
V
4
+
L
2
Γ
1
V
V
2
d
V
.


## 3.2. Nusselt Number

The nondimensional Nusselt numberN
u
=
-
∂
θ
/
∂
R
R
=
1 is obtained from the temperature profile (11) as(14)
N
u
=
2
G
c
K
-
S
π
L
2
∫
0
∞
e
-
V
2
t
Γ
1
V
d
V
+
S
G
r
+
G
c
2
L
2
i
L
K
1
i
L
K
0
i
L
+
-
i
L
K
1
-
i
L
K
0
-
i
L
+
2
S
G
r
+
G
c
π
L
2
∫
0
∞
e
-
V
2
t
V
2
cos
⁡
L
t
-
L
sin
⁡
L
t
V
4
+
L
2
Γ
1
V
V
2
d
V
.


## 3.3. Sherwood Number

The nondimensional Sherwood numberS
h
=
-
∂
ϕ
/
∂
R
R
=
1 is obtained from the concentration profile (12) as(15)
S
h
=
2
G
r
S
-
K
π
L
2
∫
0
∞
e
-
V
2
t
Γ
1
V
d
V
+
S
G
r
+
G
c
2
L
2
i
L
K
1
i
L
K
0
i
L
+
-
i
L
K
1
-
i
L
K
0
-
i
L
+
2
S
G
r
+
G
c
π
L
2
∫
0
∞
e
-
V
2
t
V
2
cos
⁡
L
t
-
L
sin
⁡
L
t
V
4
+
L
2
Γ
1
V
V
2
d
V
,where (16)
Γ
R
,
V
=
J
0
R
V
Y
0
V
-
Y
0
R
V
J
0
V
J
0
2
V
+
Y
0
2
V
,
Γ
1
V
=
J
1
V
Y
0
V
-
Y
1
V
J
0
V
J
0
2
V
+
Y
0
2
V
.


## 4. Steady State Solution

Steady state equations are obtained by neglecting the time derivative terms from (4). Solving these equations and using boundary conditions (5) we obtain the expressions of velocity, temperature, and concentration profiles as(17)
U
S
=
G
r
+
G
c
2
L
i
K
0
R
-
i
L
K
0
-
i
L
-
K
0
R
i
L
K
0
i
L
(18)
θ
S
=
G
c
K
-
S
L
2
+
S
G
r
+
G
c
2
L
2
K
0
R
-
i
L
K
0
-
i
L
+
K
0
R
i
L
K
0
i
L
,
(19)
ϕ
S
=
G
r
S
-
K
L
2
+
K
G
r
+
G
c
2
L
2
K
0
R
-
i
L
K
0
-
i
L
+
K
0
R
i
L
K
0
i
L
.


### 4.1. Skin Friction

Nondimensional skin frictionτ
S
=
-
∂
U
S
/
∂
R
R
=
1 is obtained from velocity profile (17) as(20)
τ
S
=
G
r
+
G
c
2
L
i
-
i
L
K
1
-
i
L
K
0
-
i
L
-
i
L
K
1
R
i
L
K
0
i
L
.


### 4.2. Nusselt Number

Nondimensional Nusselt numberN
u
S
=
-
∂
θ
S
/
∂
R
R
=
1 is obtained from temperature profile (18) as(21)
N
u
S
=
S
G
r
+
G
c
2
L
2
i
L
K
1
i
L
K
0
i
L
+
-
i
L
K
1
-
i
L
K
0
-
i
L
.


### 4.3. Sherwood Number

Nondimensional Sherwood numberS
h
S
=
-
∂
ϕ
S
/
∂
R
R
=
1 is obtained from concentration profile (19) as(22)
S
h
S
=
K
G
r
+
G
c
2
L
2
i
L
K
1
i
L
K
0
i
L
+
-
i
L
K
1
-
i
L
K
0
-
i
L
.It is to be noted that ast
→
∞ the expressions for unsteady velocity, temperature, concentration, skin friction, Nusselt number, and Sherwood number given by (10), (11), (12), (13), (14), and (15), respectively, approach to the corresponding expressions for steady state given by (17), (18), (19), (20), (21), and (22).


## 4.1. Skin Friction

Nondimensional skin frictionτ
S
=
-
∂
U
S
/
∂
R
R
=
1 is obtained from velocity profile (17) as(20)
τ
S
=
G
r
+
G
c
2
L
i
-
i
L
K
1
-
i
L
K
0
-
i
L
-
i
L
K
1
R
i
L
K
0
i
L
.


## 4.2. Nusselt Number

Nondimensional Nusselt numberN
u
S
=
-
∂
θ
S
/
∂
R
R
=
1 is obtained from temperature profile (18) as(21)
N
u
S
=
S
G
r
+
G
c
2
L
2
i
L
K
1
i
L
K
0
i
L
+
-
i
L
K
1
-
i
L
K
0
-
i
L
.


## 4.3. Sherwood Number

Nondimensional Sherwood numberS
h
S
=
-
∂
ϕ
S
/
∂
R
R
=
1 is obtained from concentration profile (19) as(22)
S
h
S
=
K
G
r
+
G
c
2
L
2
i
L
K
1
i
L
K
0
i
L
+
-
i
L
K
1
-
i
L
K
0
-
i
L
.It is to be noted that ast
→
∞ the expressions for unsteady velocity, temperature, concentration, skin friction, Nusselt number, and Sherwood number given by (10), (11), (12), (13), (14), and (15), respectively, approach to the corresponding expressions for steady state given by (17), (18), (19), (20), (21), and (22).


## 5. Boundary Layer Thickness

The peak vertical velocity in the steady state occurs at a nondimensional distanceR
δ, from the surface of the cylinder where R
δ can be obtained from d
U
S
/
d
R
R
=
R
δ
=
0
. This leads to the determination of R
δ from the following expression:(23)
1
+
i
K
0
i
L
K
1
R
δ
-
i
L
+
1
-
i
K
0
-
i
L
K
1
R
δ
+
i
L
=
0
.A list of calculated values of R
δ for different values S, K, Gr, and Gc is tabulated in Table 1.Table 1
Boundary layer thickness for variousS, K, Gr, and Gc.
Gr
Gc
S
K
R
δ
5
5
0.2
0
1.890311
5
5
0
0.4
1.764592
5
5
0.2
0.4
1.699095
5
5
0.2
0.6
1.655900
5
5
0.6
0.4
1.624154
8
5
0.2
0.4
1.671416
10
5
0.2
0.4
1.655901
5
8
0.2
0.4
1.648831
5
10
0.2
0.4
1.624154


## 6. Results and Discussion

In order to understand the physical insight of the problem, numerical computations for velocity, temperature, concentration, skin friction, Nusselt number, and Sherwood number are made for various physical parameters and presented in figures. Solutions of unsteady state for larger time are compared with the solutions of steady state.Velocity profiles represented by Figure2 show the effects of S and K for Gr = Gc = 5 and t
=
1.5 and Figure 3 shows the effects of Gr and Gc for S
=
0.4, K
=
0.2, and t
=
1.5. It is observed from these figures that velocity decreases with the increase in S or increase in K but increases with increase in Gr or Gc. As mentioned by Shapiro and Fedorovich [15], the fluid leads to a stably stratified flow by decreasing the fluid velocity with thermal stratification S
>
0. Remarkably, in our analysis, we have seen that, for mass stratification K
>
0, the fluid velocity still reduces. Therefore, we conclude that the inclusion of mass stratification also leads to a stably stratified flow. Figures 4 and 5, respectively, show the effects of S, K and Gr, Gc on velocity profiles against time at R
=
1.8. The time required to reach the steady state increases with decrease in stratification parameters S or K and increases with increase in the values of Gr or Gc.Figure 2
Effects of thermal and mass stratification on velocity profile at Gr = 5 and Gc = 5.Figure 3
Effects of Gr and Gc on velocity profile atS
=
0.4 and K
=
0.2.Figure 4
Effects ofS and K on velocity profile against time at Gr = 5 and Gc = 5.Figure 5
Effects of Gr and Gc on velocity profile against time atS
=
0.4 and K
=
0.2.Figure6 shows the effects of thermal and mass stratification on temperature profiles at Gr = Gc = 5 at t
=
1.2. It is observed from the figure that temperature decreases with the increase in the stratification parameter S or K. Also, temperature becomes negative in presence of dense stratified fluid. This type of behaviour has also been observed by earlier investigators (Kulkarni et al. [26], Loganathan and Ganesan [14], and Deka and Paul [24]). This is because the fluid with dense stratification near the cylinder can have temperature and concentration lower than the ambient. Figure 7 represents the effects of thermal Grashof number and mass Grashof number on temperature profiles at S
=
0.4, K
=
0.2 for t
=
1.2 and it shows that temperature decreases with increase in Gr or Gc. Figures 8 and 9, respectively, show the effects of S, K and Gr, Gc on temperature profile against time at R
=
1.8. From these figures, it is observed that initially temperature increases sharply with time but for larger time it becomes steady. Also, time required to reach the steady state increases with decrease in S, K, Gr, or Gc. It is also observed that, in absence of thermal stratification parameter, effect of mass stratification on temperature profile is negligible.Figure 6
Effects ofS and K on temperature profile at Gr = 5 and Gc = 5.Figure 7
Effects of Gr and Gc on temperature profile atS
=
0.4 and K
=
0.2.Figure 8
Effects ofS and K on temperature profile against time at Gr = Gc = 5 for R
=
1.8.Figure 9
Effects of Gr and Gc on temperature profile against time at Gr = 5, Gc = 5.Effects of thermal and mass stratification on concentration profiles are shown in Figure10 at Gr = Gc = 5 for t
=
1.5. It is observed from the figure that concentration decreases with an increase in stratification parameter S or K. Also concentration becomes negative in presence of dense stratified fluid. Figure 11 shows the effects of thermal Grashof number and mass Grashof number on concentration profiles at S
=
0.4, K
=
0.2 for t
=
1.5 and that concentration decreases with increase in Gr or Gc. Figures 12 and 13, respectively, show the effects of S, K and Gr, Gc on concentration profiles against time at R
=
1.8. It is found from these figures that initially concentration increases sharply with time but for larger time it becomes steady. Also, time required to reach the steady state increases with decrease in S, K, Gr, or Gc. It is also observed that, in absence of mass stratification parameter, effect of thermal stratification on concentration profile is negligible.Figure 10
Effects ofS and K on concentration profile at Gr = 5, Gc = 5.Figure 11
Effects of Gr and Gc on concentration profile atS
=
0.4, K
=
0.2.Figure 12
Effects ofS and K on concentration profile against time at Gr = 5, Gc = 5.Figure 13
Effects of Gr and Gc on concentration profile against time atS
=
0.4, K
=
0.2.Skin friction represented by Figures14 and 15 shows the effects of S, K and Gr, Gc, respectively. It is observed from these figures that skin friction initially decreases with time but becomes steady at large time. Skin friction increases with an increase in thermal or mass stratification but decreases with an increase in Gr or Gc.Figure 14
Effects of thermal and mass stratification on skin friction at Gr = 5, Gc = 5.Figure 15
Effects of Gr and Gc on skin friction forS
=
0.4 and K
=
0.2.Figures16 and 17, respectively, depict the effects of thermal stratification and mass stratification and thermal Grashof number and mass Grashof number on Nusselt number against time. It is found that initially Nusselt number decreases very sharply but after certain time it becomes steady. Nusselt number increases with an increase in thermal or mass stratification but decreases with an increase in thermal Grashof number or mass Grashof number. Also effect of mass stratification on Nusselt number is negligible, when thermal stratification is absent.Figure 16
Effects ofS and K on Nusselt number for Gr = 5 and Gc = 5.Figure 17
Effects of Gr and Gc on Nusselt number forS
=
0.4 and K
=
0.2.Figures18 and 19, respectively, show the effect of thermal stratification and mass stratification and thermal Grashof number and mass Grashof number on Sherwood number, that is, rate of mass transfer. It is observed that initially Sherwood number decreases sharply but becomes steady after certain time. Sherwood number increases with an increase in thermal stratification or mass stratification but decreases with an increase in thermal Grashof number or mass Grashof number. Also, effect of thermal stratification on rate of mass transfer is negligible when the fluid is free from mass stratification. A decrease in the thermal stratification parameter or the concentration stratification parameter leads to a decrease in the Nusselt number and Sherwood number and this trend becomes more pronounced as the fluid moves downstream (Figure 2). This is due to the fact that increasing the thermal and concentration stratification parameters decreases the buoyancy force due to thermal and solutal gradients and retards the flow, increasing the thermal and concentration boundary layer thickness and thus decreasing the heat and mass transfer rates between the fluid and the cylinder wall.Figure 18
Effects ofS and K on Sherwood number for Gr = 5 and Gc = 5.Figure 19
Effects of Gr and Gc on Sherwood number forS
=
0.4 and K
=
0.2.


## 7. Conclusions

(i)
Solution expressions obtained for unsteady state approaches to the solutions of steady state ast becomes large.
(ii)
The time required to reach steady state velocity increases with increase in Gr or Gc but decreases with increase inS or K, while time required to reach the steady state temperature and concentration increases with decrease in S, K, Gr, or Gc.
(iii)
Velocity increases with increase in Gr or Gc but decreases with increase inS or K.
(iv)
Temperature decreases with increase inS, K, Gr, or Gc. In presence of high stratification temperature becomes negative.
(v)
Concentration decreases with increase inS, K, Gr, or Gc. In presence of high stratification concentration becomes negative.
(vi)
Skin friction decreases with increase in Gr or Gc but increases with increase inS or K.
(vii)
Rate of heat transfer increases as Gr,K, S, or K increases.
(viii)
Rate of mass transfer increases with increase in Gr,K, S, or K.



---

*Source: 8410691-2017-01-11.xml*