Simulation of Vortex-Induced Vibrations of A Cylinder Using ANSYS CFX Rigid Body Solver

This article simulates the vortex-induced oscillations of a rigid circular cylinder with elastic support using the new ANSYS CFX rigid body solver. This solver requires no solid mesh to setup FSI (Fluid Structure Interaction) simulation. The two-way case was setup in CFX only. Specific mass of the cylinder and flow conditions were similar to previous experimental data with mass damping parameter equal to 0.04, specific mass of 1 and Reynolds number of 3800. Two dimensional simulations were setup. Both one-degree-of-freedom and two-degree-of-freedom cases were run and results were obtained for both cases with reasonable accuracy as compared with experimental results. Eight-figure XY trajectory and lock-in behavior were clearly captured. The obtained results were satisfactory.     

Validation of a finite difference method for the simulation of vortex-induced vibrations on a circular cylinder

The Karman Vortex Street generated by a circular cylinder is investigated by the numerical solution of the compressible Navier–Stokes equations in the incompressible Mach number range (Mach<0.3) using the Beam and Warming implicit scheme. The agreement with the fully incompressible projection method (Chorin, 1968) is fairly good while convergence time is very much better. The investigation suggests that the compressible Navier–Stokes equations may be used as an efficient alternative to study incompressible flows as well. Mach numbers just below 0.3 are enough to simulate incompressible flow behavior and at the same time do not cause numerical ill-conditioning in the solution. In addition, some relevant features of the vortices generated and carried by the wake of the cylinder could be fairly well captured.     

Flow induced vibrations of a sharp edge square cylinder in the wake of a circular cylinder

The response of an oscillating circular cylinder at the wake of an upstream fixed circular cylinder was classified by different researchers as galloping, wake induced galloping or wake induced vibration. Furthermore it is already known that a sharp edge square cylinder would undergo galloping if it is subjected to uniform flow. In this study the influence of the wake of a fixed circular cylinder on the response of a downstream square cylinder at different spacing ratios (S/D = 4, 8, 11) is experimentally investigated. The subject appears not to have received previous attention. The lateral displacements, lift forces and the pressure data from gauges mounted in the wake of the oscillating cylinder are recorded and analyzed. The single degree of freedom vibrating system has a low mass-damping parameter and the Reynolds number ranges from 7.7 × 10² to 3.7 × 10⁴.In contrast to that for two circular cylinders in tandem arrangement, the freely mounted downstream square cylinder displays a VIV type of response at all spacing ratios tested. There is no sign of galloping or wake induced galloping with the square cylinder. With increase at the spacing ratio the cross-flow oscillations decrease. It is shown that the vortices arriving from the upstream fixed circular cylinder play a major role on the shedding mechanism behind the downstream square cylinder and cause the square cylinder to shed vortices with frequencies above Strouhal frequency of the fixed square cylinder (St = 0.13). The VIV type of oscillations in the downstream square cylinder is most probably caused by the vortices newly generated behind the square cylinder.     

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

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

Behavior of vortex-induced vibration of a circular cylinder near a deformable wall with two degrees of freedom in steady flow

The behavior of vortex-induced vibration of a two-degree-of-freedom cylinder near a deformable wall in steady flow is investigated experimentally.The typical phenomenon of the two-degree-of-freedom cylinder’s VIV is discussed.The influences of initial gap between the cylinder and the wall on the dynamic responses of the cylinder are analyzed.The comparison is made about dynamic responses of the cylinder with one and two degrees of freedom.Experimental results show that the vibration of the cylinder near a deformable wall with a small value of initial gap-to-diameter ratios can generally be divided into two phases.The initial gap-to-diameter ratios have a noticeable influence on the occurrence of transverse vibration.The transverse maximum amplitude of the cylinder with two degrees of freedom is larger than that of the cylinder with one degree of freedom under the condition with the same values of other parameters.However,the vibration frequency of the cylinder for the two degrees of freedom case is smaller than that for the one degree of freedom case at the same value of Vr number.     

Linear irregular wave generation in a numerical wave tank

In the design of any floating or fixed marine structure, it is vital to test models in order to understand the fluid/structure interaction involved. A relatively inexpensive method, compared to physical model testing, of achieving this is to numerically model the structure and the wave conditions in a numerical wave tank. In this paper, a methodology for accurately replicating measured ocean waves in a numerical model at full scale is detailed. A Fourier analysis of the measured record allows the wave to be defined as a summation of linear waves and, therefore, Airy's linear wave theory may be used to input the wave elevation and associated water particle velocities. Furthermore, a structure is introduced into the model to display the ability of the model to accurately predict wave–structure interaction. A case study of three individual measured waves, which are recorded at the Atlantic marine energy test site, off the west coast of Ireland, is also presented. The accuracy of the model to replicate the measured waves and perform wave–structure interaction is found to be very high. Additionally, the absolute water particle velocity profile below the wave from the numerical model is compared to a filtered analytical approximation of the measured wave at a number of time-steps and is in very good agreement.     

Experimental study and numerical simulation on vortex-induced vibration of flexible riser

In order to study the effect of internal flow on vortex-induced vibration of flexible riser, the experiment on the vortex-induced vibration of flexible riser transporting fluid in the current was conducted in the physical oceanography laboratory of Ocean University of China. Considering the internal flowing fluid and external marine environment, the dynamic response of the flexible riser was measured. The corresponding numerical simulation was performed using the wake oscillatory model considering the extensibility of the riser system. Both the experiment and the numerical simulation indicated that with the increase of internal flow speed, the response amplitude increases, while the response frequency decreases.     

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

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

不同浸没长度下串列双立管涡激振动数值模拟

深海柔性立管涡激振动问题广受关注,其中多立管涡激振动较传统单立管涡激振动更为复杂。使用基于Open FOAM自主开发的深海柔性立管流固耦合求解器viv-FOAM-SJTU,对阶梯流中串列双立管涡激振动展开研究。首先针对阶梯流中串列双立管模型试验工况开展验证计算,下游立管的位移响应以及频率响应与试验结果基本吻合;而后改变立管浸没长度,分别选取浸没长度为0.4L(工况1)、0.55L(工况2)以及0.7L(工况3)的三种工况进行数值计算,其中L为立管长度。数值结果表明:3种浸没长度下,上下游立管的横流向振动均为一阶模态,且工况2和工况3中上下游立管的横流向位移均方根均大于工况1;上下游立管的顺流向主振模态在工况2和工况3中均为二阶,而在工况1中为一阶,且下游立管在工况2和工况3中的顺流向振动均表现出明显的多模态振动特性。     

不同固有频率比条件下圆柱双自由度涡激振动特性的数值研究

采用RANS方法,结合SST k-ω湍流模型,对不同顺流向与横顺流向固有频率的比值(即固有频率比,fnx/fny)条件下低质量比圆柱体的双自由度涡激振动进行了二维数值模拟。圆柱体的质量比为2.6,雷诺数范围为2 500~18 750,相应的约化速度范围为2~15,包括了经典试验中出现的整个锁定范围。通过研究发现,固有频率比是影响振动特性的重要参数,随着固有频率比的增加,响应幅值逐渐降低且向更高的约化速度偏移;在低约化速度范围内,固有频率比对顺流向和横流向振动之间的相位差以及升力频率有较大影响,从而得到各种不同偏向的8字形轨迹;最后对不同固有频率比条件下的尾涡模式进行了讨论,给出了对应不同约化速度时的尾涡模式。     

张力腿平台涡激运动模型试验与数值分析

为揭示张力腿平台涡激运动响应规律,采用数值模拟与模型试验相结合的方法研究了全水深系泊张力腿平台涡激运动响应。根据张力腿平台主尺度参数,按照几何相似制作了水池试验模型,在满足运动相似和动力相似的条件下开展了均匀流、剖面流模型试验,并将模型试验结果与数值模拟结果进行了对比,验证了数值模拟与模型试验结果的一致性。分析结果表明在均匀流作用下张力腿平台涡激运动的锁定区间在5.5<U_r<8.5,来流角对涡激运动影响较大。剖面流作用下平台涡激运动规律与均匀流作用基本一致,但剖面流造流引起的能量分散,使平台在XY平面的运动轨迹变得不规律,系泊系统提供的回复力对涡激运动幅值有抑制作用,来流角和流速对张力腿或立管模态影响明显。论文得到结论对于张力腿平台的工程设计有一定借鉴意义。     

Numerical simulation of a short flexible pipe subject to forced motion and vortex-induced vibration

A series of numerical simulations about a small scale(aspect ratio:63.2) flexible pipe undergoing forced harmonious oscillation and vortex-induced vibration(VIV) have been taken into account.The wake hydrodynamics and pipe deformation were accomplished by ANSYS MFX solution strategy designed for fluid-structure interaction(FSI) problem with well-performed LES model.The configuration of structured mesh,multi-domain design,different mesh stiffness admeasured by User Fortran ensured that the numerical task was competent to deal with large deformation related to this case.The introduction of instantaneous amplitude definition and modeless component decomposition method(Chen and Kim,2008) was helpful to reveal much more information from modal analysis.Most results from numerical simulation are generally consistent with those from model test(Choi and Hong,2000) via the comparison between them.As supplementary to model test,visualization of the vortex wake was also provided.It has been proved that the forced oscillation doesn't only excite a complicated dumbbell-like wake pattern around the outer thimble,but also results in inner flow inside the PVC pipe.The velocity of the inner flow increases with the frequency of forced oscillation.     

Study on Pure IL VIV of A Free Spanning Pipeline Under General Boundary Conditions

Pipeline spans may occur due to natural seabed irregularities or local scour of bed sediment. The pure in-line (IL) vortex-induced vibrations (VIV) analysis of the free spans is an important subject for design of pipeline in uneven seabed. The main objective of this paper is to analyze the characteristics of pure IL VIV of a free spanning pipeline under general boundary conditions. An IL wake oscillator model which can describe the coupling of pipeline structure and fluctuating drag is introduced and employed. The coupled partial differential equations of structure and wake are transformed into a set of ordinary differential equations using two-mode Galerkin method. Some case studies are presented and thoroughly discussed in order to investigate the effects of internal fluid, axial force and boundary conditions on the pure IL VIV.     

Numerical and Physical Investigation on Vortex-Induced Vibrations of Marine Risers

1 .Introduction Risers ,in deep water and strong current environments ,are prone to vortex-induced vibrations(VIV) .Vortex-inducedforces may excite the risersintheir normal modes of transverse direction.Un-der the“lock-in”condition,large resonant oscill…     

均匀流中深水系泊Truss Spar平台涡激运动试验研究

在工程中广泛应用的Truss Spar平台自诞生之初就为涡激运动问题所困扰,对其运动特性及抑制方法的机理性研究始终都在进行。采用截断系泊系统水池模型试验对设计中的Truss Spar平台在均匀流中的涡激运动响应进行了研究,分析了有无侧板平台的运动轨迹以及不同速度及流向角来流对平台运动的平衡位置的影响。在影响平台涡激运动特性的重要因素中,选取折合速度及流向角分别加以研究,分析了幅值在锁定区间中的变化规律,验证了减涡侧板对涡激运动良好的抑制作用,并对热点问题进行了探讨。为进一步研究Spar平台的涡激运动特性奠定了基础。     

Analysis and prediction of vortex-induced vibrations of variable-tension vertical risers in linearly sheared currents

Many studies have tackled the problem of vortex-induced vibrations (VIV) of a vertical riser with a constant tension and placed in uniform currents. In this study, attention is focused on the cross-flow VIV modelling, time-domain analysis and prediction of variable-tension vertical risers in linearly sheared currents. The partial-differential equation governing the riser transverse motion is based on a flexural tensioned-beam model with typical pinned-pinned supports. The hydrodynamic excitation model describing the modulation of lift force is based on a distributed van der Pol wake oscillator whose nonlinear equation is also partial-differential due to the implementation of a diffusion term. The variation of empirical wake coefficients with system parameters and the water depth-dependent Reynolds number is introduced. Based on the assumed Fourier mode shape functions obtained by accounting for the effect of non-uniform tension, the Galerkin technique is utilized to construct a low-dimensional multi-mode model governing the coupled fluid-riser interaction system due to VIV. Numerical simulations in the case of varying sheared flow profiles are carried out to systematically evaluate riser nonlinear dynamics and highlight the influence of fluid-structure parameters along with associated VIV aspects. In particular, the effects of shear and tensioned-beam (tension versus bending) parameters are underlined. Some comparisons with published experimental results and observations are qualitatively and quantitatively discussed. Overall parametric analysis and prediction results may be worthwhile for being a new benchmark against future experimental testing and/or numerical results predicted by an alternative model and methodology.     

Wave–current interaction with a vertical square cylinder

A three-dimensional numerical model is developed in this study to investigate the problem of wave–current–body interaction. The model solves the spatially averaged Navier–Stokes equations. Turbulence effects are modeled by a subgrid-scale (SGS) model using the concept of large eddy simulation (LES). The model is employed to study the wave–current interaction with a square cylinder that is mounted on the bottom and vertically pierces the free surface. The force analysis demonstrates that the presence of waves can reduce both the strength and frequency of vortex shedding induced by a uniform current due to the nonlinear wave–current interaction. The free surface elevation, strain rates of the mean flow, and eddy viscosity are found to closely correlate with the mechanism of vortex shedding. It is also shown that when the vortex shedding is neglected in the calculation such as by the potential flow approach, one may significantly underestimate the magnitude of in-line force. The energy spectral analysis reveals that there exist initiating, growing, and decaying regions for shedding vortices around the cylinder. In the vortex initiating region, both coherent and turbulent structures are nearly two-dimensional that become three-dimensional in the vortex growing region. The kinetic energy of both coherent and turbulent motions is dissipated in the vortex decaying region, within which the mean flow gradually returns back to two-dimensional.     

Cross-flow vortex-induced vibration reduction of a long flexible cylinder using 3 and 4 control rods with different configurations

Control rod is one of the common passive control methods to suppress the vortex-induced vibration (VIV) of cylindrical structures. In this paper, the experimental study is conducted to detailed understand the performance of multiple control rods in suppressing the cross-flow (CF) VIV for a long flexible cylinder. The influence of the spatial arrangement of 3 and 4 control rods on CF VIV response of the main cylinder is investigated in a towing tank. It is observed that the attack angle θ is a very significant parameter to affect the vibration response, dominant frequency and the VIV suppression efficiency of the main cylinder. Based on the suppression efficiencies analysis of VIV response in the present experimental investigation, the spatial arrangement of 3 control rods with θ = 40° and 4 control rods with θ = 30° is the best choice for suppressing the CF VIV response of the main flexible cylinder. Overall, the use of 4 control rods could reduce VIV more effectively than the application of 3 control rods.     

Experimental study of vortex-induced vibrations of a pipeline near an erodible sandy seabed

Based on similarity analyses, a series of experiments have been conducted with a newly established hydro-elastic facility to investigate the transverse vortex-induced vibrations (VIVs) of a submarine pipeline near an erodible sandy seabed under the influence of ocean currents. Typical characteristics of coupling processes between pipe vibration and soil scour in the currents have been summarized for Case I: pipe is laid above seabed and Case II: pipe is partially embedded in seabed on the basis of the experimental observations. Pipe vibration and the corresponding local scour are usually two coupled physical processes leading to an equilibrium state. The influence of initial gap-to-diameter ratio (e₀/D) on the interaction between pipe vibration and local scour has been studied. Experimental results show that the critical values of V_r for the initiation of VIVs of the pipe near an erodible sand bed get bigger with decreasing initial gap-to-diameter ratio within the examined range of e₀/D (−0.25<e₀/D<0.75). The comparison of the pipe vibrations near an erodible soil with those near a rigid boundary and under wall-free conditions indicates that the vibration amplitudes of the pipe near an erodible sand bed are close to the curve fit under wall-free conditions; nevertheless, for the same stability parameter, the maximum amplitudes for the VIV coupled with local scour increase with the increase of initial gap-to-diameter ratio.     

细长输流管内外流耦合振动特性研究

管内流动会影响输流管的振动响应,目前关于输流弹性管涡激振动方面的研究较少。基于计算流体力学(CFD)方法,开展内外流对细长输流弹性管振动特性影响的研究。首先在不考虑内流的情况下将弹性管涡激振动数值预报结果与模型试验数据进行对比,验证了数值方法的可靠性。再者考虑内外流耦合作用情况下,对不同内流流速下细长输流弹性管振动位移时—空分布、顺流向最大平均偏移、振动轨迹、内部横向涡的形成与分布等进行了对比分析。结果发现,与外流流速相比,内流流速的增加虽然难以改变弹性管的主振模态,但对沿管体的振动强度影响显著。顺流向最大偏移处管体运动轨迹发生明显的变形和跳跃。在剪切外流和均匀内流对弹性管的联合作用下,沿管跨方向模态间能量转换频繁,伴随着间歇性出现或消失的沿弹性管传播的行波组分,这主要归因于复杂的双重流固耦合系统(外流—管体,内流—管体)。在内流以附加质量力、离心力和科氏力形式的激励下,弹性管内二次流现象明显。在振动过程中,内部横向涡沿管壁生成、脱落并逐渐散布于整个横截面。