NUMERICAL SIMULATION OF UNSTEADY FLOW AROUND A CIRCULAR CYLINDER AT HIGH REYNOLDS NUMBERS

A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusion, the vortex-in-cell method is used to obtain the convection velocity, and nascent vortices are created on a cylinder to satisfy the zero-slip condition. The impulsively started flow around a circular cylinder and the separation induced by a pair of incident vortices symmetrically approaching a circular cylinder have been successfully simulated by the hybrid scheme. The impulsively started flow from rest has been computed at Reynolds numbers 3000 and 9500. Comparisons are made with those results of finite-difference method, vortex method and flow visualization. Agreement is good. The particular attention has been paid to the evolutions of flow pattern. A topological analysis has been proposed in the region of the near wake. The bulge, isolated secondary vortex, a pair of secondary vortices, ' forewake phenomenon and other patterns are simulated numerically. The separation induced by a pair of incident vortices approaching a circular cylinder has been investigated by using the same scheme. The rebounding phenomenon of the incident vortex is observed and is attributed to the effect of the secondary vortex. In particular, we have found that a tertiary vortex can be formed near the surface; this phenomenon has been verified by flow visualization reported recently.

Flow Dynamics around Two Side by Side Circular Cylinders with Alternating Movements

The flow dynamics is analyzed through two-dimensional numerical simulations around two circular cylinders arranged side by side, with 4 combinations of alternating motions. All simulations are performed for Re = 1000, amplitude of oscillation (A) equal to 3, frequency ratio (fr) of 0.5, specific rotation (α) equal to 0.5 and different values of spacing ratio (L/D). It is verified that the combination of the type of movement, together with the position of one cylinder in relation to the other, exerts significant influence on the flow dynamics, as well as on the pressure distribution around the cylinder surface and on the average values of the fluid dynamics coefficients. The smallest value of the average pressure coefficient (Cp = -3.3), is obtained for the oscillating cylinder when placed side by side with the clockwise rotation cylinder, case 3 and L/D = 1.5. On the other hand, the lowest mean drag coefficient (Cd = 1.0788), is obtained for the cylinder with counterclockwise rotation, located in the lower position in relation to oscillating cylinder in the upper position, with spacing between them of 1.5. Furthermore, it is observed that the rotation movement is more effective in reducing drag than the rotation-oscillation movement, for the studied frequency ratio.

Unsteady Numerical Simulation of Flow around 2-D Circular Cylinder for High Reynolds Numbers

In this paper, 2-D computational analyses were conducted for unsteady high Reynolds number flows around a smooth circular cylinder in the supercritical and upper-transition flow regimes, i.e. 8.21×104〈Re〈1.54×106. The calculations were performed by means of solving the 2-D Unsteady Reynolds-Averaged Navier-Stokes （URANS） equations with a k-ε turbulence model. The calculated results, produced flow structure drag and lift coefficients, as well as Strouhal numbers. The findings were in good agreement with previous published data, which also supplied us with a good understanding of the flow across cylinders of different high Reynolds numbers. Meanwhile, an effective measure was presented to control the lift force on a cylinder, which points the way to decrease the vortex induced vibration of marine structure in future.

Investigation of scale efect for the computation of turbulent flow around a circular cylinder

In order to investigate the scale effect of turbulent flow around a circular cylinder, two similarity numbers （criteria） based on turbulent kinetic and dissipation rates associ- ated with the fluctuation characteristics of turbulence wake are deduced by analyzing the Reynolds averaged NavierStokes equations （RANS）. The RNG k-s models and finite volume method are used to solve the governing equations and the second-order implicit time and upwind space discretization algorithms are used to discrete the governing equations. A numerical computation of flow parameters around a two-dimensional circular cylinder with Reynolds numbers ranging from 102 to l07 is accomplished and the result indicates that the fluctuation of turbulence flow along the center line in the wake of circular cylinder can never be changed with increasing Reynolds numbers when Re ≥ 3 × 10^6. This conclusion is useful for controlling the scale of numerical calculations and for applying model test data to engineering practice.

NUMERICAL STUDY ON THE FLOW AROUND A CIRCULAR CYLINDER WITH SURFACE SUCTION OR BLOWING USING VORTICITY-VELOCITY METHOD

A vorticity-velocity method was used to study the incompressible viscous fluid flow around a circular cylinder with surface suction or blowing. The resulted high order implicit difference equations were effeciently solved by the modified incomplete LU decomposition conjugate gradient scheme ( MILU-CG). The effects of surface suction or blowing' s position and strength on the vortex structures in the cylinder wake, as well as on the drag and lift forces at Reynoldes number Re = 100 were investigated numerically. The results show that the suction on the shoulder of the cylinder or the blowing on the rear of the cylinder can effeciently suppress the asymmetry of the vortex wake in the transverse direction and greatly reduce the lift force; the suction on the shoulder of the cylinder, when its strength is properly chosen, can reduce the drag force significantly, too.

Flow characteristics at low Reynolds number around two in-line circular cylinders with slits

This study investigated numerically the characteristics of laminar flow around two identical circular cylinders placed in tandem, with slits of the same width through their respective axis. The center to center distance between the cylinders and the slit orientation were varied to study their effects on the flow structure, lift and drag, and vortex shedding characteristics. It was found that three flow regimes could be distinguished, the transitions between which could be indicated by the sudden changes in drag and lift. Asymmetrically, configured slits destabilized the stagnant region between cylinders;whereas in-line slits connect the two cylinders to act as a single elongated bluff body, even at large cylinder separation, by stabilizing the stagnant region in between. These in turn strongly modified the transition between flow regimes. Vortex shedding was also strongly influenced by both slit configuration and cylinder separation.

Three-Dimensional Instability of Flow Around Two Circular Cylinders in Tandem Arrangement

The spatial evolution of vortices and transition to three-dimensionality in the wake of two circular cylinders in tandem arrangement have been numerically studied. An improved virtual body method developed from the virtual boundary method is used here. A Reynolds number range between 220 and 270 has been considered, and the spacing between two cylinders is selected as L/D=3 and L/D=3.5. When L/D=3, the secondary vortices of Mode-A are seen to appear at Re=240 and persist over the range of the Reynolds number of 240～270. When L/D=3.5, the similar critical Reynolds number has been found at Re=250. No obvious discontinuity has been found in the Strouhal-Reynolds number relationship, and this is different from three-dimensional flow around a single cylinder at the critical Reynolds number. The spanwise wavelength is about four times the diameter of the cylinder, and it is the characteristic wavelength for Mode-A instability. This paper can give some foremost insight into the three-dimensional instability of flow by complicated geometrical configuration.

Three-dimensional numerical simulation of the flow past a circular cylinder based on LES method

The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three dimensional N-S equations were obtained by the finite volume method. The focus of this numerical simulation was to research the characteristics of pressure distribution （drag and litt forces） and vortex tubes at high Reynolds numbers. The results of the calculations showed that the forces at every section in the spanwise direction of the cylinder were symmetrical about the middle section and smaller than the forces calculated in two dimensional cases. Moreover, the flow around the cylinder obviously presents three dimensional characteristics.

Experimental and Numerical Study of Three-Dimensional Unsteady Incompressible Flow Past Two Identical Circular Cylinders in Tandem Arrangements

The possible states in the flow past two identical cylinders in tandem arrangements are investigated. The effect of the gap （L/D = 1.5, 1.75 and 2） between the two cylinders at Reynolds number （Re = 52,639） is taken into consideration. The presence of three-dimensional flow structures was observed to include notable changes to the response of the flow as result of variation of cylinder separation. A number of planes （z/h = 0.02, 0.25, 0.5 and 0.98） were taken at 20 step times of interval 0.005 s. to cover the details of flow along the cylinders. CFD FLUENT program was used to detect the flow structure. It is observed that the gap between the two cylinders affects the flow regime, i.e., there is no distinct vortex shedding downstream of the first cylinder.

Large-eddy simulation of circular cylinder flow at subcritical Reynolds number:Turbulent wake and sound radiation

The flows past a number 3900 are simulated circular cylinder at Reynolds using large-eddy simulation （LES） and the far-field sound is calculated from the LES results. A low dissipation energy-conserving finite volume scheme is used to discretize the incompressible Navier- Stokes equations. The dynamic global coefficient version of the Vreman＇s subgrid scale （SGS） model is used to com- pute the sub-grid stresses. Curie＇s integral of Lighthill＇s acoustic analogy is used to extract the sound radiated from the cylinder. The profiles of mean velocity and turbulent fluctua- tions obtained are consistent with the previous experimental and computational results. The sound radiation at far field exhibits the characteristic of a dipole and directivity. The sound spectra display the -5/3 power law. It is shown that Vreman＇s SGS model in company with dynamic procedure is suitable for LES of turbulence generated noise.

Near Wake of a Horizontal Circular Cylinder in Stably Stratified Flows

The near wake of a circular cylinder in linearly stratified flows of finite depth was experimentally investigated by means of flow visualization and measurements of vortex shedding frequencies, at Reynolds numbers 3.5 × 103-1.2 × 104 and stratification parameters kd 0-2.0. The non-dimensional parameter kd is defined as kd = Nd/U, where N is the Brunt-Vaisala frequency, d, the diameter of the cylinder, and U, the approaching flow velocity. The study demonstrates that as kd increases from zero, the vortex shedding from a circular cylinder progressively strengthens, while the Strouhal number gradually becomes lower than that for homogeneous flow. This phenomenon can be explained by the effect of the increasingly stable stratification which enhances the two-dimensionality of the near-wake flow of the circular cylinder;the enhanced two-dimensionality of the flow strengthens the roll-up of the separated shear layer. Above a certain value of kd, however, vortex formation and shedding are strongly suppressed and the Strouhal number rises sharply. This observation is attributable to the development of stationary lee waves downstream of the circular cylinder because the lee waves strongly suppress vertical fluid motions.

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

Computation of Viscous Uniform and Shear Flow over A Circular Cylinder by A Finite Element Method

The incompressible viscous uniform and shear flow past a circular cylinder is studied. The two-dimensional Navier-Stokes equations are solved by a finite element method. The governing equations are discretized by a weighted residual method in space. The stable three-step scheme is applied to the momentum equations in the time integration. The numerical model is firstly applied to the computation of the lid-driven cavity flow for its validation. The computed results agree well with the measured data and other numerical results. Then, it is used to simulate the viscous uniform and shear flow over a circular cylinder for Reynolds numbers from 100 to 1000. The transient time interval before the vortex shedding occurs is shortened considerably by introduction of artificial perturbation. The computed Strouhal number, drag and lift coefficients agree well with the experimental data. The computation shows that the finite element model can be successfully applied to the viscous flow problem.

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.

Flow-Induced Stream-Wise Vibration of Circular Cylinders

Results from a series of studies on the stream-wise vibration of a circular cylinder verifying Japan Society of Mechanical Engineers Standard S012-1998, Guideline for Evaluation of Flow-induced Vibration of a Cylindrical Structure in a Pipe, are summarized and discussed in this paper. Experiments were carried out in a water tunnel and in a wind tunnel using a two-dimensional cylinder model elastically supported at both ends of the cylinder and a cantilevered cylinder model with a finite span length that was elastically supported at one end. These cylinder models were allowed to vibrate with one degree of freedom in the stream-wise direction. In addition, we adopted a cantilevered cylinder model that vibrated with two degrees of freedom in both the stream-wise and cross-flow directions under the same vibration conditions as an actual thermocouple well. The value of the Scruton number (structural damping parameter) was changed over a wide range, so as to evaluate the value of the critical Scruton number that suppressed vibration of the cylinder. For the two-dimensional cylinder, two different types of stream-wise excitations appeared in the reduced velocity range of approximately half of the resonance-reduced velocity. For the stream-wise vibration in the first excitation region, due to a symmetric vortex flow, the response amplitudes were sensitive to the Scruton number, while the shedding frequency of alternating vortex flow was locked-in to half of the Strouhal number of vibrating frequency of a cylinder in the second excitation region. In addition, the effects of the aspect ratio of a cantilevered cylinder on the flow-induced vibration characteristics were clarified and compared with the results of a two-dimensional cylinder. When a cantilevered circular cylinder with a finite length vibrates with one degree of freedom in the stream-wise di-rection, it is found that acylinder with a small aspect ratio has a single excitation region, whereas a cylinder with a large aspect ratio has two excitation regions. Furthermore, the vibration mechanism of a symmetric vortex flow was investigated by installing a splitter plate in the wake to prevent shedding of alternating vortices. The vibration amplitude of acylinder with a splitter plate increased surprisingly more than the amplitude of a cylinder without a splitter plate. For a cantilevered cylinder vibrating with two degrees of freedom, the Lissajous figure of vibration of the first excitation region shows the trajectories of elongated elliptical shapes, and in the second excitation region, the Lissajous trajectories draw a figure “8”. The results and information from these experimental studies proved that Standard S012-1998 provides sufficient design methods for suppressing hazardous vibrations of cylinders in liquid flows.

Combined Action of Uniform Flow and Oscillating Flow Around Marine Riser at Low Keulegan-Carpenter Number

With the increase of petroleum and gas production in deep ocean, marine risers of circular cylinder shape are widely used in the offshore oil and gas platform. In order to research the hydrodynamic performance of marine risers, the dynamic mesh technique and User-Defined Function(UDF) are used to simulate the circular cylinder motion. The motion of a transversely oscillating circular cylinder in combination of uniform flow and oscillating flow is simulated. The uniform flow and oscillating flow both are in x direction. SIMPLE algorithm is used to solve the Navier-Stokes equations. The User-Defined Function is used to control the cylinder transverse vibration and the inlet flow. The lift and drag coefficient changing with time and the map of vorticity isolines at different phase angle are obtained. Force time histories are shown for uniform flow at Reynolds number(Re) of 200 and for the combination of uniform and oscillating flows. With the increase of amplitude of oscillating flow in combined flow, the change of lift amplitude is not sensitive to the the change of cylinder oscillating frequency. Lift amplitude increases with the increase of oscillating flow amplitude in the combined flow, but there is no definite periodicity of the lift coefficient. The drag and inertia force coefficients change when the maximum velocity of the oscillating flow increases in the combined flow. The vortex shedding near the circular cylinder shows different characteristics.

Convective heat transfer from two rotating circular cylinders in tandem arrangement by using lattice Boltzmann method

The numerical investigation of the two-dimensional laminar flow past two ro- tating circular cylinders in the tandem arrangement is conducted by the lattice Boltzmann method. The numerical strategy is used for dealing with curved and moving boundaries of the second-order accuracy for velocity and temperature fields. The effects of various rotational speed ratios and gap spacing are studied with the Reynolds number of 100 and the Prandtl number of 0.71. A varied range of rotational speed ratios are investigated for four different gap spacing, i.e., 3.0, 1.5, 0.7, and 0.2. The results show that, for the first cylinder, the lift and drag coefficients for large gap spacing are similar to those for a single cylinder; for the second cylinder, the lift coefficient descends with the increase in the angular velocity for all gap spacing, while the drag coefficient ascends except for the gap spacing of 3.0. The results of the averaged periodic Nusselt number on the surface of the cylinders show that, for small distances between the cylinders and low angular velocities, conduction is a dominant mechanism of heat transfer, but for large distances and high angular velocities, convection is the main mechanism of heat transfer.

Applications of URANS on predicting unsteady turbulent separated flows

Accurate prediction of unsteady separated turbulent flows remains one of the toughest tasks and a practi cal challenge for turbulence modeling. In this paper, a 2D flow past a circular cylinder at Reynolds number 3,900 is numerically investigated by using the technique of unsteady RANS （URANS）. Some typical linear and nonlinear eddy viscosity turbulence models （LEVM and NLEVM） and a quadratic explicit algebraic stress model （EASM） are evaluated. Numerical results have shown that a high-performance cubic NLEVM, such as CLS, are superior to the others in simulating turbulent separated flows with unsteady vortex shedding.

Performance Investigation of 2D Lattice Boltzmann Simulation of Forces on a Circular Cylinder

The lattice Boltzmann method (LBM) is employed to simulate the uniform flow past a circular cylinder. The performance of the two-dimensional LBM model on the prediction of force coefficients and vortex shedding frequency is investigated. The local grid refinement technique and second-order boundary condition for curved walls are applied in the calculations. It is found that the calculated vortex shedding frequency, drag coefficient and lift coefficient are consistent with experimental results at Reynolds nu...

ON SEPARATED SHEAR LAYER OF BLUNT CIRCULAR CYLINDER

Separated shear layer of blunt circular cylinder has been experimentally investigated for the Reynolds numbers (based on the diameter) ranging from 2.8 x 10(3) to 1.0 x 10(5), with emphasis on evolution of separated shear layer, its structure and distribution of Reynolds shear stress and turbulence kinetic energy. The results demonstrate that laminar separated shear layer experiences 2 similar to 3 times vortex merging before it reattaches, and turbulence separated shear layer takes 5 similar to 6 times vortex merging. In addition, relationship between dimensionless initial frequencies of K-H instability and Reynolds numbers is identified, and reasons for the decay of turbulence kinetic energy and Reynolds shear stress in reattachment region are discussed.