On the influence of forced homogeneous-isotropic turbulence on the settling and clustering of finite-size particles

摘要

We investigate the motion of heavy particles with a diameter of several multiples of the Kolmogorov length scale in the presence of forced turbulence and gravity, resorting to interface-resolved direct numerical simulation based on an immersed boundary method. The values of the particles' relative density (1.5) and of the Galileo number (180) are such that strong wake-induced particle clustering would occur in the absence of turbulence (Uhlmann and Doychev in J Fluid Mech 752:310-348, 2014. 10.1017/jfm.2014.330). The forced turbulence in the two present cases (with Taylor-scale Reynolds number 95 and 140) would lead to mild levels of clustering in the absence of gravity (Uhlmann and Chouippe in J Fluid Mech 812:991-1023, 2017. 10.1017/jfm.2016.826). Here we detect a tendency to cluster with an intensity (quantified via the standard deviation of the distribution of Voronoi cell volumes) which is intermediate between these two limiting cases, meaning that forced background turbulence decreases the level of clustering otherwise observed under ambient settling. However, the clustering strength does not monotonically decay with the relative turbulence intensity. Various mechanisms by which coherent structures can affect particle motion are discussed. It is argued that the reduced interaction time due to particle settling through the surrounding eddy (crossing trajectories) has the effect of shifting upward the range of eddies with a time-scale matching the characteristic time-scale of the particle. In the present cases this shift might bring the particles into resonance with the energetic eddies of the turbulent spectrum. Concerning the average particle settling velocity, we find very small deviations (of the order of one percent) from the value obtained for an isolated particle in ambient fluid when defining the relative velocity as an apparent slip velocity (i.e., as the difference between the averages computed separately for the velocities of each phase). This is cons