skip to main content


Title: A direct numerical investigation of two-way interactions in a particle-laden turbulent channel flow
Understanding the two-way interactions between finite-size solid particles and a wall-bounded turbulent flow is crucial in a variety of natural and engineering applications. Previous experimental measurements and particle-resolved direct numerical simulations revealed some interesting phenomena related to particle distribution and turbulence modulation, but their in-depth analyses are largely missing. In this study, turbulent channel flows laden with neutrally buoyant finite-size spherical particles are simulated using the lattice Boltzmann method. Two particle sizes are considered, with diameters equal to 14.45 and 28.9 wall units. To understand the roles played by the particle rotation, two additional simulations with the same particle sizes but no particle rotation are also presented for comparison. Particles of both sizes are found to form clusters. Under the Stokes lubrication corrections, small particles are found to have a stronger preference to form clusters, and their clusters orientate more in the streamwise direction. As a result, small particles reduce the mean flow velocity less than large particles. Particles are also found to result in a more homogeneous distribution of turbulent kinetic energy (TKE) in the wall-normal direction, as well as a more isotropic distribution of TKE among different spatial directions. To understand these turbulence modulation phenomena, we analyse in detail the total and component-wise volume-averaged budget equations of TKE with the simulation data. This budget analysis reveals several mechanisms through which the particles modulate local and global TKE in the particle-laden turbulent channel flow.  more » « less
Award ID(s):
1706130 1513031
NSF-PAR ID:
10110914
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
875
ISSN:
0022-1120
Page Range / eLocation ID:
1096 to 1144
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Current understanding of turbulence modulation by solid particles is incomplete as making reliable predictions on the nature and level of modulation remains a challenging task. Multiple modulation mechanisms may be simultaneously induced by particles, but the lack of reliable methods to identify these mechanisms and quantify their effects hinders a complete understanding of turbulence modulation. In this work, we present a full analysis of the turbulent kinetic energy (TKE) equation for a turbulent channel flow laden with a few fixed particles near the channel walls, in order to investigate how the wall generated turbulence interacts with the particles and how, as a result, the global turbulence statistics are modified. All terms in the budget equations of total and component-wise TKEs are explicitly computed using the data from direct numerical simulations. Particles are found to modify turbulence by two competing mechanisms: the reduction of the intrinsic turbulence production associated with a reduced mean shear due to the resistance imposed by solid particles (the first mechanism), and an additional TKE production mechanism by displacing incoming fluid (the second mechanism). The distribution of TKE in the wall-normal direction is also made more homogeneous due to the significantly enhanced pressure transport of TKE. Finally, the budget analysis of component-wise TKE reveals an enhanced inter-component TKE transfer due to the presence of particles, which leads to a more isotropic distribution of TKE among three velocity components. 
    more » « less
  2. The lattice Boltzmann method is employed to conduct direct numerical simulations of turbulent open channel flows with the presence of finite-size spherical sediment particles. The uniform particles have a diameter of approximately 18 wall units and a density of ρp=2.65ρf, where ρp and ρf are the particle and fluid densities, respectively. Three low particle volume fractions ϕ=0.11%, 0.22%, and 0.44% are used to investigate the particle-turbulence interactions. Simulation results indicate that particles are found to result in a more isotropic distribution of fluid turbulent kinetic energy (TKE) among different velocity components, and a more homogeneous distribution of the fluid TKE in the wall-normal direction. Particles tend to accumulate in the near-wall region due to the settling effect and they preferentially reside in low-speed streaks. The vertical particle volume fraction profiles are self-similar when normalized by the total particle volume fractions. Moreover, several typical transport modes of the sediment particles, such as resuspension, saltation, and rolling, are captured by tracking the trajectories of particles. Finally, the vertical profiles of particle concentration are shown to be consistent with a kinetic model. 
    more » « less
  3. We conducted a spectral analysis of the turbulence kinetic energy (TKE) budget in a bubble plume using particle image velocimetry with fluorescent particles. Our findings confirmed the hypothesis of an inverse energy cascade in the bubble plume, where TKE is transferred from small to large eddies. This is attributed to direct injection of TKE by bubble passages across a wide range of scales, in contrast to canonical shear production of TKE in large scales. Turbulence dissipation was identified as the primary sink of the bubble-produced TKE and occurred at all scales. The decomposition of velocities using the critical length scale of inter-scale energy transfer allowed us to distinguish between large- and small-scale motions in the bubble plume. The large-scale turbulent fluctuations exhibited a skewed distribution and were likely associated with the return flow after bubble passage and the velocities induced by the bubble wake. The small-scale turbulent fluctuations followed a Gaussian distribution relatively well. The large-scale motions contributed to over half of the Reynolds stresses, while there were significant small-scale contributions to the normal stresses near the plume center but not to the shear stress. The large-scale motions in the vorticity field induced a street of vertically elongated vortex pairs, while the small-scale vortices exhibited similar sizes in both horizontal and vertical directions.

     
    more » « less
  4. Particle-laden flows of sedimenting solid particles or droplets in a carrier gas have strong inter-phase coupling. Even at low particle volume fractions, the two-way coupling can be significant due to the large particle to gas density ratio. In this semi-dilute regime, the slip velocity between phases leads to sustained clustering that strongly modulates the overall flow. The analysis of perturbations in homogeneous shear reveals the process by which clusters form: (i) the preferential concentration of inertial particles in the stretching regions of the flow leads to the formation of highly concentrated particle sheets, (ii) the thickness of the latter is controlled by particle-trajectory crossing, which causes a local dispersion of particles, (iii) a transverse Rayleigh–Taylor instability, aided by the shear-induced rotation of the particle sheets towards the gravity normal direction, breaks the planar structure into smaller clusters. Simulations in the Euler–Lagrange formalism are compared to Euler–Euler simulations with the two-fluid and anisotropic-Gaussian methods. It is found that the two-fluid method is unable to capture the particle dispersion due to particle-trajectory crossing and leads instead to the formation of discontinuities. These are removed with the anisotropic-Gaussian method which derives from a kinetic approach with particle-trajectory crossing in mind. 
    more » « less
  5. We present direct numerical simulation results of a wave-current boundary layer in a current-dominated flow regime (wave driven to steady current ratio of 0.34) over bumpy walls for hydraulically smooth flow conditions (wave orbital excursion to roughness ratio of 10). The turbulent, wave-current channel flow has a friction Reynolds number of $350$ and a wave Reynolds number of $351$ . At the lower boundary, a bumpy wall is introduced with a direct forcing immersed boundary method, while the top wall has a free-slip boundary condition. Despite the hydraulically smooth nature of the wave-driven flow, the phase variations of the turbulent statistics for the bumpy wall case were found to vary substantially when compared with the flat wall case. Results show that the addition of weak waves to a steady current over flat walls has a negligible effect on the turbulence or bottom drag. However, the addition of weak waves to a steady current over bumpy walls has a significant effect through enhancement of the Reynolds stress (RS) accompanied by a drag coefficient increase of $11\,\%$ relative to the steady current case. This enhancement occurs just below the top of the roughness elements during the acceleration portion of the wave cycle: Turbulent kinetic energy (TKE) is subsequently transported above the roughness elements to a maximum height of roughly twice the turbulent Stokes length. We analyse the TKE and RS budgets to understand the mechanisms behind the alterations in the turbulence properties due to the bumpy wall. The results provide a mechanistic picture of the differences between bumpy and flat walls in wave-current turbulent boundary layers and illustrate the importance of bumpy features even in weakly energetic wave conditions. 
    more » « less