低雷诺数下串列三圆柱涡激振动中的弛振现象及其影响因素

对间距比为1.2和雷诺数为100的串列三圆柱涡激振动进行数值模拟, 发现在某个折合流速之后, 三圆柱的响应均呈现为随着折合流速增大而增大的弛振现象, 平衡位置偏移、低频振动以及旋涡脱落与圆柱运动之间的时机三个因素共同决定了弛振现象的出现. 进一步的研究发现, 串列三圆柱的弛振现象仅出现在质量比不大于2.0和雷诺数不大于100的工况下. 当质量比较大时, 串列三圆柱的平衡位置固定不变, 且圆柱的振动不规律, 使得旋涡脱落与圆柱运动的时机处于变化之中. 当雷诺数较高时, 最上游圆柱的平衡位置在折合流速较大时回到初始位置, 不再参与对圆柱振动的调节, 使得圆柱的振动响应不再规律, 旋涡脱落与圆柱运动的时机也一直处于变化之中.      

弹性支撑圆柱绕流稳定性分析

基于CFD 技术,采用系统辨识方法,建立了亚临界雷诺数(Re < 47) 下绕圆柱流动的非定常气动力模型(reduced order model, ROM). 耦合结构运动方程和降阶气动力模型,建立了弹性支撑圆柱绕流的稳定性分析模型. 算例分析了亚临界雷诺数下,结构固有频率、质量比等参数以及支撑方式对弹性系统稳定性的影响. 对于单自由度横向支撑圆柱,当结构固有频率趋近流动最不稳定模态频率时,弹性系统会在一定频率范围内失稳,这种现象最低可在Re~20 时出现. 旋转自由度的释放能够进一步降低系统的稳定性,可将临界雷诺数进一步降低至18 左右. ROM 方法不仅具有很高的效率,而且清晰地指出了弹性系统失稳的根本原因:流动模态和结构模态耦合作用导致结构模态失稳所致. 因此,失稳状态下系统振荡频率锁定于结构固有频率. 基于ROM 技术预测的失稳边界与直接CFD/CSD 仿真结果吻合,证明了该方法的正确性和精度.      

不同控制角下附加圆柱对圆柱涡激振动影响

在弹性支撑的圆柱周围布置直径更小圆柱会影响剪切层发展以及旋涡脱落,进而改变其涡激振动状态.通过不同的布置形式和附加小圆柱个数可以实现对圆柱涡激振动的促进或抑制.激励更大幅值的振动可以更好地将水流动能转化为可利用的机械能或电能,抑制其振动则可以实现对海洋平台等结构物的保护.采用基于迭代的嵌入式浸入边界法对前侧对称布置两个小圆柱的圆柱涡激振动进行数值模拟研究,系统仅做横向振动,其中基于主圆柱直径的雷诺数为100,质量比为2.0,折合流速为3~11.小圆柱与主圆柱的直径比为0.125,间隙比为0.125.结果表明,在研究的控制角范围内(30°~90°),附加小圆柱可以很大程度上改变圆柱涡激振动的状态.当控制角较小(30°)时,附加小圆柱对主圆柱的振动起抑制作用;当控制角为45°~60°时,圆柱的振动分为涡振和弛振两个阶段,在弛振阶段,圆柱振幅随折合流速增加而持续增加;当控制角较大(75°~90°)时,附加小圆柱的促进作用随着控制角增加而减小.进一步地,结合一个周期内不同时刻旋涡脱落以及圆周压强分布,解释了附加小圆柱对主圆柱涡激振动的作用机制.应用能量系数对圆柱系统的进一步分析发现,弛振阶段由流体传递到主圆柱的能量系数随折合流速的增加逐渐下降,旋涡结构的改变是产生这种变化的直接原因.      

带有新型涡激振动抑制罩的圆柱体的水动力特性

通过模型实验和数值模拟计算,研究了带有涡激振动抑制罩的圆截面柱体的水动力特性.模型实验主要测试了柱体上附加谐波型和类圆锥型涡激振动抑制罩的单摆结构在不同流速下发生涡激振动的性质;数值模拟则针对谐波型和圆锥型扰动,在雷诺数Re为102到105范围内,研究其水动力参数,如阻力、升力和涡脱落频率等,随扰动波长和波动强度的变化.模型实验结果表明,在直圆柱开始发生共振的流速下,带抑制罩的柱体的振幅显著降低,而在更高流速下则显著增大.数值模拟结果表明,谐波型和圆锥型扰动具有相似的水动力特性;且在不同Re时,阻力、升力和涡脱落频率具有相似的变化规律;随波动强度的增大,阻力一般逐渐增大,升力则在多数情况下先减小而后增大,而涡脱落频率一般逐渐减小.      

低雷诺数下串列三圆柱涡激振动中的弛振现象及其影响因素

对间距比为1.2和雷诺数为100的串列三圆柱涡激振动进行数值模拟,发现在某个折合流速之后,三圆柱的响应均呈现为随着折合流速增大而增大的弛振现象,平衡位置偏移、低频振动以及旋涡脱落与圆柱运动之间的时机三个因素共同决定了弛振现象的出现.进一步的研究发现,串列三圆柱的弛振现象仅出现在质量比不大于2.0和雷诺数不大于100的工况下.当质量比较大时,串列三圆柱的平衡位置固定不变,且圆柱的振动不规律,使得旋涡脱落与圆柱运动的时机处于变化之中.当雷诺数较高时,最上游圆柱的平衡位置在折合流速较大时回到初始位置,不再参与对圆柱振动的调节,使得圆柱的振动响应不再规律,旋涡脱落与圆柱运动的时机也一直处于变化之中.     

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

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

不同控制角下附加圆柱对圆柱涡激振动影响

在弹性支撑的圆柱周围布置直径更小圆柱会影响剪切层发展以及旋涡脱落,进而改变其涡激振动状态.通过不同的布置形式和附加小圆柱个数可以实现对圆柱涡激振动的促进或抑制.激励更大幅值的振动可以更好地将水流动能转化为可利用的机械能或电能,抑制其振动则可以实现对海洋平台等结构物的保护.采用基于迭代的嵌入式浸入边界法对前侧对称布置两个小圆柱的圆柱涡激振动进行数值模拟研究,系统仅做横向振动,其中基于主圆柱直径的雷诺数为100,质量比为2.0,折合流速为3～11.小圆柱与主圆柱的直径比为0.125,间隙比为0.125.结果表明,在研究的控制角范围内(30°～90°),附加小圆柱可以很大程度上改变圆柱涡激振动的状态.当控制角较小(30°)时,附加小圆柱对主圆柱的振动起抑制作用;当控制角为45°～60°时,圆柱的振动分为涡振和弛振两个阶段,在弛振阶段,圆柱振幅随折合流速增加而持续增加;当控制角较大(75°～90°)时,附加小圆柱的促进作用随着控制角增加而减小.进一步地,结合一个周期内不同时刻旋涡脱落以及圆周压强分布,解释了附加小圆柱对主圆柱涡激振动的作用机制.应用能量系数对圆柱系统的进一步分析发现,弛振阶段由流体传递到主圆柱的能量系数随折合流速的增加逐渐下降,旋涡结构的改变是产生这种变化的直接原因.     

亚临界雷诺数下圆柱受迫振动的数值研究

在Navier-Stokes方程和k-ω湍流模型的基础上,利用流线迎风有限元方法结合ALE动网格技术对亚临界雷诺数下的圆柱受迫振动问题开展了数值模拟研究。本文的数值模型成功模拟了Re=5000条件下,圆柱发生受迫振动时尾迹区内的2S,2P和P+S尾流模式;对Re=10000情况下,无量纲振幅分别为0.3,0.4,0.5的圆柱受迫振动问题开展了数值模拟,分析了给定振幅条件下圆柱受力随振动频率的变化关系以及受迫振动的锁定区间。以上数值计算结果与Gopalkrishnan (1993)的实验结果基本符合。研究结果表明,二维数值模型能够基本正确地反映出圆柱发生受迫振动时的涡激振动特性以及有关的受力变化趋势,为今后进一步开展三维数值分析工作奠定了基础。     

质量比对D形截面柱体流致振动的影响

利用涡激振动进行海流能收集的VIVACE装置是新能源领域的研究热点.应用FLUENT软件采用k-ωSST湍流模型和Newmark-β法,通过数值模拟探究了四个质量比(2, 5, 7和10)在迎流角90°下D形截面柱体的流致振动响应,系统分析了D形截面柱体在横流向上的振动幅值、频率、平衡位置偏移量、尾涡脱落模式以及能量转化效率.所模拟的雷诺数范围为288～2880,对应的约化速度为2～20.结果表明,质量比对D形截面柱体流致振动的影响明显,质量比会改变D形截面柱体流致振动的响应分支.质量比越大, D形截面柱体进入驰振对应的约化速度越低;质量比增大, D形截面柱体平衡位置偏移量相对减小.随着约化速度的增大, D形截面柱体出现了涡激振动、涡激振动-驰振及完全驰振等响应分支.在所模拟的范围内, D形截面柱体高能量转化效率出现在涡激振动分支,而不是在驰振分支;在质量比为10且约化速度为4.5时,一级能量转化效率达到最大值44%.相关研究可为VIVACE装置的振子选型提供参考.     

冲蚀地形上悬跨海底管线涡激振动试验研究

为了给海底管线服役期间的安全评估提供科学依据,基于水弹性相似关系,对海底管线服役期间常见水流环境条件下,冲蚀地形上悬跨弹性海底管线涡激振动进行试验研究.采用粘贴于管线内壁的应变片测量振动过程中的应变.通过与平底地形上相关研究结果进行对比分析发现:冲蚀地形上悬跨海底管线涡激振动强度与振动频率均受到冲蚀地形的影响;振动强度随水流折减速度的变化比平底地形要快,在相对冲坑长(悬跨长)L/D为20与L/D为35时涡激振动强度的差别较平底地形上悬跨管线要小得多;管线垂向振动频率随冲蚀地形冲坑长度的减少而增大;冲蚀地形上悬跨管线周围的斯特劳哈尔数St在0.21左右.比相同雷诺数下平底地形上悬跨管线周围St数略大.     

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.     

The ‘void’ structure in the wake of a self-oscillating flexible circular cylinder

We discuss the experimental vortex wake of a flexible circular cylinder undergoing vortex-induced vibration at low Reynolds number and a large cylinder aspect ratio. Hydrogen bubbles formed on the cylinder track the von Karman vortex cores. They show a characteristic ‘void’ structure. We propose a vortex skeleton model that includes a pinch-off of opposite-signed cores. Voids occurred at a node in streamwise vibration when close to an antinode in transverse cylinder vibration. A vibration model predicts the ratio of shedding frequency to natural cylinder vibration frequency necessary for void formation at specific spanwise locations.     

单窄条控制件对横向振荡柱体尾流 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种新发现的旋涡脱落模式,并对各模式旋涡的形成过程作了详细描述. 在控制件位置平面上给出了各旋涡模式的存在区域,画出了旋涡脱落强度的等值线图,并发现在一个相当大的区域内,旋涡脱落强 度可减小一半以上,尾流变窄. 发现柱体大幅振荡引起的横向剪切流在旋涡生成中起关键作用. 探讨了控制件对横向剪切流的影响,分析了控制件在每种旋涡模式形成中的作用机制.      

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.     

Vortex-induced vibrations of a diamond cross-section: Sensitivity to corner sharpness

This paper studies the fluid–structure interaction of an elastically mounted square cross-section cylinder immersed in a free stream. The cross-section is mounted such that its sides are at 45° to the free stream direction, in a “diamond” configuration, and its motion is constrained to the transverse direction relative to the flow direction. Apart from the cross-section, this setup is the same as the majority of single-degree-of-freedom vortex-induced vibration studies of cylinders. Two-dimensional direct numerical simulations of this system have been performed. The Reynolds number based on the point-to-point distance of the cross-section has been fixed at Re=200). Simulations at this Reynolds number allow a direct comparison with previous results from circular cylinders, and therefore focus directly on the impact of the geometry.The sensitivity of the flow, and therefore the motion of the cylinder, to geometrical effects is considered. This is done by rounding the two side corners (those pointing across the flow) at a given radius. For well-rounded corners, the flow behaviour resembles that of a circular cylinder undergoing vortex-induced vibration. However, below a critical radius, the dynamics are considerably altered. Highly disordered and irregular wakes and body motions are observed, as well as a synchronized, periodic P+S wake mode (Williamson and Roshko, 1988), which consists of a pair of vortices on one side, and a single vortex on the other side, shed per oscillation cycle, which results in a non-zero mean lift force. A period-doubled version of this P+S wake is also presented. The spatial structure, and the spatio-temporal symmetries of each of these modes is reported. The results show that even though the spatio-temporal symmetry of the flow is unaffected by the geometry when the body is rigidly mounted (the flow always saturating to a Kármán vortex street) geometric features such as sharp corners can induce a number of spontaneous symmetry breaking bifurcations when the body is elastically mounted. Which of these various modes is observed is shown to be a function of both the corner radius and the spring stiffness, expressed through the reduced velocity.     

Finite element study of vortex‐induced cross‐flow and in‐line oscillations of a circular cylinder at low Reynolds numbers

Vortex‐induced vibrations of a circular cylinder placed in a uniform flow at Reynolds number 325 are investigated using a stabilized space–time finite element formulation. The Navier–Stokes equations for incompressible fluid flow are solved for a two‐dimensional case along with the equations of motion of the cylinder that is mounted on lightly damped spring supports. The cylinder is allowed to vibrate, both in the in‐line and in the cross‐flow directions. Results of the computations are presented for various values of the structural frequency of the oscillator, including those that are sub and superharmonics of the vortex‐shedding frequency for a stationary cylinder. In most of the cases, the trajectory of the cylinder corresponds to a Lissajou figure of 8. Lock‐in is observed for a range of values of the structural frequency. Over a certain range of structural frequency (F_s), the vortex‐shedding frequency of the oscillating cylinder does not match F_s exactly; there is a slight detuning. This phenomenon is referred to as soft‐lock‐in. Computations show that this detuning disappears when the mass of the cylinder is significantly larger than the mass of the surrounding fluid it displaces. A self‐limiting nature of the oscillator with respect to cross‐flow vibration amplitude is observed. It is believed that the detuning of the vortex‐shedding frequency from the structural frequency is a mechanism of the oscillator to self‐limit its vibration amplitude. The dependence of the unsteady solution on the spatial resolution of the finite element mesh is also investigated. Copyright © 1999 John Wiley & Sons, Ltd.     

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.     

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.     

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.