scholarly journals Small-scale dynamics of settling, bidisperse particles in turbulence

2018 ◽  
Vol 839 ◽  
pp. 594-620 ◽  
Author(s):  
Rohit Dhariwal ◽  
Andrew D. Bragg

Mixing and collisions of inertial particles at the small scales of turbulence can be investigated by considering how pairs of particles move relative to each other. In real problems the two particles will have different sizes, i.e. they are bidisperse, and the effect of gravity on their motion is often important. However, how turbulence and gravity compete to control the motion of bidisperse inertial particles is poorly understood. Motivated by this, we use direct numerical simulations (DNS) to investigate the dynamics of settling, bidisperse particles in isotropic turbulence. In agreement with previous studies, we find that without gravity (i.e. $Fr=\infty$, where $Fr$ is the Froude number), bidispersity leads to an enhancement of the relative velocities, and a suppression of their spatial clustering. For $Fr<1$, the relative velocities in the direction of gravity are enhanced by the differential settling velocities of the bidisperse particles, as expected. However, we also find that gravity can strongly enhance the relative velocities in the ‘horizontal’ directions (i.e. in the plane normal to gravity). This non-trivial behaviour occurs because fast settling particles experience rapid fluctuations in the fluid velocity field along their trajectory, leading to enhanced particle accelerations and relative velocities. Indeed, the results show that even when $Fr\ll 1$, turbulence can still play an important role, not only on the horizontal motion, but also on the vertical motion of the particles. This is related to the fact that $Fr$ only characterizes the importance of gravity compared with some typical acceleration of the fluid, yet accelerations in turbulence are highly intermittent. As a consequence, there is a significant probability for particles to be in regions of the flow where the Froude number based on the local, instantaneous fluid acceleration is ${>}1$, even though the typically defined Froude number is $\ll 1$. This could imply, for example, that extreme events in the mixing of settling, bidisperse particles are only weakly affected by gravity even when $Fr\ll 1$. We also find that gravity drastically reduces the clustering of bidisperse particles. These results are strikingly different to the monodisperse case, for which recent results have shown that when $Fr<1$, gravity strongly suppresses the relative velocities in all directions, and can enhance clustering. Finally, we consider the implications of these results for the collision rates of settling, bidisperse particles in turbulence. We find that for $Fr=0.052$, the collision kernel is almost perfectly predicted by the collision kernel for bidisperse particles settling in quiescent flow, such that the effect of turbulence may be ignored. However, for $Fr=0.3$, turbulence plays an important role, and the collisions are only dominated by gravitational settling when the difference in the particle Stokes numbers is ${\geqslant}O(1)$.

2018 ◽  
Vol 860 ◽  
pp. 465-486 ◽  
Author(s):  
Nimish Pujara ◽  
Greg A. Voth ◽  
Evan A. Variano

We examine the dynamics of slender, rigid rods in direct numerical simulation of isotropic turbulence. The focus is on the statistics of three quantities and how they vary as rod length increases from the dissipation range to the inertial range. These quantities are (i) the steady-state rod alignment with respect to the perceived velocity gradients in the surrounding flow, (ii) the rate of rod reorientation (tumbling) and (iii) the rate at which the rod end points move apart (stretching). Under the approximations of slender-body theory, the rod inertia is neglected and rods are modelled as passive particles in the flow that do not affect the fluid velocity field. We find that the average rod alignment changes qualitatively as rod length increases from the dissipation range to the inertial range. While rods in the dissipation range align most strongly with fluid vorticity, rods in the inertial range align most strongly with the most extensional eigenvector of the perceived strain-rate tensor. For rods in the inertial range, we find that the variance of rod stretching and the variance of rod tumbling both scale as $l^{-4/3}$, where $l$ is the rod length. However, when rod dynamics are compared to two-point fluid velocity statistics (structure functions), we see non-monotonic behaviour in the variance of rod tumbling due to the influence of small-scale fluid motions. Additionally, we find that the skewness of rod stretching does not show scale invariance in the inertial range, in contrast to the skewness of longitudinal fluid velocity increments as predicted by Kolmogorov’s $4/5$ law. Finally, we examine the power-law scaling exponents of higher-order moments of rod tumbling and rod stretching for rods with lengths in the inertial range and find that they show anomalous scaling. We compare these scaling exponents to predictions from Kolmogorov’s refined similarity hypotheses.


2017 ◽  
Vol 822 ◽  
pp. 640-663 ◽  
Author(s):  
J. L. G. Oliveira ◽  
C. W. M. van der Geld ◽  
J. G. M. Kuerten

Three-dimensional particle tracking velocimetry is applied to particle-laden turbulent pipe flows at a Reynolds number of 10 300, based on the bulk velocity and the pipe diameter, for developed fluid flow and not fully developed flow of inertial particles, which favours assessment of the radial migration of the inertial particles. Inertial particles with Stokes number ranging from 0.35 to 1.11, based on the particle relaxation time and the radial-dependent Kolmogorov time scale, and a ratio of the root-mean-square fluid velocity to the terminal velocity of order 1 have been used. Core peaking of the concentration of inertial particles in up-flow and wall peaking in down-flow have been found. The difference in mean particle and Eulerian mean liquid velocity is found to decrease to approximately zero near the wall in both flow directions. Although the carrier fluid has all of the characteristics of the corresponding turbulent single-phase flow, the Reynolds stress of the inertial particles is different near the wall in up-flow. These findings are explained from the preferential location of the inertial particles with the aid of direct numerical simulations with the point-particle approach.


Author(s):  
Lionel Thomas ◽  
Benoiˆt Oesterle´

The dispersion of small inertial particles moving in a homogeneous, hypothetically stationary, shear flow is investigated using both theoretical analysis and numerical simulation, under one-way coupling approximation. In the theoretical approach, the previous studies are extended to the case of homogeneous shear flow with a corresponding anisotropic spectrum. As it is impossible to obtain a closed theoretical solution without some drastic simplifications, the motion of dispersed particles is also investigated using kinematic simulation where random Fourier modes are generated according to a prescribed anisotropic spectrum with a superimposed linear mean fluid velocity profile. The combined effects of particle Stokes number and dimensionless drift velocity (magnitude and direction) are investigated by computing the statistics from Lagrangian tracking of a large number of particles in many flow field realizations, and comparison is made between the observed effects in shear flow and in isotropic turbulence.


2014 ◽  
Vol 745 ◽  
pp. 279-299 ◽  
Author(s):  
Ryo Onishi ◽  
J. C. Vassilicos

AbstractThis study investigates the collision statistics of inertial particles in inverse-cascading two-dimensional (2D) homogeneous isotropic turbulence by means of a direct numerical simulation (DNS). A collision kernel model for particles with small Stokes number ($\mathit{St}$) in 2D flows is proposed based on the model of Saffman & Turner (J. Fluid Mech., vol. 1, 1956, pp. 16–30) (ST56 model). The DNS results agree with this 2D version of the ST56 model for $\mathit{St}\lesssim 0.1$. It is then confirmed that our DNS results satisfy the 2D version of the spherical formulation of the collision kernel. The fact that the flatness factor stays around 3 in our 2D flow confirms that the present 2D turbulent flow is nearly intermittency-free. Collision statistics for $\mathit{St}= 0.1$, 0.4 and 0.6, i.e. for $\mathit{St}<1$, are obtained from the present 2D DNS and compared with those obtained from the three-dimensional (3D) DNS of Onishi et al. (J. Comput. Phys., vol. 242, 2013, pp. 809–827). We have observed that the 3D radial distribution function at contact ($g(R)$, the so-called clustering effect) decreases for $\mathit{St}= 0.4$ and 0.6 with increasing Reynolds number, while the 2D $g(R)$ does not show a significant dependence on Reynolds number. This observation supports the view that the Reynolds-number dependence of $g(R)$ observed in three dimensions is due to internal intermittency of the 3D turbulence. We have further investigated the local $\mathit{St}$, which is a function of the local flow strain rates, and proposed a plausible mechanism that can explain the Reynolds-number dependence of $g(R)$. Meanwhile, 2D stochastic simulations based on the Smoluchowski equations for $\mathit{St}\ll 1$ show that the collision growth can be predicted by the 2D ST56 model and that rare but strong events do not play a significant role in such a small-$\mathit{St}$ particle system. However, the probability density function of local $\mathit{St}$ at the sites of colliding particle pairs supports the view that powerful rare events can be important for particle growth even in the absence of internal intermittency when $\mathit{St}$ is not much smaller than unity.


2019 ◽  
Vol 864 ◽  
pp. 925-970 ◽  
Author(s):  
Alec J. Petersen ◽  
Lucia Baker ◽  
Filippo Coletti

We study experimentally the spatial distribution, settling and interaction of sub-Kolmogorov inertial particles with homogeneous turbulence. Utilizing a zero-mean-flow air turbulence chamber, we drop size-selected solid particles and study their dynamics with particle imaging and tracking velocimetry at multiple resolutions. The carrier flow is simultaneously measured by particle image velocimetry of suspended tracers, allowing the characterization of the interplay between both the dispersed and continuous phases. The turbulence Reynolds number based on the Taylor microscale ranges from $Re_{\unicode[STIX]{x1D706}}\approx 200{-}500$, while the particle Stokes number based on the Kolmogorov scale varies between $St_{\unicode[STIX]{x1D702}}=O(1)$ and $O(10)$. Clustering is confirmed to be most intense for $St_{\unicode[STIX]{x1D702}}\approx 1$, but it extends over larger scales for heavier particles. Individual clusters form a hierarchy of self-similar, fractal-like objects, preferentially aligned with gravity and with sizes that can reach the integral scale of the turbulence. Remarkably, the settling velocity of $St_{\unicode[STIX]{x1D702}}\approx 1$ particles can be several times larger than the still-air terminal velocity, and the clusters can fall even faster. This is caused by downward fluid fluctuations preferentially sweeping the particles, and we propose that this mechanism is influenced by both large and small scales of the turbulence. The particle–fluid slip velocities show large variance, and both the instantaneous particle Reynolds number and drag coefficient can greatly differ from their nominal values. Finally, for sufficient loadings, the particles generally augment the small-scale fluid velocity fluctuations, which however may account for a limited fraction of the turbulent kinetic energy.


2016 ◽  
Vol 796 ◽  
pp. 659-711 ◽  
Author(s):  
Peter J. Ireland ◽  
Andrew D. Bragg ◽  
Lance R. Collins

In Part 1 of this study (Ireland et al., J. Fluid Mech., vol. 796, 2016, pp. 617–658), we analysed the motion of inertial particles in isotropic turbulence in the absence of gravity using direct numerical simulation (DNS). Here, in Part 2, we introduce gravity and study its effect on single-particle and particle-pair dynamics over a wide range of flow Reynolds numbers, Froude numbers and particle Stokes numbers. The overall goal of this study is to explore the mechanisms affecting particle collisions, and to thereby improve our understanding of droplet interactions in atmospheric clouds. We find that the dynamics of heavy particles falling under gravity can be artificially influenced by the finite domain size and the periodic boundary conditions, and we therefore perform our simulations on larger domains to reduce these effects. We first study single-particle statistics that influence the relative positions and velocities of inertial particles. We see that gravity causes particles to sample the flow more uniformly and reduces the time particles can spend interacting with the underlying turbulence. We also find that gravity tends to increase inertial particle accelerations, and we introduce a model to explain that effect. We then analyse the particle relative velocities and radial distribution functions (RDFs), which are generally seen to be independent of Reynolds number for low and moderate Kolmogorov-scale Stokes numbers $St$. We see that gravity causes particle relative velocities to decrease by reducing the degree of preferential sampling and the importance of path-history interactions, and that the relative velocities have higher scaling exponents with gravity. We observe that gravity has a non-trivial effect on clustering, acting to decrease clustering at low $St$ and to increase clustering at high $St$. By considering the effect of gravity on the clustering mechanisms described in the theory of Zaichik & Alipchenkov (New J. Phys., vol. 11, 2009, 103018), we provide an explanation for this non-trivial effect of gravity. We also show that when the effects of gravity are accounted for in the theory of Zaichik & Alipchenkov (2009), the results compare favourably with DNS. The relative velocities and RDFs exhibit considerable anisotropy at small separations, and this anisotropy is quantified using spherical harmonic functions. We use the relative velocities and the RDFs to compute the particle collision kernels, and find that the collision kernel remains as it was for the case without gravity, namely nearly independent of Reynolds number for low and moderate $St$. We conclude by discussing practical implications of the results for the cloud physics and turbulence communities and by suggesting possible avenues for future research.


2000 ◽  
Author(s):  
Daniel Huilier

Abstract A Lagrangian approach is developed to describe particle’s dispersion in a stationary, homogeneous and isotropic turbulent flow. Obviously, the particles’ dispersion is influenced by the fluid velocity fluctuations, which are classically simulated by a Monte Carlo process or Markov chains. However, some studies have shown the restrictions of these methods generating the fluid turbulent velocity and have suggested improvements to ensure that the Lagrangian model accounts for the three main effects governing the dispersion in gas-particle flows, namely the inertia, crossing trajectories and continuity effects. The first aim of this paper is to present an improved Lagrangian model which integrates the spatio-temporal characteristics of the fluid turbulence experienced by the particle. The agreement between the numerical results obtained and the analytical expressions derived by Wang and Stock (1993) will be very satisfying. Another interest is to investigate the role of the traditionally-neglected and troublesome added mass and history terms in numerical studies when long time dispersion of inertial particles is the primary concern. Indeed, we will observe that for a large range of values of the ratio of particle to fluid density, these non-stationary forces have statistically no influence on the characteristics of the turbulent particle dispersion and can be safely omitted.


2002 ◽  
Vol 02 (02) ◽  
pp. 295-310 ◽  
Author(s):  
H. SIGURGEIRSSON ◽  
A. M. STUART

The motion of an inertial particle in a Gaussian random field is studied. This is a model for the phenomenon of preferential concentration, whereby inertial particles in a turbulent flow can correlate significantly. Mathematically the motion is described by Newton's second law for a particle on a 2-D torus, with force proportional to the difference between a background fluid velocity and the particle velocity itself. The fluid velocity is defined through a linear stochastic PDE of Ornstein–Uhlenbeck type. The properties of the model are studied in terms of the covariance of the noise which drives the stochastic PDE. Sufficient conditions are found for almost sure existence and uniqueness of particle paths, and for a random dynamical system with a global random attractor. The random attractor is illustrated by means of a numerical experiment, and the relevance of the random attractor for the understanding of particle distributions is highlighted.


2019 ◽  
Vol 4 (10) ◽  
Author(s):  
Mohamad Ibrahim Cheikh ◽  
James Chen ◽  
Mingjun Wei

Sign in / Sign up

Export Citation Format

Share Document