An official website of the United States government
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
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
- PAR ID:
- 10110914
- 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
-
-
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
-
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
-
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
-
null (Ed.)Abstract The performance of different versions of the discrete random walk models in turbulent flows with nonuniform normal root-mean-square (RMS) velocity fluctuations and turbulence time scales were carefully investigated. The OpenFOAM v2−f low Reynolds number turbulence model was used for evaluating the fully developed streamwise velocity and the wall-normal RMS velocity fluctuations profiles in a turbulent channel flow. The results were then used in an in-house matlab particle tracking code, including the drag and Brownian forces, and the trajectories of randomly injected point-particles with diameters ranging from 10 nm to 30 μm were evaluated under the one-way coupling assumption. The distributions and deposition velocities of fluid-tracer and finite-size particles were evaluated using the conventional-discrete random walk (DRW) model, the modified-DRW model including the velocity gradient drift correction, and the new improved-DRW model including the velocity and time gradient drift terms. It was shown that the conventional-DRW model leads to superfluous migration of fluid-point particles toward the wall and erroneous particle deposition rate. The concentration profiles of tracer particles obtained by using the modified-DRW model still are not uniform. However, it was shown that the new improved-DRW model with the velocity and time scale drift corrections leads to uniform distributions for fluid-point particles and reasonable concentration profiles for finite-size heavy particles. In addition, good agreement was found between the estimated deposition velocities of different size particles by the new improved-DRW model with the available data.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2019.509
