Research Article  Open Access
Ashish Paul, Rudra Kanta Deka, "Unsteady Natural Convection Flow past an Infinite Cylinder with Thermal and Mass Stratification", International Journal of Engineering Mathematics, vol. 2017, Article ID 8410691, 13 pages, 2017. https://doi.org/10.1155/2017/8410691
Unsteady Natural Convection Flow past an Infinite Cylinder with Thermal and Mass Stratification
Abstract
This paper presents an analytical solution of unsteady onedimensional 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.
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 twodimensional unsteady natural convective flow past semiinfinite 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 onedimensional 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 highporosity 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, semiinfinite vertical cylinder with constant temperature and mass diffusion in a thermally stratified medium by employing a finite difference scheme of CrankNicolson type.
Shapiro and Fedorovish [15] presented analytical solution of onedimensional 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 nonNewtonian 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 onedimensional 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 with constant temperature and concentration in presence of thermal and mass stratification. The axis of the cylinder is taken vertically upward along the axis of the cylinder and the radial coordinate 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 () 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 from the cylinder. Behind the wave there must be a dependence on the vertical coordinate, . The basic premise in this work is that convective effects will begin at a position, , as soon as fluid which was initially located at the leading edge rises to this position, regardless of the distance, , away from the cylinder at which it first arrives. Since the surface temperature and concentration above the leading edge are uniform with , the temperature of the fluid and concentration may be assumed to be independent of . In addition, the vertical velocity, , must be independent of 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 only. At time , the uniform temperature () and concentration () 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 and termed as thermal stratification and mass stratification, respectively. Here, and 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 selfstratifications. Then following Boussinesq’s approximation, the onedimensional equations for momentum, energy, and concentration are as follows:with initial and boundary conditions as
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 for Boussinesq flow of liquids or gases, for a perfect gas with pressure work term retained], and is height. Thus is the combination of thermal stratification and compression, where stands for thermal stratification and for compression. In an adiabatic environment, . In that case and () 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 (stable), (neutral), and (unstable). It should be noted, however, that the compression work term is generally quite small (since = 9.8 m s^{−2}, = 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 is termed as mass stratification. We shall also study the simultaneous effect of mass stratification in addition to thermal stratification.
Introducing the nondimensional quantitiesthe governing equations (1) reduce toand the corresponding initial and boundary conditions in nondimensional form are
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 becomewhere 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 , , and are presented in Appendix)
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: where . It is to be noted that and in expressions , , 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 is obtained from the velocity profile (10) as
3.2. Nusselt Number
The nondimensional Nusselt number is obtained from the temperature profile (11) as
3.3. Sherwood Number
The nondimensional Sherwood number is obtained from the concentration profile (12) aswhere
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
4.1. Skin Friction
Nondimensional skin friction is obtained from velocity profile (17) as
4.2. Nusselt Number
Nondimensional Nusselt number is obtained from temperature profile (18) as
4.3. Sherwood Number
Nondimensional Sherwood number is obtained from concentration profile (19) as
It is to be noted that as 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 , from the surface of the cylinder where can be obtained from This leads to the determination of from the following expression:A list of calculated values of for different values , , Gr, and Gc is tabulated in Table 1.

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 and for Gr = Gc = 5 and and Figure 3 shows the effects of Gr and Gc for , , and . It is observed from these figures that velocity decreases with the increase in or increase in 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 . Remarkably, in our analysis, we have seen that, for mass stratification , 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 , and Gr, Gc on velocity profiles against time at . The time required to reach the steady state increases with decrease in stratification parameters or and increases with increase in the values of Gr or Gc.
Figure 6 shows the effects of thermal and mass stratification on temperature profiles at Gr = Gc = 5 at . It is observed from the figure that temperature decreases with the increase in the stratification parameter or . 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 , for and it shows that temperature decreases with increase in Gr or Gc. Figures 8 and 9, respectively, show the effects of , and Gr, Gc on temperature profile against time at . 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 , , 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 thermal and mass stratification on concentration profiles are shown in Figure 10 at Gr = Gc = 5 for . It is observed from the figure that concentration decreases with an increase in stratification parameter or . 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 , for and that concentration decreases with increase in Gr or Gc. Figures 12 and 13, respectively, show the effects of , and Gr, Gc on concentration profiles against time at . 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 , , Gr, or Gc. It is also observed that, in absence of mass stratification parameter, effect of thermal stratification on concentration profile is negligible.
Skin friction represented by Figures 14 and 15 shows the effects of , 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.
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.
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.
7. Conclusions
(i)Solution expressions obtained for unsteady state approaches to the solutions of steady state as becomes large.(ii)The time required to reach steady state velocity increases with increase in Gr or Gc but decreases with increase in or , while time required to reach the steady state temperature and concentration increases with decrease in , , Gr, or Gc.(iii)Velocity increases with increase in Gr or Gc but decreases with increase in or .(iv)Temperature decreases with increase in , , Gr, or Gc. In presence of high stratification temperature becomes negative.(v)Concentration decreases with increase in , , 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 in or .(vii)Rate of heat transfer increases as Gr, , , or increases.(viii)Rate of mass transfer increases with increase in Gr, , , or .
Appendix
Laplace transform of (4) with initial conditions in (5) and with is (6), (7), and (8) with transformed boundary conditionswhere , , and are Laplace transforms of , , and , respectively, and is Laplace transform parameter.
Upon applying (7) and (8), we eliminate and from (6) that give rise to where , , .
Now to determine from (A.2), we rewrite (A.2) as with
Equation (A.3) is the modified Bessel equation of order zero and its solution is where and are arbitrary constants. In view of (A.5) and using the following identities (see Carslaw and Jaeger [25]) the expression for is obtained after applying variation of parameter technique aswhere and are arbitrary constants.
Now, since and are unbounded as , we set and so (A.7) is then read asNow, using (A.8) we have from (6) Using boundary conditions (A.1) in (A.8) and (A.9) and solving we haveThus we haveSimilarly, substituting the expression of in (7) and (8) and solving, the expressions of and can be obtained as
Nomenclature
Concentration  
Concentration at the surface of the cylinder  
Mass diffusion coefficient  
Thermal Grashof number  
Mass Grashof number  
Acceleration due to gravity  
Bessel function of first kind and order zero  
Bessel function of first kind and order one  
Dimensionless mass stratification parameter  
Modified Bessel function of second kind and order zero  
Modified Bessel function of second kind and order one  
Nusselt number  
Prandtl number  
Radial coordinate measured from the axis of the cylinder  
Radius of the cylinder  
Dimensionless radial coordinate  
Dimensionless thermal stratification parameter  
Schmidt number  
Sherwood number  
Time  
Dimensionless time  
Temperature  
Temperature at the surface of the cylinder  
component of velocity  
Dimensionless velocity  
Bessel function of second kind and order zero  
Bessel function of second kind and order one  
Thermal diffusivity  
Kinematic viscosity  
Volumetric coefficient of thermal expansion  
Volumetric coefficient of expansion with concentration  
Dimensionless temperature  
Dimensionless concentration  
Thermal stratification parameter  
Mass stratification parameter . 
Competing Interests
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.
References
 B. Gebhart and L. Pera, “The nature of vertical natural convection flows resulting from the combined buoyancy effects of thermal and mass diffusion,” International Journal of Heat and Mass Transfer, vol. 14, no. 12, pp. 2025–2050, 1971. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 F. A. Bottemanne, “Experimental results of pure and simultaneous heat and mass transfer by free convection about a vertical cylinder for Pr=0.71 and Sc=0.63,” Applied Scientific Research, vol. 25, no. 5, pp. 372–382, 1972. View at: Publisher Site  Google Scholar
 T. S. Chen and C. F. Yuh, “Combined heat and mass transfer in natural convection along a vertical cylinder,” International Journal of Heat and Mass Transfer, vol. 23, no. 4, pp. 451–461, 1980. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 K. Velusamy and V. K. Garg, “Transient natural convection over a heat generating vertical cylinder,” International Journal of Heat and Mass Transfer, vol. 35, no. 5, pp. 1293–1306, 1992. View at: Publisher Site  Google Scholar
 P. Ganesan and H. P. Rani, “Transient natural convection along vertical cylinder with heat and mass transfer,” Heat and Mass Transfer, vol. 33, no. 56, pp. 449–455, 1998. View at: Publisher Site  Google Scholar
 P. Ganesan and P. Loganathan, “Effects of mass transfer and flow past a moving vertical cylinder with constant heat flux,” Acta Mechanica, vol. 150, no. 34, pp. 179–190, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. Ganesan and P. Loganathan, “Unsteady natural convective flow past a moving vertical cylinder with heat and mass transfer,” Heat and Mass Transfer, vol. 37, no. 1, pp. 59–65, 2001. View at: Publisher Site  Google Scholar
 P. Ganesan and P. Loganathan, “Heat and mass flux effects on a moving vertical cylinder with chemically reactive species diffusion,” Journal of Engineering Physics and Thermophysics, vol. 75, no. 4, pp. 899–909, 2002. View at: Publisher Site  Google Scholar
 H. P. Rani, “Transient natural convection along a vertical cylinder with variable surface temperature and mass diffusion,” Heat and Mass Transfer, vol. 40, no. 12, pp. 67–73, 2003. View at: Publisher Site  Google Scholar
 M. Devakar, D. Sreenivasu, and B. Shankar, “Analytical solutions of couple stress fluid flows with slip boundary conditions,” Alexandria Engineering Journal, vol. 53, no. 3, pp. 723–730, 2014. View at: Publisher Site  Google Scholar
 R. K. Deka, A. Paul, and A. Chaliha, “Transient free convection flow past vertical cylinder with constant heat flux and mass transfer,” Ain Shams Engineering Journal, 2015. View at: Publisher Site  Google Scholar
 H. S. Takhar, A. J. Chamkha, and G. Nath, “Natural convection flow from a continuously moving vertical surface immersed in a thermally stratified medium,” Heat and Mass Transfer/Waerme und Stoffuebertragung, vol. 38, no. 1, pp. 17–24, 2001. View at: Publisher Site  Google Scholar
 H. S. Takhar, A. J. Chamkha, and G. Nath, “Natural convection on a vertical cylinder embedded in a thermally stratified highporosity medium,” International Journal of Thermal Sciences, vol. 41, no. 1, pp. 83–93, 2002. View at: Publisher Site  Google Scholar
 P. Loganathan and P. Ganesan, “Numerical study of doublediffusive, free convective flow past a moving vertical cylinder,” Journal of Engineering Physics and Thermophysics, vol. 79, no. 1, pp. 73–78, 2006. View at: Publisher Site  Google Scholar
 A. Shapiro and E. Fedorovich, “Unsteady convectively driven flow along a vertical plate immersed in a stably stratified fluid,” Journal of Fluid Mechanics, vol. 498, pp. 333–352, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Shapiro and E. Fedorovich, “Prandtl number dependence of unsteady natural convection along a vertical plate in a stably stratified fluid,” International Journal of Heat and Mass Transfer, vol. 47, no. 22, pp. 4911–4927, 2004. View at: Publisher Site  Google Scholar
 C.Y. Cheng, “Combined heat and mass transfer in natural convection flow from a vertical wavy surface in a powerlaw fluid saturated porous medium with thermal and mass stratification,” International Communications in Heat and Mass Transfer, vol. 36, no. 4, pp. 351–356, 2009. View at: Publisher Site  Google Scholar
 D. Srinivasacharya and G. S. Reddy, “Effect of double stratification on free convection in a powerlaw fluid saturated porous medium,” Fluid Dynamics and Materials Processing, vol. 9, no. 3, pp. 291–305, 2013. View at: Google Scholar  MathSciNet
 D. Srinivasacharya and C. RamReddy, “Natural convection heat and mass transfer in a micropolar fluid with thermal and mass stratification,” International Journal for Computational Methods in Engineering Science and Mechanics, vol. 14, no. 5, pp. 401–413, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 B. V. Rathish Kumar and S. V. S. S. N. V. G. Krishna Murthy, “Double diffusive free convection induced by verticalwavy surface in a doubly stratified Darcy porous medium under the influence of Soret and Dufour effect,” Journal of Porous Media, vol. 15, no. 9, pp. 877–890, 2012. View at: Publisher Site  Google Scholar
 M. Neagu, “Free convective heat and mass transfer induced by a constant heat and mass fluxes vertical wavy wall in a nonDarcy double stratified porous medium,” International Journal of Heat and Mass Transfer, vol. 54, no. 1112, pp. 2310–2318, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. K. Deka and A. Paul, “Transient free convection flow past an infinite moving vertical cylinder in a stably stratified fluid,” Journal of Heat Transfer, vol. 134, no. 4, Article ID 042503, 8 pages, 2012. View at: Publisher Site  Google Scholar
 R. K. Deka and A. Paul, “Transient free convection flow past an infinite vertical cylinder with thermal stratification,” Journal of Mechanical Science and Technology, vol. 26, no. 8, pp. 2229–2237, 2012. View at: Publisher Site  Google Scholar
 R. K. Deka and A. Paul, “Convectively driven flow past an infinite moving vertical cylinder with thermal and mass stratification,” Pramana—Journal of Physics, vol. 81, no. 4, pp. 641–665, 2013. View at: Publisher Site  Google Scholar
 H. S. Carslaw and J. C. Jaeger, Operational Methods in Applied Mathematics, Oxford University Press, London, UK, 2nd edition, 1948.
 A. K. Kulkarni, H. R. Jacobs, and J. J. Hwang, “Similarity solution for natural convection flow over an isothermal vertical wall immersed in thermally stratified medium,” International Journal of Heat and Mass Transfer, vol. 30, no. 4, pp. 691–698, 1987. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Ashish Paul and Rudra Kanta Deka. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.