Erratum: Interface-resolved simulations of small inertial particles in turbulent channel flow (Journal of Fluid Mechanics (2020)883 (A54) DOI: 10.1017/jfm.2019.918)

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 |U_s| 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 ϵ = √|ω|ν/|U_s|, 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/U_b, 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 F_l presented above, with J = 1, predicts very well all the observables in figure 1. We have therefore accidentally discovered that the expression (2) for F_l 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.