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.     

Flow-induced vibrations of a tethered circular cylinder

One of the most basic examples of fluid-structure interaction is provided by a tethered body in a fluid flow. The tendency of a tethered buoy to oscillate when excited by waves is a well-known phenomenon; however, it has only recently been found that a submerged buoy will act in a similar fashion when exposed to a uniform flow at moderate Reynolds numbers, with a transverse peak-to-peak amplitude of approximately two diameters over a wide range of velocities. This paper presents results for the related problem of two-dimensional simulations of the flow past a tethered cylinder. The coupled Navier–Stokes equations and the equations of motion of the cylinder are solved using a spectral-element method. The response of the tethered cylinder system was found to be strongly influenced by the mean layover angle as this parameter determined if the oscillations would be dominated by in-line oscillations, transverse oscillations or a combination of the two. Three branches of oscillation are noted, an in-line branch, a transition branch and a transverse branch. Within the transition branch, the cylinder oscillates at the shedding frequency and modulates the drag force such that the drag signal is dominated by the lift frequency. It is found that the mean amplitude response is greatest at high reduced velocities, i.e., when the cylinder is oscillating predominantly transverse to the fluid flow. Furthermore, the oscillation frequency is synchronized to the vortex shedding frequency of a stationary cylinder, except at very high reduced velocities. Visualizations of the pressure and vorticity in the wake reveal the mechanisms behind the motion of the cylinder.     

Vortex shedding from a square cylinder in presence of a moving wall

A numerical study on the flow past a square cylinder placed parallel to a wall, which is moving at the speed of the far field has been made. Flow has been investigated in the laminar Reynolds number (based on the cylinder length) range. We have studied the flow field for different values of the cylinder to wall separation length. The governing unsteady Navier–Stokes equations are discretized through the finite volume method on a staggered grid system. A SIMPLE type of algorithm has been used to compute the discretized equations iteratively. A shear layer of negative vortex generates along the surface of the wall, which influences the vortex shedding behind the cylinder. The flow‐field is distinct from the flow in presence of a stationary wall. An alternate vortex shedding occurs for all values of gap height in the unsteady regime of the flow. The strong positive vortex pushes the negative vortex upwards in the wake. The gap flow in the undersurface of the cylinder is strong and the velocity profile overshoots. The cylinder experiences a downward force for certain values of the Reynolds number and gap height. The drag and lift are higher at lower values of the Reynolds number. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical investigation of confined wakes behind a square cylinder in a channel

The structure of confined wakes behind a square cylinder in a channel is investigated via the numerical solution of the unsteady Navier–Stokes equations. Vortex shedding behind the cylinder induces periodicity in the flow field. Details of the phenomenon are simulated through numerical flow visualization. The unsteady periodic wake can be characterized by the Strouhal number, which varies with the Reynolds number and the blockage ratio of the channel. The periodicity of the flow is, however, damped in the downstream region of a long duct. This damping may be attributed to the influence of side walls on the flow structure.     

URANS vs. experiments of flow induced motions of multiple circular cylinders with passive turbulence control

Two-dimensional Unsteady Reynolds-Average Navier–Stokes equations with the Spalart–Allmaras turbulence model are used to simulate the flow induced motions of multiple circular cylinders with passive turbulence control (PTC) in steady uniform flow. Four configurations with 1, 2, 3, and 4 cylinders in tandem are simulated and studied at a series of Reynolds numbers in the range of 30 000<Re<120 000. Simulation results are verified by experimental data measured in the Marine Renewable Energy Laboratory. Good agreement was observed between the values of vorticity, amplitude ratio, and frequency ratio predicted by numerical simulations and experimental measurements. The amplitude and frequency response show the initial and upper branches in vortex induced vibration (VIV), transition from VIV to galloping, and galloping branch for all PTC-cylinders. The maximum amplitude of 2.9 diameters for the first cylinder is achieved at Re=104 356 in the numerical results. Compared with the first cylinder, the VIV initial branch starts at higher Re for the downstream cylinders due to the presence of the upstream cylinder(s). 2P and 2P+2S vortex patterns are observed at Re=62 049 and Re=90 254 for the single PTC-cylinder. Furthermore, the shed vortices of the downstream cylinders are strongly disrupted and modified by the vortices shed from the upstream one in the cases of multiple PTC-cylinders.     

Modes of vortex formation and transition to three-dimensionality in the wake of a freely vibrating cylinder

The present work is aimed to give some insight into the relation between vortex shedding modes and transition to three-dimensionality in the wake of a freely vibrating cylinder by establishing a numerical model and analyzing the relevant results of two- and three-dimensional simulations. The compressible flow past an elastically-mounted cylinder is solved by using the immersed boundary method (IB method). The cylinder is free to vibrate in the transverse direction with zero structure damping. The response of displacement amplitude is studied with the variation of reduced velocity. Whether P+S mode exists in three-dimensional flow and the occurrence of 2P mode is caused by flow transition from two-dimensional to three-dimensional are problems of concern. Both 2P and P+S wake modes are observed in two- and three-dimensional simulations. The numerical results indicate that the flow transition from two-dimensional to three-dimensional is coupled with the cylinder vibration in the synchronization/lock-in regime. The wake formation given by three-dimensional simulations suggests that the P+S mode might exist in reality when the flow is reverted to two-dimensional by vortex induced vibration (VIV) at Re=300–350. When Reynolds number increases to 425, the wake formation undergoes transition to three-dimensionality and 2P mode is observed. The effect of mass ratio on the flow transition to three-dimensionality is studied. The relationship between wake modes and aerodynamic forces is discussed.     

Reduction in drag and vortex shedding frequency through porous sheath around a circular cylinder

A numerical study on the laminar vortex shedding and wake flow due to a porous‐wrapped solid circular cylinder has been made in this paper. The cylinder is horizontally placed, and is subjected to a uniform cross flow. The aim is to control the vortex shedding and drag force through a thin porous wrapper around a solid cylinder. The flow field is investigated for a wide range of Reynolds number in the laminar regime. The flow in the porous zone is governed by the Darcy–Brinkman–Forchheimer extended model and the Navier–Stokes equations in the fluid region. A control volume approach is adopted for computation of the governing equations along with a second‐order upwind scheme, which is used to discretize the convective terms inside the fluid region. The inclusion of a thin porous wrapper produces a significant reduction in drag and damps the oscillation compared with a solid cylinder. Dependence of Strouhal number and drag coefficient on porous layer thickness at different Reynolds number is analyzed. The dependence of Strouhal number and drag on the permeability of the medium is also examined. Copyright © 2009 John Wiley & Sons, Ltd.     

Direct numerical simulation of three-dimensional flow past a yawed circular cylinder of infinite length

Direct numerical simulation of flow past a stationary circular cylinder at yaw angles (α) in the range of 0–60° was conducted at Reynolds number of 1000. The three-dimensional (3-D) Navier–Stokes equations were solved using the Petrov–Galerkin finite element method. The transition of the flow from 2-D to 3-D was studied. The phenomena that were observed in flow visualization, such as the streamwise vortices, the vortex dislocation and the instability of the shear layer, were reproduced numerically. The effects of the yaw angle on wake structures, vortex shedding frequency and hydrodynamic forces of the cylinder were investigated. It was found that the Strouhal number at different yaw angles (α) follows the independence principle. The mean drag coefficient agrees well with the independence principle. It slightly increases with the increase of α and reaches a maximum value at α=60°, which is about 10% larger than that when α=0°. The root-mean-square (r.m.s.) values of the lift coefficient are noticeably dependent on α.     

Application of wake oscillators to two-dimensional vortex-induced vibrations of circular cylinders in oscillatory flows

A nonlinear time-domain simulation model for predicting two-dimensional vortex-induced vibration (VIV) of a flexibly mounted circular cylinder in planar and oscillatory flow is presented. This model is based on the utilization of van der Pol wake oscillators, being unconventional since wake oscillators have typically been applied to steady flow VIV predictions. The time-varying relative flow–cylinder velocities and accelerations are accounted for in deriving the coupled hydrodynamic lift, drag and inertia forces leading to the cylinder cross-flow and in-line oscillations. The system fluid–structure interaction equations explicitly contain the time-dependent and hybrid trigonometric terms. Depending on the Keulegan–Carpenter number (KC) incorporating the flow maximum velocity and excitation frequency, the model calibration is performed, entailing a set of empirical coefficients and expressions as a function of KC and mass ratio. Parametric investigations in cases of varying KC, reduced flow velocity, cylinder-to-flow frequency ratio and mass ratio are carried out, capturing some qualitative features of oscillatory flow VIV and exploring the effects of system parameters on response prediction characteristics. The model dependence of hydrodynamic coefficients on the Reynolds number is studied. Discrepancies and limitations versus advantages of the present model with different feasible solution scenarios are illuminated to inform the implementation of wake oscillators as a computationally efficient prediction model for VIV in oscillatory flows.     

Numerical simulation of the unsteady flow over an elliptic cylinder at different orientations

A numerical method is developed for investigating the two‐dimensional unsteady viscous flow over an inclined elliptic cylinder placed in a uniform stream of infinite extent. The direction of the free stream is normal to the cylinder axis and the flow field unsteadiness arises from two effects, the first is due to the flow field development following the start of the motion and the second is due to vortex shedding in the wake region. The time‐dependent flow is governed by the full conservation equations of mass and momentum with no boundary layer approximations. The parameters involved are the cylinder axis ratio, Reynolds number and the angle of attack. The investigation covers a Reynolds number range up to 5000. The minor–major axis ratio of the elliptic cylinder ranges between 0.5 and 0.6, and the angle of attack ranges between 0° and 90°. A series truncation method based on Fourier series is used to reduce the governing Navier–Stokes equations to two coupled infinite sets of second‐order differential equations. These equations are approximated by retaining only a finite number of terms and are then solved by approximating the derivatives using central differences. The results reveal an unusual phenomenon of negative lift occurring shortly after the start of motion. Various comparisons are made with previous theoretical and experimental results, including flow visualizations, to validate the solution methodology. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical simulation of vortex-induced vibration of a circular cylinder at low mass-damping using RANS code

Fundamental research on vortex-induced vibration (VIV) of a circular cylinder is still needed to build more rational VIV analysis tools for slender marine structures. Numerical results are presented for the response of an elastically mounted rigid cylinder at low mass damping constrained to oscillate transversely to a free stream. A two-dimensional Reynolds-averaged Navier–Stokes (RANS) code equipped with the SST k−ω turbulence model is applied for the numerical calculations. The numerical results are compared in detail with recent experimental and computational work. The Reynolds-averaging procedure erases the random disturbances in the vortex shedding process, so that the comparison between experimental data and the numerical results obtained by RANS codes may reveal some random characteristics of the VIV response. How random disturbance affects the observation in the experiments is discussed in this paper and the issues influencing the appearance of the upper branch in experiments are especially investigated. The absence of the upper branch in RANS simulations is explained in depth on account of discrepancies, which exist between experiments and RANS simulations. In addition, the formation of the 2P vortex shedding mode and its transition through the lock-in region are well reproduced in this investigation.     

Wake modes of a cylinder undergoing free streamwise vortex-induced vibrations

Simultaneous measurements of the response of a circular cylinder experiencing vortex-induced vibrations (VIVs) in the streamwise direction and the resulting wake field were performed for a range of reduced velocities using time-resolved Particle-Image Velocimetry in the Reynolds number range 450–3700. The dominant vortex shedding mode was identified using phase-averaged vorticity fields. The cylinder response amplitude was characterised by two response branches, separated by a low amplitude region at resonance, as has been previously reported in the literature. During the first response branch the wake exhibited not only the symmetric S-I mode, but also the alternate A-II mode at slightly higher reduced velocities. For both modes, the vortices were observed to be shed at the cylinder response frequency, but rearranged downstream into a more stable structure in which the velocity fluctuations were no longer synchronised to the cylinder motion. A special case of the A-II mode, referred to as the SA mode, was found to dominate in the second response branch and the low amplitude region, while the far wake and the cylinder motion were synchronised (lock-in). A change in the timing of the vortex shedding with respect to the cylinder motion was observed between the low amplitude region and the second response branch. This is likely to correspond to a change in the fluid forcing and levels of excitation, and may explain the variation in the cylinder amplitude observed in this region. Lock-in and the second response branch were found to coincide with a contraction of the wake and an increase in strength of the shed vortices. This work reveals the inherent differences between the extensively studied case of transverse-only VIV and the streamwise-only case, which is crucial if the wealth of information available on transverse VIV is to be extended to the more practical two degree-of-freedom case.     

OBLIQUE VORTEX SHEDDING BEHIND TAPERED CYLINDERS

The vortex shedding in the wake behind linearly tapered circular cylinders has been considered for the two taper ratios 75:1 and 100:1. The Reynolds number based on the velocity of the incoming flow and the largest diameter was in the range from 130 to 180. The low Reynolds number assured that laminar flow prevailed in the entire flow field. The full unsteady three-dimensional Navier–Stokes equations were solved numerically with the view of exploring the rather complex vortex shedding phenomena caused by the variation of the natural shedding frequency along the span of the cylinder. The accurate computer simulations showed that this variation gave rise to discrete shedding cells, each with its own characteristic frequency and inclined with respect to the axis of the cylinder. Flow visualizations revealed that vortex dislocation and splitting took place in the numerically simulated flow fields. The computer simulations compared surprisingly well with the extensive laboratory experiments reported by Piccirillo & Van Atta in 1993 for a range of comparable conditions; this has enabled detailed analyses of other flow variables (notably pressure and vorticity) than those readily accessible in a physical experiment. However, distinct differences in the vortex dynamics are observed in some of the cases.     

Control of vortex shedding behind circular cylinder for flows at low Reynolds numbers

It has been observed by researchers in the past that vortex shedding behind circular cylinders can be altered, and in some cases suppressed, over a limited range of Reynolds numbers by proper placement of a second, much smaller, ‘control’ cylinder in the near wake of the main cylinder. Results are presented for numerical computations of some such situations. A stabilized finite element method is employed to solve the incompressible Navier–Stokes equations in the primitive variables formulation. At low Reynolds numbers, for certain relative positions of the main and control cylinder, the vortex shedding from the main cylinder is completely suppressed. Excellent agreement is observed between the present computations and experimental findings of other researchers. In an effort to explain the mechanism of control of vortex shedding, the streamwise variation of the pressure coefficient close to the shear layer of the main cylinder is compared for various cases, with and without the control cylinder. In the cases where the vortex shedding is suppressed, it is observed that the control cylinder provides a local favorable pressure gradient in the wake region, thereby stabilizing the shear layer locally. Copyright © 2001 John Wiley & Sons, Ltd.     

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.     

Vortex Induced Vibrations of a Pair of Cylinders at Reynolds Number 1000

Vortex induced vibrations of two equal-sized cylinders in tandem and staggered arrangement placed in uniform incompressible flow is studied. A stabilized finite element formulation is utilized to solve the governing equations. The Reynolds number for these 2D simulations is 1000. The cylinders are separated by 5.5 times the cylinder diameter in the streamwise direction. For the staggered arrangement, the cross-flow spacing between the two cylinders is 0.7 times the cylinder diameter. In this arrangement, the downstream cylinder lies in the wake of the upstream one and therefore experiences an unsteady inflow. The wake looses its temporal periodicity, beyond a few diameters downstream of the front cylinder. The upstream cylinder responds as an isolated single cylinder while the downstream one undergoes disorganized motion. Soft-lock-in is observed in almost all the cases.     

Suppression of vortex shedding for flow around a circular cylinder using optimal control

Adjoint formulation is employed for the optimal control of flow around a rotating cylinder, governed by the unsteady Navier–Stokes equations. The main objective consists of suppressing Karman vortex shedding in the wake of the cylinder by controlling the angular velocity of the rotating body, which can be constant in time or time‐dependent. Since the numerical control problem is ill‐posed, regularization is employed. An empirical logarithmic law relating the regularization coefficient to the Reynolds number was derived for 60?Re?140. Optimal values of the angular velocity of the cylinder are obtained for Reynolds numbers ranging from Re=60 to Re=1000. The results obtained by the computational optimal control method agree with previously obtained experimental and numerical observations. A significant reduction of the amplitude of the variation of the drag coefficient is obtained for the optimized values of the rotation rate. Copyright © 2002 John Wiley & Sons, Ltd.     

The onset of three‐dimensionality in an oscillating flow past a fixed circular cylinder

A computational study of the sinusoidally oscillating flow past a fixed circular cylinder has been performed to examine the viscous transition from 2D to 3D in the wake of the cylinder, i.e., the Honji instability. The primitive variables form of the Navier–Stokes equations was used with the discretization in the form of a combined finite‐difference/spectral‐method. Numerical results were obtained at a frequency parameter value of 196 and for a range of Keulegan–Carpenter (KC) numbers from 1 to 4. The calculations agreed quite well with experimental results. At KC=1, the wake is 2D; spanwise structures begin to appear at KC=2; separation has occurred by KC=3.2; and the vortex structures have become chaotic by KC=4. Calculated values of the force coefficients also agree extremely well with the experimental values. Copyright © 1999 John Wiley & Sons, Ltd.     

Global stability analysis of axisymmetric boundary layer over a circular cylinder

This paper presents a linear global stability analysis of the incompressible axisymmetric boundary layer on a circular cylinder. The base flow is parallel to the axis of the cylinder at inflow boundary. The pressure gradient is zero in the streamwise direction. The base flow velocity profile is fully non-parallel and non-similar in nature. The boundary layer grows continuously in the spatial directions. Linearized Navier–Stokes (LNS) equations are derived for the disturbance flow quantities in the cylindrical polar coordinates. The LNS equations along with homogeneous boundary conditions forms a generalized eigenvalues problem. Since the base flow is axisymmetric, the disturbances are periodic in azimuthal direction. Chebyshev spectral collocation method and Arnoldi’s iterative algorithm is used for the solution of the general eigenvalues problem. The global temporal modes are computed for the range of Reynolds numbers and different azimuthal wave numbers. The largest imaginary part of the computed eigenmodes is negative, and hence, the flow is temporally stable. The spatial structure of the eigenmodes shows that the disturbance amplitudes grow in size and magnitude while they are moving towards downstream. The global modes of axisymmetric boundary layer are more stable than that of 2D flat-plate boundary layer at low Reynolds number. However, at higher Reynolds number they approach 2D flat-plate boundary layer. Thus, the damping effect of transverse curvature is significant at low Reynolds number. The wave-like nature of the disturbance amplitudes is found in the streamwise direction for the least stable eigenmodes.     

On the study of vortex-induced vibration of a cylinder with helical strakes

While the effect of helical strakes on suppression of Vortex-Induced Vibrations (VIV) has been studied extensively, the mechanism of VIV mitigation using helical strakes is much less well documented in the literature. In the present study, a rigid circular cylinder of diameter d=80 mm attached with three-strand helical strakes of dimensions of 10d in pitch and 0.12d in height was tested in a wind tunnel. It was found that the helical strakes can reduce VIV by about 98%. Unlike the bare cylinder, which experiences lock-in over the reduced velocity in the range of 5-8.5, the straked cylinder does not show any lock-in region. In exploring the mechanism of VIV reduction by helical strakes, measurements in stationary bare and straked cylinder wakes using both a single X-probe at four different Reynolds numbers, i.e. Re=10 240, 20 430, 30 610 and 40 800, and two X-probes with variable separations in the spanwise direction at Re=20 430 were conducted. It was found that vortices shed from the straked cylinder are weakened significantly. The dominate frequency varies by about 30% over the range of x/d=10-40 in the streamwise direction while that differs by about 37.2% of the averaged peak frequency over a length of 3.125d in the spanwise direction. The latter is supported by the phase difference between the velocity signals measured at two locations separated in the spanwise direction. The correlation length of the vortex structures in the bare cylinder wake is much larger than that obtained in the straked cylinder wake. As a result, the straked cylinder wake agrees more closely with isotropy than the bare cylinder wake. Flow visualization on the plane perpendicular to the cylinder axis at Reynolds number of about 300 reveals small-scale vortices in the shear layers of the straked cylinder wake. However, these vortices do not roll up and interact with each other to form the well-organized Karman-type vortices. Flow visualization on the plane parallel to the cylinder axis shows vortex dislocation and swirling flow, which should be responsible for the variations of the peak frequency in the streamwise as well as spanwise directions.