Vortex Induced Vibrations of a Pair of Cylinders at Reynolds Number 1000

Vortex induced vibrations of two equal-sized cylinders in tandem and staggered arrangement placed in uniform incompressible flow is studied. A stabilized finite element formulation is utilized to solve the governing equations. The Reynolds number for these 2D simulations is 1000. The cylinders are separated by 5.5 times the cylinder diameter in the streamwise direction. For the staggered arrangement, the cross-flow spacing between the two cylinders is 0.7 times the cylinder diameter. In this arrangement, the downstream cylinder lies in the wake of the upstream one and therefore experiences an unsteady inflow. The wake looses its temporal periodicity, beyond a few diameters downstream of the front cylinder. The upstream cylinder responds as an isolated single cylinder while the downstream one undergoes disorganized motion. Soft-lock-in is observed in almost all the cases.     

串列双圆柱绕流问题的数值模拟

本文运用有限体积方法,对绕串列放置的双圆柱的二维不可压缩流动进行了数值计算。为研究两圆柱不同间距对圆柱相互作用和尾流特征的影响,选取间距比Ｌ／Ｄ（Ｌ为两圆柱中心间的距离,Ｄ为圆柱直径）在１．５～５．０之间每隔０．５共八个有代表性的间距进行了计算模拟。计算均在Ｒｅ＝２００条件下进行。计算结果表明：对该绕流问题,流动特征在很大程度上取决于间距的大小。且间距存在一临界值,间距比从小于临界值变化到大于临界     

Wake development in staggered short cylinder arrays within a channel

Staggered arrays of short cylinders, known as pin?Cfins, are commonly used as a heat exchange method in many applications such as cooling electronic equipment and cooling the trailing edge of gas turbine airfoils. This study investigates the near wake flow as it develops through arrays of staggered pin fins. The height-to-diameter ratio was unity while the transverse spacing was kept constant at two cylinder diameters. The streamwise spacing was varied between 3.46 and 1.73 cylinder diameters. For each geometric arrangement, experiments were conducted at Reynolds numbers of 3.0e3 and 2.0e4 based on cylinder diameter and velocity through the minimum flow area of the array. Time-resolved flowfield measurements provided insight into the dependence of row position, Reynolds number, and streamwise spacing. Decreasing streamwise spacing resulted in increased Strouhal number as the near wake length scales were confined. In the first row of the bundle, low Reynolds number flows were mainly shear-layer-driven while high Reynolds number flows were dominated by periodic vortex shedding. The level of velocity fluctuations increased for cases having stronger vortex shedding. The effect of streamwise spacing was most apparent in the reduction of velocity fluctuations in the wake when the spacing between rows was reduced from 2.60 diameters to 2.16 diameters.     

THE FLOW INTERACTION BETWEEN A STATIONARY CYLINDER AND A DOWNSTREAM FLEXIBLE CYLINDER

This paper deals with an experimental investigation of the dynamic response of a long flexible circular cylinder in the wake of a stationary geometrically similar cylinder. The system has a low damping ratio. The tandem and staggered arrangements have been considered. The separations investigated are typical of transmission line bundles and range from 7 to 25 diameters and Reynolds number ranges from 5 000 to 27 000. For the tandem cylinders, the results indicate that (a) the dynamic response of the downstream cylinder is no longer hysteretic, (b) the synchronization onset is at higher reduced velocities and (c) the synchronization region is wider than that of an isolated cylinder; this region is twice as large for a separation of 10 diameters and decreases with increasing spacing. Despite the absence of the hysteresis loop and its associated discontinuities, two vortex patterns are present around the resonance velocity. A small decrease of the maximum amplitude of oscillation with increasing spacing has also been observed, as well as a reduction of the velocity for which this maximum occurs. For separations of 7 and 8·5 diameters, the leeward cylinder exhibits a combination of vortex-induced and wake-galloping oscillations. For the staggered arrangement, the increase of the stagger angle β reduces the synchronization onset velocity and the synchronization region and beyond β=20°, the downstream cylinder behaves as an isolated one.     

Wake-induced vibration of tandem and staggered cylinders with two degrees of freedom

The wake-induced vibration (WIV) of two staggered cylinder with two degrees of freedom (2-dof) has been investigated by experiments in a water channel for Reynolds number between 2000 and 25 000. The streamwise separation was fixed to 4 diameters and the lateral separation varied between 0 and 3 diameters for tandem and staggered configurations. Results are presented in the form of trajectories of motion and dynamic response curves of displacements, frequencies and force coefficients. Excitation caused by the WIV mechanism is found to get weaker as the initial position of the downstream cylinder is increased from the centreline of the wake (tandem arrangement) towards the sides. For a lateral separation of 3 diameters wake interference was already found to be negligible. Evidence of a type of wake-stiffness concept is also observed to occur for 2-dof WIV in tandem arrangement, especially for higher reduced velocities. A similar mechanism may also be occurring for staggered arrangements around the centreline.     

Finite-element analysis of fluid flow and heat transfer for staggered bundles of cylinders in cross flow

The two-dimensional Navier-Stokes equations and the energy equation governing steady laminar incompressible flow are solved by a penalty finite-element model for flow across finite depth, five-row deep, staggered bundles of cylinders. Pitch to diameter ratios of 1·5 and 2·0 are considered for cylinders in equilateral triangular and square arrangements. Reynolds numbers studied range from 100 to 400, and a Prandtl number of 0·7 is used. Velocity vector fields, streamline patterns, vorticity, pressure and temperature contours, local and average Nusselt numbers, pressure and shear stress distributions around the cylinder walls and drag coefficients are presented. The results obtained agree well with available experimental and numerical data.     

Forced‐convection heat transfer from tandem square cylinders in cross flow at low Reynolds numbers

This paper presents the results of a numerical study on the flow characteristics and heat transfer over two equal square cylinders in a tandem arrangement. Spacing between the cylinders is five widths of the cylinder and the Reynolds number ranges from 1 to 200, Pr=0.71. Both steady and unsteady incompressible laminar flow in the 2D regime are performed with a finite volume code based on the SIMPLEC algorithm and non‐staggered grid. A study of the effects of spatial resolution and blockage on the results is provided. In this study, the instantaneous and mean streamlines, vorticity and isotherm patterns for different Reynolds numbers are presented and discussed. In addition, the global quantities such as pressure and viscous drag coefficients, RMS lift and drag coefficients, recirculation length, Strouhal number and Nusselt number are determined and discussed for various Reynolds numbers. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical predictions for unsteady viscous flow past an array of cylinders

The unsteady two-dimensional flow around an array of circular cylinders submerged in a uniform onset flow is analysed. The fluid is taken to be viscous and incompressible. The array of cylinders consists of two horizontal rows extending to infinity in the upstream and downstream directions. The centre-to-centre distance between adjacent cylinders is fixed at three diameters, and the rows are staggered. Advantage is taken of spatially periodic boundary conditions in the flow direction. This reduces the computational domain to a rectangular region surrounding a single circular cylinder. Two cases, for Reynolds numbers of 1000 and 10,000, are presented.     

The effect of arrangement of two circular cylinders on the maximum efficiency of Vortex-Induced Vibration power using a Scale-Adaptive Simulation model

The complex behaviour of an unsteady flow around two circular cylinders in tandem is of interest for many civil engineering applications across a wide range of aerospace, mechanical and marine applications. The present paper analyses Vortex-Induced Vibration (VIV) for the flow around two circular cylinders. It has been shown that the amount of kinetic energy which can be captured by VIV is a function of the arrangement of the two cylinders. The upstream cylinder is fixed while the downstream is mounted elastically with one degree of freedom normal to the mean flow direction. The efficiency of the VIV power obtained from downstream cylinder is compared for different arrangement of the cylinders. For this purpose, the longitudinal and lateral distances between the cylinders were varied and the Reynolds number was kept constant. Scale-Adaptive Simulation (SAS) and Shear Stress Transport (SST) CFD models are utilized to analyse the validity of the SAS turbulence model. The results indicate that both turbulence models predict the flow characteristics around the cylinders with reasonable precision; however, the predictions from SAS were more accurate compared to the SST. Based on this comparison, SAS model was chosen as a tool to analyse the VIV response of the downstream cylinder. The location of the downstream cylinder has been altered in the wake of upstream one in order to obtain the optimum efficiency of the VIV power. The results reveal that the arrangement of the cylinders can significantly change the efficiency. It is also observed that cylinders offset from one another show a higher efficiency compared to cylinders with their centres aligned.     

CONTROL OF FLOW PAST BLUFF BODIES USING ROTATING CONTROL CYLINDERS

Computational results for control of flow past a circular cylinder using small rotating cylinders are presented. A well-proven stabilized finite-element method, that has been applied to various flow problems earlier, is utilized to solve the incompressible Navier–Stokes equations in the primitive variables formulation. The formulation is first applied to study flow past an isolated rotating cylinder. Excellent match with experimental results, reported earlier, is observed. It is found that in purely two-dimensional flows, very high lift coefficients can be realized. However, it is observed, via three-dimensional Navier–Stokes simulations, that the end-effects and centrifugal instabilities along the cylinder span lead to a loss of lift and increase in drag. The aspect ratio of the cylinder plays an important role. The flow past a bluff body with two rotating control cylinders is studied using 2-D numerical simulations. The effect of the Reynolds number is studied by carrying out simulations for Re=10²and 10⁴. Finite element meshes with an adequate number of grid points are employed to resolve the flow in the gap between the main and control cylinders. Two values of the gap are considered: 0·01D and 0·075 D, where D is the diameter of the main cylinder. It is observed that when the control cylinders rotate at high speed, such that the tip speed is 5 times the free-stream speed, the flow at Re=100 achieves a steady state. For Re=10⁴, even though the flow remains unsteady, the wake is highly organized and narrower compared to the one without control. The results are in good agreement with the flow-visualization studies conducted by other researchers for bluff bodies using similar control concepts. In all the cases, a significant reduction in the overall drag coefficient and the unsteady aerodynamic forces acting on the main cylinder is observed. Results are also presented for the power requirements of the system for translation and rotation. It is found that the coefficient of power required for the rotation of control cylinders is significant for Re=100 but negligible for Re=10⁴flow. The size of the gap is found to be more critical for the Re=10⁴flows. This study brings out the relevance of the gap as a design parameter for such flow control devices.     

The evolution of starting vortex and secondary separation on an elliptic cylinder

Flow visualization was used to investigate experimentally the evolution process from symmetrical shedding to staggered shedding of the starting vortex and the phenomenon of secondary separation on an elliptic cylinder at moderate Reynolds numbers. The vortex structure of the flow separation was studied. The temporal variation of separation angle and length of wake vortex were given. The photographs and experimental results provided basis for further investigation of the complicated feature of the starting process of unsteady separated flows around an elliptic cylinder. The project supported by the National Natural Science Foundation of China.     

Heat transfer and flow behavior around four staggered elliptic cylinders in cross flow

Experimental and numerical studies were carried out to investigate forced convection heat transfer and flow features around the downstream elliptic cylinder in four staggered cylinders in cross flow. The elliptic cylinders examined had an axis ratio (b/c) of 1:2, and they were arranged with zero angle of attack to the upstream flow. The present heat transfer measurements were obtained by heating only the downstream elliptic cylinder (test cylinder) under the condition of constant heat flux. The testing fluid was air and the Reynolds number based on the major axis length (c) was ranged from 4,000 to 45,570. The tested longitudinal spacing ratio (S_x/c) and the transversal spacing ratio (S_y/b) were in the ranges of 1.5 ≤ S_x/c ≤ 4.0 and 1.5 ≤ S_y/b ≤ 4.0, respectively. The air flow pattern and temperature fields around the four staggered elliptic cylinders were predicted by using CFD software package. Also, a flow visualization study was made to show the flow features around the elliptic cylinders. It was observed that Nu_m of the downstream elliptic cylinder in four staggered cylinders was higher than that of three in-line cylinders for all tested spacing ratios and Reynolds numbers except for Re = 4,000. It was clear that, at lower Reynolds number values (Re < 14,100), the average Nusselt number of the downstream elliptic cylinder in three staggered arrangement was higher than that of the downstream cylinder in four staggered arrangement for all tested spacing ratios. On the other hand, at Re > 14,100, the tested elliptic cylinder in four staggered arrangement had the higher values of the average Nusselt number. Moreover, in four staggered arrangement, the maximum average Nusselt number enhancement ratio (average Nusselt number of the tested downstream cylinder/average Nusselt number of a single elliptic cylinder) was found to be about 2.0, and was obtained for spacing ratios of S_x/c = 2.5, S_y/b = 2.5 and at Re = 32,000. Finally, the average Nusselt number of the tested cylinder in four staggered arrangement was correlated in terms of Reynolds number and cylinder spacing ratios.     

On the origin of wake-induced vibration in two tandem circular cylinders at low Reynolds number

We numerically investigate flow-induced vibrations of circular cylinders arranged in a tandem configuration at low Reynolds number. Results on the coupled force dynamics are presented for an isolated cylinder and a pair of rigid cylinders in a tandem configuration where the downstream cylinder is elastically mounted and free to vibrate transversely. Contrary to turbulent flows at high Reynolds number, low frequency component with respect to shedding frequency is absent in laminar flows. Appearance and disappearance of the vorticity regions due to reverse flow on the aft part of the vibrating cylinder is characterized by a higher harmonic in transverse load, which is nearly three times of the shedding frequency. We next analyze the significance of pressure and viscous forces in the composition of lift and their phase relations with respect to the structural velocity. For both the isolated and tandem vibrating cylinders, the pressure force supplies energy to the moving cylinder, whereas the viscous force dissipates the energy. Close to the excitation frequency ratio of one, the ratio of transverse viscous force to pressure force is found to be maximum. In addition, movement of stagnation point plays a major role on the force dynamics of both configurations. In the case of isolated cylinder, displacement of the stagnation point is nearly in-phase with the velocity. During vortex-body interaction, the phase difference between the transverse pressure force and velocity and the location of stagnation point determines the loads acting on the cylinder. When the transverse pressure force is in-phase with velocity, the stagnation point moves to higher suction region of the cylinder. In the case of the tandem cylinder arrangement, upstream vortex shifts the stagnation point on the downstream cylinder to the low suction region. Thus a larger lift force is observed for the downstream cylinder as compared to the vibrating isolated cylinder. Phase difference between the transverse load and the velocity of the downstream cylinder determines the extent of upstream wake interaction with the downstream cylinder. When the cylinder velocity is in-phase with the transverse pressure load component, interaction of wake vortex with the downstream cylinder is lower compared to other cases considered in this study. We extend our parametric study of tandem cylinders for the longitudinal center-to-center spacing ranging from 4 to 10 diameter.     

Vortex shedding from two finite circular cylinders in a staggered configuration

Wind tunnel experiments were conducted to measure the vortex shedding frequencies for two circular cylinders of finite height arranged in a staggered configuration. The cylinders were mounted normal to a ground plane and were partially immersed in a flat-plate turbulent boundary layer. The Reynolds number based on the cylinder diameter was Re_D=2.4×10⁴, the cylinder aspect ratio was AR=9, the boundary layer thickness relative to the cylinder height was δ/H=0.4, the centre-to-centre pitch ratio was varied from P/D=1.125 to 5, and the incidence angle was incremented in small steps from α=0° to 90°. The Strouhal numbers were obtained behind the upstream and downstream cylinders using hot-wire anemometry. From the behaviour of the Strouhal number data obtained at the mid-height position, the staggered configuration could be broadly classified by the pitch ratio as closely spaced (P/D<1.5), moderately spaced (1.5?P/D?3), or widely spaced (P/D>3). The closely spaced staggered finite cylinders were characterized by the same Strouhal number measured behind both cylinders, an indication of single bluff-body behaviour. Moderately spaced staggered finite cylinders were characterized by two Strouhal numbers at most incidence angles. Widely spaced staggered cylinders were characterized by a single Strouhal number for both cylinders, indicative of synchronized vortex shedding from both cylinders at all incidence angles. For selected staggered configurations representative of closely spaced, moderately spaced, or widely spaced behaviour, Strouhal number measurements were also made along the vertical lengths of the cylinders, from the ground plane to the free end. The power spectra showed that for certain cylinder arrangements, because of the influences of the cylinder–wall junction and free-end flow fields, the Strouhal numbers and flow patterns change along the cylinder.     

结构振动对湍流近尾迹的影响

研究了圆柱绕流中流体与结构的相互作用,侧重结构振动对湍流尾迹的影响,用激光测振仪测量圆柱在升力方向的位移;用热线和LDA（二维）测量湍流的近尾迹,通过变化自由流的速度和圆柱体直径（特征尺寸）来变化雷诺数,用两个振动特性不同的（一个相对刚性,一个相对弹）圆柱来产生尾迹,研究固体结构振动对湍流近尾迹的平均速度场和湍流场的影响,结果表明,结构自由振动对湍流近尾迹场影响明显,该影响随雷诺数的变化不明显。     

基于浸入边界-格子Boltzmann通量求解法的椭圆柱流动特性分析

基于浸入边界-格子Boltzmann通量求解法,开展了雷诺数Re=100不同几何参数下单椭圆柱及串列双椭圆柱绕流流场与受力特性对比研究。结果表明,随长短轴比值的增加,单椭圆柱绕流阻力系数先减小后缓慢上升,最大升力系数则随长短轴比值的增大而减小;尾迹流动状态从周期性脱落涡到稳定对称涡。间距是影响串列圆柱及椭圆柱流场流动状态的主要因素,间距较小时,串列圆柱绕流呈周期性脱落涡状态,而椭圆柱则为稳定流动;随着间距增加,上下游圆柱及椭圆柱尾迹均出现卡门涡街现象,且串列椭圆柱临界间距大于串列圆柱。串列椭圆柱阻力的变化规律与圆柱的基本相同,上游平均阻力大于下游阻力;上游椭圆柱阻力随着间距的变大先减小,下游随间距的变大而增加,当间距达到临界间距时上下游阻力跃升,随后出现小幅度波动再逐渐增加,并趋近于相同长短轴比值下单柱体绕流的阻力。     

Cross-flow past a pair of nearly in-line cylinders with the upstream cylinder subjected to a transverse harmonic oscillation

An experimental investigation is presented for the cross-flow past a pair of staggered circular cylinders, with the upstream cylinder subject to forced harmonic oscillation transverse to the flow direction. Experiments were conducted in a water tunnel with Reynolds numbers, based on upstream velocity, U, and cylinder diameter, D, in the range 1440⩽Re⩽1680. The longitudinal separation between cylinder centres is L/D=2.0, with a transverse separation (for the mean position of the upstream cylinder) of T/D=0.17; the magnitude of the harmonic oscillation is 0.44D peak-to-peak and the nondimensional frequency range of the excitation is 0.05⩽f_eD/U⩽0.44. Flow visualization of the wake-formation region and hot-film measurements of the wake spectra are used to investigate the wake-formation process. An earlier study showed that stationary cylinders in this nearly in-line configuration straddle two very different flow regimes, the so-called shear-layer reattachment (SLR) and induced separation (IS) regimes. The present study, demonstrates that oscillation of the upstream cylinder causes considerable modification of the flow patterns around the cylinders. In particular, the wake experiences strong periodicities at the frequency of the oscillating cylinder; in addition to the usual fundamental lock-in, both sub- and superharmonic resonances are obtained. It is also observed that, although the flow exhibits regions of SLR and IS for excitation frequencies below the fundamental lock-in, for frequencies above the lock-in range the flow no longer resembles either of these flow regimes and vortices are formed in the gap between the cylinders.     

Unsteady mixed convection heat transfer from tandem square cylinders in cross flow at low Reynolds numbers

A two-dimensional numerical study is carried out to understand the influence of cross buoyancy on the vortex shedding processes behind two equal isothermal square cylinders placed in a tandem arrangement at low Reynolds numbers. The spacing between the cylinders is fixed with five widths of the cylinder dimension. The flow is considered in an unbounded medium, however, fictitious confining boundaries are chosen to make the problem computationally feasible. Numerical calculations are performed by using a finite volume method based on the PISO algorithm in a collocated grid system. The range of Reynolds number is chosen to be 50–150. The flow is unsteady laminar and two-dimensional in this Reynolds number range. The mixed convection effect is studied for Richardson number range of 0–2 and the Prandtl number is chosen constant as 0.71. The effect of superimposed thermal buoyancy on flow and isotherm patterns are presented and discussed. The global flow and heat transfer quantities such as overall drag and lift coefficients, local and surface average Nusselt numbers and Strouhal number are calculated and discussed for various Reynolds and Richardson numbers.     

On the flow behaviour of confined finite-length wavy cylinders

Confined aspect-ratio of 6 wavy cylinders with a mean blockage-ratio of 0.5 were studied using time-resolved particle-image velocimetry at a sub-critical Reynolds number of 2700. Wavelengths and wave amplitudes of 2–4 and 0.1–0.3 mean diameters respectively were investigated. Results show that vortices are generally shed from the wavy cylinder and channel walls regularly, reminiscent of the unsteady symmetric flow configuration in confined non-wavy cylinders. Furthermore, vortex formation lengths for confined wavy cylinders are generally shorter than their unconfined counterparts, though their variations with respect to geometrical changes remain consistent with unconfined flow conditions. Gross cross-stream flow behaviour does not differ significantly between confined and unconfined wavy cylinders, indicating that finite-length effects are independent of the present confinement. Confined wavy cylinder wake regions are more sensitive towards geometrical changes and a combination of small wavelength and large wave amplitude leads to significant suppression of coherent cylinder and wall vortex-shedding. This is supported by phase-averaged flow reconstructions derived from Proper Orthogonal Decomposition analysis. Lastly, larger wave amplitudes lead to redistributions of dominant flow energy further downstream and to higher mode numbers, which suggests a causal link to the formation of stronger and more coherent streamwise vortices.