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.     

Experimental study of vortex-induced vibrations of a cylinder near a rigid plane boundary in steady flow

In this study, the vortex-induced vibrations of a cylinder near a rigid plane boundary in a steady flow are studied experimentally. The phenomenon of vortex-induced vibrations of the cylinder near the rigid plane boundary is reproduced in the flume. The vortex shedding frequency and mode are also measured by the methods of hot film velocimeter and hydrogen bubbles. A parametric study is carded out to investigate the influences of reduced velocity, gap-to-diameter ratio, stability parameter and mass ratio on the amplitude and frequency responses of the cylinder. Experimental results indicate： （1） the Strouhal number （St） is around 0.2 for the stationary cylinder near a plane boundary in the sub-criti- cal flow regime; （2） with increasing gap-to-diameter ratio （eo/D）, the amplitude ratio （A/D） gets larger but frequency ratio （f/fn） has a slight variation for the case of larger values of eo/D（eo/D 〉 0.66 in this study）; （3） there is a clear difference of amplitude and frequency responses of the cylin- derbetween the larger gap-to-diameter ratios （e0/D 〉 0.66） and the smaller ones （e0/D 〈 0.3）; （4） the vibration of the cylinder is easier to occur and the range of vibration in terms of Vr number becomes more extensive with decrease of the stability parameter, but the frequency response is affected slightly by the stability parameter; （5） with decreasing mass ratio, the width of the lock-in ranges in terms of Vr and the frequency ratio （f/fn） become larger.     

Interference of vortex-induced vibration and transverse galloping for a rectangular cylinder

The phenomenon of interference between vortex-induced vibration (VIV) and galloping in the transverse degree of freedom was studied in the wind tunnel in the case of a spring-mounted slender rectangular cylinder with a side ratio of 1.5 having the short side perpendicular to the flow. The tests were carried out in a wide Scruton number range, starting from low values and increasing it in small steps by using eddy-current viscous dampers. This study helped understanding the dynamics of the interaction between the two excitation mechanisms and clearly highlighted the transition through four regimes of VIV-galloping interference. It was found that a high value of the mass-damping parameter is required to decouple the ranges of excitation of vortex-induced vibration and galloping completely, and for the quasi-steady theory to predict the galloping critical wind speed correctly. This conclusion is also relevant from the engineering point of view, as it means that structures and structural elements with ordinary mass-damping properties can exhibit sustained vibrations in flow speed ranges where no excitation is predicted by classical theories of vortex-induced vibration and galloping. Although most of the experimental tests were conducted in smooth flow at zero angle of attack, the paper also discusses the sensitivity of the results to a small variation of the mean flow incidence and to the presence of a low-intensity free-stream turbulence.     

Flow characteristics around a circular cylinder subjected to vortex-induced vibration near a plane boundary

Although vortex-induced vibration (VIV) has been extensively studied, much of existing literature deals with uniform flow in the absence of a boundary. The VIV flow field of a structure close to a boundary generally remains unexplored, but it can have important engineering implications, such as pipeline scour if the boundary is an erodible seabed. In this paper, laboratory experiments are performed to investigate the flow characteristics of an elastically mounted circular cylinder undergoing VIV, and a rigid plane boundary is considered to simplify the problem. The initial gap-to-diameter ratio is fixed at 0.8, and six different reduced velocities are considered. The velocity field is measured using a high resolution particle image velocimetry (PIV) system, which has several advantages over traditional PIV systems, including high sampling rate and the ability to mitigate scatter of laser light near the boundary, allowing accurate measurements at the viscous sublayer. This paper presents the vibration amplitude and oscillation frequency for different V_r; in addition, the mean velocity field, turbulence characteristics, vortex behavior, gap flow velocity, and normal/shear stresses on the boundary were measured/calculated, leading to new insights on the flow field behavior.     

Investigation of scale effect for the computation of turbulent flow around a circular cylinder

In order to investigate the scale effect of turbulent flow around a circular cylinder, two similarity numbers （criteria） based on turbulent kinetic and dissipation rates associ- ated with the fluctuation characteristics of turbulence wake are deduced by analyzing the Reynolds averaged NavierStokes equations （RANS）. The RNG k-s models and finite volume method are used to solve the governing equations and the second-order implicit time and upwind space discretization algorithms are used to discrete the governing equations. A numerical computation of flow parameters around a two-dimensional circular cylinder with Reynolds numbers ranging from 102 to l07 is accomplished and the result indicates that the fluctuation of turbulence flow along the center line in the wake of circular cylinder can never be changed with increasing Reynolds numbers when Re ≥ 3 × 10^6. This conclusion is useful for controlling the scale of numerical calculations and for applying model test data to engineering practice.     

不同剪切率来流作用下柔性圆柱涡激振动数值模拟

采用浸入边界法对细长柔性圆柱在线性剪切流条件下的涡激振动进行三维数值模拟。细长柔性圆柱振动采用三维索模型模拟,其两端铰接,质量比为6,长细比为50,无量纲顶张力为496。来流为线性剪切流,剪切率从0到0.024变化,最大雷诺数为250。研究发现:剪切流作用下柔性立管横流向振动表现为驻波模式,而顺流向振动表现为行波-驻波混合模式。随着剪切率增大,振动频谱呈现多频响应,振动能量逐渐向低频转移。阻力系数平均值随着展向变化,脉动阻力系数和升力系数的均方根值均表现为“双峰”模式。流固能量传递系数沿立管轴向的分布表明,振动激励区集中于高流速区,而振动阻尼区多位于低流速区。剪切率较小时,圆柱的泻涡为平行交叉模式;剪切率较大时,圆柱的泻涡为倾斜泻涡模式,且由于泻涡频率沿立管轴向变化,尾流发生涡裂现象,形成泻涡频率不同的胞格结构。      

Numerical study on the suppression of the vortex-induced vibration of an elastically mounted cylinder by a traveling wave wall

In the present paper, the commercial CFD code “Fluent” was employed to perform 2-D simulations of an entire process that included the flow around a fixed circular cylinder, the oscillating cylinder (vortex-induced vibration, VIV) and the oscillating cylinder subjected to shape control by a traveling wave wall (TWW) method. The study mainly focused on using the TWW control method to suppress the VIV of an elastically supported circular cylinder with two degrees of freedom at a low Reynolds number of 200. The cross flow (CF) and the inline flow (IL) displacements, the centroid motion trajectories and the lift and drag forces of the cylinder that changed with the frequency ratios were analyzed in detail. The results indicate that a series of small-scale vortices will be formed in the troughs of the traveling wave located on the rear part of the circular cylinder; these vortices can effectively control the flow separation from the cylinder surface, eliminate the oscillating wake and suppress the VIV of the cylinder. A TWW starting at the initial time or at some time halfway through the time interval can significantly suppress the CF and IL vibrations of the cylinder and can remarkably decrease the fluctuations of the lift coefficients and the average values of the drag coefficients; however, it will simultaneously dramatically increase the fluctuations of the drag coefficients.     

Three-dimensional numerical simulation of vortex-induced vibration of an elastically mounted rigid circular cylinder in steady current

Vortex-induced vibration (VIV) of an elastically mounted rigid circular cylinder in steady current is investigated by solving the three-dimensional Navier–Stokes equations. The cylinder is allowed to vibrate only in the cross-flow direction. The aim of this study is to investigate the variation of the vortex shedding flow in the axial direction of the cylinder and to study the transition of the flow from two-dimensional (2D) to three-dimensional (3D) for VIV of a cylinder. Simulations are carried out for a constant mass ratio of 2, the Reynolds numbers ranging from 150 to 1000 and the reduced velocities ranging from 2 to 12. The three-dimensionality of the flow is found to be the strongest in the upper branch of the VIV response and weakest in the initial branch. The 2S and 2P vortex shedding modes are found to coexist along the cylinder span in the upper branch, leading to strong variations of the lift coefficient in the axial direction of the cylinder. The difference between the flow transition from 2D to 3D in the VIV lock-in regime and that in the wake of a stationary cylinder is identified. The transition mode B found in the wake of a stationary cylinder is also found in the wake of a vibrating cylinder. The critical Reynolds number for flow transition from 2D to 3D of a cylinder undergoing cross-flow VIV at a reduced velocity of 6 is found to be greater than that for a stationary cylinder. For a constant reduced velocity of 6, the wake flow changes from 2D to 3D as the Reynolds number is increased from 250 to 300. Some 2D numerical simulations are performed and it is found that the 2D Navier–Stokes (NS) equations are not able to predict the VIV in the turbulent flow regime, while the 2D Reynolds-averaged Navier–Stokes (RANS) equations improve the results.     

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.     

NUMERICAL STUDY ON THE FLOW AROUND A CIRCULAR CYLINDER WITH SURFACE SUCTION OR BLOWING USING VORTICITY-VELOCITY METHOD

A vorticity-velocity method was used to study the incompressible viscous fluid flow around a circular cylinder with surface suction or blowing. The resulted high order implicit difference equations were effeciently solved by the modified incomplete LU decomposition conjugate gradient scheme (MILU-CG). The effects of surface suction or blowing ' s position and strength on the vortex structures in the cylinder wake, as well as on the drag and lift forces at Reynoldes number Re = 100 were investigated numerically. The results show that the suction on the shoulder of the cylinder or the blowing on the rear of the cylinder can effeciently suppress the asymmetry of the vortex wake in the transverse direction and greatly reduce the lift force; the suction on the shoulder of the cylinder, when its strength is properly chosen , can reduce the drag force significantly , too.     

The effect of end conditions on the vortex-induced vibration of cylinders

In the present investigation we study the effect of end conditions on the vortex-induced vibration of an elastically mounted rigid cylinder. This work was triggered by some initial controlled vibration experiments which showed that spanwise end conditions can have a large effect on measured fluid forces on a cylinder, and this suggested that some of the disparity amongst previous free vibration studies may possibly be attributed to differences in end conditions. In the principal experiments here, we are concerned with a vertical cylinder piercing the clean free surface of a water channel, and attached to a carriage system mounted atop the channel. The upper end of the submerged cylinder is thus the free surface, while the lower end is manipulated to yield three different conditions, namely: an attached endplate; an unattached endplate fixed to the channel floor (with a variable gap between cylinder and plate); and a condition of no endplate at all. Interestingly, we find that the free vibration response for the attached and unattached endplate cases were nearly identical. One expectation was that the case without an endplate would lead to a flow around the end of the body, modifying the vortex dynamics, and thereby reducing the correlation of the induced fluid forces on the body. Surprisingly, over the entire response plot, the vibration amplitude is markedly higher in the absence of an endplate, with the exception of the peak amplitude, which remains nearly unchanged. Unexpectedly, the vibrations become much more steady at flow velocities in the vicinity of the peak response, if the endplate is removed. In a further set of experiments, we undertake controlled vibration, where we vary the gap between cylinder and endplate. We discover a large discontinuous jump in the magnitude of fluid excitation, when the gap exceeds 15% of a diameter. For larger gaps, the fluid excitation becomes independent of the gap size, effectively equivalent to having no plate at all. This study is consistent with some of the disparity between the character of vibration response plots in previous studies, if one takes into account the particular end conditions chosen in those studies.     

On the study of vortex-induced vibration of a cylinder with helical strakes

While the effect of helical strakes on suppression of Vortex-Induced Vibrations (VIV) has been studied extensively, the mechanism of VIV mitigation using helical strakes is much less well documented in the literature. In the present study, a rigid circular cylinder of diameter d=80 mm attached with three-strand helical strakes of dimensions of 10d in pitch and 0.12d in height was tested in a wind tunnel. It was found that the helical strakes can reduce VIV by about 98%. Unlike the bare cylinder, which experiences lock-in over the reduced velocity in the range of 5-8.5, the straked cylinder does not show any lock-in region. In exploring the mechanism of VIV reduction by helical strakes, measurements in stationary bare and straked cylinder wakes using both a single X-probe at four different Reynolds numbers, i.e. Re=10 240, 20 430, 30 610 and 40 800, and two X-probes with variable separations in the spanwise direction at Re=20 430 were conducted. It was found that vortices shed from the straked cylinder are weakened significantly. The dominate frequency varies by about 30% over the range of x/d=10-40 in the streamwise direction while that differs by about 37.2% of the averaged peak frequency over a length of 3.125d in the spanwise direction. The latter is supported by the phase difference between the velocity signals measured at two locations separated in the spanwise direction. The correlation length of the vortex structures in the bare cylinder wake is much larger than that obtained in the straked cylinder wake. As a result, the straked cylinder wake agrees more closely with isotropy than the bare cylinder wake. Flow visualization on the plane perpendicular to the cylinder axis at Reynolds number of about 300 reveals small-scale vortices in the shear layers of the straked cylinder wake. However, these vortices do not roll up and interact with each other to form the well-organized Karman-type vortices. Flow visualization on the plane parallel to the cylinder axis shows vortex dislocation and swirling flow, which should be responsible for the variations of the peak frequency in the streamwise as well as spanwise directions.     

Suppression of vortex shedding for flow around a circular cylinder using optimal control

Adjoint formulation is employed for the optimal control of flow around a rotating cylinder, governed by the unsteady Navier–Stokes equations. The main objective consists of suppressing Karman vortex shedding in the wake of the cylinder by controlling the angular velocity of the rotating body, which can be constant in time or time‐dependent. Since the numerical control problem is ill‐posed, regularization is employed. An empirical logarithmic law relating the regularization coefficient to the Reynolds number was derived for 60?Re?140. Optimal values of the angular velocity of the cylinder are obtained for Reynolds numbers ranging from Re=60 to Re=1000. The results obtained by the computational optimal control method agree with previously obtained experimental and numerical observations. A significant reduction of the amplitude of the variation of the drag coefficient is obtained for the optimized values of the rotation rate. Copyright © 2002 John Wiley & Sons, Ltd.     

Active vortex-induced vibration control of a circular cylinder at low Reynolds numbers using an adaptive fuzzy sliding mode controller

An adaptive fuzzy sliding mode control (AFSMC) scheme is applied to actively suppress the two-dimensional vortex-induced vibrations (VIV) of an elastically mounted circular cylinder, free to move in in-line and cross-flow directions. Laminar flow regime at Re=90, low non-dimensional mass with equal natural frequencies in both directions, and zero structural damping coefficients, are considered. The natural oscillator frequency is matched with the vortex shedding frequency of a stationary cylinder at Re=100. The strongly coupled unsteady fluid/cylinder interactions are captured by implementing the moving mesh technology through integration of an in-house developed User Define Function (UDF) into the main code of the commercial CFD solver Fluent. The AFSMC approach comprises of a fuzzy system designed to mimic an ideal sliding-mode controller, and a robust controller intended to compensate for the difference between the fuzzy controller and the ideal one. The fuzzy system parameters as well as the uncertainty bound of the robust controller are adaptively tuned online. A collaborative simulation scheme is realized by coupling the control model implemented in Matlab/Simulink to the plant model constructed in Fluent, aiming at determination of the transverse control force required for complete suppression of the cylinder streamwise and cross-flow oscillations. The simulation results demonstrate the high performance and effectiveness of the adopted control algorithm in attenuating the 2D-VIV of the elastic cylinder over a certain flow velocity range. Also, the enhanced transient performance of the AFSM control strategy in comparison with a conventional PID control law is demonstrated. Furthermore, the effect of control action on the time evolution of vortex shedding from the cylinder is discussed. In particular, it is observed that the coalesced vortices in the far wake region of the uncontrolled cylinder, featuring the C(2S)-type vortex shedding characteristic mode, are ultimately forced to switch to the classical von Kármán vortex street of 2S-type mode, displaying wake vortices of moderately weaker strengths very similar to those of the stationary cylinder. Lastly, robustness of AFSMC is verified against relatively large structural uncertainties as well as with respect to a moderate deviation in the uniform inlet flow velocity.     

串列双圆柱绕流下游圆柱两自由度涡致振动研究

数值研究了串列双圆柱绕流下游圆柱两自由度涡致振动问题，研究发现：(1) 双自由度的圆柱振幅峰值及出现振峰的频率比都比单自由度的大；(2) 尾流圆柱中的升力远大于均匀来流的，而阻力却相反；(3) 下游圆柱的位移响应对于频率比的变化没有均匀来流中的"敏感"；(4) 尾流中，在频率比1.16和0.87之间，出现了明显的"拍"现象，即圆柱的振幅响应包含不同的频率，而在均匀来流中，并无明显的"拍"现象. 采用ALE方法，计算网格采用H-O非交错网格系统，结合分块耦合方法. N-S方程的对流项和扩散项分别采用三阶迎风紧致格式和四阶中心紧致格式离散. 圆柱振动采用弹簧柱体阻尼器模型，柱体的振动方程采用龙格-库塔法求解. 通过模拟柱体和流体之间的非线性耦合作用，成功地捕捉到了"拍"和"相位开关"等现象.     

Two-degree-of-freedom vortex-induced vibrations of a circular cylinder at Re=3900

The vortex-induced vibrations of an elastically mounted circular cylinder are investigated on the basis of direct numerical simulations. The body is free to move in the in-line and cross-flow directions. The natural frequencies of the oscillator are the same in both directions. The Reynolds number, based on the free stream velocity and cylinder diameter, is set to 3900 and kept constant in all simulations. The behavior of the coupled flow-structure system is analyzed over a wide range of the reduced velocity (inverse of the natural frequency) encompassing the lock-in range, i.e. where body motion and flow unsteadiness are synchronized. The statistics of the structural responses and forces are in agreement with prior experimental results. Large-amplitude vibrations develop in both directions. The in-line and cross-flow oscillations are close to harmonic; they exhibit a frequency ratio of 2 and a variable phase difference across the lock-in range. Distinct trends are noted in the force-displacement phasing mechanisms in the two directions: a phase difference jump associated with a sign change of the effective added mass and a vibration frequency crossing the natural frequency is observed in the cross-flow direction, while no phase difference jump occurs in the in-line direction. Higher harmonic components arise in the force spectra; their contributions become predominant when the cylinder oscillates close to the natural frequency. The force higher harmonics are found to impact the transfer of energy between the flow and the moving body, in particular, by causing the emergence of new harmonics in the energy transfer spectrum.     

基于重叠网格与结构网格的圆柱绕流数值模拟

为研究重叠网格与结构网格在圆柱绕流数值模拟中的区别,以二维圆柱为例,利用有限元分析软件ANSYS 19.2中的DM与Mesh建立模型并划分重叠网格,利用ANSYS 19.2中的ICEM建立模型并划分结构网格。采用FLUENT 19.2中laminar模型模拟分析系统中的平均升力系数、平均阻力系数、斯特劳哈尔数St等流体动力特性。通过改变流体流速得到两种不同网格下各6组雷诺数Re,这6组雷诺数在60~160之间。结果表明:结构网格与重叠网格的St都随着Re的增加而增加,但相同雷诺数下重叠网格对应的St数值更大,St的增长速度更快;重叠网格与结构网格的平均升力系数与平均阻力系数随着Re的增加趋于稳定的速度都加快,但结构网格的平均升力系数与平均阻力系数趋于稳定的速度更快,且两种网格的平均升力系数与平均阻力系数趋于稳定速度的差距逐渐缩小,当Re=160时,两种网格的平均升力系数与平均阻力系数趋于稳定的速度几乎相同;当雷诺数在60~160之间时,采用重叠网格计算出来的斯特劳哈尔数比结构网格更加接近理论值;从升力功率谱密度分布曲线中可以看出,随着雷诺数的增加,两种网格下的频率逐渐变大,并且相同雷诺数下重叠网格的频率比结构网格大。     

Numerical simulation of two-degree-of-freedom vortex-induced vibration of a circular cylinder close to a plane boundary

Two-degree-of-freedom vortex-induced vibrations (VIV) of a circular cylinder close to a plane boundary are investigated numerically. The Reynolds-Averaged Navier-Stokes (RANS) equations are solved using the Arbitrary Lagrangian Eulerian (ALE) scheme with a k-ω turbulence model closure. The numerical model is validated against experimental data of VIV of a cylinder in uniform flow and VIV of a cylinder close to a plane boundary at low mass ratios. The numerical results of the vibration mode, vibration amplitude and frequency agree well with the experimental data. VIV of a circular cylinder close to a plane boundary is simulated with a mass ratio of 2.6 and gap ratios of e/D=0.002 and 0.3 (gap ratio is defined as the ratio of gap between the cylinder and the bed (e) to cylinder diameter (D)). Simulations are carried out for reduced velocities ranging from 1 to 15 and Reynolds numbers ranging from 1000 to 15 000. It is found that vortex-induced vibrations occur even if the initial gap ratio is as small as e/D=0.002, although reported research indicated that vortex shedding behind a fixed circular cylinder is suppressed at small gap ratios (e/D<0.3 or 0.2). It was also found that vibration amplitudes are dependant on the bouncing back coefficient when the cylinder hits the plane boundary. Three vortex shedding modes are identified according to the numerical results: (i) single-vortex mode where the vortices are only shed from the top of the cylinder; (ii) vortex-shedding-after-bounce-back mode; (iii) vortex-shedding-before-bounce-back mode. It was found that the vortex shedding mode depends on the reduced velocity.