MODELING THE UNSTEADY LIFT AND DRAG ON A FINITE-LENGTH CIRCULAR CYLINDER IN CROSS-FLOW

Semi-empirical models for unsteady lift and drag are developed to predict the spectral features of the unsteady forces on a finite-length, right circular cylinder in cross-flow. In general, the models consist of two parts; the spatial variation of r.m.s wall pressure on the cylinder, and the correlation lengths which describe the spatial extent of the correlation of the unsteady wall pressures. Experiments were conducted in a low noise wind tunnel as a function of cylinder diameter Reynolds number (19 200<Re<32 000) and the Strouhal number (0·05< St<3·33), to measure the statistics of the unsteady wall pressures on a model cylinder. These results are incorporated into the theoretical models, and predictions of the spectral characteristics of the lift and drag are made. The r.m.s. wall pressures on the cylindrical surface are found to have the largest amplitude near the cylinder end-cap, and on the rearward portion of the cylinder body. The high levels in these locations are attributed to the separated flow region over the end-cap. The circumferential and axial length-scales decrease exponentially with Strouhal number. Both length-scales exhibit maxima near the Strouhal shedding frequency of St=0·21. The axial length-scales are found to depend on the measurement reference location due to the three-dimensional flow and separated flow region near the end-cap. The unsteady lift and drag predictions using the models developed in this work agree well with previously measured unsteady force data measured on inertial hydrophones exposed to flow. The broadband unsteady lift is found to be greater than the broadband unsteady drag by nominally 3dB.     

Numerical investigation on vortex-induced vibration of an elastically mounted circular cylinder at low Reynolds number using the fictitious domain method

A direct numerical simulation of two-dimensional (2D) flow past an elastically mounted circular cylinder at low Reynolds number using the fictitious domain method had been undertaken. The cylinder motion was modelled by a two degree-of-freedom mass–spring–damper system. The computing code was verified against a benchmark problem in which flow past a stationary circular cylinder is simulated. Then, analyses of vortex-induced vibration (VIV) responses, drag and lift forces and the phase and vortex structures were carried out. Results show that the cylinder's non-dimensional cross-flow response amplitude reaches its summit of 0.572 in the ‘lock-in’ regime. The ‘2S’, instead of the ‘2P’, vortex shedding mode is dominated in the ‘lower’ branch for this 2D low-Re VIV. A secondary oscillation is observed in the lift force when ‘lock-in’ occurs. It is shown that this secondary component changes the phase, offset the energy input by the primary component and thus reduces the cylinder responses. Effects of the Skop–Griffin parameter on cylinder responses were also investigated.     

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.     

VORTEX SHEDDING RESONANCE FROM A ROTATIONALLY OSCILLATING CYLINDER

Vortex shedding resonance of a circular cylinder wake to a forced rotational oscillation has been investigated experimentally by measuring the velocity fluctuations in the wake, pressure distributions over the cylinder surface, and visualizing the flow field with respect to cylinder oscillations. The vortex shedding resonance occurs near the natural shedding frequency at small amplitude of cylinder oscillations, while the peak resonance frequency shifts to a lower value with an increase in oscillation amplitude. The drag and lift forces acting on the cylinder at fixed forcing Strouhal number indicate that the phase lag of fluid forces to the cylinder oscillations increases with an increase in oscillation amplitude, supporting the variation of resonance frequency with oscillation amplitude. The comparative study of the measured pressure distributions and the simultaneous flow visualizations with respect to cylinder rotation shows the mechanisms of phase lag, which is due to the strengthened vortex formation and the modification of the surface pressure distributions.     

Numerical study of the suppression mechanism of vortex-induced vibration by symmetric Lorentz forces

In this paper, the electro-magnetic control of vortex-induced vibration (VIV) of a circular cylinder is investigated numerically based on the stream function–vorticity equations in the exponential–polar coordinates attached on the moving cylinder for Re=150. The effects of the instantaneous wake geometries and the corresponding cylinder motion on the hydrodynamic forces for one entire period of vortex shedding are discussed using a drag–lift phase diagram. The drag–lift diagram is composed of the upper and lower closed curves due to the contributions of the vortex shedding but is magnified, translated and turned under the action of the cylinder motion. The Lorentz force for controlling the vibration cylinder is classified into the field Lorentz force and the wall Lorentz force. The symmetric field Lorentz force will symmetrize the flow passing over the cylinder and decreases the lift oscillation, which, in turn, suppresses the VIV, whereas the wall Lorentz force has no effect on the lift. The cylinder vibration increases as the work performed by the lift dominates the energy transfer. Otherwise, the cylinder vibration decreases. If the net transferred energy per motion is equal to zero, the cylinder will vibrate steadily or be fixed.     

带有柔性薄板三维方柱流固耦合的数值模拟研究

采用双向流固耦合方法,对带有柔性薄板三维方柱的流场变化特性进行了研究。通过对比单方柱,分析了带有柔性薄板三维方柱阻力系数、升力系数以及斯特劳哈尔数的变化规律。研究表明,在方柱尾流区域附加一柔性薄板可以使其阻力系数降低34.6%,同时其变化幅值大大减小;其升力系数的均方根减小84.8%,流场脉动大幅度减小;斯特劳哈尔数降低79.5%。研究结果表明,在三维方柱后设置柔性薄板可有效抑制涡脱落,从而改善三维方柱的尾流特性。     

棱柱绕流流动的数值模拟

利用数值方法对长宽比为1／3, 1和3的棱柱绕流在雷诺数为100的非稳态流动特性进行了分析和研究。采用有限体积法对棱柱绕流的二维流动N-S方程进行离散求解,分析和研究了非稳态的棱柱绕流流场,升力系数,阻力系数和涡动特性,数值模拟的结果与相关文献的数据比较吻合。通过上述研究能够为了解棱柱绕流的非稳态流动特性提供有力的帮助。而对棱柱三维流动的模拟分析和对雷诺数的变化对棱柱流动特性的影响进行研究,将为掌握棱柱绕流的工程特性打下基础。     

An experimental study of fluctuating lift on a square-section cylinder oscillating transversely in a uniform stream

Experiments on a square-section cylinder fixed and forced to oscillate transversely in a uniform stream were conducted in a water tank. The Reynolds number of the experiments is in the range of 3·10³ to 10⁴, the amplitude to side length ratioA/D is up to 0.7 and the range of reduced velocity is 4.5<V _r <12. This study aims at investigating the lock-in phenomenon, the fluctuating lift and the phase shift between fluctuating lift and displacement of the oscillating cylinder. The problems on the aeroelastic instability relating to present experimental results have been discussed. The flow visualization clearly shows that there are drastic changes of vortex-shedding from cylinder at the resonance point and the upper end of the lock-in range. The results of the flow visualization give better understanding of the physical mechanism of the phase shift. Project supported by National Natural Science Foundation of China     

Modal analysis of measurements from a large-scale VIV model test of a riser in linearly sheared flow

Large-scale model testing of a tensioned steel riser in well-defined sheared current was performed at Hanøytangen outside Bergen, Norway in 1997. The length of the model was 90 m and the diameter was 3 cm. The aim of the present work is to look into this information and try to improve the understanding of vortex-induced vibrations (VIV) for cases with very high order of responding modes, and in particular to study if and under which circumstances the riser motions would be single-mode or multi-mode. The measurement system consisted of 29 biaxial gauges for bending moment. The signals are processed to yield curvature and displacement and further to identify modes of vibration. A modal approach is used successfully employing a combination of signal filtering and least-squares fitting of precalculated mode-shapes. As a part of the modal analysis, it is demonstrated that the equally spaced instrumentation limited the maximum mode number to be extracted to be equal to the number of instrumentation locations. This imposed a constraint on the analysis of in-line (IL) vibration, which occurs at higher frequencies and involves higher modes than cross-flow (CF). The analysis has shown that in general the riser response was irregular (i.e. broad-banded) and that the degree of irregularity increases with the flow speed. In some tests distinct spectral peaks could be seen, corresponding to a dominating mode. No occurrences of single-mode (lock-in) were seen. The IL response is more broad-banded than the CF response and contains higher frequencies. The average value of the displacement r.m.s over the length of the riser is computed to indicate the magnitude of VIV motion during one test. In the CF direction the average displacement is typically 1/4 of the diameter, almost independent of the flow speed. For the IL direction the values are in the range 0.05–0.08 of the diameter. The peak frequency taken from the spectra of the CF displacement at riser midpoint show approximately to be equal to the Strouhal frequency. The peak frequency in IL direction was typically twice the Strouhal frequency.     

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.     

Frequency effects on lift and drag for flow past an oscillating cylinder

A transversely oscillating cylinder in a uniform flow is modeled to investigate frequency effects of flow-induced wake on lift and drag of the cylinder. Specifically, verified unsteady fluid dynamic simulations using an immersed-boundary method in a fixed Cartesian grid predict the flow structure around the cylinder and reveal how the integration of surface pressure and shear distributions provides lift and drag on the oscillating cylinder. In this study, frequency ranges to be considered are both near and away from the natural frequency of wake vortex shedding. Subsequently, the effects of frequency lock-in, superposition and demultiplication on lift and drag are discussed based on the spectral analysis of time histories of lift and drag.     

NUMERICAL STUDY OF THE FLOW PAST A CYLINDER EXCITED TRANSVERSELY TO THE INCIDENT STREAM. PART 1: LOCK-IN ZONE,HYDRODYNAMIC FORCES AND WAKE GEOMETRY

The numerical study of the flow past a circular cylinder forced to oscillate transversely to the incident stream is presented herein, at a fixed Reynolds number equal to 106. The finite element technique was favoured for the solution of the Navier–Stokes equations, in the formulation where the stream function and the vorticity are the field variables. The cylinder oscillation frequency ranged between 0·80 and 1·20 of the natural vortex-shedding frequency, and the oscillation amplitude extended up to 50% of the cylinder diameter. Since the resolution of the characteristics of synchronized wakes is the focus of the study, the first task is the determination of the boundary of the lock-in region. The computation revealed that, when the cylinder oscillation frequency exceeds the frequency of the natural shedding of vortices, the flow is not absolutely periodic at subsequent cycles but a quasiperiodic flow pattern occurs, which creates difficulty in the determination of the lock-in boundary. The time histories of the drag and lift forces for various oscillation parameters are presented, while the vorticity contours were favoured for the numerical flow visualization. The hydrodynamic forces, the phase angle between the lift force and the cylinder displacement, and the parameters of the wake geometry when steady state was reached, are presented in cumulative diagrams. These diagrams indicate the effect of the oscillation parameters on the hydrodynamic forces and on the wake geometry.     

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

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

UNSTEADY LIFT FORCE ON A TOWED SPHERE

An experimental effort to characterize the broadband flow-induced lift forces on a spherical body that is towed underwater is described. The body itself is the transducer which is comprised of a small geophone encased in a near-neutrally-buoyant sphere, 7·62 cm in diameter. The research described in this paper quantifies the flow-induced unsteady lift force signal as a function of the sphere diameter Reynolds number (7620<Re<34 290) and the Strouhal number (1·5<St<30). It is found that the broadband flow-induced unsteady lift forces are proportional to the product of an area and the dynamic pressure of the flow, as expected. These data are compared to similar data measured previously on a finite-length, right-circular cylinder in cross flow. This comparison indicates that the cylindrical body creates more unsteady side force than does the spherical one, particularly at the lower end of the Strouhal number range.     

Vortex shedding from a circular cylinder near a moving wall

Flow over a circular cylinder near a moving plane wall is simulated numerically. The influence of the moving wall on the vortex shedding from the cylinder is demonstrated and the corresponding mechanism is illustrated using instability theory. A critical gap ratio between the circular cylinder and the moving wall is defined, and a precise method for determining the critical gap ratio is proposed. The drag and lift forces and the pressure coefficient are presented as a function of the gap ratio. The scaling of the Strouhal number is discussed.     

Frequency lock-in phenomenon for oscillating airfoils in buffeting flows

Navier-Stokes based computer simulations are conducted to determine the aerodynamic flow field response that is observed for a NACA0012 airfoil that undergoes prescribed harmonic oscillation in transonic buffeting flows, and also in pre-buffet flow conditions. Shock buffet is the term for the self-sustained shock oscillations that are observed for certain combinations of Mach number and steady mean flow angle of attack even in the absence of structural motion. The shock buffet frequencies are typically on the order of the elastic structural frequencies, and therefore may be a contributor to transonic aeroelastic response phenomena, including limit-cycle oscillations. Numerical simulations indicate that the pre-shock-buffet flow natural frequency increases with mean angle of attack, while the flow damping decreases and approaches zero at the onset of buffet. Airfoil harmonic heave motions are prescribed to study the interaction between the flow fields induced by the shock buffet and airfoil motion, respectively. At pre-shock-buffet conditions the flow response is predominantly at the airfoil motion frequency, with some smaller response at multiplies of this frequency. At shock buffet conditions, a key effect of prescribed airfoil motions on the buffeting flow is to create the possibility of a lock-in phenomenon, in which the shock buffet frequency is synchronized to the prescribed airfoil motion frequency for certain combinations of airfoil motion frequencies and amplitudes. Aerodynamic gain-phase models for the lock-in region, as well as for the pre-shock-buffet conditions are suggested, and also a possible relationship between the lock-in mechanism and limit-cycle oscillation is discussed.     

ENERGY DISTRIBUTION IN MODES IN THE WAKE OF A FINITE-LENGTH CYLINDER BEFORE AND AFTER TRANSITION

The flow in the wake of a finite-length cylinder has been studied experimentally both before and after the transition to turbulence. This instability occurs at Reynolds numbers around 180–190. One end of the cylinder was fixed to the bottom of the test-section of a wind tunnel, whilst the other terminated in the open flow (free end). For these boundary conditions four main frequency modes within the wake can be identified. These are a centre-cell mode at the Strouhal frequency, end-cell modes with a frequency below the Strouhal frequency, a mode exhibiting the difference frequency between the centre-cell and end-cell modes, and a low-frequency mode (appearing only after transition to turbulence). In this work the energy content of these four modes has been determined throughout the wake, both before and after transition to turbulence. For three of the modes, the energy content is the same before and after transition, whereas the low-frequency mode exhibits energy two to four orders of magnitude greater after transition than before. Hence it is clear that the additional turbulence energy appearing in the wake after transition is located predominantly in this low-frequency mode. The appearance of this low-frequency mode is characterized by the simultaneous appearance of a peak in the power spectra of the velocity fluctuations centred about zero frequency (but with finite width). Consequently, the appearance of this zero frequency peak can be taken as the signature of the onset of turbulence. By considering the downstream growth rates of this low-frequency mode, evidence is presented which suggests that transition to turbulence may occur as a result of wake transition in the downstream central plane of the cylinder.     

Numerical simulation of shear effect on vortex shedding behind a square cylinder

This paper is concerned with the numerical simulation of the flow structure around a square cylinder in a uniform shear flow. The calculations were conducted by solving the unsteady 2D Navier–Stokes equations with a finite difference method. The effect of the shear parameter K of the approaching flow on the vortex-shedding Strouhal number and the force coefficients acting on the square cylinder is investigated in the range K=0·0–0·25 at various Reynolds numbers from 500 to 1500. The computational results are compared with some existing experimental data and previous studies. The effect of shear rate on the Strouhal number and the force acting on the cylinder has a tendency to reduce the oscillation. The Strouhal number, mean drag and amplitude of the fluctuating force tend to decrease as the shear rate increases, but show no significant change at low shear rate. Increasing the Reynolds number decreases the Strouhal number and increases the force acting on the cylinder. At high shear rate the shedding frequencies of the fluctuating drag and lift coefficients are identical. © 1997 John Wiley & Sons, Ltd.     

Mode competition in streamwise-only vortex induced vibrations

Time-resolved Particle-Image Velocimetry (PIV) has been used to study mode competition and transient behaviour in the wake of a cylinder experiencing Vortex-Induced Vibrations (VIV) in the streamwise direction. The cylinder response regime contained two branches, occurring above and below the onset of synchronisation between the wake and the cylinder motion (lock-in). During the first branch, the wake exhibited both the S-I mode (in which two vortices are shed simultaneously per vibration cycle) and the alternate A-II mode (similar to the well known von Kármán vortex street). An extended PIV data set acquired in this region revealed mode switching between the S-I and A-II modes. A criterion based on Proper-Orthogonal Decomposition was developed to identify which mode was dominant as a function of time. The A-II mode was found to be dominant for over 90% of the instantaneous fields examined, while the S-I mode appeared to be more unstable.Symmetrically shed vortices were found to rearrange downstream into an alternate structure in which the wake was no longer synchronised to the cylinder motion. The dominant frequency of transverse velocity fluctuations was measured throughout the wake in order to study the effects of this breakdown in more detail. For the majority of the wake, the fluctuations occurred at the Strouhal frequency, while in a region in the near wake the fluctuations occurred at the frequency of the cylinder motion. It is thought that during the first response branch vortices are formed at the cylinder response frequency, but tend to quickly rearrange downstream into an alternate structure which is no longer synchronised to the cylinder motion. As a result, the fluctuating drag will be synchronised to the structural motion, and is capable of providing positive energy transfer in the apparent absence of lock-in. Finally, the spatial dependence of the frequency of velocity fluctuations throughout the wake is used to explain some of the conflicting results in the literature regarding streamwise VIV, and the implications for the general study of VIV are discussed.     

Application of Differing Forcing Function Models on Simulated Flow Past an Oscillating Cylinder in a Uniform Low Reynolds Number Flow

A Computational Fluid Dynamics (CFD) model is presented for the uniform viscous two dimensional flow past an oscillating cylinder at low Reynolds number. Numerical simulations are made to study the effect of differing forced induced oscillation mechanisms with a large range of cylinder forcing frequencies. In the first case sinusoidal velocity slip boundary conditions are adopted for the cylinder surface to simulate cylinder oscillation. The implication suggests that no modification or additional term need to be added to the Navier-Stokes equations. In the second case this time extra body force terms which are assumed to account for velocity effects due to cylinder movement are included in the Navier-Stokes equations with the imposition of same boundary conditions. Drag and lift coefficients are extracted from present numerical results and other detailed computations of these coefficients are made at a Reynolds number of 80 and an amplitude-to diameter ratio 0.14. The results are found to be in agreement with each other at low force driving frequencies below and near lock-in. However, differences are found at higher frequencies above lock-in. Agreement are also found with experimental results at some frequency ranges.