Features of the supercritical transition in the wake of a circular cylinder

The features of the wake behind a uniform circular cylinder atRe=200, which is just beyond the critical Reynolds number of 3-D transition, are investigated in detail by direct numerical simulations by solving 3-D incompressible Navier-Stokes equations using mixed spectral-spectral-element method. The high-order splitting algorithm based on the mixed stiffly stable scheme is employed in the time discretization. Due to the nonlinear evolution of the secondary instability of the wake, the spanwise modes with different wavelengths emerge. The spanwise characteristic length determines the transition features and global properties of the wake. The existence of the spanwise phase difference of the primary vortices shedding is confirmed by Fourier analysis of the time series of the spanwise vorticity and attributed to the dominant spanwise mode. The spatial energy distributions of various modes and the velocity profiles in the near wake are obtained. The numerical results indicate that the near wake is in 3-D quasi-periodic laminar state with transitional behaviors at this supercritical Reynolds number. The project supported by the State Key Fundamental Research Project of “Large Scale Scientific Computation Research” (G199903281)     

Three-dimensional numerical simulations of cross-flow around four cylinders in an in-line square configuration

This paper presents a numerical study of three-dimensional (3-D) laminar flow around four circular cylinders in an in-line square configuration. The investigation focuses on effects of spacing ratio (L/D) and aspect ratio (H/D) on 3-D flow characteristics, and the force and pressure coefficients of the cylinders. Extensive 3-D numerical simulations were performed at Reynolds number of 200 for L/D from 1.6 to 5.0 at H/D=16 and H/D from 6 to 20 at L/D=3.5. The results show that the 3-D numerical simulations have remedied the inadequacy of 2-D simulations and the results are in excellent agreement with the experimental results. The relation between 3-D flow patterns and pressure characteristics around the four cylinders is examined and discussed. The critical spacing ratio for flow pattern transformation was found to be L/D=3.5 for H/D=16, while a bistable wake pattern was observed at L/D=1.6 for the same aspect ratio. Moreover, a transformation of flow pattern from a stable shielding flow pattern to a vortex shedding flow pattern near the middle spanwise positions of the cylinders was observed and was found to be dependent on the aspect ratio, spacing ratio, and end wall conditions. Due to the highly 3-D nature of the flows, different flow patterns coexist over different spanwise positions of the cylinders even for the same aspect ratio. It is concluded that spacing ratio, aspect ratio, and the no-slip end wall condition have important combined effects on free shear layer development of the cylinders and hence have significant effects on the pressure field and force characteristics of the four cylinders with different spacing ratios and aspect ratios.     

Transition phenomena in the wake of an inclined square cylinder

The division of flow regimes in a square cylinder wake at various angles of attack (α) is studied. This study provides evidence of the existence of modes A and B instabilities in the wake of an inclined square cylinder. The critical Reynolds numbers for the inception of these instability modes were identified through the determination of discontinuities in the Strouhal number versus Reynolds number curves. The spectra and time traces of wake streamwise velocity were observed to display three distinct patterns in different flow regimes. Streamwise vortices with different wavelengths at various Reynolds numbers were visualized. A PIV technique was employed to quantitatively measure the parameters of wake vortices. The wavelengths of the streamwise vortices in the modes A and B regimes were measured by using the auto-correlation method. From the present investigation, the square cylinder wake at various angles of attack undergoes a similar transition path to that of a circular cylinder, although various quantitative parameters measured which include the critical Reynolds numbers, spanwise wavelength of secondary vortices, and the circulation and vorticity of wake vortices all show an α dependence.     

A three-dimensional simulation of a steady approach flow past a circular cylinder at low Reynolds number

The three-dimensional (3D) unsteady viscous wake of a circular cylinder exposed to a steady approach flow is calculated using a fractional-step finite-difference/spectral-element method. The calculated flow fields at Reynolds numbers of 100 (2D) and 200 (3D) are examined in detail. The flow field at Re = 100 is 2D as expected, while the flow field at Re = 200 has distinct 3D features, with spanwise wavelengths of about 3.75 cylinder diameters. The calculated results produce drag and lift coefficients and Strouhal numbers that agree extremely well with the experimental values. These 3D values at Re = 200 are in better agreement with experimental values than the results of a 2D calculation at Re = 200, which is expected. © 1998 John Wiley & Sons, Ltd.     

Numerical Simulation of the Flow Around an Infinitely Long Circular Cylinder in the Transition Regime

The development of three-dimensional structures and the succeeding transition to turbulence occurs in the wake of a circular cylinder at Reynolds numbers 190≤Re≤330. This regime is investigated numerically by means of a spectral element method. Earlier numerical works aimed mainly at reproducing characteristic wake patterns observed in experiments. Small sizes of computational domains and short integration times were chosen to save computational resources. Consequently, the quantitative results show a considerable scatter. Within this work, a step by step approach to highly accurate direct numerical simulations is described. Thorough studies of the effect of resolution and blockage are performed in the laminar, two-dimensional regime, resulting in Reynolds number relationships that exactly reproduce experimental data. Based on these results, a stability analysis is performed to obtain wavelengths that are unstable against spanwise perturbations and the critical Reynolds number for the onset of three-dimensionality. The most unstable wavelengths of the “mode A” and “mode B” instabilities and its multiples are used as periodicity length for direct numerical simulations. Effects of integration time, resolution in streamwise as well as spanwise directions, and periodicity length on the flow quantities are studied. Numerically obtained Reynolds number relationships of Strouhal number and base-pressure coefficients that fit accurately within measured results are given for the first time. Curves for drag and lift coefficients are provided and compared with previous numerical studies. Furthermore, physical interpretations of the wake transition are discussed. Since the separation of physical features and effects of experimental arrangements are frequently an open question, our numerical results are able to supply a contribution to the understanding of the physics of cylinder flow. Received 12 September 2000 and accepted 26 June 2001     

LES-based Study of the Roughness Effects on the Wake of a Circular Cylinder from Subcritical to Transcritical Reynolds Numbers

This paper investigates the effects of surface roughness on the flow past a circular cylinder at subcritical to transcritical Reynolds numbers. Large eddy simulations of the flow for sand grain roughness of size k/D = 0.02 are performed (D is the cylinder diameter). Results show that surface roughness triggers the transition to turbulence in the boundary layer at all Reynolds numbers, thus leading to an early separation caused by the increased momentum deficit, especially at transcritical Reynolds numbers. Even at subcritical Reynolds numbers, boundary layer instabilities are triggered in the roughness sublayer and eventually lead to the transition to turbulence. The early separation at transcritical Reynolds numbers leads to a wake topology similar to that of the subcritical regime, resulting in an increased drag coefficient and lower Strouhal number. Turbulent statistics in the wake are also affected by roughness; the Reynolds stresses are larger due to the increased turbulent kinetic energy production in the boundary layer and separated shear layers close to the cylinder shoulders.     

Effects of a splitter plate on the near wake of a circular cylinder in 2 and 3-dimensional flow configurations

 Results of flow visualization, hot wire, and base pressure measurements were conducted for an investigation of the near wake of a circular cylinder at subcritical Reynolds numbers between 2700 to 46000. A base mounted splitter plate allowed for the modification of the formation region characteristics without disrupting the normal Kármán shedding. The results provide an explanation for the non-linearity in the relationship between shedding frequency and splitter plate length and extend the previous investigations of Roshko (1954), Gerrard (1966) and Apelt et al. (1973). In addition to the nominal 2-D configurations, a sinuous trailing edge splitter plate, cylinder taper, and shear flow were incorporated to study the effects of mild 3-dimensionality. A strong spanwise coherence was found to exist in the formation region. A superposition principle was discovered which showed that certain 3-D geometry and flow configurations could be combined to produce a nominal 2-D wake. Received: 26 June 1996 / Accepted: 7 January 1997     

Momentum and heat transport in a finite-length cylinder wake

This paper reports an experimental study of turbulent momentum and heat transport in the wake of a wall-mounted finite-length square cylinder, with its length-to-width ratio L/d = 3–7. The cylinder was slightly heated so that heat produced could be considered as a passive scalar. A moveable three-wire probe (a combination of an X-wire and a cold wire) was used to measure velocity and temperature fluctuations at a Reynolds number of 7,300 based on d and the free-stream velocity. Measurements were performed at 10 and 20d downstream of the cylinder at various spanwise locations. Results indicate that L/d has a pronounced effect on Reynolds stresses, temperature variance and heat fluxes. The downwash flow from the free end of the cylinder acts to suppress spanwise vortices and, along with the upwash flow from the cylinder base, makes the finite-length cylinder wake highly three-dimensional. Reynolds stresses, especially the lateral normal stress, are significantly reduced as a result of suppressed spanwise vortices at a small L/d. The downwash flow acts to separate the two rows of spanwise vortices further apart from the wake centerline, resulting in a twin-peak distribution in temperature variance. While the downwash flow entrains high-speed fluid into the wake, responsible for a small deficit in the time-averaged streamwise velocity near the free end, it does not alter appreciably the distribution of time-averaged temperature. It has been found that the downwash flow gives rise to a counter-gradient transport of momentum about the central region of the wake near the free end of the cylinder, though such a counter-gradient transport does not occur for heat transport.     

Simulation of turbulent flows around a circular cylinder using nonlinear eddy-viscosity modelling: steady and oscillatory ambient flows

Steady incident flow past a circular cylinder for sub- to supercritical Reynolds number has been simulated as an unsteady Reynolds-averaged Navier–Stokes (RANS) equation problem using nonlinear eddy-viscosity modelling assuming two-dimensional flow. The model of Craft et al. (Int. J. Heat Fluid Flow 17 (1996) 108), with adjustment of the coefficients of the ‘cubic’ terms, predicts the drag crisis at a Reynolds number of about 2×10⁵ due to the onset of turbulence upstream of separation and associated changes in Strouhal number and separation positions. Slightly above this value, at critical Reynolds numbers, drag is overestimated because attached separation bubbles are not simulated. These do not occur at supercritical Reynolds numbers and drag coefficient, Strouhal number and separation positions are in approximate agreement with experimental measurements (which show considerable scatter). Fluctuating lift predictions are similar to sectional values measured experimentally for subcritical Reynolds numbers but corresponding measurements have not been made at supercritical Reynolds numbers. For oscillatory ambient flow, in-line forces, as defined by drag and inertia coefficients, have been compared with the experimental values of Sarpkaya (J. Fluid Mech. 165 (1986) 61) for values of the frequency parameter, β=D²/νT, equal to 1035 and 11240 and Keulegan–Carpenter numbers, KC=U₀T/D, between 0.2 and 15 (D is cylinder diameter, ν is kinematic viscosity, T is oscillation period, and U₀ is the amplitude of oscillating velocity). Variations with KC are qualitatively reproduced and magnitudes show best agreement when there is separation with a large-scale wake, for which the turbulence model is intended. Lift coefficients, frequency and transverse vortex shedding patterns for β=1035 are consistent with available experimental information for β≈250−500. For β=11240, it is predicted that separation is delayed due to more prominent turbulence effects, reducing drag and lift coefficients and causing the wake to be more in line with the flow direction than transverse to it. While these oscillatory flows are highly complex, attached separation bubbles are unlikely and the flows probably two dimensional.     

Numerical investigations into effects of three-dimensional wake patterns on unsteady aerodynamic characteristics of a circular cylinder at Re=1.3×105

In contrast with a wide range of applications concerning flows around a circular cylinder at upper subcritical Reynolds numbers (Re), there is no systematic understanding about the fundamentals of so-called random flow patterns, and their effects on intermittent modulations in the time history of pressure or force, and the decrease in their spanwise correlations. This paper employed the large-eddy simulation (LES) technique to predict flows past a circular cylinder at Re=1.3×10⁵ and to provide images based on flow visualization that can clarify the physical mechanism responsible for these outcomes. A reasonably sufficient spanwise length was adopted for the numerical model by taking into consideration the effect of aspect ratios (the spanwise length to the diameter). We found that even at such high Res, a three-dimensional pattern of vortical field is present in the wake resulting in total force modulation and weak spanwise correlation, e.g., obvious oblique shedding. The whole development process of the three-dimensional wake is exhibited as a universal. The results revealed that local phase variations in primary vortex shedding are the starting points of three-dimensional wake patterns, which are induced by the “irregular” streamwise vortex. The three-dimensional near wake following local phase variations is associated with a successive evolution composed of certain stages in order. Quantitative analyses based on the time series of sectional lift coefficients show that intermittent increase in primary shedding periods and sectional lift streak divisions are closely related to local phase variations and vortex division in the development process of the three-dimensional pattern. In addition to the phase difference along the span, the three-dimensional pattern also weakens vortex shedding in cross sections perpendicular to the axis of the cylinder, resulting in modulation of the sectional lift coefficient.     

Separated flow structures around a cylindrical obstacle in a narrow channel

A detailed experimental study is performed on the separated flow structures around a low aspect-ratio circular cylinder (pin-fin) in a practical configuration of liquid cooling channel. Distinctive features of the present arrangement are the confinement of the cylinder at both ends, water flow at low Reynolds numbers (Re = 800, 1800, 2800), very high core flow turbulence and undeveloped boundary layers at the position of the obstacle. The horseshoe vortex system at the junctions between the cylinder and the confining walls and the near wake region behind the obstacle are deeply investigated by means of Particle Image Velocimetry (PIV). Upstream of the cylinder, the horseshoe vortex system turns out to be perturbed by vorticity bursts from the incoming boundary layers, leading to aperiodical vortex oscillations at Re = 800 or to break-away and secondary vorticity eruptions at the higher Reynolds numbers. The flow structures in the near wake show a complex three-dimensional behaviour associated with a peculiar mechanism of spanwise mass transport. High levels of free-stream turbulence trigger an early instabilization of the shear layers and strong Bloor–Gerrard vortices are observed even at Re = 800. Coalescence of these vortices and intense spanwise flow inhibit the alternate primary vortex shedding for time periods whose length and frequency increase as the Reynolds number is reduced. The inhibition of alternate vortex shedding for long time periods is finally related to the very large wake characteristic lengths and to the low velocity fluctuations observed especially at the lowest Reynolds number.     

Flow past a circular cylinder over a free surface: Interaction between the near wake and the free surface deformation

Air-flow around a circular cylinder placed above a free surface and liquid flow under the free surface were investigated experimentally in a wind/wave tunnel. The cylinder spanned the tunnel test-section and was oriented normal to the freestream direction. The main objective of this study was to investigate the interaction of the cylinder wake with the free surface. The flow structure was analyzed for various gap widths, H, between the cylinder and the free surface using a digital particle image velocimetry (PIV) system with a spatial resolution of 2048×2048 pixels. The Reynolds number based on the cylinder diameter was 3.3×10³. For each experimental condition, 400 instantaneous velocity fields were measured and ensemble-averaged to obtain spatial distributions of the mean velocity and turbulence statistics. The results showed that the cylinder near-wake inclined upward due to the influence of the free surface elevation. Vortices were shed, even at a small gap ratio of H/D=0.25, where D is the cylinder diameter. Strong jet-like flow appeared in the gap beneath the cylinder. At a gap ratio of H/D=0.50, the jet flow exhibited a quasi-periodic vibration with a period of 2–3 s. The free surface deformation was caused by the pressure difference in the air-flow immediately above it. As the gap ratio increased, the inclination angle of the wake and the height of the free surface elevation decreased gradually. The liquid flow under the free surface followed a convective flow motion, and the range of the convection depended on the gap width between the cylinder and the free surface.     

Passive control of wake flow by two small control cylinders at Reynolds number 80

Passive control of the wake behind a circular cylinder in uniform flow is studied by numerical simulation at Re_D=80. Two small control cylinders are placed symmetrically along the separating shear layers at various stream locations. In the present study, the detailed flow mechanisms that lead to a significant reduction in the fluctuating lift but maintain the shedding vortex street are clearly revealed. When the stream locations lie within 0.8≤X_C/D≤3.0, the alternate shedding vortex street remains behind the control cylinders. In this case, the symmetric standing eddies immediately behind the main cylinder and the downstream delay of the shedding vortex street are the two primary mechanisms that lead to a 70–80% reduction of the fluctuating lift on the main cylinder. Furthermore, the total drag of all the cylinders still has a maximum 5% reduction. This benefit is primarily attributed to the significant reduction of the pressure drag on the main cylinder. Within X_C/D>3.0, the symmetry of the standing eddy breaks down and the staggered vortex street is similar to that behind a single cylinder at the same Reynolds number. In the latter case, the mean pressure drag and the fluctuating lift coefficients on the main cylinder will recover to the values of a single cylinder.     

Anisotropic Organised Eddy Simulation for the prediction of non-equilibrium turbulent flows around bodies

The unsteady turbulent flow around bodies at high Reynolds number is predicted by an anisotropic eddy-viscosity model in the context of the Organised Eddy Simulation (OES). A tensorial eddy-viscosity concept is developed to reinforce turbulent stress anisotropy, that is a crucial characteristic of non-equilibrium turbulence in the near-region. The theoretical aspects of the modelling are investigated by means of a phase-averaged PIV in the flow around a circular cylinder at Reynolds number 1.4×10⁵. A pronounced stress–strain misalignment is quantified in the near-wake region of the detached flow, that is well captured by a tensorial eddy-viscosity concept. This is achieved by modelling the turbulence stress anisotropy tensor by its projection onto the principal matrices of the strain-rate tensor. Additional transport equations for the projection coefficients are derived from a second-order moment closure scheme. The modification of the turbulence length scale yielded by OES is used in the Detached Eddy Simulation hybrid approach. The detached turbulent flows around a NACA0012 airfoil (2-D) and a circular cylinder (3-D) are studied at Reynolds numbers 10⁵ and 1.4×10⁵, respectively. The results compared to experimental ones emphasise the predictive capabilities of the OES approach concerning the flow physics capture for turbulent unsteady flows around bodies at high Reynolds numbers.     

Numerical solution for the flow of viscoelastic fluids around an inclined circular cylinder.

Finite difference solutions have been obtained by the perturbation method to investigate the influence of shear thinning and elasticity on the flow around an inclined circular cylinder of finite length in a uniform flow. In this numerical analysis a generalized upper-convected Maxwell model, in which the viscosity changes according to the Cross model, has been used.The local flow over the cylinder is only slightly deflected. However, in the wake flow behind the cylinder the particle path is remarkably influenced by the axial flow and rapidly flows up parallel to the cylinder's axis. Then it gradually rejoins direction of the incoming flow. It is found that viscoelastic fluids are prone to flow axially in the vicinity of the cylinder. The numerical predictions generally agree with the flow visualization results.The numerical solutions also demonstrate that elasticity has a strong effect on the velocity profile especially around both ends of the cylinder; elasticity increases the asymmetric profiles of both circumferential velocity and axial velocity with respect to equal to 90° and decreases a difference in the circumferential velocity between the windward end and the leeward end.For non-Newtonian fluids, the length of the wake flow is influenced by not only the Reynolds number but also the cylinder diameter and it is larger for the cylinder with the smaller diameter at the same Reynolds number.Partly presented at the 9th Australasian Fluid Mechanics Conference, University of Auckland, New Zealand, 8–12 December, 1986     

Effect of Large Scale 3-D Structures on the Flow Around a Heated Cylinder at Low Reynolds Number

This work presents results of flow around a heated circular cylinder in mixed convection regime and demonstrates that Prandtl number and angle of attack of the incoming flow have a large influence on the characterisation of the flow transition from 2-D to 3-D. Previous studies show that heat transfer can enhance the formation of large 3-D structures in the wake of the cylinder for Reynolds numbers between 75 and 127 and a Richardson number larger than 0.35. This transitional mode is generally identified as “mode E”. In this work, we compare the results for water-based flow (large Prandtl number) with the ones for air-based flows (low Prandtl number). The comparison is carried out at two Reynolds numbers (100 and 150) and at a fixed Richardson number of 1. It shows that at the low Reynolds number of 100 the low Prandtl number flow does not enter into transition. This is caused by the impairment of the baroclinic vorticity production provoked by the spanwise temperature gradient. At low Prandtl number temperature gradients are less steep. For an air-based flow at Reynolds number 150, several Richardson numbers have been simulated. In this situation, the flow enters into transition and exhibits the characteristics of “mode E”, with the development of Λ-shaped structures in the near wake and mushroom-like structures in the far wake. It is also observed that the transition is delayed at Richardson number of 0.5. Simulations are also carried to investigate the effect of the angle of attack on the incoming flow on the development of large coherent structures. When the angle of attack is positive, the development of the wake tends to return to a more bi-dimensional configuration, where large scale coherent structures are impaired. In contrast, when the angle of attack is negative, large scale tri-dimensional structures dominate the flow in the wake, but with a very chaotic behaviour and the regular pattern of zero angle of attack is destroyed. The different behaviour of the flow with the variation of the angle of attack is also related to the baroclinic vorticity production, where new terms appear in the equations, leading to a positive effect of the vorticity production in case of a negative angle of attack and the opposite for a positive angle of attack.     

Effect of the domain spanwise periodic length on the flow around a circular cylinder

We assess the effect of the choice of spanwise periodic length on simulations of the flow around a fixed circular cylinder. The Reynolds number is set to 400 because, at this value, both lift coefficient and shedding frequency show significant drop due to three-dimensional flow structures. From the analysis of the three-dimensionalities of the wake and of the integral quantities such as Strouhal number, RMS of lift coefficient and energy contained in the three-dimensional portion of the flow we obtain an estimate of the minimum spanwise length to satisfactorily represent the flow. Furthermore, we observe a distinct wake behavior when the spanwise length is approximately the mode B instability wavelength.     

Numerical study of the effects occurring near a circular cylinder in stratified fluid flows with short buoyancy periods

Flow around a two-dimensional circular cylinder of a stratified fluid with periods buoyancy T_b = 25.2 and 6.28 sec is studied numerically over a wide range of Reynolds and Froude numbers. It is found that in the presence of a perturbation ahead of a cylinder which moves downstream with increasing Reynolds number, the salinity isolines have the shape of a semi-circular comb with sharp teeth. The shape change of the attached waves and the occurrence of fluid layers of different densities in the cylinder wake are studied. In flows with a buoyancy period T_b = 6.28 sec at Reynolds numbers Re < 60, stagnant zones are found in the cylinder wake, and at Re > 60, these zones are absent.     

Experiments on the flow past long circular cylinders in a shear flow

This paper describes an experimental investigation of the flow past circular cylinders, with the mean flow perpendicular to the cylinder axis, at conditions that yield a strong three-dimensional behaviour. The experiments were carried out in the subcritical regime. Long cylinders with end plates were subjected to shear flow with a linear velocity profile in the spanwise direction generated by means of a curved gauze. It was concluded that spanwise cellular structures of vortex shedding emerged in the wake, more clearly for some boundary conditions than others. These structures are characterised by a portion of spanwise length, a cell, having constant shedding frequency over a time average, which implies that there were no vortex dislocations inside that cell during that time. These features were studied using flow visualisation and hot-film anemometry. Spectra of the local shedding frequency are shown, revealing the effect of the shear parameter (=0.02 and 0.04) and aspect ratio L/D (=20.6 and 8) on the stability and geometry of the cells at several Reynolds numbers in the range of 3.13×10³Re_m1.25×10⁴.     

Effect of self-issuing jets along the span on the near-wake of a square cylinder

The study herein focuses on the vortex shedding characteristics and near-wake vorticity patterns of a square cylinder having self-issuing jets through holes along its span. Three different values of spacing between the consecutive holes λ with respect to the cylinder diameter D, i.e., λ/D = 1.5, 3 and 4 are studied experimentally via Digital Particle Image Velocimetry for the Reynolds number range extending from 200 to 1,000. It has been observed that the three-dimensionality of the wake flow depends on the spacing between the holes and Re number. For sufficiently low Reynolds numbers, the jet flows issuing from the holes yield a non-uniform distribution of mean flow characteristics like the shedding frequency and the formation length of vortices along the span of the cylinder when the spacing between jets along centerline is close to wavelength of the naturally existing three-dimensional wake instability. Additionally, for Re number up to 500, the self-issuing jets emanating from the holes show an indirect interaction with shear layers originating from upper and lower separation lines of the cylinder. However, for higher Re numbers of 750 and 1,000, they directly interact with and modify the vortices forming from the cylinder.