Suspension flow past a cylinder: particle interactions with recirculating wakes

2014 ◽  
Vol 760 ◽  
Author(s):  
Hamed Haddadi ◽  
Shahab Shojaei-Zadeh ◽  
Kevin Connington ◽  
Jeffrey F. Morris

AbstractExperimental observations of the flow of a suspension of solid fraction ${\it\phi}\approx 0.084$ over a circular cylindrical post in a shallow microchannel (depth smaller than the cylinder radius) find that the recirculating wake behind the obstacle at moderate Reynolds numbers is depleted or devoid of particles. Particles injected into the wake exit to regain the depleted state. By numerical simulation of the discrete particle motion, the basis for the depletion behind the cylinder is studied; rather than a shallow channel, the numerical simulations consider a periodic domain, mimicking the flow past an infinite cylinder. The Reynolds number is defined, using the average axial velocity ${\bar{U}}$, diameter of the obstacle $D$ and the kinematic viscosity of the suspension ${\it\nu}$, as $Re={\bar{U}}D/{\it\nu}$, and is studied for $Re<30$ in the simulation – conditions for which the pure fluid exhibits an extended steady closed-streamline (recirculating) wake behind the cylinder; unsteadiness is found to be suppressed by the channel walls in the experiments, allowing steady flow at a larger $Re$ than expected for an infinite cylinder (up to at least $Re=300$). The simulations use the lattice-Boltzmann method to determine the motion of the fluid and neutrally buoyant particles. The trajectory of a single particle (small relative to the cylinder) shows migration to a limit cycle inside the wake. With an increase of the number of particles in the wake alone (no particles in the free stream), particles can escape the wake due to velocity fluctuations. Simulation of the flow of suspensions of ${\it\phi}=0.04,0.06$ and 0.08 demonstrates that there is particle exchange between the wake and the free stream; the net flux of particles out of the wake leads to a particle-depleted wake, qualitatively very similar to the experimental observation.

2009 ◽  
Vol 61 (6) ◽  
pp. 658-682 ◽  
Author(s):  
A. P. Singh ◽  
A. K. De ◽  
V. K. Carpenter ◽  
V. Eswaran ◽  
K. Muralidhar

1994 ◽  
Vol 258 ◽  
pp. 287-316 ◽  
Author(s):  
C. Norberg

The investigation is concentrated on two important quantities – the Strouhal number and the mean base suction coefficient, both measured at the mid-span position. Reynolds numbers from about 50 to 4 × 104 were investigated. Different aspect ratios, at low blockage ratios, were achieved by varying the distance between circular end plates (end plate diameter ratios between 10 and 30). It was not possible, by using these end plates in uniform flow and at very large aspect ratios, to produce parallel shedding all over the laminar shedding regime. However, parallel shedding at around mid-span was observed throughout this regime in cases when there was a slight but symmetrical increase in the free-stream velocity towards both ends of the cylinder. At higher Re, the results at different aspect ratios were compared with those of a ‘quasi-infinite cylinder’ and the required aspect ratio to reach conditions independent of this parameter, within the experimental uncertainties, are given. For instance, aspect ratios as large as L/D = 60–70 were needed in the range Re ≈ 4 × 103–104. With the smallest relative end plate diameter and for aspect ratios smaller than 7, a bi-stable flow switching between regular vortex shedding and ‘irregular flow’ was found at intermediate Reynolds number ranges in the subcritical regime (Re ≈ 2 × 103).


2017 ◽  
Vol 14 (01) ◽  
pp. 1750002 ◽  
Author(s):  
Yi-Kun Wei ◽  
Xu-Qu Hu

Two-dimensional simulations of channel flow past an array of cylinders are carried out at high Reynolds numbers. Considering the thickness fluctuating effect on the equation of motion, a modified lattice Boltzmann method (LBM) is proposed. Special attention is paid to investigate the thickness fluctuations and vortex shedding mechanisms between 11 cylinders. Results for the velocity and vorticity differences are provided, as well as for the energy density and enstrophy spectra. The numerical results coincide very well with some published experimental data that was obtained by turbulent soap films. The spectra extracted from the velocity and vorticity fields are displayed from simulations, along with the thickness fluctuation spectrum H(k). Our results show that the statistics of thickness fluctuations resemble closely those of a passive scalar in turbulent flows.


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.


2021 ◽  
Vol 42 (7) ◽  
pp. 1715-1727
Author(s):  
Alexander Mazo ◽  
Valeriy Molochnikov ◽  
Evgeniy Kalinin ◽  
Anton Paereliy ◽  
Nickolay Dushin ◽  
...  

2010 ◽  
Vol 668 ◽  
pp. 304-334 ◽  
Author(s):  
V. B. L. BOPPANA ◽  
J. S. B. GAJJAR

The effect of blockage on the onset of instability in the two-dimensional uniform flow past a cascade of cylinders is investigated. The same techniques as those described in Gajjar & Azzam (J. Fluid Mech., vol. 520, 2004, p. 51) are used to tackle the generalized eigenvalue problem arising from a global stability analysis of the linearized disturbance equations. Results have been obtained for the various mode classes, and our results show that for the odd–even modes, which correspond to anti-phase oscillatory motion about the midplane between the cylinders and are the modes most extensively studied in the literature, the effect of blockage has a marginal influence on the critical Reynolds numbers for instability. This is in sharp contrast to results cited in many studies with a fully developed inlet flow past a cylinder placed between confining walls. We are also able to find other unstable modes and in particular for low blockage ratios, the odd–odd modes which correspond to the in-phase oscillatory motion about the midplane between the cylinders are the first to become unstable as compared with the odd–even modes, and with much lower frequencies.


2002 ◽  
Vol 461 ◽  
pp. 365-386 ◽  
Author(s):  
DONGJOO KIM ◽  
HAECHEON CHOI

Numerical simulations are conducted for laminar flow past a sphere rotating in the streamwise direction, in order to investigate the effect of the rotation on the characteristics of flow over the sphere. The Reynolds numbers considered are Re = 100, 250 and 300 based on the free-stream velocity and sphere diameter, and the rotational speeds are in the range of 0 [les ] ω* [les ] 1, where ω* is the maximum azimuthal velocity on the sphere surface normalized by the free-stream velocity. At ω* = 0 (without rotation), the flow past the sphere is steady axisymmetric, steady planar-symmetric, and unsteady planar-symmetric, respectively, at Re = 100, 250 and 300. Thus, the time-averaged lift forces exerted on the stationary sphere are not zero at Re = 250 and 300. When the rotational speed increases, the time-averaged drag force increases for the Reynolds numbers investigated, whereas the time-averaged lift force is zero for all ω* > 0. On the other hand, the lift force fluctuations show a non-monotonic behaviour with respect to the rotational speed. At Re = 100, the flow past the sphere is steady axisymmetric for all the rotational speeds considered and thus the lift force fluctuation is zero. At Re = 250 and 300, however, the flows are unsteady with rotation and the lift force fluctuations first decrease and then increase with increasing rotational speed, showing a local minimum at a specific rotational speed. The vortical structures behind the sphere are also significantly modified by the rotation. For example, at Re = 300, the flows become ‘frozen’ at ω* = 0.5 and 0.6, i.e. the vortical structures in the wake simply rotate without temporal variation of their strength and the magnitude of the instantaneous lift force is constant in time. It is shown that the flow becomes frozen at higher rotational speed with increasing Reynolds number. The rotation speed of the vortical structures is shown to be slower than that of the sphere.


Sign in / Sign up

Export Citation Format

Share Document