scispace - formally typeset
Search or ask a question

Showing papers on "Vortex shedding published in 2012"


Journal ArticleDOI
Yan Bao1, Cheng Huang1, Dai Zhou1, Jiahuang Tu1, Zhaolong Han1 
TL;DR: In this article, the effect of the natural frequency ratio on the characteristics of vortex-induced vibration (VIV) responses, including wake frequencies, orbital trajectories, response amplitudes, hydrodynamic forces and wake mode patterns, was investigated.

159 citations


Journal ArticleDOI
TL;DR: In this article, the influence of angle of attack of a square section cylinder on the cylinder's flow-induced vibration was examined, where the direction of the vibration is transverse to the oncoming flow.
Abstract: This study examines the influence of angle of attack of a square section cylinder on the cylinder’s flow-induced vibration, where the direction of the vibration is transverse to the oncoming flow. Our experiments, which traversed the velocity–angle of attack parameter space in considerable breadth and depth, show that a low-mass ratio body can undergo combinations of both vortex-induced vibration and galloping. When the body has an angle of attack that makes it symmetric to the flow, such as when it assumes the square or diamond orientation, the two mechanisms remain independent. However, when symmetry is lost we find a mixed mode response with a new branch of vortex-induced oscillations that exceeds the amplitudes resulting from the two phenomena independently. The oscillations of this higher branch have amplitudes larger than the ‘upper branch’ of vortex-induced vibrations and at half the frequency. For velocities above this resonant region, the frequency splits into two diverging branches. Analysis of the amplitude response reveals that the transition between galloping and vortex-induced vibrations occurs over a narrow range of angle of incidence. Despite the rich set of states found in the parameter space the vortex shedding modes remain very similar to those found previously in vortex-induced vibration.

156 citations


Journal ArticleDOI
TL;DR: In this paper, the flow over a circular cylinder at Reynolds number 3900 and Mach number 0.2 was predicted numerically using the technique of large-eddy simulation using an O-type curvilinear grid of size of 300 × 300 × 64.
Abstract: The flow over a circular cylinder at Reynolds number 3900 and Mach number 0.2 was predicted numerically using the technique of large-eddy simulation. The computations were carried out with an O-type curvilinear grid of size of 300 × 300 × 64. The numerical simulations were performed using a second-order finite-volume method with central-difference schemes for the approximation of convective terms. A conventional Smagorinsky and a dynamic k-equation eddy viscosity sub-grid scale models were applied. The integration time interval for data sampling was extended up to 150 vortex shedding periods for the purpose of obtaining a fully converged mean flow field. The present numerical results were found to be in good agreement with existing experimental data and previously obtained large-eddy simulation results. This gives an indication on the adequacy and accuracy of the selected large-eddy simulation technique implemented in the OpenFOAM toolbox.

150 citations


Journal ArticleDOI
TL;DR: Ultrahigh-speed video imaging is used in combination with high-resolution numerical simulations to show how this ejecta gives way to irregular splashing, and at higher Reynolds numbers, its base becomes unstable, shedding vortex rings into the liquid from the free surface in an axisymmetric von Kármán vortex street.
Abstract: The splashing of a drop impacting onto a liquid pool produces a range of different sized microdroplets. At high impact velocities, the most significant source of these droplets is a thin liquid jet emerging at the start of the impact from the neck that connects the drop to the pool. We use ultrahigh-speed video imaging in combination with high-resolution numerical simulations to show how this ejecta gives way to irregular splashing. At higher Reynolds numbers, its base becomes unstable, shedding vortex rings into the liquid from the free surface in an axisymmetric von Karman vortex street, thus breaking the ejecta sheet as it forms.

138 citations


Journal ArticleDOI
TL;DR: In this paper, a cantilevered piezoelectric beam is excited in a heating, ventilation and air conditioning (HVAC) flow by the interactions between an aerodynamic fin attached at the end of the piezolectric cantilever and the vortex shedding downstream from a bluff body placed in the air flow ahead of the fin/cantilever assembly.
Abstract: A cantilevered piezoelectric beam is excited in a heating, ventilation and air conditioning (HVAC) flow. This excitation is amplified by the interactions between (a) an aerodynamic fin attached at the end of the piezoelectric cantilever and (b) the vortex shedding downstream from a bluff body placed in the air flow ahead of the fin/cantilever assembly. The positioning of small weights along the fin enables tuning of the energy harvester to operate at resonance for flow velocities from 2 to 5 m s−1, which are characteristic of HVAC ducts. In a 15 cm diameter air duct, power generation of 200 μW for a flow speed of 2.5 m s−1 and power generation of 3 mW for a flow speed of 5 m s−1 was achieved. These power outputs are sufficient to power a wireless sensor node for HVAC monitoring systems or other sensors for smart building technology.

131 citations


Journal ArticleDOI
TL;DR: It is reported that a distributed pressure sensing mechanism has the capability to discriminate KáRMán vortex streets from uniform flows, and determine the orientation and position of the platform with respect to the incoming flow and the centre axis of the Kármán vortex street.
Abstract: With the overall goal being a better understanding of the sensing environment from the local perspective of a situated agent, we studied uniform flows and Karman vortex streets in a frame of reference relevant to a fish or swimming robot. We visualized each flow regime with digital particle image velocimetry and then took local measurements using a rigid body with laterally distributed parallel pressure sensor arrays. Time and frequency domain methods were used to characterize hydrodynamically relevant scenarios in steady and unsteady flows for control applications. Here we report that a distributed pressure sensing mechanism has the capability to discriminate Karman vortex streets from uniform flows, and determine the orientation and position of the platform with respect to the incoming flow and the centre axis of the Karman vortex street. It also enables the computation of hydrodynamic features which may be relevant for a robot while interacting with the flow, such as vortex shedding frequency, vortex travelling speed and downstream distance between vortices. A Karman vortex street was distinguished in this study from uniform flows by analysing the magnitude of fluctuations present in the sensor measurements and the number of sensors detecting the same dominant frequency. In the Karman vortex street the turbulence intensity was 30% higher than that in the uniform flow and the sensors collectively sensed the vortex shedding frequency as the dominant frequency. The position and orientation of the sensor platform were determined via a comparative analysis between laterally distributed sensor arrays; the vortex travelling speed was estimated via a cross-correlation analysis among the sensors.

117 citations


Journal ArticleDOI
TL;DR: In this article, the authors studied the forces experienced by a falling sphere in both cavity-forming and non-cavity-forming cases during free-surface water entry of spheres of varying masses, diameters, and surface treatments.
Abstract: We present a study of the forces during free-surface water entry of spheres of varying masses, diameters, and surface treatments. Previous studies have shown that the formation of a subsurface air cavity by a falling sphere is conditional upon impact speed and surface treatment. This study focuses on the forces experienced by the sphere in both cavity-forming and non-cavity-forming cases. Unsteady force estimates require accurate determination of the deceleration for both high and low mass ratios, especially as inertial and hydrodynamic effects approach equality. Using high-speed imaging, high-speed particle image velocimetry, and numerical simulation, we examine the nature of the forces in each case. The effect of mass ratio is shown, where a lighter sphere undergoes larger decelerations and more dramatic trajectory changes. In the non-cavity-forming cases, the forces are modulated by the growth and shedding of a strong, ring-like vortex structure. In the cavity-forming cases, little vorticity is shed by the sphere, and the forces are modulated by the unsteady pressure required for the opening and closing of the air cavity. A data-driven boundary-element-type method is developed to accurately describe the unsteady forces using cavity shape data from experiments.

112 citations


Journal ArticleDOI
01 Sep 2012
TL;DR: In this paper, the effect of Reynolds number, surface roughness, freestream turbulence, proximity and wake interference on the radiated noise was studied on single and multiple rod configurations.
Abstract: Acoustic measurements were performed on single and multiple rod configurations to study the effect of Reynolds number, surface roughness, freestream turbulence, proximity and wake interference on the radiated noise. The Reynolds number ranged from 3.8 × 103 to 105. Directivity measurements were performed to determine how well the dipole assumption for the radiation of vortex shedding noise holds for the different model configurations tested. The dependence of the peak Sound Pressure Level on velocity was also examined. Several concepts for the reduction of the noise radiating from cylindrical rods were tested. It was shown that wire wraps and collar distributions could be used to significantly reduce the noise radiating from rods in tandem configurations.

111 citations


Journal ArticleDOI
TL;DR: In this article, the authors studied flexural vibrations of a cantilever beam with thin rectangular cross section submerged in a quiescent viscous fluid and undergoing oscillations whose amplitude is comparable with its width.

106 citations


Journal ArticleDOI
TL;DR: In this article, experimental results on the unsteady behaviour of the wake of a modelled wind turbine in an atmospheric boundary layer (ABL) wind tunnel are presented. But the results are limited to the case of a single turbine.

100 citations


Journal ArticleDOI
TL;DR: In this paper, a wind tunnel experiment on the flow around a circular cylinder with diameter D (40mm) attached with ten splitter plates freely rotatable around the cylinder axis has been carried out with different ratios of length to cylinder diameter (L / D ) from 0.5 to 6.0, in a range of Reynolds number from 3×10 4 to 6×104.

Journal ArticleDOI
TL;DR: In this article, a uniform flow vibrates a linear array of four cylinders affixed to piezoelectric energy transducers, and the energy generated by downstream cylinders can exceed that of leading cylinders by more than an order of magnitude.

Journal ArticleDOI
TL;DR: In this article, a computational fluid dynamics simulation was performed for a small-scale, high solidity H-type Darrieus vertical axis wind turbine, which has a large stationary domain and smaller rotating subdomain connected by a sliding mesh interface.
Abstract: A computational fluid dynamics simulation was performed for a small-scale, high solidity (σ = 0.48) H-type Darrieus vertical axis wind turbine. Two-dimensional unsteady Reynolds-averaged Navier–Stokes equations were solved for the turbine numerical model, which has a large stationary domain and smaller rotating subdomain connected by a sliding mesh interface. The simulation results were first validated against steady-state airfoil data. The model was then used to solve for three rotating blades with constant ambient flow velocity (Re = 360,000) over numerous blade speed ratios. The high solidity and the associated low blade speed ratio and rotational speed of the turbine result in complex flow–blade interaction mechanisms. These include dynamic stall resulting in vortex shedding, vortex impingement on the source blade and significant flow momentum extraction causing reduced power production from the downstream blade pass. Copyright © 2011 John Wiley & Sons, Ltd.

Journal ArticleDOI
TL;DR: In this paper, the authors investigated high Reynolds-number flow over a rotating cylinder by two-dimensional numerical computations, and the results revealed stabilization of the acting forces at high spin rates, thus indicating a flowfield with suppressed vortex shedding activity, as it is expected, in accordance with theoretical considerations in previous studies.

Journal ArticleDOI
TL;DR: In this article, a new electromagnetic energy harvester for harnessing energy from vibration induced by Karman vortex street is proposed, which converts flow energy into electrical energy by fluid flow, vortex shedding from a bluff body and electromagnetic induction.

Journal ArticleDOI
TL;DR: In this article, a numerical study of the two-dimensional and three-dimensional unsteady flow over two square cylinders arranged in an in-line configuration for Reynolds numbers from 40 to 1000 and a gap spacing of 4D, where D is the cross-sectional dimension of the cylinders.
Abstract: SUMMARY This paper describes a numerical study of the two-dimensional and three-dimensional unsteady flow over two square cylinders arranged in an in-line configuration for Reynolds numbers from 40 to 1000 and a gap spacing of 4D, where D is the cross-sectional dimension of the cylinders. The effect of the cylinder spacing, in the range G = 0.3D to 12D, was also studied for selected Reynolds numbers, that is, Re = 130, 150 and 500. An incompressible finite volume code with a collocated grid arrangement was employed to carry out the flow simulations. Instantaneous and time-averaged and spanwise-averaged vorticity, pressure, and streamlines are computed and compared for different Reynolds numbers and gap spacings. The time averaged global quantities such as the Strouhal number, the mean and the RMS values of the drag force, the base suction pressure, the lift force and the pressure coefficient are also calculated and compared with the results of a single cylinder. Three major regimes are distinguished according to the normalized gap spacing between cylinders, that is, the single slender-body regime (G < 0.5), the reattach regime (G < 4) and co-shedding or binary vortex regime (G ≥4). Hysteresis with different vortex patterns is observed in a certain range of the gap spacings and also for the onset of the vortex shedding. Copyright © 2011 John Wiley & Sons, Ltd.

Journal ArticleDOI
TL;DR: In this paper, a second-order characteristic-based split finite element (SDFE) algorithm was used to solve the Navier-Stokes equations for an inline cylinder array consisting of six square cylinders at a Reynolds number of 100.

Journal ArticleDOI
TL;DR: In this article, the authors investigated the feasibility of using open-cell metal foam covering layer to reduce the low Mach number aerodynamic noise generated by the flow around a circular cylinder which is the typical section of pantographs.

Journal ArticleDOI
TL;DR: In this paper, the root-to-tip vortex formation of a trapezoidal flat-plate fin undergoing rotation from rest at a 90° angle of attack and Reynolds numbers of O(103) is investigated.
Abstract: We investigate experimentally the unsteady, three-dimensional vortex formation of low-aspect-ratio, trapezoidal flat-plate fins undergoing rotation from rest at a 90° angle of attack and Reynolds numbers of O(103). The objectives are to characterize the unsteady three-dimensional vortex structure, examine vortex saturation, and understand the effects of the root-to-tip flow for different velocity programs. The experiments are conducted in a water tank facility, and the diagnostic tools are dye flow visualization and digital particle image velocimetry. The dye visualizations show that the low-aspect-ratio plate produces symmetric ring-like vortices comprised mainly of tip-edge vorticity. They also indicate the presence of the root-to-tip velocity. For large rotational amplitudes, the primary ring-like vortex sheds and a secondary ring-like vortex is generated while the plate is still in motion, indicating saturation of the leading vortex. The time-varying vortex circulation in the flow symmetry plane provides quantitative evidence of vortex saturation. The phenomenon of saturation is observed for several plate velocity programs. The temporal development of the vortex circulation is often complex, which prevents an objective determination of an exact saturation time. This is the result of an interaction between the developing vortex and the root-to-tip flow, which breaks apart the vortex. However, it is possible to define a range of time during which the vortex reaches saturation. A formation-parameter definition is investigated and is found to reasonably predict the state corresponding to the pinch-off of the initial tip vortex across the velocity programs tested. This event is the lower bound on the saturation time range.

Journal ArticleDOI
TL;DR: In this article, the influence of three-dimensional effects in hovering flapping flight was analyzed using numerical simulations at a Reynolds number of 1000 to compare two types of flapping kinematics whose plunging phase is characterized by either a rectilinear translation or a revolving motion.
Abstract: This paper aims at understanding the influence of three-dimensional effects in hovering flapping flight. Numerical simulations at a Reynolds number of 1000 are performed to compare two types of flapping kinematics whose plunging phase is characterized by either a rectilinear translation or a revolving motion. In this way, we are able to isolate the three-dimensional effects induced by the free end condition from that induced by the spanwise incident velocity gradient (and the associated implicit Coriolis and centrifugal effects). In the rectilinear translation case, the analysis of the wake and of the aerodynamic loads reveals that the wingspan can be compartmented into three distinct regions whether it is predominantly subjected to an unstable two-dimensional flow, a stable three-dimensional flow or both two-dimensional and three-dimensional effects. It is found that this partitioning exhibits common features for three different aspect ratios of the wing. In conjunction with the previous results of Ringuette, Milano & Gharib (J. Fluid Mech., vol. 581, 2007, pp. 453–468), this suggests that the influence of the tip vortex over the wingspan is driven by a characteristic length scale. In addition, this length scale matches the position of the connecting point between leading and tip vortices observed in the revolving case, providing insight into the connecting process. In both translating and revolving cases, leading edge vortex attachment and strong spanwise velocities are found to be strongly correlated phenomena. Spanwise velocities (that mostly confine at the periphery of the vortices), together with downward velocities, do not only affect the leading edge vortex but also act as an inhibitor for the trailing edge vortex growth. As a consequence, cross-wake interactions between leading and trailing edge vortices are locally limited, hence contributing to flow stabilization.

Journal ArticleDOI
TL;DR: In this article, the effects of gap-width on the aerodynamic characteristics of a twin-deck bridge in a nominally smooth flow were investigated, with a specific goal of studying the vortex shedding mechanisms of the bridge deck.

Journal ArticleDOI
TL;DR: In this article, the shedding process in the near wake of a surface-mounted, square cross-section cylinder of height-to-width aspect ratio 4 at a Reynolds number of 12,000 based on free-stream velocity and the obstacle width was investigated.
Abstract: The shedding process in the near wake of a surface-mounted, square cross-section cylinder of height-to-width aspect ratio 4 at a Reynolds number of 12,000 based on free-stream velocity and the obstacle width was investigated. The boundary layer thickness was 0.18 obstacle heights based on 99% free-stream velocity. The study is performed using planar high frame-rate particle image velocimetry synchronized with pressure measurements and hot-wire anemometry. Spatial cross-correlation, instantaneous phase relationships, and phase-averaged velocity data are reported. Two dominant vortex-shedding regimes are observed. During intervals of high-amplitude pressure fluctuations on the obstacle side faces, alternate formation and shedding of vortices is observed (regime A) similar to the von Karman process. Regime B is characterized by two co-existing vortices in the obstacle lee throughout the shedding cycle and is observed within low-amplitude pressure fluctuation intervals. Despite the coexisting vortices in the base region, opposite sign vorticity is still shed out-of-phase downstream of this vortex pair giving rise to a staggered arrangement of counter-rotating vortices downstream. While the probability of occurrence of Regime B increases toward the free end, the amplitude modulation remains coherent along the obstacle height. Conditionally phase-averaged reconstructions of the flow field are consistent with the spatial distribution of the phase relationships and their probability density function. Earlier observations are reconciled showing that the symmetric shedding of vortices is a rare occurrence.


Journal ArticleDOI
TL;DR: In this paper, a wind tunnel study was carried out with the objective of exploring the effect of shape on vortex induced vibrations of trapezoidal box girder bridge decks, and it was demonstrated that it is possible to derive a virtually vibration free deck and that the angle between the horizontal bottom plate and the inclined side panes is an important parameter in achieving this goal.

Journal ArticleDOI
TL;DR: In this article, an experimental analysis of path and shape oscillations of an air bubble of diameter d rising in water at high Reynolds number in a vertical Hele-Shaw cell of width h is presented.
Abstract: We report an experimental analysis of path and shape oscillations of an air bubble of diameter d rising in water at high Reynolds number in a vertical Hele-Shaw cell of width h. Liquid velocity perturbations induced by the relative movement have also been investigated to analyze the coupling between the bubble motion and the wake dynamics. The confinement ratio h/d is lower than unity so that the bubble is flattened in between the walls of the cell. As the bubble diameter is increased, the Archimedes and the Bond numbers increase within 10 6 Ar 6 104 and 6 × 10−3 6 Bo 6 140. Mean shapes become more and more elongated. They first evolve from in-plane circles to ellipses, then to complicated shapes without fore-aft symmetry and finally to semi-circular capped bubbles. The scaling law Re = 0.5Ar is however valid for a large range of Ar, indicating that the liquid films between the bubble and the walls do no contribute significantly to the drag force exerted on the bubble. The coupling between wake dynamics, bubble path and shape oscillations evolves and a succession of contrasted regimes of oscillations is observed. The rectilinear bubble motion becomes unstable from a critical value Ar1 through an Hopf bifurcation while the bubble shape is still circular. The amplitude of path oscillations first grows as Ar increases above Ar1 but then surprisingly decreases beyond a second Archimedes number Ar2. This phenomenon, observed for steady ellipsoidal shape with moderate eccentricity, can be explained by the rapid attenuation of bubble wakes caused by the confinement. Shape oscillations around a significantly elongated mean shape starts for Ar > Ar3. The wake structure progressively evolves due to changes in the bubble shape. After the break-up of the fore-aft symmetry, a fourth regime involving complicated shape oscillations is then observed for Ar > Ar4. Vortex shedding disappears and unsteady attached vortices coupled to shape oscillations trigger path oscillations of moderate amplitude. Path and shape oscillations finally decrease when Ar is further increased. For Ar > Ar5, capped bubbles followed by a steady wake rise on a straight path.

Journal ArticleDOI
TL;DR: In this article, large-eddy simulations of a coherent counter-rotating vortex pair in different environments are performed and it is revealed that the vortex bursting phenomenon, known from photos of aircraft contrails or smoke visualization, is caused by collisions of secondary vortical structures traveling along the vortex tube which expel material from the vortex but do not result in a sudden decay of circulation or an abrupt change of vortex core structure.
Abstract: Large-eddy simulations of a coherent counter-rotating vortex pair in different environments are performed. The environmental background is characterized by varying turbulence intensities and stable temperature stratifications. Turbulent exchange processes between the vortices, the vortex oval, and the environment, as well as the material redistribution processes along the vortex tubes are investigated employing passive tracers that are superimposed to the initial vortex flow field. It is revealed that the vortex bursting phenomenon, known from photos of aircraft contrails or smoke visualization, is caused by collisions of secondary vortical structures traveling along the vortex tube which expel material from the vortex but do not result in a sudden decay of circulation or an abrupt change of vortex core structure. In neutrally stratified and weakly turbulent conditions, vortex reconnection triggers traveling helical vorticity structures which is followed by their collision. A long-lived vortex ring links once again establishing stable double rings. Key phenomena observed in the simulations are supported by photographs of contrails. The vertical and lateral extents of the detrained passive tracer strongly depend on environmental conditions where the sensitivity of detrainment rates on initial tracer distributions appears to be low.

Journal ArticleDOI
TL;DR: In this article, the authors presented the results of numerical investigations of bridge aeroelasticity, in particular static coefficients and instationary flutter derivatives for a symmetric bridge deck section using the Unsteady Reynolds-Averaged Navier-Stokes (URANS) method.

Journal ArticleDOI
TL;DR: In this paper, the authors investigated the two-degree-of-freedom Vortex-Induced Vibration (VIV) of four rigidly coupled circular cylinders in a square configuration and showed that the approaching angle has significant effects on the response of the four-cylinder system.

Journal ArticleDOI
TL;DR: In this paper, the authors investigate global instability and vortex shedding in the separated laminar boundary layer beneath internal solitary waves (ISWs) of depression in a two-layer stratified fluid by performing high-resolution two-dimensional direct numerical simulations.
Abstract: We investigate global instability and vortex shedding in the separated laminar boundary layer beneath internal solitary waves (ISWs) of depression in a two-layer stratified fluid by performing high-resolution two-dimensional direct numerical simulations. The simulations were conducted with waves propagating over a flat bottom and shoaling over relatively mild and steep slopes. Over a flat bottom, the potential for vortex shedding is shown to be directly dependent on wave amplitude, for a particular stratification, owing to increase of the adverse pressure gradient ( for leftward propagating waves) beneath the trailing edge of larger amplitude waves. The generated eddies can ascend from the bottom boundary to as high as 33 % of the total depth in two-dimensional simulations. Over sloping boundaries, global instability occurs beneath all waves as they steepen. For the slopes considered, vortex shedding begins before wave breaking and the vortices, shed from the bottom boundary, can reach the pycnocline, modifying the wave breaking mechanism. Combining the results over flat and sloping boundaries, a unified criterion for vortex shedding in arbitrary two-layer continuous stratifications is proposed, which depends on the momentum-thickness Reynolds number and the non-dimensionalized ISW-induced pressure gradient at the point of separation. The criterion is generalized to a form that may be readily computed from field data and compared to published laboratory experiments and field observations. During vortex shedding events, the bed shear stress, vertical velocity and near-bed Reynolds stress were elevated, in agreement with laboratory observations during re-suspension events, indicating that boundary layer instability is an important mechanism leading to sediment re-suspension.

Journal ArticleDOI
TL;DR: In this paper, a flow regime map is proposed as a function of spacing and Reynolds number for six in-line square cylinders, where is the surface-to-surface distance between two cylinders, is the size of the cylinder and is the Reynolds number.
Abstract: The flow around six in-line square cylinders has been studied numerically and experimentally for and , where is the surface-to-surface distance between two cylinders, is the size of the cylinder and is the Reynolds number The effect of spacing on the flow regimes is initially studied numerically at for which a synchronous flow regime is observed for , while quasi-periodic-I, quasi-periodic-II and chaotic regimes occur between , and , respectively These regimes have been confirmed via particle-image-velocimetry-based experiments A flow regime map is proposed as a function of spacing and Reynolds number The flow is predominantly quasi-periodic-II or chaotic at higher Reynolds numbers The quasi-periodic and chaotic nature of the flow is due to the wake interference effect of the upstream cylinders which becomes more severe at higher Reynolds numbers The appearance of flow regimes is opposite to that for a row of cylinders The Strouhal number for vortex shedding is the same for all the cylinders, especially for synchronous and quasi-periodic-I flow regimes The mean drag () experienced by the cylinders is less than that for an isolated cylinder, irrespective of the spacing The first cylinder is relatively insensitive to the presence of downstream cylinders and the is almost constant at 12 The for the second and third cylinders may be negative, with the value of increasing monotonically with spacing The changes in root mean square lift coefficient are consistent with changes in Interestingly, the instantaneous lift force can be larger than the instantaneous drag force on the cylinders These results should help improve understanding of flow around multiple bluff bodies