Application of wake oscillators to two-dimensional vortex-induced vibrations of circular cylinders in oscillatory flows

A nonlinear time-domain simulation model for predicting two-dimensional vortex-induced vibration (VIV) of a flexibly mounted circular cylinder in planar and oscillatory flow is presented. This model is based on the utilization of van der Pol wake oscillators, being unconventional since wake oscillators have typically been applied to steady flow VIV predictions. The time-varying relative flow–cylinder velocities and accelerations are accounted for in deriving the coupled hydrodynamic lift, drag and inertia forces leading to the cylinder cross-flow and in-line oscillations. The system fluid–structure interaction equations explicitly contain the time-dependent and hybrid trigonometric terms. Depending on the Keulegan–Carpenter number (KC) incorporating the flow maximum velocity and excitation frequency, the model calibration is performed, entailing a set of empirical coefficients and expressions as a function of KC and mass ratio. Parametric investigations in cases of varying KC, reduced flow velocity, cylinder-to-flow frequency ratio and mass ratio are carried out, capturing some qualitative features of oscillatory flow VIV and exploring the effects of system parameters on response prediction characteristics. The model dependence of hydrodynamic coefficients on the Reynolds number is studied. Discrepancies and limitations versus advantages of the present model with different feasible solution scenarios are illuminated to inform the implementation of wake oscillators as a computationally efficient prediction model for VIV in oscillatory flows.     

柔性圆柱涡激振动流体力系数识别及其特性

涡激振动是诱发海洋立管、浮式平台系泊缆和海底悬跨管道等柔性圆柱结构疲劳损伤的重要因素.目前,海洋工程中用于柔性圆柱涡激振动预报的流体力系数主要来源刚性圆柱横流向受迫振动的实验数据,存在一定缺陷和误差.本文综合考虑横流向与顺流向振动耦合作用,建立了柔性圆柱涡激振动流体力模型,运用有限元法和最小二乘法确定升力系数、脉动阻力系数和附加质量系数.为了准确识别柔性圆柱涡激振动流体力系数,设计并开展了拖曳水池模型实验,实验用柔性圆柱模型的质量比为1.82,长径比为195.5.通过与刚性圆柱流体力系数对比,深入分析了柔性圆柱流体力系数的特性.结果表明:柔性圆柱在一阶模态控制区,流体力系数随约化速度变化趋势与刚性圆柱大致相似;二阶模态控制区,升力系数和脉动阻力系数显著增大;附加质量系数在响应频率较低时与振动位移的相关性增强;当响应频率较低时,振动位移较大区域为能量耗散区,当响应频率较高时,振动位移较大区域为能量输入区.     

Time domain simulation of vortex-induced vibrations in stationary and oscillating flows

This paper focuses on the further development of a previously published semi-empirical method for time domain simulation of vortex-induced vibrations (VIV). A new hydrodynamic damping formulation is given, and the necessary coefficients are found from experimental data. It is shown that the new model predicts the observed hydrodynamic damping in still water and for cross-flow oscillations in stationary incoming flow with high accuracy. Next, the excitation force model, which is one component of the total hydrodynamic force model, is optimized by simulating the VIV response of an elastic cylinder in a series of experiments with stationary flow. The optimization is performed by repeating the simulations until the best possible agreement with the experiments is found. The optimized model is then applied to simulate the cross-flow VIV of an elastic cylinder in oscillating flow, without introducing any changes to the hydrodynamic force modeling. By comparison with experiment, it is shown that the model predicts the frequency content, mode and amplitude of vibration with a high level of realism, and the amplitude modulations occurring at high Keulegan–Carpenter numbers are well captured. The model is also utilized to investigate the effect of increasing the maximum reduced velocity and the mass ratio of the elastic cylinder in oscillating flow. Simulations show that complex response patterns with multiple modes and frequencies appear when the maximum reduced velocity is increased. If, however, the mass ratio is increased by a factor of 5, a single mode dominates. This illustrates that, in oscillating flows, the mass ratio is important in determining the mode participation at high maximum reduced velocities.     

正弦振荡来流下柔性立管涡激振动发展过程

在风浪流的作用下,海洋浮式结构物将带动悬链线立管在水中作周期性往复运动,从而在立管运动方向上产生相对振荡来流,这种振荡来流将激励立管悬垂段发生“间歇性” 的涡激振动. 在海洋工程水池中对一个4m 长的立管微段进行模型试验研究,以探索相对振荡来流作用下立管涡激振动产生的机理及其发展的物理过程. 试验通过振荡装置带动模型作正弦运动来模拟不同最大约化速度U_Rmax、不同KC（Keulegan-Carpenternumber）的相对振荡来流,利用光纤应变片测量立管涡激振动响应. 结合模态分析方法处理试验数据得到位移响应时历,继而提出相对振荡来流下柔性立管涡激振动发展的3 个阶段:建立阶段、锁定阶段以及衰减阶段. 并进一步总结了最大约化速度U_Rmax,KC 对涡激振动发展过程的影响规律. 最终获得不同最大约化速度U_Rmax下,涡激振动各发展阶段随KC 所占时间分布比例图.      

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.     

Suppression of vortex-induced vibration of a circular cylinder using suction-based flow control

In the present study, a flow control method is employed to mitigate vortex-induced vibration (VIV) of a circular cylinder by using a suction flow method. The VIV of a circular cylinder was first reproduced in a wind tunnel by using a spring–mass system. The time evolution of the cylinder oscillation and the time histograms of the surface pressures of 119 taps in four sections of the circular cylinder model were measured during the wind tunnel experiments. Four steady suction flow rates were used to investigate the effectiveness of the suction control method to suppress VIV of the circular cylinder. The vibration responses, the mean and fluctuating pressure coefficients, and the resultant aerodynamic force coefficients of the circular cylinder under the suction flow control are analyzed. The measurement results indicate clearly that the steady suction flow control method exhibits excellent control effectiveness and can distinctly suppress the VIV by dramatically reducing the amplitudes of cylinder vibrations, fluctuating pressure coefficients and lift coefficients of the circular cylinder model. By comparing the test cases with different suction flow rates, it is found that there exists an optimal suction flow rate for the maximum VIV control. The cases with higher suction flow rates do not necessarily behave better than those with lower suction flow rates. With the experimental setting used in the present study, the suction flow control method is found to behave better for VIV suppression when the ratio of the suction flow velocity to the oncoming flow velocity is less than one.     

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.     

A harmonic balance technique for the reduced-order computation of vortex-induced vibration

We present a harmonic balance (HB) method to model frequency lock-in effect during vortex-induced vibration (VIV) of elastically mounted circular cylinder and a flexible riser section in a freestream uniform flow. The fluid flow and structure are coupled by a fixed-point iteration process through a frequency updating algorithm. By minimizing the structural residual in the standard least-square norm, the convergence of HB-based fixed-point algorithm is achieved for a range of reduced velocity. To begin with, the HB solver is first assessed for a periodic unsteady flow around a stationary circular cylinder. A freely vibrating circular cylinder is then adopted for the reduced-order computation of VIV at low Reynolds numbers of Re=100 and 180 with one- and two-degrees-of-freedom. The coupled VIV dynamics and the frequency lock-in phenomenon are accurately captured. The results show that the HB solver is able to predict the amplitude of vibration, frequency and forces comparable to its time domain counterpart, while providing a significant reduction with regard to overall computational cost. The proposed new scheme is then demonstrated for a fully-coupled three dimensional (3D) analysis of a linear-elastic riser section undergoing vortex-induced vibration in the lock-in range. The results reveal the 3D effects through isosurfaces of streamwise vorticity blobs distributed over the span of flexible riser section. In comparison to time domain results, the 3D flow-structure interactions are accurately predicted while providing a similar speed up rate that of 2D simulations. This further corroborates that the HB solver can be extended to 3D flow-structure dynamics without compromising efficiency and accuracy.     

剪切流下发生涡激振动的柔性立管阻力特性研究

利用缩尺模型试验的方法研究了线性剪切流下涡激振动发生时柔性立管的阻力特性.文中基于光纤光栅应变传感器测得的模型应变信息,采用梁复杂弯曲理论计算了立管的平均阻力,继而分析了阻力系数沿管长方向和雷诺数的分布特性以及涡激振动对阻力系数的放大效应,并提出了用于估算柔性立管发生涡激振动时阻力系数的经验公式.结果表明:涡激振动对阻力系数有放大效应,使得立管局部阻力系数高达3.2;平均阻力系数在1.0×104到1.2×105的雷诺数区间内的值为1.3~2.0,并随雷诺数的增大而减小.本文提出的经验公式可准确估算高雷诺数下涡激振动发生时柔性立管的阻力系数,此经验公式考虑了流速、涡激振动主导模态以及主导频率对阻力系数的影响.      

On shear flow single mode lock-in with both cross-flow and in-line lock-in mechanisms

A long flexible cylinder exposed to ocean currents is known to undergo vortex-induced vibration (VIV). In a spatially sheared flow the response of a riser to VIV can vary from single mode lock-in to multimodal. A new experimental facility was designed and built to investigate the above-mentioned areas. The facility consisted of a long flexible cylinder in either a uniform or a simplified vertically sheared flow. The instrumentation consisted of direct local fluid force measurement at two locations on the cylinder as well as accelerometers spaced along the cylinder axis. The simplified shear flow was a 2-slab flow, with each slab having uniform velocity. Test conditions included forcing the cylinder simultaneously at resonance in both regions to investigate modal competition issues and multimodal response patterns. Resonant VIV excitation of two different modes simultaneously, was conducted which revealed single mode lock-in of the higher frequency through an unexpected mechanism. The higher frequency mode's damping region underwent in-line excitation at four times the predicted shedding frequency that provided a power-in effect to support the dominant mode's cross-flow response.     

Instability analysis of modulated Taylor vortices

This study investigates the instability analysis of modulated Taylor vortices flow by utilising a numerical method. Based on the consideration that the outer cylinder is fixed and the inner cylinder rotates at a non-zero averaged speed under varying modulated amplitudes and frequencies, the flow is converted from one-dimension Couette flow to Taylor vortices. When the modulated amplitude is greater than 1 and the rotation speed of the inner cylinder exceeds the threshold value for one-dimensional flow, the flow will be more stable at intermediate and high frequencies. When the modulated amplitude is sufficiently large and the inner cylinder rotates at medium frequency, subharmonic flow arises.     

近壁面柱体涡激振动的迟滞效应

柱体涡激振动是典型的流固耦合问题,其响应规律标识码在升速流动和远离壁面条件下获得的. 而自然环境流动通常不断经历升速和降速过程,近壁面柱体的涡激振动可呈现与远离标识码体不同的响应特征. 本研究结合大型波流水槽,设计了具有微结构阻尼的柱体涡激振动装置. 基于量纲分析,开展系列水槽标识码验,通过同步测量柱体涡激振动位移时程和绕流流场变化,研究了升降流速作用下柱体涡激振动触发和停振的临界速度(即上临标识码临界速度)变化规律,探究了近壁面柱体涡激振动迟滞效应. 采用自下向上激光扫射的 PIV 流场测量系统,对比分析了固定柱体标识码振动柱体的绕流特征. 实验观测表明,近壁面柱体涡激振动触发的临界速度呈现随壁面间距比减小而逐渐减小的变化趋势;但标识码速条件下的涡激振动停振所对应的下临界速度却明显小于升速时的涡激振动触发所对应的上临界速度. 采用上临界与下临界约标识码差值可定量表征涡激振动迟滞程度,研究发现该值随着柱体间距比减小呈线性增大趋势. 涡激振动迟滞现象通常伴随振幅阶跃标识码阶跃值则随着间距比减小而非线性减小.      

Numerical investigation of two-degree-of-freedom vortex-induced vibration of a circular cylinder in oscillatory flow

Two-degree-of-freedom (2dof) vortex-induced vibration (VIV) of a circular cylinder in oscillatory flow is investigated numerically. The direction of the oscillatory flow is perpendicular to the spanwise direction of the circular cylinder. Simulations are carried out for the Keulegan–Carpenter (KC) numbers of 10, 20 and 40 and the Reynolds numbers ranging from 308 to 9240. The ratio of the Reynolds number to the reduced velocity is 308. At KC=10, the amplitude of the primary frequency component is much larger than those of other frequency components. Most vibrations for KC=20 and 40 have multiple frequencies. The primary frequency of the response in the cross-flow direction decreases with the increasing reduced velocity, except when the reduced velocity is very small. Because the calculated primary frequencies of the response in the cross-flow direction are multiple of the oscillatory flow frequency in most of the calculated cases, the responses are classified into single-frequency mode, double-frequency mode, triple frequency mode, etc. If the reduced velocity is in the range where the VIV is transiting from one mode to another, the vibration is very irregular.For each KC number the range of the reduced velocity can be divided into a cross-flow-in-phase regime (low V_r), where the response and the hydrodynamic force in the cross-flow direction synchronize, and a cross-flow-anti-phase regime (high V_r), where the response and the hydrodynamic force in the cross-flow direction are in anti-phase with each other. The boundary values of V_r between the cross-flow-in-phase and the cross-flow-anti-phase regimes are 7, 9 and 11 for KC=10, 20 and 40, respectively. For KC=20, another cross-flow-anti-phase regime is found between 15≤V_r≤19. Similarly the in-line-in-phase and the in-line-anti-phase regimes are also identified for the response in the in-line direction. It is found that the boundary value of V_r between the in-line-in-phase and the in-line-anti-phase regimes is greater than that in the cross-flow direction. They are 14 and 26 for KC=10 and 20, respectively. Maximum amplitude occurs at the boundary value of the reduced velocity between in-phase regime and anti-phase regime in both the x- and the y-directions.     

两端铰接的细长柔性圆柱体涡激振动响应特性数值研究

目前细长柔性圆柱体涡激振动响应的研究方法主要包括实验方法、计算流体动力学方法以及半经验模型方法. 鉴于实验方法研究成本较高、计算流体动力学方法计算时间较长,本文基于尾流振子模型对线性剪切来流下两端铰接的细长柔性圆柱体涡激振动响应特性进行了半经验模型方法研究. 先建立了柔性圆柱体结构振子以及尾流振子之间的耦合模型,紧接着基于二阶精度中心差分格式对耦合模型先离散后迭代进行求解. 对不同剪切参数下柔性圆柱体涡激振动响应的振动波长、振动频率、振动位移以及响应频率随时间的变化特性等参数进行了分析. 分析结果表明:圆柱体的涡激振动响应由驻波和行波混合组成. 当无量纲弯曲刚度较小时,在圆柱体两端附近,驻波占主导;而在圆柱体中间段附近,行波占主导. 当无量纲弯曲刚度较大时,在圆柱体整个长度区间上均为驻波占主导. 随着剪切参数的增大,振动位移以及振动波长均逐渐减小,而振动频率和频率带宽均逐渐增大.      

Throughflow and g-jitter effects on binary fluid saturated porous medium

A non-autonomous complex Ginzburg-Landau equation (CGLE) for the finite amplitude of convection is derived, and a method is presented here to determine the amplitude of this convection with a weakly nonlinear thermal instability for an oscillatory mode under throughflow and gravity modulation. Only infinitesimal disturbances are considered. The disturbances in velocity, temperature, and solutal fields are treated by a perturbation expansion in powers of the amplitude of the applied gravity field. Throughflow can stabilize or destabilize the system for stress free and isothermal boundary conditions. The Nusselt and Sherwood numbers are obtained numerically to present the results of heat and mass transfer. It is found that throughflow and gravity modulation can be used alternately to heat and mass transfer. Further, oscillatory flow, rather than stationary flow, enhances heat and mass transfer.     

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.     

Parametric excitations and lock-in of flexible hydrofoils in two-phase flows

This work introduces a reduced-order method to study the parametric excitations and lock-in of flexible hydrofoils caused by unsteady two-phase (cavitating) flow. The reduced-order method is based on a 1-DOF structural model coupled with a van der Pol wake oscillator with empirically derived relations for the variation in lift, cavity-length, and cavity-shedding frequency as a function of a non-dimensional cavitation parameter. The results are compared with several available data from both numerical simulations and experimental measurements. The frequency content of both the predicted and measured vibrations suggested that, in addition to the primary cavity-shedding frequency and the hydrofoil natural frequencies, unsteady two-phase flows may excite additional modulated frequencies due to time-varying fluid-added mass effects. The results show that these frequency modulations might cause the flexible hydrofoil to undergo higher-order resonances, as well as parametric resonances. While the maximum deformations for the primary and higher-order resonances were observed to damp out, parametric resonances might persist even with realistic fluid damping coefficients (4–12%). It was observed that with higher effective foil flexibility, the cavity-shedding frequencies may be significantly modified from the rigid foil trends, and may instead lock-in with the system natural frequencies.     

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.     

剪切流场中含内流立管横向涡激振动特性

立管是海洋工程中输送油气或其他矿产资源的必备结构, 外部洋流引起的立管涡激振动影响着立管的疲劳寿命, 危害深海资源开发. 本文基于欧拉?伯努利梁方程, 结合半经验时域水动力模型, 建立剪切流与内流耦合作用下海洋立管涡激振动预报模型, 运用有限元方法和Newmark-β逐步积分法求解方程, 首先将数值模拟结果与实验数据进行对比, 验证模型正确性. 然后, 运用此模型, 对剪切流作用下含内流的顶张立管在不同内流速度和密度下的横向涡激振动响应特性进行研究, 主要分析了立管的横向振动模态、振动频率以及均方根位移等涡激振动参数随内流速度和密度等参数的变化规律. 结果表明, 在剪切流场中, 含内流海洋立管在横向上表现出多模态多频率的涡激振动;立管横向振动的最大均方根位移随内流速度和密度的增大而增大, 特别是当内流速度较大时, 横向最大均方根位移增大明显;立管横向振动的主导频率随内流速度和密度的增大而减小, 并且内流密度的增大同样会引起模态转换和频率转换.      

An experimental investigation of one- and two-degree of freedom VIV of cylinders

In the paper,an experiment investigation was conducted for one-and two-degree of freedom vortex-induced vibration(VIV) of a horizontally-oriented cylinder with diameter of 11 cm and length of 120 cm.In the experiment,the spring constants in the cross-flow and in-line flow directions were regulated to change the natural vibration frequency of the model system.It was found that,in the one-degree of freedom VIV experiment,a "double peak" phenomenon was observed in its amplitude within the range of the reduced velocities tested,moreover,a "2T" wake appeared in the vicinity of the second peak.In the two-degree of freedom VIV experiment,the trajectory of cylinder exhibited a reverse "C" shape,i.e.,a "new moon" shape.Through analysis of these data,it appears that,besides the non-dimensional in-line and cross-flow natural vibration frequency ratios,the absolute value of the natural vibration frequency of cylinder is also one of the important parameters affecting its VIV behavior.