Shear flow induced vibrations of long slender cylinders with a wake oscillator model

A time domain model is presented to study the vibrations of long slender cylinders placed in shear flow. Long slender cylinders such as risers and tension legs are widely used in the field of ocean engineering. They are subjected to vortex-induced vibrations（VIV） when placed within a transverse incident flow. A three dimensional model coupled with wake oscillators is formulated to describe the response of the slender cylinder in cross-flow and in-line directions. The wake oscillators are distributed along the cylinder and the vortex-shedding frequency is derived from the local current velocity. A non-linear fiuid force model is accounted for the coupled effect between cross-flow and in-line vibrations. The comparisons with the published experimental data show that the dynamic features of VIV of long slender cylinder placed in shear flow can be obtained by the proposed model,such as the spanwise average displacement,vibration frequency,dominant mode and the combination of standing and traveling waves. The simulation in a uniform flow is also conducted and the result is compared with the case of nonuniform flow. It is concluded that the flow shear characteristic has significantly changed the cylinder vibration behavior.

A Novel Wake Oscillator Model for Vortex-Induced Vibrations Prediction of A Cylinder Considering the Influence of Reynolds Number

It is well known that the Reynolds number has a significant effect on the vortex-induced vibrations(VIV) of cylinders. In this paper, a novel in-line(IL) and cross-flow(CF) coupling VIV prediction model for circular cylinders has been proposed, in which the influence of the Reynolds number was comprehensively considered. The Strouhal number linked with the vortex shedding frequency was calculated through a function of the Reynolds number. The coefficient of the mean drag force was fitted as a new piecewise function of the Reynolds number, and its amplification resulted from the CF VIV was also taken into account. The oscillating drag and lift forces were modelled with classical van der Pol wake oscillators and their empirical parameters were determined based on the lock-in boundaries and the peak-amplitude formulas. A new peak-amplitude formula for the IL VIV was developed under the resonance condition with respect to the mass-damping ratio and the Reynolds number. When compared with the results from the experiments and some other prediction models, the present model could give good estimations on the vibration amplitudes and frequencies of the VIV both for elastically-mounted rigid and long flexible cylinders. The present model considering the influence of the Reynolds number could generally provide better results than that neglecting the effect of the Reynolds number.

A NEW WAKE OSCILLATOR MODEL FOR PREDICTING VORTEX INDUCED VIBRATION OF A CIRCULAR CYLINDER

This article proposes a new wake oscillator model for vortex induced vibrations of an elastically supported rigid circular cylinder in a uniform current. The near wake dynamics related with the fluctuating nature of vortex shedding is modeled based on the classical van der Pol equation, combined with the equation for the oscillatory motion of the body. An appropriate approach is developed to estimate the empirical parameters in the wake oscillator model. The present predicted results are compared to the experimental data and previous wake oscillator model results. Good agreement with experimental results is found.

NONLINEAR FLUID DAMPING IN STRUCTURE-WAKE OSCILLATORS IN MODELING VORTEX-INDUCED VIBRATIONS

A Nonlinear Fluid Damping （NFD） in the form of the square-velocity is applied in the response analysis of Vortex-Induced Vibrations （VIV）. Its nonlinear hydrodynamic effects on the coupled wake and structure oscillators are investigated. A comparison between the coupled systems with the linear and nonlinear fluid dampings and experiments shows that the NFD model can well describe response characteristics, such as the amplification of body displacement at lock-in and frequency lock-in, both at high and low mass ratios. Particularly, the predicted peak amplitude of the body in the Griffin plot is in good agreement with experimental data and empirical equation, indicating the significant effect of the NFD on the structure motion.

Application of time–frequency entropy from wake oscillation to gas–liquid flow pattern identification

Gas–liquid two-phase flow abounds in industrial processes and facilities. Identification of its flow pattern plays an essential role in the field of multiphase flow measurement. A bluff body was introduced in this study to recognize gas–liquid flow patterns by inducing fluid oscillation that enlarged differences between each flow pattern. Experiments with air–water mixtures were carried out in horizontal pipelines at ambient temperature and atmospheric pressure. Differential pressure signals from the bluff-body wake were obtained in bubble, bubble/plug transitional, plug, slug, and annular flows. Utilizing the adaptive ensemble empirical mode decomposition method and the Hilbert transform, the time–frequency entropy S of the differential pressure signals was obtained. By combining S and other flow parameters, such as the volumetric void fraction β, the dryness x, the ratio of density φ and the modified fluid coefficient ψ, a new flow pattern map was constructed which adopted S（1–x）φ and （1–β）ψ as the vertical and horizontal coordinates, respectively. The overall rate of classification of the map was verified to be 92.9% by the experimental data. It provides an effective and simple solution to the gas–liquid flow pattern identification problems.

Coupled Cross-Flow/in-Line Vortex-Induced Vibration Responses of a Catenary-Type Riser Subjected to Uniform Flows

A three-dimensional numerical scheme was developed to investigate the vortex-induced vibration(VIV)of a catenary-type riser(CTR)in the in-line(IL)and cross-flow(CF)directions.By using the vector form intrinsic finite element method,the CTR was discretized into a finite number of spatial particles whose motions satisfy Newton’s second law.The Van der Pol oscillator was used to simulate the effect of vortex shedding.The coupling equations of structural vibration and wake oscillator were solved using an explicit central differential algorithm.The numerical model was verified with the published results.The VIV characteristics of the CTR subjected to uniform flows,including displacement,frequency,standing wave,traveling wave,motion trajectory,and energy transfer,were studied comprehensively.The numerical results revealed that the multimode property occurs in the CF-and IL-direction VIV responses of the CTR.An increase in the flow velocity has slight effects on the maximum VIV displacement.Due to structural nonlin-earity,the double-frequency relationship in the CF and IL directions is rarely captured.Therefore,the vibration trajectories display the shape of an inclined elliptical orbit.Moreover,the negative energy region is inconspicuous under the excitation of the uniform flow.

Nonlinear Coupled in-Line and Cross-Flow Vortex-Induced Vibration Analysis of Top Tensioned Riser

The fluctuating furces of the fluid exerted on the top terrsioned riser （＇FIR） in the in-line and cross-flow directions are both modeled by van del Pol wake oscillator model and the nonlinear coupled dynamics of the in-line and cross-flow vortex-induced vibrations （VIV） of the riser are analyzed in time domain in this papar. The numencal shnulation results of the riser＇s in-line and cross-flow displacements and curvatures are compared with experimental measurements and the comparison shows the validity of this method in modeling some main features of the riser＇s VIV. Finally, the effects of the riser＇s top tensions and internal flow velocities on the coupled vibrations of the riser are investigated.

Nonlinear Numerical Analysis of Vortex-Induced Vibration of A Three-Dimensional Deepwater Steep Wave Riser with Large Deformation Features

The cross-flow(CF)vortex-induced vibration(VIV)of a deepwater steep wave riser(SWR)subjected to uniform or shear flow loads is investigated numerically.The model is based on a three-dimensional(3D)nonlinear elastic rod theory coupled with a wake oscillator model.In this numerical simulation,the nonlinear motion equations of the riser with large deformation features are established in a global coordinate system to avoid the transformation between global and local coordinate systems,and are discretized with the time-domain finite element method(FEM).A wakeoscillator model is employed to study the vortex shedding,and the lift force generated by the wake flow is described in a van der Pol equation.A Newmark-βiterative scheme is used to solve their coupling equation for the VIV response of the SWR.The developed model is validated against the existing experimental results for the VIV response of the top-tension riser(TTR).Then,the numerical simulations are executed to determine VIV characteristics of the SWR.The effects of both flow velocity and the spanwise length of the flow field on the drag coefficient in the inline(IL)direction and the lift coefficient in the CF direction are investigated systematically.The results indicate that compared with TTR,the low frequency and multi-modal vibration are the main components of the SWR due to the large deformation and flexible characteristics.For shear flow,the multi-frequency resonance dominates the VIV response of the SWR,especially at the hang-off segment.

Prediction of Streamwise Flow-Induced Vibration of A Circular Cylinder in the First Instability Range

The streamwise flow-induced vibration of a circular cylinder with symmetric vortex shedding in the first instability range is investigated, and a wake oscillator model for the dynamic response prediction is proposed. An approach is applied to calibrate the empirical parameters in the present model; the numerical and experimental results are compared to validate the proposed model. It can be found that the present prediction model is accurate and sufficiently simple to be easily applied in practice.

Random vortex induced vibration response of suspended flexible cable to fluctuating wind

A popular dynamical model for the vortex induced vibration(VIV)of a suspended flexible cable consists of two coupled equations.The first equation is a partial differential equation governing the cable vibration.The second equation is a wake oscillator that models the lift coefficient acting on the cable.The incoming wind acting on the cable is usually assumed as the uniform wind with a constant velocity,which makes the VIV model be a deterministic one.In the real world,however,the wind velocity is randomly fluctuant and makes the VIV of a suspended flexible cable be treated as a random vibration.In the present paper,the deterministic VIV model of a suspended flexible cable is modified to a random one by introducing the fluctuating wind.Using the normal mode approach,the random VIV system is transformed into an infinite-dimensional modal vibration system.Depending on whether a modal frequency is close to the aeolian frequency or not,the corresponding modal vibration is characterized as a resonant vibration or a non-resonant vibration.By applying the stochastic averaging method of quasi Hamiltonian systems,the response of modal vibrations in the case of resonance or non-resonance can be analytically predicted.Then,the random VIV response of the whole cable can be approximately calculated by superimposing the response of the most influential modal vibrations.Some numerical simulation results confirm the obtained analytical results.It is found that the intensity of the resonant modal vibration is much higher than that of the non-resonant modal vibration.Thus,the analytical results of the resonant modal vibration can be used as a rough estimation for the whole response of a cable.

Nonlinear Analysis of Bidirectional Vortex-Induced Vibration of A Deepwater Steep Wave Riser Subjected to Oblique Currents

An improved three-dimensional(3D)time-domain couple model is established in this paper to simulate the bidirectional vortex-induced vibration(VIV)of a deepwater steep wave riser(SWR)subjected to oblique currents.In this model,the nonlinear motion equations of the riser are established in the global coordinate system based on the slender rod theory with the finite element method.Van der Pol equations are used to describe the lift forces induced by the x-and y-direction current components,respectively.The coupled equations at each time step are solved by a Newmark-βiterative scheme for the SWR VIV.The present model is verified by comparison with the published experimental results for a top-tension riser.Then,a series of simulations are executed to determine the influences of the oblique angle/velocity of the current,different top-end positions and the length of the buoyancy segment on the VIV displacement,oscillating frequency as well as hydrodynamic coefficients of the SWR.The results demonstrate that there exists a coupled resonant VIV corresponding to x-direction and y-direction,respectively.However,the effective frequency is almost identical between the vibrations at the hang-off segment along x and y directions.The addition of the buoyancy modules in the middle of the SWR has a beneficial impact on the lift force of three segments and simultaneously limits the VIV response,especially at the decline segment and the hang-off segments.Additionally,the incident current direction significantly affects the motion trajectory of the SWR which mainly includes the fusiform and rectangle shapes.

PREDICTION OF VORTEX-INDUCED VIBRATION OF LONG FLEXIBLE CYLINDERS MODELED BY A COUPLED NONLINEAR OSCILLATOR:INTEGRAL TRANSFORM SOLUTION

The Generalized Integral Transform Technique （GITT） was applied to predict dynamic response of Vortex-Induced Vibration （VIV） of a long flexible cylinder. A nonlinear wake oscillator model was used to represent the cross-flow force acting on the cylinder, leading to a coupled system of second-order Partial Differential Equations （PDEs） in temporal variable. The GITT approach was used to transform the system of PDEs to a system of Ordinary Differential Equations （ODEs）, which was numerically solved by using the Adams-Moulton and Gear method （DIVPAG） developed by the International Mathematics and Statistics Library （IMSL）. Numerical results were presented for comparison to those given by the finite difference method and experimental results, allowing a critical evaluation of the technique performance. The influence of variation of mean axial tension induced by elongation of flexible cylinder was evaluated, which was shown to be not negligible in numerical simulation of VIV of a long flexible cylinder.

Numerical Investigation of the Superimposed Effects on Stator Wake Oscillation in an Axial-radial Combined Compressor

The superimposed influences of the blade rows in a multistage compressor are important because different matches of upstream and downstream blades can result in significant differences in the stator wake oscillation. Numerical investigation of the axial stator wake oscillation, which is affected upstream by the axial rotor and downstream by the radial rotor, was performed in an axial-radial combined compressor. Many configurations with different blade numbers and locations, which influence axial stator wake oscillation were investigated. When rotors have equal blade numbers, the axial stator wake oscillates periodically versus time within time T(moving blade passing 1/3 revolution). In contrast, stator wake oscillates irregularly within T when rotors have different blade numbers. A model-split subtraction method is presented in order to separate the influences of the individual blade rows on the wake oscillation of the axial stator. Analysis from the rotor-stator configuration showed that the unsteady flow angle fluctuation response is caused by the upstream rotor. For the rotor-stator-rotor configuration, the unsteady flow angle fluctuations are influenced by upand downstream blade rows. With the model-split subtraction method, the upand downstream influences on the flow angle fluctuation could be clearly separated and quantified. Low amplitudes could be observed when the influences from upand downstream moving rows were superimposed with the "positive peaknegative peak" type wave. Clocking investigations were carried out to change the relative superimposed phase of influences from the surrounding blade rows in order to modulate the amplitudes of the axial stator wake oscillation. However, the amplitudes did not reach the maximum when they were superimposed with "positive peak-positive peak" type wave due to the impact of the interaction between the two moving blade rows.

Unsteady Flow Variability Driven by Rotor-stator Interaction at Rotor Exit

Numerical investigation of the unsteady flow variability driven by rotorstator interaction in a transonic axial compressor is performed. Two models with close and far axial gap between rotor and stator rows are studied in the simulation. Particular attention is attached to the analysis of mechanisms involved in driving rotor wake oscillation, rotor wake skewing and flow angle fluctuation at rotor exit. The results show that smaller axial gap is favorable to enhance the interaction in the region between two adjacent rows, and the fluctuation of the static pressure difference between two sides of rotor wake is improved by potential field from down stator, which is the driving force for rotor wake oscillation. The interaction between rotor and stator is weakened by increasing axial distance, rotor wake shifts to suction side of rotor blade with 5%-10% of rotor pitch, the absolute value of flow angle at rotor exit is less than that in the case of close interspace for every time step, and the fluctuation amplitude is also decreased.