During the last decade, many approaches for resolved-particle simulation (RPS) have been developed for numerical studies of finite-size particle-laden turbulent flows. In this paper, three RPS approaches are compared for a particle-laden decaying turbulence case. These methods are, the Volume-of-Fluid Lagrangian method, based on the viscosity penalty method (VoF-Lag); a direct forcing Immersed Boundary Method, based on a regularized delta function approach for the fluid/solid coupling (IBM); and the Bounce Back scheme developed for Lattice Boltzmann method (LBM-BB). The physics and the numerical performances of the methods are analyzed. Modulation of turbulence is observed for all the methods, with a faster decay of turbulent kinetic energy compared to the single-phase case. Lagrangian particle statistics, such as the velocity probability density function and the velocity autocorrelation function, show minor differences among the three methods. However, major differences between the codes are observed in the evolution of the particle kinetic energy. These differences are related to the treatment of the initial condition when the particles are inserted in an initially single-phase turbulence. The averaged particle/fluid slip velocity is also analyzed, showing similar behavior as compared to the results referred in the literature. The computational performances of the different methods differ significantly. The VoF-Lag method appears to be computationally most expensive. Indeed, this method is not adapted to turbulent cases. The IBM and LBM-BB implementations show very good scaling.

We use interface-resolved numerical simulations to study finite-size effects in turbulent channel flow of neutrally buoyant spheres. Two cases with particle sizes differing by a factor of two, at the same solid volume fraction of 20 % and bulk Reynolds number are considered. These are complemented with two reference single-phase flows: the unladen case, and the flow of a Newtonian fluid with the effective suspension viscosity of the same mixture in the laminar regime. As recently highlighted in Costa et al. (Phys. Rev. Lett., vol. 117, 2016, 134501), a particle–wall layer is responsible for deviations of the mesoscale-averaged statistics from what is observed in the continuum limit where the suspension is modelled as a Newtonian fluid with (higher) effective viscosity. Here we investigate in detail the fluid and particle dynamics inside this layer and in the bulk. In the particle–wall layer, the near-wall inhomogeneity has an influence on the suspension microstructure over a distance proportional to the particle size. In this layer, particles have a significant (apparent) slip velocity that is reflected in the distribution of wall shear stresses. This is characterized by extreme events (both much higher and much lower than the mean). Based on these observations we provide a scaling for the particle-to-fluid apparent slip velocity as a function of the flow parameters. We also extend the scaling laws in Costa et al. (Phys. Rev. Lett., vol. 117, 2016, 134501) to second-order Eulerian statistics in the homogeneous suspension region away from the wall. The results show that finite-size effects in the bulk of the channel become important for larger particles, while negligible for lower-order statistics and smaller particles. Finally, we study the particle dynamics along the wall-normal direction. Our results suggest that single-point dispersion is dominated by particle–turbulence (and not particle–particle) interactions, while differences in two-point dispersion and collisional dynamics are consistent with a picture of shear-driven interactions.

We present interface-resolved numerical simulations of turbulent channel flow laden with non-spherical rigid and neutrally-buoyant particles. We first focus on the case of oblate particles of aspect ratio 1/3 at volume fractions up to 15% and show that the turbulent drag is decreasing when increasing the particle volume fraction although the effective viscosity of the suspension actually increases. We relate the observed drag reduction to turbulence attenuation and to particle migration away from the near-wall region. Particles tend to align parallel to the wall with rotation rates significantly lower than those reported for spheres. In the second part of the study, we examine the effect of the particle slenderness on the observed drag reduction and show that the drag increases for flatter particles.

The gravity-driven motion of rigid particles in a viscous fluid is relevant in many natural and industrial processes, yet this has mainly been investigated for spherical particles. We therefore consider the sedimentation of non-spherical (spheroidal) isolated and particle pairs in a viscous fluid via numerical simulations using the Immersed Boundary Method. The simulations performed here show that the critical Galileo number for the onset of secondary motions decreases as the spheroid aspect ratio departs from 1. Above this critical threshold, oblate particles perform a zigzagging motion whereas prolate particles rotate around the vertical axis while having their broad side facing the falling direction. Instabilities of the vortices in the wake follow when farther increasing the Galileo number. We also study the drafting-kissing-tumbling associated with the settling of particle pairs. We find that the interaction time increases significantly for non-spherical particles and, more interestingly, spheroidal particles are attracted from larger lateral displacements. This has important implications for the estimation of collision kernels and can result in increasing clustering in suspensions of sedimenting spheroids.

The macroscopic behavior of dense suspensions of neutrally buoyant spheres in turbulent plane channel flow is examined. We show that particles larger than the smallest turbulence scales cause the suspension to deviate from the continuum limit in which its dynamics is well described by an effective suspension viscosity. This deviation is caused by the formation of a particle layer close to the wall with significant slip velocity. By assuming two distinct transport mechanisms in the near-wall layer and the turbulence in the bulk, we define an effective wall location such that the flow in the bulk can still be accurately described by an effective suspension viscosity. We thus propose scaling laws for the mean velocity profile of the suspension flow, together with a master equation able to predict the increase in drag as a function of the particle size and volume fraction.