scholarly journals In Search of Random Uncorrelated Particle Motion (RUM) in a Simple Random Flow Field

Author(s):  
Michael W. Reeks ◽  
Luca Fabbro ◽  
Alfredo Soldati

DNS studies of dispersed particle motion in isotropic homogeneous turbulence [1] have revealed the existence of a component of random uncorrelated motion (RUM) dependent on the particle inertia τp (normalised particle response time or Stoke number). This paper reports the presence of RUM in a simple linear random smoothly varying flow field of counter rotating vortices where the two-particle velocity correlation was measured as a function of spatial separation. Values of the correlation less than one for zero separation indicated the presence of RUM. In terms of Stokes number, the motion of the particles in one direction corresponds to either a heavily damped (τp < 0.25) or lightly damped (τp > 0.25) harmonic oscillator. In the lightly damped case the particles overshoot the stagnation lines of the flow and are projected from one vortex to another (the so-called sling-shot effect). It is shown that RUM occurs only when τp > 0.25, increasing monotonically with increasing Stokes number. Calculations of the particle pair separation distribution function show that equilibrium of the particle concentration field is never reached, the concentration at zero separation increasing monotonically with time. This is consistent with the calculated negative values of the average Liapounov exponent (finite compressibility) of the particle velocity field.

Author(s):  
Michael W. Reeks

Comparisons are made between the Advection-Diffusion Equation (ADE) approach for particle transport and the two fluid model approach based on the PDF method. In principal the ADE approach offers a simpler way of calculating the inertial deposition of particles in a turbulent boundary layer than that based on the PDF approach. However the ADE equations that have recently been used are only strictly valid for a simple Gaussian process when particle inertia is small. Using a prescribed but in general non-Gaussian random particle velocity field, it is shown that the net particle mass flux contains an extra drift term to that from the mean velocity of the particle velocity field, associated with the compressibility of the velocity field. Furthermore the diffusive flux in general depends not only upon the gradient of the mean concentration (true only for a Gaussian random flow field) but also upon higher order derivatives whose relative contribution depends on diffusion coefficients Dijk... etc. These coefficients depend upon the statistical moments associated with random displacements and compressibility of the particle flow field along particle trajectories which in turn depend upon particle inertia. In contrast the PDF approach offers the advantage of using a simple gradient (Gaussian) approximation in particle phase space which can lead to a non-Gaussian spatial dispersion process when particle inertia is important. Conditions based on the particle mean free path are derived for which a simple ADE is appropriate. Some of the features of particle transport in an inhomogeneous turbulent flow are illustrated by examining particle dispersion in a random flow field composed of pairs of counter rotating vortices which has an rms velocity which increase linearly from a stagnation point.


2012 ◽  
Vol 694 ◽  
pp. 426-463 ◽  
Author(s):  
Marco Martins Afonso ◽  
Andrea Mazzino ◽  
Paolo Muratore-Ginanneschi

AbstractThe large-scale/long-time transport of inertial particles of arbitrary mass density under gravity is investigated by means of a formal multiple-scale perturbative expansion in the scale-separation parameter between the carrier flow and the particle concentration field. The resulting large-scale equation for the particle concentration is determined, and is found to be diffusive with a positive definite eddy diffusivity. The calculation of the latter tensor is reduced to the resolution of an auxiliary differential problem, consisting of a coupled set of two differential equations in a $(6+ 1)$-dimensional coordinate system (three space coordinates plus three velocity coordinates plus time). Although expensive, numerical methods can be exploited to obtain the eddy diffusivity, for any desirable non-perturbative limit (e.g. arbitrary Stokes and Froude numbers). The aforementioned large-scale equation is then specialized to deal with two different relevant perturbative limits: (i) vanishing of both Stokes time and sedimenting particle velocity; (ii) vanishing Stokes time and finite sedimenting particle velocity. Both asymptotics lead to a greatly simplified auxiliary differential problem, now involving only space coordinates and thus easily tackled by standard numerical techniques. Explicit, exact expressions for the eddy diffusivities have been calculated, for both asymptotics, for the class of parallel flows, both static and time-dependent. This allows us to investigate analytically the role of gravity and inertia on the diffusion process by varying relevant features of the carrier flow, such as the form of its temporal correlation function. Our results exclude a universal role played by gravity and inertia on the diffusive behaviour: regimes of both enhanced and reduced diffusion may exist, depending on the detailed structure of the carrier flow.


2011 ◽  
Vol 686 ◽  
pp. 338-351 ◽  
Author(s):  
Elena Meneguz ◽  
Michael W. Reeks

AbstractA full Lagrangian method (FLM) is used in direct numerical simulations (DNS) of incompressible homogeneous isotropic and statistically stationary turbulent flow to measure the statistical properties of the segregation of small inertial particles advected with Stokes drag by the flow. Qualitative good agreement is observed with previous kinematic simulations (KS) (IJzermans, Meneguz & Reeks, J. Fluid Mech., vol. 653, 2010, pp. 99–136): in particular, the existence of singularities in the particle concentration field and a threshold value for the particle Stokes number $\mathit{St}$ above which the net compressibility of the particle concentration changes sign (from compression to dilation). A further KS analysis is carried out by examining the distribution in time of the compression of an elemental volume of particles, which shows that it is close to Gaussian as far as the third and fourth moments but non-Gaussian (within the uncertainties of the measurements) for higher-order moments when the contribution of singularities in the tails of the distribution increasingly dominates the statistics. Measurements of the rate of occurrence of singularities show that it reaches a maximum at $\mathit{St}\ensuremath{\sim} 1$, with the distribution of times between singularities following a Poisson process. Following the approach used by Fevrier, Simonin & Squires (J. Fluid Mech., vol. 553, 2005, pp. 1–46), we also measured the random uncorrelated motion (RUM) and mesoscopic components of the compression for $\mathit{St}= 1$ and show that the non-Gaussian highly intermittent part of the distribution of the compression is associated with the RUM component and ultimately with the occurrence of singularities. This result is consistent with the formation of caustics (Wilkinson et al. Phys. Fluids, vol. 19, 2007, p. 113303), where the activation of singularities precedes the crossing of trajectories (RUM).


2000 ◽  
Vol 421 ◽  
pp. 185-227 ◽  
Author(s):  
E. MEIBURG ◽  
E. WALLNER ◽  
A. PAGELLA ◽  
A. RIAZ ◽  
C. HÄRTEL ◽  
...  

The two-way coupling mechanisms in particle-laden mixing layers are investigated, with and without particle settling, and with an emphasis on the resulting modifications to the fluid vorticity field. The governing equations are interpreted with respect to the production and cancellation of vorticity. These mechanisms are shown to be related to the misalignment of the concentration gradient and the slip velocity, as well as to the difference in fluid and particle vorticities. Preliminary insight into the physics is obtained from an analysis of the unidirectional base flow. For this model problem, the conditions are established under which the particle velocity remains a single-valued function of space for all times. The resulting simplified set of two-way-coupled equations governing the vorticity of the fluid and particulate phases, respectively, is solved numerically. The formation of a decaying travelling wave solution is demonstrated over a wide range of parameters. Interestingly, the downward propagation of the fluid vorticity field is not accomplished through convection, but rather by the production and loss of vorticity on opposite sides of the mixing layer. For moderate settling velocities, the simulation results reveal an optimal coupling mechanism between the fluid and particle vorticities at intermediate values of the mass loading parameter. For large settling velocities and intermediate mass loadings, more than one local maximum is seen to evolve in the vorticity field. A scaling law for the downward propagation rates of the vorticity fronts is derived.Two-dimensional particle-laden mixing layers are investigated by means of a mixed Lagrangian–Eulerian approach which is based on the vorticity variable. For uniformly seeded mixing layers, the simulations confirm some of the features observed by Druzhinin (1995b) for the model problem of a two-way-coupled particle-laden Stuart vortex, as well as by Dimas & Kiger (1998) in a linear stability analysis. For small values of the Stokes number, a mild destabilization of the mixing layer is observed. At moderate and large Stokes numbers, on the other hand, the transport of vorticity from the braids into the core of the evolving Kelvin–Helmholtz vortices is seen to be slowed by the two-way coupling effects. As a result, the particle ejection from the vortex cores is weakened. For constant mass loadings, the two-way coupling effects are strongest at intermediate Stokes number values. For moderately large Stokes numbers, the formation of two bands of high particle concentration is observed in the braids, which reflects the multi-valued nature of the particle velocity field. For mixing layers in which only one stream is seeded, the particle concentration gradient across the mixing layer leads to strong vorticity production and loss, which results in an effective net motion of the vortex in the flow direction of the seeded stream. Under particle settling, the vortex propagates downward as well. For the parameter range explored here, its settling velocity agrees well with the scaling law derived from the unidirectional flow analysis.


2010 ◽  
Vol 653 ◽  
pp. 99-136 ◽  
Author(s):  
RUTGER H. A. IJZERMANS ◽  
ELENA MENEGUZ ◽  
MICHAEL W. REEKS

The results presented here are part of a long-term study in which we analyse the segregation of inertial particles in turbulent flows using the so called full Lagrangian method (FLM) to evaluate the ‘compressibility’ of the particle phase along a particle trajectory. In the present work, particles are advected by Stokes drag in a random flow field consisting of counter-rotating vortices and in a flow field composed of 200 random Fourier modes. Both flows are incompressible and, like turbulence, have structure and a distribution of scales with finite lifetime. The compressibility is obtained by first calculating the deformation tensor Jij associated with an infinitesimally small volume of particles following the trajectory of an individual particle. The fraction of the initial volume occupied by the particles centred around a position x at time t is denoted by |J|, where J ≡ det(Jij) and Jij ≡ ∂xi(x0, t)/∂x0,j, x0 denoting the initial position of the particle. The quantity d〈ln|J|〉/dt is shown to be equal to the particle averaged compressibility of the particle velocity field 〈∇ · v〉, which gives a measure of the rate-of-change of the total volume occupied by the particle phase as a continuum. In both flow fields the compressibility of the particle velocity field is shown to decrease continuously if the Stokes number St (the dimensionless particle relaxation time) is below a threshold value Stcr, indicating that the segregation of particles continues indefinitely. We show analytically and numerically that the long-time limit of 〈∇ · v〉 for sufficiently small values of St is proportional to St2 in the flow field composed of random Fourier modes, and to St in the flow field consisting of counter-rotating vortices. If St > Stcr, however, the particles are ‘mixed’. The level of mixing can be quantified by the degree of random uncorrelated motion (RUM) of particles which is a measure of the decorrelation of the velocities of two nearby particles. RUM is zero for fluid particles and increases rapidly with the Stokes number if St > Stcr, approaching unity for St ≫ 1. The spatial averages of the higher-order moments of the particle number density are shown to diverge with time indicating that the spatial distribution of particles may be very intermittent, being associated with non-zero values of RUM and the occurrence of singularities in the particle velocity field. Our results are consistent with previous observations of the radial distribution function in Chun et al. (J. Fluid Mech., vol. 536, 2005, p. 219).


Fluids ◽  
2021 ◽  
Vol 6 (1) ◽  
pp. 27
Author(s):  
J. Barry Greenberg ◽  
David Katoshevski

A theoretical investigation of the influence of a standing wave flow-field on the dynamics of a laminar two-dimensional spray diffusion flame is presented for the first time. The mathematical analysis permits mild slip between the droplets and their host surroundings. For the liquid phase, the use of a small Stokes number as the perturbation parameater enables a solution of the governing equations to be developed. Influence of the standing wave flow-field on droplet grouping is described by a specially constructed modification of the vaporization Damkohler number. Instantaneous flame front shapes are found via a solution for the usual Schwab–Zeldovitch parameter. Numerical results obtained from the analytical solution uncover the strong bearing that droplet grouping, induced by the standing wave flow-field, can have on flame height, shape, and type (over- or under-ventilated) and on the existence of multiple flame fronts.


Author(s):  
Johannes Gradl ◽  
Florian Schwertfirm ◽  
Hans-Christoph Schwarzer ◽  
Hans-Joachim Schmid ◽  
Michael Manhart ◽  
...  

Mixing and consequently fluid dynamic is a key parameter to tailor the particle size distribution (PSD) in nanoparticle precipitation. Due to fast and intensive mixing a static T-mixer configuration is capable for synthesizing continuously nanoparticles. The flow and concentration field of the applied mixer is investigated experimentally at different flow rates by Particle Image Velocimetry (PIV) and Laser Induced Fluorescence (LIF). Due to the PIV measurements the flow field in the mixer was characterized qualitatively and the mixing process itself is quantified by the subsequent LIF-measurements. A special feature of the LIF set up is to detect structures in the flow field, which are smaller than the Batchelor length. Thereby a detailed insight into the mixing process in a static T-Mixer is given. In this study a CFD-based approach using Direct Numerical Simulation (DNS) in combination with the solid formation kinetics solving population balance equations (PBE) is applied, using barium sulfate as modeling material. A Lagrangian Particle Tracking strategy is used to couple the flow field information with a micro mixing model and with the classical theory of nucleation. We found that the DNS-PBE approach including macro and micro mixing, combined with the population balance is capable of predicting the full PSD in nanoparticle precipitation for different operating parameters. Additionally to the resulting PSD, this approach delivers a 3D-information about all running subprocesses in the mixer, i.e. supersaturation built-up or nucleation, which is visualized for different process variables.


Author(s):  
Alireza Dastan ◽  
Omid Abouali

In this paper pressure drop and particle deposition in a microchannel with a hydraulic diameter of 225 micrometer is investigated numerically. Several hundred micron length fibers caught at the entrance of the channels making a “fiber web” also is modeled in this research. Governing equations for the flow field are solved with an Eulerian approach while the equations of particle motion in the flow are solved by a Lagrangian approach. Assuming the symmetry in the domain, one channel and the corresponding plenum are studied in the computational domain. For studying the effects of fibers in the flow, two fiber webs with four and six solid fibers are studied. The increase of pressure drop in the microchannel because of the entrance fiber web is computed and discussed. Also deposition and collection of the particles with various diameters at the fiber webs are also presented.


Sign in / Sign up

Export Citation Format

Share Document