Direct numerical simulation is used to investigate effects of turbulent flow in the confined geometry of a facecentred cubic porous unit cell on the transport, clustering and deposition of fine particles at different Stokes numbers ( đť‘†đť‘ˇ=0.01,0.1,0.5,1,2 ) and at a pore Reynolds number of 500. Particles are advanced using oneway coupling and the collision of particles with pore walls is modelled as perfectly elastic with specular reflection. Tools for studying inertial particle dynamics and clustering developed for homogeneous flows are adapted to take into account the embedded, curved geometry of the pore walls. The pattern and dynamics of clustering are investigated using the volume change of Voronoi tesselation in time to analyse the divergence and convergence of the particles. Similar to the case of homogeneous, isotropic turbulence, the cluster formation is present at large volumes, while cluster destruction is prominent at small volumes and these effects are amplified with the Stokes number. However, unlike homogeneous, isotropic turbulence, the formation of a large number of very small volumes was observed at all Stokes numbers and attributed to the collision of particles with the pore wall. Multiscale wavelet analysis of the particle number density indicates that the peak of the energy density spectrum, representative of enhanced particle clustering, shifts towards larger scales with an increase in the Stokes number. Scaledependent skewness and flatness quantify the intermittent void and cluster distribution, with cluster formation observed at small scales for all Stokes numbers, and void regions at large scales for large Stokes numbers.
more »
« less
Clustering of inertial particles in turbulent flow through a porous unit cell
Direct numerical simulation is used to investigate effects of turbulent flow in the confined geometry of a facecentred cubic porous unit cell on the transport, clustering and deposition of fine particles at different Stokes numbers ( $St = 0.01, 0.1, 0.5, 1, 2$ ) and at a pore Reynolds number of 500. Particles are advanced using oneway coupling and the collision of particles with pore walls is modelled as perfectly elastic with specular reflection. Tools for studying inertial particle dynamics and clustering developed for homogeneous flows are adapted to take into account the embedded, curved geometry of the pore walls. The pattern and dynamics of clustering are investigated using the volume change of Voronoi tesselation in time to analyse the divergence and convergence of the particles. Similar to the case of homogeneous, isotropic turbulence, the cluster formation is present at large volumes, while cluster destruction is prominent at small volumes and these effects are amplified with the Stokes number. However, unlike homogeneous, isotropic turbulence, the formation of a large number of very small volumes was observed at all Stokes numbers and attributed to the collision of particles with the pore wall. Multiscale wavelet analysis of the particle number density indicates that the peak of the energy density spectrum, representative of enhanced particle clustering, shifts towards larger scales with an increase in the Stokes number. Scaledependent skewness and flatness quantify the intermittent void and cluster distribution, with cluster formation observed at small scales for all Stokes numbers, and void regions at large scales for large Stokes numbers.
more »
« less
 Award ID(s):
 2053248
 NSFPAR ID:
 10324166
 Date Published:
 Journal Name:
 Journal of Fluid Mechanics
 Volume:
 937
 ISSN:
 00221120
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this


The relative velocities and positions of monodisperse highinertia particle pairs in isotropic turbulence are studied using direct numerical simulations (DNS), as well as Langevin simulations (LS) based on a probability density function (PDF) kinetic model for pair relative motion. In a prior study (Rani etÂ al. , J. Fluid Mech. , vol. 756, 2014, pp. 870â€“902), the authors developed a stochastic theory that involved deriving closures in the limit of high Stokes number for the diffusivity tensor in the PDF equation for monodisperse particle pairs. The diffusivity contained the time integral of the Eulerian twotime correlation of fluid relative velocities seen by pairs that are nearly stationary. The twotime correlation was analytically resolved through the approximation that the temporal change in the fluid relative velocities seen by a pair occurs principally due to the advection of smaller eddies past the pair by largescale eddies. Accordingly, two diffusivity expressions were obtained based on whether the pair centre of mass remained fixed during flow time scales, or moved in response to integralscale eddies. In the current study, a quantitative analysis of the (Rani etÂ al. 2014) stochastic theory is performed through a comparison of the pair statistics obtained using LS with those from DNS. LS consist of evolving the Langevin equations for pair separation and relative velocity, which is statistically equivalent to solving the classical Fokkerâ€“Planck form of the pair PDF equation. Langevin simulations of particlepair dispersion were performed using three closure forms of the diffusivity â€“ i.e. the one containing the time integral of the Eulerian twotime correlation of the seen fluid relative velocities and the two analytical diffusivity expressions. In the first closure form, the twotime correlation was computed using DNS of forced isotropic turbulence laden with stationary particles. The two analytical closure forms have the advantage that they can be evaluated using a model for the turbulence energy spectrum that closely matched the DNS spectrum. The three diffusivities are analysed to quantify the effects of the approximations made in deriving them. Pair relativemotion statistics obtained from the three sets of Langevin simulations are compared with the results from the DNS of (moving) particleladen forced isotropic turbulence for $St_{\unicode[STIX]{x1D702}}=10,20,40,80$ and $Re_{\unicode[STIX]{x1D706}}=76,131$ . Here, $St_{\unicode[STIX]{x1D702}}$ is the particle Stokes number based on the Kolmogorov time scale and $Re_{\unicode[STIX]{x1D706}}$ Â is the Taylor microscale Reynolds number. Statistics such as the radial distribution function (RDF), the variance and kurtosis of particlepair relative velocities and the particle collision kernel were computed using both Langevin and DNS runs, and compared. The RDFs from the stochastic runs were in good agreement with those from the DNS. Also computed were the PDFs $\unicode[STIX]{x1D6FA}(Ur)$ and $\unicode[STIX]{x1D6FA}(U_{r}r)$ of relative velocity $U$ and of the radial component of relative velocity $U_{r}$ respectively, both PDFs conditioned on separation $r$ . The first closure form, involving the Eulerian twotime correlation of fluid relative velocities, showed the best agreement with the DNS results for the PDFs.more » « less

null (Ed.)We investigate the dynamics of cohesive particles in homogeneous isotropic turbulence, based on oneway coupled simulations that include Stokes drag, lubrication, cohesive and direct contact forces. We observe a transient flocculation phase, followed by a statistically steady equilibrium phase. We analyse the temporal evolution of floc size and shape due to aggregation, breakage and deformation. Larger turbulent shear and weaker cohesive forces yield smaller elongated flocs. Flocculation proceeds most rapidly when the fluid and particle time scales are balanced and a suitably defined Stokes number is $O(1)$ . During the transient stage, cohesive forces of intermediate strength produce flocs of the largest size, as they are strong enough to cause aggregation, but not so strong as to pull the floc into a compact shape. Small Stokes numbers and weak turbulence delay the onset of the equilibrium stage. During equilibrium, stronger cohesive forces yield flocs of larger size. The equilibrium floc size distribution exhibits a preferred size that depends on the cohesive number. We observe that flocs are generally elongated by turbulent stresses before breakage. Flocs of size close to the Kolmogorov length scale preferentially align themselves with the intermediate strain direction and the vorticity vector. Flocs of smaller size tend to align themselves with the extensional strain direction. More generally, flocs are aligned with the strongest Lagrangian stretching direction. The Kolmogorov scale is seen to limit floc growth. We propose a new flocculation model with a variable fractal dimension that predicts the temporal evolution of the floc size and shape.more » « less

null (Ed.)Particle pair statistics of inertial particles having average Stokes numbers of 2.1 and 14 are measured in isotropic turbulence at a Reynolds number of ReÎ» = 240. The radial distribution function (RDF) and mean relative approach velocity are obtained at small separation distances using 2frame stereoscopic particle tracking velocimetry (stereoPTV). At small separation distance, the RDF varies by an order of magnitude in the range of Stokes numbers investigated. However, the mean relative approach velocity is found to have a weak dependence on Stokes number. The results are shown to have high accuracy when compared to analogous monoPTV datasets, and can be used to provide a more reliable estimate of the interparticle collision rate. The main limitation of the measurement is observed at separation distances less than the laser sheet thickness, where the technique tended to underestimate the mean relative approach velocity.more » « less

null (Ed.)The effect of turbulence on snow precipitation is not incorporated into present weather forecasting models. Here we show evidence that turbulence is in fact a key influence on both fall speed and spatial distribution of settling snow. We consider three snowfall events under vastly different levels of atmospheric turbulence. We characterize the size and morphology of the snow particles, and we simultaneously image their velocity, acceleration and relative concentration over vertical planes approximately $30\ \textrm {m}^2$ in area. We find that turbulencedriven settling enhancement explains otherwise contradictory trends between the particle size and velocity. The estimates of the Stokes number and the correlation between vertical velocity and local concentration are consistent with the view that the enhanced settling is rooted in the preferential sweeping mechanism. When the snow vertical velocity is large compared to the characteristic turbulence velocity, the crossing trajectories effect results in strong accelerations. When the conditions of preferential sweeping are met, the concentration field is highly nonuniform and clustering appears over a wide range of scales. These clusters, identified for the first time in a naturally occurring flow, display the signature features seen in canonical settings: powerlaw size distribution, fractallike shape, vertical elongation and large fall speed that increases with the cluster size. These findings demonstrate that the fundamental phenomenology of particleladen turbulence can be leveraged towards a better predictive understanding of snow precipitation and ground snow accumulation. They also demonstrate how environmental flows can be used to investigate dispersed multiphase flows at Reynolds numbers not accessible in laboratory experiments or numerical simulations.more » « less