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 GEOMETRIC THEORY IN INVESTIGATION ON MODE LOCALIZATION AND FREQUENCY LOCI VEERING PHENOMENA

According to the mapping theory in complex plane,the geometric features of eigen fre-quency loci of systems undergoing free vibrations are investigated.It is concluded that the phenomenaof curve coalescence and veering can be described in a unified manner from the singularities of map-ping from the complex parameter plane onto the complex frequency plane.The formation of a branchpoint in the parameter space is the foundation of explaining localization and veering phenornena.Bythe use of condensation to reduce the dimension of a system,the scope of application of the geometrictheory is widely expanded.The theory is applied to examples to verify the validity of the proposed ap-proach.The present work is an improvement and extension of recent work by M.S.Traintafyllou etal.,     

ELECTRICALLY FORCED THICKNESS-SHEAR VIBRATIONS OF QUARTZ PLATE WITH NONLINEAR COUPLING TO EXTENSION

We study electrically forced nonlinear thickness-shear vibrations of a quartz plate resonator with relatively large amplitude. It is shown that thickness-shear is nonlinearly coupled to extension due to the well-known Poynting effect in nonlinear elasticity. This coupling is relatively strong when the resonant frequency of the extensional mode is about twice the resonant frequency of the thickness-shear mode. This happens when the plate length/thickness ratio assumes certain values. With this nonlinear coupling, the thickness-shear motion is no longer sinusoidal. Coupling to extension also affects energy trapping which is related to device mounting. When damping is 0.01, nonlinear coupling causes a frequency shift of the order of 10^-6 which is not insignificant,and an amplitude change of the order of 10^-8. The effects are expected to be stronger under real damping of 10^-5 or larger. To avoid nonlinear coupling to extension, certain values of the aspect ratio of the plate should be avoided.     

Effects of tension on vortex-induced vibration(ⅥⅤ) responses of a long tensioned cylinder in uniform flows

The effects of tension on vortex-induced vibration(ⅥⅤ) responses for a tension-dominated long cylinder with an aspect ratio of 550 in uniform flows are experimentally investigated in this paper. The results show that elevated tension suppresses fluctuations of maximum displacement with respect to flow velocity and makes chaotic ⅥⅤ more likely to appear. With respect to periodic ⅥⅤ, if elevated tension is applied, the dominant vibration frequency in the in-line(IL) direction will switch from a fundamental vibration frequency to twice the value of the fundamental vibration frequency, which results in a ratio of the dominant vibration frequency in the IL direction to that in the cross-flow direction of 2.0. The suppression of the elevated tension in the fluctuation of the maximum displacement causes the axial tension to become an active control parameter for the ⅥⅤ maximum displacement of a tension-dominated long riser or tether of an engineering structure in deep oceans. However,the axial tension must be optimized before being used since the high dominant vibration frequency due to the elevated tension may unfavorably affect the fatigue life of the riser or tether.     

Scale-adaptive simulation of flow past wavy cylinders at a subcritical Reynolds number

The Xu & Yan scale-adaptive simulation (XYSAS) model is employed to simulate the flows past wavy cylinders at Reynolds number 8 × 10 3.This approach yields results in good agreement with experimental measurements.The mean flow field and near wake vortex structure are replicated and compared with that of a corresponding circular cylinder.The effects of wavelength ratios λ/D m from 3 to 7,together with the amplitude ratios a /D m of 0.091 and 0.25,are fully investigated.Owing to the wavy configuration,a maximum reduction of Strouhal number and root-meansquare (r.m.s) fluctuating lift coefficients are up to 50% and 92%,respectively,which means the vortex induced vibration (VIV) could be effectively alleviated at certain larger values of λ/D m and a /D m.Also,the drag coefficients can be reduced by 30%.It is found that the flow field presents contrary patterns with the increase of λ/D m.The free shear layer becomes much more stable and rolls up into mature vortex only further downstream when λ/D m falls in the range of 5-7.The amplitude ratio a /D m greatly changes the separation line,and subsequently influences the wake structures.     

Effects of tension on vortex-induced vibration (VIV) responses of a long tensioned cylinder in uniform flows

The effects of tension on vortex-induced vibra-tion (VIV) responses for a tension-dominated long cylinder with an aspect ratio of 550 in uniform flows are experimen-tally investigated in this paper. The results show that elevated tension suppresses fluctuations of maximum displacement with respect to flow velocity and makes chaotic VIV more likely to appear. With respect to periodic VIV, if elevated tension is applied, the dominant vibration frequency in the in-line (IL) direction will switch from a fundamental vibra-tion frequency to twice the value of the fundamental vibration frequency, which results in a ratio of the dominant vibration frequency in the IL direction to that in the cross-flow direc-tion of 2.0. The suppression of the elevated tension in the fluctuation of the maximum displacement causes the axial tension to become an active control parameter for the VIV maximum displacement of a tension-dominated long riser or tether of an engineering structure in deep oceans. However, the axial tension must be optimized before being used since the high dominant vibration frequency due to the elevated tension may unfavorably affect the fatigue life of the riser or tether.     

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.     

An investigation into natural vibrations of fluid–structure interaction systems subject to Sommerfeld radiation condition

A fluid-structure interaction system subject to Sommerfeld＇s condition is defined as a Sommerfeld system which is divided into three categories： Fluid Sommerfeld （FS） System, Solid Sommerfeld （SS） System and Fluid Solid Sommerfeld （FSS） System of which Sommerfeld conditions are imposed on a fluid boundary only, a solid boundary only and both fluid and solid boundaries, respectively. This paper follows the previous initial results claimed by simple examples to further mathematically investigate the natural vibrations of generalized Sommerfeld systems. A new parameter representing the speed of radiation wave for generalized 3-D problems with more complicated boundary conditions is introduced into the Sommerfeld condition which allows investigation of the natural vibrations of a Sommerfeld system involving both free surface and compressible waves. The mathematical demonstrations and selected examples confirm and reveal the natural behaviour of generalized Sommerfeld systems defined above. These generalized conclusions can be used in theoretical or engineering analysis of the vibrations of various Sommerfeld systems in engineering.     

Regimes of vortex shedding from an in-line oscillating circular cylinder in the uniform flow

A finite volume method for the time dependent viscous incompressible flow around an in-line oscillating circular cylinder at Reynolds number of 200, 855 is presented in this paper. The Navier-Stokes equations in a finite volume form are solved with a moving grid system, based on a time dependent coordinate transformation. To investigate the vortex-shedding characteristics behind the circular cylinder and the effects of Reynolds number and other non-dimensional parameters such as reduced amplitude and reduced frequency， several numerical schemes have been tested with different amplitude and frequency close to Sto and a harmonic at each Reynolds number. Present numerical results indicate several types of vortex shedding mode which is known mainly depending on the reduced frequency and also the reduced amplitude, which is called synchronization or lock-on.     

Linear stability of a fluid channel with a porous layer in the center简

We perform a Poiseuille flow in a channel linear stability analysis of a inserted with one porous layer in the centre, and focus mainly on the effect of porous filling ratio. The spectral collocation technique is adopted to solve the coupled linear stability problem. We investigate the effect of permeability, σ, with fixed porous filling ratio ψ = 1/3 and then the effect of change in porous filling ratio. As shown in the paper, with increasing σ, almost each eigenvalue on the upper left branch has two subbranches at ψ = 1/3. The channel flow with one porous layer inserted at its middle （ψ = 1/3） is more stable than the structure of two porous layers at upper and bottom walls with the same parameters. By decreasing the filling ratio ψ, the modes on the upper left branch are almost in pairs and move in opposite directions, especially one of the two unstable modes moves back to a stable mode, while the other becomes more instable. It is concluded that there are at most two unstable modes with decreasing filling ratio ψ. By analyzing the relation between ψ and the maximum imaginary part of the streamwise phase speed, Cimax, we find that increasing Re has a destabilizing effect and the effect is more obvious for small Re, where ψ a remarkable drop in Cimax can be observed. The most unstable mode is more sensitive at small filling ratio ψ, and decreasing ψ can not always increase the linear stability. There is a maximum value of Cimax which appears at a small porous filling ratio when Re is larger than 2 000. And the value of filling ratio 0 corresponding to the maximum value of Cimax in the most unstable state is increased with in- creasing Re. There is a critical value of porous filling ratio （= 0.24） for Re = 500; the structure will become stable as ψ grows to surpass the threshold of 0.24; When porous filling ratio ψ increases from 0.4 to 0.6, there is hardly any changes in the values of Cimax. We have also observed that the critical Reynolds number is especially sensitive for small ψ where the fastest drop is observed, and there may be a wide range in which the porous filling ratio has less effect on the stability （ψ ranges from 0.2 to 0.6 at σ = 0.002）. At larger permeability, σ, the critical Reynolds number tends to converge no matter what the value of porous filling ratio is.     

Numerical simulation of two-degree-of-freedom vortex-induced vibration of a circular cylinder close to a plane boundary

Two-degree-of-freedom vortex-induced vibrations (VIV) of a circular cylinder close to a plane boundary are investigated numerically. The Reynolds-Averaged Navier-Stokes (RANS) equations are solved using the Arbitrary Lagrangian Eulerian (ALE) scheme with a k-ω turbulence model closure. The numerical model is validated against experimental data of VIV of a cylinder in uniform flow and VIV of a cylinder close to a plane boundary at low mass ratios. The numerical results of the vibration mode, vibration amplitude and frequency agree well with the experimental data. VIV of a circular cylinder close to a plane boundary is simulated with a mass ratio of 2.6 and gap ratios of e/D=0.002 and 0.3 (gap ratio is defined as the ratio of gap between the cylinder and the bed (e) to cylinder diameter (D)). Simulations are carried out for reduced velocities ranging from 1 to 15 and Reynolds numbers ranging from 1000 to 15 000. It is found that vortex-induced vibrations occur even if the initial gap ratio is as small as e/D=0.002, although reported research indicated that vortex shedding behind a fixed circular cylinder is suppressed at small gap ratios (e/D<0.3 or 0.2). It was also found that vibration amplitudes are dependant on the bouncing back coefficient when the cylinder hits the plane boundary. Three vortex shedding modes are identified according to the numerical results: (i) single-vortex mode where the vortices are only shed from the top of the cylinder; (ii) vortex-shedding-after-bounce-back mode; (iii) vortex-shedding-before-bounce-back mode. It was found that the vortex shedding mode depends on the reduced velocity.     

An immersed interface method for flow past circular cylinder in the vicinity of a plane moving wall

Two‐dimensional flows past a stationary circular cylinder near a plane boundary are numerically simulated using an immersed interface method with second‐order accuracy. Instead of a fixed wall, a moving wall with no‐slip boundary is considered to avoid the complex involvement of the boundary layer and to focus only on the shear‐free wall proximity effects for investigating the force dynamics and flow fields. To analyze the convergence and accuracy of our implementation, numerical studies have been first performed on a simple test problem of rotational flow, where the second order of convergence is confirmed through numerical experiments and an optimal range of relative grid‐match ratio of Lagrangian to Eulerian grid sizes has been recommended. By comparing the force quantities and the Strouhal number, the accuracy of this method has been demonstrated on the flow past a stationary isolated cylinder. The cylinder is then put in proximity to the wall to investigate the shear‐free wall proximity effects in the low Reynolds number regime (20≤Re≤200). The gap ratio, e/D, where e denotes the gap between the cylinder and the moving wall and D denotes the diameter of the cylinder, is taken from 0.10 to 2.00 to determine the critical gap ratio, (e/D)_critical, for the alternate vortex shedding, where the fluid forces, flow fields and the streamwise velocity profiles are studied. One of the key findings is that the (e/D)_critical for the alternate vortex shedding decreases as the Reynolds number increases. We also find that, in this low Reynolds number regime, the mean drag coefficient increases and peaks at e/D = 0.5 with the increase of e/D and keeps decreasing gently from e/D = 0.5 to e/D = 2.0, while the mean lift coefficient decreases monotonically with the increase of e/D. New correlations are then proposed for computing force coefficients as a function of Re and e/D for a cylinder in the vicinity of a moving plane wall. Copyright © 2015 John Wiley & Sons, Ltd.     

An experimental investigation of flow past a dual step cylinder

Flow development in the wake of a dual step cylinder has been investigated experimentally using Laser Doppler Velocimetry and flow visualization. The dual step cylinder model is comprised of a large diameter cylinder (D) mounted at the mid-span of a small diameter cylinder (d). The experiments have been performed for a Reynolds number (Re _D ) of 1,050, a diameter ratio (D/d) of 2, and a range of large cylinder aspect ratios (L/D). The results show that the flow development is highly dependent on L/D. The following four distinct flow regimes can be identified based on vortex dynamics in the wake of the large cylinder: (1) for L/D ≥ 15, three vortex shedding cells form in the wake of the large cylinder, one central cell bounded by two cells of lower frequency, (2) for 8 < L/D ≤ 14, a single vortex shedding cell forms in the wake of the large cylinder, (3) for 2 < L/D ≤ 6, vortex shedding from the large cylinder is highly three-dimensional. When spanwise vortices are shed, they deform substantially and attain a hairpin shape in the near wake, (4) for 0.2 ≤ L/D ≤ 1, the large cylinder induces vortex dislocations between small cylinder vortices. The results show that for Regimes I to III, on the average, the frequency of vortex shedding in the large cylinder wake decreases with L/D, which is accompanied by a decrease in coherence of the shed vortices. In Regime IV, small cylinder vortices connect across the large cylinder wake, but these connections are interrupted by vortex dislocations. With decreasing L/D, the frequency of dislocations decreases and the dominant frequency in the large cylinder wake increases toward the small cylinder shedding frequency.     

Flow characteristics around a circular cylinder subjected to vortex-induced vibration near a plane boundary

Although vortex-induced vibration (VIV) has been extensively studied, much of existing literature deals with uniform flow in the absence of a boundary. The VIV flow field of a structure close to a boundary generally remains unexplored, but it can have important engineering implications, such as pipeline scour if the boundary is an erodible seabed. In this paper, laboratory experiments are performed to investigate the flow characteristics of an elastically mounted circular cylinder undergoing VIV, and a rigid plane boundary is considered to simplify the problem. The initial gap-to-diameter ratio is fixed at 0.8, and six different reduced velocities are considered. The velocity field is measured using a high resolution particle image velocimetry (PIV) system, which has several advantages over traditional PIV systems, including high sampling rate and the ability to mitigate scatter of laser light near the boundary, allowing accurate measurements at the viscous sublayer. This paper presents the vibration amplitude and oscillation frequency for different V_r; in addition, the mean velocity field, turbulence characteristics, vortex behavior, gap flow velocity, and normal/shear stresses on the boundary were measured/calculated, leading to new insights on the flow field behavior.     

Shear layers of a circular cylinder with rotary oscillation

The behavior of the separated shear layers and the near wake of a circular cylinder with small-amplitude rotary oscillations (Ω₁ = 0.05−0.15 for f _f/f _o ≤ 1.25) were investigated experimentally at Re = 3,700. Measurements of an unforced cylinder were also made for 2,000 ≤ Re ≤ 10,000 to better understand the effects of rotary oscillations. The results show that the shear-layer vortices formed closer to the cylinder and the distance separating them was found to decrease with cylinder oscillations. The shear-layer frequency, however, increased with increasing forcing frequency f _f. The formation-region length l _f decreased significantly with increasing f _f while decreased to a lesser extent with increasing normalized oscillation amplitude Ω₁. The shear layer also diffused to a length L _d larger than that of an unforced cylinder, while the l _f-L _d-Strouhal frequency offsetting mechanism was generally maintained. The near wake was of lower momentum compared to an unforced cylinder, and the transverse velocity fluctuations associated with the unforced vortex-shedding frequency f _o always presented a local peak at f _f/f _o = 0.5, regardless of Ω₁ tested.     

Lock-on characteristics behind two side-by-side cylinders of diameter ratio two at small gap ratio

The lock-on characteristics, the detailed interactions and downstream evolutions of the wakes behind side-by-side cylinders of unequal diameter (D/d?=?2), spaced by a gap ratio 0.75 (G/D?=?0.75), are investigated at Reynolds number 600 by the dye flow visualization, laser Doppler anemometry (LDA) and particle image velocimeter (PIV) velocity measurements. The lock-on frequency bands are studied by LDA and PIV at Reynolds number 2,000. The D, d and G are the diameters of the large, the small cylinders and the net gap between two cylinders, respectively. Periodic excitations, in form of rotary oscillation about the cylinder center, are applied to the large cylinder with the same amplitude. It is found that while the large cylinder is excited, two lock-on frequency bands of the wake behind the large cylinder are detected. These two lock-on frequency bands correspond to the primary and the one-third sub-harmonic lock-on of the wake behind large cylinder, respectively. These two lock-on frequency bands distribute symmetrically about the fundamental and the third superharmonic of the natural shedding frequency behind a single cylinder at the same Reynolds number. The left-shifted frequency band (1.8?≤?f _e /f _os ?≤?2.0) is not considered as a locked-on frequency band because the phase difference between two excitation frequencies across f _e /f _os ?=?2.0 vary significantly. While the wake behind the large cylinder is locked-on at f _e /3 (or f _os ), the gap flow becomes unbiased and the frequency of the wake behind small cylinder remains around the natural shedding frequency. Thus, the frequency band of 3.0?≤?f _e /f _os ?≤?3.22 is also not locked-on because the phase difference in the narrow wake excited at f _e /f _os ?=?2.93 and 3.07 changes significantly. Note f _e and f _os denote the excitation frequency and the natural shedding frequency behind a single large cylinder, respectively.     

Lock-on characteristics of a cavity shear layer

This investigation focuses on defining the lock-on regions of a cavity shear layer subject to local periodic excitations. A circular cylinder of small diameter (d=4 mm), located very close to the upstream edge of cavity, is used to generate the local periodic excitations in the form of oscillatory rotation about its center with various angular amplitudes (Δθ) and frequencies (f_e). All the experiments were conducted in a recirculating water channel at three different Reynolds numbers that are based on the momentum thickness at the upstream edge of cavity (Re_θ0=152, 216 and 278). The LDV system and the laser-sheet technique are employed to perform the quantitative velocity measurements and the qualitative flow visualization, respectively. For cavity flows at three Reynolds numbers studied, the resonant lock-on is found to be the primary lock-on region within the range of frequency ratio (f_e/f₀=0.28–2.0). Here f₀ denotes the natural instability frequency of an unexcited cavity shear layer. The frequency bandwidth of resonant lock-on region does increase with increasing excitation amplitudes (Δθ). While the excitation amplitudes are smaller than 5° (Δθ5°), the resonant lock-on region, at Reynolds numbers 216 and 278, distributes asymmetrically about f_e/f₀=1.0 and biases to the high frequency (or large f_e/f₀) side. However, the sidewise expansion of resonant lock-on region is enlarged and the degree of asymmetric distribution is alleviated at large excitation amplitudes (Δθ>5°). The amount of sidewise expansion of the resonant lock-on region biased toward the high-frequency side is more significant at the lowest Reynolds number (152) than those at two higher Reynolds numbers (216 and 278). Besides, there exists a sub-harmonic lock-on region only at the lowest Reynolds number 152. The existence of a sub-harmonic lock-on region clearly reveals that the differential equation governing the self-excited oscillation within a cavity contains the quadratic nonlinear term. Further, at the lowest Reynolds number (152), the sidewise expansion of the sub-harmonic lock-on region is much narrower than that of the resonant lock-on region.     

Stress boundary layers in the viscoelastic flow past a cylinder in a channel： limiting solutions

The flow past a cylinder in a channel with the aspect ratio of 2:1 for the upper convected Maxwell (UCM) fluid and the Oldroyd-B fluid with the viscosity ratio of 0.59 is studied by using the Galerkin/Least-square finite element method and a p-adaptive refinement algorithm. A posteriori error estimation indicates that the stress-gradient error dominates the total error. As the Deborah number, D_e, approaches 0.8 for the UCM fluid and 0.9 for the Oldroyd-B fluid, strong stress boundary layers near the rear stagnation point are forming, which are characterized by jumps of the stress-profiles on the cylinder wall and plane of symmetry, huge stress gradients and rapid decay of the gradients across narrow thicknesses. The origin of the huge stress-gradients can be traced to the purely elongational flow behind the rear stagnation point, where the position at which the elongation rate is of 1/2D_e approaches the rear stagnation point as the Deborah number approaches the critical values. These observations imply that the cylinder problem for the UCM and Oldroyd-B fluids may have physical limiting Deborah numbers of 0.8 and 0.9, respectively.The project supported by the National Natural Science Foundation of China (50335010 and 20274041) and the MOLDFLOW Comp. Australia.     

Vortex dynamics of a cylinder wake in proximity to a wall

Large-eddy simulations (LES) are used to investigate the modifications of wake dynamics and turbulence characteristics behind a circular cylinder placed near a wall for varying gap-to-diameter (G/D) ratios (where G signifies the gap between the wall and the cylinder, and D the cylinder diameter). The three-dimensional (3-D), time-dependent, incompressible Navier–Stokes equations with a dynamic subgrid-scale model are solved using a symmetry-preserving finite-difference scheme of second-order spatial and temporal accuracy. The immersed boundary (IB) method is employed to impose the no-slip boundary condition on the cylinder surface. Flow visualizations along with turbulence statistics are presented to gain insight into the flow structures that are due to interaction between the shear layers and the approaching boundary layer. Apart from the vortex shedding mechanism, the paper illustrates the physics involving the shear layer transition, stretching, breakdown and turbulence generation, either qualitatively or quantitatively, in the presence of a wall for a Reynolds number of Re=1440 (based on D and the inlet free-stream velocity U_∞).     

Flow structure of wake behind a rotationally oscillating circular cylinder

Flow around a circular cylinder oscillating rotationally with a relatively high forcing frequency has been investigated experimentally. The dominant parameters affecting this experiment are the Reynolds number (Re), oscillation amplitude (θ_A), and frequency ratio F_R=f_f/f_n, where f_f is the forcing frequency and f_n is the natural frequency of vortex shedding. Experiments were carried out under conditions of Re=4.14×10³, 0°θ_A60° and 0.0F_R2.0. Rotational oscillation of the cylinder significantly modified the flow structure in the near-wake. Depending on the frequency ratio F_R, the cylinder wake showed five different flow regimes, each with a distinct wake structure. The vortex formation length and the vortex shedding frequency were greatly changed before and after the lock-on regime where vortices shed at the same frequency as the forcing frequency. The lock-on phenomenon always occurred at F_R=1.0 and the frequency range of the lock-on regime expanded with increasing oscillation amplitude θ_A. In addition, the drag coefficient was reduced when the frequency ratio F_R was less than 1.0 (F_R<1.0) while fixing the oscillation amplitude at θ_A=30°. When the oscillation amplitude θ_A was used as a control parameter at a fixed frequency ratio F_R=1.0 (lock-on regime), the drag reduction effect was observed at all oscillation amplitudes except for the case of θ_A=30°. This type of active flow control method can be used effectively in aerodynamic applications while optimizing the forcing parameters.