We perform interface-resolved simulations to study the modulation of statistically steady-state homogeneous shear turbulence by neutrally buoyant finite-size particles. We consider two shapes, spheres and oblates, and various solid volume fractions, up to 20%. The results show that a statistically steady state is not exclusive to single-phase homogeneous shear turbulence as the production and dissipation rates of the turbulent kinetic energy are also statistically in balance in particle-laden cases. The turbulent kinetic energy shows a non-monotonic behaviour with increasing solid volume fraction: increasing turbulence attenuation up to a certain concentration of solid particles and then enhancement of the turbulent kinetic energy at higher concentrations. This behaviour is observed at lower volume fractions for oblate particles than for spheres. The attenuation of the turbulence activity at lower volume fractions is explained through the enhancement of the dissipation rate close to the surface of particles. At higher volume fractions, however, particle pair interactions induce regions of high Reynolds shear stress, resulting in the enhancement of the turbulence activity. We show that the oblate particles of the considered size have larger rotational rates than spheres with no preferential orientation. This is in contrast to previous studies in wall-bounded flows where preferential orientation close to the wall and reduced rotation rates result in turbulence attenuation and thus drag reduction. Our results shed some light on the effect of rigid particles, smaller than the near-wall turbulent structures but still comparable to the viscous length scale, on the dynamics of the equilibrium logarithmic layer in wall-bounded flows.

We present results of interface-resolved simulations of heat transfer in suspensions of finite-size neutrally-buoyant spherical particles for solid volume fractions up to 35% and bulk Reynolds numbers from 500 to 5600. An Immersed Boundary-Volume of Fluid method is used to solve the energy equation in the fluid and solid phase.

We relate the heat transfer to the regimes of particle motion previously identified, i.e. a viscous regime at low volume fractions and low Reynolds number, particle-laden turbulence at high Reynolds and moderate volume fraction and particulate regime at high volume fractions. We show that in the viscous dominated regime, the heat transfer is mainly due to thermal diffusion with enhancement due to the particle-induced fluctuations. In the turbulent-like regime, we observe the largest enhancement of the global heat transfer, dominated by the turbulent heat flux. In the particulate shear-thickening regime, however, the heat transfer enhancement decreases as mixing is quenched by the particle migration towards the channel core. As a result, a compact loosely-packed core region forms and the contribution of thermal diffusion to the total heat transfer becomes significant once again. The global heat transfer becomes, in these flows at volume fractions larger than 25%, lower than in single phase turbulence.

We use interface-resolved direct numerical simulations to study the dynamics of a single sediment particle in a turbulent open channel flow over a fixed porous bed. The relative strength of the gravitational acceleration, quantified by the Galileo number, is varied so as to reproduce the different modes of sediment transport - resuspension, saltation and rolling. The results show that the sediment dynamics at lower Galileo numbers (i.e. resuspension and saltation) are mainly governed by the mean flow. Here, the regime of motion can be predicted by the ratio between the gravity and the shear-induced boundary force. In these cases, the sediment particle rapidly takes off when exposed to the flow, and proceeds with an oscillatory motion. Increasing the Galileo number, the frequency of these oscillations increases and their amplitude decreases, until the transport mode switches from resuspension to saltation. In this case, the sediment travels by short successive collisions with the bed. Further increasing the Galileo number, the particle rolls without detaching from the bed. Differently from the previous modes, the motion is triggered by extreme turbulent events, and the particle response depends on the specific initial conditions, at fixed Reynolds number. The results reveal that close to the onset of sediment motion, only turbulent sweeps can effectively trigger the particle motion by increasing the stagnation pressure upstream. We show that for the parameters in this study, a criterion based on the streamwise flow-induced force can successfully predict the incipient movement.