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.
University Grants Commission1. 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 radius r0 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 (T0) and concentration (C0) 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 dT∞′/dz and dC∞′/dz 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∂∂rr∂u∂r+gβT′-T∞′+gβ∗C′-C∞′,∂T′∂t′=αr∂∂rr∂T′∂r-γu,∂C′∂t′=Dr∂∂rr∂C′∂r-ξuwith initial and boundary conditions as(2)t′≤0:u=0,T′=T∞′C′=C∞′∀rt′>0:u=0,T′=T0,C′=C0atr=r0u⟶0,T′⟶T∞′,C′⟶T∞′asr⟶∞.
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 [≡dT∞′(z)/dz for Boussinesq flow of liquids or gases, ≡dT∞′(z)/dz+g/Cp for a perfect gas with pressure work term retained], and z is height. Thus γ is the combination of thermal stratification and compression, where dT∞′(z)/dz stands for thermal stratification and g/Cp for compression. In an adiabatic environment, dT∞′z/dz+g/Cp=0. In that case dT∞′(z)/dz=-g/Cp and (-g/Cp) 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}, Cp = 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 ξ=dC∞′(z)/dz 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=rr0,U=ur0ν,t=t′νr02,θ=T′-T∞′zT0-T∞′z,ϕ=C′-C∞′zC0-C∞′z,Pr=να,Sc=νDGr=gβr03T0-T∞′zν2,Gc=gβ∗r03C0-C∞′zν2,S=γr0T0-T∞′z,K=ξr0C0-C∞′zthe governing equations (1) reduce to(4)∂U∂t=1R∂∂RR∂U∂R+Grθ+Gcϕ,∂θ∂t=1Pr1R∂∂RR∂θ∂R-SU,∂ϕ∂t=1Sc1R∂∂RR∂ϕ∂R-KUand the corresponding initial and boundary conditions in nondimensional form are(5)t≤0:U=0,θ=0,ϕ=0∀Rt>0:U=0,θ=1,ϕ=1atR=1U⟶0,θ⟶0,ϕ⟶0asR⟶∞.
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 Pr=Sc=1 become(6)d2U¯dR2+1RdU¯dR-pU¯+Grθ¯+Gcϕ¯=0,(7)d2θ¯dR2+1Rdθ¯dR-pθ¯-SU¯=0,(8)d2ϕ¯dR2+1Rdϕ¯dR-pϕ¯-KU¯=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¯=Gr+Gc2iLK0Rp-iLpK0p-iL-K0Rp+iLpK0p+iL,θ¯=GcK-SL2K0RppK0p+SGr+Gc2L2K0Rp-iLpK0p-iL+K0Rp+iLpK0p+iL,ϕ¯=GrS-KL2K0RppK0p+KGr+Gc2L2K0Rp-iLpK0p-iL+K0Rp+iLpK0p+iL.
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=Gr+Gc2LiK0R-iLK0-iL-K0RiLK0iL+2Gr+GcπL∫0∞e-V2tV2sinLt+LcosLtV4+L2ΓR,VVdV(11)θ=GcK-SL21+2π∫0∞e-V2tΓR,VdVV+SGr+Gc2L2K0R-iLK0-iL+K0RiLK0iL+2SGr+GcπL2∫0∞e-V2tV2cosLt-LsinLtV4+L2ΓR,VVdV,(12)ϕ=GrS-KL21+2π∫0∞e-V2tΓR,VdVV+KGr+Gc2L2K0R-iLK0-iL+K0RiLK0iL+2KGr+GcπL2∫0∞e-V2tV2cosLt-LsinLtV4+L2ΓR,VVdV,where L2=SGr+KGc. 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/∂RR=1 is obtained from the velocity profile (10) as(13)τ=Gr+Gc2Li-iLK1-iLK0-iL-iLK1iLK0iL+2Gr+GcπL∫0∞e-V2tV2sinLt+LcosLtV4+L2Γ1VV2dV.
3.2. Nusselt Number
The nondimensional Nusselt number Nu=-∂θ/∂RR=1 is obtained from the temperature profile (11) as(14)Nu=2GcK-SπL2∫0∞e-V2tΓ1VdV+SGr+Gc2L2iLK1iLK0iL+-iLK1-iLK0-iL+2SGr+GcπL2∫0∞e-V2tV2cosLt-LsinLtV4+L2Γ1VV2dV.
3.3. Sherwood Number
The nondimensional Sherwood number Sh=-∂ϕ/∂RR=1 is obtained from the concentration profile (12) as(15)Sh=2GrS-KπL2∫0∞e-V2tΓ1VdV+SGr+Gc2L2iLK1iLK0iL+-iLK1-iLK0-iL+2SGr+GcπL2∫0∞e-V2tV2cosLt-LsinLtV4+L2Γ1VV2dV,where (16)ΓR,V=J0RVY0V-Y0RVJ0VJ02V+Y02V,Γ1V=J1VY0V-Y1VJ0VJ02V+Y02V.
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)US=Gr+Gc2LiK0R-iLK0-iL-K0RiLK0iL(18)θS=GcK-SL2+SGr+Gc2L2K0R-iLK0-iL+K0RiLK0iL,(19)ϕS=GrS-KL2+KGr+Gc2L2K0R-iLK0-iL+K0RiLK0iL.
4.1. Skin Friction
Nondimensional skin friction τS=-∂US/∂RR=1 is obtained from velocity profile (17) as(20)τS=Gr+Gc2Li-iLK1-iLK0-iL-iLK1RiLK0iL.
4.2. Nusselt Number
Nondimensional Nusselt number NuS=-∂θS/∂RR=1 is obtained from temperature profile (18) as(21)NuS=SGr+Gc2L2iLK1iLK0iL+-iLK1-iLK0-iL.
4.3. Sherwood Number
Nondimensional Sherwood number ShS=-∂ϕS/∂RR=1 is obtained from concentration profile (19) as(22)ShS=KGr+Gc2L2iLK1iLK0iL+-iLK1-iLK0-iL.
It is to be noted that as t→∞ 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 distance Rδ, from the surface of the cylinder where Rδ can be obtained from dUS/dRR=Rδ=0. This leads to the determination of Rδ from the following expression:(23)1+iK0iLK1Rδ-iL+1-iK0-iLK1Rδ+iL=0.A list of calculated values of Rδ for different values S, K, Gr, and Gc is tabulated in Table 1.
Boundary layer thickness for various S, 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 Figure 2 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.
Effects of thermal and mass stratification on velocity profile at Gr = 5 and Gc = 5.
Effects of Gr and Gc on velocity profile at S=0.4 and K=0.2.
Effects of S and K on velocity profile against time at Gr = 5 and Gc = 5.
Effects of Gr and Gc on velocity profile against time at S=0.4 and K=0.2.
Figure 6 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.
Effects of S and K on temperature profile at Gr = 5 and Gc = 5.
Effects of Gr and Gc on temperature profile at S=0.4 and K=0.2.
Effects of S and K on temperature profile against time at Gr = Gc = 5 for R=1.8.
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 Figure 10 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.
Effects of S and K on concentration profile at Gr = 5, Gc = 5.
Effects of Gr and Gc on concentration profile at S=0.4, K=0.2.
Effects of S and K on concentration profile against time at Gr = 5, Gc = 5.
Effects of Gr and Gc on concentration profile against time at S=0.4, K=0.2.
Skin friction represented by Figures 14 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.
Effects of thermal and mass stratification on skin friction at Gr = 5, Gc = 5.
Effects of Gr and Gc on skin friction for S=0.4 and K=0.2.
Figures 16 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.
Effects of S and K on Nusselt number for Gr = 5 and Gc = 5.
Effects of Gr and Gc on Nusselt number for S=0.4 and K=0.2.
Figures 18 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.
Effects of S and K on Sherwood number for Gr = 5 and Gc = 5.
Effects of Gr and Gc on Sherwood number for S=0.4 and K=0.2.
7. Conclusions
Solution expressions obtained for unsteady state approaches to the solutions of steady state as t becomes large.
The time required to reach steady state velocity increases with increase in Gr or Gc but decreases with increase in S or K, while time required to reach the steady state temperature and concentration increases with decrease in S, K, Gr, or Gc.
Velocity increases with increase in Gr or Gc but decreases with increase in S or K.
Temperature decreases with increase in S, K, Gr, or Gc. In presence of high stratification temperature becomes negative.
Concentration decreases with increase in S, K, Gr, or Gc. In presence of high stratification concentration becomes negative.
Skin friction decreases with increase in Gr or Gc but increases with increase in S or K.
Rate of heat transfer increases as Gr, K, S, or K increases.
Rate of mass transfer increases with increase in Gr, K, S, or K.
Appendix
Laplace transform of (4) with initial conditions in (5) and with Pr=Sc=1 is (6), (7), and (8) with transformed boundary conditions(A.1)U¯=0,θ¯=1pϕ¯=1patR=1,U¯⟶0,θ¯⟶0ϕ¯⟶0asR⟶∞,where U¯p,R, T¯p,R, and C¯p,R are Laplace transforms of Ut,R, Tt,R, and Ct,R, respectively, and p is Laplace transform parameter.
Upon applying (7) and (8), we eliminate θ¯ and ϕ¯ from (6) that give rise to (A.2)d2dR2+1RddR-a2d2dR2+1RddR-b2U¯=0,where a2=p+iL, b2=p-iL, L2=SGr+KSc.
Now to determine U¯ from (A.2), we rewrite (A.2) as (A.3)d2dR2+1RddR-a2W=0with (A.4)d2dR2+1RddR-b2U¯=W.
Equation (A.3) is the modified Bessel equation of order zero and its solution is (A.5)W=C1I0aR+C2K0aR,where C1 and C2 are arbitrary constants. In view of (A.5) and using the following identities (see Carslaw and Jaeger [25]) (A.6)∫RI0aRK0bRdR=Ra2-b2aI1aRK0bR+bI0aRK1bR,∫RK0aRK0bRdR=Ra2-b2bK0aRK1bR-aK0bRK1aR,∫RI0aRI0bRdR=Ra2-b2aI0bRI1aR-bI0aRI1bR,I0axK1ax+I1axK0ax=1axthe expression for U¯ is obtained after applying variation of parameter technique as(A.7)U¯=C1I0aR2iL+C2K0aR2iL+C3I0bR+C4K0bR,where C3 and C4 are arbitrary constants.
Now, since I0aR and I0bR are unbounded as R→∞, we set C1=C1=0 and so (A.7) is then read as(A.8)U¯=C2K0aR2iL+C4K0bR.Now, using (A.8) we have from (6) (A.9)C2K0aR-2iLC4K0bR+2Grθ¯+2Gcϕ¯=0.Using boundary conditions (A.1) in (A.8) and (A.9) and solving we have(A.10)C2=-Gr+GcpK0p+iL,C4=Gr+Gc2piLK0p-iL.Thus we have(A.11)U¯=Gr+Gc2iLK0Rp-iLpK0p-iL-K0Rp+iLpK0p+iL.Similarly, substituting the expression of U¯ in (7) and (8) and solving, the expressions of θ¯ and ϕ¯ can be obtained as (A.12)θ¯=GcK-SL2K0RppK0p+SGr+Gc2L2K0Rp-iLpK0p-iL+K0Rp+iLpK0p+iL,ϕ¯=GrS-KL2K0RppK0p+KGr+Gc2L2K0Rp-iLpK0p-iL+K0Rp+iLpK0p+iL.
NomenclatureC:
Concentration
C0:
Concentration at the surface of the cylinder
D:
Mass diffusion coefficient
Gr:
Thermal Grashof number
Gc:
Mass Grashof number
g:
Acceleration due to gravity
J0:
Bessel function of first kind and order zero
J1:
Bessel function of first kind and order one
K:
Dimensionless mass stratification parameter
K0:
Modified Bessel function of second kind and order zero
K1:
Modified Bessel function of second kind and order one
Nu:
Nusselt number
Pr:
Prandtl number
r:
Radial coordinate measured from the axis of the cylinder
r0:
Radius of the cylinder
R:
Dimensionless radial coordinate
S:
Dimensionless thermal stratification parameter
Sc:
Schmidt number
Sh:
Sherwood number
t′:
Time
t:
Dimensionless time
T:
Temperature
T0:
Temperature at the surface of the cylinder
u:
z-component of velocity
U:
Dimensionless velocity
Y0:
Bessel function of second kind and order zero
Y1:
Bessel function of second kind and order one
α:
Thermal diffusivity
ν:
Kinematic viscosity
β:
Volumetric coefficient of thermal expansion
β∗:
Volumetric coefficient of expansion with concentration
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
One of the authors acknowledges the support from the University Grants Commission, New Delhi, India.
GebhartB.PeraL.The nature of vertical natural convection flows resulting from the combined buoyancy effects of thermal and mass diffusionBottemanneF. A.Experimental results of pure and simultaneous heat and mass transfer by free convection about a vertical cylinder for Pr=0.71 and Sc=0.63ChenT. S.YuhC. F.Combined heat and mass transfer in natural convection along a vertical cylinderVelusamyK.GargV. K.Transient natural convection over a heat generating vertical cylinderGanesanP.RaniH. P.Transient natural convection along vertical cylinder with heat and mass transferGanesanP.LoganathanP.Effects of mass transfer and flow past a moving vertical cylinder with constant heat fluxGanesanP.LoganathanP.Unsteady natural convective flow past a moving vertical cylinder with heat and mass transferGanesanP.LoganathanP.Heat and mass flux effects on a moving vertical cylinder with chemically reactive species diffusionRaniH. P.Transient natural convection along a vertical cylinder with variable surface temperature and mass diffusionDevakarM.SreenivasuD.ShankarB.Analytical solutions of couple stress fluid flows with slip boundary conditionsDekaR. K.PaulA.ChalihaA.Transient free convection flow past vertical cylinder with constant heat flux and mass transferTakharH. S.ChamkhaA. J.NathG.Natural convection flow from a continuously moving vertical surface immersed in a thermally stratified mediumTakharH. S.ChamkhaA. J.NathG.Natural convection on a vertical cylinder embedded in a thermally stratified high-porosity mediumLoganathanP.GanesanP.Numerical study of double-diffusive, free convective flow past a moving vertical cylinderShapiroA.FedorovichE.Unsteady convectively driven flow along a vertical plate immersed in a stably stratified fluidShapiroA.FedorovichE.Prandtl number dependence of unsteady natural convection along a vertical plate in a stably stratified fluidChengC.-Y.Combined heat and mass transfer in natural convection flow from a vertical wavy surface in a power-law fluid saturated porous medium with thermal and mass stratificationSrinivasacharyaD.ReddyG. S.Effect of double stratification on free convection in a power-law fluid saturated porous mediumSrinivasacharyaD.RamReddyC.Natural convection heat and mass transfer in a micropolar fluid with thermal and mass stratificationRathish KumarB. V.Krishna MurthyS. V. S. S. N. V. G.Double diffusive free convection induced by verticalwavy surface in a doubly stratified Darcy porous medium under the influence of Soret and Dufour effectNeaguM.Free convective heat and mass transfer induced by a constant heat and mass fluxes vertical wavy wall in a non-Darcy double stratified porous mediumDekaR. K.PaulA.Transient free convection flow past an infinite moving vertical cylinder in a stably stratified fluidDekaR. K.PaulA.Transient free convection flow past an infinite vertical cylinder with thermal stratificationDekaR. K.PaulA.Convectively driven flow past an infinite moving vertical cylinder with thermal and mass stratificationCarslawH. S.JaegerJ. C.KulkarniA. K.JacobsH. R.HwangJ. J.Similarity solution for natural convection flow over an isothermal vertical wall immersed in thermally stratified medium