Turbulence attenuation by large neutrally buoyant particles

Turbulence modulation by inertial-range-size, neutrally-buoyant particles is investigated experimentally in a von K\'arm\'an flow. Increasing the particle volume fraction $\Phi_\mathrm{v}$, maintaining constant impellers Reynolds number attenuates the fluid turbulence. The inertial-range energy transfer rate decreases as $\propto\Phi_\mathrm{v}^{2/3}$, suggesting that only particles located on a surface affect the flow. Small-scale turbulent properties, such as structure functions or acceleration distribution, are unchanged. Finally, measurements hint at the existence of a transition between two different regimes occurring when the average distance between large particles is of the order of the thickness of their boundary layers.Comment: 7 pages, 4 figure

Rotational intermittency and turbulence induced lift experienced by large particles in a turbulent flow

The motion of a large, neutrally buoyant, particle, freely advected by a turbulent flow is determined experimentally. We demonstrate that both the translational and angular accelerations exhibit very wide probability distributions, a manifestation of intermittency. The orientation of the angular velocity with respect to the trajectory, as well as the translational acceleration conditioned on the spinning velocity provide evidence of a lift force acting on the particle.Comment: 4 page, 4 figure

Droplet size distribution in homogeneous isotropic turbulence

We study the physics of droplet breakup in a statistically stationary homogeneous and isotropic turbulent flow by means of high resolution numerical investigations based on the multicomponent lattice Boltzmann method. We verified the validity of the criterion proposed by Hinze (1955) for droplet breakup and we measured the full probability distribution function (pdf) of droplets radii at different Reynolds numbers and for different volume fraction. By means of a Lagrangian tracking we could follow individual droplets along their trajectories, define a local Weber number based on the velocity gradients and study its cross-correlation with droplet deformation.Comment: 10 pages, 6 figure

Velocity gradients statistics along particle trajectories in turbulent flows: the refined similarity hypothesis in the Lagrangian frame

We present an investigation of the statistics of velocity gradient related quantities, in particluar energy dissipation rate and enstrophy, along the trajectories of fluid tracers and of heavy/light particles advected by a homogeneous and isotropic turbulent flow. The Refined Similarity Hypothesis (RSH) proposed by Kolmogorov and Oboukhov in 1962 is rephrased in the Lagrangian context and then tested along the particle trajectories. The study is performed on state-of-the-art numerical data resulting from numerical simulations up to Re~400 with 2048^3 collocation points. When particles have small inertia, we show that the Lagrangian formulation of the RSH is well verified for time lags larger than the typical response time of the particle. In contrast, in the large inertia limit when the particle response time approaches the integral-time-scale of the flow, particles behave nearly ballistic, and the Eulerian formulation of RSH holds in the inertial-range.Comment: 7 pages, 7 figures; Physical Review E (accepted Dec 7, 2009

Simultaneous 3D measurement of the translation and rotation of finite size particles and the flow field in a fully developed turbulent water flow

We report a novel experimental technique that measures simultaneously in three dimensions the trajectories, the translation, and the rotation of finite size inertial particles together with the turbulent flow. The flow field is analyzed by tracking the temporal evolution of small fluorescent tracer particles. The inertial particles consist of a super-absorbent polymer that renders them index and density matched with water and thus invisible. The particles are marked by inserting at various locations tracer particles into the polymer. Translation and rotation, as well as the flow field around the particle are recovered dynamically from the analysis of the marker and tracer particle trajectories. We apply this technique to study the dynamics of inertial particles much larger in size (Rp/{\eta} \approx 100) than the Kolmogorov length scale {\eta} in a von K\'arm\'an swirling water flow (R{\lambda} \approx 400). We show, using the mixed (particle/fluid) Eulerian second order velocity structure function, that the interaction zone between the particle and the flow develops in a spherical shell of width 2Rp around the particle of radius Rp. This we interpret as an indication of a wake induced by the particle. This measurement technique has many additional advantages that will make it useful to address other problems such as particle collisions, dynamics of non-spherical solid objects, or even of wet granular matter.Comment: 18 pages, 7 figures, submitted to "Measurement Science and Technology" special issue on "Advances in 3D velocimetry

Effect of energy spectrum law on clustering patterns for inertial particles subjected to gravity in Kinematic Simulation

We study the clustering of inertial particles using a periodic kinematic simulation. Particles clustering is observed for different pairs of Stokes number and Froude number and different spectral power laws (1.4 6 p 6 2.1). The main focus is to identify and then quantify the effect of p on the clustering attractor - by attractor we mean the set of points in the physical space where the particles settle when time tends to infinity. It is observed that spectral power laws can have a dramatic effect on the attractor shape. In particular, we observed a new attractor type which was not present in previous studies for Kolmogorov spectra (p = 5/3)

Convergent Chaos

Chaos is widely understood as being a consequence of sensitive dependence upon initial conditions. This is the result of an instability in phase space, which separates trajectories exponentially. Here, we demonstrate that this criterion should be refined. Despite their overall intrinsic instability, trajectories may be very strongly convergent in phase space over extremely long periods, as revealed by our investigation of a simple chaotic system (a realistic model for small bodies in a turbulent flow). We establish that this strong convergence is a multi-facetted phenomenon, in which the clustering is intense, widespread and balanced by lacunarity of other regions. Power laws, indicative of scale-free features, characterize the distribution of particles in the system. We use large-deviation and extreme-value statistics to explain the effect. Our results show that the interpretation of the 'butterfly effect' needs to be carefully qualified. We argue that the combination of mixing and clustering processes makes our specific model relevant to understanding the evolution of simple organisms. Lastly, this notion of convergent chaos, which implies the existence of conditions for which uncertainties are unexpectedly small, may also be relevant to the valuation of insurance and futures contracts

Analysis of the correlation dimension for inertial particles

We obtain an implicit equation for the correlation dimension which describes cluster- ing of inertial particles in a complex flow onto a fractal measure. Our general equation involves a propagator of a nonlinear stochastic process in which the velocity gradient of the fluid appears as additive noise. When the long-time limit of the propagator is considered our equation reduces to an existing large-deviation formalism from which it is difficult to extract concrete results. In the short-time limit, however, our equation reduces to a solvability condition on a partial differential equation. In the case where the inertial particles are much denser than the fluid, we show how this approach leads to a perturbative expansion of the correlation dimension, for which the coefficients can be obtained exactly and in principle to any order. We derive the perturbation series for the correlation dimension of inertial particles suspended in three-dimensional spatially smooth random flows with white-noise time correlations, obtaining the first 33 non-zero coefficients exactly

PDF model based on Langevin equation for polydispersed two-phase flows applied to a bluff-body gas-solid flow,

The aim of the paper is to discuss the main characteristics of a complete theoretical and numerical model for turbulent polydispersed two-phase flows, pointing out some specific issues. The theoretical details of the model have already been presented [Minier and Peirano, Physics Reports, Vol. 352/1-3, 2001 ]. Consequently, the present work is mainly focused on complementary aspects, that are often overlooked and that require particular attention. In particular, the following points are analysed : the necessity to add an extra term in the equation for the velocity of the fluid seen in the case of twoway coupling, the theoretical and numerical evaluations of particle averages and the fulfilment of the particle mass-continuity constraint. The theoretical model is developed within the PDF formalism. The important-physical choice of the state vector variables is first discussed and the model is then expressed as a stochastic differential equation (SDE) written in continuous time (Langevin equations) for the velocity of the fluid seen. The interests and limitations of Langevin equations, compared to the single-phase case, are reviewed. From the numerical point of view, the model corresponds to an hybrid Eulerian/Lagrangian approach where the fluid and particle phases are simulated by different methods. Important aspects of the Monte Carlo particle/mesh numerical method are emphasised. Finally, the complete model is validated and its performance is assessed by simulating a bluff-body case with an important recirculation zone and in which two-way coupling is noticeable.Comment: 23 pages, 10 figure

An Euler Solver Based on Locally Adaptive Discrete Velocities

A new discrete-velocity model is presented to solve the three-dimensional Euler equations. The velocities in the model are of an adaptive nature---both the origin of the discrete-velocity space and the magnitudes of the discrete-velocities are dependent on the local flow--- and are used in a finite volume context. The numerical implementation of the model follows the near-equilibrium flow method of Nadiga and Pullin [1] and results in a scheme which is second order in space (in the smooth regions and between first and second order at discontinuities) and second order in time. (The three-dimensional code is included.) For one choice of the scaling between the magnitude of the discrete-velocities and the local internal energy of the flow, the method reduces to a flux-splitting scheme based on characteristics. As a preliminary exercise, the result of the Sod shock-tube simulation is compared to the exact solution.Comment: 17 pages including 2 figures and CMFortran code listing. All in one postscript file (adv.ps) compressed and uuencoded (adv.uu). Name mail file `adv.uu'. Edit so that `#!/bin/csh -f' is the first line of adv.uu On a unix machine say `csh adv.uu'. On a non-unix machine: uudecode adv.uu; uncompress adv.tar.Z; tar -xvf adv.ta