<figure><img src="image_1.png" /></figure>

HAL
open science

# On the relevance of lift force modelling in turbulent wall flows with small inertial particles

Wei Gao, Pengyu Shi, Matteo Parsani, Pedro Costa

► To cite this version:

Wei Gao, Pengyu Shi, Matteo Parsani, Pedro Costa. On the relevance of lift force modelling in turbulent wall flows with small inertial particles. Journal of Fluid Mechanics, 2024, 988, pp.A47. 10.1017/jfm.2024.474. hal-04684938

HAL Id: hal-04684938

https://ut3-toulouseinp.hal.science/hal-04684938v1

Submitted on 3 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# On the relevance of lift force modelling in turbulent wall flows with small inertial particles

Wei Gao<sup>1,2,†</sup>, Pengyu Shi<sup>3,4</sup>, Matteo Parsani<sup>1,2</sup> and Pedro Costa<sup>5</sup>

<sup>1</sup>Mechanical Engineering, Physical Science and Engineering Division, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia

<sup>2</sup>Applied Mathematics and Computational Science, Computer Electrical and Mathematical Science and Engineering Division, Extreme Computing Research Center, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia

<sup>3</sup>Institut de Mécanique des Fluides de Toulouse (IMFT), Université de Toulouse, CNRS, Toulouse, France

<sup>4</sup>Helmholtz-Zentrum Dresden – Rossendorf, Institute of Fluid Dynamics, 01328 Dresden, Germany

<sup>5</sup>Process and Energy Department, TU Delft, Leegewaterstraat 39, 2628 CB Delft, The Netherlands

(Received 1 December 2023; revised 5 April 2024; accepted 10 May 2024)

In particle-laden turbulent wall flows, lift forces can influence the near-wall turbulence. This has been observed recently in particle-resolved simulations, which, however, are too expensive to be used in upscaled models. Instead, point-particle simulations have been the method of choice to simulate the dynamics of these flows during the last decades. While this approach is simpler, cheaper and physically sound for small inertial particles in turbulence, some issues remain. In the present work, we address challenges associated with lift force modelling in turbulent wall flows and the impact of lift forces in the near-wall flow. We performed direct numerical simulations of small inertial point particles in turbulent channel flow for fixed Stokes number and mass loading while varying the particle size. Our results show that the particle dynamics in the buffer region, causing the apparent particle-to-fluid slip velocity to vanish, raises major challenges for modelling lift forces accurately. While our results confirm that lift forces have little influence on particle dynamics for sufficiently small particle sizes, for inner-scaled diameters of order one and beyond, lift forces become quite important near the wall. The different particle dynamics under lift forces results in the modulation of streamwise momentum transport in the near-wall region. We analyse this lift-induced turbulence modulation for different lift force models, and the results indicate that realistic models are critical for particle-modelled simulations to correctly predict turbulence modulation by particles in the near-wall region.

**Key words:** multiphase and particle-laden flows

† Email address for correspondence: wei.gao@kaust.edu.sa
# 1. Introduction

Wall-bounded turbulent flows laden with inertial particles abound in environmental and industrial contexts, such as the transport of particulate matter in the atmosphere, sediment transport in rivers, the separation of fine particles within industrial cyclones, and fluidized bed reactors. The chaotic and multiscale nature of the fluid turbulence coupled with the particle dynamics results in fascinating phenomena, which, however, are still challenging to understand and model (Balachandar & Eaton 2010; Brandt & Coletti 2022).

In the past decades, many studies have been devoted to the dynamics of turbulent wall transport of small inertial particles, driven by their prevalent nature and rich physics. Traditionally, these systems have been classified by how the dispersed phase influences the overall flow behaviour (Elghobashi 1994). Precisely, one-way coupling (1WC) corresponds to very small particle loadings, so small that their influence on the dynamics of the suspending fluid turbulence is negligible; two-way coupling (2WC) denotes flow regimes where mass loading is high enough such that the particles modify flow observably, but particle-particle interactions are negligible; finally, regimes where both particle-particle and particle-fluid interactions influence the flow dynamics are grouped into the four-way coupling category.

In first-principles, fully resolved direct numerical simulations (DNS) of a particle-laden flow, the flow around each particle needs to be resolved explicitly (Balachandar & Eaton 2010; Maxey 2017) (so-called particle-resolved DNS, PR-DNS). While this approach is free from modelling assumptions for the dispersed phase dynamics, it is computationally expensive due to the explicit imposition of no-slip and no-penetration boundary conditions at the surface of many particles moving in a turbulent medium. This is incredibly challenging when particles are tiny and there is a scale separation between the particle size and the smallest (Kolmogorov) turbulence scale due to extreme resolution requirements. Fortunately, in this case, one may be able to resort to the so-called point-particle approximation (PP-DNS), where interphase coupling is considered to be localized to a point.

Indeed, the point-particle approximation has been the method of choice for simulating the particle dynamics in turbulent wall flows. In these cases, it is assumed that the local properties of an undisturbed flow at the particle position drive the dispersed-phase dynamics (Gatignol 1983; Maxey & Riley 1983). In the case of highly inertial particles (i.e. large particle-to-fluid density ratios), the so-called Maxey-Riley-Gatignol equations simplify to a drag force term (Arcen, Tanière & Oesterlé 2006), which nevertheless yields highly non-trivial particle dynamics, even in isotropic turbulence (Toschi & Bodenschatz 2009). In wall-bounded turbulent flows, the inhomogeneous turbulence results in even richer dynamics, with the particle distribution driven by the interplay between small-scale clustering, turbophoresis, and the interaction between the particles and near-wall turbulence structures (Reeks 1983; Soldati & Marchioli 2009; Sardina et al. 2012), resulting in very inhomogeneous particle concentrations peaking at the wall, with strong preferential sampling of low-speed regions, as reproduced in a plethora of numerical studies, such as Fessler, Kulick & Eaton (1994), Uijttewaal & Oliemans (1996), Marchioli et al. (2003, 2008), Kuerten (2006), Soldati & Marchioli (2009), Bernardini, Pirozzoli & Orlandi (2013) and Jie et al. (2022).

When 2WC effects are important, point-particle DNS must describe the back-reaction of the dispersed phase in the flow. This is a major challenge, as the point-particle dynamics is driven by a local undisturbed fluid velocity, while the local flow field is being disturbed by the particles (see e.g. Gualtieri et al. 2015). The classical approach, known as the particle-in-cell method, was developed by Crowe, Sharma & Stock (1977)
Relevance of lift in particulate turbulent wall flows

and is widely used – even in the present study – but it requires a sufficiently high number of particles per grid cell, and cannot reproduce simple benchmarks of a sedimenting sphere in a quiescent medium. Indeed, approaches for a consistent treatment are being actively investigated; see Gualtieri *et al.* (2015), Horwitz & Mani (2016, 2020), Ireland & Desjardins (2017) and Horwitz *et al.* (2022). Investigations of particle-laden turbulent flows in the two-way-coupling regime are found in e.g. Vreman *et al.* (2009), Zhao, Andersson & Gillissen (2010), Kuerten (2016), Richter & Sullivan (2014), Capecelatro, Desjardins & Fox (2018), Wang & Richter (2019) and Battista *et al.* (2023).

While employing the point-particle approximation for small inertial particles is
physically sound, validating the fidelity of the approximation in 1WC and 2WC conditions
remains a challenge. Experimental data are available (Eaton & Fessler 1994; Kaftori,
Hetsroni & Banerjee 1995). Still, there are few parameter-matched numerical studies
due to limitations in terms of Reynolds number and well-controlled experiments, and
only in recent years have efforts in this direction started to appear (Wang *et al.* 2019).
Fortunately, PR-DNS of small inertial particles in turbulence has become possible, thanks
to the continuous growth in available computer power and development of efficient
numerical methods, with the first direct comparisons between point-particle models
and particle-resolved simulations starting to appear for forced homogeneous isotropic
turbulence and decaying HIT with moving particles (Schneiders, Meinke & Schröder
2017a,b; Fröhlich *et al.* 2018; Mehrabadi *et al.* 2018), and turbulent channel flow (Horne
& Mahesh 2019; Costa, Brandt & Picano 2020a).

These results from PR-DNS of particles in wall-bounded turbulence confirm that a sole drag force may not suffice to accurately reproduce the particle dynamics even for relatively small particles with a large density ratio under 1WC conditions. Where the shear rate is high near the wall, lift forces are also important. This has already been suggested in early works using PP-DNS, since the work by McLaughlin (1989) (see also Botto et al. 2005). In this regard, while Arcen et al. (2006) reported that the lift force has a negligible impact on the dispersed phase statistics, Marchioli, Picciotto & Soldati (2007) and Shin et al. (2022) found that the inclusion of the lift force in PP-DNS can lead to weakened near-wall particle accumulation in upward and horizontal channels, respectively. McLaughlin (1989) showed that the inclusion of the lift force resulted in a higher deposition rate; Wang et al. (1997) found that neglecting the lift force resulted in a slight reduction in the deposition rate.

While these findings are not necessarily contradictory, as there are some variations in
the governing parameters in the different studies and lift force models, many questions
remain elusive. (1) Which form of lift force model is appropriate for reproducing with high
fidelity the dynamics of small inertial particles in a turbulent wall flow? (2) Under which
conditions do lift forces matter in the particle dynamics? (3) What are the consequences
of choosing an inaccurate lift force model in dispersed phase dynamics and near-wall
turbulence modulation?

In the recent direct comparison between PR-DNS and PP-DNS for small particles in turbulent channel flow by Costa et al. (2020a) (see also Costa, Brandt & Picano 2020b), it was shown that lift models are vital for reproducing the near-wall particle dynamics, at least for inertial particles with a size of the order of one viscous wall unit. Surprisingly, a modified Saffman lift model predicted the particle dynamics perfectly, with the Saffman force scaled by a normalized shear rate, while conventional lift models yield poorer predictions (Costa et al. 2020b). In addition to reducing the near-wall concentration peak, it was also seen that lift force causes a large increase in correlated streamwise-wall-normal particle velocity fluctuations. This quantity has dramatic drag-increasing effects for sufficiently high mass loading, as it modulates the streamwise momentum budget; see
Costa, Brandt & Picano (2021). This direct link between lift force and drag increase makes their accurate modelling crucial.

In the present work, we address the three questions presented above by performing DNS of turbulent channel flow laden with small inertial particles, using the point-particle approximation, with three different lift models, from the classical Saffman lift model to the one that perfectly predicts previous PR-DNS data. We consider different inner-scaled particle sizes $D^+ \sim 1-0.1$ in 1WC and 2WC conditions for a Stokes number that is known to feature strong wall accumulation and preferential concentration in low-speed regions. Our analysis shows that currently available lift models are bound to fail near the wall for small inertial particles due to a vanishing particle-to-fluid slip velocity. We then use 2WC simulations to illustrate how different lift force models can result in qualitatively different turbulence modulation (i.e. turbulence attenuation versus turbulence enhancement). Near-wall accumulation is still significantly reduced for the smallest particle size ($D^+ = 0.1$). Still, the lift-induced increase of correlated velocity fluctuations near the wall is negligible, thus lift force has few consequences in turbulence modulation.

This paper is organized as follows. Next, in § 2, the numerical method, lift force models and computational set-ups are described. Then in § 3, the effects of lift force models on the particle dynamics and near-wall accumulation are investigated using 1WC DNS, based on which we try to explore the reason for the failure of conventional lift force models. Following this, we demonstrate qualitatively the effect of lift force models on turbulence modulation and momentum transfer with 2WC DNS. Finally, conclusions are drawn in § 4.

## 2. Methodology

### 2.1. Governing equations and numerical method

The fluid phase is governed by the incompressible Navier-Stokes equations,

$$ \left. \begin{array}{l} \nabla \cdot \mathbf{u} = 0, \\ \rho(\partial_t \mathbf{u} + (\mathbf{u} \cdot \nabla)\mathbf{u}) = -\nabla p + \mu \nabla^2 \mathbf{u} + \mathbf{f}, \end{array} \right\} \qquad (2.1) $$

where $\mathbf{u}$ denotes the fluid velocity vector, $p$ is the fluid pressure, $\mu$ is the dynamic viscosity, $\rho$ is the fluid density, and $\mathbf{f}$ is the particle feedback force to the fluid phase in the case of 2WC point-particle simulations, here computed by using a standard particle-in-cell approach, which spreads the particle force to the nearest Eulerian grid points with a linear kernel (Boivin, Simonin & Squires 1998; Lee & Lee 2019; Zhang, Cui & Zheng 2023).

We consider $x$, $y$ and $z$ as the streamwise, wall-normal and spanwise directions, respectively. No-slip and no-penetration boundary conditions are specified at the domain walls, while periodicity is applied along the streamwise and spanwise directions. These equations are discretized in space using a pseudo-spectral approach along $x$ and $z$ and second-order finite differences along $y$. Wray's low-storage third-order Runge-Kutta scheme is employed for time marching (see Wray 1990).

Spherical particles with density $\rho_p \gg \rho$ and particle diameter $D$ in the absence of gravity are tracked with a standard Lagrangian point-particle method, with their dynamics governed by

$$ m_p \dot{U}_p = F_d + F_l, \quad \dot{X}_p = U_p, \qquad (2.2a,b) $$

with the Schiller-Naumann drag force

$$ F_d = -3\pi\mu DU_s(1 + 0.15 Re_p^{0.687}), \qquad (2.3) $$
where $U_p$, $X_p$, $F_d$ and $F_l$ are the particle velocity, position, drag and lift forces, $m_p$ is the particle mass, $U_s = U_p - u|_{x=X_p}$ is the local slip velocity evaluated at the particle location, and $Re_p = |U_s|D/\nu$ is the particle Reynolds number. Since the particle-to-fluid density ratio is high, other dynamic effects such as added mass, fluid acceleration, and the Basset history force are negligible. While in most practical scenarios of particle-laden wall transport, gravity is important at high density ratios, we neglect it in the present work. We do this to isolate the interplay between more intricate wall accumulation mechanisms (e.g. turbophoresis) and lift forces, in a flow that could otherwise feature significant settling effects. Finally, for the sake of simplicity, a perfectly elastic hard-sphere rebound is employed for particle-wall interactions.

### 2.1.1. Lift force models

We consider standard shear-induced lift force models where lift force $F_l$ acts perpendicular to the local shear sampled by the particle and points along $\omega \times U_s$, where $\omega = (\nabla \times u)|_{x=X_p}$ is the undisturbed flow vorticity evaluated at the particle position. For convenience, let us express the lift force in terms of the dimensionless lift coefficient $C_L$, defined as

$$C_L = \frac{F_l}{\frac{1}{8}\pi\rho D^2 |U_s|^2} \cdot \frac{\omega \times U_s}{|\omega \times U_s|} \qquad (2.4)$$

Assuming that the Oseen length $\ell_u = \nu/|U_s|$ is much larger than the Saffman length $\ell_\omega = \sqrt{\nu/|\omega|}$, i.e. $\varepsilon \gg 1$ with $\varepsilon = \sqrt{|\omega|\nu/|U_s|}$, an explicit lift solution can be derived:

$$C_L = \frac{18}{\pi^2} \varepsilon J(\varepsilon), \qquad (2.5)$$

with

$$J(\varepsilon \gg 1) = J^\infty = 2.255. \qquad (2.6)$$

Hereafter, the lift expression (2.5) together with the $J$ function by (2.6) will be referred to as the Saffman model, which is expected to be valid in the double limits $Re_p \to 0$ and $Re_\omega \to 0$ provided that $\varepsilon \gg 1$, where $Re_\omega = |\omega|D^2/\nu$ is the shear Reynolds number.

Still in the double limits $Re_p \to 0$ and $Re_\omega \to 0$, the $J$ function in (2.5) at finite $\varepsilon$ turns out to be a volume integral in Fourier space. Its value cannot be put in closed form, but was estimated numerically in Asmolov (1989), McLaughlin (1991) and more recently, Shi & Rzehak (2019). Based on these numerical data, various empirical correlations of $J(\varepsilon)$ were proposed (see Shi & Rzehak (2019) for a comprehensive review), with the most commonly used being the one proposed by Mei (1992), i.e.

$$J(\varepsilon) = 0.3 J^\infty \left( 1 + \tanh\left[\frac{5}{2}(\log_{10}\varepsilon + 0.191)\right] \right) \left( \frac{2}{3} + \tanh(6\varepsilon - 1.92) \right). \qquad (2.7)$$

Hereafter, the lift expression (2.5) together with the $J$ function by (2.7) will be referred to as the Mei model, which is expected to be valid irrespective of $\varepsilon$, provided that $Re_p$ and $Re_\omega$ are small.

The flow is assumed to be unbounded for the two lift models introduced above. In particle-laden channel flows, near-wall accumulation of particles is often observed, as discussed above. Typically, the peak in the particle concentration appears within the viscous sublayer, i.e. for $y^+ \le 5$, wherein the wall effect on the lift is crucial if $D^+ \ge 1$ (Balachandar & Eaton 2010; Shi et al. 2021). In this context, the explicit lift solution can be derived if the wall lies in the ‘inner region’ of the flow disturbance in the low-$Re_p$ limit,
specifically, if the separation between the particle and the wall, $ℓ$, is much smaller than the inertial lengths ($ℓ \ll \min(ℓ_u, ℓ_ω)$). The corresponding lift solutions (see Shi & Rzehak (2020) for a comprehensive review) take the general form

$$C_L = A + B Sr + C Sr^2, \qquad (2.8)$$

where $Sr = |ω|D/|U_s|$ is the normalized shear rate, and A, B and C are pre-factors that are independent of $Sr$. Hereafter, the terms proportional to $Sr$ and $Sr^2$ will be referred to as the linear and quadratic inner contributions $C_{L,in}^{linear}$ and $C_{L,in}^{quad}$, respectively.

Situations where the wall is in the 'outer region' of the flow disturbance ($ℓ \gg \max(ℓ_u, ℓ_ω)$), while still in the limit $Re_p \to 0$, were considered in Asmolov (1990), McLaughlin (1993) and Takemura, Magnaudet & Dimitrakopoulos (2009). Assuming that $Sr = Re_ω^{1/2} \varepsilon \le O(1)$ (i.e. the same condition for the shear rate as considered in Saffman 1965), the outer-region lift approaches the Saffman solution (2.5), namely $C_{L,out}^{linear} \propto \varepsilon$. Together with the scaling in the inner region that $C_{L,in}^{linear} \propto Sr$, it appears that the inner-region lift contribution transitions into the outer region by a pre-factor $\varepsilon/Sr$. If one assumes that the quadratic contribution follows the same transition, then it may be speculated that $C_{L,in}^{quad}$ scales as $\varepsilon Sr$ in the outer region, i.e. $C_{L,out}^{quad} \propto \varepsilon Sr$. For the particle-laden channel flow,

$$Sr \approx \frac{(u_{\tau}^2/v)D}{u_s} = \left(\frac{u_s}{u_{\tau}}\right)^{-1} D^+, \qquad (2.9)$$

where $u_s$ is the streamwise mean slip velocity, namely $u_s = \langle U_s \cdot e_x \rangle$, with $e_x$ the unit vector along the streamwise direction, and $u_τ$ is the conventional wall friction velocity. As will be demonstrated in figure 5 in § 3.1, $u_s$ might change its sign in the buffer layer, leading to extremely large values of $Sr$ in the near-wall region.

The analysis above implies that there is also a quadratic lift contribution $C_L^{quad} \propto \varepsilon Sr$, which might dominate the lift generation in the inner wall region. This is consistent with the form of the lift force proposed in Costa et al. (2020a,b), where

$$C_L = \frac{18}{\pi^2} \varepsilon Sr J^\infty, \qquad (2.10)$$

with $J^\infty = 2.255$ according to Saffman (1965). Despite its simplified form, this correlation aligns well with their PR-DNS results. Hereafter, the lift model by (2.10) will be referred to as the CBP model.

Finally, it should be noted that in practice, to avoid singularities in the numerical calculation, we implemented a slight variant of the lift coefficient as described in (2.4), where the quotient on the right-hand side to define the unit vector is modified to $ω \times U_s/|ω| |U_s|$. In practice, the results are not very sensitive to this choice, as $ω \times U_s$ is nearly aligned with the wall-normal direction.

## 2.2. Validation

The present code has demonstrated successful applications in wall-bounded particle-laden flows, including open channel flow (Wang & Richter 2019; Gao, Samtaney & Richter 2023) and planar Couette flow (Richter & Sullivan 2014; Richter 2015) loaded with inertial particles, but the lift force is neglected. Hence, for completeness, we present the validation of the 1WC PP-DNS code, including the lift force models. Specifically, we reproduce
<figure><img src="image_4.png" /><figcaption>Figure 1. (a) Normalized particle volume fraction $\phi/\Phi$. Outer-scaled second-order moments of particle velocity: (b) streamwise velocity root mean square, (c) wall-normal velocity root mean square, (d) Reynolds shear stress profile. The points denote the corresponding lift-included 1WC PP-DNS from Costa et al. (2020b). Here, $U_b$ is the flow bulk velocity.</figcaption></figure>

in figure 1 the PP-DNS of particle-laden turbulent channel flow in the 1WC regime using these three lift models, as reported in Costa et al. (2020a,b) in the 1WC regime. This corresponds to case CL1 in table 1. The agreement is excellent, which validates our implementation.

## 2.3. Computational parameters

We perform point-particle DNS of channel flow at $Re_τ = 180$ in a computational domain ($L_x, L_y, L_z$) = (6h, 2h, 3h), where h is the half-channel height. The domain is discretized on ($N_x, N_y, N_z$) = (160, 320, 160) grid points; the grid is slightly stretched to refine the near-wall resolution, corresponding to grid spacing ($\Delta_x^+, \Delta_y^+, \Delta_z^+$) = (6.75, 0.9, 3.375), with the conventional '+' superscript denoting viscous wall scaling. Note that the grid spacing near the wall is comparable to that of the largest particle sizes. All simulations are carried out at a constant time-marching step fixed at $\Delta t^+ = 0.1$ (normalized by viscous unit $v/u_τ^2$), which corresponds to CFL $\approx 0.4$. The DNS results confirm this time step to be sufficiently small that particles and fluid elements could not pass through a grid cell per time step (Zheng, Feng & Wang 2021). The total simulation time is $T \approx 450h/u_τ$ (approximately 80 000 viscous time scales $v/u_τ^2$), which is long enough to guarantee converged statistics (Sardina et al. 2012).

We consider three different set-ups for varying particle diameter, while $St^+ = \tau_p/(v/u_τ^2) = 50$ is fixed by varying the particle density. These parameters are described
<table><caption>Table 1. Physical parameters used in the present point-particle DNS campaign: $N_p$ denotes the total number of particles, and $D^+$ is the inner-scaled particle diameter. The particle density ratio was adjusted to keep the same value of the Stokes number $St^+ = \tau_p/(\nu/u_t^2) = 50$ and bulk solid mass fraction $\Psi_m = 0.337$ for all cases.</caption><thead><tr><th>Case</th><th>D<sup>+</sup></th><th>ρ<sub>p</sub>/ρ</th><th>N<sub>p</sub></th></tr></thead><tbody><tr><td>CPF</td><td></td><td></td><td></td></tr><tr><td>CL1</td><td>3</td><td>100</td><td>5 × 10<sup>4</sup></td></tr><tr><td>CM1</td><td>1</td><td>900</td><td>1.5 × 10<sup>5</sup></td></tr><tr><td>CS1</td><td>0.1</td><td>90 000</td><td>1.5 × 10<sup>6</sup></td></tr></tbody></table>

in table 1. This target particle Stokes number was chosen since it is known to feature highly inhomogeneous particle distributions in wall turbulence; see Sardina et al. (2012). Moreover, it should be noted that a value $D^+ = 3$ corresponds to about one Kolmogorov length scale in the channel bulk, an order of magnitude that is often investigated in the literature using PP-DNS (see e.g. Bernardini 2014; Motoori, Wong & Goto 2022; Zhang et al. 2022). We performed 1WC and 2WC simulations at a fixed mass fraction for different particle sizes and lift force models. Naturally, in the 1WC simulations, the mass fraction is not a governing parameter, and the high number of particles ensures statistical convergence of the results. In the 2WC PP-DNS, the bulk mass fraction is fixed to $\Psi_m \approx 30\%$, to ensure significant turbulence modulation and highlight the effects of lift force model choice in the flow statistics.

# 3. Results and discussion

## 3.1. Particle dynamics

We start by describing the dynamics of the dispersed phase under the effect of lift force models using a simple 1WC approach. This allows us to measure observables sampled by the particle that dictate the validity region of the lift models (e.g. $Sr$, $\epsilon$ and $Re_p$). As we will see, current lift models are bound to fail in the buffer layer.

Figure 2 shows the profiles of the normalized mean particle concentration for all particle sizes. When the lift force is not taken into account, the preferential accumulation of particles in the wall region is significant, which has been observed widely (Picano, Sardina & Casciola 2009; Sardina et al. 2012; Lee & Lee 2015; Gao et al. 2023; Gualtieri et al. 2023), and termed as turbophoresis (Reeks 1983; Johnson, Bassenne & Moin 2020). However, the particles show less wall accumulation when subjected to lift forces, and this effect is more pronounced the larger the particle size. The Saffman model yields a much smaller peak, while the other two models feature smaller deviations, with the highest peak corresponding to the Mei model. When the particles are very small ($D^+ = 0.1$), the differences in the lift force models become almost indistinguishable, with all the models robustly showing the same reduction of wall accumulation. Yet the difference in accumulation in the absence of lift is clear.

To qualitatively describe the lift force effect on the spatial localization of wall accumulation, figure 3 shows instantaneous snapshots of streamwise velocity fluctuations in a wall-parallel plane within the buffer layer ($y^+ = 12$) and in a plane of constant streamwise location ($x^+ = 540$), along with the corresponding particle positions for the case with largest particles, CL1. As expected, the particles show strong spatial localization and inhomogeneous distribution, with larger local density occurring in elongated clusters
Relevance of lift in particulate turbulent wall flows

<figure><img src="image_5.png" /><figcaption>Figure 2. Local particle volume fraction $\phi$, normalized by the bulk value $\Phi$, as a function of the wall-normal distance (1WC): (a) CL1, (b) CM1, (c) CS1. The insets show the normalized particle volume fraction versus the inner-scaled wall distance.</figcaption></figure>

(Sardina et al. 2012). Effects of the lift force model on the gross characteristics of the particle distribution are readily discernible. One notable manifestation is the difference in particle numbers in the wall-parallel (x–z) plane. Consistent with figure 2, the lift-free model yields the highest number of particles, while the Saffman model yields the lowest. Conversely, the Mei and CBP models fall somewhere between these two. This indicates the weakening accumulation of particles near the wall by the lift force, consistent with observations by Marchioli et al. (2007) and Shin et al. (2022). The other manifestation concerns small-scale particle clustering. At first glance, lift-free, Mei and CBP simulations show a more significant tendency to over-sample low-speed regions than the Saffman one. This tendency is checked in figure 4, showing the probability density functions (p.d.f.s) of streamwise fluid velocity fluctuations sampled by the particles within the $y^+ \le 12$ region. For larger particle sizes, particles under the Saffman lift force sample a less focused range of streamwise velocity fluctuations than other models for larger particle sizes. In particular, for $D^+ = 3$ and 1, the CBP and Mei cases show approximately the same tendency to sample low-speed regions, differently than the Saffman model and comparable to the lift-free case for $D^+ = 3$. Preferential accumulation gets more focused with decreasing particle size, with the lift-free case producing the strongest preferential accumulation for $D^+ = 1$ and 0.1.

Figure 5 shows the inner-scaled profiles of streamwise mean particle velocity $\langle u_p \rangle$ and local slip velocity $\langle u_s \rangle$. Compared with the lift-free particles with $D^+ \ge 1$, lift forces increase $\langle u_p \rangle$ in the viscous sublayer and buffer layer, while the differences between cases become negligible beyond $y^+ \sim 30$. Once more, particles under the Saffman lift force show the highest variation with respect to the lift-free case, yielding the highest $\langle u_p \rangle$; particles under the Mei lift force yield the lowest, and the ones under the CBP lift force fall in between. Consistently, for the smallest particle sizes, the overall modification of $\langle u_p \rangle$ by lift forces is quite small. This tendency of lift forces becoming less important with decreasing particle size is expected – as we will discuss later in more detail (§ 3.3), the relative importance of lift to drag forces close to the wall approximately scales with the inner-scaled particle diameter $D^+$.

At the wall, while the fluid velocity must vanish, particles can feature a mean apparent slip where they flow faster than the fluid due to their inertia (Zhao, Marchioli & Andersson 2012). This is shown clearly in the particle mean slip (figures 5d–f). At larger wall distances, the tendency of particles to over-sample slower-than-average fluid velocity (Kiger & Pan 2002; Baker & Coletti 2021) is also clearly reflected in the mean particle
<figure><img src="image_6.png" /><figcaption>Figure 3. Inner-scaled instantaneous streamwise velocity fluctuation and the corresponding instantaneous snapshots of particle locations in the x-z plane ($y^+ = 12$, left-hand column) and y-z plane ($x^+ = 540$, right-hand column) (case CL1, IWC). The wall-parallel plane shows all particles with inner-scaled wall normal position $Y_p^+ \le 12$. (a) Lift-free, (b) Saffman model, (c) Mei model, (d) CBP model.</figcaption></figure>

slip velocity. A similar trend can be observed in the corresponding streamwise mean slip velocity profiles; however, the difference in $\langle u_s \rangle$ is larger than $\langle u_p \rangle$, and this is attributed to the negative apparent slip velocity ($\langle u|_{x=X_p} \rangle - \langle u \rangle < 0$), which reflects preferential sampling of slower-than-average fluid. Interestingly, lift forces increase the apparent slip velocity near the wall. This can be understood from the sign of the lift force in this region, which points towards the wall due to the positive slip velocity. Particles sampling higher-momentum regions will be driven towards the wall by the lift force, and their inertia results in a higher mean slip.

The sign of $\langle u_s \rangle$ changes in approximately the same location of the buffer layer, irrespective of the particle size and lift force model. This mechanism seems relatively
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>Figure 4. Probability density functions (p.d.f.s) of inner-scaled streamwise fluid velocity fluctuations conditioned at the particle location within the $y^+ \le 12$ region (1WC): (a) CL1, (b) CM1, (c) CS1.</figcaption></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /><figcaption>Figure 5. Inner-scaled streamwise mean particle velocity and slip velocity profiles (1WC): (a,d) CL1, (b,e) CM1, (c,f) CS1. The horizontal dashed lines denote $\langle u_s \rangle = 0$.</figcaption></figure>

robust and has been observed in numerous other studies with inertial particles at different Stokes numbers and different flow Reynolds numbers (see e.g. Mortensen et al. 2008; Zhao et al. 2012; Wang & Richter 2020; Gao et al. 2023). It will be shown that this change in sign makes lift force modelling in wall turbulence extremely challenging (see discussions in § 3.3).

Figure 6 shows the x-y component inner-scaled velocity covariances for the particle and fluid phases, respectively. The spanwise velocity variances are hardly affected by lift forces, and thus not shown. As expected, away from the wall where the mean shear is low, the particle statistics closely follow those of the lift-free cases. Near the wall, instead, the lift force enhances the streamwise and wall-normal particle velocity variances ($\langle u_p^2 \rangle$ and $\langle v_p^2 \rangle$) and Reynolds stress ($\langle u'_p v'_p \rangle$), even for the smallest particle size ($D^+ = 0.1$). Enhancement of wall-normal velocity fluctuations near the wall has indeed been reported in previous numerical and experimental studies (Fong, Amili & Coletti 2019; Costa et al. 2020a). Indeed, as also shown in Costa et al. (2020a), this enhancement of velocity fluctuations near the wall cannot be reproduced when lift forces are not considered in
<figure><img src="image_16.png" /><figcaption>Figure 6. Second-order moments of mean particle (solid lines) and fluid (dashed lines) velocity (1WC). (a–c) Inner-scaled streamwise and (d–f) wall-normal velocity variances, and (g–i) Reynolds shear stress profiles, where (a,d,g), (b,e,h) and (c,f,i) denote cases CL1, CM1 and CS1, respectively.</figcaption></figure>

the particle dynamics. It seems that lift force is very effective at correlating streamwise and wall-normal velocity fluctuations near the wall – under high shear rate, small variations in streamwise velocity naturally translate to changes in vertical acceleration through the lift forces, which in turn induce wall-normal velocity fluctuations. Streamwise velocity fluctuations will also be amplified naturally, since streamwise and wall-normal velocity fluctuations are correlated through the mean shear. As we will see in § 3.4, this amplification of near-wall velocity fluctuations by lift forces has major consequences on turbulent modulation and drag changes in the flow at high particle mass loading.

## 3.2. *Mechanism for near-wall accumulation*

To understand the role of the lift force in altering the near-wall particle distributions, it is essential to examine three primary mechanisms responsible for wall-normal inertial particle transport. The first is the lift-induced migration, which is the direct consequence of the inertial lift force. The second is turbophoresis (Reeks 1983), causing particle migration from regions of higher to lower turbulence intensity. In wall-bounded flows, turbophoresis results in strong particle accumulation in the viscous sublayer, where turbulent fluctuations vanish. Note, however, that in most cases, the peak in the particle concentration appears not at the wall but at a distance of *O*(*D*) from the wall (Marchioli & Soldati 2002; Costa et al. 2020a), owing to the wall-particle collisions as well as hydrodynamic wall-particle
Relevance of lift in particulate turbulent wall flows

<figure><img src="image_17.png" /><figcaption>Figure 7. Illustration of the interplay between different mechanisms responsible for the wall-normal particle transport.</figcaption></figure>

interactions (Goldman, Cox & Brenner 1967; Vasseur & Cox 1977; Zeng, Balachandar & Fischer 2005). Finally, there is a biased sampling effect, due to the tendency of particles to sample ejection regions in the buffer layer featuring negative streamwise velocity and positive (i.e. repelling) wall-normal velocity fluctuations (Marchioli & Soldati 2002; Sardina et al. 2012). As a result, biased sampling usually provides a net slow particle drift away from the wall. The interaction among the three mechanisms is sketched in figure 7 and will be elaborated below.

Lift-induced migration is directed along $\omega \times U_s$, which in the case of a turbulent wall flow is predominantly in the wall-normal direction. According to figures 5(d–f), particles lead the liquid flow for $y^+ \lesssim 10$, and lag behind the flow further away from the wall. Hence the lift-induced migration is towards the wall in the viscous sublayer, and acts cooperatively with the turbophoresis effect. The way these two mechanisms interact in the viscous sublayer may seem surprising at first glance, as figure 2 clearly indicates a suppression of the near-wall particle accumulation by lift force. The reasoning behind this suppression may be understood by noting the following two issues. First, the lift force is towards the channel centre within the buffer layer and beyond $y^+ \approx 10$. The resulting outward migration, together with the strong turbulence intensity in the buffer layer, tends to entrain particles in this region outwards, leading to a net decrease in the concentration in the region $y^+ \lesssim 10$ as revealed by figure 2. Second, the lift force influences the wall-normal particle transport by turbophoresis. This phenomenon is linked to the particle wall-normal velocity variance $\langle v_p^2 \rangle$ (Reeks 1983; Johnson et al. 2020). Indeed, in the limit of vanishing particle Reynolds number, the particle wall-normal momentum balance at steady state yields a turbophoresis pseudo-force proportional to $d\langle v_p^2 \rangle/dy$, which creates the migration of particles down gradients in particle wall-normal velocity variance (Sikovsky 2014; Johnson et al. 2020). Recall figures 6(d–f), where the variation of $\langle v_p^2 \rangle$ with the wall distance is depicted. Clearly, the lift force tends to increase the ‘apparent inertia’ of the particles, enhancing the deviation between the particle velocity variance and that of the fluid within the viscous sublayer. This leads to an attenuation in the corresponding turbophoresis pseudo-force. In particular, for the two cases where $D^+ \ge 1$ (figures 6d,e), the prediction with the Saffman lift model yields
$d\langle v_p^2 \rangle/dy < 0$ in the viscous sublayer, indicating that the turbophoresis pseudo-force is outwards, driving the particles away from the wall. Consequently, the lift force competes with the turbophoresis pseudo-force in these two cases, leading to a highly flattened near-wall peak in the particle fraction, as seen in figures 2(a,b). In contrast, they cooperate in predictions using the Mei and the CBP models, as figure 2(a) shows. There, the lift force compensates the attenuation in turbophoresis, making the negative slopes in $\phi(y^+)$ from results using Mei and the CBP models closely follow the lift-free case (highlighted in the figure insets). Note, however, that the magnitude of the Mei lift (CBP lift) is proportional to $Re_\omega^{1/2} (Re_\omega)$, which approaches $D^+(D^{+2})$ in the viscous sublayer. Consequently, for the two cases where $D^+ \le 1$, the lift-induced migration is weak and incapable of compensating the still-pronounced attenuation in turbophoresis revealed by figure 6(e,f). This is why in these two cases the slopes in $\phi(y^+)$ from results employing the Mei and CBP models are smaller in magnitude than in the lift-free case. Still, the difference in accumulation with respect to the lift-free case for smaller sizes is clear, due to near-wall particles residing for long times near the wall while being subjected to small but persistent lift forces.

Given the discussion above, it is worth examining in more detail the link between the chosen lift model and the corresponding modulation in the turbophoresis pseudo-force, particularly to what extent this modulation relates to the lift force magnitude. Combining the turbophoresis pseudo-force (figures 6d–f) and lift force as shown in figure 8, a monotonic increase in the suppression of the turbophoresis pseudo-force with increasing lift is observed at $D^+ \le 1$, while no such trend is evident at $D^+ = 3$. In the latter case, the suppression of the turbophoresis pseudo-force is most pronounced in the prediction using the Saffman model, even though the magnitude of the corresponding lift force is only half that of the prediction using the CBP model for $y^+ \lesssim 5$. Hence the differences in peak concentration are directly linked not to the lift force magnitudes, but rather to the modulation in the turbophoresis pseudo-force induced by these models. In the present work, the turbophoresis pseudo-force is always suppressed by the presence of lift, leading to a decrease in the peak concentration. Hence the precise relationship between the extent of this suppression and the lift force magnitude remains unclear, making it challenging to directly connect the chosen lift models with the changes in the predicted concentration profiles. Finally, it would be worth investigating these dynamics at lower Stokes numbers, as both turbophoresis and preferential sweeping will be less pronounced, and the nature of the interaction between lift and turbulence may change drastically. Given the increased issues with lift force modelling at lower Stokes number, discussed in the next subsection, addressing this question should be the object of a future, dedicated study.

To gain insights into biased sampling of sweeps/ejection events under lift forces, we carried out a quadrant analysis for the fluid flow experienced by the particles. In the $u'-v'$ plane, with positive $v'$ directed away from the wall, the flow experienced by the particle is categorized into four types of events: first quadrant events (Q1), characterized by outward motion of high-speed fluid, with $u' > 0$ and $v' > 0$; second quadrant events (Q2), characterized by outward motion of low-speed fluid, with $u' < 0$ and $v' > 0$, which are usually called ejections; third quadrant events (Q3), characterized by inward motion of low-speed fluid, with $u' < 0$ and $v' < 0$; and finally, fourth quadrant events (Q4), which represent motions of high-speed fluid towards the wall, with $u' > 0$ and $v' < 0$, and are usually called sweeps. Figure 9 shows the joint p.d.f.s in the $u'-v'$ plane of the fluid fluctuation seen by ascending ($v_p > 0$) and descending ($v_p < 0$) particles within the $y^+ \le 12$ region. The first observation is that a larger proportion of particles tends
Relevance of lift in particulate turbulent wall flows

<figure><img src="image_18.png" /><figcaption>Figure 8. The inner-scaled lift-induced acceleration (i.e. normalized by $u_t^3/\nu$) profiles along the wall-normal direction (1WC): (a) CL1, (b) CM1, (c) CS1. The horizontal dashed lines denote $\langle F_{l,y} \rangle = 0$.</figcaption></figure>

to move outwards (ascending), irrespective of the presence of lift force. The proportion of ascending particles without lift force decreases as the particle size decreases, whereas particles subject to lift force exhibit an inverse trend. Second, the lift-free ascending particles are highly likely to sample ejections (Q2), while the lift-free descending ones are prone to sample sweeps (Q4) for particles with $D^+ \ge 1$. This agrees with experimental observations in channel flow laden with finite-size particles (Kiger & Pan 2002; Yu, Vinkovic & Buffat 2016; Baker & Coletti 2021). However, the particles under the effect of lift forces tend not to over-sample sweeps, irrespective of the wall-normal particle velocity sign, which indicates that the wall-normal velocity of the particles is less correlated with the fluid one. This also confirms that inertial particles under the effect of lift forces are not prone to slowly drive towards the wall in low-speed regions, because the competing wall-repelling effect of lift forces seems to dominate their dynamics. In close inspection of the high-intensity region of the joint p.d.f.s on the $u'$-axis, we note that particles under the Saffman lift force tend to sample a wider range of negative $u'$ regions, compared with particles under the Mei and CBP lift forces. This aligns with the p.d.f.s of $u'$ as shown in figure 4. Third, irrespective of the presence of the lift force or wall-normal particle velocity sign, the smallest particles ($D^+ = 0.1$) over-sample ejection (Q2) and inward motions of low-speed fluid regions (Q3). Sweep events, bringing high-momentum fluid towards the wall, are experienced only weakly by particles. Overall, the percentage of particles sampling Q2 events only slightly surpasses that for Q3, suggesting that biased sampling may play a minor role in driving high-inertia particles ($St^+ = 50$) to depart from the inner wall region. This observation aligns with numerical results from Marchioli & Soldati (2002), where biased sampling is found to be more efficient for transferring particles with smaller inertia.

### 3.3. Emerging challenges in lift force modelling

The aforementioned results allow us to address several issues concerning the effects of lift force in turbulent wall flows, which may have been overlooked previously in PP-DNS. Ideally, such discussions would be based on results from PR-DNS. However, in this instance, we will extrapolate from PP-DNS results, in conjunction with the CBP model, as an analogue to PR-DNS outcomes. This approach is justified by findings from Costa et al. (2020a,b), which demonstrate that PP-DNS are capable of satisfactorily reproducing the corresponding PR-DNS results for the most challenging case considered in this study, with $D^+ = 3$.
<figure><img src="image_19.png" /><figcaption>Figure 9. Joint p.d.f.s of the streamwise and wall-normal fluid fluctuating velocities seen by ascending ($v_p > 0$) and descending ($v_p < 0$) particles within the $y^+ \le 12$ region (1WC): (a) CL1, (b) CM1, (c) CS1. The numbers in the images denote the proportion of ascending particles.</figcaption></figure>
Let us first address the circumstances under which the effects of lift force could be considered negligible for large inertial particles, and thereby disregarded in PP-DNS. Given that biased sampling does not play a key role in the wall-normal transport of inertial particles, a criterion based on the ratio of lift to the turbophoresis pseudo-force seems sufficient to describe the relative significance of the lift force. In this context, the dimensionless particle size $D^+$ seems a promising candidate, for the following two reasons. First, $D^{+2}$ approaches the shear Reynolds number $Re_\omega$ in the viscous sublayer; the latter directly measures the magnitude of the shear lift irrespective of the choice of lift models (Saffman 1965; McLaughlin 1991; Mei 1992; Shi & Rzehak 2019; Costa et al. 2020b). Second, the turbophoresis effects seem to decay with increasing $D^+$. This is clearly revealed by comparing the slopes of $\langle v_p^2 \rangle(y^+)$ in figures 6(d-f), which correspond to results with increasing $D^+$. Gluing these two aspects together, it seems that the force ratio of lift to turbophoresis approximately scales as $D^{+k}$ with $k > 0$. Although it is not feasible to determine the exact value of $k$, the results summarized in § 3.1 made it clear that by no means should the lift force be neglected for $D^+ \ge 1$. On the other hand, it is reasonable for the lift force to be neglected for $D^+ \le 0.1$, except for the under-prediction of near-wall accumulation discussed above (recall figure 2).

We should now address a major issue in commonly applied lift force models for PP-DNS of turbulent wall flows. As outlined in § 2.1, the shear lift models originating from the pioneering work of Saffman (1965) are applicable rigorously only in the double limits $Re_p \to 0$ and $Re_\omega \to 0$. In the viscous sublayer, $Re_p \approx D^+(\langle u_s \rangle/u_\tau) \approx D^+$ for $D^+ \ge 1$ according to figure 5, whereas $Re_\omega \approx D^{+2}$. Consequently, neither of the two limits is satisfied for $D^+ \ge 1$. In addition to the double limits above, there is a constraint concerning the dimensionless shear rate $Sr$. This is made clear by noting that $Sr = Re_\omega^{1/2} \epsilon$ and hence $Sr \le O(1)$ even in the strong shear limit $\epsilon \ge 1$ considered in Saffman (1965). This is why in the wall-bounded situation, the ‘outer-region’ lift solutions (Asmolov 1990; McLaughlin 1993; Takemura et al. 2009) following the methodology of Saffman (1965) match the corresponding ‘inner-region’ solutions (Cherukat & McLaughlin 1994; Magnaudet, Takagi & Legendre 2003; Shi & Rzehak 2020) only for $Sr = O(1)$. For particle-laden channel flows with small inertial particles (Mortensen et al. 2008; Zhao et al. 2012; Wang & Richter 2020; Gao et al. 2023), this constraint is more seriously violated than that for the double limits above, as the slip velocity reverses at $y^+ \approx O(10)$ (see e.g. figures 5(d-f) in the present work), making $Sr \to \infty$ for particles within the buffer layer. This issue is highlighted in figure 10, where we show the p.d.f.s of the inverse of $Sr$ (horizontal axis) as functions of the wall-normal distance $y^+$ (vertical axis). Apparently, a large portion of particles experience a strong shear where $Sr \ge 1$, in particular for the two cases where $D^+ \ge 1$. Finally, it is worth noting that the probability of a particle experiencing very high values of $Sr$ near the wall will increase with decreasing values of $St^+$ as, for the same local shear rate, particles will experience an ever smaller slip velocity. This exacerbates the violation of the $Sr = O(1)$ constraint, meaning that current lift force models are even less suitable for turbulent channel flow laden with particles at lower Stokes numbers.

Additionally, the discussion above focuses on the lift force arising from the shear in the ambient flow. In wall-bounded flows, as considered in this work, the lift force may deviate significantly from its unbounded counterpart if the particle is close to the wall. In brief, the presence of the wall leads to a repulsive transverse force in the absence of the ambient shear (Vasseur & Cox 1977; Zeng et al. 2005), and in the presence of an ambient shear, a
<figure><img src="image_20.png" /><figcaption>Figure 10. The p.d.f.s of $u_s/(\omega_z D)$ along the wall-normal direction (1WC DNS with the CBP model): (a) CL1, (b) CM1, (c) CS1. The vertical dashed lines denote $u_s/(\omega_z D) = \pm 1$.</figcaption></figure>

suppression of the shear-induced lift (McLaughlin 1993; Magnaudet et al. 2003; Takemura et al. 2009). Although analytical solutions for the combined lift force are achievable in some asymptotic limits (see e.g. Wang et al. (1997) and, more recently, Shi & Rzehak (2020) for an overview of these solutions), there is so far no satisfactory way to glue these solutions together to achieve a general treatment of the wall effect (Ekanayake et al. 2020; Ekanayake, Berry & Harvie 2021). Here, instead of trying to make any further progress in achieving such a general solution, it is more feasible to discuss under which condition this intricate wall effect may be neglected. Within the $O(D)$-depth layer close to the wall, wall-particle collisions govern the wall-normal position of the peak particle concentration. Hence the hydrodynamic wall-particle interaction may play a significant role only beyond $y^+ \approx D^+$. On the other hand, it is known that (McLaughlin 1993) the amount of reduction of the shear-induced lift due to the wall decreases as $(\ell/\ell_\omega)^{-5/3}$, and reduces to only 20% of its magnitude in the limit $\ell/\ell_\omega \to 0$ at $\ell/\ell_\omega \approx 3$. Hence the wall effect may be considered negligible if the position of the peak particle concentration, which is of $O(D)$ from the wall, is larger than approximately $3\ell_\omega$, which yields approximately $D^+ \ge 3$. This finding, together with the previous one where (shear-induced) lift force is non-negligible for $D^+ \ge 1$, indicates that such a wall effect can be disregarded except in the overlap regime where $1 \le D^+ \le 3$.

In addition to the ambient shear, particle rotation is known to produce a lift force as well (Rubinow & Keller 1961). This lift contribution is often neglected in prior PP-DNS and not considered in the present work. Nevertheless, it might be worth discussing its relevant importance in the following. Assuming that the particle inertia is small enough for the particle to always stay in the torque-free state, particle rotation leads to a lift contribution $C_L^\Omega = 0.5 \text{ Sr}$ in the double limits considered by Saffman (1965). The ratio between this spin-induced lift contribution and that induced by ambient shear is proportional to $Re_\omega^{1/2}$, indicating that $C_L^\Omega$ is merely a component of the second-order lift contribution from the ambient shear. This is confirmed by the recent work of Candelier et al. (2023), where second-order lift contributions (with respect to $Re_\omega^{1/2}$) were obtained using matched asymptotic expansions. In addition to the spin-induced contribution $C_L^\Omega = 0.5 \text{ Sr}$, they also obtained a second-order contribution from the ambient shear, $C_L^{(2)} = -0.505 \text{ Sr}$, which counterbalances the spin-induced contribution. In other words, the effect of particle rotation on the lift is negligible in the torque-free state. This conclusion is valid irrespective
Relevance of lift in particulate turbulent wall flows

of the presence of the wall (Cherukat & McLaughlin 1994; Magnaudet et al. 2003; Shi & Rzehak 2020). However, large-inertia particles with a large density ratio as considered in this work might not always stay torque-free. Specifically, for particles with density ratio $ρ_p/ρ \ge 100$, the relaxation time in response to torque is comparable with that to slip (see e.g. Bagchi & Balachandar 2002). This finding, together with the pronounced slip in the viscous sublayer revealed by figures 5(a,b), indicates that particles with $D^+ \ge 1$ may experience a strong angular acceleration and possibly a non-negligible lift contribution from forced rotation. Investigating the effects of rotation by accounting for the particles’ conservation of angular momentum would therefore be interesting but is outside the scope of the present work due to the higher complexity associated with drag and lift closures for angular momentum.

Finally, it is worth noting that while the CBP lift force model used in Costa et al. (2020a,b) appears to provide reasonable predictions in some cases, it is not valid universally. Specifically, it does not align with any of the inner-region solutions (Shi & Rzehak 2020), which are applicable for particles extremely close to the wall (e.g. within the viscous sublayer), nor does it approach the unbounded solution (Saffman 1965; McLaughlin 1991) for particles located within the outer layer of the channel flow. The agreement between the PP-DNS utilizing this lift force model and PR-DNS suggests that the model offers a satisfactory blending of results in the two limits for the specific cases explored in Costa et al. (2020a) (see table 1). A preliminary attempt to establish a more general lift force model can be found in the work of Ekanayake et al. (2020, 2021), where a lift model satisfying the above conditions was proposed in the limit of small shear Reynolds numbers, i.e. $Re_ω \ll 1$. However, the applicability of this model is limited to $D^+ \ll 1$ because in the viscous layer, $Re_ω \approx (u_τ^2/ν)D^2/ν = (D^+)^2$. While this provides some insights, its practicality may be limited since the lift force no longer plays a pivotal role for $D^+ < 1$, as demonstrated here. Hence further work extending Ekanayake et al. (2020, 2021) to finite $Re_ω$ (hence $D^+$) would be significant.

3.4. Lift-induced turbulence modulation

Finally in this section, we use 2WC point-particle DNS to illustrate that lift force models can have tremendous consequences in basic integral quantities such as the overall drag. To achieve this, we use the same parameters as in table 1, which fixes the total mass fraction and particle Stokes number while varying the particle diameter. We should reiterate that high-fidelity 2WC point-particle DNS are highly challenging and are currently being investigated actively. Still, the simple particle-in-cell 2WC DNS employed here suffice to illustrate the impact of lift force model choice on important turbulence metrics. In what follows, we illustrate the effects of lift force on the overall drag, and analyse the streamwise momentum fluxes that contribute to it.

Figure 11 shows the friction coefficient for the 2WC DNS, defined as $C_f = 2u_τ^2/U_b^2$. While all small-inertia particles have a drag-increasing effect, the difference between lift force models is highly amplified as the particle size increases. At $D^+ = 3$, in particular, particles under the Saffman lift force show a drag reduction compared to the other particle sizes, almost reaching an overall drag decrease. Instead, the other lift models show a monotonic drag increase.

To better understand and quantify this difference, we follow the analysis in Costa et al. (2021) and investigate the stress budget in the 2WC limit of vanishing volume fraction,
<figure><img src="image_21.png" /><figcaption>Figure 11. The wall friction coefficient $C_f$ versus the particle diameter $D^+$ (2WC). The horizontal dashed line denotes the mean $C_f$ of the particle-free case.</figcaption></figure>

while keeping the mass fraction finite:

$$ \tau = \rho u_{\tau}^{2} \left(1 - \frac{y}{h}\right) \approx \underbrace{\left(\mu \frac{d\langle u \rangle}{dy} - \rho \langle u'v' \rangle\right)}_{\tau_{\nu} + \tau_{T,f}} - \underbrace{\langle \psi \rangle \langle u'_{p}v'_{p} \rangle}_{\tau_{T,p}}, \quad (3.1) $$

where $\langle \psi \rangle = \rho_p \langle \phi \rangle$ is the mean local mass fraction. The terms on the right-hand-side of this equation are the viscous stress, fluid Reynolds stress and particle Reynolds stress (hereafter $\tau_\nu$, $\tau_{T,f}$ and $\tau_{T,p}$). This equation boils down to the stress budget in the single-phase limit ($\tau_\nu$ and $\tau_{T,f}$), corrected by an additional turbulent particle momentum flux $\tau_{T,p}$ that involves the local mass fraction and the particle counterpart of the Reynolds shear stress. This term relates to the transfer of streamwise momentum by correlated particle velocity fluctuations, and shows the importance of high particle mass fraction near the wall when combined with correlated streamwise-wall-normal particle velocity fluctuations. We have seen in § 3.1 that while the former decreases in the presence of lift force, the latter may be increased by it.

The profiles of stress terms and the corresponding contributions of each momentum transfer mechanism to the overall friction coefficient, normalized by the overall friction coefficient of the particle-free case, are shown in figure 12. In the left-hand plots, the normalized shear stress profiles are shown, while on the right-hand side we show the so-called FIK (Fukagata-Iwamoto-Kasagi) identity employed to a particle suspension (Fukagata, Iwamoto & Kasagi 2002; Yu et al. 2021). In the case of a turbulent plane channel flow, this identity decomposes the mean friction coefficient ($C_f$) into laminar and turbulent contributions. Here, the latter can be further decomposed into a contribution due to the fluid phase and one due to the particles. Each contribution is given by the weighted integral $C_{f,i} = (1/h^2) \int_0^h 6(h-y)\tau_i dy/(\rho U_b^2)$, with $\tau_i$ being one of the stress mechanisms described above. Note that the residual of this stress budget is virtually zero, which validates the approximation in (3.1).

As anticipated from figure 11, all 2WC DNS show increased drag compared to the unladen flow. Yet, remarkably, different lift models show qualitatively different turbulence modulation for the largest particle size ($D^+ = 3$). Particles under the Saffman lift force
Relevance of lift in particulate turbulent wall flows

<figure><img src="image_22.png" /><figcaption>Figure 12. (a,c,e) Budget of streamwise momentum and (b,d,f) contribution of each stress contribution $\tau_i$ to the mean wall friction coefficient $C_f = 2\tau_w/(\rho U_b^2)$, i.e. $(1/\hbar^2) \int_0^h 6(h-y)\tau_i dy/(\rho U_b^2)$ normalized by that of the particle-free case: (a,b) CL1, (c,d) CM1, (e,f) CS1. Different line colours in (a,c,e) denote different cases, which match the legend in (c).</figcaption></figure>

have a major net turbulence attenuation effect, which results in an almost drag-reducing flow with respect to the unladen case. Conversely, particles under the Mei and CBP lift forces have a net drag increase, caused mostly by an increase in the 2WC stress term $\psi \langle u'_p v'_p \rangle$. This results from the amplification of particle Reynolds shear stresses $\langle u'_p v'_p \rangle$ discussed in § 3.1, which compensate for the decrease in local mass fraction $\psi$. Naturally, as the particle size decreases, the differences between lift models become less pronounced but still noticeable at all particle sizes, with particles under the Saffman lift force consistently showing the highest fluid turbulence attenuation. For the smallest
particle size $D^+ = 0.1$, the two turbulent stress (fluid and particle) terms still show noticeable but in practice negligible differences (recall figure 11).

Hence it appears that the most important dynamical effects of lift forces in 2WC conditions can be seen as finite-size effects that directly modify single-point moments of particle velocity for $D^+ \gtrsim 1$ (recall figure 6). However, the lift force still affects near-wall turbulence modulation when particles are relatively small, as they strongly affect the mean particle concentration (recall figure 2), directly linked to turbulence modulation via the 2WC stress term.

4. Conclusions

We have investigated the importance of lift forces in the dynamics of particle-laden turbulent channel flow. Three relevant shear-induced lift force models are analysed using PP-DNS for varying particle size and fixed Stokes number. Our results confirm that shear-induced lift forces strongly influence particle dynamics. We have analysed the main mechanisms for near-wall particle accumulation under lift forces: lift-induced migration, turbophoresis and biased sampling. Very close to the wall, where the particle-to-fluid slip velocity is positive and particle fluctuations are vanishing, lift forces and turbophoresis cooperate in driving particles towards the wall. Particle-wall interactions and biased sampling balance this effect. Further away from the wall, lift forces and turbophoresis compete, with the former having a wall-repelling effect. Concerning biased sampling, we observe that for larger particle sizes ($D^+ \ge 1$), the lift-free ascending particles prefer to sample ejections, while the descending ones tend to sample sweeps. Conversely, particles under lift forces tend to sample ejection regions while experiencing a marginal occurrence of sweep events, regardless of the employed lift model or wall-normal particle velocity sign. This is again a consequence of the competing wall-repelling effect of lift force. For the smallest particles ($D^+ = 0.1$), all cases exhibit over-sampling in ejection and inward motions of low-speed fluid regions (i.e. Q3 events). Overall, the sampling percentage of the ejections events only slightly surpasses that of Q3, suggesting that biased sampling plays a minor direct role in near-wall particle transport. However, the biased sampling of low-speed regions exposes particles for long times to persistent lift forces, which results in significantly reduced particle accumulation, even for small particle sizes.

These different dynamics have major consequences for turbulence modulation due to a modification of the near-wall values of (1) the local particle concentration, and (2) the correlated streamwise-wall-normal particle velocity fluctuations. While the latter is modified only for sufficiently large particle size and therefore can be seen as a finite-size effect, the former is still non-negligible for relatively small particle sizes, and can still result in small but visible changes in overall drag. Indeed, correlated particle velocity fluctuations show little sensitivity to lift forces for sufficiently small values of inner-scaled particle diameter ($D^+ \lesssim 0.1$). Conversely, for small particle sizes, the modification of near-wall particle accumulation due to lift showed little sensitivity to the flavour of shear-induced lift model. Hence to account accurately for lift-induced turbulence modulation with $D^+ = O(1)$ particle sizes, reliable shear-induced lift models are needed that correctly reproduce the particle turbulent momentum flux. This could be seen as a necessary correction for finite-size effects, in the same spirit as the Faxén correction in the Maxey–Riley–Gatignol equation.

While such correction is important for reliable particle-modelled simulations of turbulent wall flows, our analysis has shown that current lift force models are bound to lose their validity in the near-wall region for small inertial particles. The reason is that
Relevance of lift in particulate turbulent wall flows

the well-known tendency of particles to over-sample low-speed regions in the buffer layer, along with their tendency to flow faster than the fluid in the viscous sublayer, results in a region near the wall of vanishing particle-to-fluid slip velocity. Naturally, this issue will be even more pronounced at lower Stokes numbers, where the range of very low slip velocities in regions of high shear is wider. Unfortunately, current shear-induced lift force models are available only for non-dimensional shear rates $Sr$ of order unity, which robustly tends to infinity within the buffer layer. Hence the condition $Sr \le 1$ required by these models is seriously violated in wall turbulence with small inertial particles.

**Acknowledgements.** The Cray XC40 Shaheen II at KAUST was used for all simulations reported. P.C. thanks F. Picano from Padova University for insightful discussions.

**Funding.** This research was partially supported by the KAUST Office of Sponsored Research (OSR) under award no. OSR-2019-CCF-3666 and under baseline research funds of M.P.

**Declaration of interests.** The authors report no conflict of interest.

Author ORCIDs.

© Wei Gao https://orcid.org/0000-0001-7313-0058;

© Pengyu Shi https://orcid.org/0000-0001-6402-4720;

© Matteo Parsani https://orcid.org/0000-0001-7300-1280;

© Pedro Costa https://orcid.org/0000-0001-7010-1040.

REFERENCES

ARCEN, B., TANIÈRE, A. & OESTERLÉ, B. 2006 On the influence of near-wall forces in particle-laden channel flows. *Intl J. Multiphase Flow* **32**, 1326–1339.

ASMOLOV, E.S. 1989 Lift force exerted on a spherical particle in a laminar boundary layer. *Fluid Dyn.* **24** (5), 710–714.

ASMOLOV, E.S. 1990 Dynamics of a spherical particle in a laminar boundary layer. *Fluid Dyn.* **25** (6), 886–890.

BAGCHI, P. & BALACHANDAR, S. 2002 Effect of free rotation on the motion of a solid sphere in linear shear flow at moderate Re. *Phys. Fluids* **14** (8), 2719–2737.

BAKER, L.J. & COLETTI, F. 2021 Particle–fluid–wall interaction of inertial spherical particles in a turbulent boundary layer. *J. Fluid Mech.* **908**, A39.

BALACHANDAR, S. & EATON, J.K. 2010 Turbulent dispersed multiphase flow. *Annu. Rev. Fluid Mech.* **42**, 111–133.

BATTISTA, F., GUALTIERI, P., MOLLICONE, J.-P., SALVADORE, F. & CASCIOLA, C.M. 2023 Drag increase and turbulence augmentation in two-way coupled particle-laden wall-bounded flows. *Phys. Fluids* **35**, 045133.

BERNARDINI, M. 2014 Reynolds number scaling of inertial particle statistics in turbulent channel flows. *J. Fluid Mech.* **758**, R1.

BERNARDINI, M., PIROZZOLI, S. & ORLANDI, P. 2013 The effect of large-scale turbulent structures on particle dispersion in wall-bounded flows. *Intl J. Multiphase Flow* **51**, 55–64.

BOIVIN, M., SIMONIN, O. & SQUIRES, K.D. 1998 Direct numerical simulation of turbulence modulation by particles in isotropic turbulence. *J. Fluid Mech.* **375**, 235–263.

BOTTO, L., NARAYANAN, C., FULGOSI, M. & LAKEHAL, D. 2005 Effect of near-wall turbulence enhancement on the mechanisms of particle deposition. *Intl J. Multiphase Flow* **31**, 940–956.

BRANDT, L. & COLETTI, F. 2022 Particle-laden turbulence: progress and perspectives. *Annu. Rev. Fluid Mech.* **54**, 159–189.

CANDELIER, F., MEHADDI, R., MEHLIG, B. & MAGNAUDET, J. 2023 Second-order inertial forces and torques on a sphere in a viscous steady linear flow. *J. Fluid Mech.* **954**, A25.

CAPECELATRO, J., DESJARDINS, O. & FOX, R.O. 2018 On the transition between turbulence regimes in particle-laden channel flows. *J. Fluid Mech.* **845**, 499–519.

CHERUKAT, P. & MCLAUGHLIN, J.B. 1994 The inertial lift on a rigid sphere in a linear shear flow field near a flat wall. *J. Fluid Mech.* **263**, 1–18.
COSTA, P., BRANDT, L. & PICANO, F. 2020a Interface-resolved simulations of small inertial particles in turbulent channel flow. *J. Fluid Mech.* **883**, A54.

COSTA, P., BRANDT, L. & PICANO, F. 2020b Interface-resolved simulations of small inertial particles in turbulent channel flow – Corrigendum. *J. Fluid Mech.* **891**, E2.

COSTA, P., BRANDT, L. & PICANO, F. 2021 Near-wall turbulence modulation by small inertial particles. *J. Fluid Mech.* **922**, A9.

CROWE, C.T., SHARMA, M.P. & STOCK, D.E. 1977 The particle-source-in cell (PSI-CELL) model for gas-droplet flows. *J. Fluids Engng* **99** (2), 325–332.

EATON, J.K. & FESSLER, J.R. 1994 Preferential concentration of particles by turbulence. *Intl J. Multiphase Flow* **20**, 169–209.

EKANAYAKE, N.I.K., BERRY, J.D. & HARVIE, D.J.E. 2021 Lift and drag forces acting on a particle moving in the presence of slip and shear near a wall. *J. Fluid Mech.* **915**, A103.

EKANAYAKE, N.I.K., BERRY, J.D., STICKLAND, A.D., DUNSTAN, D.E., MUIR, I.L., DOWER, S.K. & HARVIE, D.J.E. 2020 Lift and drag forces acting on a particle moving with zero slip in a linear shear flow near a wall. *J. Fluid Mech.* **904**, A6.

ELGHOBASHI, S. 1994 On predicting particle-laden turbulent flows. *Appl. Sci. Res.* **52**, 309–329.

FESSLER, J.R., KULICK, J.D. & EATON, J.K. 1994 Preferential concentration of heavy particles in a turbulent channel flow. *Phys. Fluids* **6** (11), 3742–3749.

FONG, K.O., AMILI, O. & COLETTI, F. 2019 Velocity and spatial distribution of inertial particles in a turbulent channel flow. *J. Fluid Mech.* **872**, 367–406.

FRÖHLICH, K., SCHNEIDERS, L., MEINKE, M. & SCHRÖDER, W. 2018 Validation of Lagrangian two-way coupled point-particle models in large-eddy simulations. *Flow Turbul. Combust.* **101** (2), 317–341.

FUKAGATA, K., IWAMOTO, K. & KASAGI, N. 2002 Contribution of Reynolds stress distribution to the skin friction in wall-bounded flows. *Phys. Fluids* **14** (11), L73–L76.

GAO, W., SAMTANEY, R. & RICHTER, D.H. 2023 Direct numerical simulation of particle-laden flow in an open channel at $Re_τ = 5186$. *J. Fluid Mech.* **957**, A3.

GATIGNOL, R. 1983 The Faxén formulas for a rigid particle in an unsteady non-uniform Stokes-flow. *J. Méc. Théor. Appl.* **2** (2), 143–160.

GOLDMAN, A.J., COX, R.G. & BRENNER, H. 1967 Slow viscous motion of a sphere parallel to a plane wall – I. Motion through a quiescent fluid. *Chem. Engng Sci.* **22** (4), 637–651.

GUALTIERI, P., BATTISTA, F., SALVADORE, F. & CASCIOLA, C.M. 2023 Effect of Stokes number and particle-to-fluid density ratio on turbulence modification in channel flows. *J. Fluid Mech.* **974**, A26.

GUALTIERI, P., PICANO, F., SARDINA, G. & CASCIOLA, C.M. 2015 Exact regularized point particle method for multiphase flows in the two-way coupling regime. *J. Fluid Mech.* **773**, 520–561.

HORNE, W.J. & MAHESH, K. 2019 A massively-parallel, unstructured overset method to simulate moving bodies in turbulent flows. *J. Comput. Phys.* **397**, 108790.

HORWITZ, J.A.K., IACCARINO, G., EATON, J.K. & MANI, A. 2022 The discrete Green’s function paradigm for two-way coupled Euler–Lagrange simulation. *J. Fluid Mech.* **931**, A3.

HORWITZ, J.A.K. & MANI, A. 2016 Accurate calculation of Stokes drag for point-particle tracking in two-way coupled flows. *J. Comput. Phys.* **318**, 85–109.

HORWITZ, J.A.K. & MANI, A. 2020 Two-way coupled particle–turbulence interaction: effect of numerics and resolution on fluid and particle statistics. *Phys. Rev. Fluids* **5**, 104302.

IRELAND, P.J. & DESJARDINS, O. 2017 Improving particle drag predictions in Euler–Lagrange simulations with two-way coupling. *J. Comput. Phys.* **338**, 405–430.

JIE, Y., CUI, Z., XU, C. & ZHAO, L. 2022 On the existence and formation of multi-scale particle streaks in turbulent channel flows. *J. Fluid Mech.* **935**, A18.

JOHNSON, P.L., BASSENNE, M. & MOIN, P. 2020 Turbophoresis of small inertial particles: theoretical considerations and application to wall-modelled large-eddy simulations. *J. Fluid Mech.* **883**, A27.

KAFTORI, D., HETSRONI, G. & BANERJEE, S. 1995 Particle behavior in the turbulent boundary layer. I. Motion, deposition, and entrainment. *Phys. Fluids* **7** (5), 1095–1106.

KIGER, K.T. & PAN, C. 2002 Suspension and turbulence modification effects of solid particulates on a horizontal turbulent channel flow. *J. Turbul.* **3**, N19.

KUERTEN, J.G.M. 2006 Subgrid modeling in particle-laden channel flow. *Phys. Fluids* **18** (2), 025108.

KUERTEN, J.G.M. 2016 Point-particle DNS and LES of particle-laden turbulent flow – a state-of-the-art review. *Flow Turbul. Combust.* **97**, 689–713.

LEE, J. & LEE, C. 2015 Modification of particle-laden near-wall turbulence: effect of Stokes number. *Phys. Fluids* **27**, 023303.

LEE, J. & LEE, C. 2019 The effect of wall-normal gravity on particle-laden near-wall turbulence. *J. Fluid Mech.* **873**, 475–507.
Relevance of lift in particulate turbulent wall flows

MAGNAUDET, J., TAKAGI, S. & LEGENDRE, D. 2003 Drag, deformation and lateral migration of a buoyant drop moving near a wall. *J. Fluid Mech.* **476**, 115–157.

MARCHIOLI, C., GIUSTI, A., SALVETTI, M.V. & SOLDATI, A. 2003 Direct numerical simulation of particle wall transfer and deposition in upward turbulent pipe flow. *Intl J. Multiphase Flow* **29** (6), 1017–1038.

MARCHIOLI, C., PICCIOTTO, M. & SOLDATI, A. 2007 Influence of gravity and lift on particle velocity statistics and transfer rates in turbulent vertical channel flow. *Intl J. Multiphase Flow* **33**, 227–251.

MARCHIOLI, C. & SOLDATI, A. 2002 Mechanisms for particle transfer and segregation in a turbulent boundary layer. *J. Fluid Mech.* **468**, 283–315.

MARCHIOLI, C., SOLDATI, A., KUERTEN, J.G.M., ARCEN, B., TANIERE, A., GOLDENSOPH, G., SQUIRES, K.D., CARGNELUTTI, M.F. & PORTELA, L.M. 2008 Statistics of particle dispersion in direct numerical simulations of wall-bounded turbulence: results of an international collaborative benchmark test. *Intl J. Multiphase Flow* **34** (9), 879–893.

MAXEY, M. 2017 Simulation methods for particulate flows and concentrated suspensions. *Annu. Rev. Fluid Mech.* **49**, 171–193.

MAXEY, M.R. & RILEY, J.J. 1983 Equation of motion for a small rigid sphere in a nonuniform flow. *Phys. Fluids* **26** (4), 883–889.

MCLAUGHLIN, J.B. 1989 Aerosol particle deposition in numerically simulated channel flow. *Phys. Fluids A* **1** (7), 1211–1224.

MCLAUGHLIN, J.B. 1991 Inertial migration of a small sphere in linear shear flows. *J. Fluid Mech.* **224**, 261–274.

MCLAUGHLIN, J.B. 1993 The lift on a small sphere in wall-bounded linear shear flows. *J. Fluid Mech.* **246**, 249–265.

MEHRABADI, M., HORWITZ, J.A.K., SUBRAMANIAM, S. & MANI, A. 2018 A direct comparison of particle-resolved and point-particle methods in decaying turbulence. *J. Fluid Mech.* **850**, 336–369.

MEI, R. 1992 An approximate expression for the shear lift force on a spherical particle at finite Reynolds number. *Intl J. Multiphase Flow* **18** (1), 145–147.

MORTENSEN, P.H., ANDERSSON, H.I., GILLISSEN, J.J.J. & BOERSMA, B.J. 2008 Dynamics of prolate ellipsoidal particles in a turbulent channel flow. *Phys. Fluids* **20**, 093302.

MOTOORI, Y., WONG, C. & GOTO, S. 2022 Role of the hierarchy of coherent structures in the transport of heavy small particles in turbulent channel flow. *J. Fluid Mech.* **942**, A3.

PICANO, F., SARDINA, G. & CASCIOLA, C.M. 2009 Spatial development of particle-laden turbulent pipe flow. *Phys. Fluids* **21**, 093305.

REEKS, M.W. 1983 The transport of discrete particles in inhomogeneous turbulence. *J. Aerosol. Sci.* **14** (6), 729–739.

RICHTER, D.H. 2015 Turbulence modification by inertial particles and its influence on the spectral energy budget in planar Couette flow. *Phys. Fluids* **27**, 063304.

RICHTER, D.H. & SULLIVAN, P.P. 2014 Modification of near-wall coherent structures by inertial particles. *Phys. Fluids* **26**, 103304.

RUBINOW, S.I. & KELLER, J.B. 1961 The transverse force on a spinning sphere moving in a viscous fluid. *J. Fluid Mech.* **11** (3), 447–459.

SAFFMAN, P.G. 1965 The lift on a small sphere in a slow shear flow. *J. Fluid Mech.* **22** (2), 385–400.

SARDINA, G., SCHLATTER, P., BRANDT, L., PICANO, F. & CASCIOLA, C.M. 2012 Wall accumulation and spatial localization in particle-laden wall flows. *J. Fluid Mech.* **699**, 50–78.

SCHNEIDERS, L., MEINKE, M. & SCHRÖDER, W. 2017a Direct particle-fluid simulation of Kolmogorov-length-scale size particles in decaying isotropic turbulence. *J. Fluid Mech.* **819**, 188–227.

SCHNEIDERS, L., MEINKE, M. & SCHRÖDER, W. 2017b On the accuracy of Lagrangian point-mass models for heavy non-spherical particles in isotropic turbulence. *Fuel* **201**, 2–14.

SHI, P. & RZEHAK, R. 2019 Lift forces on solid spherical particles in unbounded flows. *Chem. Engng Sci.* **208**, 115145.

SHI, P. & RZEHAK, R. 2020 Lift forces on solid spherical particles in wall-bounded flows. *Chem. Engng Sci.* **211**, 115264.

SHI, P., RZEHAK, R., LUCAS, D. & MAGNAUDET, J. 2021 Drag and lift forces on a rigid sphere immersed in a wall-bounded linear shear flow. *Phys. Rev. Fluids* **6**, 104309.

SHIN, H.H., PORTELA, L.M., SCHAERER, C.E. & MANGIAVACCHI, N. 2022 Dynamics of suspended sediment transport: a direct numerical simulation study. *Intl J. Multiphase Flow* **155**, 104165.

SIKOVSKY, D.P. 2014 Singularity of inertial particle concentration in the viscous sublayer of wall-bounded turbulent flows. *Flow Turbul. Combust.* **92**, 41–64.

SOLDATI, A. & MARCHIOLI, C. 2009 Physics and modelling of turbulent particle deposition and entrainment: review of a systematic study. *Intl J. Multiphase Flow* **35**, 827–839.
TAKEMURA, F., MAGNAUDET, J. & DIMITRAKOPOULOS, P. 2009 Migration and deformation of bubbles rising in a wall-bounded shear flow at finite Reynolds number. *J. Fluid Mech.* **634**, 463–486.

TOSCHI, F. & BODENSCHATZ, E. 2009 Lagrangian properties of particles in turbulence. *Annu. Rev. Fluid Mech.* **41**, 375–404.

UIJTTEWAAL, W.S.J. & OLIEMANS, R.V.A. 1996 Particle dispersion and deposition in direct numerical and large eddy simulations of vertical pipe flows. *Phys. Fluids* **8** (10), 2590–2604.

VASSEUR, P. & COX, R.G. 1977 The lateral migration of spherical particles sedimenting in a stagnant bounded fluid. *J. Fluid Mech.* **80** (3), 561–591.

VREMAN, B., GEURTS, B.J., DEEN, N.G., KUIPERS, J.A.M. & KUERTEN, J.G.M. 2009 Two- and four-way coupled Euler–Lagrangian large-eddy simulation of turbulent particle-laden channel flow. *Flow Turbul. Combust.* **82** (1), 47–71.

WANG, G., FONG, K.O., COLETTI, F., CAPECELATRO, J. & RICHTER, D.H. 2019 Inertial particle velocity and distribution in vertical turbulent channel flow: a numerical and experimental comparison. *Intl J. Multiphase Flow* **120**, 103105.

WANG, G. & RICHTER, D.H. 2019 Two mechanisms of modulation of very-large-scale motions by inertial particles in open channel flow. *J. Fluid Mech.* **868**, 538–559.

WANG, G. & RICHTER, D.H. 2020 Multiscale interaction of inertial particles with turbulent motions in open channel flow. *Phys. Rev. Fluids* **5**, 044307.

WANG, Q., SQUIRES, K.D., CHEN, M. & MCLAUGHLIN, J.B. 1997 On the role of the lift force in turbulence simulations of particle deposition. *Intl J. Multiphase Flow* **23** (4), 749–763.

WRAY, A.A. 1990 Minimal-storage time advancement schemes for spectral methods. *Tech. Rep.* NASA Ames Research Center.

YU, W., VINKOVIC, I. & BUFFAT, M. 2016 Finite-size particles in turbulent channel flow: quadrant analysis and acceleration statistics. *J. Turbul.* **17** (11), 1048–1071.

YU, Z., XIA, Y., GUO, Y. & LIN, J. 2021 Modulation of turbulence intensity by heavy finite-size particles in upward channel flow. *J. Fluid Mech.* **913**, A3.

ZENG, L., BALACHANDAR, S. & FISCHER, P. 2005 Wall-induced forces on a rigid sphere at finite Reynolds number. *J. Fluid Mech.* **536**, 1–25.

ZHANG, H., CUI, Y. & ZHENG, X. 2023 How electrostatic forces affect particle behaviour in turbulent channel flows. *J. Fluid Mech.* **967**, A8.

ZHANG, X., NATHAN, G.J., TIAN, Z.F. & CHIN, R.C. 2022 The dominant mechanisms for each regime of secondary flows in horizontal particle-laden pipe flows. *J. Fluid Mech.* **949**, A10.

ZHAO, L., ANDERSSON, H.I. & GILLISSEN, J.J.J. 2010 Turbulence modulation and drag reduction by spherical particles. *Phys. Fluids* **22**, 081702.

ZHAO, L., MARCHIOLI, C. & ANDERSSON, H.I. 2012 Stokes number effects on particle slip velocity in wall-bounded turbulence and implications for dispersion models. *Phys. Fluids* **24**, 021705.

ZHENG, X., FENG, S. & WANG, P. 2021 Modulation of turbulence by saltating particles on erodible bed surface. *J. Fluid Mech.* **918**, A16.