A numerical investigation of the near-wake structure in the variable frequency forced oscillation of a circular cylinder

The near-wake structure of a uniform flow past a circular cylinder undergoing a constant-amplitude transverse forced oscillation is studied numerically using a 2-D large eddy simulation (LES) calculation with a Reynolds number range from 500 to 8000. Two effects are considered: First, a comparison is made between the wake structures of periodic and nonperiodic forced oscillations of the cylinder. This was done to emphasize the importance of wake-structure differences of a periodic forced oscillation and a self-excited oscillation of a circular cylinder with the latter being characterized as a nonperiodic forced oscillation. The nonperiodic constant-amplitude forced oscillations were obtained by modulating the frequency of the periodically oscillating cylinder. The differences in the vortex-shedding behavior were made evident by analyzing the vorticity field in the entire wake domain. Second, the effect of changes in the moderate values of the Reynolds number for constant and variable frequency oscillation was investigated. Significant effects on the vortex-shedding patterns in the near wake were observed for both aspects of this study.     

Resonance in flow past oscillating cylinder under subcritical conditions

Flow around an oscillating cylinder in a subcritical region are numerically studied with a lattice Boltzmann method(LBM). The effects of the Reynolds number,oscillation amplitude and frequency on the vortex wake modes and hydrodynamics forces on the cylinder surface are systematically investigated. Special attention is paid to the phenomenon of resonance induced by the cylinder oscillation. The results demonstrate that vortex shedding can be excited extensively under subcritical conditions, and the response region of vibration frequency broadens with increasing Reynolds number and oscillation amplitude. Two distinct types of vortex shedding regimes are observed. The first type of vortex shedding regime(VSR I) is excited at low frequencies close to the intrinsic frequency of flow, and the second type of vortex shedding regime(VSR II)occurs at high frequencies with the Reynolds number close to the critical value. In the VSR I, a pair of alternately rotating vortices are shed in the wake per oscillation cycle,and lock-in/synchronization occurs, while in the VSR II, two alternately rotating vortices are shed for several oscillation cycles, and the vortex shedding frequency is close to that of a stationary cylinder under the critical condition. The excitation mechanisms of the two types of vortex shedding modes are analyzed separately.     

Hydrodynamics of flow over a transversely oscillating circular cylinder beneath a free surface

In this paper, hydrodynamic force coefficients and wake vortex structures of uniform flow over a transversely oscillating circular cylinder beneath a free surface were numerically investigated by an adaptive Cartesian cut-cell/level-set method. At a fixed Reynolds number, 100, a series of simulations covering three Froude numbers, two submergence depths, and three oscillation amplitudes were performed over a wide range of oscillation frequency. Results show that, for a deeply submerged cylinder with sufficiently large oscillation amplitudes, both the lift amplitude jump and the lift phase sharp drop exist, not accompanied by significant changes of vortex shedding timing. The near-cylinder vortex structure changes when the lift amplitude jump occurs. For a cylinder oscillating beneath a free surface, larger oscillation amplitude or submergence depth causes higher time-averaged drag for frequency ratio (=oscillation frequency/natural vortex shedding frequency) greater than 1.25. All near-free-surface cases exhibit negative time-averaged lift the magnitude of which increases with decreasing submergence depth. In contrast to a deeply submerged cylinder, occurrences of beating in the temporal variation of lift are fewer for a cylinder oscillating beneath a free surface, especially for small submergence depth. For the highest Froude number investigated, the lift frequency is locked to the cylinder oscillation frequency for frequency ratios higher than one. The vortex shedding mode tends to be double-row for deep and single-row for shallow submergence. Proximity to the free surface would change or destroy the near-cylinder vortex structure characteristic of deep-submergence cases. The lift amplitude jump is smoother for smaller submergence depth. Similar to deep-submergence cases, the vortex shedding frequency is not necessarily the same as the primary-mode frequency of the lift coefficient. The frequency of the induced free surface wave is exactly the cylinder oscillation frequency. The trends of wave length variation with the Froude number and frequency ratio agree with those predicted by the linear theory of small-amplitude free surface waves.     

流向振荡圆柱绕流的格子Boltzmann方法模拟

用一种新近发展起来的格子Boltzmann方法(LBM)在相对较小的雷诺数(Re \le 200)条件下模拟了不可压缩的流向振荡圆柱绕流问题, 考查了涡脱落模态和升阻力特性. 通过模拟, 在近尾流区发现了实验研究中已经发现的对称/反对称的涡脱落模态, 包括有些传统数值方法未发现的模态. 研究了频率锁定区域的范围及其与振幅的关系, 发现振幅越大, 发生锁定的频率区域越宽. 此外还对升阻力进行了定量意义的模拟,研究了振荡频率和振幅与升阻力的关系.      

单窄条控制件对横向振荡柱体尾流 2P模式 旋涡脱落的改变

横向强迫振荡柱体尾流控制是柱体涡激振动控制的基础,在海洋、土木等工程中具有重要意义. 横向强迫振荡柱体尾流中存在一种锁频旋涡脱落模式,即在一个振荡周期内柱体上、下侧各脱落旋转方向相反的一对涡,称为2P模式. 本文将相对宽度b/D=0.32的窄条控制件置于横向强迫振荡柱体下游,对振幅比A/D=1.25, 无量纲振频f_e D/V_∞=0.22,雷诺数Re=1 200的2P模式旋涡脱落进行干扰,并通过改变控制件位置,研究旋涡的变化规律. 采用二维大涡模拟和实验验证方法进行研究,在控制件位置范围0.8≤X/D≤3.2, 0.4≤Y/D≤3.2内,得到了2P, 2S, P+S和另外6种新发现的旋涡脱落模式,并对各模式旋涡的形成过程作了详细描述. 在控制件位置平面上给出了各旋涡模式的存在区域,画出了旋涡脱落强度的等值线图,并发现在一个相当大的区域内,旋涡脱落强 度可减小一半以上,尾流变窄. 发现柱体大幅振荡引起的横向剪切流在旋涡生成中起关键作用. 探讨了控制件对横向剪切流的影响,分析了控制件在每种旋涡模式形成中的作用机制.      

A numerical study on the flow patterns of two oscillating cylinders

Flow around two oscillating cylinders in a side-by-side arrangement at Reynolds number (Re)=185 is simulated using the immersed boundary method. The purpose of this study is to investigate the combined effects of the gap between the two cylinders and their oscillation in the flow. The cylinders oscillate transversely to a uniform cross-flow with a prescribed sinusoidal function in the opposite direction, with the oscillation amplitude equal to 20% of the cylinder diameter. The gap between the two cylinders and the oscillating frequency are chosen as major variables for the parametric study to investigate their influence on the flow pattern. The ratio of mean gap distance between the two oscillating cylinders to the cylinder diameter is chosen to be 0.6, 1.0, 1.4, and 1.8, and the ratio of oscillating frequencies to the natural vortex shedding frequency of a fixed cylinder is 0.8, 1.0, and 1.2. Wake patterns and the drag and lift coefficients are described and compared with those from a single oscillating cylinder and two stationary cylinders. The wake patterns of two oscillating cylinders can be explained by flow mechanisms of two stationary cylinders, a single oscillating cylinder, and their combinations, and are in agreement with classifications of flow over two stationary cylinders presented in previous studies. In the case of two oscillating cylinders, the modulation phenomenon appears from a lower excitation frequency than in a single oscillating cylinder. Generally, oscillating cylinders have higher drag and root-mean-square (r.m.s.) values of drag coefficients than stationary cylinders.     

均匀来流中旋转圆柱黏性绕流的数值研究

从不可压非定常Ｎ－Ｓ方程出发，首次数值求解了均匀来流中圆柱作周向旋转振荡的黏性绕流问题。探讨了旋转角速度振幅、振荡频率及Ｒｅ数等因素对流场结构及其非定常演化过程的影响，并根据计算结果，给出了在旋转振动频率。速度振幅平面内流场涡结构的分区图。     

均匀来流中旋转圆柱黏性绕流的数值研究

从不可压非定常Ｎ－Ｓ方程出发，首次数值求解了均匀来流中圆柱作周向旋转振荡的黏性绕流问题。探讨了旋转角速度振幅、振荡频率及Ｒｅ数等因素对流场结构及其非定常演化过程的影响，并根据计算结果，给出了在旋转振动频率。速度振幅平面内流场涡结构的分区图。     

Numerical investigation of the turbulent energy budget in the wake of freely oscillating elastically mounted cylinder at low reduced velocities

We present a numerical study of the turbulent kinetic energy budget in the wake of cylinders undergoing Vortex-Induced Vibration (VIV). We show three-dimensional Large Eddy Simulations (LES) of an elastically mounted circular cylinder in the synchronization regime at Reynolds number of Re=8000. The Immersed Boundary Method (IBM) is used to account for the presence of the cylinder. The flow field in the wake is decomposed using the triple decomposition splitting the flow variables in mean, coherent and stochastic components. The energy transfer between these scales of motions are then studied and the results of the free oscillation are compared to those of a forced oscillation. The turbulent kinetic energy budget shows that the maximum amplitude of VIV is defined by the ability of the mean flow to feed energy to the coherent structures in the wake. At amplitudes above this maximum amplitude, the energy of the coherent structures needs to be fed additionally by small scale, stochastic energy in form of backscatter to sustain its motion. Furthermore, we demonstrate that the maximum amplitude of the VIV is defined by the integral length scale of the turbulence in the wake.     

Flow interference between a stationary cylinder and an elastically mounted cylinder arranged in proximity

The flow interference between two circular cylinders, one stationary and the other free to oscillate in the transverse direction, are studied numerically at Re=150. The incompressible Navier–Stokes equation in two space dimensions, an assumption that is expected to be valid at the considered Re, is solved by the characteristic-based-split (CBS) finite element method using the T4/C3 MINI triangular element. The center-to-center spacing between the two cylinders is fixed at 4D, where D is the cylinder diameter. The angle of incident flow, α, with respect to the line through the two cylinder centers, varies within the range from α=0° to 90°. For the elastically mounted cylinder, the reduced mass considered is M_r=2.0; the structural damping coefficient is assigned to be zero, which encourages high amplitude oscillations. For each α, the computations are conducted for a wide range of reduced velocities, U_r. The flow interference is examined by scrutinizing (i) the frequency characteristics of the vortex shedding and oscillation; (ii) the dynamic response of the oscillating cylinder, including the amplitude of displacement, the drag and lift force characteristics and the phase relationship between the lift and the displacement series; and (iii) the flow response in terms of the instantaneous vorticity field. It was found that the flow interference type is significantly affected by the angle of the incident flow. As the cylinder is oscillated outside of the region of the wake behind the stationary cylinder (α≥30°), it behaves similarly to its isolated counterpart. In contrast, if the cylinder is partially or entirely submerged within the upstream wake (α<30°), then both the flow and body responses are substantially modified due to the vigorous interaction between the upstream wake and the oscillating cylinder; the response therefore belongs to the wake-induced regime. The U_r range associated with the higher amplitude response is significantly shifted toward a higher U_r. The maximum vibration amplitude builds up to a significantly higher level, even increasing the U_r far beyond the resonance regime. In general, the wake flow associated with the wake-induced vibration (WIV) regime appears to be more unperiodic than does that corresponding to the vortex-induced oscillation regime. It is also revealed that both the vortex-cylinder and the shear layer-cylinder interaction mechanisms are responsible for the characteristics of the responses of oscillating cylinder. The larger momentum required for the higher oscillation amplitude is obtained from the duration of the energy transfer from the fluid to the cylinder, which is ascribed to the phase lag between the lift force and the cylinder displacement.     

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.     

System identification of the kinematics of an oscillating cylinder using wake velocities

A classical problem in vortex-induced vibration is to know the flow field past an oscillating cylinder. In this paper we use system theory to identify the oscillatory behaviour of a circular cylinder from flow variables in the wake. We use numerical simulations (CFD) of the flow past a cylinder oscillating in the cross-flow direction at different oscillation frequencies and amplitudes to construct a transfer function that relates the displacement of the cylinder and the resulting flow field. This transfer function can then be inverted to ‘predict’ the displacement of the cylinder given the flow field (as determined by simulations or measurements). We investigate this technique in the so-called lock-in region, where the vortex shedding frequency is synchronised with the oscillation frequency of the cylinder.     

Numerical Study of the Flow Behind a Rotary Oscillating Circular Cylinder

A numerical study is performed of flow behind a rotationally oscillating circular cylinder in a uniform flow by solving the two-dimensional incompressible Navier-Stokes equations. The flow behavior in lock-on regime and the timing of vortex formation from the oscillating cylinder are studied. When the frequency of excitation of the cylinder is in the vicinity of the natural vortex formation frequency, a lock-on vortex formation regime appears. As the excitation frequency being increased relative to the natural frequency the initially formed vorticity concentration switches to the opposite side of the cylinder. The effects of oscillating frequency and amplitude on the vortex structures formed in the near wake of the cylinder are also investigated. Based on the present calculated results, some complicated vortex patterns are identified and are consistent with the previous experimental visualizations.     

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.     

Flow past a transversely oscillating square cylinder in free stream at low Reynolds numbers

This paper reports simulation results for free‐stream flow past an oscillating square cylinder at Re=100 and 150, for oscillating‐to‐natural‐shedding frequency ratios of 0.5?f_r?3.0 at a fixed oscillation amplitude of 0.2 of the cylinder width. The transformed governing equations are solved in a non‐inertial frame of reference using the finite volume technique. The ‘lock‐in’ phenomena, where the vortex shedding becomes one with the oscillation frequency, is observed near the natural shedding frequency (f_r≈1). Beyond the synchronization band, downstream recovery of the wake to its stationary (natural) state (frequency) is observed in cross‐stream velocity spectra. At higher forcing frequencies, a phase lag between the immediate and the far wake results in a shear layer having multi‐polar vortices. A ‘Vortex‐switch’ accompanied by a change in the direction of energy transfer is identified at the ‘lock‐in’ boundaries. The variation of aerodynamic forces is noticed to be different in the lock‐in regime. The velocity phase portrait in the far wake revealed a chaotic state of flow at higher excitation though a single (natural) frequency appears in the spectra. Copyright © 2008 John Wiley & Sons, Ltd.     

The study of self-sustained oscillating plane jet flow impinging upon a small cylinder

 Experimental studies of a plane jet impinging upon a small circular cylinder are conducted by hot-wire measurements. The cylinder is located on the jet centerline within the potential-core region. The jet–cylinder interactions on the instability shear layer frequency, the cylinder wake shedding frequency, and the induced self-sustained oscillation phenomenon are carefully investigated. Test data indicate that the self-sustained flow oscillation is mainly generated by the resonant effect of the flow between the jet exit and the cylinder. Its resonant frequency is found to vary linearly and exhibits jump-stage pattern as a function of the distance between the jet exit and the cylinder. The feedback mechanism and the hydrodynamic instability theorem are proposed to predict correctly the frequency jump position, wave number and the convection speed of the self-sustained oscillating flow for different jet exit velocities. Received: 15 July 1998/Accepted: 9 December 1998     

A NUMERICAL SIMULATION OF VORTEX SHEDDING FROM AN OSCILLATING CIRCULAR CYLINDER

Vortex shedding from an oscillating circular cylinder is studied by numerical solutions of the two-dimensional unsteady Navier–Stokes equations. A physically consistent method is used for the reconstruction of velocity fluxes which arise from discrete equations for the mass and momentum balances. This method ensures a second-order accuracy. Two phenomena are investigated and, in both cases, the cylinder oscillation is forced. The first is the flow induced by the harmonic in-line oscillation of cylinder in water at rest. The Reynolds number is equal to 100 and the Keulegan–Carpenter number is equal to 5. A comparison of phase-averaged velocity vectors between measurements and predictions is presented. Applying the widely used model of Morison to the computed in-line force history, the drag and the inertia coefficients are calculated and compared for different grid levels. Using these to reproduce the force functions, deviations from those originally computed are revealed. The second problem is an investigation of a transversely oscillating cylinder in a uniform flow at fixed Reynolds number equal to 185. The cylinder oscillation frequency ranges between 0·80 and 1·20 of the natural vortex-shedding frequency, and the oscillation amplitude is 20% of the cylinder diameter. As the frequency of excitation of the cylinder increases relative to the inherent vortex formation frequency, the initially formed concentration of vorticity moves closer to the cylinder until a limiting position is reached. At this point, the vorticity concentration abruptly switches to the opposite side of the cylinder. This process induces distinct changes of the topology of the corresponding streamline patterns.     

Characterization of long time fluctuations of forces exerted on an oscillating circular cylinder at KC=10

Flow dynamics, in-line and transverse forces exerted on an oscillating circular cylinder in a fluid initially at rest are studied by numerical resolution of the two-dimensional Navier-Stokes equations. The Keulegan-Carpenter number is held constant at KC=10 and Re is increased from 40 to 500. For the different flow regimes, links between flow spatio-temporal symmetries and force histories are established. Besides simulations of long duration show that in two ranges of Re, forces exhibit low frequency fluctuations compared to the cylinder oscillation frequency. Such observations have been only mentioned in the literature and are more deeply examined here. In both ranges, force fluctuations correspond to oscillations of the front and rear stagnation points on the cylinder surface. However, they occur in flow regimes whose basic patterns (V-shaped mode or diagonal mode) have different symmetry features, inducing two distinct behaviors. For 80≤Re≤100, fluctuations are related to a spectral broadening of the harmonics and to a permutation between three vortex patterns (V-shaped, transverse and oblique modes). In the second range 150≤Re≤280, amplitude fluctuations are correlated to the appearance of low frequency peaks interacting with harmonics of the cylinder frequency. Fluctuations are then a combination of a wavy fluctuation and an amplitude modulation. The carrier frequency corresponding to the wavy fluctuation depends on Re and is related to a fluid characteristic time; the modulation frequency is independent of Re and equal to 1/4 of the cylinder oscillation frequency.     

Characteristics of Oscillations in Supersonic Open Cavity Flows

Characteristics of oscillations in supersonic open cavity flows are investigated numerically using hybrid RANS/LES (Reynolds-Averaged Navier-Stokes/Large Eddy Simulation) method. The oscillation regimes and feedback mechanisms for the supersonic cavity flows are identified and analyzed. The calculation captures a mixed shear-layer/wake oscillation mode in the flow of Ma = 1.75, where these two modes occur alternately. The shear-layer mode and wake mode are driven by vortex convection-acoustic feedback and absolute instability, respectively. In particular, the results indicate that the feedback-acoustic-wave in the shear-layer mode is probably generated by the reflection of the downstream-traveling pressure wave, associated with the shed vortex in the shear layer, on the aft wall. The cavity flow of Ma = 2.52 is then simulated to see the influence of Mach number. It is found that the increase of Mach number may decrease the amplitude of the fluctuations in the shear layer, inhibiting the transition to wake mode. Furthermore, the influence of upstream injection is also studied, where the results show that the injection only weakens the oscillations and faintly shifts the resonant frequencies.     

Flow-induced vibrations of a tethered circular cylinder

One of the most basic examples of fluid-structure interaction is provided by a tethered body in a fluid flow. The tendency of a tethered buoy to oscillate when excited by waves is a well-known phenomenon; however, it has only recently been found that a submerged buoy will act in a similar fashion when exposed to a uniform flow at moderate Reynolds numbers, with a transverse peak-to-peak amplitude of approximately two diameters over a wide range of velocities. This paper presents results for the related problem of two-dimensional simulations of the flow past a tethered cylinder. The coupled Navier–Stokes equations and the equations of motion of the cylinder are solved using a spectral-element method. The response of the tethered cylinder system was found to be strongly influenced by the mean layover angle as this parameter determined if the oscillations would be dominated by in-line oscillations, transverse oscillations or a combination of the two. Three branches of oscillation are noted, an in-line branch, a transition branch and a transverse branch. Within the transition branch, the cylinder oscillates at the shedding frequency and modulates the drag force such that the drag signal is dominated by the lift frequency. It is found that the mean amplitude response is greatest at high reduced velocities, i.e., when the cylinder is oscillating predominantly transverse to the fluid flow. Furthermore, the oscillation frequency is synchronized to the vortex shedding frequency of a stationary cylinder, except at very high reduced velocities. Visualizations of the pressure and vorticity in the wake reveal the mechanisms behind the motion of the cylinder.