Numerical solutions of flow past a circular cylinder at Reynolds numbers up to 160

1998 ◽  
Vol 12 (6) ◽  
pp. 1200-1205 ◽  
Author(s):  
Jeongyoung Park ◽  
Kiyoung Kwon ◽  
Haecheon Choi
1980 ◽  
Vol 98 (4) ◽  
pp. 819-855 ◽  
Author(s):  
Bengt Fornberg

Numerical solutions have been obtained for steady viscous flow past a circular cylinder at Reynolds numbers up to 300. A new technique is proposed for the boundary condition at large distances and an iteration scheme has been developed, based on Newton's method, which circumvents the numerical difficulties previously encountered around and beyond a Reynolds number of 100. Some new trends are observed in the solution shortly before a Reynolds number of 300. As vorticity starts to recirculate back from the end of the wake region, this region becomes wider and shorter. Other flow quantities like position of separation point, drag, pressure and vorticity distributions on the body surface appear to be quite unaffected by this reversal of trends.


1997 ◽  
Vol 345 ◽  
pp. 101-131
Author(s):  
M. D. KUNKA ◽  
M. R. FOSTER

Because of the importance of oscillatory components in the oncoming flow at certain oceanic topographic features, we investigate the oscillatory flow past a circular cylinder in an homogeneous rotating fluid. When the oncoming flow is non-reversing, and for relatively low-frequency oscillations, the modifications to the equivalent steady flow arise principally in the ‘quarter layer’ on the surface of the cylinder. An incipient-separation criterion is found as a limitation on the magnitude of the Rossby number, as in the steady-flow case. We present exact solutions for a number of asymptotic cases, at both large frequency and small nonlinearity. We also report numerical solutions of the nonlinear quarter-layer equation for a range of parameters, obtained by a temporal integration. Near the rear stagnation point of the cylinder, we find a generalized velocity ‘plateau’ similar to that of the steady-flow problem, in which all harmonics of the free-stream oscillation may be present. Further, we determine that, for certain initial conditions, the boundary-layer flow develops a finite-time singularity in the neighbourhood of the rear stagnation point.


2013 ◽  
Vol 730 ◽  
pp. 145-161 ◽  
Author(s):  
Qiqi Wang ◽  
Jun-Hui Gao

AbstractThis paper analyses the adjoint solution of the Navier–Stokes equation. We focus on flow across a circular cylinder at three Reynolds numbers, ${\mathit{Re}}_{D} = 20, 100$ and $500$. The quantity of interest in the adjoint formulation is the drag on the cylinder. We use classical fluid mechanics approaches to analyse the adjoint solution, which is a vector field similar to a flow field. Production and dissipation of kinetic energy of the adjoint field is discussed. We also derive the evolution of circulation of the adjoint field along a closed material contour. These analytical results are used to explain three numerical solutions of the adjoint equations presented in this paper. The adjoint solution at ${\mathit{Re}}_{D} = 20$, a viscous steady state flow, exhibits a downstream suction and an upstream jet, the opposite of the expected behaviour of a flow field. The adjoint solution at ${\mathit{Re}}_{D} = 100$, a periodic two-dimensional unsteady flow, exhibits periodic, bean-shaped circulation in the near-wake region. The adjoint solution at ${\mathit{Re}}_{D} = 500$, a turbulent three-dimensional unsteady flow, has complex dynamics created by the shear layer in the near wake. The magnitude of the adjoint solution increases exponentially at the rate of the first Lyapunov exponent. These numerical results correlate well with the theoretical analysis presented in this paper.


1990 ◽  
Vol 220 ◽  
pp. 459-484 ◽  
Author(s):  
H. M. Badr ◽  
M. Coutanceau ◽  
S. C. R. Dennis ◽  
C. Ménard

The unsteady flow past a circular cylinder which starts translating and rotating impulsively from rest in a viscous fluid is investigated both theoretically and experimentally in the Reynolds number range 103 [les ] R [les ] 104 and for rotational to translational surface speed ratios between 0.5 and 3. The theoretical study is based on numerical solutions of the two-dimensional unsteady Navier–Stokes equations while the experimental investigation is based on visualization of the flow using very fine suspended particles. The object of the study is to examine the effect of increase of rotation on the flow structure. There is excellent agreement between the numerical and experimental results for all speed ratios considered, except in the case of the highest rotation rate. Here three-dimensional effects become more pronounced in the experiments and the laminar flow breaks down, while the calculated flow starts to approach a steady state. For lower rotation rates a periodic structure of vortex evolution and shedding develops in the calculations which is repeated exactly as time advances. Another feature of the calculations is the discrepancy in the lift and drag forces at high Reynolds numbers resulting from solving the boundary-layer limit of the equations of motion rather than the full Navier–Stokes equations. Typical results are given for selected values of the Reynolds number and rotation rate.


1957 ◽  
Vol 2 (3) ◽  
pp. 237-262 ◽  
Author(s):  
Ian Proudman ◽  
J. R. A. Pearson

This paper is concerned with the problem of obtaining higher approximations to the flow past a sphere and a circular cylinder than those represented by the well-known solutions of Stokes and Oseen. Since the perturbation theory arising from the consideration of small non-zero Reynolds numbers is a singular one, the problem is largely that of devising suitable techniques for taking this singularity into account when expanding the solution for small Reynolds numbers.The technique adopted is as follows. Separate, locally valid (in general), expansions of the stream function are developed for the regions close to, and far from, the obstacle. Reasons are presented for believing that these ‘Stokes’ and ‘Oseen’ expansions are, respectively, of the forms $\Sigma \;f_n(R) \psi_n(r, \theta)$ and $\Sigma \; F_n(R) \Psi_n(R_r, \theta)$ where (r, θ) are spherical or cylindrical polar coordinates made dimensionless with the radius of the obstacle, R is the Reynolds number, and $f_{(n+1)}|f_n$ and $F_{n+1}|F_n$ vanish with R. Substitution of these expansions in the Navier-Stokes equation then yields a set of differential equations for the coefficients ψn and Ψn, but only one set of physical boundary conditions is applicable to each expansion (the no-slip conditions for the Stokes expansion, and the uniform-stream condition for the Oseen expansion) so that unique solutions cannot be derived immediately. However, the fact that the two expansions are (in principle) both derived from the same exact solution leads to a ‘matching’ procedure which yields further boundary conditions for each expansion. It is thus possible to determine alternately successive terms in each expansion.The leading terms of the expansions are shown to be closely related to the original solutions of Stokes and Oseen, and detailed results for some further terms are obtained.


2007 ◽  
Vol 592 ◽  
pp. 89-115 ◽  
Author(s):  
A. MILIOU ◽  
A. DE VECCHI ◽  
S. J. SHERWIN ◽  
J. M. R. GRAHAM

Three-dimensional spectral/hp computations have been performed to study the fundamental mechanisms of vortex shedding in the wake of curved circular cylinders at Reynolds numbers of 100 and 500. The basic shape of the body is a circular cylinder whose centreline sweeps through a quarter section of a ring and the inflow direction lies on the plane of curvature of the quarter ring: the free stream is then parallel to the geometry considered and the part of the ring that is exposed to it will be referred to as the ‘leading edge’. Different configurations were investigated with respect to the leading-edge orientation. In the case of a convex-shaped geometry, the stagnation face is the outer surface of the ring: this case exhibited fully three-dimensional wake dynamics, with the vortex shedding in the upper part of the body driving the lower end at one dominant shedding frequency for the whole cylinder span. The vortex-shedding mechanism was therefore not governed by the variation of local normal Reynolds numbers dictated by the curved shape of the leading edge. A second set of simulations were conducted with the free stream directed towards the inside of the ring, in the so-called concave-shaped geometry. No vortex shedding was detected in this configuration: it is suggested that the strong axial flow due to the body's curvature and the subsequent production of streamwise vorticity plays a key role in suppressing the wake dynamics expected in the case of flow past a straight cylinder. The stabilizing mechanism stemming from the concave curved geometry was still found to govern the wake behaviour even when a vertical extension was added to the top of the concave ring, thereby displacing the numerical symmetry boundary condition at this point away from the top of the deformed cylinder. In this case, however, the axial flow from the deformed cylinder was drawn into the wake of vertical extension, weakening the shedding process expected from a straight cylinder at these Reynolds numbers. These considerations highlight the importance of investigating flow past curved cylinders using a full three-dimensional approach, which can properly take into account the role of axial velocity components without the limiting assumptions of a sectional analysis, as is commonly used in industrial practice. Finally, towing-tank flow visualizations were also conducted and found to be in qualitative agreement with the computational findings.


Sign in / Sign up

Export Citation Format

Share Document