skip to main content


Title: Flow modulation by a few fixed spherical particles in a turbulent channel flow
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
Award ID(s):
1706130 1513031
NSF-PAR ID:
10189096
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
884
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  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. null (Ed.)
    Abstract Turbulence driven by wind and waves controls the transport of heat, momentum, and matter in the ocean surface boundary layer (OSBL). For realistic ocean conditions, winds and waves are often neither aligned nor constant, for example, when winds turn rapidly. Based on a Large Eddy Simulation (LES) method, which captures shear-driven turbulence (ST) and Langmuir turbulence (LT) driven by the Craik-Leibovich vortex force, we investigate the OSBL response to abruptly turning winds. We design idealized LES experiments, whose winds are initially constant to equilibrate OSBL turbulence before abruptly turning 90° either cyclonically or anticyclonically. The transient Stokes drift for LT is estimated from a spectral wave model. The OSBL response includes three successive stages that follow the change in direction. During stage 1, turbulent kinetic energy (TKE) decreases due to reduced TKE production. Stage 2 is characterized by TKE increasing with TKE shear production recovering and exceeding TKE dissipation. Transient TKE levels may exceed their stationary values due to inertial resonance and non-equilibrium turbulence. Turbulence relaxes to its equilibrium state at stage 3, but LT still adjusts due to slowly developing waves. During stages 1 and 2, greatly misaligned wind and waves lead to Eulerian TKE production exceeding Stokes TKE production. A Reynolds stress budget analysis and Reynolds-averaged Navier-Stokes equation models indicate that Stokes production furthermore drives the OSBL response. The Coriolis effects result in asymmetrical OSBL responses to wind turning directions. Our results suggest that transient wind conditions play a key role in understanding realistic OSBL dynamics. 
    more » « less
  5. null (Ed.)
    Dimensional analysis suggests that the dissipation length scale ( $\ell _{{\it\epsilon}}=u_{\star }^{3}/{\it\epsilon}$ ) is the appropriate scale for the shear-production range of the second-order streamwise structure function in neutrally stratified turbulent shear flows near solid boundaries, including smooth- and rough-wall boundary layers and shear layers above canopies (e.g. crops, forests and cities). These flows have two major characteristics in common: (i) a single velocity scale, i.e. the friction velocity ( $u_{\star }$ ) and (ii) the presence of large eddies that scale with an external length scale much larger than the local integral length scale. No assumptions are made about the local integral scale, which is shown to be proportional to $\ell _{{\it\epsilon}}$ for the scaling analysis to be consistent with Kolmogorov’s result for the inertial subrange. Here ${\it\epsilon}$ is the rate of dissipation of turbulent kinetic energy (TKE) that represents the rate of energy cascade in the inertial subrange. The scaling yields a log-law dependence of the second-order streamwise structure function on ( $r/\ell _{{\it\epsilon}}$ ), where $r$ is the streamwise spatial separation. This scaling law is confirmed by large-eddy simulation (LES) results in the roughness sublayer above a model canopy, where the imbalance between local production and dissipation of TKE is much greater than in the inertial layer of wall turbulence and the local integral scale is affected by two external length scales. Parameters estimated for the log-law dependence on ( $r/\ell _{{\it\epsilon}}$ ) are in reasonable agreement with those reported for the inertial layer of wall turbulence. This leads to two important conclusions. Firstly, the validity of the $\ell _{{\it\epsilon}}$ -scaling is extended to shear flows with a much greater imbalance between production and dissipation, indicating possible universality of the shear-production range in flows near solid boundaries. Secondly, from a modelling perspective, $\ell _{{\it\epsilon}}$ is the appropriate scale to characterize turbulence in shear flows with multiple externally imposed length scales. 
    more » « less