Inertial Particles In Isotropic Turbulence Research Articles

Detailed characterization of extreme clustering at near-contact scales in isotropic turbulence

Recent measurements of inertial particles in isotropic turbulence (Hammond & Meng, J. Fluid Mech., vol. 921, 2021, A16) revealed surprising extreme clustering of particles at near-contact separations $(r)$ , whereby the radial distribution function, $g(r)$ , grows from $O(10)$ to $O(10^3)$ with a $(r/a)^{-6}$ scaling (where $a$ is the particle radius), and a surprising upturn of the mean inward particle-pair relative velocity (MIRV). Hydrodynamic interactions (HIs) were proposed to explain the extreme clustering, but despite predicting the correct scaling $(r/a)^{-6}$ , the HI theory underpredicted $g(r)$ by at least two orders of magnitude (Bragg et al., J. Fluid Mech., vol. 933, 2022, A31). To further understand the extreme clustering phenomenon and the relevance of HI, we characterize $g(r)$ and particle-pair kinematics for Stokes numbers $0.07 \leq St \leq 3.68$ in a homogeneous isotropic turbulence chamber using three-dimensional (3-D) particle tracking resolved to near–contact. A drift–diffusion equation governing $g(r)$ is presented to investigate the kinematic mechanisms of particle pairs. Measurements in all 24 conditions show that when $r/a\lessapprox 20$ , extreme clustering consistently occurs, scaling as $g(r) \sim (r/a)^{-k}$ with $4.5 \leq k \leq 7.6$ , which increases with $St$ . Here $g(r)$ varies with $St$ , particle size, density and polydispersity in ways that HI cannot explain. The extreme clustering region features an inward drift contributed by particle-pair turbophoresis and an inward radial relative acceleration. The latter indicates an interparticle attractive force at these separations that HI also cannot explain. The MIRV turns upward when approaching the extreme clustering region, opposite to direct numerical simulation predictions. These observations further support our previous assessment that extreme clustering arises from particle–particle interactions, but HI is not the main mechanism.

Hydrodynamic interactions and extreme particle clustering in turbulence

Expanding recent observations by Hammond & Meng (J. Fluid Mech., vol. 921, 2021, A16), we present a range of detailed experimental data of the radial distribution function (r.d.f.) of inertial particles in isotropic turbulence for different Stokes number,$St$, showing that the r.d.f. grows explosively with decreasing separationr, exhibiting$r^{-6}$scaling as the collision radius is approached, regardless of$St$or particle radius$a$. To understand such explosive clustering, we correct a number of errors in the theory by Yavuzet al.(Phys. Rev. Lett., vol. 120, 2018, 244504) based on hydrodynamic interactions between pairs of small, weakly inertial particles. A comparison between the corrected theory and the experiment shows that the theory by Yavuzet al.underpredicts the r.d.f. by orders of magnitude. To explain this discrepancy, we explore several alternative mechanisms for this discrepancy that were not included in the theory and show that none of them are likely the explanation. This suggests new, yet-to-be-identified physical mechanisms are at play, requiring further investigation and new theories.

Local analysis of the clustering, velocities, and accelerations of particles settling in turbulence

Using 3D Vorono\text{\"i} analysis, we explore the local dynamics of small, settling, inertial particles in isotropic turbulence using Direct Numerical Simulations (DNS). We independently vary the Taylor Reynolds number $R_\lambda \in[90,398]$, Froude number $Fr\equiv a_\eta/g\in[0.052,\infty]$ (where $a_\eta$ is the Kolmogorov acceleration, and $g$ is the acceleration due to gravity), and Kolmogorov scale Stokes number $St\equiv\tau_p/\tau_\eta\in[0,3]$. In agreement with previous results using global measures of particle clustering, such as the Radial Distribution Function (RDF), we find that for small Vorono\text{\"i} volumes (corresponding to the most clustered particles), the behavior is strongly dependent upon $St$ and $Fr$, but only weakly dependent upon $R_\lambda$, unless $St>1$. However, larger Vorono\text{\"i} volumes (void regions) exhibit a much stronger dependence on $R_\lambda$, even when $St\leq 1$, and we show that this, rather than the behavior at small volumes, is the cause of the sensitivity of the standard deviation of the Vorono\text{\"i} volumes that has been previously reported. We also show that the largest contribution to the particle settling velocities is associated with increasingly larger Vorono\text{\"i} volumes as the settling parameter $Sv\equiv St/Fr$ is increased. Our local analysis of the acceleration statistics of settling inertial particles shows that clustered particles experience a net acceleration in the direction of gravity, while particles in void regions experience the opposite. The particle acceleration variance, however, is a convex function of the Vorono\text{\"i} volumes, with or without gravity, which seems to indicate a non-trivial relationship between the Vorono\text{\"i} volumes and the sizes of the turbulent flow scales. Results for the variance of the fluid acceleration at the inertial particle "..."

Influence of Reynolds number on the motion of settling, bidisperse inertial particles in turbulence

Using Direct Numerical Simulations (DNS), we examine the effects of Taylor Reynolds number, $R_\lambda$, and Froude number, $Fr$, on the motion of settling, bidisperse inertial particles in isotropic turbulence. Particle accelerations play a key role in the relative motion of bidisperse particles, and we find that reducing $Fr$ leads to an enhancement of the accelerations, but a suppression of their intermittency. For Stokes numbers $St>1$, the effect of $R_\lambda$ on the accelerations is enhanced by gravity, since settling causes the particle accelerations to be affected by a larger range of flow scales. The results for the Probability Density Function (PDF) of the particle relative velocities show that for bidisperse particles, decreasing $Fr$ leads to an enhancement of their relative velocities in both the vertical (parallel to gravity) and horizontal directions. Importantly, our results show that even when the particles are settling very fast, turbulence continues to play a key role in their vertical relative velocities, and increasingly so as $R_\lambda$ is increased. This occurs because although the settling velocity may be much larger than typical velocities of the turbulence, due to intermittency, there are significant regions of the flow where the turbulence contribution to the particle motion is of the same order as that from gravitational settling...

Enhanced and suppressed multiscale dispersion of bidisperse inertial particles due to gravity

Using Direct Numerical Simulations (DNS), we investigate how gravity modifies the multiscale dispersion of bidisperse inertial particles in isotropic turbulence. The DNS has a Taylor Reynolds number $R_\lambda=398$, and we simulate Stokes numbers (based on the Kolmogorov timescale) in the range $St\leq 3$ , and consider Froude numbers $Fr = 0.052$ and $\infty$, corresponding to strong gravity and no gravity, respectively. The degree of bidispersity is quantified by the difference in the Stokes number of the particles $|\Delta St|$. We first consider the mean-square separation of bidisperse particle-pairs and find that without gravity (i.e. $Fr = \infty$), bidispersity leads to an enhancement of the the mean-square separation over a significant range of scales. When $|\Delta St|\geq O(1)$, the relative dispersion is further enhanced by gravity due to the large difference in the settling velocities of the two particles. However, when $|\Delta St|\ll1$, gravity suppresses the relative dispersion as the settling velocity contribution is small, and gravity suppresses the non-local contribution to the particle dynamics. In order to gain further insights, we consider separately the relative dispersion in the vertical (parallel to gravity) and horizontal directions. As expected, the vertical relative dispersion can be strongly enhanced by gravity due to differences in the settling velocities of the two particles. However, a key finding of our study is that gravity can also significantly enhance the horizontal relative dispersion. This non-trivial effect occurs because fast settling particles experience rapid fluctuations in the fluid velocity field along their trajectory, leading to enhanced particle accelerations and relative velocities. For sufficiently large initial particle separations...

Wavelet multiresolution analysis of particle-laden turbulence

The preferential concentration of inertial particles in isotropic turbulence is studied using wavelets, enabling the joint scale-space analysis of the flow field of the carrier phase, the concentration field of the dispersed phase, and interphase conditioned statistics.

Analysis of the forward and backward in time pair-separation probability density functions for inertial particles in isotropic turbulence

In this paper we investigate, using theory and direct numerical simulations (DNS), the forward in time (FIT) and backward in time (BIT) probability density functions (PDFs) of the separation of inertial particle pairs in isotropic turbulence. In agreement with our earlier study (Bragg et al., Phys. Fluids, vol. 28, 2016, 013305), where we compared the FIT and BIT mean-square separations, we find that inertial particles separate much faster BIT than FIT, with the strength of the irreversibility depending upon the final/initial separation of the particle pair and their Stokes number $St$. However, we also find that the irreversibility shows up in subtle ways in the behaviour of the full PDF that it does not in the mean-square separation. In the theory, we derive new predictions, including a prediction for the BIT/FIT PDF for $St\geqslant O(1)$, and for final/initial separations in the dissipation regime. The prediction shows how caustics in the particle relative velocities in the dissipation range affect the scaling of the pair-separation PDF, leading to a PDF with an algebraically decaying tail. The predicted functional behaviour of the PDFs is universal, in that it does not depend upon the level of intermittency in the underlying turbulence. We also analyse the pair-separation PDFs for fluid particles at short times, and construct theoretical predictions using the multifractal formalism to describe the fluid relative velocity distributions. The theoretical and numerical results both suggest that the extreme events in the inertial particle-pair dispersion at the small scales are dominated by their non-local interaction with the turbulent velocity field, rather than due to the strong dissipation range intermittency of the turbulence itself. In fact, our theoretical results predict that for final/initial separations in the dissipation range, when $St\gtrsim 1$, the tails of the pair-separation PDFs decay faster as the Taylor Reynolds number $Re_{\unicode[STIX]{x1D706}}$ is increased, the opposite of what would be expected for fluid particles.

Developments and difficulties in predicting the relative velocities of inertial particles at the small-scales of turbulence

In this paper, we consider the development of theoretical models to predict the relative velocities of inertial particles in isotropic turbulence. In particular, we use our recently developed theory for the backward-in-time (BIT) relative dispersion of inertial particles in turbulence [Bragg et al., “Forward and backward in time dispersion of fluid and inertial particles in isotropic turbulence,” Phys. Fluids 28, 013305 (2016)] to develop the theoretical model by Pan and Padoan [“Relative velocity of inertial particles in turbulent flows,” J. Fluid Mech. 661, 73 (2010)]. We focus on the most difficult regime to model, the dissipation range, and find that the modified Pan and Padoan model (that uses the BIT dispersion theory) can lead to significantly improved predictions for the relative velocities, when compared with the Direct Numerical Simulation (DNS) data. However, when the particle separation distance, r, is less than the Kolmogorov lengthscale, η, the modified model overpredicts the DNS data. We explain how these overpredictions arise from two assumptions in the BIT dispersion theory that are in general not satisfied when the final separation of the BIT dispersing particles is <η. We then demonstrate the failure of both the original and modified versions of the Pan and Padoan model to predict the correct scale-invariant forms for the inertial particle relative velocity structure functions in the dissipation regime. It is shown how this failure, which is also present in other models, is associated with our current inability to correctly predict not only the quantitative but also the qualitative behavior of the radial distribution function in the dissipation range when St=O(1).

The effect of Reynolds number on inertial particle dynamics in isotropic turbulence. Part 2. Simulations with gravitational effects

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.

Forward and backward in time dispersion of fluid and inertial particles in isotropic turbulence

In this paper, we investigate both theoretically and numerically the Forward-In-Time (FIT) and Backward-In-Time (BIT) dispersion of fluid and inertial particle-pairs in isotropic turbulence. Fluid particles are known to separate faster BIT than FIT in three-dimensional turbulence, and we find that inertial particles do the same. However, we find that the irreversibility in the inertial particle dispersion is in general much stronger than that for fluid particles. For example, the ratio of the BIT to FIT mean-square separation can be up to an order of magnitude larger for the inertial particles than for the fluid particles. We also find that for both the inertial and fluid particles, the irreversibility becomes stronger as the scale of their separation decreases. Regarding the physical mechanism for the irreversibility, we argue that whereas the irreversibility of fluid particle-pair dispersion can be understood in terms of a directional bias arising from the energy transfer process in turbulence, inertial particles experience an additional source of irreversibility arising from the non-local contribution to their velocity dynamics, a contribution that vanishes in the limit St → 0, where St is the particle Stokes number. For each given initial (final, in the BIT case) separation, r0, there is an optimum value of St for which the dispersion irreversibility is strongest, as such particles are optimally affected by both sources of irreversibility. We derive analytical expressions for the BIT, mean-square separation of inertial particles and compare the predictions with numerical data obtained from a Reλ ≈ 582 (where Reλ is the Taylor Reynolds number) Direct Numerical Simulation (DNS) of particle-laden isotropic turbulent flow. The small-time theory, which in the dissipation range is valid for times ≤max[Stτη, τη] (where τη is the Kolmogorov time scale), is in excellent agreement with the DNS. The theory for long-times is in good agreement with the DNS provided that St is small enough so that the inertial particle motion at long-times may be considered as a perturbation about the fluid particle motion, a condition that would in fact be satisfied for arbitrary St at sufficiently long-times in the limit Reλ → ∞.

Settling regimes of inertial particles in isotropic turbulence

AbstractWe investigate the settling speeds and root mean square (r.m.s.) velocities of inertial particles in isotropic turbulence with gravity using experiments with water droplets in air turbulence from 32 loudspeaker jets and direct numerical simulations (DNS). The dependence on particle inertia, gravity and the scales of both the smallest and largest turbulent eddies is investigated. We isolate the mechanisms of turbulence settling modification and find that the reduced settling speeds of large particles in experiments are due to nonlinear drag effects. We demonstrate using DNS that reduced settling speeds with linear drag (e.g. see Nielsen, J. Sedim. Petrol., vol. 63, 1993, pp. 835–838) only arise in artificial flows that, by design, eliminate preferential sweeping by the eddies. Gravity and inertia both reduce the particle r.m.s. velocities and falling particles are more responsive to vertical than to horizontal fluctuations. The model by Wang & Stock (J. Atmos. Sci., vol. 50, 1993, pp. 1897–1913) captures these trends.

New insights from comparing statistical theories for inertial particles in turbulence: II. Relative velocities

Part I of this two-part series compared two theories for the radial distribution function (RDF), a statistical measure of the clustering of inertial particles in isotropic turbulence. In Part II, we will contrast three theoretical models for the relative velocities of inertial particles in isotropic turbulence, one by Zaichik et al (2009 New. J. Phys. 11 103018), the second by Pan et al (2010 J. Fluid Mech. 661 73) and the third by Gustavsson et al (2011 Phys. Rev. E. 84 045304). We find that in general they describe the relative velocities in qualitatively similar ways, capturing the influence of the non-local dynamics on the formation of caustics and non-smooth scaling behavior in the dissipation range. We then compare the theoretical predictions with direct numerical simulation data and find that although they capture the qualitative behavior of the data consistently, they differ quantitatively, and we discuss the possible sources of error in each of the theories. Finally, we consider how the Zaichik et al theory predicts that the formation of caustics modifies the form of the RDF, and show that the theory describes the behavior of the RDF for particles whose response time scales with the inertial range time scales of the turbulence.

New insights from comparing statistical theories for inertial particles in turbulence: I. Spatial distribution of particles

In this paper, we contrast two theoretical models for the spatial clustering of inertial particles in isotropic turbulence, one by Chun et al (2005 J. Fluid Mech. 536 219) and the other by Zaichik et al (2007 Phys. Fluids 19 113308). Although their predictions for the radial distribution function are similar in the regime , they appear to describe the physical mechanism responsible for the clustering in quite different ways. We demonstrate why the theories generate such similar results in the regime by showing that the clustering mechanism in the Chun et al theory captures the leading order effects of the clustering mechanism in the Zaichik et al theory for . However, outside of this regime, the similarity between the predictions of the theories breaks down, and we consider the sources of the differences as well as the physical meaning and implications of the differences. Using DNS data we then show that the clustering mechanism described by the Zaichik et al theory accurately describes the clustering up to , and we identify a possible source of error for some of the slight quantitative discrepancies at larger St. We then compare these theories with others in the literature and attempt to reconcile as many of the physical explanations for clustering as we can. Finally, we consider the relationship between clustering in isotropic turbulence and the near-wall accumulation of inertial particles in a turbulent boundary layer, and how they scale with the Stokes number in the weak inertia limit.

Dispersion and temperature statistics of inertial particles in isotropic turbulence

The dispersion and temperature distribution of inertial particles are important in many turbulent, multiphase flow problems. In order to understand these better, direct numerical simulations (DNSs) are performed for inertial particles in a fluid with a constant temperature gradient and whose motion is either statistically stationary or decaying, isotropic turbulence. It is found that, for long times, the dispersion of inertial particles is the greatest when the Stokes number, Stη=τp/τη, is of order 1, where τp and τη are, respectively, the particle response time and the flow Kolmogorov time scale. A similar result is found for the long time behavior of the time rate of change of the mean-square particle temperature fluctuations, d⟨Tp′2⟩/dt. To understand the DNS results better, an evolution equation for ⟨Tp′2⟩, along with the short and long time limits, is derived analytically from the thermal energy equation for inertial particles.

Lagrangian measurements of inertial particle accelerations in a turbulent boundary layer

Two-dimensional Lagrangian acceleration statistics of inertial particles in a turbulent boundary layer with free-stream turbulence are determined by means of a particle tracking technique using a high-speed camera moving along the side of the wind tunnel at the mean flow speed. The boundary layer is formed above a flat plate placed horizontally in the tunnel, and water droplets are fed into the flow using two different methods: sprays placed downstream from an active grid, and tubes fed into the boundary layer from humidifiers. For the flow conditions studied, the sprays produce Stokes numbers varying from 0.47 to 1.2, and the humidifiers produce Stokes numbers varying from 0.035 to 0.25, where the low and high values refer to the outer boundary layer edge and the near-wall region, respectively. The Froude number is approximately 1.0 for the sprays and 0.25 for the humidifiers, with a small variation within the boundary layer. The free-stream turbulence is varied by operating the grid in the active mode as well as a passive mode (the latter behaves as a conventional grid). The boundary layer momentum-thickness Reynolds numbers are 840 and 725 for the active and passive grid respectively. At the outer edge of the boundary layer, where the shear is weak, the acceleration probability density functions are similar to those previously observed in isotropic turbulence for inertial particles. As the boundary layer plate is approached, the tails of the probability density functions narrow, become negatively skewed, and their peak occurs at negative accelerations (decelerations in the streamwise direction). The mean deceleration and its root mean square (r.m.s.) increase to large values close to the plate. These effects are more pronounced at higher Stokes number. In the vertical direction, there is a slight downward mean deceleration and its r.m.s., which is lower in magnitude than that of the streamwise component, peaks in the buffer region. Although there are free-stream turbulence effects, and the complex boundary layer structure plays an important role, a simple model suggests that the acceleration behaviour is dominated by shear, gravity and inertia. The results are contrasted with inertial particles in isotropic turbulence and with fluid particle acceleration statistics in a boundary layer. The background velocity field is documented by means of hot-wire anemometry and laser Doppler velocimetry measurements. These appear to be the first Lagrangian acceleration measurements of inertial particles in a shear flow.

Relative dispersion of two inertial particles in turbulent flow

A two-point statistical model based on the kinetic equation for the probability density function (PDF) of relative velocity is used to perform analysis of relative dispersion of inertial particles in isotropic turbulence. It is demonstrated that, in the case of highly inertial particles, the Richardson law for dispersion may be valid only in the case of very high values of the Reynolds number.

Collision rates of bidisperse inertial particles in isotropic turbulence

This paper presents two statistical models for predicting collision rates of bidisperse heavy particles suspended in homogeneous isotropic turbulence. One of the models is based on the assumption that the joint fluid-particle velocity distribution function is Gaussian. The other model stems from a kinetic equation for the two-point probability density function of the velocity distributions of two particles. The validity of these models is tested against DNS data by Zhou, Wexler, and Wang [J. Fluid Mech. 433, 77 (2001)]. Comparisons between the model predictions and DNS results demonstrate encouraging agreement.