Experimental and high-order LES analysis of the flow in near-wall region of a square cylinder

A coupled experimental/numerical analysis of turbulent flow past a square cylinder is performed at the ERCOFTAC Reynolds number Re = U_∞D/ν = 21,400, where U_∞ is the inflow velocity and D the cylinder height. Complementary Laser Doppler Velocimetry (LDV) and high-order large-eddy simulations (LES) approaches, based on a spectral vanishing technique (SVV-LES), provide a comprehensive data base including both instantaneous data and post-processed statistics. Beyond these results, an achievement of the paper is to investigate the coherent structures developing on the sides and in the wake of the cylinder with a special focus on the flow features in the near-wall region. The flow is found to separate at the leading edge of the cylinder with the occurence of three-dimensional Kelvin-Helmholtz (KH) pairings localized in the separating shear layer. The interaction between these KH vortical structures and Von Kármán vortex shedding (VK) in the near wake is discussed based on both visualisations and frequency analysis. In particular, signatures of VK and KH vortical structures are found on velocity time samples.     

Thermal characteristics of the wake shear layers from a slightly heated circular cylinder

In this paper, we investigate the thermal characteristics of wake shear layers generated by a slightly heated circular cylinder. Measurements of the fluctuating temperature were made in the region x/d = 0.6 to x/d = 3 (where x is the downstream distance from the cylinder axis and d is the cylinder diameter) using a single cold-wire probe. The Reynolds number Re was varied in the range 2,600–8,600. For Re = 5,500, simultaneous measurements were made with a rake of 16 cold wires, aligned in the direction of the mean shear, at x/d = 2 and 3. The results indicate that the passive temperature can be an effective marker of various instabilities of the wake shear layers, including the Kelvin–Helmholtz (KH) instability. The temperature data have confirmed the approximate Re ^m dependence of the KH instability frequency (f _KH) with different values of m over different ranges of Re, as reported previously in the literature. However, it is found that this power-law dependence is not exact, and a third-order polynomial dependence appears to fit the data well over the full range of Re. Importantly, it is found that the wake shear-layer instabilities can be grouped into three categories: (1) one with frequencies much smaller than the Bénard–Kármán-vortex shedding frequency, (2) one associated with the vortex shedding and (3) one related to the KH instability. The low-frequency shear-layer instabilities from both sides of the cylinder are in-phase, in contrast to the anti-phase high-frequency KH instabilities. Finally, the observed streamwise decrease in the mean KH frequency provides strong support for the occurrence of vortex pairing in wake shear layers from a circular cylinder, thus implying that both the wake shear layer and a mixing layer develop in similar fashion.     

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.     

Time-averaged topological flow patterns and their influence on vortex shedding of a square cylinder in crossflow at incidence

Flow characteristics around the square cylinder and their influence on the wake properties are studied. Time-averaged flow patterns on the surfaces of square cylinder in a cross-stream at incidence are experimentally probed by surface-oil flow technique and analyzed by flow topology for Reynolds numbers between 3.9×10⁴ and 9.4×10⁴ as the incidence angle changes from 0° to 45°. Vortex shedding characteristics are measured by a single-wire hot-wire anemometer for Reynolds numbers between 5×10³ and 1.2×10⁵. The effects of topological flow patterns on the wake properties then are revealed and discussed. Flows around the square cylinder are identified as three categories: the subcritical, supercritical, and wedge flows according to the prominently different features of the topological flow patterns. The Strouhal number of vortex shedding, turbulence in the wake, and wake width present drastically different behaviors in different characteristic flow regimes. A critical incidence angle of 15° separates the subcritical and supercritical regimes. At the critical incidence angle the wake width and shear-layer turbulence present minimum values. The minimum wake width appearing at the critical incidence angle, which leads to the maximum Strouhal number, is due to the reattachment of one of the separated boundary layer to the lateral face of the square cylinder. If the Strouhal numbers are calculated based on the wake width instead of the cross-stream projection width of cylinder, the data in the subcritical and supercritical regimes are well correlated into two groups, which would approach constants at high Reynolds numbers.     

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.     

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.     

Wake transition for surface mounted rectangular cylinder due to incoming shear flow

Three-dimensional unsteady wake characteristics have been investigated numerically in flow past surface mounted finite-height rectangular cylinder using Open Source Field Operation and Manipulation. Effect of impinging shear (shear intensity, K) on transitional characteristics of wake flow has been studied using iso-Q surfaces for Reynolds number (Re) in the range from 150 to 250. Various flow regimes, such as steady flow, symmetric and asymmetric modes of vortex shedding have been identified based on the values of Re and K for different side ratios (SR) of the cylinder. Unsteady wake oscillations have been analyzed using time signal of transverse velocity component in the wake. These signals have been decomposed into different component signals using Hilbert-Huang transformation (HHT). Variation of frequency and energy density with time of the decomposed signals has been presented in the form of Hilbert spectra. Effects of Re, SR and K on wake oscillation frequency have been illustrated in the form of marginal spectra. Time-delay reconstructions and Poincare sections have been examined to study periodic and aperiodic nature of the wake flow. Non-stationarity associated with the wake fluctuation is quantified in terms of degree of stationarity. Symmetric and asymmetric modes have been confirmed using singular value decomposition of the vorticity field and presented using dynamic modes. Growth rate and frequency of the modes corresponding to symmetric shedding are found to be lower than those for asymmetric shedding. In addition, variation in mean drag coefficient has been reported with change in Re and K for each value of SR.     

On vortex shedding from low aspect ratio dual step cylinders

A dual-step cylinder is comprised of two cylinders of different diameters. A large diameter cylinder (D) with low aspect ratio (L/D) is attached to the mid-span of a small diameter cylinder (d). The present study investigates the effect of Reynolds number (Re_D) and L/D on dual step cylinder wake development for D/d=2, 0.2≤L/D≤3, and two Reynolds numbers, Re_D=1050 and 2100. Experiments have been performed in a water flume facility utilizing flow visualization, Laser Doppler Velocimetry (LDV), and Particle Image Velocimetry (PIV). The results show that vortex shedding occurs from both the large and small diameter cylinders for 1≤L/D≤3 at Re_D=2100 and 2≤L/D≤3 at Re_D=1050. At these conditions, large cylinder vortices predominantly form vortex loops in the wake and small cylinder vortices form half-loop vortex connections. At lower aspect ratios, vortex shedding from the large cylinder ceases, with the dominant frequency in the large cylinder wake attributed to the passage of vortex filaments connecting small cylinder vortices. At these lower aspect ratios, the presence of the large cylinder induces periodic vortex dislocations. Increasing L/D increases the frequency of occurrence of vortex dislocations and decreases the dominant frequency in the large cylinder wake. The identified changes in wake topology are related to substantial variations in the location of boundary layer separation on the large cylinder, and, consequently, changes in the size of the vortex formation region. The results also show that the Reynolds number has a substantial effect on wake vortex shedding frequency, which is more profound than that expected for a uniform cylinder.     

Flow around a cylinder surrounded by a permeable cylinder in shallow water

The change in flow characteristics downstream of a circular cylinder (inner cylinder) surrounded by an outer permeable cylinder was investigated in shallow water using particle image velocimetry technique. The diameter of the inner cylinder and the water height were kept constant during the experiments as d?=?50?mm and h _w ?=?25?mm, respectively. The depth-averaged free-stream velocity was also kept constant as U?=?170?mm/s which corresponded to a Reynolds number of Re_d?=?8,500 based on the inner cylinder diameter. In order to examine the effect of diameter and porosity of the outer cylinder on flow characteristics of the inner cylinder, five different outer cylinder diameters (D?=?60, 70, 80, 90 and 100?mm) and four different porosities (???=?0.4, 0.5, 0.6 and 0.7) were used. It was shown that both porosity and outer cylinder diameter had a substantial effect on the flow characteristics downstream of the circular cylinder. Turbulent statistics clearly demonstrated that in comparison with the bare cylinder (natural case), turbulent kinetic energy and Reynolds stresses decreased remarkably when an outer cylinder was placed around the inner cylinder. Thereby, the interaction of shear layers of the inner cylinder has been successfully prevented by the presence of outer cylinder. It was suggested by referring to the results that the outer cylinder having 1.6????D/d????2.0 and 0.4????D/d????0.6 should be preferred to have a better flow control in the near wake since the peak magnitude of turbulent kinetic energy was considerably low in comparison with the natural case and it was nearly constant for these mentioned porosities ??, and outer cylinder to inner cylinder diameter ratios D/d.     

The Effect of Reynolds Number on the Passive Scalar Field in the Turbulent Wake of a Circular Cylinder

In this paper we investigate by experiments the effect of Reynolds number on a passive scalar (temperature) field in the turbulent wake of a slightly heated circular cylinder. The Reynolds number defined by Re≡U _∞ d/ν (see Nomenclature) is varied from Re= 1200 to Re= 8600. Temperature differential above ambient is chosen to be the passive scalar quantity. Present measurements are conducted using a cold wire (0.63 μm) probe. Results obtained suggest that Reynolds number in general has significant influence on the scalar mixing characteristics in the entire wake flow. Specifically, as Re increases, the mean scalar spreads out more rapidly, the scalar fluctuation intensity increases; however, its variance decays at a lower rate with downstream distance. It is also found that an increase of Re accelerates the streamwise evolution of the scalar probability density function from highly non-Gaussianity to near Gaussianity along the wake centreline. This reflects the reduction in length of the Karman-vortex street caused by an increase of Re.     

Large eddy simulation of bluff body wake in planar shear flow

Large eddy simulation of planar shear flow past a square cylinder has been investigated. Dynamic Smagorinsky model has been used to model subgrid scale stress. The shear parameter, K, namely the nondimensional streamwise velocity gradient in the lateral direction, is 0.0, 0.1 and 0.2. Reynolds number based on the centerline velocity is fixed at Re=21400. The time and span‐averaged velocity components, pressure coefficient, Reynolds stresses for uniform are in good agreement with the literature. In shear flow the calculated flow structure and mean velocity components are shown to be markedly different from those of the uniform flow. With increasing shear parameter, the cylinder wake is dominated by clockwise vortices. Both the velocity components in shear flow are compared with respective components in uniform flow. Comparison of normal and shear stresses between shear and no shear case have also been presented. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

Wake flow pattern modified by small control cylinders at low Reynolds number

Passive wake control behind a circular cylinder in uniform flow is studied by numerical simulation for Re_D ranging from 80 to 300. Two small control cylinders, with diameter d/D=1/8, are placed at x/D=0.5 and y/D=±0.6. Unlike the 1990 results of Strykowski and Sreenivasan, in the present study, the vortex street behind the main cylinder still exists but the fluctuating lift and the form drag on the main cylinder reduces significantly and monotonously as the Reynolds number increases from 80 to 300. Obstruction of the control cylinders to the incoming flow deflects part of the fluid to pass through the gap between the main and control cylinders, forming two symmetric streams. These streams not only eliminate the flow separation along the rear surface of the main cylinder, they also merge toward the wake centerline to create an advancing momentum in the immediate near-wake region. These two effects significantly reduce the wake width behind the main cylinder and lead to monotonous decrease of the form drag as the Reynolds number increases. As the Reynolds number gets higher, a large amount of the downstream advancing momentum significantly delays the vortex formation farther downstream, leading to a more symmetric flow structure in the near-wake region of the main cylinder. As the Reynolds number increases from 80 to 300, both increasing symmetry of the flow structure in the near-wake and significant delay of the vortex formation are the main reasons for the fluctuating lift to decrease monotonously.     

Flow past a cylinder: shear layer instability and drag crisis

Flow past a circular cylinder for Re=100 to 10⁷ is studied numerically by solving the unsteady incompressible two‐dimensional Navier–Stokes equations via a stabilized finite element formulation. It is well known that beyond Re ～ 200 the flow develops significant three‐dimensional features. Therefore, two‐dimensional computations are expected to fall well short of predicting the flow accurately at high Re. It is fairly well accepted that the shear layer instability is primarily a two‐dimensional phenomenon. The frequency of the shear layer vortices, from the present computations, agree quite well with the Re^0.67 variation observed by other researchers from experimental measurements. The main objective of this paper is to investigate a possible relationship between the drag crisis (sudden loss of drag at Re ～ 2 × 10⁵) and the instability of the separated shear layer. As Re is increased the transition point of shear layer, beyond which it is unstable, moves upstream. At the critical Reynolds number the transition point is located very close to the point of flow separation. As a result, the shear layer eddies cause mixing of the flow in the boundary layer. This energizes the boundary layer and leads to its reattachment. The delay in flow separation is associated with narrowing of wake, increase in Reynolds shear stress near the shoulder of the cylinder and a significant reduction in the drag and base suction coefficients. The spatial and temporal power spectra for the kinetic energy of the Re=10⁶ flow are computed. As in two‐dimensional isotropic turbulence, E(k) varies as k^?5/3 for wavenumbers higher than energy injection scale and as k^?3 for lower wavenumbers. The present computations suggest that the shear layer vortices play a major role in the transition of boundary layer from laminar to turbulent state. Copyright © 2004 John Wiley & Sons, Ltd.     

Influence of the orientation of a square cylinder on the wake properties

An experimental investigation of flow around a square cylinder placed at various angles with respect to the approach fluid velocity is reported. The focus of the study is toward examining the sensitivity of the wake properties to the cylinder orientation and Reynolds number. Angles of incidence in the range of 0-60° and Reynolds numbers of 1340, 4990, and 9980 have been considered. Velocity measurements have been carried out using an X-wire hotwire anemometer. The Strouhal number and the drag coefficient of the cylinder have been computed from the wake measurements. Utilizing the velocity traces at distinct probe locations in the near and the far wake, statistical properties such as the RMS velocities and the spectra have been obtained. Results obtained in the present work revealed that for a cylinder with zero inclination, flow separates from the corners on the face exposed to the incoming flow. For inclinations greater than zero, the points of separation on the cylinder move downstream and the wake size increases, but the separated shear layer rolls up over a shorter distance. These factors lead to a reduced drag coefficient and a higher Strouhal number. The center-line recovery of the time-averaged velocity and the decay rates of velocity fluctuations depend on the Reynolds number. A marginal effect of the cylinder orientation is also seen.     

Non-unique turbulent boundary layer flows having a moderately large velocity defect: a rational extension of the classical asymptotic theory

The classical analysis of turbulent boundary layers in the limit of large Reynolds number Re is characterised by an asymptotically small velocity defect with respect to the external irrotational flow. As an extension of the classical theory, it is shown in the present work that the defect may become moderately large and, in the most general case, independent of Re but still remain small compared to the external streamwise velocity for non-zero pressure gradient boundary layers. That wake-type flow turns out to be characterised by large values of the Rotta–Clauser parameter, serving as an appropriate measure for the defect and hence as a second perturbation parameter besides Re. Most important, it is demonstrated that also this case can be addressed by rigorous asymptotic analysis, which is essentially independent of the choice of a specific Reynolds stress closure. As a salient result of this procedure, transition from the classical small defect to a pronounced wake flow is found to be accompanied by quasi-equilibrium flow, described by a distinguished limit that involves the wall shear stress. This situation is associated with double-valued solutions of the boundary layer equations and an unconventional weak Re-dependence of the external bulk flow—a phenomenon seen to agree well with previous semi-empirical studies and early experimental observations. Numerical computations of the boundary layer flow for various values of Re reproduce these analytical findings with satisfactory agreement.     

TR-PIV measurement of the wake behind a grooved cylinder at low Reynolds number

A comparative study of the wakes behind cylinders with grooved and smooth surfaces was performed with a view to understand the wake characteristics associated with the adult Saguaro cacti. A low-speed recirculation water channel was established for the experiment; the Reynolds number, based on the free-stream velocity and cylinder diameter (D), was kept at Re_D=1500. State-of-the-art time-resolved particle image velocimetry (TR-PIV) was employed to measure a total of 20 480 realizations of the wake field at a frame rate of 250 Hz, enabling a comprehensive view of the time- and phase-averaged wake pattern. In comparison to the wake behind the smooth cylinder, the length of the recirculation zone behind the grooved cylinder was extended by nearly 18.2%, yet the longitudinal velocity fluctuation intensity was considerably weakened. A global view of the peaked spectrum of the longitudinal velocity component revealed that the intermediate region for the grooved cylinder, which approximately corresponds to the transition region where the shear layer vortices interact, merge and shed before the formation of the Karman-like vortex street, was much wider than that for the smooth one. The unsteady events near St=0.3-0.4 were detected in the intermediate region behind the grooved cylinder, but no such events were found in the smooth cylinder system. Although the formation of the Karman-like vortex street was delayed by about 0.6D downstream for the grooved cylinder, no prominent difference in the vortex street region was found in the far wake for both cylinders. The Proper Orthogonal Decomposition (POD) method was used extensively to decompose the vector and swirling strength fields, which gave a close-up view of the vortices in the near wake. The first two POD modes of the swirling strength clarified the spatio-temporal characteristics of the shear layer vortices behind the grooved cylinder. The small-scale vortices superimposed on the shear layers behind the grooved cylinder were found to be generated and convected downstream in the same phase, which would significantly reduce the fluctuating force on the cylinder surface.     

Flow past a square cylinder at low Reynolds numbers

Results are presented for the flow past a stationary square cylinder at zero incidence for Reynolds number, Re ? 150. A stabilized finite‐element formulation is employed to discretize the equations of incompressible fluid flow in two‐dimensions. For the first time, values of the laminar separation Reynolds number, Re_s, and separation angle, θ_s, at Re_s are predicted. Also, the variation of θ_s with Re is presented. It is found that the steady separation initiates at Re = 1.15. Contrary to the popular belief that separation originates at the rear sharp corners, it is found to originate from the base point, i.e. θ_s=180^° at Re = Re_s. For Re > 5, θ_s approaches the limit of 135 ^°. The length of the separation bubble increases approximately linearly with increasing Re. The drag coefficient varies as Re^?0.66. Flow characteristics at Re ? 40 are also presented for elliptical cylinders of aspect ratios 0.2, 0.5, 0.8 and 1 (circle) having the same characteristic dimension as the square and major axis oriented normal to the free‐stream. Compared with a circular cylinder, the flow separates at a much lower Re from a square cylinder leading to the formation of a bigger wake (larger bubble length and width). Consequently, at a given Re, the drag on a square cylinder is more than the drag of a circular cylinder. This suggests that a cylinder with square section is more bluff than the one with circular section. Among all the cylinder shapes studied, the square cylinder with sharp corners generates the largest amount of drag. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical simulation of Reynolds number effects on velocity shear flow around a circular cylinder

Three-dimensional Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) are performed to investigate the shear effects on flow around a circular cylinder at Reynolds numbers of Re=60–1000. The shear parameter, β, which is based on the velocity gradient, cylinder diameter and upstream mean velocity at the center plane of the cylinder, varies from 0 to 0.30. Variations of Strouhal number, drag and lift coefficients, and unsteady wake structures with shear parameter are studied, along with their dependence on Reynolds number. The presented simulation provides detailed information for the flow field around a circular cylinder in shear flow. This study shows that the Strouhal number exhibits no significant variation with shear parameter. The stagnation point moves to the high-velocity side almost linearly with shear parameter, and this result mainly influences the aerodynamic forces acting on a circular cylinder in shear flow. Both the Reynolds number and shear parameter influence the movement of the stagnation point and separation point. Mode A wake instability is suppressed into parallel vortex shedding mode at a certain shear parameter. The lift force increases with increasing shear parameter and acts from the high-velocity side to the low-velocity side. In addition, a simple method to estimate the lift force coefficient in shear flow is provided.     

Role of the shear layer instability in the near wake behavior of two side-by-side circular cylinders

Wakes, and their interaction behind two parallel cylinders lying in a plane perpendicular to the flow, have been investigated experimentally in the sub-critical Reynolds number regime. The experiments were performed in a water channel using laser Doppler velocimetry. The gap between the two cylinders was less than the cylinder diameter, a geometry referred to as strong interaction configuration. In this case the blockage is strong and a gap-jet appears between the cylinders. Two flow regimes of the near wake region have been identified: one below a critical Reynolds number Re _c ]1000;1700[, where the gap jet is stably deflected to one side and the double near-wake becomes asymmetric; the other, above Re _c, where the gap-jet deflection is unstable and a random flopping phenomenon takes place. When Re<Re _c, two different Strouhal numbers are identified, related to the Kármán vortex shedding behind each cylinder. When Re>Re _c, a third frequency appears in the near wake, related to the development of Kelvin-Helmholtz vortices in the separated shear layer of the cylinders [Prasad A, Williamson CHK (1997) J Fluid Mech 333:375]. The observed flopping behavior is attributed to the birth of these Kelvin-Helmholtz instabilities and their intermittent nature. Further downstream, beyond about five cylinder diameters, the random flopping flow phenomena disappear while a slightly asymmetric single wake persists. It is characterized by a Strouhal number St=0.13, a value that one would normally measure behind a single cylinder of twice its diameter.