Simulation of Vortex-Shedding Flow About a Circular Cylinder at High Reynolds Numbers

1990 ◽  
Vol 112 (2) ◽  
pp. 155-161 ◽  
Author(s):  
Charles C. S. Song ◽  
Mingshun Yuan

Vortex shedding over a circular cylinder is modeled based on the weakly compressible flow equations with a simple subgrid scale turbulence model and a simple hybrid boundary condition. An explicit finite volume method is used. A subcritical and a supercritical case are computed. It is shown that the large-scale vortex-shedding phenomenon, the primary vortices, and the related oscillatory lift and drag can be calculated fairly well with a grid system coarser than the boundary layer thickness. The secondary vortices and the related higher frequency oscillations are also calculated by using somewhat finer grids.

1970 ◽  
Vol 4 (1) ◽  
pp. 27-42 ◽  
Author(s):  
Md Mahbubar Rahman ◽  
Md. Mashud Karim ◽  
Md Abdul Alim

The dynamic characteristics of the pressure and velocity fields of unsteady incompressible laminar and turbulent wakes behind a circular cylinder are investigated numerically and analyzed physically. The governing equations, written in the velocity pressure formulation are solved using 2-D finite volume method. The initial mechanism for vortex shedding is demonstrated and unsteady body forces are evaluated. The turbulent flow for Re = 1000 & 3900 are simulated using k-? standard, k-? Realizable and k-? SST turbulence models. The capabilities of these turbulence models to compute lift and drag coefficients are also verified. The frequencies of the drag and lift oscillations obtained theoretically agree well with the experimental results. The pressure and drag coefficients for different Reynolds numbers were also computed and compared with experimental and other numerical results. Due to faster convergence, 2-D finite volume method is found very much prospective for turbulent flow as well as laminar flow.Keywords: Viscous unsteady flow, laminar & turbulent flow, finite volume method, circular cylinder.DOI: 10.3329/jname.v4i1.914Journal of Naval Architecture and Marine Engineering 4(2007) 27-42


Author(s):  
Shufan Hu ◽  
Chen Zhang ◽  
Hong Liu ◽  
Fuxin Wang

This paper studied the unsteadiness of vortex motion produced by a three-dimensional wing section with horn/ridge ice contamination. Using improved delayed detached eddy simulation method, multi-scale vortex and their associated flow structures were successfully captured. Results have shown a diversity of unsteadiness scales at different time series, including shear layer instability, vortex pairing, co-rotating and breaking up. Proper orthogonal decomposition was then introduced to extract the characteristic vortex shedding modes with scheduling the eigenvalues λi from large to small. The dominate and secondary proper orthogonal decomposition modes under horn ice condition were displayed, which could be illustrated as fluctuations near recirculation zone, and large-scale vortex shedding/reattaching motion, respectively. The proper orthogonal decomposition modal characteristics for ridge ice showed that vortex scales varied from large to small. The trajectory of large-scale vortex reattaching and co-rotating exist simultaneously with the pressure peak and recover, which also verified the association of proper orthogonal decomposition modes with different scales of vortices. Future works would be presented on demonstration of the complex structures and the dynamic features in such flow.


Author(s):  
Antoine Placzek ◽  
Jean-Franc¸ois Sigrist ◽  
Aziz Hamdouni

The numerical simulation of the flow past a circular cylinder forced to oscillate transversely to the incident stream is presented here for a fixed Reynolds number equal to 100. The 2D Navier-Stokes equations are solved with a classical Finite Volume Method with an industrial CFD code which has been coupled with a user subroutine to obtain an explicit staggered procedure providing the cylinder displacement. A preliminary work is conducted in order to check the computation of the wake characteristics for Reynolds numbers smaller than 150. The Strouhal frequency fS, the lift and drag coefficients CL and CD are thus controlled among other parameters. The simulations are then performed with forced oscillations f0 for different frequency rations F = f0/fS in [0.50–1.50] and an amplitude A varying between 0.25 and 1.25. The wake characteristics are analysed using the time series of the fluctuating aerodynamic coefficients and their FFT. The frequency content is then linked to the shape of the phase portrait and to the vortex shedding mode. By choosing interesting couples (A,F), different vortex shedding modes have been observed, which are similar to those of the Williamson-Roshko map.


1966 ◽  
Vol 25 (3) ◽  
pp. 481-494 ◽  
Author(s):  
B. J. Vickery

The paper presents the results of measurements of fluctuating lift and drag on a long square cylinder. The measurements include the correlation of lift along the cylinder and the distribution of fluctuating pressure on a cross-section. The magnitude of the fluctuating lift was found to be considerably greater than that for a circular cross-section and the spanwise correlation much stronger.It was found that the presence of large-scale turbulence in the stream had a marked influence on both the steady and the fluctuating forces. The most significant changes were at small angles of attack (%alpha; < 10°) and included a reduction in base suction and a decrease in fluctuating lift of about 50%.


1976 ◽  
Vol 78 (3) ◽  
pp. 561-576 ◽  
Author(s):  
A. Richter ◽  
E. Naudascher

The fluctuating lift and drag acting on a long, rigidly supported circular cylinder placed symmetrically in a narrow rectangular duct were investigated for various blockage percentages over a wide range of Reynolds numbers around the critical value. The data obtained permit a full assessment of the effect of confinement on the mean-drag coefficient, the root-mean-square values of both the drag and the lift fluctuations, the Strouhal number of the dominant vortex shedding, and the Reynolds number marking transition from laminar to turbulent flow separation. Besides experimental information on a subject on which little is known so far, the paper provides a basis for the deduction of better correction procedures concerning the effects of blockage.


This paper reviews how Kolmogorov postulated for the first time the existence of a steady statistical state for small-scale turbulence, and its defining parameters of dissipation rate and kinematic viscosity. Thence he made quantitative predictions of the statistics by extending previous methods of dimensional scaling to multiscale random processes. We present theoretical arguments and experimental evidence to indicate when the small-scale motions might tend to a universal form (paradoxically not necessarily in uniform flows when the large scales are gaussian and isotropic), and discuss the implications for the kinematics and dynamics of the fact that there must be singularities in the velocity field associated with the - 5/3 inertial range spectrum. These may be particular forms of eddy or ‘eigenstructure’ such as spiral vortices, which may not be unique to turbulent flows. Also, they tend to lead to the notable spiral contours of scalars in turbulence, whose self-similar structure enables the ‘box-counting’ technique to be used to measure the ‘capacity’ D K of the contours themselves or of their intersections with lines, D' K . Although the capacity, a term invented by Kolmogorov (and studied thoroughly by Kolmogorov & Tikhomirov), is like the exponent 2 p of a spectrum in being a measure of the distribution of length scales ( D' K being related to 2 p in the limit of very high Reynolds numbers), the capacity is also different in that experimentally it can be evaluated at local regions within a flow and at lower values of the Reynolds number. Thus Kolmogorov & Tikhomirov provide the basis for a more widely applicable measure of the self-similar structure of turbulence. Finally, we also review how Kolmogorov’s concept of the universal spatial structure of the small scales, together with appropriate additional physical hypotheses, enables other aspects of turbulence to be understood at these scales; in particular the general forms of the temporal statistics such as the high-frequency (inertial range) spectra in eulerian and lagrangian frames of reference, and the perturbations to the small scales caused by non-isotropic, non-gaussian and inhomogeneous large-scale motions.


Author(s):  
Kun Luo ◽  
Jianren Fan ◽  
Kefa Cen

A direct numerical simulation technique combined with a two-way coupling method was developed to study a gas–solid turbulent jet with a moderately high Reynolds number. The flow was weakly compressible and spatially developing. A high-resolution solver was performed for the gas phase flow-field and the Lagrangian method was used to trace particles. The modulations on flow structures and other turbulent characteristics by particles at different Stokes numbers were investigated. It is found that the particles at Stokes numbers of 0.01 and 50 can advance the development of the large-scale vortex structures and make the turbulence intensity profiles wider and lower, but the particles at a Stokes number of 1 delay the evolution of the large-scale vortex structures and decrease the turbulence intensities. The jet velocity half-width and the decay of the streamwise mean velocity in the jet centreline are reduced by all particles, in which particles at a Stokes number of 0.01 result in a larger reduction of the velocity half-width and particles at a Stokes number of 1 lead to a larger reduction of the streamwise mean velocity decay. All particles decrease the vorticity thickness, but increase the fluid momentum thickness. In addition, the two-way coupled particle distribution is more uniform than that of the one-way coupled case.


2017 ◽  
Vol 21 (6 Part A) ◽  
pp. 2553-2563 ◽  
Author(s):  
Yidan Song ◽  
Rui Zhu ◽  
Terrence Simon ◽  
Gongnan Xie

The flow over four square cylinders in an in-line, square arrangement was numerically investigated by using the finite volume method with CFD techniques. The working fluid is an incompressible ideal gas. The length of the sides of the array, L, is equal. The analysis is carried out for a Reynolds number of 300, with center-to-center distance ratios, L/D, ranging from 1.5 to 8.0. To fully understand the flow mechanism, details in terms of lift and drag coefficients and Strouhal numbers of the unsteady wake frequencies are analyzed, and the vortex shedding patterns around the four square cylinders are described. It is concluded that L/D has important effects on the drag and lift coefficients, vortex shedding frequencies, and flow field characteristics.


Sign in / Sign up

Export Citation Format

Share Document