A Parametric Study of Turbulent Flow Past a Circular Cylinder Using Large Eddy Simulation

2015 ◽  
Vol 137 (9) ◽  
Author(s):  
W. Sidebottom ◽  
A. Ooi ◽  
D. Jones

Flow over a circular cylinder at a Reynolds number of 3900 is investigated using large eddy simulations (LES) to assess the affect of four numerical parameters on the resulting flow-field. These parameters are subgrid scale (SGS) turbulence models, wall models, discretization of the advective terms in the governing equations, and grid resolution. A finite volume method is employed to solve the incompressible Navier–Stokes equations (NSE) on a structured grid. Results are compared to the experiments of Ong and Wallace (1996, “The Velocity Field of the Turbulent Very Near Wake of a Circular Cylinder,” Exp. Fluids, 20(6), pp. 441–453) and Lourenco and Shih (1993, “Characteristics of the Plane Turbulent Near Wake of a Circular Cylinder: A Particle Image Velocimetry Study,” private communication (taken from Ref. [2]); and the numerical results of Beaudan and Moin (1994, “Numerical Experiments on the Flow Past a Circular Cylinder at Sub-Critical Reynolds Number,” Technical Report No. TF-62), Kravchenko and Moin (2000, “Numerical Studies of Flow Over a Circular Cylinder at ReD = 3900,” Phys. Fluids, 12(2), pp. 403–417), and Breuer (1998, “Numerical and Modelling Influences on Large Eddy Simulations for the Flow Past a Circular Cylinder,” Int. J. Heat Fluid Flow, 19(5), pp. 512–521). It is concluded that the effect of the SGS models is not significant; results with and without a wall model are inconsistent; nondissipative discretization schemes, such as central finite difference methods, are preferred over dissipative methods, such as upwind finite difference methods; and it is necessary to properly resolve the boundary layer in the vicinity of the cylinder in order to accurately model the complex flow phenomena in the cylinder wake. These conclusions are based on the analysis of bulk flow parameters and the distribution of mean and fluctuating quantities throughout the domain. In general, results show good agreement with the experimental and numerical data used for comparison.

2011 ◽  
Vol 47 (1) ◽  
pp. 44-50 ◽  
Author(s):  
Stephen Wornom ◽  
Hilde Ouvrard ◽  
Maria Vittoria Salvetti ◽  
Bruno Koobus ◽  
Alain Dervieux

2013 ◽  
Vol 736 ◽  
pp. 414-443 ◽  
Author(s):  
Y. Ueda ◽  
T. Kida ◽  
M. Iguchi

AbstractThe long-time viscous flow about two identical rotating circular cylinders in a side-by-side arrangement is investigated using an adaptive numerical scheme based on the vortex method. The Stokes solution of the steady flow about the two-cylinder cluster produces a uniform stream in the far field, which is the so-called Jeffery’s paradox. The present work first addresses the validation of the vortex method for a low-Reynolds-number computation. The unsteady flow past an abruptly started purely rotating circular cylinder is therefore computed and compared with an exact solution to the Navier–Stokes equations. The steady state is then found to be obtained for $t\gg 1$ with ${\mathit{Re}}_{\omega } {r}^{2} \ll t$, where the characteristic length and velocity are respectively normalized with the radius ${a}_{1} $ of the circular cylinder and the circumferential velocity ${\Omega }_{1} {a}_{1} $. Then, the influence of the Reynolds number ${\mathit{Re}}_{\omega } = { a}_{1}^{2} {\Omega }_{1} / \nu $ about the two-cylinder cluster is investigated in the range $0. 125\leqslant {\mathit{Re}}_{\omega } \leqslant 40$. The convection influence forms a pair of circulations (called self-induced closed streamlines) ahead of the cylinders to alter the symmetry of the streamline whereas the low-Reynolds-number computation (${\mathit{Re}}_{\omega } = 0. 125$) reaches the steady regime in a proper inner domain. The self-induced closed streamline is formed at far field due to the boundary condition being zero at infinity. When the two-cylinder cluster is immersed in a uniform flow, which is equivalent to Jeffery’s solution, the streamline behaves like excellent Jeffery’s flow at ${\mathit{Re}}_{\omega } = 1. 25$ (although the drag force is almost zero). On the other hand, the influence of the gap spacing between the cylinders is also investigated and it is shown that there are two kinds of flow regimes including Jeffery’s flow. At a proper distance from the cylinders, the self-induced far-field velocity, which is almost equivalent to Jeffery’s solution, is successfully observed in a two-cylinder arrangement.


2012 ◽  
Vol 232 ◽  
pp. 471-476 ◽  
Author(s):  
Rui Zhao ◽  
Chao Yan

The flow past a circular cylinder at a subcritical Reynolds number 3900 was simulated by the method of detached-eddy simulation (DES). The objective of this present work is not to investigate the physical phenomena of the flow but to study modeling as well as numerical aspects which influence the quality of DES solutions in detail. Firstly, four typical spanwise lengths are chosen and the results are systematically compared. The trend of DES results along the span increment is different from previous large-eddy simulation (LES) investigation. A wider spanwise length does not necessary improve the results. Then, the influence of mesh resolution is studied and found that both too coarse and over refined grids will deteriorate the performance of DES. Finally, different orders of numerical schemes are applied in the inviscid fluxes and the viscous terms. The discrepancies among different schemes are found tiny. However, the instantaneous flow structures produced by 5th order WENO with 4th order central differencing scheme are more abundant than the others. That is, for the time-averaged quantities, the second-order accurate schemes are effective enough, whereas the higher-order accurate methods are needed to resolve the transient characteristics of the flow.


1991 ◽  
Vol 113 (3) ◽  
pp. 377-383 ◽  
Author(s):  
Xuegeng Wang ◽  
Charles Dalton

A finite-difference study of the sinusoidally oscillating flow past a fixed circular cylinder is made using vorticity and stream function as the dependent variables. Calculations are performed for conditions which lead to both a symmetric wake and an unsymmetric wake. The Reynolds number ranges from 100 to 3000 and the Keulegan-Carpenter number ranges from 1 to 12. A hybrid differencing scheme is introduced to provide a stable for large values of the parameters. Good comparison to flow visualization results and calculated force coefficients is found. The results are given a physical interpretation for the various vortex patterns observed.


Sign in / Sign up

Export Citation Format

Share Document