Behaviour of macroscopic rigid spheres in lid-driven cavity flow

Experiments are conducted to investigate the behaviour of macroscopic rigid particles suspended in a fully three-dimensional viscous flow. The flow considered takes place in a closed cubic cavity, steadily driven along its upper face by a translating lid. Navier–Stokes computations are first performed to characterize the fluid flow, and the resulting kinematic template is checked using laser-illuminated micro-particles. Nearly neutrally buoyant rigid spheres are then inserted in the cavity, and their three-dimensional motions are tracked using stereoscopic imaging. The measured macro-particle motions are compared with those of simulated passive tracers, and their responses to changes in experimental conditions are examined. Although steric effects are observed to hinder passage through narrow throats of the flow field, macro-particle trajectories are otherwise found to align closely with passive tracer paths. The macro-particle orbits, however, are not evenly distributed within the cavity, and cluster closer to the periphery as the Reynolds and Stokes numbers increase. With support from observations of particle rotations relative to the fluid, we interpret this behaviour as resulting from weak forces pulling the macroscopic spheres towards preferential paths, similar to the Segré–Silberberg effect in Poiseuille flow.     

Effects of shear thinning on migration of neutrally buoyant particles in pressure driven flow of Newtonian and viscoelastic fluids

The pattern of cross stream migration of neutrally buoyant particles in a pressure driven flow depends strongly on the properties of the suspending fluid. These migration effects have been studied by direct numerical simulation in planar flow. Shear thinning has a large effect when the inertia or elasticity is large, but only a small effect when they are small. At moderate Reynolds numbers, shear thinning causes particles to migrate away from the centerline, creating a particle-free zone in the core of the channel, which increases with the amount of shear thinning. In a viscoelastic fluid with shear thinning, particles migrate either toward the centerline or toward the walls, creating an annular particle-free zone at intermediate radii. The simulations also give rise to precise determination of slip velocity distributions in the various cases studied.     

Direct simulation of the motion of neutrally buoyant balls in a three-dimensional Poiseuille flow

In a previous article the authors introduced a Lagrange multiplier based fictitious domain method. Their goal in the present article is to apply a generalization of the above method to: (i) the numerical simulation of the motion of neutrally buoyant particles in a three-dimensional Poiseuille flow; (ii) study – via direct numerical simulations – the migration of neutrally buoyant balls in the tube Poiseuille flow of an incompressible Newtonian viscous fluid. Simulations made with one and several particles show that, as expected, the Segré–Silberberg effect takes place. To cite this article: T.-W. Pan, R. Glowinski, C. R. Mecanique 333 (2005).     

Time-resolved PIV measurements of the interaction of polystyrene beads with near-wall-coherent structures in a turbulent channel flow

Time-resolved PIV measurements were performed in a dilute particle-laden flow tracking near-neutrally buoyant polystyrene beads and the velocity field of a near wall turbulent boundary layer. Data were taken in a vertical light sheet aligned in the streamwise direction at the center of a horizontal, closed loop, transparent square water channel facility. In addition, low speed measurements were performed characterizing the effects of the dispersed phase on mean and turbulence flow quantities. Reynolds shear stress slightly differed from clear water conditions whereas fluid mean and rms values were not affected. A case study for several beads revealed a clear relation between their movement and near-wall coherent structures. Several structures having 2D vorticity signatures of near-wall hairpin vortices and hairpin packets, directly affected bead movement. A statistical analysis showed that the mean streamwise velocity of ascending beads lagged behind the mean fluid velocity and bead rms values were higher than fluid ones. Particle Reynolds numbers based on the magnitude of the instantaneous relative velocity vector peaked near the wall; values not exceeding 100, too low for vortex shedding to occur. Quadrant analysis showed a clear preference for ascending beads to reside in ejections while for descending beads the preference for sweeps was less.     

Three-dimensional magnetic resonance velocimetry measurements of turbulence quantities in complex flow

A magnetic resonance velocimetry (MRV) experimental technique based on magnetic resonance imaging and capable of measuring the turbulent Reynolds stresses in a 3D flow domain is described. Results are presented in backward facing step flow in a square channel with a Reynolds number of 48,000 based on step height and freestream velocity at the step. MRV results are compared to particle image velocimetry (PIV) measurements in the centerplane containing the streamwise and cross-stream axes. MRV and PIV mean velocity measurements show excellent agreement. MRV measurements for Reynolds normal stresses compare to within ±20% of the PIV results while results for the turbulent shear are less accurate.     

Flow over periodic hills: an experimental study

Two-dimensional flow over periodically arranged hills was investigated experimentally in a water channel. Two-dimensional particle image velocimetry (PIV) and one-dimensional laser Doppler anemometry (LDA) measurements were undertaken at four Reynolds numbers ( \text5,600 ￡ Re ￡ \text37,000\text{5,600} \le Re \le \text{37,000}). Two-dimensional PIV field measurements were thoroughly validated by means of point-by-point 1D LDA measurements at certain positions of the flow. A detailed study of the periodicity and the homogeneity was undertaken, which demonstrates that the flow can be regarded as two-dimensional and periodic for Re 3 \text10,000Re \ge \text{10,000}. We found a decreasing reattachment length with increasing Reynolds number. This is connected to a higher momentum in the near-wall zone close to flow separation which comes from the velocity speed up above the obstacle. This leads to a velocity overshoot directly above the hill crest which increases with Reynolds number as the inner layer depth decreases. The flow speed up above that layer is independent of the Reynolds number which supports the assumption of inviscid flow disturbance in the outer layer usually made in asymptotic theory for flow over small hills.     

Low-Reynolds-number motion of a deformable drop between two parallel plane walls

The motion of a three-dimensional deformable drop between two parallel plane walls in a low-Reynolds-number Poiseuille flow is examined using a boundary-integral algorithm that employs the Green’s function for the domain between two infinite plane walls, which incorporates the wall effects without discretization of the walls. We have developed an economical calculation scheme that allows long-time dynamical simulations, so that both transient and steady-state shapes and velocities are obtained. Results are presented for neutrally buoyant drops having various viscosity, size, deformability, and channel position. For nearly spherical drops, the decrease in translational velocity relative to the undisturbed fluid velocity at the drop center increases with drop size, proximity of the drop to one or both walls, and drop-to-medium viscosity ratio. When deformable drops are initially placed off the centerline of flow, lateral migration towards the channel center is observed, where the drops obtain steady shapes and translational velocities for subcritical capillary numbers. With increasing capillary number, the drops become more deformed and have larger steady velocities due to larger drop-to-wall clearances. Non-monotonic behavior for the lateral migration velocities with increasing viscosity ratio is observed. Simulation results for large drops with non-deformed spherical diameters exceeding the channel height are also presented.     

Apparent wall slip velocity measurements in free surface flow of concentrated suspensions

In this work we have experimentally measured the apparent wall slip velocity in open channel flow of neutrally buoyant suspension of non-colloidal particles. The free surface velocity profile was measured using the tool of particle imaging velocimetry (PIV) for two different channels made of plane and rough walls. The rough walled channel prevents wall slip, whereas the plane wall showed significant wall slip due to formation of slip layer. By comparing the velocity profiles from these two cases we were able to determine the apparent wall slip velocity. This method allows characterization of wall slip in suspension of large sized particles which cannot be performed in conventional rheometers. Experiments were carried out for concentrated suspensions of various particle volume concentrations and for two different sizes of particles. It was observed that wall slip velocity increases with particle size and concentration but decreases with increase in the viscosity of suspending fluid. The apparent wall slip velocity coefficients are in qualitative agreement with the earlier measurements. The effect of wall slip on free surface corrugation was also studied by analyzing the power spectral density (PSD) of the refracted light from the free surface. Our results indicate that free surface corrugation is a bulk flow response and it does not arise from boundary problem such as development of slip layer.     

Free-surface flow of concentrated suspensions

Free-surface flows of concentrated suspensions exhibit many interesting phenomena such as particle segregation and surface corrugation. In this work the flow structures associated with free-surface has been studied experimentally. The free-surface velocity for neutrally buoyant suspension of uniform spheres in a gravity driven inclined channel flow was determined by particle imaging velocimetry (PIV) technique. Experiments were carried out for concentrated suspensions with particle fractions ? ranging from 0.40 to 0.50. The measured velocities show blunted profile in the channel. The blunting of the velocity profile increases with the particle concentration. The rms velocity fluctuations measured at the free-surface progressively increase with particle fraction ? and are linear in shear rate γ. The surface roughness were characterized by analyzing the power spectral density of the refracted light from the free-surface. The characteristics observed are in support of earlier findings.     

Numerical simulation of the motion of neutrally buoyant particles in the plane Poiseuille flow of a Newtonian fluid

In this note we discuss the generalization of a Lagrange multiplier based fictitious domain method to the simulation of migration of neutrally buoyant particles in plane Poiseuille flow of a Newtonian fluid. The migration away from the center of the channel is believed to be an effect of the curvature of the velocity profile. We found this effect is not weakened by the presence of many particles, but instead by the collisions among the particles. Experiments and simulations show that the particles concentrate in the central region where the shear rate is low.     

A new facility for time-resolved PIV measurements in rotating channels

A new facility to measure the time evolution of 2D velocity fields in a rotating channel is presented, and the accuracy is discussed in detail. Measurements are made by means of a time-resolved PIV system composed of a continuous laser diode, coupled by a fiber optics cable to a laser plane optical module, and a CMOS high-speed camera. Both the PIV system and divergent channel are fixed on a 2.5 m rotating disk. This allows a direct measurement of the relative velocity of flows with Reynolds numbers between 3 × 10³ and 3 × 10⁴ and Rotation numbers between 0.0 and 0.52. These values correspond to the flow conditions in small radial impellers and can be independently adjusted by a change of the relative flow velocity and RPM. It is shown that this new facility allows high signal-to-noise ratios, and that the direct acquisition of the data in a rotating frame drastically reduces the measurement error. The accuracy and high spatial and temporal resolution of the measurements allow a detailed analysis of boundary layer characteristics in stationary and rotating conditions.     

PIV–LES analysis of channel flow rotating about the streamwise axis

The flow field of a channel rotating about the streamwise axis is analyzed experimentally and numerically. The current investigations were carried out at a bulk velocity based Reynolds number of Re_m = 2850 and a friction velocity based Reynolds number of Re_τ = 180, respectively. Particle-image velocimetry (PIV) measurements are compared with large-eddy simulation data to show earlier direct numerical simulation findings to generate too large a reverse flow region in the center region of the spanwise flow. The development of the mean spanwise velocity distribution and the influence of the rotation on the turbulent properties, i.e., the Reynolds stresses and the two-point correlations of the flow, are confirmed in both investigations. The rotation primarily influences those components of the Reynolds shear stresses, which contain the spanwise velocity component. The size of the correlation areas and thus the length scales of the flow generally grow in all three coordinate directions leading to longer structures. Furthermore, experimental results of the same channel flow at a significantly lower bulk Reynolds number of Re_{m, l} = 665, i.e., a laminar flow in a non-rotating channel, are introduced. The experiments show the low Reynolds number flow to become turbulent under rotation and to develop the same characteristics as the high Reynolds number flow.     

Experimental investigation of turbulent natural convection flow in a channel

This paper presents the results of velocity measurements of natural convection in symmetrically heated vertical channel using the particle image velocimetry (PIV) system. Velocity measurements were conducted at three different sections on the horizontal plane to validate the flow two-dimensionality and at three different heights in the vertical plane to establish vertical mean velocity profiles. The results indicate a considerable influence of the Rayleigh number and aspect ratio on the mean velocity profile. The results also indicate significant diffusion rates of the vertical mean velocity component and normal Reynolds stresses towards the center of the channel.     

The lateral migration of solid particles in a laminar flow near a plane

The lateral inertial migration of a solid spherical particle suspended in a laminar flow over a vertical wall is considered theoretically. Formulae for the migration velocity are obtained for both neutrally buoyant and non-neutrally buoyant particles and also for the case of zero flow over the wall. Situations in which the particle is either free to rotate or prevented from rotating are considered. Such results are found to agree qualitatively with known experimental data.     

Pipe flow of suspensions

This study shows that fully developed pipe flow of a particulate suspension is defined by four dimensionless parameters of particle-fluid interactions in addition to the Reynolds number. Effects accounted for include the Magnus effect due to fluid shear, electrostatic repulsion due to electric charges on the particles, and Brownian or turbulent diffusion. In the case of a laminar liquid-solid suspension, electrostatic effect is negligible but shear effect is prominent. Solution of the basic equations gives the density distribution of particles with a peak at the center (Einstein, Jeffery) or at other radii between the center and the pipe wall (Segré et al) depending on the magnitudes of the various flow parameters. In the case of a turbulent gas-solid suspension, the Magnus effect is significant only within the thickness of the laminar sublayer. However, charges induced on the particles by the impact of particles at the wall produce a higher density at the wall than at the center of the pipe. The velocity distribution of particles is characterized by a slip velocity at the wall and a lag in velocity in the core from the fluid phase. These results are verified by earlier measurements.     

An automated single-particle tracker: application to characterization of non-azimuthal motion in Couette flows at low Reynolds number

We describe an experiment that allows us to record 3-dimensional trajectories of single particles in Couette shear flows, at low Reynolds number. The core of the apparatus is a Couette cell with transparent contra-rotating cylinders. Fluorescent spherical particles are used as tracers. A single tracer is imaged onto a webcam, equipped with a home-made autofocus system. For a given average shear rate, tracking of an individual tracer is performed automatically by driving the amount of contra-rotation between both cylinders and the position of the webcam. The performance of the tracker is illustrated through examples of trajectories of neutrally buoyant tracers in a Newtonian fluid. The setup is mostly aimed at characterizing complex flows in non-colloidal concentrated suspensions and wet granular materials. We show examples of 3d trajectories in a dense suspension of 200 μm spherical grains, revealing details of the short-scale diffusive-like particle motion, together with flow localization and large-scale non-azimuthal flow patterns.     

PIV measurements in a weakly separating and reattaching turbulent boundary layer

A high Reynolds number flat plate turbulent boundary layer was studied in a wind-tunnel experiment using particle image velocimetry (PIV). The flow is subjected to an adverse pressure gradient (APG) which is designed such that the boundary layer separates and reattaches, forming a weak separation bubble. With PIV we are able to get a more complete picture of this complex flow phenomenon. The view of a separation bubble being composed of large scale coherent regions of instantaneous backflow occurring randomly in a three-dimensional manner in space and time is verified by the present PIV measurements. The PIV database was used to test the applicability of various velocity scalings around the separation bubble. We found that the mean velocity profiles in the outer part of the boundary layer, and to some extent also the Reynolds shear-stress, are self-similar when using a velocity scale based on the local pressure gradient. The same can be said for the so called Perry–Schofield scaling, which suggests that the two velocity scales are connected. This can also be interpreted as an experimental evidence of the claimed relation between the latter velocity scale and the maximum Reynolds shear-stress.     

Experimental measurement of flow past cavities of different shapes

Experiments are carried out in order to investigate the flow structure past a rectangular, triangular and semi-circular cavity of length-to-depth ratio of 2:1 using the Particle Image Velocimetry (PIV) technique. The experiments are performed in a large scale water channel with three different upstream velocities resulting in Reynolds numbers of 1230, 1460 and 1700, based on inflow momentum thickness, for each cavity type. Contours of constant averaged streamwise and transverse components of velocity, contours of constant averaged vorticity, Reynolds stress and streamline plots for each cavity type for the aforementioned three Reynolds numbers are presented. In addition, streamwise velocity, Reynolds stress and turbulence intensity values are compared for all cavity types. Effect of cavity shape on flow structure within the cavity is discussed in detail. Moreover, spectrum of instantaneous streamwise velocity fluctuations in shear layer near the downstream of the leading corner and the upstream of the trailing corner of the cavities are obtained and it was found that no organized oscillations are present in the flow; rectangular and triangular cavities have the largest amplitudes while semi-circular cavity has the smallest. Calculated turbulence intensities also reveal that the maximum turbulence intensities occur at cavity lid in the centerline section and rectangular and triangular cavities have larger turbulence intensity compared to semi-circular cavity.     

The circular cylinder in simple shear at moderate Reynolds numbers: An experimental study

The flow around a neutrally buoyant circular cylinder in simple shear was studied using particle-image velocimetry. Steady experimental results are presented for both torque-free and fixed cylinders at shear-based Reynolds numbers between 8 and 25 at two different flow confinement values. These data, which are to our knowledge the first experimental results for this Reynolds number range, show that the normalized rotation rate decays much more slowly with Reynolds number than previous theoretical predictions for unconfined flow, with the rate of decay decreasing as confinement increases. A new simple model is proposed to explain these effects. Given physical and technological limitations upon steady and unconfined flow experiments and simulations, the effect of confinement on a cylinder in simple shear merits further study. Received: 17 February 2000/Accepted: 22 June 2000     

Cross-stream forces and velocities of fixed and freely suspended particles in viscoelastic Poiseuille flow: Perturbation and numerical analyses

The cross-stream migration of a circular particles (or infinitely long cylinder) in two dimensional, inertia-less viscoelastic pressure-driven flows is examined through complementary finite element simulations and second-order fluid perturbation analyses for small Deborah number (De), where De is defined as the fluid relaxation time divided by the characteristic flow time. A neutrally buoyant, freely suspended particle migrates toward the center of the channel for all particle sizes and cross-stream positions due to the coupled effects of the linear and quadratic variations of the imposed velocity. A particle that is held at a fixed position, in contrast, experiences a cross-stream force directed toward the wall as a result of the coupled effects of the local shear flow and the flow relative to the particle.