Influence of thermal buoyancy on vortex shedding behind a circular cylinder in parallel flow
Introduction
One of the most interesting phenomena observed in viscous flows past bluff bodies is that of vortex shedding, with the wake structure predominantly called a von-Kármán vortex street. The characteristics of vortex shedding behind circular cylinders or similar two-dimensional obstacles have been extensively discussed; relevant reviews were provided by Bergere and Wille [1], and later by Williamson [2]. This phenomenon has significant implications for numerous thermal applications such as chimney stacks, cooling towers, fuel rods of nuclear reactors, and heat exchanger tubes, where a circular cylinder exists as one of the fundamental structural components. For applications, it must be recognised that buoyancy forces can influence the thermal and flow fields considerably under moderate Reynolds numbers and with a relatively high temperature discrepancy between the cylinder surface and flowing fluid. Indeed, the situation becomes physically complex when the heat transfer strongly influences the nature of the wake behind the hot cylinder owing to the temperature-induced buoyancy forces.
In mixed convective flows past cylinders, the heat transfer characteristics depend mainly on the Reynolds, Grashof or Richardson, and Prandtl numbers, and the angle between the forced flow and buoyancy force direction. These flows can be classified on the basis of the direction of the main flow to the buoyancy-induced flow into two primary forms, horizontal cross flows and vertical flows. In particular, the wake flow behaviour is completely different for the situation of a vertical flow under the influence of vertical buoyancy compared to that for a horizontal cross-flow configuration. The behaviour behind a circular cylinder within mixed convective horizontal cross flows has been reasonably studied, including in Biswas and Sandip [3], Ankit and Amit [4], Jiansheng and Yunjian [5], Jiansheng and Chen [6], Kumar and Sameen [7], Sumit and Dominic [8], Vivien et al. [9], Aleksyuk and Osiptsov [10], Amit et al. [11], Sonal et al. [12], and Zheng et al. [13].
Moreover, in a vertical flow configuration, the forced flow and buoyancy force can be in the same (parallel-flow) or opposing (counter-flow) directions. In fact, several numerical and experimental studies have been conducted previously for the situation wherein, the forced flow is pointed vertically upward (aiding-flow), which is the case considered in this study. Oosthuizen and Madan [14] conducted pioneering experimental work on the problem of mixed convective parallel flow for Reynolds numbers over the range 100 300, and proposed an empirical relationship for the Nusselt number. Jackson and Yen [15] also established a correlation that was determined to be suitable for high Grashof numbers. Other authors, for example, Acrivos [16], Joshi and Sukhatme [17], Nakai and Okazaki [18], Sparrow and Lee [19], Merkin [20], Badr [21,22], and Syakila et al. [23] numerically studied the steady-state condition of combined convection from a circular cylinder under boundary-layer flows using approximate methods. However, the aforementioned investigations were directed at studying only steady flows, disregarding the periodic flow behaviour and determining the critical Grashof number where the flow abruptly transforms from steady to unsteady.
The transient behaviour of mixed convective flow about a circular cylinder was investigated numerically by Jain and Lohar [24] and Noto and Matsumoto [25]. They revealed that the increase in the cylinder temperature increases the vortex shedding frequency. These studies, however, were performed while the flow was developing; it is not clear if they apply to a fully developed flow. Following this, Farouk and Guceri [26] analysed numerically mixed convective flow about a heated cylinder positioned between vertical parallel insulated plates for a low Reynolds number of 6.20. However, their study focused only on the effect of varying the separation distance between the plates. Ho et al. [27] reported a numerical investigation on buoyancy-aided convective flow around a cylinder placed inside an open vertical insulated duct for Reynolds numbers of 20, 40, and 60. Vortex shedding was not observed at a Reynolds number of 60, presumably owing to the considerable blockage caused by the narrow channel wall spacing. Therefore, they investigated the combined influence of the buoyancy force and channel wall spacing on the steady wake and heat transfer enhancement. It was determined that a considerable augmentation in heat transfer can be obtained through blockage variation.
Singh et al. [28] identified unsteady (full-periodic) flow and temperature fields about a hot/cold cylinder within an adiabatic vertical duct for Richardson numbers over the range −1 1, and at a fixed blockage ratio of 0.25 and Reynolds number of 100. The effect on the cylinder heating/cooling was related to the positive/negative effect of the buoyancy forces. It was demonstrated that the shedding incident persists to characterise the flow for Richardson numbers approaching 0.15; however, the aiding-buoyancy alters the dynamics entirely by stopping the shedding thereafter. They predicted that over the range −1 0.15, the increase in the magnitude of Richardson number tends to marginally decrease the average Nusselt number. However, beyond the critical value of 0.15, the increase in Richardson number serves to enhance the average Nusselt number. A similar problem was also examined by Gandikota et al. [29]. A comprehensive investigation was performed for several Reynolds numbers over the range 50 150, for Richardson numbers over the range −0.5 0.5, and for two blockage ratios of 0.02 and 0.25. For = 100, their results indicated that wake unsteadiness stops entirely at critical Richardson numbers of 0.15 and 0.18 at the above mentioned blockage parameters, with a higher critical value for the lower blockage ratio. They indicated that the proximity of the side wall restricts the fluid flow neighbouring the cylinder, which stabilises against the tendency to shed vortices, and consequently lower cylinder heating is necessitated to stabilise the flow. Remarkably, the average Nusselt number was found to persistently increase with Richardson number over its entire negative and positive ranges in both channels with a lesser gradient below the critical value.
The current literature reveals that only two studies, Chang and Sa [30] and Hatanaka and Kawahara [31], have investigated full-periodic parallel flows (aiding-flows) over a circular cylinder in the absence of side walls (unconfined flows). They analysed the behaviour of near-wake vortices behind a hot and/or cold circular cylinder for Richardson numbers over the range −1 1 and for a constant Reynolds number of 100. These studies used different numerical methods; however, both predicted that the cooling of the cylinder surface strengthens the shear layer and the periodic flow becomes more activated. However, heating was found to suppress the frequent flow into a stable twin-eddy structure at a critical Richardson number of 0.15, identifying a collapse of the Kármán vortex street. Hatanaka and Kawahara [31] did not document the heat transfer enhancement/diminution, as their focus was only on the changing dynamic and thermal patterns. Surprisingly, Chang and Sa [30] reported that the trend of the average Nusselt number demonstrated interesting and sensitive features in the unsteady flow region before the suppression of the vortex shedding. Thus, in the negative range of the Richardson numbers −1 0, the average Nusselt number increases notably; however, in the positive range of the Richardson numbers 0 0.15, it decreases marginally to compose a localised minimum prior to the discontinuity at = 0.15. Moreover, in the region of steady flow with twin stationary vortices connected to the cylinder, for 0.15, the average Nusselt number was demonstrated to increase again.
In the aforementioned investigations of Singh et al. [28], Gandikota et al. [29], and Chang and Sa [30], there are apparently conflicting conclusions concerning the reliance of the wall convective heat transfer on the Richardson number in the unsteady flow regime, where the flow structure is dynamically altered and characterised by the vortex shedding. Therefore, the first objective of the current study is to investigate thoroughly the effect of this parameter, and provide a deeper understanding of the augmentation or diminution of the convective heat transfer with the Richardson number, with special attention near its critical value. Secondly, the role of the Reynolds number was ignored in these investigations. Importantly, a higher critical Richardson number could occur at a higher Reynolds number, or vice versa. This could be an indication that the critical Richardson number is not a universal property for the suppression of vortex shedding; rather, it could be Reynolds number dependent. Furthermore, after the stage of shedding suppression, stable twin vortices are formed and attached to the aft of the cylinder. A complete picture and detailed information regarding the dependency of the Nusselt number on the decaying and then vanishing of these dual vortices by the heating effect have not been provided satisfactorily. For these reasons, the present study further investigates and analyses the transient buoyant upward flows over a heated unconfined circular cylinder to better realise the connection between the temperature-induced buoyancy forces and vortex shedding patterns changing in the combined forced and free convection regime, destruction of stable wakes by the heating influence, and the convection heat transfer characteristics, for different Reynolds numbers.
Section snippets
Physical problem and governing equations
The schematic configuration of the physical problem considered in the present investigation accompanied by the Cartesian coordinate system employed are displayed in Fig. 1. The system is comprised of a stationary unconfined horizontal circular cylinder of diameter D heated to and maintained at a fixed and steady temperature . The cylinder is placed within a vertical upward laminar cold free flow stream of temperature and velocity . The oncoming cold flow interacts with the
Numerical method
Artificial boundaries are positioned at the ends of the computational domain forcing the blockage ratio of the cylinder diameter to the domain width equal to 1/30. The downstream and upstream of the computational domain from the cylinder centre are fixed in this investigation at 45D and 15D, respectively. These values are selected to minimise the impact of the inlet and outlet boundaries on the hydrodynamic/thermal fields in the proximity of the cylinder. An existing in-house solver was used
Results and discussion
In the current computational investigation, laminar mixed convection heat transfer out of a heated horizontal unconfined circular cylinder is investigated within a parallel flow regime, such that the buoyancy forces and forced flow are in the same orientation, i.e. the forced flow is vertically upward directed. In this investigation, the Reynolds number is over the range 20 150, at several values of Grashof number () for each case; hence, the Richardson number is changed
Conclusion
A numerical study was performed to investigate the problem of mixed forced and natural convection from a heated circular cylinder subject to a vertical upward laminar flow of a viscous incompressible fluid. The influence of assisting thermal buoyancy on the flow and thermal fields in general, and on the collapse of vortex shedding, next to the cylinder in particular, for different ranges of Reynolds number, and Grashof or Richardson numbers, was considered. The results confirm that the natural
Declaration of competing interest
None.
Acknowledgement
This research was supported in part by the Monash eResearch Centre and eSolutions-Research Support Services through the use of the MonARCH HPC Cluster.
References (50)
- et al.
Effect of thermal buoyancy on vortex shedding past a circular cylinder in cross-flow at low Reynolds numbers
Int. J. Heat Mass Tran.
(2009) - et al.
Direct numerical simulation of energy separation effect in the near wake behind a circular cylinder
Int. J. Heat Mass Tran.
(2018) On the combined effect of forced and free convection heat transfer in laminar boundary layer flows
Chem. Eng. Sci.
(1966)- et al.
Analysis of mixed convection about a horizontal cylinder
Int. J. Heat Mass Tran.
(1976) Mixed convection from a horizontal circular cylinder
Int. J. Heat Mass Tran.
(1977)Laminar combined convection from a horizontal cylinder-parallel and contra flow regimes
Int. J. Heat Mass Tran.
(1984)- et al.
Hydrodynamics of a particle impact on a wall
Appl. Math. Model.
(2006) Flow patterns around heart valves: a numerical method
J. Comput. Phys.
(1972)- et al.
Three-dimensional instabilities in the wake of a circular cylinder
Exp. Therm. Fluid Sci.
(1996) - et al.
Analysis of forced convection heat transfer from a circular cylinder embedded in a porous medium
Int. J. Therm. Sci.
(2012)
Forced convection from a circular cylinder in pulsating flow with and without the presence of porous media
Int. J. Heat Mass Tran.
Bnard convection from a circular cylinder in a packed bed
Int. Commun. Heat Mass Tran.
Numerical study of heat transfer characteristics on tandem cylinders, inline and staggered tube banks in cross-flow of air
Int. Commun. Heat Mass Tran.
Periodic flow phenomena
Annu. Rev. Fluid Mech.
Vortex dynamics in the cylinder wake
Annu. Rev. Fluid Mech.
Effect of a circular cylinder on separated forced convection at a backward-facing step
Int. J. Therm. Sci.
Heat and fluid flow characteristics of a rectangular channel with a small diameter circular cylinder as vortex generator
Int. J. Therm. Sci.
Heat transfer and flow characteristics of a rectangular channel with a small circular cylinder having slit-vent vortex generator
Int. J. Therm. Sci.
Flow past a moderately heated horizontal cylinder at low Reynolds number
J. Aero. Eng.
Effects of thermal induced buoyancy forces on the vortex shedding of a circular cylinder
Int. Commun. Heat Mass Tran.
Numerical study of mixed convection around a heated circular cylinder
Appl. Mech. Mater.
Cross-buoyancy mixed convection from a heated cylinder placed asymmetrically in a channel
Int. Commun. Heat Mass Tran.
Effect of channel confinement on wake dynamics and forced convective heat transfer past a blunt headed cylinder
Int. J. Therm. Sci.
Airfoil-shaped self-agitator for convective heat transfer enhancement
Int. J. Therm. Sci.
Combined convection heat transfer from horizontal cylinders in air
Trans. ASME J. Heat Trans.
Cited by (14)
Mixed convective transport around tandem circular cylinders in an unconfined medium
2024, Nuclear Engineering and DesignNumerical analysis of flow shedding over an obstacle at low Reynolds number
2023, International Journal of Heat and Fluid FlowAn experimental investigation of the flow characteristics and heat transfer properties of a heated and an unheated vertically positioned circular cylinder at different Reynolds numbers
2023, International Journal of Heat and Mass TransferCitation Excerpt :They depend on the angle between the forced flow and the buoyancy direction. In other words, the Richardson number represents the ratio between thermally induced buoyancy and inertial force [5–10]. Understanding how and when the transition from laminar flow structure to turbulent flow occurs in natural convection boundary layers also plays a critical role in developing a control strategy to increase heat and mass transfer [11].
Study on the fluid flow and heat transfer characteristics of a horizontal elliptical cylinder under thermal buoyancy effect
2022, International Journal of Heat and Mass TransferCitation Excerpt :It was observed that increasing the Grashof number eliminates the vortex shedding, although the drag coefficient also increases. Moreover, Al-Sumaily et al. [36] studied the thermal buoyancy effect on the vortex shedding of a circular cylinder for the Reynolds numbers between 20 and 150. According to the studies mentioned above, the effect of thermal buoyancy on the flow past an elliptical cylinder has yet to be determined.
Thermal Features of mixed convection from an inline periodic array of circular cylinders in non-Newtonian power-law fluids
2022, Case Studies in Thermal Engineering