Flow-excited acoustic resonance of two side-by-side cylinders in cross-flow

The aeroacoustic response of two side-by-side circular cylinders in cross-flow is investigated experimentally. In order to investigate the effect of the gap between the cylinders on the acoustic resonance mechanism, six spacing ratios between the cylinders, in the range of T/D=1.25–3, have been investigated, where D is the diameter of the cylinders and T the centre-to-centre distance between them. Special attention is given to the intermediate spacing ratio range, which exhibits bistable flow regimes in the absence of resonance. During the tests, the acoustic cross-modes of the duct housing the cylinders are self-excited. For the intermediate spacing ratios, T/D=1.25, 1.35, 1.46 and 1.75, two distinct vortex-shedding frequencies at the off-resonance conditions are observed. These are associated with the wide and narrow wakes of the cylinders, as described in the literature. In this case, acoustic resonances occur at a Strouhal number, which is between those observed before the onset of resonance. The acoustic resonance synchronizes vortex shedding in the two wakes and thereby eliminates the bistable flow phenomenon. For large spacing ratios, T/D=2.5 and 3, vortex shedding occurs at a single Strouhal number at which the acoustic resonance is excited.     

Effect of acoustic resonance on the dynamic lift forces acting on two tandem cylinders in cross-flow

Direct measurements of the dynamic lift force acting on two tandem cylinders in cross-flow are performed in the presence and absence of acoustic resonance. The dynamic lift force is measured because it represents the integrated effect of the unsteady wake and therefore it is directly related to the dipole sound source generated by vortex shedding from the cylinder. Three spacing ratios inside the proximity interference region, L/D=1.75, 2.5 and 3 are considered. During the tests, the first transverse acoustic mode of the duct housing the cylinders is self-excited. In the absence of acoustic resonance, the measured dynamic lift coefficients agree with those reported in the literature. When the acoustic resonance is initiated, a drastic increase in the dynamic lift coefficient is observed, especially for the downstream cylinder. This can be associated with abrupt changes in the phase between the lift forces and the acoustic pressure. The dynamic lift forces on both cylinders are also decomposed into in-phase and out-of-phase components, with respect to the resonant sound pressure. The lift force components for the downstream cylinder are found to be dominant. Moreover, the out-of-phase component of the lift force on the downstream cylinder is found to become negative over two different ranges of flow velocity and to virtually vanish between these two ranges. Acoustic resonance of the first mode is therefore excited over two ranges of flow velocity separated by a non-resonant range near the velocity of frequency coincidence. It is therefore concluded that the occurrence of acoustic resonance is controlled by the out-of-phase lift component of the downstream cylinder, whereas the effect of the in-phase lift component is confined to causing small changes in the acoustic resonance frequency.     

Vortex shedding and acoustic resonance of single and tandem finned cylinders

The effect of fins on vortex shedding and acoustic resonance is investigated for isolated and two tandem cylinders exposed to cross-flow in a rectangular duct. Three spacing ratios between the tandem cylinders (S/D_e=1.5, 2 and 3) are tested for a Reynolds number range from 1.6×10⁴ to 1.1×10⁵. Measurements of sound pressure as well as mean and fluctuating velocities are performed for bare and finned cylinders with three different fin densities. The effect of fins on the sound pressure generated before the onset of acoustic resonance as well as during the pre-coincidence and coincidence resonance is found to be rather complex and depends on the spacing ratio between cylinders, the fin density and the nature of the flow-sound interaction mechanism.For isolated cylinders, the fins reduce the strength of vortex shedding only slightly, but strongly attenuate the radiated sound before and during the occurrence of acoustic resonance. This suggests that the influence of the fins on correlation length is stronger than on velocity fluctuations. In contrast to isolated cylinders, the fins in the tandem cylinder case enhance the vortex shedding process at off-resonant conditions, except for the large spacing case which exhibits a reversed effect at high Reynolds numbers. Regarding the acoustic resonance of the tandem cylinders, the fins promote the onset of the coincidence resonance, but increasing the fin density drastically weakens the intensity of this resonance. The fins are also found to suppress the pre-coincidence resonance for the tandem cylinders with small spacing ratios (S/D_e=1.5, 2 and 2), but for the largest spacing case (S/D_e=3), they are found to have minor effects on the sound pressure and the lock-in range of the pre-coincidence resonance.     

FLUID BEHAVIOUR OF SIDE-BY-SIDE CIRCULAR CYLINDERS IN STEADY CROSS-FLOW

The flow field for two and three circular cylinders of equal diameter D arranged in a side-by-side configuration in steady cross-flow was investigated using flow visualization, hot-film anemometry and particle image velocimetry (PIV), for centre-to-centre pitch ratios from T/D=1·0 to 6·0, and Reynolds numbers from Re=500 to 3000. For two-cylinder arrangements, three basic flow patterns were observed: single bluff-body vortex shedding at small T/D , biased flow with synchronized vortex shedding at intermediateT /D , and symmetric flow with synchronized vortex shedding at larger T/D . For three-cylinder arrangements, either single bluff-body behaviour or an asymmetric biased flow pattern could be observed at smallT /D , whereas a symmetric-biased flow pattern was found at intermediate T/D . Instantaneous representations of the in-plane vorticity field obtained from the PIV technique revealed some variation in these basic flow patterns at given T/D and Re.     

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.     

Characteristics and suppression of flow-induced vibrations of two side-by-side circular cylinders

A free-vibration experiment was conducted to examine flow-induced vibration (FIV) characteristics of two identical circular cylinders in side-by-side arrangements at spacing ratio T^⁎ (=T/D)=0.1–3.2, covering all possible flow regimes, where T is the gap spacing between the cylinders and D is the cylinder diameter. Each of the cylinders was two-dimensional, spring mounted, and allowed to vibrate independently in the cross-flow direction. Furthermore, an attempt to suppress flow-induced vibrations was undertaken by attaching flexible sheets at the rear stagnation lines of the cylinders. Based on the vibration responses of the two cylinders, four vibration patterns I, II, III and IV are identified at 0.1≤T^⁎<0.2, 0.2≤T^⁎≤0.9, 0.9<T^⁎<2.1 and 2.1≤T^⁎≤3.2, respectively. Pattern I is characterized by the two cylinders vibrating inphase^, with the maximum amplitudes occurring at the same reduced velocity Ur=10.47 almost two times that (Ur=5.25) for an isolated cylinder. Pattern II features no vibration generated for either cylinder. Pattern III exemplifies the occurrence of the maximum vibration amplitude of a cylinder at a smaller Ur and that of the other cylinder at a higher Ur compared to its counterpart in an isolated cylinder. Pattern IV represents each cylinder response resembling an isolated cylinder response; the vibrations of the two cylinders are, however, coupled inphase or antiphase. Linking maximum vibration amplitudes to fluctuating lift forces acting on fixed cylinders reveals that fluid–structure interactions between two fixed cylinders and between two elastic cylinders are not the same, even though vibration is not significant. As such, while two fixed cylinders generate narrow and wide wakes at 0.2≤T^⁎<1.7, two elastic cylinders do the same for a longer range of T^⁎ (0.2≤T^⁎<2.1). The flexible sheets effectively suppress FIV of the two cylinders in patterns III and IV, and reduce the vibration amplitude in pattern I. For the effectively controlled cases (patterns III and IV), the flexible sheet of each cylinder folds into a semicircle at the base, forming two free edges.     

Measurements of the dynamic lift force acting on a circular cylinder in cross-flow and exposed to acoustic resonance

A piezoelectric transducer is developed to perform direct measurements of the dynamic lift force acting on a circular cylinder in cross-flow, in the presence and absence of acoustic resonance. Details of the force transducer design are presented in the paper. The dynamic lift force is measured for a single cylinder with two different diameters, D=12.7 and 15.8 mm. During the tests, the first transverse acoustic mode of the duct housing the cylinder is self-excited. The fluctuating pressure on the top wall of the duct is measured simultaneously with the dynamic lift force. In the absence of acoustic resonance, the measured dynamic lift coefficients agree favorably with those reported in the literature. However, when the acoustic resonance is initiated, the dynamic lift experiences a drastic increase in amplitude associated with abrupt changes in the phase between the lift force and the acoustic pressure. A methodology to extract the hydrodynamic lift component from the total lift measured during acoustic resonance is also proposed. The hydrodynamic lift force is then decomposed into in-phase and out-of-phase components, with respect to the resonant sound pressure. This decomposition procedure provides new insights into the nature of the aeroacoustic sources in the cylinder wake. The proposed methodology, together with the test results provide a general design approach to assess the increase in the dynamic fluid loading on bluff bodies in cross-flow due to the excitation of acoustic resonance.     

基于声场精细积分算法的潜艇流激噪声预报

建立了Suboff潜艇流场、结构场和声场耦合模型,对其水下航行时流激引起的艇体低频结构振动及水下辐射噪声进行了研究。采用CFD方法对潜艇三维流场进行非定常计算,以此得到艇体表面的脉动力,并通过数据映射将该脉动力加载到潜艇有限元模型进行响应分析。以潜艇壳体响应结果为边界条件,采用声场精细积分算法对潜艇流激引起的水下辐射噪声进行预报,该算法可对Helmholtz边界积分方程(HBIE)中的弱奇异积分系数及近场积分系数进行自适应求解,当细化步数达到6时,积分以10-4的残差单调收敛。结果表明,流激作用下,潜艇水下辐射声功率及声辐射效率皆出现明显峰值,其中流激的线谱成份引起艇体受激强迫振动,辐射较强水声,宽带谱成份则引起艇体结构共振,可辐射更强的水声。因此,降低潜艇的流激噪声,需同时优化艇体型线及内部结构以改善伴流,降低湍流强度,重点是避免流激共振。     

Thermal study of an array of inline horizontal cylinders below a nearly adiabatic ceiling

Steady state two-dimensional free convection heat transfer from a horizontal, isothermal cylinder in a horizontal array of cylinders consists of three isothermal cylinders, located underneath a nearly adiabatic ceiling is studied experimentally. A Mach–Zehnder interferometer is used to determine thermal field and smoke test is made to visualize flow field. Effects of the cylinders spacing to its diameter (S/D), and cylinder distance from ceiling to its diameter (L/D) on heat transfer from the centered cylinder are investigated for Rayleigh numbers from 1500 to 6000. Experiments are performed for an inline array configuration of horizontal cylinders of diameters D = 13 mm. Results indicate that due to the nearly adiabatic ceiling and neighboring cylinders, thermal plume resulted from the centered cylinder separates from cylinder surface even for high L/D values and forming recirculation regions. By decreasing the space ratio S/D, the recirculation flow strength increases. Also, by decreasing S/D, boundary layers of neighboring cylinders combine and form a developing flow between cylinders. The strength of developing flow depends on the cylinders Rayleigh number and S/D ratio. Due to the developing flow between cylinders, the vortex flow on the top of the centered cylinder appears for all L/D ratios and this vortex influences the value of local Nusselt number distribution around the cylinder.Variation of average Nusselt number of the centered cylinder depends highly on L/D and the trend with S/D depends on the value of Rayleigh number.     

Vortex-shedding from single and tandem cylinders in the presence of applied sound

A single cylinder and two tandem cylinder configurations with longitudinal pitch ratios L/D=1.75 and 2.5 were rigidly mounted in an open circuit wind tunnel and a standing acoustic pressure wave was imposed so that the acoustic particle velocity was normal to both the cylinder axis and the mean flow velocity. The effect of sound on the vortex-shedding was investigated for various amplitudes by means of pressure taps on the cylinders and wake hot-wire probes. These tests show that applied sound can entrain and shift the natural vortex-shedding frequency to the frequency of excitation and produce nonlinearities in the wake. The lock-in envelope for the tandem cylinders is considerably larger than for the single cylinder. The lock-in range for the smaller tandem cylinder spacing was broader still than either the single cylinder, or the L/D=2.5 tandem cylinder case. The pressure and hot-wire measurements show for the single cylinder, and tandem cylinder configuration with pitch ratio L/D=2.5, that there was a phase jump near the coincidence of the vortex-shedding frequency and the excitation frequency, while there was no jump for the pitch ratio of 1.75. As well, the applied sound field was also noted to induce vortex-shedding in the gap for the L/D=2.5 case, while no vortex-shedding was noted for the smaller pitch ratio.     

Phenomenology of Kármán vortex streets in oscillatory flow

Vortex wakes of circular cylinders at low Reynolds numbers have been investigated. Sound waves are superimposed on the flow in mean flow direction. In this configuration the Kármán vortices are shed at the sound frequency or at subharmonics of the sound frequency. The Karman vortex street is treated as a nonlinear self-excited flow oscillator with forced oscillations. Using a flow visualization technique a variety of wake structures has been identified as a function of sound frequency and sound amplitude, but independent of the Reynolds number. The superimposed sound influences the distribution of circulation and accordingly the shedding mechanism. Primary vortex and secondary vortex are shed simultaneously from one side of the cylinder. The alternate vortex shedding is arranged spatially and temporally. Structures along the vortex axes are revealed.Parts of this paper have been presented at IUTAM Congress 1984, Lyngby, Denmark, and at ICNM Conference 1985, Shanghai, P. R. China     

Features of flow-induced forces deduced from wavelet analysis

In the present study, the effect of Reynolds number (Re) on flow interference between two side-by-side stationary cylinders and the associated flow-induced forces are investigated using finite element method and wavelet analysis. The pitch ratio chosen is T/D=1.7, where T is the separation distance measured between cylinder centers and D is the diameter, and Re, based on the free-stream velocity and the diameter of the cylinder, is varied within the laminar flow regime, i.e., 60<Re<200. The method of continuous wavelet transform is used to analyze time-variant features of flow-induced forces in the time–frequency domain. Flow patterns in the form of vorticity plots are presented to demonstrate the underlying physics. It is found that flow interference initially occurs in the inner vortices shed from the two cylinders, and extends to the outer vortices with increasing Re. The flow behind two cylinders undergoes three regimes: Regime I—unbiased gap flow, Regime II—stable biased gap flow, and Regime III—unstable gap flow. Flow-induced forces show significant variations when the flow transits from one regime to another. In particular, during the transition from Regimes II to III, the forces not only increase by amplitude, but also change their nature from deterministic to random, and show some nonstationary features. This is shown to be caused by the amalgamation of inner and outer vortices behind the two cylinders when the flow interference extends from inner vortices to outer vortices. Whenever possible, the present results are compared with experimental measurements and theoretical predictions. The numerical simulations are consistent with these other results.     

Large eddy simulations of the flow past two side-by-side circular cylinders

A LES Large Eddy Simulation is performed to study the flow past two side-by-side circular cylinders at a Reynolds number of 5800, based on the free-stream velocity and the cylinders diameter. The centre-to-centre transverse pitch ratio T/D is varied from 1.5 to 3. Both cylinders are slightly heated and the small amount of heat can be treated as a passive scalar. The numerical simulations are in good agreement with experimental observations.     

Global aerodynamic instability of twin cylinders in cross flow

This paper comprises an in-depth physical discussion of the flow-induced vibration of two circular cylinders in view of the time-mean lift force on stationary cylinders and interaction mechanisms. The gap-spacing ratio T/D is varied from 0.1 to 5 and the attack angle α from 0° to 180° where T is the gap width between the cylinders and D is the diameter of a cylinder. Mechanisms of interaction between two cylinders are discussed based on time-mean lift, fluctuating lift, flow structures and flow-induced responses. The whole regime is classified into seven interaction regimes, i.e., no interaction regime; boundary layer and cylinder interaction regime; shear-layer/wake and cylinder interaction regime; shear-layer and shear-layer interaction regime; vortex and cylinder interaction regime; vortex and shear-layer interaction regime; and vortex and vortex interaction regime. Though a single non-interfering circular cylinder does not correspond to a galloping following quasi-steady galloping theory, two circular cylinders experience violent galloping vibration due to shear-layer/wake and cylinder interaction as well as boundary layer and cylinder interaction. A larger magnitude of fluctuating lift communicates to a larger amplitude vortex excitation.     

Flow around four cylinders arranged in a square configuration

This paper presents an experimental study of the flow around four circular cylinders arranged in a square configuration. The Reynolds number was fixed at Re=8000, the pitch-to-diameter ratio between adjacent cylinders was varied from P/D=2 to 5 and the incidence angle was changed from α=0° (in-line square configuration) to 45° (diamond configuration) at an interval of 7.5°. The flow field was measured using digital Particle Image Velocimetry (PIV) to examine the vortex shedding characteristics of the cylinders, together with direct measurement of fluid dynamic forces (lift and drag) on each cylinder using a piezoelectric load cell. Depending on the pitch ratio, the flow could be broadly classified as shielding regime (P/D≤2), shear layer reattachment regime (2.5≤P/D≤3.5) and vortex impinging regime (P/D≥4). However, this classification is valid only in the case that the cylinder array is arranged nearly in-line with the free stream (α≈0°), because the flow is also sensitive to α. As α increases from 0° to 45°, each cylinder experiences a transition of vortex shedding pattern from a one-frequency mode to a two-frequency mode. The flow interference among the cylinders is complicated, which could be non-synchronous, quasi-periodic or synchronized with a definite phase relationship with other cylinders depending on the combined value of α and P/D. The change in vortex pattern is also reflected by some integral parameters of the flow such as force coefficients, power spectra and Strouhal numbers.     

Numerical simulation of cross-flow around four cylinders in an in-line square configuration

Successful numerical simulations can reveal important flow characteristics and information which are extremely difficult to obtain experimentally. Two- and three-dimensional (3-D) numerical simulations of cross-flow around four cylinders in an in-line square configuration are performed using a finite-volume method. For 2-D studies, the Reynolds numbers (Re) are chosen to be Re=100 and 200 and the spacing ratio L/D is set at 1.6, 2.5, 3.5, 4.0 and 5.0. For the 3-D investigation, the simulation is only performed at a Re=200, a spacing ratio L/D=4.0 and an aspect ratio H/D=16. The 2-D studies reveal three distinct flow patterns: (I) a stable shielding flow; (II) a wiggling shielding flow and (III) a vortex shedding flow. A transformation of the flow pattern from (I) to (II) at Re=100 will increase the amplitude of the maximum fluctuating pressure on the downstream cylinder surface by 4–12 times, while a transformation of the flow pattern from (II) to (III) will enhance the maximum fluctuating pressure amplitude by 2–3 times. There is a large discrepancy between 2-D simulation and flow visualization results at L/D=4.0 and Re=200. A probable cause could be the strong 3-D effect at the ends of the cylinder at low H/D. It was found that, for an in-line square configuration at critical L/D and when H/D is lower than a certain value, 3-D effects are very significant at the ends of the cylinders. In such cases, a time-consuming 3-D numerical simulation will have to be performed if full replication of the flow phenomenon were to be achieved.     

The role of shape and relative submergence on the structure of wakes of low-aspect-ratio wall-mounted bodies

The effects of shape and relative submergence (the ratio of flow depth to obstacle height, d/H) are investigated on the wakes around four different low-aspect-ratio wall-mounted obstacles at Re _H  = 17,800: semi-ellipsoids with the major axes of the base ellipses aligned in the streamwise and transverse directions, and two cylinders with aspect ratios matching the ellipsoids (H/D = 0.89 and 0.67, where D is the maximum transverse dimension). Particle Image Velocimetry was used to interrogate the flow. Streamwise features observed in the mean wake include counter-rotating distributions of vorticity inducing downwash (tip structures), upwash (base structures), and horseshoe vortices. In particular, the relatively subtle change in geometry produced by the rotation of the ellipsoid from the streamwise to the transverse orientation results in a striking modification of the mean streamwise vorticity distribution in the wake. Tip structures are dominant in the former case, while base structures are dominant in the latter. A vortex skeleton model of the wake is proposed in which arch vortex structures, shed from the obstacle, are deformed by the competing mechanisms of Biot-Savart self-induction and the external shear flow. The selection of tip or base structures in the ellipsoid wakes is caused by tilting of the arch structures either upstream or downstream, respectively, which is governed by ellipsoid curvature. An inverse relationship was observed between the relative submergence and the strength of the base structures for the ellipsoids, with a dominant base structure observed for d/H = 1 in both cases. These results demonstrate a means by which to achieve significant modifications to flow structure and thereby also to transport mechanisms in the flow. Therefore, this work provides insight into the modeling and control of flow over wall-mounted bodies.     

Sources of acoustic resonance generated by flow around a long rectangular plate in a duct

This paper describes a numerical investigation of the acoustic resonance occurring at specific flow speeds when a long two-dimensional rectangular plate is placed on the centre-line of a duct. The flow Mach number is sufficiently small that the flow and acoustic fields can be modelled separately; however, the effect of the acoustic field in modifying the flow field is accounted for and, in turn, the flow field solution determines the time-dependent source distribution for the acoustic model. This allows the range of flow speeds, or equivalently the associated Strouhal numbers, where resonance is possible to be predicted. It is shown that the Strouhal number based on plate chord (or length) displays a stepping behaviour as the plate length is increased. The main source or sink region where energy is transferred between the acoustic and flow fields is shown to be immediately downstream of the trailing edge of the plate. Visualizations of the numerical solution show that the timing when the vortices enter this region relative to the phase of the acoustic cycle is crucial in determining if resonance can occur and is the cause of the observed stepwise increase. Comparison is made with previous physical experiments.     

Some useful hybrid approaches for predicting aerodynamic noise

In recent years, several numerical studies have shown the feasibility of Direct Noise Computation (DNC) where the turbulent flow and the radiated acoustic field are obtained simultaneously by solving the compressible Navier–Stokes equations. The acoustic radiation obtained by DNC can be used as reference solution to investigate hybrid methods in which the sound field is usually calculated as a by-product of the flow field obtained by a more conventional Navier–Stokes solver. A hybrid approach is indeed of practical interest when only the non-acoustic part of the aerodynamic field is available. In this review, some acoustic analogies or hybrid approaches are revisited in the light of CAA. To cite this article: C. Bailly et al., C. R. Mecanique 333 (2005).     

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.