Kolmogorov-size particles in homogeneous and isotropic turbulence

TL;DR

This study compares detailed and simplified simulations of particles in turbulence, revealing how particle density affects flow behavior, energy spectra, and clustering, with implications for modeling particle-laden flows.

Contribution

It provides a comprehensive comparison between interface-resolved and point-particle simulations, highlighting their accuracy and limitations in modeling particle-turbulence interactions.

Findings

Energy spectrum follows -5/3 and -4 scaling laws.

Particles favor axial strain and vortex compression events.

Point-particle simulations underpredict clustering at higher densities.

Abstract

We investigate the fluid-solid interaction of suspensions of Kolmogorov-size spherical particles moving in homogeneous isotropic turbulence at a microscale Reynolds number of $Re_\lambda \approx 140$. Two volume fractions are considered, $10^{-5}$ and $10^{-3}$, and the solid-to-fluid density ratio is set to $5$ and $100$. We present a comparison between interface-resolved (PR-DNS) and one-way-coupled point-particle (PP-DNS) direct numerical simulations. We find that the modulated energy spectrum shows the classical $-5/3$ Kolmogorov scaling in the inertial range of scales and a $-4$ scaling at smaller scales, with the latter resulting from a balance between the energy injected by the particles and the viscous dissipation, in an otherwise smooth flow. An analysis of the small-scale flow topology shows that the particles mainly favour events with axial strain and vortex compression. The dynamics of the particles and their collective motion studied for PR-DNS are used to assess the validity of the PP-DNS. We find that the PP-DNS predicts fairly well both the Lagrangian and Eulerian statistics of the particles motion for the low-density case, while some discrepancies are observed for the high-density case. Also, the PP-DNS is found to underpredict the level of clustering of the suspension compared to the PR-DNS, with a larger difference for the high-density case.