Equation (2.10) in Costa, Brandt & Picano (2020) for the lift force model used in the point-particle direct numerical simulations (DNS), and which is derived from the classical lift force of Saffman (1965), (Equation presented) does not correspond to the force model actually used in the point-particle DNS with lift force presented in the manuscript. Instead, the following equation was used: (Equation presented) which replaces the first occurrence of the term |Us| on the right-hand-side of (1) with |ω|D. We recall that two cases were considered in the manuscript depending on the value of J in the lift force equation: J = 1 in the case denoted PP-Saffman; and J given by (Equation presented) with ϵ = √|ω|ν/|Us|, in the case denoted PP-McLaughlin. Also, equation (2.13) of the manuscript - describing the perfectly elastic hard-sphere rebound - is incorrect; the term D/2 should be D: (Equation presented) Despite the lapse in the manuscript, equation (4) was implemented correctly (Costa et al. 2020). The results from the point-particle DNS with the model reported in (2.10) of Costa et al. (2020) ((1) above) differ from those reported in the manuscript, and are shown (Figure presented) in figure 1 (cf. figures 7 and 8 of Costa et al. (2020)). The statistics presented here have been collected in the fully developed state from 600 samples over a time interval of 250h/Ub, which ensured statistical convergence of the results. The results from the point-particle cases presented in the original manuscript are also reproduced here with this (higher) statistical sampling, and show very minor differences with respect to figures 7 and 8 of Costa et al. (2020). In light of these results, the conclusions drawn from the results in the last section of § 3 of the manuscript must be therefore reformulated: (i) The Saffman lift model does not correctly predict the near-wall statistics of the interface-resolved DNS very close to the wall, including the near-wall concentration peak. (ii) The equation proposed by Mei (1992) that fits the model of McLaughlin (1991) shows results similar to those reported in the original manuscript for this model. That is, it predicts well the near-wall concentration peak, and fails to predict the other observables near the wall. (iii) Equation (2) for Fl presented above, with J = 1, predicts very well all the observables in figure 1. We have therefore accidentally discovered that the expression (2) for Fl predicts the observed particle statistics very well. Still, the reason for the strikingly good agreement remains elusive to us. We hope that this result can be further exploited for the improvement lift force models for point-particle simulations of wall-bounded turbulent flows.@en

We consider suspensions of finite-size neutrally buoyant rigid spherical particles in channel flow and investigate the relevance of different momentum transfer mechanisms and the relation between the local particle dynamics and the bulk flow properties in the highly inertial regime. Interface-resolved simulations are performed in the range of Reynolds numbers and solid volume fractions. The Lagrangian particle statistics show that pair interactions are highly inhomogeneous and dependent on the distance from the wall: in their vicinity, the underlying mean shear drives the pair interactions, while a high degree of isotropy, dictated by more frequent collisions, characterizes the core region. Analysis of the momentum balance reveals that while the particle-induced stresses govern the dynamics in dense conditions, and moderate Reynolds numbers, <![CDATA[$Re, the turbulent stresses take over at higher Reynolds numbers. This behaviour is associated with a reduced particle migration toward the channel core, which decreases the importance of the particle-induced stress and increases the turbulent activity. Our results indicate that Reynolds stresses and the associated velocity fluctuations, characteristics of near-wall turbulence, prevail at high inertia over the resistance to deformation presented by the particles for volume fractions lower than 30 %. @en

We use interface-resolved simulations to study near-wall turbulence modulation by small inertial particles, much denser than the fluid, in dilute/semi-dilute conditions. We considered three bulk solid mass fractions, and, with only the latter two showing turbulence modulation. The increase of the drag is strong at, but mild in the densest case. Two distinct regimes of turbulence modulation emerge: for smaller mass fractions, the turbulence statistics are weakly affected and the near-wall particle accumulation increases the drag so the flow appears as a single-phase flow at slightly higher Reynolds number. Conversely, at higher mass fractions, the particles modulate the turbulent dynamics over the entire flow, and the interphase coupling becomes more complex. In this case, fluid Reynolds stresses are attenuated, but the inertial particle dynamics near the wall increases the drag via correlated velocity fluctuations, leading to an overall drag increase. Hence, we conclude that, although particles at high mass fractions reduce the fluid turbulent drag, the solid phase inertial dynamics still increases the overall drag. However, inspection of the streamwise momentum budget in the two-way coupling limit of vanishing volume fraction, but finite mass fraction, indicates that this trend could reverse at even higher particle load. @en

We present a direct comparison between interface-resolved and one-way-coupled point-particle direct numerical simulations (DNS) of gravity-free turbulent channel flow laden with small inertial particles, with high particle-to-fluid density ratio and diameter of approximately three viscous units. The most dilute flow considered, solid volume fraction , shows the particle feedback on the flow to be negligible, whereas differences with respect to the unladen case, notably a drag increase of approximately 10 %, are found for a volume fraction. This is attributed to a dense layer of particles at the wall, caused by turbophoresis, flowing with large particle-to-fluid apparent slip velocity. The most dilute case is therefore taken as the benchmark for assessing the validity of a widely used point-particle model, where the particle dynamics results only from inertial and nonlinear drag forces. In the bulk of the channel, the first- and second-order moments of the particle velocity from the point-particle DNS agree well with those from the interface-resolved DNS. Close to the wall, however, most of the statistics show major qualitative differences. We show that this difference originates from the strong shear-induced lift force acting on the particles in the near-wall region. This mechanism is well captured by the lift force model due to Saffman (J. Fluid Mech., vol. 22 (2), 1965, pp. 385-400), while other widely used, more elaborate, approaches aiming at extending the lift model for a wider range of particle Reynolds numbers can actually underpredict the magnitude of the near-wall particle velocity fluctuations for the cases analysed here.@en

We study turbulent channel flow of a binary mixture of finite-sized neutrally buoyant rigid particles by means of interface-resolved direct numerical simulations. We fix the bulk Reynolds number and total solid volume fraction, Reb=5600 and Φ=20% , and vary the relative fraction of small and large particles. The binary mixture consists of particles of two different sizes, 2h/dl=20 and 2h/ds=30 where h is the half-channel height and dl and ds the diameters of the large and small particles. While the particulate flow statistics exhibit a significant alteration of the mean velocity profile and turbulent fluctuations with respect to the unladen flow, the differences between the mono-disperse and bi-disperse cases are small. However, we observe a clear segregation of small particles at the wall in binary mixtures, which affects the dynamics of the near-wall region and thus the overall drag. This results in a higher drag in suspensions with a larger number of large particles. As regards bi-disperse effects on the particle dynamics, a non-monotonic variation of the particle dispersion in the spanwise (homogeneous) direction is observed when increasing the percentage of small/large particles. Finally, we note that particles of the same size tend to cluster more at contact whereas the dynamics of the large particles gives the highest collision kernels due to a higher approaching speed.@en

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.@en

Dense suspensions are usually investigated in the laminar limit where inertial effects are insignificant. In this regime, the main effect of the suspended phase is to alter the rheological behavior of the flow which always displays higher effective viscosity with respect to the carrier fluid. When the flow rate is high enough, i.e. at high Reynolds number, the flow may become turbulent and the interaction between solid and liquid phase modifies the turbulent dynamics that we know in single-phase fluids. In the present work, we study turbulent channel flows laden with finite-size particles at high volume fraction (F = 0:2) by means of Direct Numerical Simulations. A direct-forcing Immersed Boundary Method has been adopted to couple liquid and solid phases. The two-phase simulations have been performed fixing the bulk Reynolds number at Reb = Ub 2h=n = 12000 (Ub bulk velocity, h channel half-width and n the fluid kinematic viscosity). The particle size is relatively large with respect to the viscous length, i.e. 10 and 20 times, but smaller than large scales. We will present a detailed comparison of the statistical behavior of the particle-laden flow and the corresponding single-phase flow. The presence of the solid phase strongly alters the wall turbulence dynamics and its effect cannot be accounted only considering the higher rheological effective viscosity.@en

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.@en

Suspensions of solid particles are frequently found in applications and environmental flows. Several studies concern the rheological properties of suspensions in laminar flows, but much less is known of turbulent suspensions. The present work fills this gap providing DNS data on dense suspensions of neutrally-buoyant rigid sphere in a turbulent channel flow at the bulk Reynolds number of Re = U0h/ν = 2800. We show that considering volume fractions Φ ≤ 0.1 the turbulent flow is similar to the unladen case with higher turbulence intensities. On the contrary, the flow behavior strongly changes at Φ = 0.2where the turbulence appears to be attenuated.@en

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.@en

Inertial regimes in a channel flow of suspension of finite-size neutrally buoyant particles are studied for a wide range of Reynolds numbers: 500 Re 5000, and particle volume fractions: 0 0:3. The flow is classified in three different regimes according to the phase-averaged stress budget across the channel [2]. The laminar viscous regime at low Re and - where the viscous stress is the dominating term in the budget, the turbulent regime at high Re and relatively low where the momentum is mainly transferred by the action of the Reynolds stress and the inertial shear-thickening regime where the particle stress contributes the most to the significant enhancement of the wall shear stress. Particle distribution and dispersion properties provide additional evidence for the existence of the three different regimes.@en

Dense suspensions are usually investigated in the laminar limit where inertial effects are insignificant. When the flow rate is high enough, i.e. at high Reynolds number, the flow may become turbulent and the interaction between solid and liquid phases modifies the turbulence we know in single-phase fluids. In the present work, we study turbulent channel flows laden with finite-size particles at high volume fraction by means of Direct Numerical Simulations. A direct-forcing Immersed Boundary Method has been adopted to couple liquid and solid phases. We will show that the turbulence is attenuated in dense cases, even though the overall drag is increased because of the particle contribution to the total stress.@en