圆柱涡激振动研究进展与展望

圆柱涡激振动广泛存在于机械工程、海洋工程等诸多领域,研究者们取得了许多研究成果,而系统综述圆柱涡激振动的论文距今已近10 a,因此,有必要对近10 a的研究进展进行系统分析。文中系统地总结了近10 a圆柱涡激振动研究成果,阐述了圆柱涡激振动的尾流模态和其对应的响应分支之间的因果关系,分析了影响圆柱涡激振动的关键因素(如质量比、阻尼比和雷诺数)对涡激振动响应的影响,介绍了圆柱涡激振动最大响应振幅的曲线拟合公式及其局限性,最后对圆柱涡激振动的研究方向提出建议。     

分离盘控制圆柱涡激振动的数值模拟研究

圆柱涡激振动问题一直以来备受关注,分离盘作为涡激振动抑制装置得到广泛研究。分离盘长度L与圆柱直径D之比L/D是影响抑制效果的主要因素。运用有限体积法结合RANS方程与一定的湍流模式离散和求解流场,通过编写自定义程序,使用动网格模拟结构物的运动带来的流域边界的变化,针对弹性支撑的圆柱及附加长度为0.5 D的分离盘模型,在约化速度Ur为2.5~13的情况下,对涡激振动及其抑制进行研究。结果表明:分离盘可以抑制甚至消除圆柱涡激振动,99%以上的振幅被抑制;锁定区始点被推后,锁定区变窄;附加分离盘的圆柱阻力和升力被抑制;其斯特鲁哈数(St)稍高于单圆柱St但差别不大。     

单自由度圆柱涡激振动特性及其影响因素分析

对二维刚性圆柱涡激振动进行了系统理论和数值模拟研究,对比了不同因素下圆柱振幅、锁定区间、升阻力系数等关键参数和特征,探究了影响圆柱涡激振动特性的因素。通过线性化分析推导出了锁定区间柱体响应的频率比、振幅与流体荷载和系统物性参数之间的定性关系。为检验理论分析结果,对相关经典实验的结果进行了整理,并用RANS模型开展了相应的数值模拟。理论分析、实验和数值模拟三者的结论相符,可以确认以下结论:当圆柱处在锁定状态时,其振幅大小由组合因子决定。组合因子越小,圆柱的振幅越大,反之则越小,但是组合因子对振幅的影响幅度有限;柱体锁定区间范围的大小由质量比决定,质量比越小,锁定范围越大,反之则越小,质量比对锁定区间范围的大小有显著影响,影响的主要范围是在约化速度较大一端。     

并列双圆柱涡激振动的经验性模型研究

建立一种新的预报并列双圆柱涡激振动响应的经验性模型,根据特定间距比条件下旋涡脱落频率出现分支的现象,提出以两个具有不同固有频率的尾流振子来共同描述结构的近壁尾涡动力特性,同时两个振子均满足van der Pol方程,进而得到结构振子和流体振子的耦合方程组。使用该模型分别对中高质量比和低质量比的并列圆柱涡激振动问题进行数值计算,结果表明,结构的位移响应和最大振动幅值等变化规律与实验结果趋势一致,数值基本吻合。     

不同运动自由度组合串列双圆柱涡激振动数值模拟

基于自主研发的紧致插值曲线CIP (Constrained Interpolation Profile)方法数学模型,对均匀来流条件时不同运动自由度组合下的串列双圆柱涡激振动问题开展二维数值模拟。模型针对雷诺数Re=100,质量比m*=2的串列双圆柱涡激振动问题,选取上、下游圆柱不同运动自由度组合工况进行模拟。重点分析圆柱的升阻力系数、运动位移随折合速度Ur变化的响应。研究表明:当上游圆柱双自由度运动时,随着下游圆柱运动自由度的增加,下游圆柱对上游圆柱涡激振动响应的影响减弱;当下游圆柱双自由度运动时,随着上游圆柱运动自由度的增加,上游圆柱对下游圆柱涡激振动响应的影响变强。研究结果表明圆柱运动自由度组合形式对串列双圆柱涡激振动的影响不可忽略。     

考虑流固耦合时的海底管道悬跨段非线性动力分析

通过对管道的涡激振动试验,提出了考虑流固耦合的非线性涡激升力表达式,并用该式进行了海底管线悬跨段非线性动力响应时程分析。对考虑流固耦合与未考虑流固耦合情况下得到的管道动力响应时程进行对比,算例表明:当管外流场流速与管道顺流向振动速度值较接近时,不考虑流固耦合时的计算结果明显小于考虑流固耦合时的计算结果。分析认为,在管外流场流速与管道顺流向振动速度值较接近的情况下,管道的涡激振动计算宜采用非线性涡激力模型。     

涡激振动潮流能转换装置获能实验研究

涡激振动潮流能转换装置是适用于低流速、浅水地区的海洋潮流能资源开发的一种新型装置。本文对涡激振动装置的工作机理进行研究,利用涡激振动流体力分解模型,建立了涡激振动方程和预测模型并对预测模型进行了验证。建立涡激振动的数学耦合模型,结合理论分析和预测模型设计了模型实验样机,进行水槽实验,测试流速、弹性系数等参数对涡激振动响应的影响,实验数据验证了数学耦合模型的合理性。     

不同排列方式和间隙比下双圆柱流致振动特性研究

采用基于ALE的流固耦合分析方法,将圆柱体振动模型简化为平面双自由度质量-弹簧-阻尼系统,进行了计算参数无关性分析。对不同θ角和间隙比T下的双圆柱流致振动进行数值模拟,并与单圆柱振动计算结果进行比较,得到了不同θ角和间隙比T下中心圆柱的升阻力系数、流场分布和运动轨迹变化情况。通过分析发现小间隙比下由于两圆柱间的相互干扰,改变了中心圆柱的尾涡结构,导致了升力和阻力系数随不同θ角而发生变化;小间隙比下由于两圆柱间的相互干扰影响,均使中心圆柱的轨迹在不同θ角下发生了较大的变化,顺流向的振动频率与单圆柱振动时相比出现了多倍频,导致其运动轨迹发生了变化。     

悬跨海底管道疲劳寿命分析研究

海底管道在服役期间由于各种原因会在某些管段形成悬跨。这些悬跨在海流力作用下,将产生涡激振动。这种涡激振动最终可能导致管道疲劳失效。管道在海流力作用下发生的涡激振动是管道振动和漩涡尾流振动耦合的结果。在建立管道振动模型和Matteoluca尾流振子模型基础上,对管道涡激振动动力响应特性进行分析。依据Miner线性损伤累积理论,采用S—N曲线法分析计算管道疲劳寿命。最后,针对海洋油气开发与生产,提出延长海底管道疲劳寿命的方法和措施。     

二维圆柱涡激振动数值模拟中湍流模式适用性的探讨

采用RANS结合4种湍流模型对低质量比单自由度涡激振动进行数值模拟,对比分析其对该问题的适用性。用四阶Runge-Kutta法离散运动方程,基于动网格技术处理圆柱振荡引起的网格运动,并对壁面条件的处理进行了细致分析。根据Williamson水槽实验,从振幅比、频率比响应、水动力系数及相位突变、尾涡模式等方面对比分析了4种模型的性能和表现,结果表明Standardκ-ω模型与实验差异较大;目前本类研究中运用较少的Realizableκ-ε模型也是可以适用于涡激振动计算的;κ-ω系列模型得到的最大振幅(0.55D)小于κ-ε系列模型的结果(0.87D);SSTκ-ω模型以及2种κ-ε模型都反映出锁定区振动频率fex与泄涡频率fst分离,其中SSTκ-ω模型较为接近实验结果;尾涡模式上,SSTκ-ω模型在各个分支与实验结果一致;总之,各种模型针对不同物理现象各有优势和缺陷。     

Numerical Investigation on Vortex-Induced Rotations of A Triangular Cylinder Using An Immersed Boundary Method

A numerical study of vortex-induced rotations(VIRs) of an equivalent triangular cylinder, which is free to rotate in the azimuthal direction in a uniform flow, is presented. Based on an immersed boundary method, the numerical model is established, and is verified through the benchmark problem of flow past a freely rotating rectangular body.The computation is performed for a fixed reduced mass of m~*=2.0 and the structural stiffness and damping ratio are set to zero. The effects of Reynolds number(Re=25-180) on the characteristics of VIR are studied. It is found that the dynamic response of the triangular cylinder exhibits four distinct modes with increasing Re: a rest position,periodic rotational oscillation, random rotation and autorotation. For the rotational oscillation mode, the cylinder undergoes a periodic vibration around an equilibrium position with one side facing the incoming flow. Since the rotation effect, the outset of vortex shedding from cylinder shifts to a much lower Reynolds number. Further increase in Re leads to 2 P and P+S vortex shedding modes besides the typical 2 S pattern. Our simulation results also elucidate that the free rotation significantly changes the drag and lift forces. Inspired by these facts, the effect of free rotation on flow-induced vibration of a triangular cylinder in the in-line and transverse directions is investigated. The results show that when the translational vibration is coupled with rotation, the triangular cylinder presents a galloping response instead of vortex-induced vibration(VIV).     

The validity of the independence principle applied to the vortex-induced vibration of an inclined cylinder in steady flow

The validity of the independence principle applied to the vortex-induced vibration (VIV) of an inclined cylinder in steady flow is investigated by conducting numerical simulations. In order to create a perfect end-effect-free condition, periodic boundary condition is applied on the two end boundaries that are perpendicular to the cylinder. It is found that the response amplitude and frequency for an inclination angle of α = 45° agree well with their counterparts for α = 0°. The numerical results demonstrated the validity of the independence principle in the case of vortex-induced vibration, which has not been demonstrated by laboratory tests due to the difficulty in avoiding the end effects.     

Numerical evaluation of passive control of VIV by small control rods

Flow past a circular cylinder with multiple small control rods is studied by numerical simulation for Re_D ranging from 1161.3 to 6387.1. The Reynolds-Averaged-Navier–Stokes (RANS) equations and shear stress transport (SST) k − ω turbulence model are used to calculate the vortex field, while a fourth-order Runge–Kutta method is employed for evaluating the structure dynamics of the cylinder group. Comparisons with experimental results demonstrate the validation of this method. This study is concerned with the vortex induced vibration (VIV) suppression efficacy of small control rods placed around a main cylinder. The effects of control rod number, diameter ratio, spacing ratio and Reynolds number on the hydrodynamics and vibration responses of the main cylinder are investigated. The reduced percents of in-line and cross-flow amplitudes and the increased percents of the whole cross-sectional area of cylinders and the drag coefficient are used to give a comprehensive evaluation. Results of simulation indicate that placing small rods with appropriate number at appropriate locations can achieve good suppression effectiveness at a wide range of Reynolds number. The numerical result for the case with nine control rods, diameter ratio of 0.15 and spacing ratio of 0.6 shows the best suppression effect among the cases investigated in this study.     

Vibration of two elastically mounted cylinders of different diameters in oscillatory flow

Vibration of two elastically mounted cylinders in an oscillatory flow at a Keulegan-Carpenter number of 10 is simulated numerically. The two cylinders are rigidly connected with each other and are allowed to vibrate in the cross-flow direction only. The aim of this paper is to identify the effects of the orientation of the cylinders and the gap between the cylinders on the vibration. The two-dimensional Reynolds-Averaged Navier-Stokes equations are solved to predict the flow and the cylinder vibration is predicted using the equation of motion. When the two cylinders are in a tandem arrangement, a combined single pair flow regime and attached pair flow regime are observed as reduced velocity exceeds 10 and this combined regime and the single pair regime occurs intermittently. Periodic vibration is found when the two cylinders are in a staggered arrangement with a 45° flow attack angle. When the two cylinders are in a side-by-side arrangement, a new single vortex regime is observed. This single vortex remains attached to the cylinder surface and rotates around the cylinder. The intermittent switch between this single vortex regime and the single pair regime are observed.     

Vortex-Induced Vibrations of A Long Flexible Cylinder in Linear and Exponential Shear Flows

A numerical study based on a wake oscillator model was conducted to determine the response performance of vortex-induced vibration(VIV) on a long flexible cylinder with pinned-pinned boundary conditions subjected to linear and exponential shear flows. The coupling equations of a structural vibration model and wake oscillator model were solved using a standard central finite difference method of the second order. The VIV response characteristics including the structural displacement, structural frequency, structural wavenumber, standing wave behavior,travelling wave behavior, structural velocity, lift force coefficient and transferred energy from the fluid to the structure with different flow profiles were compared. The numerical results show that the VIV displacement is a combination of standing waves and travelling waves. For linear shear flow, standing waves and travelling waves dominate the VIV response within the low-velocity and high-velocity zones, respectively. The negative values of the transferred energy only occur within the low-velocity zone. However, for exponential shear flow, travelling waves dominate the VIV response and the negative energy occurs along the entire length of the cylinder.     

Reynolds and Mass-Damping Effect on Prediction of the Peak Amplitude of A Freely Vibrating Cylinder

The Reynolds effect and mass-damping effect on the peak amplitude of a freely vibrating cylinder is studied by using forced oscillating data from Gopalkrishnan＇ s research in 1993, in which all experimental cases were carried out at a fixed Reynolds and the tested cylinder was recognized as a body that had no mass and damping. However, the Reynolds and roass-damping are the very important parameters for the peak amplitude of a freely vibrating cylinder. In the present study, a function F is introduced to connect the forced oscillation and free vibration. Firstly the peak amplitude AG^＊ can be obtained from the function F using forced oscillation data of Gopalkrishnan＇ s experimental at Re = 10^4, and then the Reynolds effect is taken into account in the function f（Re）, while the mass-damping effect is considered in the function K（ α ）, where a is the mass-damping ratio. So the peak amplitude of a freely vibrating cylinder can be predicted by the expression： A ^＊ = K（ α ）f（ Re ）AG^＊ . It is found that the peak transverse amplitudes predicted by the above equation agree very well with many recent experimental data under both high and low Reynolds conditions while roass-damping varies. Furthermore, it is seen that the Reynolds number does have a great effect on the peak amplitude of a freely vibrating cylinder. The present idea in this paper can be applied as an update in the empirical models that also use forced oscillation data to predict the vortex induced vibration （VIV） response of a long riser in the frequency domain.     

Seismic response of submerged floating tunnel tether

A mathematical equation for vibration of submerged floating tunnel tether under the effects of earthquake and parametric excitation is presented.Multi-step Galerkin method is used to simplify this equation and the fourth-order Runge-Kuta integration method is used for numerical analysis.Finally,vibration response of submerged floating tunnel tether subjected to earthquake and parametric excitation is analyzed in a few numerical examples.The results show that the vibration response of tether varies with the seismic wave type;the steady maximum mid-span displacement of tether subjected to seismic wave keeps constant when parametric resonance takes place;the transient maximum mid-span displacement of tether is related to the peak value of input seismic wave acceleration.     

Use of Helical Strakes for FIV Suppression of Two Inclined Flexible Cylinders in A Side-by-Side Arrangement

The experimental studies on flow-induced vibrations(FIV) reduction of two side-by-side flexible cylinders inclined at 45° by using the helical strakes were carried out in a towing tank. The main aim of the experiment is to check whether the helical strakes with a pitch of 17.5 D and a height of 0.25 D, which is considered as the most effective vibration suppression device for the isolated cylinder undergoing vortex-shedding, still perform very well to reduce FIV of two inclined flexible cylinders in a side-by-side arrangement. The vibration of two identical inclined cylinders with a mass ratio of 1.90 and an aspect ratio of 350 was tested in the experiment. The center-to-center distance between the two cylinders was 3.0 D. The uniform flow was simulated by towing the cylinder models along the tank.The towing velocity varied from 0.05 to 1.0 m/s with an interval of 0.05 m/s. The maximum Reynolds number can be up to 1.6×104. Three cases were experimentally studied in this paper, including two side-by-side inclined smooth cylinders, only one smooth cylinder fitted with helical strakes in the two side-by-side inclined cylinders system and both two cylinders attached with helical strakes. The variations of displacement amplitude, dominant frequency, FIV suppression efficiency and dominant mode for the two side-by-side inclined cylinders with reduced velocity were shown and discussed.