Chemical Reaction Effect on Transient Free Convective Flow past an Infinite Moving Vertical Cylinder

Abstract

An analysis is performed to study the heat and mass transfer on the flow past an infinite moving vertical cylinder, in the presence of first-order chemical reaction. The closed-form solutions of the dimensionless governing partial differential equations are obtained in terms of Bessel's functions and modified Bessel's functions by the Laplace transform technique. The transient velocity profiles, temperature profiles, and concentration profiles are studied for various sets of physical parameters, namely, the chemical reaction parameter, Prandtl number, Schmidt number, thermal Grashof number, mass Grashof number, and time. The skin friction, Nusselt number, and Sherwood number are also obtained and presented in graphs. It is observed that in presence of as well as increase in chemical reaction the flow velocity decreases. Also, in presence of destructive chemical reaction the concentration profile and Sherwood number tend to the steady state at large time.

1. Introduction

Combined heat and mass transfer in natural convection flows along moving vertical cylinders has got considerable attention in the last few decades because of its wide engineering and industrial application such as in hot rolling, hot extrusion, nuclear reactor cooling system, and underground energy system. Sparrow and Gregg [1] first studied the heat transfer from vertical cylinders. Goldstein and Briggs [2] presented an analysis of the transient free convective flow past vertical flat plate and circular cylinder for the unit and variable Prandtl number by employing the Laplace transform technique. Ishak [3] studied mixed convection over a vertical cylinder in presence of heat flux, while Lien et al. [4] studied the effects of free convection and mass transfer on the flow past an impulsively moving infinite vertical circular cylinder.

But, in nature, the presence of pure air or water is rather impossible. It is always possible that either foreign mass is present naturally in air or water, or foreign masses are mixed with air or water. A simple example is the naturally available water vapor that causes the flow of air. In many chemical engineering processes, there is chemical reaction between a foreign mass and the fluid in which the cylinder is moving. The rate of reaction, which is directly proportional to the concentration, is termed as first order chemical reaction. The effect of a chemical reaction depends on whether the reaction is homogeneous or heterogeneous. Chambré and Young [5] have analyzed a first order chemical reaction in the neighborhood of a horizontal plate. Das et al. [6] presented an analytical investigation of the effect of a homogeneous first order chemical reaction on the flow past an impulsively started infinite vertical plate with uniform heat flux and mass transfer by the usual Laplace transform technique. Muthucumaraswamy and Ganesan [7] studied the effect of first-order chemical reaction on flow past an impulsively started vertical plate with uniform heat and mass flux by employing an implicit finite-difference scheme of the Crank-Nicolson type. Gangadhar et al. [8] studied hydromagnetic mixed convection heat and mass transfer over a vertical plate in the presence of magnetic field and chemical reaction. They have transformed the governing partial differential equations to a system of nonlinear differential equations and obtained the solution by using the Runge-Kutta-fourth order method with a shooting technique.

Ganesan and Rani [9] presented a numerical solution of transient convection flow over a vertical cylinder under the combined buoyancy effects of heat and mass transfer in presence of first order chemical reaction. Again, Ganesan and Loganathan [10] investigated the combined effects of heat and mass transfer on a moving vertical cylinder under the assumption that chemically reactive species are emitted from the surface of the cylinder and are diffused into the fluid. Recently, Rani et al. [11] and Machireddy [12] studied the effect of chemical reaction on flow past over a vertical cylinder numerically under different physical situations.

However, closed-form solutions on natural convection along a moving vertical cylinder are very limited in the literature, which are very much necessary for validating numerical models. This may be due to the complicated mathematical approach involved in the analyses of such models. In the present paper, we have studied the combined buoyancy effects of heat and mass transfer from an infinite moving vertical cylinder in the presence of first order chemical reaction. The closed-form solutions of the system of nondimensional partial differential equations are obtained in terms of Bessel’s functions and modified Bessel’s functions by the Laplace transform technique. The effects of various physical parameters in the flow problem are discussed graphically.

2. Mathematical Analysis

Consider the unsteady, laminar, and incompressible viscous flow over an infinite vertical cylinder of radius 𝑟₀. Initially, it is assumed that the cylinder and fluid are at the same temperature 𝑇^′_∞ and concentration 𝐶^′_∞. At time 𝑡^′≥0, the cylinder starts to move in the vertical direction with constant velocity 𝑢₀, and the temperature and concentration near the cylinder are raised to constant temperature 𝑇^′_𝑤 and concentration 𝐶^′_𝑤, respectively. Here, the 𝑥-axis is taken vertically upward along the axis of the cylinder, and the radial coordinate 𝑟 is taken normal to the cylinder. The viscous dissipation is assumed to be negligible in the energy equation, as the motion is due to free convection only. It is also assumed that there exists a homogeneous first-order chemical reaction between the fluid and species concentration. Also, it is assumed that the level of species concentration is very low, and, hence, heat generated during chemical reaction can be neglected. Under these assumptions, the governing boundary-layer equations for momentum, energy, and concentration for free convective flow with Boussinesq’s approximation are as follows: 𝜕𝑢𝜕𝑡^′=𝜈𝑟𝜕𝑟𝜕𝑟𝜕𝑢𝑇𝜕𝑟+𝑔𝛽^′−𝑇^′_∞+𝑔𝛽^∗𝐶^′−𝐶^′_∞,𝜕𝑇^′𝜕𝑡^′=𝛼𝑟𝜕𝑟𝜕𝑟𝜕𝑇^′,𝜕𝑟𝜕𝐶^′𝜕𝑡^′=𝐷𝑟𝜕𝑟𝜕𝑟𝜕𝐶^′𝜕𝑟−𝑘₁𝐶^′, (1) with the following initial and boundary conditions: 𝑡^′≤0∶𝑢=0,𝑇^′=𝑇^′_∞,𝐶^′=𝐶^′_∞𝑡∀𝑟,_′>0∶𝑢=𝑢₀,𝑇^′=𝑇^′_𝑤,𝐶^′=𝐶^′_𝑤at𝑟=𝑟₀,𝑢⟶0,𝑇^′⟶𝑇^′_∞,𝐶^′⟶𝐶^′_∞as𝑟⟶∞. (2) All the physical variables are mentioned in the Nomenclature section. In order to transform the governing equations into nondimensional form, we introduce the following dimensionless quantities: 𝑟𝑅=𝑟₀𝑢,𝑈=𝑢₀𝑇,𝑇=^′−𝑇^′_∞𝑇^′_𝑤−𝑇^′_∞,𝐶𝐶=^′−𝐶^′_∞𝐶^′_𝑤−𝐶^′_∞,𝜈Pr=𝛼𝜈,Sc=𝐷𝑡,𝑡=^′𝜈𝑟²₀,Gr=𝑔𝛽𝑟²₀𝑇^′_𝑤−𝑇^′_∞𝑢₀𝜈,Gc=𝑔𝛽^∗𝑟²₀𝐶^′_𝑤−𝐶^′_∞𝑢₀𝜈𝑘,𝐾=¹𝑟²₀𝜈, (3) which reduces the governing equations (1) into 𝜕𝑈=𝜕𝜕𝑡²𝑈𝜕𝑅²+1𝑅𝜕𝑈𝜕𝑅+Gr𝑇+Gc𝐶,𝜕𝑇=1𝜕𝑡𝜕Pr²𝑇𝜕𝑅²+1𝑅𝜕𝑇,𝜕𝑅𝜕𝐶=1𝜕𝑡𝜕Sc²𝐶𝜕𝑅²+1𝑅𝜕𝐶𝜕𝑅−𝐾𝐶, (4) with the following initial and boundary conditions: 𝑡≤0∶𝑈=0,𝑇=0,𝐶=0∀𝑅,𝑡>0∶𝑈=1,𝑇=1,𝐶=1at𝑅=1,𝑈⟶0,𝑇⟶0,𝐶⟶0as𝑅⟶∞. (5)

3. Solution Technique

To solve the governing nondimensional unsteady equations (4) subject to initial and boundary conditions (5), we apply the Laplace transform technique.

The Laplace transformation of (4) subject to initial conditions in (5) gives the following: 𝑑²𝑈𝑑𝑅²+1𝑅𝑑𝑈𝑑𝑅−𝑝𝑈+Gr𝑇+Gc𝑑𝐶=0,₂𝑇𝑑𝑅²+1𝑅𝑑𝑇𝑑𝑅−𝑝Pr𝑑𝑇=0,₂𝐶𝑑𝑅²+1𝑅𝑑𝐶𝑑𝑅−Sc𝑝+𝐾𝐶=0, (6)(7)(8) where 𝑝 is the parameter of the Laplace transformation defined by 𝐿{𝑓(𝑡)}=𝐹(𝑝), with 𝐿 being the Laplace operator, and 𝑈(𝑅,𝑝), 𝑇(𝑅,𝑝), and 𝐶(𝑅,𝑝) are the Laplace transforms of 𝑈(𝑅,𝑡), 𝑇(𝑅,𝑡), and 𝐶(𝑅,𝑡) respectively.

Solutions of (7) and (8) subject to the transformed boundary conditions (5) give the following: 𝐾𝑇=⁰𝑅𝑝Pr𝑝𝐾₀,𝑝Pr𝐾𝐶=⁰𝑅Sc𝑝+𝐾𝑝𝐾₀.Sc𝑝+𝐾 (9)

Substituting the expressions of 𝑇 and 𝐶 in (6) and solving, the expression of 𝑈 can be obtained as following: 𝐾𝑈=⁰𝑅√𝑝𝑝𝐾₀√𝑝+Gc𝑝𝑝()×𝐾1−Sc−𝐾₀𝑅Sc𝑝+𝐾𝐾₀−𝐾Sc𝑝+𝐾⁰𝑅√𝑝𝐾₀√𝑝+Gr𝑝²()𝐾1−Pr⁰𝑅𝑝Pr𝐾₀−𝐾𝑝Pr⁰𝑅√𝑝𝐾₀√𝑝. (10)

Using the theorem of the inverse Laplace transform on (9)-(10), it can be readily shown that (cf. Carslaw and Jaeger [13]) 𝐶=𝑒^−𝐾𝑡12𝜋𝑖^𝛾+𝑖∞_{𝛾−𝑖∞}𝑒^𝑝𝑡𝐾₀𝑅Sc𝑝𝑝−𝐾𝐾₀1Sc𝑝𝑑𝑝,𝑈=×2𝜋𝑖_{𝛾+𝑖∞𝛾−𝑖∞}𝑒^𝑝𝑡𝐾₀𝑅√𝑝𝑝𝐾₀√𝑝+Gc𝑝𝑝()×𝐾1−Sc−𝐾₀𝑅Sc𝑝+𝐾𝐾₀−𝐾Sc𝑝+𝐾⁰𝑅√𝑝𝐾₀√𝑝+Gr𝑝²()𝐾1−Pr⁰𝑅𝑝Pr𝐾₀−𝐾𝑝Pr₀𝑅√𝑝𝐾₀√𝑝1𝑑𝑝,𝑇=2𝜋𝑖^𝛾+𝑖∞_{𝛾−𝑖∞}𝑒^𝑝𝑡𝐾₀𝑅𝑝Pr𝑝𝐾₀𝑝Pr𝑑𝑝. (11)(12)(13)

Now, the integrand of (11) has a branch point at 𝑝=0 and a simple pole at 𝑝=𝐾. Also, 𝐾₀(Sc𝑝) do not have zero at any point in the real and imaginary axes, if the branch cut is made along the negative real axis. To obtain 𝐶(𝑡,𝑅) from 𝐶(𝑝,𝑅), we use the adjoining contour shown in Figure 1. Accordingly, the line integral in (11) may be replaced by the limit of the sum of the integrals over FE, ED, DC, CB, and BA as 𝑆₁→∞ and 𝑆₀→0. It can be easily shown that the values of the integral (11) along the paths DC, BA, and FE approache zero as 𝑆₁→∞ and 𝑆₀→0.

Figure 1 

Contour for the inverse integral.

Along the paths CB and ED, we choose 𝑝=𝑉²𝑒^𝑖𝜋/Sc and 𝑝=𝑉²𝑒^−𝑖𝜋/Sc, respectively. Then we consider the following:

The line integral along CB=𝑒^−𝐾𝑡1𝜋𝑖^∞₀𝑒^{−𝑉2𝑡/Sc}𝐽₀(𝑅𝑉)−𝑖𝑌₀(𝑅𝑉)𝑉²+Sc𝐾𝐽₀(𝑉)−𝑖𝑌₀(𝑉)𝑉𝑑𝑉, (14)

and the line integral along ED=𝑒^−𝐾𝑡1𝜋𝑖^∞₀𝑒^{−𝑉2𝑡/Sc}𝐽₀(𝑅𝑉)+𝑖𝑌₀(𝑅𝑉)𝑉²+Sc𝐾𝐽₀(𝑉)+𝑖𝑌₀(𝑉)𝑉𝑑𝑉. (15)

Now, the residue of the integrand of (11) at 𝑝=𝐾 is 𝑒^𝐾𝑡𝐾₀𝑅√𝐾Sc𝐾₀√𝐾Sc. (16)

Thus, using (14), (15), and (16) and applying theory of residues, we have the Laplace inverse of 𝐶, which is the concentration given as 𝐶=𝑒^𝐾𝑡𝐾₀𝑅√𝐾Sc𝐾₀√+𝐾Sc2𝑒^−𝐾𝑡𝜋^∞₀𝑒^{−𝑉2𝑡/Sc}𝑉²Γ()().+Sc𝐾𝑅,𝑉𝑉𝑑𝑉for𝐾≠0 (17) In absence of chemical reaction (𝐾=0), it can be readily shown that 2𝐶=1+𝜋^∞₀𝑒^{−𝑉2𝑡/Sc}Γ(𝑅,𝑉)𝑑𝑉𝑉, (18) where Γ()=𝐽𝑅,𝑉⁰(𝑅𝑉)𝑌₀(𝑉)−𝑌₀(𝑅𝑉)𝐽₀(𝑉)𝐽²₀(𝑉)+𝑌²₀(𝑉). (19)

Similarly, we obtain the expressions for velocity and temperature from (12) and (13) for 𝐾≠0 and 𝐾=0, respectively, as follows: 2𝑈=1+𝜋^∞₀𝑒^−𝑉2𝑡Γ(𝑅,𝑉)𝑑𝑉𝑉+Gc𝐾𝐾Sc1−⁰𝑅√𝐾Sc𝐾₀√+𝐾Sc2PrGr(Pr−1)𝜋^∞₀1−𝑒^{−𝑡𝑉2/Pr}×Γ𝑉𝑅,√(Pr−Γ𝑅,𝑉)𝑑𝑉𝑉³−2Gc𝜋^∞₀𝑉𝑒^{−(𝑉2/Sc+𝐾)𝑡}𝑉²+𝐾Sc𝑉²)−𝑒/Sc+𝐾(1−Sc+𝐾Sc_−𝑉2𝑡𝑉𝑉²()Γ()21−Sc+𝐾Sc𝑅,𝑉𝑑𝑉,𝑈=1+𝜋^∞₀𝑒^−𝑉2𝑡Γ(𝑅,𝑉)𝑑𝑉𝑉+2PrGr(Pr−1)𝜋^∞₀1−𝑒^{−𝑉2𝑡/Pr}×Γ𝑉𝑅,√(Pr−Γ𝑅,𝑉)𝑑𝑉𝑉³+2GcSc(Sc−1)𝜋^∞₀1−𝑒^{−𝑉2𝑡/Sc}×Γ𝑉𝑅,√(Sc−Γ𝑅,𝑉)𝑑𝑉𝑉³,2𝑇=1+𝜋^∞₀𝑒^{−𝑡𝑉2/Pr}Γ(𝑅,𝑉)𝑑𝑉𝑉. (20)

Knowing the velocity, temperature, and concentration field deduced above, we now determine the skin friction (𝜏), Nusselt number (Nu), and Sherwood number (Sh). In nondimensional form, the skin friction and Sherwood number for 𝐾≠0 and 𝐾=0, respectively, are obtained as follows: 𝜏=−𝜕𝑈𝜕𝑅_𝑅=1=2𝜋^∞₀𝑒^−𝑉2𝑡Γ₁(𝑉)𝑑𝑉−Gc√𝐾𝐾Sc¹√𝐾Sc𝐾₀√+𝐾Sc2PrGr()𝜋×Pr−1_∞01−𝑒^{−𝑉2𝑡/Pr}1√ΓPr¹𝑉√Pr−Γ₁(𝑉)𝑑𝑉𝑉²+2Gc𝜋^∞₀𝑉²𝑒^{−𝑉2𝑡/Sc−𝐾𝑡}𝑉²+𝐾Sc𝑉²()−𝑒/Sc+𝐾1−Sc+𝐾Sc_−𝑉2𝑡𝑉²()Γ1−Sc+𝐾Sc¹(𝑉)(),2𝑑𝑉𝐾≠0𝜏=𝜋^∞₀𝑒^−𝑉2𝑡Γ₁(𝑉)+𝑑𝑉2PrGr(Pr−1)𝜋^∞₀1−𝑒^{−𝑉2𝑡/Pr}×1√ΓPr¹𝑉√Pr−Γ₁(𝑉)𝑑𝑉𝑉²+2GcSc()𝜋×Sc−1_∞01−𝑒^{−𝑉2𝑡/Sc}1√ΓSc¹𝑉√Sc−Γ₁(𝑉)𝑑𝑉𝑉²(𝐾=0), (21) where Γ₁(𝑉)=𝐽₁(𝑉)𝑌₀(𝑉)−𝑌₁(𝑉)𝐽₀(𝑉)𝐽²₀(𝑉)+𝑌²₀(𝑉),Nu=−𝜕𝑇𝜕𝑅_𝑅=1=2𝜋^∞₀𝑒^{−𝑉2𝑡/Pr}Γ₁(𝑉)𝑑𝑉,Sh=−𝜕𝐶𝜕𝑅_𝑅=1=√𝐾Sc𝐾₁√𝐾Sc𝐾₀√+𝐾Sc2𝑒^−𝐾𝑡𝜋^∞₀𝑒^{−𝑉2𝑡/Sc}𝑉²Γ+𝐾Sc¹(𝑉)𝑉²(),2𝑑𝑉𝐾≠0Sh=𝜋^∞₀𝑒^{−(𝑉2𝑡)/Sc}Γ₁(𝑉)𝑑𝑉(𝐾=0). (22)

4. Results and Discussions

In order to understand the physical insight of the problem, numerical computations for the velocity, temperature, concentration, skin friction, Nusselt number, and Sherwood number are made for various physical parameters and presented in Figures 2–11. Since water and air are the most commonly occurring fluids in nature, we basically restricted our observations to Pr=0.71 (air) and Pr=7 (water) only.

Figure 2 

Velocity profiles for different 𝐾 and Pr at Gr = 5, Gc = 5, 𝑡=0.6, and Sc = 0.6.

Figure 3 

Effect of Sc on velocity profiles for Pr = 0.71 and 7 at 𝐾=0.6, Gr = 5, Gc = 5, and 𝑡=0.6.

Figure 4 

Effects of Gr and Gc on velocity profiles at 𝐾=0.5, Sc = 0.6, 𝑡=0.6, and Pr = 0.71.

Figure 5 

Effect of time on velocity profiles at 𝐾=0.5, Gr = 4, Gc = 4, Sc = 0.6, and Pr = 0.71.

Figure 6 

Temperature profiles for different values of the Prandtl number and time.

Figure 7 

Effects of 𝐾 and Sc on concentration profiles at 𝑡=0.4.

Figure 8 

Concentration profiles with respect to time for different values of 𝐾 and Sc at 𝑅=1.5.

Figure 9 

Effects of 𝐾 and Pr on skin friction at Gr = 2, Gc = 2, and Sc = 0.6.

Figure 10 

Effects of Gr and Gc on skin friction at 𝐾=0.5, Pr = 0.71, and Sc = 0.6 and 0.96.

Figure 11 

Effects of Pr on the Nusselt number with respect to time.

Figures 2 and 3 depict the effects of the chemical reaction and Schmidt number on velocity profiles for Pr=0.71 (air) and Pr=7 (water) at Gr=5, Gc=5, and 𝑡=0.6. For the case of 𝐾>0, that is, for destructive reaction, increasing values of 𝐾 lead to a fall in velocity profiles. For generative reaction, 𝐾<0, a fall in velocity is observed for increasing 𝐾. This is due to the fact that as 𝐾<0, the last term in the momentum equation becomes positive and plays a crucial role. Also, velocity profile decreases very fast with generative chemical reaction in comparison with destructive chemical reaction for both air and water. The Schmidt number Sc embodies the ratio of the momentum diffusivity to the species diffusivity. It physically relates the relative thickness of the hydrodynamic boundary layer and concentration boundary layer. It is observed that as the Schmidt number increases, the transient velocity decreases. Effects of the thermal Grashof number (Gr) and mass Grashof number (Gc) at 𝐾=0.5, Sc = 0.6, 𝑡=0.6, and Pr=0.71 are presented in Figure 4. It is observed that an increase in the thermal Grashof number or mass Grashof number leads to a rapid increase in the fluid velocity. The thermal Grashof number signifies the relative effect of the buoyancy force to the hydrodynamic viscous force. The positive values of Gr correspond to cooling of the cylinder by natural convection. Heat is, therefore, conducted away from the vertical cylinder into the fluid which increases temperature and thereby enhances the buoyancy force. It is observed that the transient velocity accelerates due to enhancement in the thermal buoyancy force. The solutal Grashof number Gc defines the ratio of the species buoyancy force to the viscous hydrodynamic force. It is noticed that the transient velocity increases considerably with a rise in the species buoyancy force. In both cases, it is interesting to note that as Gr or Gc increases, there is rapid rise in the velocity near the surface of vertical cylinder, and then it descends smoothly to the free stream velocity. Velocity profiles for different values of time at Gr = 4, Gc = 4, 𝐾=0.5, Sc = 0.6, and Pr = 0.71 are plotted in Figure 5. It is noticed that fluid velocity increases with increasing time.

The transient temperature profiles for different Pr and 𝑡 are plotted in Figure 6. It is noticed that the temperature decreases with an increase in Pr and increases with an increase in time 𝑡. The transient concentration profiles for different Schmidt number and chemical reaction parameter are plotted in Figure 7. It is observed from the figure that the concentration decreases with the Schmidt number. For destructive reaction (𝐾>0), as the reaction parameter increases, the concentration profiles decrease. This is due to the fact that as 𝐾>0, the last term in the mass diffusion equation becomes positive for generative reaction (𝐾<0); the opposite effect is observed.

Figure 8 depicts the concentration profiles for various values of the chemical reaction parameter and Schmidt number with respect to time. It is clear from the figure that, initially, concentration profiles increase with time, but, for larger time, they become steady for destructive chemical reaction or even without chemical reaction. It is interesting to see that for generative chemical reaction (𝐾<0), concentration profile increases unboundedly for all values of the Schmidt number.

Figure 9 demonstrates the skin friction against time for various values of the chemical reaction parameter 𝐾 and Prandtl number Pr at Gr = 2, Gc = 2, and Sc = 0.6, which reflects that skin friction increases with the destructive chemical reaction (𝐾>0) or Prandtl number. It is noticeable that for generative chemical reaction (𝐾<0), the skin friction shows an opposite character and decreases as time progresses. Effects of the thermal Grashof number, mass Grashof number, and Schmidt number are plotted in Figure 10. It is observed that an increase in the value of the Schmidt number increases skin friction, whereas it decreases with the increase in the value of the thermal Grashof number or mass Grashof number.

Figure 11 depicts the effects of the Prandtl number on rate of heat transfer with respect to time, which shows that, initially, the Nusselt number decreases sharply with time and tends to steady at larger time. Also, the Nusselt number increases with increase in the Prandtl number.

Effects of the chemical reaction parameter and Schmidt number on the Sherwood number are depicted in Figure 12. It can be observed from the figure that, initially, the Sherwood number falls sharply with time but tends to steady state for destructive chemical reaction or no chemical reaction, whereas, for generative chemical reaction (𝐾<0), the Sherwood number decreases unboundedly. Also, the Sherwood number increases with increased values of the Schmidt number.

Figure 12 

Effects of 𝐾 and Sc on the Sherwood number with respect to time.

5. Conclusions

In this paper, the authors have analyzed the effect of first order chemical reaction on unsteady free convection flow past an infinite vertical cylinder with combined heat and mass transfer. The closed-form solutions of the governing partial differential equations are obtained by the Laplace transform technique. Effects of various physical parameters involved in the flow problem are also analyzed and discussed with graphical results. Based on their observations and discussions above, they summarize their conclusions as follows.(i)Transient velocity decreases with chemical reaction. Velocity profile increases with increase in buoyancy force parameters but decreases with increase in the Prandtl number or Schmidt number. (ii)Transient temperature profiles increase with the Prandtl number. (iii)Transient concentration profile decreases with the Schmidt number and destructive chemical reaction but increases with generative chemical reaction. (iv)Skin friction increases with the destructive chemical reaction or Schmidt number but decreases with the thermal and mass Grashof number.(v)The Nusselt number increases with the Prandtl number. (vi)The Sherwood number increases with increase in destructive chemical reaction and tends to steady state at larger time, while, for generative chemical reaction, there is a decrease in the Sherwood number.

Nomenclature

𝐶′:	Concentration
𝐶:	Dimensionless concentration
𝐷:	Mass diffusion coefficient
Gr:	The thermal Grashof number
Gc:	The mass Grashof number
𝑔:	Acceleration due to gravity
𝐽0:	Bessel’s function of the first kind and order zero
𝐽1:	Bessel’s function of the first kind and order one
𝑘1:	Chemical reaction parameter
𝐾:	Dimensionless chemical reaction parameter
𝐾0:	Modified Bessel’s function of the second kind and order zero
𝐾1:	Modified Bessel’s function of the second kind and order one
Nu:	The Nusselt number
Pr:	The Prandtl number
𝑟:	Radial coordinate measured from the axis of the cylinder
𝑅:	Dimensionless radial coordinate
Sc:	The Schmidt number
Sh:	The Sherwood number
𝑡′:	Time
𝑡:	Dimensionless time
𝑇′:	Temperature of fluid
𝑇:	Dimensionless fluid temperature
𝑢:	Vertical velocity
𝑈:	Dimensionless vertical velocity
𝑌0:	Bessel’s function of the second kind and order zero
𝑌1:	Bessel’s function of the second kind and order one.

Greek Symbols

𝛼:	Thermal diffusivity of fluid
𝜈:	Kinematic viscosity
𝛽:	Volumetric coefficient of thermal expansion
𝛽∗:	Volumetric coefficient of expansion with concentration.