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

Mohammadreza Momenifar and Andrew D. Bragg ∗ Department of Civil and Environmental Engineering,

Duke University, Durham, North Carolina 27708, USA

## Abstract

Using 3D Vorono¨ ı 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 λ ∈ [90 , 398], Froude number Fr ≡ a /g η ∈ [0 052 . , ∞ ] (where a η is the Kolmogorov acceleration, and g is the acceleration due to gravity), and Kolmogorov scale Stokes number St ≡ τ /τ p η ∈ [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¨ ı volumes (corresponding to the most clustered particles), the behavior is strongly dependent upon St and Fr , but only weakly dependent upon R λ , unless St &gt; 1. However, larger Vorono¨ ı volumes (void regions) exhibit a much stronger dependence on R λ , even when St ≤ 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¨ ı volumes that has been previously reported. We also show that the largest contribution to the particle settling velocities is associated with increasingly larger Vorono¨ ı volumes as the settling parameter Sv ≡ 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¨ ı volumes, with or without gravity, which seems to indicate a non-trivial relationship between the Vorono¨ ı volumes and the sizes of the turbulent flow scales. Results for the variance of the fluid acceleration at the inertial particle positions are of the order of the square of the Kolmogorov acceleration and depend only weakly on Vorono¨ ı volumes. These results call into question the 'sweep-stick' mechanism for particle clustering in turbulence which would lead one to expect that clustered particles reside in the special regions where the fluid acceleration is zero (or at least small).

We then consider the properties of particles in clusters, which are regions of connected Vorono¨ ı cells whose volume is less than a certain threshold. The results show self-similarity of the clusters, and that the statistics of the cluster volumes depends only weakly on St , with a stronger dependance on Fr and R λ . Finally, we compare the average settling velocities of all particles in the flow with those in clusters, and show that those in the clusters settle much faster, in agreement with previous work. However, we also find that this difference grows significantly with increasing R λ and exhibits a non-monotonic dependence on Fr . The kinetic energy of the particles, however, are almost the

same for particles whether they are in clusters or not.

∗ andrew.bragg@duke.edu

## I. INTRODUCTION

Turbulent multiphase flows have been the subject of a numerous studies for many years because of their broad range of applications in nature and industrial contexts such as distribution of greenhouse gases [1], cloud formation [2], volcanic eruptions [3], contaminant transport [4], heat exchangers [5], boiling heat transfer [6], combustion processes [7], air-lift pumps [8], cavitating flows [9], sedimentation process in settling tanks [10], thermosyphons [11] and fluidized bed reactors [12] to name a few. The present work explores the class of multi-phase flows involving the motion of small, settling inertial particles in turbulent flows.

A well-known and striking effect of particle inertia is that it leads to the spatial clustering of the particles in turbulent flows [4, 13]. This has important implications for the rate at which particles collide in turbulent flows [14, 15, 16, 17]. It is also important since regions with a high concentration of particles have the potential to significantly modify the underlying turbulent flow in those same regions, through momentum coupling with the fluid [18]. Several different methods and techniques have been used to analyze the particle clustering, such as box-counting ([19]), Vorono¨ ı analysis ([20]), Lyapunov exponents ([21]), Minkowski functionals ([22]) and radial distribution functions (RDFs) ([14, 17, 23]). These methods can provide different insights and perspectives into the spatial clustering, and two of the most commonly used methods are the RDF and Vorono¨ ı analysis.

The RDF, defined by the probability of finding a particle at a given distance from a reference particle, is unaffected by variation in the number of particles in the domain (except for statistical noise) and is directly related to the particle collision kernel ([14]). However, since it is constructed using spatial averages, the RDF cannot be used to consider the properties of individual clusters of particles, and only provides global information on the clustering. In Vorono¨ ı analysis, the domain is divided into cells associated with each individual particle, where any point in a given cell is closer to that particle than any other. This approach gives local insight into the particle clustering, with the local particle concentration field given by the inverse of the volume of the Vorono¨ ı cell. Using this technique enables characterization of the topology, kinematics and dynamics of individual clusters in the turbulent flow. However, a disadvantage is that unlike the RDF, the results of Vorono¨ ı analysis are sensitive to the number of particles in the flow, at least quantitatively ([24]). As discussed in [25], since the RDF and Vorono¨ ı analysis provide different information on the particle cluster-

ing (e.g. global vs. local information), they can both be used to provide a more complete understanding of the problem.

Recently, extensive investigations into the effects of the flow Taylor Reynolds number, R λ , and the Froude number Fr (quantifying the effect of gravity) on the particle clustering have been performed using the RDF computed from Direct Numerical Simulations (DNS) [17, 23, 26]. These studies have shown that when the Kolmogorov-scale Stokes number of the particle, St , is glyph[lessorsimilar] 1, there is virtually no effect of R λ on the clustering, whereas the clustering increases with increasing R λ when St &gt; 1. Concerning the effect of gravity, it was found that for St glyph[lessorsimilar] 1, gravity leads to a reduction in the clustering, while for St &gt; 1 it enhances the clustering. In contrast, it appears that such a systematic study of the parametric dependence of clustering as quantified by a Vorono¨ ı analysis has not been undertaken. We now summarize what has been done on this.

The experimental study of [27] considered the clustering of small water droplets in isotropic turbulence using Vorono¨ ı analysis. They examined the range 200 ≤ R λ ≤ 400 and 2 ≤ St ≤ 10, and reported that the level of clustering was enhanced with increasing R λ . The experimental study of [28] considered 170 ≤ R λ ≤ 450, 0 1 . ≤ St ≤ 5 and found that the clustering is strongly enhanced with increasing R λ , but with weak dependence on St . These studies focused on the regime of large Fr , where the effects of gravity on the particle motion are weak. The recent experimental study of [25] also used Vorono¨ ı analysis to quantify the clustering and explored the range 200 ≤ R λ ≤ 500, 0 37 . ≤ St ≤ 20 8, and . 0 73 . ≤ Fr ≤ 3 7. . They found that the clusters exhibit a range of sizes (up to the integral scale of the flow), posses self-similarity, and are elongated and aligned with the gravity direction [25]. They also found that the clustering is enhanced with increasing R λ , and in agreement with the RDF analysis in [17, 23, 26], that gravity can enhance the clustering of the inertial particles. A limitation of these experimental results is that their Vorono¨ ı analysis is only two-dimensional, and it is known that at least in the context of the RDF, two-dimensional analysis of particle clustering can differ significantly from the full threedimensional results [29]. Moreover, and perhaps more importantly, in these experiments it was not possible to independently vary R λ and Fr , and so a systematic understanding of their individual effects on the clustering as quantified by Vorono¨ ı analysis is lacking. DNS investigations into the effect of R λ on the Vorono¨ ı analysis of particle clustering have been undertaken, and in agreement with the experimental results, they show that the clustering

is enhanced with increasing R λ [24]. However, only the limited range 75 ≤ R λ ≤ 180 was considered, and they did not consider the effect of gravity.

In addition to analysis of the local particle concentration using Vorono¨ ı analysis, the DNS study of [30] explored the idea of coherent clusters of inertial particles, with Fr = 0 1 . in a low Reynolds number flow, R λ = 65. They defined a coherent cluster as a group of connected Vorono¨ ı cells whose volumes are smaller than a certain threshold. Moreover, they specified that the total volume of the cluster of connected cells must exceed 8 η 3 (where η is the Kolmogorov length scale) in order to be considered a coherent cluster. They found that the coherent clusters exhibit a tendency to align with the direction of the local vorticity vector, and that the cluster volume depends strongly on St . Furthermore, in the presence of gravity, coherent clusters were found to be elongated and aligned with the gravity direction, as also found experimentally in [25]. Since they only considered a single set of values for R λ and Fr , the effects of these parameters on the coherent clusters is not known.

In addition to using Vorono¨ ı analysis to quantify the particle clustering, has also been used to understand the settling rates of particles in turbulence and how this is related to the local particle concentration and turbulence properties [19, 31, 32, 33, 34]. Again, however, the effects of R λ and Fr on the results was not systematically explored. It has also been recently shown theoretically and numerically how the settling of particles in turbulence is affected by a range of flow scales that depends on St, R λ , and Fr [35]. It is therefore of interest to explore the effect of St, R λ , and Fr on the particle settle speeds using Vorono¨ ı analysis in order to gain further insights into the multiscale nature of the problem.

A recent striking finding concerning the effect of gravitational settling is that it can dramatically enhance the fluctuating accelerations of inertial particles in turbulence [23, 36], which in turn can profoundly affect collision and mixing rates of settling particles in turbulence [26, 37, 38]. An analysis of how these enhanced accelerations might be related to the local particle concentration has not yet been undertaken, but which would provide insights concerning the physical mechanism responsible for the effect, and also the scales at which it occurs.

In summary of this brief review, it is clear that there is a need to systematically explore the effects of R λ , and Fr on particle clustering and geometry of clusters as analyzed using Vorono¨ ı analysis, and to understand how this might differ from the perspective gained using the RDF analysis in [17, 23, 26]. Moreover, there is a need to explore how the en-

hanced fluctuating particle accelerations due to gravitational settling might be connected to the properties of the local particle concentration in order to gain further insights into the mechanisms responsible for this effect. In the present work we address these issues.

## II. COMPUTATIONAL DETAILS

Since one of the main goals of this work is to understand how the R , Fr λ dependence of particle clustering as determined by Vorono¨ ı analysis might differ from that based on an RDF analysis in [17, 23, 26], we therefore consider the same description for the particles and fluid as that in [17, 23, 26]. In particular, we focus on the motion of a dilute suspension of one-way coupled, small, heavy, spherical inertial particles whose motion is governed by a simplified version of the equation of Maxey &amp; Riley ([39])

$$\ddot { x } ^ { p } ( t ) \equiv \dot { v } ^ { p } ( t ) = \frac { u ( x ^ { p } ( t ), t ) - v ^ { p } ( t ) } { \tau _ { p } } + g,$$

where x p ( ) t , v p ( ) are the particle position and velocity vectors, t u x ( p ( ) t , t ) is the fluid velocity at the particle position, g is the gravitational acceleration, and τ p is the particle response time. Furthermore, fluid particles are tracked by solving ˙ x p ( ) t ≡ u x ( p ( ) t , t ).

| Parameter       | DNS1      | DNS2       | DNS3       | DNS4      |
|-----------------|-----------|------------|------------|-----------|
| N               | 128       | 128        | 1024       | 512       |
| R λ             | 93        | 94         | 90         | 224       |
| Fr              | ∞         | 0.3        | 0.052      | ∞         |
| L               | 2 π       | 2 π        | 16 π       | 2 π       |
| ν               | 0.005     | 0.005      | 0.005      | 0.0008289 |
| glyph[epsilon1] | 0.324     | 0.332      | 0.257      | 0.253     |
| l               | 1.48      | 1.49       | 1.47       | 1.40      |
| l/η             | 59.6      | 60.4       | 55.6       | 204       |
| u ′             | 0.984     | 0.996      | 0.912      | 0.915     |
| u ′ /u η        | 4.91      | 4.92       | 4.82       | 7.60      |
| T L             | 1.51      | 1.50       | 1.61       | 1.53      |
| T L /τ η        | 12.14     | 12.24      | 11.52      | 26.8      |
| κ max η         | 1.5       | 1.48       | 1.61       | 1.66      |
| N p             | 262,144   | 262,144    | 16,777,216 | 2,097,152 |
| Parameter       | DNS5      | DNS6       | DNS7       | DNS8      |
| N               | 512       | 1024       | 1024       | 1024      |
| R λ             | 237       | 230        | 398        | 398       |
| Fr              | 0.3       | 0.052      | ∞          | 0.052     |
| L               | 2 π       | 4 π        | 2 π        | 2 π       |
| ν               | 0.0008289 | 0.0008289  | 0.0003     | 0.0003    |
| glyph[epsilon1] | 0.2842    | 0.239      | 0.223      | 0.223     |
| l               | 1.43      | 1.49       | 1.45       | 1.45      |
| l/η             | 214       | 213        | 436        | 436       |
| u ′             | 0.966     | 0.914      | 0.915      | 0.915     |
| u ′ /u η        | 7.82      | 7.7        | 10.1       | 10.1      |
| T L             | 1.48      | 1.63       | 1.58       | 1.58      |
| T L /τ η        | 27.36     | 27.66      | 43.0       | 43.0      |
| κ max η         | 1.62      | 1.68       | 1.60       | 1.60      |
| N p             | 2,097,152 | 16,777,216 | 2,097,152  | 2,097,152 |

TABLE I: Simulation parameters for the DNS study of isotropic turbulence (arbitrary units). N is the number of grid points in each direction, R λ ≡ u λ/ν ′ is the Taylor micro-scale Reynolds number ( R λ ≡ √ 15 Re for homogeneous and isotropic flows), λ ≡ u / ′ 〈 ( ∇ u ) 2 〉 1 2 / is the Taylor micro-scale, L is the box size, ν is the fluid kinematic viscosity, glyph[epsilon1] ≡ 2 ν ∫ κ max 0 κ E κ 2 ( )d κ is the mean turbulent kinetic energy dissipation rate, l ≡ 3 π/ (2 k ) ∫ κ max 0 E κ /κ ( ) d κ is the integral length scale, η ≡ ν 3 4 / /glyph[epsilon1] 1 4 / is the Kolmogorov length scale, u ′ ≡ √ (2 k/ 3) is the fluid r.m.s. fluctuating velocity, k is the turbulent kinetic energy, u η is the Kolmogorov velocity scale, T L ≡ l/u ′ is the large-eddy turnover time, τ η ≡ √ ( ν/glyph[epsilon1] ) is the Kolmogorov time scale, κ max = √ 2 N/ 3 is the maximum resolved wavenumber, and N p is the number of particles per Stokes number.

The relevant non-dimensional parameters in this problem are the Stokes number ( St ≡ τ /τ p η ), the settling parameter ( Sv ≡ τ g/u p η , where the numerator represents particle's terminal velocity in the laminar flow) and the Froude number ( Fr ≡ a /g η = glyph[epsilon1] 3 4 / / ν ( 1 4 / g )) which characterize the particle's inertia, the effect of gravity on inertial particles and the effect of gravity on the flow, respectively. Here a η is the Kolmogorov acceleration, τ η is the Kolmogorov timescale, u η is the Kolmogorov velocity scale, and glyph[epsilon1] is the mean turbulent kinetic energy dissipation rate.

Direct Numerical Simulations (DNS) of the incompressible Navier-Stokes equation are performed on a triply periodic cube of length L , using a pseudo-spectral method on a uniform mesh with N 3 grid points. A deterministic forcing scheme is used to generate statistically stationary, isotropic turbulence. Details of numerical setup and simulations can be found in [17, 23, 40], and information on the simulations are provided in Table I.

As mentioned in the introduction, Vorono¨ ı analysis can be sensitive to the number of particles used in the analysis. Therefore, when comparing results from different DNS with different R λ and Fr we choose the number of particles such that the average distance between the particles is the same in each study, r ≡ L /N 1 3 / p = 7 9 . η . This number is comparable to that in the studies of [20, 24, 41]. A smaller value for r would be desirable, especially for performing the coherent cluster analysis [30]. However, the computational expense for the higher R λ cases we consider make it unfeasible to consider significantly smaller values of r .

Using our DNS data, we perform a 3D Vorono¨ ı analysis of the particles, and then compute a range of statistical quantities using the information from the analysis. In order to address the issue of open Vorono¨ ı cells at the boundaries of the domain, the particles across the boundaries are periodically repeated to ensure that all the cells are closed and the sum of Vorono¨ ı volumes is equal to the domain volume ([42, 43, 44]). Therefore, the Vorono¨ ı volumes associated with the particles adjacent to the boundaries are ill-defined and so are ignored in our analysis.

## III. RESULTS AND DISCUSSION

## A. Vorono¨ ı volume distributions

In figure 1, the Probability Distribution Function (PDF) of the Vorono¨ ı cell volumes (normalized by the mean Vorono¨ ı cells volume) for different St and with different Fr and R λ combinations are shown. The small/large Vorono¨ ı volumes represent the high/low concentration regions of the flow. For St = 0 the particles are randomly distributed, and their PDF follows that of a random Poisson process (RPP), described by [45]. Compared to the St = 0 case, the PDFs for St &gt; 0 are significantly higher at small volumes, indicating an enhanced probability of finding regions with high local particle concentration, i.e. clustering. In the experimental study of [46] they observed a crossover, such that at sufficiently small volumes, the PDFs of the Vorono¨ ı volumes for inertial particles drops below that for fluid particles. They conjectured that this effect was due to interactions between the particles when they are sufficiently close. In our results we do not observe this crossover, and neither did the numerical studies of [24, 30]. While this may be due to our neglect of the effect of particle interactions, we note that the crossover was also not observed in the recent experimental work of [25].

Concerning the effect of R λ , the results in figure 1 indicate that with respect to the behavior at small volumes, the effect of R λ depends on both St and Fr . When St &lt; 3, the effect of R λ is quite weak, sometimes leading to a slight increase of the PDF at small volumes that then saturates as R λ is increased (e.g. for the case St = 0 5, . Fr = 0 052), . while in other cases it leads to a slight decrease of the PDF at small volumes (e.g. for the case St = 0 7, . Fr = 0 3). . However, for St = 3, the sensitivity of the Vorono¨ ı volume PDF to R λ becomes much more apparent, with the clustering becoming stronger as R λ is increased. This dependence of the clustering on R λ is similar to that found in [17, 23] where the RDF was used to analyze the clustering. Arguments in [17, 23, 47] suggest that this behavior arises because unless St is sufficiently large, particles in the dissipation range are not able to remember their interaction with the inerital-range turbulence along their pathhistory. In this case, their motion is dominated by the dissipate range dynamics of the flow, and their dissipation range motion is not affected by the changing size of the inertial range as R λ is increased. Nevertheless, even if their dissipation range clustering is not affected

by the inertial range, the clustering could still be affected by the strong intermittency of the dissipation range turbulence [48]. In [17, 49] it was argued that because the RDF is a low-order measure of the particle phase-space dynamics, it is not very sensitive to the flow intermittency that as most apparent in the high-order moments. However, the PDF of the Vorono¨ ı volumes is not confined to low-order information, and could in principle be sensitive to intermittent fluctuations in the particle motion, yet our results imply that the particle clustering is only weakly affected by the flow intermittency. One possible explanation for this surprising behavior is that due to the preferential sampling of the flow by the inertial particles, the inertial particles avoid the intermittent regions of the flow [17, 23, 50], and as a result their spatial clustering is not affected by the increased intermittency in the underlying turbulence as R λ is increased.

Concerning the effect of Fr , the results in figure 1 indicate that with respect to the behavior at small volumes, the clustering initially becomes stronger as Fr is reduced from ∞ (no gravity case) to 0 3, but then becomes weaker as . Fr is further reduced. Therefore, gravity has a non-monotonic effect on the clustering, and whether it enhances or reduces the clustering depends on St and Fr .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-d4c56d8b-a5c7-4299-b69b-67f5b90b6499.jpg)

(a) St = 0

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-d2f3c47f-52cf-4d53-b0ba-aa1e8383362f.jpg)

(c) St = 0 5 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-f5e0b1ae-5083-42f9-8c68-252cdc67f406.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-46999664-0075-4333-999d-1727b92726f1.jpg)

(b) St = 0 2 .

(d) St = 0 7 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-1b4e7253-84c4-4230-842d-c360faa17872.jpg)

(e) St = 1

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-a1d4b569-3531-476a-94e6-deb01fd260c6.jpg)

(f) St = 3

FIG. 1: PDF of the Vorono¨ ı volumes (normalized by the mean volume) at different cases of Fr and R λ combinations for (a) St = 0, (b) St = 0 2, (c) . St = 0 5, (d) . St = 0 7,(e) . St = 1, and (f) St = 3. Different colors represents different cases. The dashed line represents the Random Poisson Distribution.

This effect of gravity on clustering agrees with the conclusions based on the RDF analysis in [23], although in that analysis the non-monotonic dependance of the clustering on Fr was not so clear since only the values Fr = ∞ , 0 052 were considered, and not the intermediate . value Fr = 0 3. . However, such non-monotonic dependence was observed using an RDF analysis in our previous study [26]. Physical arguments for the effect of gravity on clustering can be found in [23].

To obtain more insight into the R λ dependence, in figure 2 we plot the standard deviation (s.d.) of the Vorono¨ ı volumes as a function of St . In agreement with the experimental study of [46], we observe that for R λ = 90, the s.d. increases with increasing St until it reaches a peak around St ≈ 2, after which it decreases. Also in agreement with their results we find that the s.d. becomes larger and the peak shifts to higher values of St as R λ is increased. This is essentially due to the fact that as R λ is increased, more flow scales are present in the turbulence at which the particles can cluster. However, for any finite St , the growth of the s.d. with increasing R λ must saturate, since at sufficiently large scales the particles behave as fluid particles and do not cluster [51]. The results in figure 2 reveal a stronger sensitivity to R λ than is implied by the results in figure 1 which may indicate that the sensitivity of the s.d. to R λ is mainly in the behavior of the PDF of Vorono¨ ı cells with large volume, whose behavior is not clear in a loglog plot. To consider this, in figure 3 we plot the PDF of Vorono¨ ı volumes in a log-lin plot. Comparing the results with those in figure 1 confirms that it is the Vorono¨ ı cells with larger volumes that exhibit the strongest sensitivity to R λ and that the sensitivity to R λ is generally stronger in the presence of gravity.

The data in figure 1 shows that the PDF for inertial particles intersects that for fluid particles at two points, denoted by ϑ C and ϑ v , respectively. These points are often used to define thresholds for detecting clusters and voids, with volumes less than ϑ C denoting clustered regions, and volumes greater than ϑ v denoting void regions. In view of these definitions, it is then apparent that the R λ dependence of the s.d. of the Vorono¨ ı volumes comes mainly from the R λ of the void regions, and not the clustered regions. In this way our conclusions are consistent with the findings in [17, 23, 26] based on the RDF that the small-scale clustering is quite insensitive to R λ when St glyph[lessorsimilar] 1, and the strong R λ dependence of the s.d. of the Vorono¨ ı volumes does not contradict this.

A related point is that the results in figure 2 indicate that gravity suppresses the s.d. of the Vorono¨ ı volumes, while for the same cases (e.g. St = 3 , R λ = 398), figure 1 indicates

FIG. 2: Standard deviation of Vorono¨ ı cells normalized by their mean as a function of St .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-6412c33f-3433-49e0-9022-583d9f321fe9.jpg)

that gravity enhances the small-scale clustering. This is again because the s.d. is dominated by the behavior of the void regions rather than the clustered regions, and the void regions tend to be suppressed by gravity. This shows again the non-trivial effect of gravity on the particle clustering; while the probability of clustered and void regions both increase with increasing St , for a given St , decreasing Fr can enhance the probability of clustered regions while lowering the probability of void regions. These points illustrate that using the s.d. Vorono¨ ı volumes as a measure of clustering should be used with caution, since its properties are dominated by void, rather than clustered regions of the flow.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-6bd15c86-8979-4628-bc2e-63ee96146da7.jpg)

(a) St = 0

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-a3ece822-f142-45c9-b8e3-a3bf4b40abbb.jpg)

(c) St = 0 5 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-48b8f18a-2bd4-46bf-84b2-79d0d5f18c8e.jpg)

(e)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-ff37c946-1713-4060-ab5a-e674552509a4.jpg)

(b) St = 0 2 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-84c3a473-6445-460e-ae14-92fc6e5c8d15.jpg)

(d) St = 0 7 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-3bb7879e-c939-4340-86d4-36a6168d8e4d.jpg)

St

= 1

FIG. 3: PDF of the Vorono¨ ı volumes (normalized by the mean volume) at different cases of Fr and R λ combinations for (a) St = 0, (b) St = 0 2, (c) . St = 0 5, (d) . St = 0 7,(e) . St = 1, and (f) St = 3. Different colors represents different cases. The dashed line represents the Random Poisson Distribution.

In a number of previous studies it has been reported that the PDF of Vorono¨ ı volumes follows a log-normal distribution [20, 25, 27, 28, 31]. In figure 4 the centered and normalized PDFs of the logarithm of the Vorono¨ ı volumes are presented. The results show that for St ≥ 0 5 a Gaussian distribution can approximate the shape of . log ϑ ( ) in the interval ± 2 σ log ϑ ( ) , while large deviations from the Gaussian PDF are found outside of this range. In particular, the tails of the PDF are heavier than a Gaussian distribution, as also observed in [31]. However, the PDF seems to approach a Gaussian shape as St increases. For smaller St the Gaussian PDF is not followed, which is to be expected since in the limit St → 0 the PDF follows that for a RPP, a shown earlier. Therefore, the log-normal shape for the PDF of Vorono¨ ı volumes is only a reasonable approximation in certain regimes. Departures from the log-normal behavior exhibit different (and not entirely clear) dependencies on R λ and Fr for the left and right tails of the PDF, an explanation for which is not clear.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-f451d5f0-0d8e-4cd6-88cc-1e5da2692ce2.jpg)

(a) St = 0

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-163fec34-f619-4393-be9c-976c0fc9be3a.jpg)

(c) St = 0 5 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-bc7c3cbe-5f7e-4067-b381-08be2e748d0e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-2c2632af-379a-44ad-b99f-b6fd31bea802.jpg)

(b) St = 0 2 .

(d) St = 0 7 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-abdf969e-7541-4745-805d-73cfee624995.jpg)

(e) St = 1

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-fa16618f-0c55-4352-8d75-9ef28fbf4b74.jpg)

(f) St = 3

FIG. 4: PDF of the log-normal distribution of Vorono¨ ı volumes (centered at the mean and normalized by standard deviation) at different cases of Fr and R λ combinations for (a) St = 0, (b) St = 0 2, (c) . St = 0 5, (d) . St = 0 7,(e) . St = 1, and (f) St = 3. Different colors represents different cases and dashed line denotes the Gaussian distribution.

## B. Particle settling velocities

We now turn to consider the average velocity of particles conditioned on the Vorono¨ ı volumes to explore how the particle settling velocities depend upon the local particle concentration. The results in figure 5 clearly illustrate the correlation between the particle settling velocity and the local concentration (inverse of the Vorono¨ ı volumes), as also observed over smaller ranges of R λ and Fr in the DNS studies of [30, 32, 52] as well as the experimental study of [19]. The recent experiments of [25] also considered this quantity under a wide range of particle loadings, to examine the influence of two-way coupling between the particles and fluid on the particle settling velocities and its relation to the local concentration. Their results at low particle loadings are similar to ours, suggesting, as expected, that the one-way coupled assumption yields physically realistic results in that regime.

The results show that the Fr for which the settling velocity reaches the largest value depends upon St , which is simply because in the particle equation of motion, the important non-dimensional number is the settling number Sv ≡ St/Fr . However, for St &gt; 0 2, . the settling enhancement is strongest for intermediate values of Fr , in agreement with the DNS study of [31]. For Fr = 0 3, the average settling velocity enhances almost monotonically with . increasing particle concentration (the decrease observed at small volumes for St = 0 2 may . be due to noise). However, for Fr = 0 052, there is a local concentration at which the mean . settling velocity is maximum, and then decreases for higher concentrations. If we assume that larger Vorono¨ ı volumes are associated with larger flow scales, then this is consistent with the arguments in [35] that as Fr is decreased, the scales of the flow responsible for enhancing the particle settling velocity increase. The results also show that increasing R λ can significantly increase the settling velocities, especially for St ≥ 0 5. . This is explained by the recent theoretical analysis of [35] that shows that as R λ is increased, more and more energetic flow scales are introduced to the flow that are able to contribute to the particle settling velocity enhancement through the preferential sweeping mechanism described below.

Concerning the physical mechanism underlying the correlation between the settling velocity and concentration/clustering, since our DNS are one-way coupled the effect cannot be explained in terms of the particles modifying the local flow field, but in terms of the way the particles passively sample the underlying flow field. [4] provided a theoretical analysis and on the basis of this proposed that in the regime St glyph[lessmuch] 1, the particles are preferentially swept

around the downward (direction of gravity) moving side of vortices in the flow, and that as a result of this preferential sweeping, their average settling velocity exceeds the Stokes settling velocity. This is the so-called 'preferential sweeping mechanism' [53]. This theoretical analysis was recently extended by [35] to arbitrary St , wherein the basic mechanism is the same as that proposed by Maxey, except that now the size of the vortices around which the particles are preferentially swept is shown to depend essentially on St, Fr and R λ , and is a multiscale preferential sweeping mechanism. However, it was pointed out in [35] that unlike the regime St glyph[lessmuch] 1, in the regime St ≥ O (1) the mechanism generating the clustering is distinct from the preferential sweeping mechanism. Indeed, for St ≥ O (1) the clustering is generated by a non-local mechanism and not the preferential sampling of the local flow field [47, 51]. As a result, while the results in figure 5 for St glyph[lessmuch] 1 can be directly explained in terms of the preferential sweeping mechanism, the same does not apply for St ≥ O (1). A direct and unambiguous test of the preferential sweeping mechanism would be to condition the average particle velocity not on the concentration, but on some measure of the local preferential sampling of the flow.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-31fe4e2b-e693-4afc-8ed6-52fe41a74e19.jpg)

(a) St = 0 2 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-cca05347-a935-4d37-a4f1-524e11f50ada.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-e1bae086-44ea-48e4-a3b2-9ad65ea0ea4d.jpg)

(b)

St

= 0 5

.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-a4271a2e-a19e-4d6f-992a-c5f5defb25d1.jpg)

(c)

St

= 0 7

.

(d)

St

= 1

FIG. 5: Average velocity of particles conditioned on the PDF of the Vorono¨ ı volumes at different cases of Fr and R λ combinations for (a) St = 0, (b) St = 0 2, (c) . St = 0 5, (d) . St = 0 7, and (e) . St = 1. Different colors represents different cases.

## C. Local Analysis of Particle Accelerations

In [23, 36] it was shown that settling can significantly enhance the accelerations of inertial particles in turbulent flows, and in [26] it was also shown that the variance of the particle accelerations generally increase with increasing R λ and/or with decreasing Fr . To gain further insight into the behavior of the particle accelerations, we perform a local analysis by conditioning the particle acceleration behavior on the local Vorono¨ ı volume, to understand how the accelerations are linked to the concentration field.

Figures 6 and 7 show results for the particle acceleration ( a p ( ) t ≡ ˙ v p ( )) conditioned on t the local Vorono¨ ı volume, in the direction parallel and perpendicular to gravity, respectively. According to these results, in the absence of gravity ( Fr = ∞ ), the particle accelerations are zero (to within statistical noise) in both directions, independent of St and R λ . In the presence of gravity, this same behavior is observed for the accelerations in the direction normal to gravity, however, in the direction of gravity the quantity becomes finite for finite St . In particular, 〈 a p x 〉 ϑ becomes negative at small volumes and positive for large volumes. Moreover, 〈 a p x 〉 ϑ becomes increasingly negative as R λ is increased. Some of this behavior may be understood by considering the following identity

$$\langle a _ { x } ^ { p } \rangle \equiv \int _ { 0 } ^ { \infty } \langle a _ { x } ^ { p } \rangle _ { \vartheta } P ( \vartheta ) \, d \vartheta.$$

Since P ϑ ( ) ≥ 0, and 〈 a p x 〉 = 0 for statistically stationary, homogeneous turbulence, then the only way for (2) to be satisfied is if either 〈 a p x 〉 ϑ = 0 ∀ ϑ , or else its sign varies with ϑ . The latter situation is what we observe in the presence of gravity, with particles in high concentration regions (small ϑ ) exhibiting a net acceleration in the direction of gravity even though the net acceleration of all the particles is zero, i.e. 〈 a p x 〉 = 0. However, we are not able to explain why 〈 a p x 〉 ϑ is negative, rather than positive at small ϑ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-93db96ac-19c8-4d1c-990b-1be94626aabc.jpg)

(a) St = 0

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-c51d70d2-b76c-4fc0-98f3-03306177326a.jpg)

(c) St = 0 5 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-0bfbe16d-10c4-4bd6-9abe-d66f3b427e11.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-02efcc5c-c745-4f38-82d1-87e7c533410f.jpg)

(b) St = 0 2 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-d3d13341-6b59-4f64-8f1f-e1cab0c083a7.jpg)

(d) St = 0 7 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-12273f69-99c7-4eaa-8d9b-8a6488d7aeb2.jpg)

(e)

St

= 1

(f) St = 3

FIG. 6: Average acceleration of particles, in the gravity direction, conditioned on the PDF of the Vorono¨ ı volumes at different cases of Fr and R λ combinations for (a) St = 0, (b) St = 0 2, (c) . St = 0 5, (d) . St = 0 7,(e) . St = 1, and (f) St = 3. Different colors represents different cases.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-a8556e62-8b8d-42ab-a993-70f51b6c8883.jpg)

(a) St = 0

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-0530ac22-8a8b-405a-a43c-73ee226ae3b1.jpg)

(c) St = 0 5 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-2a7f7e66-8f1a-4e6c-acf7-fe8eaf524dfe.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-b07800ab-9f5d-480b-97a6-d6a0018d7ebe.jpg)

(b) St = 0 2 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-f0698a34-5316-4d34-8499-bcef96e4409d.jpg)

(d) St = 0 7 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-00676fde-a338-4358-938a-1b0528125b66.jpg)

(e)

St

= 1

(f)

St

= 3

FIG. 7: Average acceleration of particles, in the plane normal to the gravity direction, conditioned on the PDF of the Vorono¨ ı volumes at different cases of Fr and R λ combinations for (a) St = 0, (b) St = 0 2, (c) . St = 0 5, (d) . St = 0 7,(e) . St = 1, and (f) St = 3. Different colors represents different cases.

Figures 8 and 9 show results for the second moment of the particle accelerations conditioned on the Vorono¨ ı volumes, in the direction parallel and perpendicular to gravity, respectively. The results reveal a surprising nonlinear dependence on ϑ , that is present irrespective of Fr, R λ , that disappears for sufficiently small or large St . If one adopts a standard view that the Vorono¨ ı volumes are connected to the different range of scales in a turbulent flow, with small volumes associated with small turbulent flow scales etc, then this would explain why the acceleration variances increase with decreasing ϑ , since accelerations are largest at the smallest scales of a turbulent flow. However, this explanation cannot account for why at sufficiently large ϑ , the acceleration variances increase with increasing ϑ . The results in figures 8 and 9 therefore perhaps indicate that there is not a simple relationship between the size of the Vorono¨ ı volumes and the size of the turbulent flow scales. Furthermore, the strong dependence of this behavior on St indicates that the nonlinear dependence of the acceleration variances on ϑ is related to the spatial clustering of the particles.

This nonlinear dependence of the acceleration variances on ϑ is not qualitatively affected by Fr, R λ , however, the results show that decreasing Fr and/or increasing R λ lead to an enhancement of the particle accelerations at all ϑ (to within statistical noise). This would seem to imply, in agreement with the physical arguments in [23], that the enhancement of the inertial particle accelerations due to gravitational settling are not fundamentally associated with a change in how gravity causes the particles to cluster. Rather, it is caused by the fact that gravity causes the particles to experience rapid changes in the fluid velocity along their trajectory, which in turn leads to large particle accelerations. This effect, described more fully in [37], would operate even if the particles were uniformly distributed throughout the flow. Indeed, our results in figures 8 and 9 for St = 3 reveal strong enhancements of the particle accelerations due to gravity, even though the acceleration variances are almost independent of ϑ .

In figure 10 we show results for the variances of the fluid acceleration ( a F ) at the inertial particle positions, conditioned on ϑ , and for the direction of gravity. These results also exhibit, in general, a convex functional shape with respect to ϑ . They also show that as Fr is decreased, the accelerations become stronger, even for cases where the St = 0 are independent of Fr (in general the fluid particle accelerations should be independent of Fr , however since the different Fr cases correspond to different simulations, the results can depend on Fr , for example due to different domain sizes in the simulations). Since

the fluid acceleration measured at the inertial particle positions can only differ from the fluid acceleration at fluid particle positions due to preferential sampling of the flow by the particles, this indicates that the settling causes the inertial particles to sample regions of the flow exhibiting strong accelerations.

The results in figure 10 also present a challenge to the 'sweep-stick' mechanism that has been proposed to explain particle clustering in turbulent flows [54]. The sweep-stick mechanism states that inertial particles stick to points in the flow where the fluid acceleration is zero, and since these acceleration stagnation points are themselves clustered, this therefore explains the clustering of the inertial particles. In [51] a number of theoretical arguments were given that call into question the validity of this proposed mechanism for particle clustering, supported by results from DNS. If the sweep-stick mechanism were correct, we should expect to find that in regions where the particles are clustered (i.e. where ϑ is small), the fluid acceleration at the particle position should be zero, or at least very small compared with a η . The results in figure 10 show, however, that in the regions where the particles are clustered, the fluid accelerations at the inertial particle positions are large, i.e. their magnitudes are ≥ O ( a η ). Moreover, over a significant range of ϑ , the accelerations are almost independent of ϑ . These results therefore strongly call into question the validity of the sweep-stick mechanism. We refer the reader to [51] for a detailed discussion concerning theoretical issues with the sweep-stick mechanism.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-af5d8861-0a72-4e6c-bac3-2ee7c4a67d1f.jpg)

(a) St = 0

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-937837da-fe8f-499a-a625-6da464667d50.jpg)

(c) St = 0 5 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-4f89c4e8-64e3-4bb7-82d4-ee1be722d822.jpg)

(e)

St

= 1

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-351fdf84-2e93-4742-a2dd-44dadb106c75.jpg)

(b) St = 0 2 .

(d) St = 0 7 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-8fb13d0c-b361-46a0-856c-8b4679776c3c.jpg)

FIG. 8: Second moment of particle accelerations, in the gravity direction, conditioned on the PDF of the Vorono¨ ı volumes at different cases of Fr and R λ combinations for (a) St = 0, (b) St = 0 2, (c) . St = 0 5, (d) . St = 0 7,(e) . St = 1, and (f) St = 3. Different colors represents different cases.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-e56361aa-d535-40b3-bdd3-0b05e6112012.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-2fad2d46-52ff-4616-b38d-88020aa5d49d.jpg)

(a) St = 0

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-9f08cde5-2d76-4d85-b2c7-22e4241f1df1.jpg)

(c) St = 0 5 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-05019707-fb02-4fb7-ab5c-725f75b6165a.jpg)

(e)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-47d11834-3706-464b-b40e-9416fff3c7fe.jpg)

(b) St = 0 2 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-34ea2134-7e8b-4bfc-b487-e232fcf62221.jpg)

(d) St = 0 7 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-974db9c0-654c-49c7-94b6-1eb308d4874c.jpg)

St

= 1

FIG. 9: Second moment of particle accelerations, in the plane normal to the gravity direction, conditioned on the PDF of the Vorono¨ ı volumes at different cases of Fr and R λ combinations for (a) St = 0, (b) St = 0 2, (c) . St = 0 5, (d) . St = 0 7,(e) . St = 1, and (f) St = 3. Different colors represents different cases.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-ff4debec-3e21-44a4-9ed2-5635c4b1d822.jpg)

(a) St = 0

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-4630a225-39c2-47b8-8094-10c59a53c40b.jpg)

(c) St = 0 5 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-142b41d4-08d4-4566-b4b2-b08f82e70f00.jpg)

(e) St = 1

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-ec1c65a0-4d36-4dbb-8661-6aa63a8ee498.jpg)

(b) St = 0 2 .

(d) St = 0 7 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-a6823ea5-115f-4b57-acd9-81462094a4ca.jpg)

FIG. 10: Second moment of fluid accelerations at the particle position, in the gravity direction, conditioned on the PDF of the Vorono¨ ı volumes at different cases of Fr and R λ combinations for (a) St = 0, (b) St = 0 2, (c) . St = 0 5, (d) . St = 0 7,(e) . St = 1, and (f) St = 3. Different colors represents different cases.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-f5084229-cf85-4e3a-b069-2814c160aa81.jpg)

We now turn our attention to the results for clusters of particles, which are defined as connected Vorono¨ ı cells that have volumes ≤ ϑ C [30]. Figure 11 shows the PDF of cluster volumes normalized by the Kolmogorov scale. The results show that cluster sizes are distributed over a wide range of scales, from the dissipation up to integral range scales. In agreement with previous studies ([20, 27, 30, 32, 52]), we find that the right hand side of the PDF is well described as a power law with exponent -2, which implies self-similarity of the clusters, in line with results in [20, 55]. This slope was found to be ∼ -5 / 3 in the experimental study of [28].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-6bb98801-b516-4f61-8ff9-9e211af224fd.jpg)

(a) St = 0 4 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-58187b37-7371-4e97-9fe5-6bfdf14e494c.jpg)

(b) St = 0 7 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-179087b0-de89-4d5e-bbc1-c56f273fa0f1.jpg)

(c) St = 1

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-7e842138-dbc5-450a-8e8e-f13466096a93.jpg)

(d) St = 1 5 .

(e) St = 3

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-617d7c4a-fc1e-4824-b0e0-259040d2d5be.jpg)

FIG. 11: PDF of the coherent cluster volumes (normalized by the Kolmogorov scale) at different cases of Fr and R λ combinations for (a) St = 0 4, (b) . St = 0 7, (c) . St = 1, (d) St = 1 5, and (e) . St = 3. Different colors represents different cases. The dashed line represents a power law with exponent -2.

The results in Figure 11 reveal a much weaker dependence on St, R λ and Fr than the PDF of the Vorono¨ ı volumes considered earlier. Figure 12 shows the same results in a log-lin plot to emphasize the behavior for larger clusters, and again we find a weak dependence on St, R λ and Fr . One possible reason for the weak dependence is related to the fact that given that the average inter-particle distance is 7 9 . η , the clusters are quite large and so the effects of inertia are weaker than they would be for smaller clusters. Indeed, the study of [30] had a much smaller average inter-particle distance of ≈ 2 2 . η , and in their results the strongest effects of the particle inertia were found for the smallest clusters, as expected.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-8047a73c-935c-48b1-90a4-3efed2ba3494.jpg)

(a) St = 0 4 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-a021c8df-a995-4978-8fe7-6869d83d2378.jpg)

(c)

St

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-4d033e13-56ce-4683-9ce5-644928d1c707.jpg)

(b) St = 0 7 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-84383a95-f74e-45ee-a93a-3731cf38c41e.jpg)

= 1

(d) St = 1 5 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-31627e95-a6a5-4af1-8ca3-49c6f7f5ec75.jpg)

FIG. 12: PDF of the coherent cluster volumes (normalized by the Kolmogorov scale) at different cases of Fr and R λ combinations for (a) St = 0 4, (b) . St = 0 7, (c) . St = 1, (d) St = 1 5, and (e) . St = 3. Different colors represents different cases. The dashed line represents a power law with exponent -2.

To obtain clearer quantitative insight, the standard deviation (s.d.) of the coherent cluster volumes is shown in the figure 13. The results reveal a quite weak dependence on St , which was also reported in the experimental study of [28]. This is in striking contrast to the results for the s.d. of the Vorono¨ ı volumes shown in Figure 2, for which a very strong dependence on St was observed. However, the effect of Fr and R λ on the s.d. of the coherent cluster volumes is much stronger. In particular, the s.d. increases strongly with increasing R λ , while it increases only slightly in going from Fr = ∞ to 0.3, and then increases substantially in going from Fr = 0 3 to 0.052. .

FIG. 13: Standard deviation of Vorono¨ ı cells of clusters as a function of St .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-4244732d-c78c-4763-a72b-442e924bf891.jpg)

In figure 14, the average percentage of particles within clusters compared with the number in the whole domain is shown. In agreement with [30], we observe that this percentage increases with increasing St , though for R λ this reaches a maximum at St ≈ 2 and then reduces. We expect that a maximum would also be observed for the higher R λ , but at larger St (a maximum must occur since in the limits St → 0 and St → ∞ the percentage must go to zero). The results show that generally, increasing R λ leads to an increase of the percentage, while the dependence on Fr is non-monotonic, with the percentage generally being largest for the intermediate value Fr = 0 3. . It is important to note, however, that the percentages we observe are much larger than those observed by [30]. This difference is likely due to the fact that in their simulations the average inter-particle distance was 2 2 . η , while in ours it has the significantly larger value 7 9 . η . If this indeed the explanation, then it shows that the significance of the particle clusters, in terms of the number of particles contained within them as a fraction of the total number in the flow, depends essentially on the particle loading in the flow.

FIG. 14: Average percentage of total particles in clusters to the total number of particles in the whole domain, as a function of St and for different cases

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-384123f3-bc36-4969-bbc5-215d6ac66795.jpg)

Figure 15 shows results for the average cluster size as a function of the number of particles in the clusters, which provides information on the relationship between the coherent cluster sizes and the particle concentrations within them. (For brevity we only show the results for St = 1, but the behavior is the same for the other St cases). This relationship is found to be approximately linear, and is independent of all the control parameters over the range we have explored. Such a linear relationship, and its lack of sensitivity to St and Fr , were already reported in [30] for R λ = 65, however, our results confirm that this quantity is also independent of R λ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-dd476667-d693-4e58-872e-4fb702158691.jpg)

(a) Re = 90 &amp; Fr = ∞

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-ee03b07f-f69e-406d-8a44-fdc43468a9d5.jpg)

(c) Re = 224 &amp; Fr = ∞

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-bbc3d292-4a3c-4e6d-9efd-baa6b0cd0306.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-9474e65d-606c-470a-87ef-4a8fd8888234.jpg)

(b) Re = 90 &amp; Fr = 0 052 .

(d) Re = 224 &amp; Fr = 0 052 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-ae989a25-61ee-4c19-a7ea-7a0ff4919521.jpg)

(e) Re = 398 &amp; Fr = ∞

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-81cd58a7-7911-4f3a-a1e5-a97b8b58a68f.jpg)

(f) Re = 398 &amp; Fr = 0 052 .

FIG. 15: The cluster volume normalized by the Kolmogorov scale versus particle counts in clusters normalized by the whole particles in the domain, for St = 1.

FIG. 16: Difference between the average velocity of particles within the clusters and the whole particles in the domain.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-241a1566-1332-48a8-af0c-94b87d5ea5a4.jpg)

Finally, we turn to consider the velocity statistics of the particles in coherent clusters as compared with those in the entire domain. In figure 16 we show results for the difference in the average velocity of particles in clusters to that based on all particles in the flow. In agreement with those of [30], our results show that the mean settling velocity of particles in clusters can be significantly larger than than based upon averaging over all particles in the flow. Our results also show that this difference increases with increasing R λ , with a nonmonotonic dependence on Fr (since the results are zero for Fr = ∞ ). The average settling velocity of all particles is known to increase with increasing R λ , which is principally due to the enhanced range of scales for the particles to interact with as R λ grows [35]. However, the R λ dependence of the results in figure 16 would seem to be due to something more subtle to this, since all particles in the flow, whether clustered or not, would be subject to the enhanced range of scales available for the particles to interact with as R λ grows.

The results for the particle fluctuating kinetic energy K ( St ) ≡ (1 / 2) 〈‖ v p ( ) t -〈 v p ( ) t 〉‖ 〉 2 are shown in figure 17. The results show that while the settling velocity of particles in coherent clusters is significantly different from that of all the particles in the flow, there is only a small difference for their kinetic energy. In the absence of gravity, the particles in coherent clusters have slightly less kinetic energy compared with the average based on all particles in the flow. The most likely explanation for this is that the difference is caused by the preferential sampling of the turbulent flow, which is known to be suppressed in the presence of gravity [23]. As explained in [17], the decrease of the kinetic energy with increasing St is due to the inertial particles filtering out fluctuations in the turbulent flow, while the increase of the kinetic energy with increasing R λ (for a fixed St and Fr ) is due to a reduction of the filtering effect owing to the increased timescale separation in the turbulent flow as R λ is increased. Finally, the decrease of the kinetic energy as Fr is decreased (for a fixed St and R λ ) is because the filtering effect is enhanced due to the timescale of the fluid velocity seen by the particle becoming shorter as the particles settle faster.

FIG. 17: Average kinetic energy of particles within the clusters (solid lines) and the whole particles in the domain (dashed lines).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00341v1-figure-539e177c-0f84-4194-b842-b480884b2319.jpg)

## IV. CONCLUSIONS

In this work, we have used DNS and three-dimensional Vorono¨ ı analysis to explore the local distribution, settling velocity and acceleration of inertial particles (in the point particle limit) in statistically stationary, isotropic turbulence. In these simulations, we independently varied the Taylor Reynolds number R λ ∈ [90 , 398], Froude number Fr ≡ a /g η ∈ [0 052 . , ∞ ] (where a η is the Kolmogorov acceleration, and g is the acceleration due to gravity), and Kolmogorov scale Stokes number St ≡ τ /τ p η ∈ [0 , 3]. Independently varying these parameters was not possible in previous experimental investigations on the problem. The average inter-particle distance was held fixed in all the DNS for consistency when comparing the Vorono¨ ı analysis results for different R λ .

In agreement with previous results using the Radial Distribution Function (RDF) to quantify particle clustering [17, 23], which is a global measure, we find that for small Vorono¨ ı volumes (corresponding to the most clustered particles), the behavior is strongly dependent upon St and Fr , but only weakly dependent upon R λ , unless St &gt; 1. However, larger Vorono¨ ı volumes (void regions) exhibit a much stronger dependence on R λ , even when St ≤ 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¨ ı volumes that has been previously reported. Our results also show that the standard deviation of the Vorono¨ ı volumes is dominated by the behavior of the void regions, rather than the clustered regions. As a result, the results show that the standard deviation of the Vorono¨ ı volumes is reduced by gravity at all St , even though the Probability Density Function (PDF) results show that gravity can enhance the small-scale clustering, as was also observed using the RDF analysis in [23]. This is because gravity seems to always suppress the void regions. This highlights that using the standard deviation of Vorono¨ ı volumes as a measure of particle clustering (which has been done in many previous studies) can be somewhat misleading, since its properties are dominated by void, rather than small-scale clustered regions of the flow. We also show that the PDF of Vorono¨ ı volumes exhibits a quasi-lognormal behavior over a certain range of volumes, as previously reported. However, we show that this is only approximate, and does not describe well Vorono¨ ı volumes that are sufficiently large or small. The validity of the log-normal assumption also depends in a rather complex way on R λ and Fr .

We find that the average settling velocities conditioned on the Vorono¨ ı volumes exhibits a non-monotonic dependence on the Vorono¨ ı volumes, with the largest contribution to the particle settling velocities being associated with increasingly larger Vorono¨ ı volumes as the settling parameter Sv ≡ St/Fr is increased. This non-monotonic behavior does not seem to have been previously reported, likely because the Sv considered were not sufficiently large. It can be explained by the recent work of [35] who show theoretically and numerically that as Sv is increased, the scales responsible for the particle settling velocities shift to larger scales.

Even though the globally averaged particle acceleration is zero in the system we are considering, 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. In the direction normal to gravity, and in the absence of gravity, the average particle accelerations are independent of the Vorono¨ ı volumes. The particle acceleration variance, however, is a convex function of the Vorono¨ ı volumes, with or without gravity, which seems to indicate a non-trivial relationship between the Vorono¨ ı volumes and the sizes of the turbulent flow scales. Results for the variance of the fluid acceleration at the inertial particle positions are of the order of the square of the Kolmogorov acceleration and depend only weakly on Vorono¨ ı volumes. These results call

into question the 'sweep-stick' mechanism for particle clustering in turbulence which would lead one to expect that clustered particles reside in regions where the fluid acceleration is zero [54].

We then consider the properties of particles in clusters, which are regions of connected Vorono¨ ı cells whose volume is less than a certain threshold. The results for the PDFs of the cluster volumes reveal significant self-similarity of the clusters. The standard deviation of the cluster volumes provides further insight and shows that the statistics of the cluster volumes depends only weakly on St , with a stronger dependance on Fr and R λ . The weak dependence upon St may however be due to the fact that given the average inter-particle distance in our simulations, the cluster sizes are comparable to scales of the flow where one might expect the effects of particle inertia to be weak anyway. Finally, we compared the average settling velocities of all particles in the flow with those in clusters, and showed that those in the clusters settle much faster, in agreement with previous work. However, we also find that this difference grows significantly with increasing R λ and exhibits a non-monotonic dependence on Fr . The kinetic energy of the particles, on the other hand, is almost the same for particles in and not in the clusters, especially in the presence of gravity.

## V. ACKNOWLEDGMENTS

This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1548562 [56]. Specifically, the Comet cluster was used under allocation CTS170009. The Comet cluster and the Duke Computing Cluster (DCC) were used extensively for the computations in this study. The authors would like to express their sincere thanks and gratitude to Tom Milledge and Roberto E. Gonzalez for their generous helps with computational challenges of this work. We also would like to thank Yuanqing Liu for all his input during the discussion on the coding issues.

[1] Eric A D'Asaro. Turbulence in the upper-ocean mixed layer. Annual review of marine science , 6:101-115, 2014.

- [2] Kenneth V Beard and Harry T Ochs III. Warm-rain initiation: An overview of microphysical mechanisms. Journal of Applied Meteorology , 32(4):608-625, 1993.
- [3] T Esposti Ongaro, Carlo Cavazzoni, Giovanni Erbacci, Augusto Neri, and Maria-Vittoria Salvetti. A parallel multiphase flow code for the 3d simulation of explosive volcanic eruptions. Parallel Computing , 33(7-8):541-560, 2007.
- [4] M. R. Maxey. The gravitational settling of aerosol particles in homogeneous turbulence and random flow fields. J. Fluid Mech. , 174:441-465, 1987.
- [5] Vincenzo Bianco, Oronzio Manca, and Sergio Nardini. Numerical investigation on nanofluids turbulent convection heat transfer inside a circular tube. International Journal of Thermal Sciences , 50(3):341-349, 2011.
- [6] MR Momenifar, MA Akhavan-Behabadi, M Nasr, and P Hanafizadeh. Effect of lubricating oil on flow boiling characteristics of r-600a/oil inside a horizontal smooth tube. Applied Thermal Engineering , 91:62-72, 2015.
- [7] Kenneth Kuan-yun Kuo and Ragini Acharya. Fundamentals of turbulent and multiphase combustion . John Wiley &amp; Sons, 2012.
- [8] Pedram Hanafizadeh, Mohammadreza Momenifar, A Nouri Geimassi, and S Ghanbarzadeh. Void fraction and wake analysis of a gas-liquid two-phase cross-flow. Multiphase Science and Technology , 26(4), 2014.
- [9] Mohammad Hossein Arabnejad, Ali Amini, Mohamed Farhat, and Rickard E Bensow. Numerical and experimental investigation of shedding mechanisms from leading-edge cavitation. International Journal of Multiphase Flow , 2019.
- [10] Saeid Vahidifar, Mohammad Reza Saffarian, and Ebrahim Hajidavalloo. Introducing the theory of successful settling in order to evaluate and optimize the sedimentation tanks. Meccanica , 53(14):3477-3493, 2018.
- [11] P Hanafizadeh, MR Mmomenifar, M Bastankhah, and MH Saidi. Thermosyphon analysis using the particle image velocimetry method (piv). MECHANICAL ENGINEERING SHARIF (SHARIF: MECHANICAL ENGINEERING) , 32-3(1):11-20, 2016.
- [12] HT Bi, N Ellis, IA Abba, and JR Grace. A state-of-the-art review of gas-solid turbulent fluidization. Chemical Engineering Science , 55(21):4789-4825, 2000.
- [13] J. Bec, L. Biferale, M. Cencini, A. S. Lanotte, S. Musacchio, and F. Toschi. Heavy particle concentration in turbulence at dissipative and inertial scales. Phys. Rev. Lett. , 98:084502,

2007.

- [14] S. Sundaram and L. R. Collins. Collision statistics in an isotropic, particle-laden turbulent suspension I. Direct numerical simulations. J. Fluid Mech. , 335:75-109, 1997.
- [15] Federico Toschi and Eberhard Bodenschatz. Lagrangian properties of particles in turbulence. Annual review of fluid mechanics , 41:375-404, 2009.
- [16] Raymond A Shaw. Particle-turbulence interactions in atmospheric clouds. Annual Review of Fluid Mechanics , 35(1):183-227, 2003.
- [17] Peter J Ireland, Andrew D Bragg, and Lance R Collins. The effect of reynolds number on inertial particle dynamics in isotropic turbulence. part 1. simulations without gravitational effects. Journal of Fluid Mechanics , 796:617-658, 2016.
- [18] S. Balachandar and J. K. Eaton. Turbulent dispersed multiphase flow. Annu. Rev. Fluid Mech. , 42:111-133, 2010.
- [19] Alberto Aliseda, Alain Cartellier, F Hainaux, and Juan C Lasheras. Effect of preferential concentration on the settling velocity of heavy particles in homogeneous isotropic turbulence. Journal of Fluid Mechanics , 468:77-105, 2002.
- [20] Romain Monchaux, Micka¨ el Bourgoin, and Alain Cartellier. Preferential concentration of heavy particles: a vorono¨ ı analysis. Physics of Fluids , 22(10):103304, 2010.
- [21] J´ er´mie e Bec, Luca Biferale, Guido Boffetta, Massimo Cencini, Stefano Musacchio, and Federico Toschi. Lyapunov exponents of heavy particles in turbulence. Physics of Fluids , 18(9):091702, 2006.
- [22] Enrico Calzavarini, Martin Kerscher, Detlef Lohse, and Federico Toschi. Dimensionality and morphology of particle and bubble clusters in turbulent flow. Journal of fluid mechanics , 607:13-24, 2008.
- [23] Peter J Ireland, Andrew D Bragg, and Lance R Collins. The effect of reynolds number on inertial particle dynamics in isotropic turbulence. part 2. simulations with gravitational effects. Journal of Fluid Mechanics , 796:659-711, 2016.
- [24] Yoshiyuki Tagawa, Julin Martnez Mercado, Vivek N. Prakash, Enrico Calzavarini, Chao Sun, and Detlef Lohse. Three-dimensional lagrangian vorono analysis for clustering of particles and bubbles in turbulence. Journal of Fluid Mechanics , 693:201215, 2012.
- [25] Alec J Petersen, Lucia Baker, and Filippo Coletti. Experimental study of inertial particles clustering and settling in homogeneous turbulence. Journal of Fluid Mechanics , 864:925-970,

2019.

- [26] Mohammadreza Momenifar, Rohit Dhariwal, and Andrew D Bragg. Influence of reynolds number on the motion of settling, bidisperse inertial particles in turbulence. Physical Review Fluids , 4(5):054301, 2019.
- [27] Martin Obligado, Tomas Teitelbaum, Alain Cartellier, Pablo Mininni, and Micka¨l Bourgoin. e Preferential concentration of heavy particles in turbulence. Journal of Turbulence , 15(5):293310, 2014.
- [28] Sholpan Sumbekova, Alain Cartellier, Alberto Aliseda, and Mickael Bourgoin. Preferential concentration of inertial sub-kolmogorov particles: the roles of mass loading of particles, stokes numbers, and reynolds numbers. Physical Review Fluids , 2(2):024302, 2017.
- [29] GRETCHEN L. Holtzer and LANCE R. Collins. Relationship between the intrinsic radial distribution function for an isotropic field of particles and lower-dimensional measurements. Journal of Fluid Mechanics , 459:93?102, 2002.
- [30] Lucia Baker, Ari Frankel, Ali Mani, and Filippo Coletti. Coherent clusters of inertial particles in homogeneous turbulence. Journal of Fluid Mechanics , 833:364-398, 2017.
- [31] A Dejoan and R Monchaux. Preferential concentration and settling of heavy particles in homogeneous turbulence. Physics of Fluids , 25(1):013301, 2013.
- [32] R Monchaux and A Dejoan. Settling velocity and preferential concentration of heavy particles under two-way coupling effects in homogeneous turbulence. Physical Review Fluids , 2(10):104302, 2017.
- [33] Thorsten Bosse, Leonhard Kleiser, and Eckart Meiburg. Small particles in homogeneous turbulence: settling velocity enhancement by two-way coupling. Physics of Fluids , 18(2):027102, 2006.
- [34] Anne Dejoan. Dns experiments on the settling of heavy particles in homogeneous turbulence: two-way coupling and reynolds number effects. In Journal of Physics: Conference Series , volume 333, page 012006. IOP Publishing, 2011.
- [35] Josin Tom and Andrew D. Bragg. Multiscale preferential sweeping of particles settling in turbulence. Journal of Fluid Mechanics , 871:244?270, 2019.
- [36] H Parishani, O Ayala, B Rosa, L-P Wang, and WW Grabowski. Effects of gravity on the acceleration and pair statistics of inertial particles in homogeneous isotropic turbulence. Physics of Fluids , 27(3):033304, 2015.

- [37] Rohit Dhariwal and Andrew D Bragg. Small-scale dynamics of settling, bidisperse particles in turbulence. Journal of Fluid Mechanics , 839:594-620, 2018.
- [38] Rohit Dhariwal and Andrew D. Bragg. Enhanced and suppressed multiscale dispersion of bidisperse inertial particles due to gravity. Phys. Rev. Fluids , 4:034302, Mar 2019.
- [39] Martin R Maxey and James J Riley. Equation of motion for a small rigid sphere in a nonuniform flow. The Physics of Fluids , 26(4):883-889, 1983.
- [40] Peter J Ireland, T Vaithianathan, Parvez S Sukheswalla, Baidurja Ray, and Lance R Collins. Highly parallel particle-laden flow solver for turbulence research. Computers &amp; Fluids , 76:170177, 2013.
- [41] Romain Monchaux. Measuring concentration with vorono¨ ı diagrams: the study of possible biases. New Journal of Physics , 14(9):095013, 2012.
- [42] Christopher Nilsen, Helge I Andersson, and Lihao Zhao. A vorono¨ ı analysis of preferential concentration in a vertical channel flow. Physics of Fluids , 25(11):115108, 2013.
- [43] B Rabencov and R van Hout. Voronoi analysis of beads suspended in a turbulent square channel flow. International Journal of Multiphase Flow , (68):10-13, 2015.
- [44] Wenjun Yuan, Lihao Zhao, Helge I Andersson, and Jianqiang Deng. Three-dimensional vorono¨ ı analysis of preferential concentration of spheroidal particles in wall turbulence. Physics of Fluids , 30(6):063304, 2018.
- [45] J´ arai-Szab´ Ferenc and Zolt´n N´da. On the size distribution of poisson voronoi cells. o a e Physica A: Statistical Mechanics and its Applications , 385(2):518-526, 2007.
- [46] Martin Obligado, Mahrane Missaoui, Romain Monchaux, Alain Cartellier, and Micka¨l Boure goin. Reynolds number influence on preferential concentration of heavy particles in turbulent flows. In Journal of Physics: Conference Series , volume 318, page 052015. IOP Publishing, 2011.
- [47] A. D. Bragg, P. J. Ireland, and L. R. Collins. On the relationship between the non-local clustering mechanism and preferential concentration. Journal of Fluid Mechanics , 780:327343, 10 2015.
- [48] Uriel Frisch. Turbulence: The Legacy of A. N. Kolmogorov . Cambridge University Press, 1995.
- [49] A.D. Bragg and L.R. Collins. New insights from comparing statistical theories for inertial particles in turbulence: I. spatial distribution of particles. New J. Phys. , 16:055013, 2014.

- [50] J. Bec, L. Biferale, G. Boffetta, A. Celani, M. Cencini, A. S. Lanotte, S. Musacchio, and F. Toschi. Acceleration statistics of heavy particles in turbulence. J. Fluid Mech. , 550:349358, 2006.
- [51] Andrew D. Bragg, Peter J. Ireland, and Lance R. Collins. Mechanisms for the clustering of inertial particles in the inertial range of isotropic turbulence. Phys. Rev. E , 92:023029, Aug 2015.
- [52] Ari Frankel, Hadi Pouransari, Filippo Coletti, and Ali Mani. Settling of heated particles in homogeneous turbulence. Journal of Fluid Mechanics , 792:869-893, 2016.
- [53] L. P. Wang and M. R. Maxey. Settling velocity and concentration distribution of heavy particles in homogeneous isotropic turbulence. J. Fluid Mech. , 256:27-68, 1993.
- [54] S. W. Coleman and J. C. Vassilicos. A unified sweep-stick mechanism to explain particle clustering in two- and three-dimensional homogeneous, isotropic turbulence. Physics of Fluids , 21(11):113301, 2009.
- [55] Hiroshi Yoshimoto and Susumu Goto. Self-similar clustering of inertial particles in homogeneous turbulence. Journal of Fluid Mechanics , 577:275-286, 2007.
- [56] J. Towns, T. Cockerill, M. Dahan, I. Foster, K. Gaither, A. Grimshaw, V. Hazlewood, S. Lathrop, D. Lifka, G. D. Peterson, R. Roskies, J. R. Scott, and N. Wilkins-Diehr. Xsede: Accelerating scientific discovery. Computing in Science &amp; Engineering , 16(5):62-74, Sept.-Oct. 2014.
Mohammadreza Momenifar Andrew D. Bragg
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
"..."