skip to main content


Title: Particle Resolved DNS Study of Turbulence Effects on Hyporheic Mixing in Randomly Packed Sediment Beds
Pore-resolved direct numerical simulations (DNS) are used to investigate the interactions between stream-water flow turbulence and groundwater flow through a porous sediment bed in the hyporheic zone. Two permeability Reynolds numbers (2.56 and 5.17), representative of aquatic systems and representing ratio of permeability to viscous length scales, were simulated to understand its influence on the momentum exchange at the sediment-water interface (SWI). A doubleaveraging methodology is used to compute the Reynolds stresses, form-induced stresses, and pressure fluctuations. It is observed that both shear layer and turbulent shear stress penetration increases with ReK. Reynolds and form-induced bed-normal stresses increase with ReK. The peak values of the form-induced stresses for the lower (2.56) and higher (5.17) ReK happen within the top layer of the sediment bed. The sum of turbulent and form-induced pressure fluctuations, analyzed at their respective zero-displacement planes, are statistically similar and can be well approximated by a t location-scale distribution fit providing with a model that could potentially be used to impose boundary conditions at the SWI in reach scale simulations.  more » « less
Award ID(s):
2053248
NSF-PAR ID:
10412944
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
12th International Symposium on Turbulence and Shear Flow Phenomena (TSFP12), Osaka, Japan, July 19--22, 2022
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We quantify the lake dynamics, near‐bottom turbulence, flux of dissolved oxygen (DO) across the sediment‐water interface (SWI) and their interactions during oxygenation in two lakes. Field observations show that the lake dynamics were modified by the bubble plumes, showing enhanced mixing in the near‐field of the plumes. The interaction of the bubble‐induced flow with the internal density structure resulted in downwelling of warm water into the hypolimnion in the far‐field of the plumes. Within the bottom boundary layer (BBL), both lakes show weak oscillating flows primarily induced by seiching. The vertical profile of mean velocity within 0.4 m above the bed follows a logarithmic scaling. One lake shows a larger drag coefficient than those in stationary BBLs, where the classic law‐of‐the‐wall is valid. The injection of oxygen elevated the water column DO and hence, altered the DO flux across the SWI. The gas transfer velocity is driven by turbulence and is correlated with the bottom shear velocity. The thickness of the diffusive boundary layer was found to be consistent with the Batchelor length scale. The dynamics of the surface renewal time follow a log‐normal distribution, and the turbulent integral time scale is comparable to the surface renewal time. The analyses suggest that the effect of bubble plumes on the BBL turbulence is limited and that the canonical scales of turbulence emerge for the time‐average statistics, validating the turbulence scaling of gas transfer velocity in low‐energy lakes.

     
    more » « less
  2. This study discusses turbulent suspension flows of non-Brownian, non-colloidal, neutrally buoyant and rigid spherical particles in a Newtonian fluid over porous media with particles too large to penetrate and move through the porous layer. We consider suspension flows with the solid volume fraction ${{\varPhi _b}}$ ranging from 0 to 0.2, and different wall permeabilities, while porosity is constant at 0.6. Direct numerical simulations with an immersed boundary method are employed to resolve the particles and flow phase, with the volume-averaged Navier–Stokes equations modelling the flow within the porous layer. The results show that in the presence of particles in the free-flow region, the mean velocity and the concentration profiles are altered with increasing porous layer permeability because of the variations in the slip velocity and wall-normal fluctuations at the suspension-porous interface. Furthermore, we show that variations in the stress condition at the interface significantly affect the particle near-wall dynamics and migration toward the channel core, thereby inducing large modulations of the overall flow drag. At the highest volume fraction investigated here, ${{\varPhi _b}}= 0.2$ , the velocity fluctuations and the Reynolds shear stress are found to decrease, and the overall drag increases due to the increase in the particle-induced stresses. 
    more » « less
  3. null (Ed.)
    Abstract Mangrove swamps are extremely productive ecosystems providing many ecological services in coastal regions. The hydrodynamic interactions of mangrove roots and water flow have been proposed as a key element to mitigate erosion. Several studies reveal that precise prediction of the morphological evolution of coastal areas, in the face of global warming and the consequent sea-level rise, requires an understanding of interactions between root porosity (the fraction of the volume of void space over the total volume), water flows, and sediment transport. Water flows around the mangrove prop roots create a complex energetic process that mixes up sediments and generates a depositional region posterior to the roots. In this work, we investigated the boundary layer behind permeable arrays of cylinders (patch) that represent the mangrove roots to explore the impact of patch porosity on the onset of sediment transport. The flow measurements were performed in a vertical plane along the water depth downstream of the mangrove root models. A high-resolution Particle Image Velocimetry (PIV) was used in a flume to observe the impact of porosity on the mean flow, velocity derivatives, skin friction coefficient, and production of turbulent kinetic energy for Reynolds number of 2500 (based on patch diameter length-scale). Here, we proposed a predictive model for critical velocity for incipient motion that takes into account the mangrove roots porosity and the near-bed turbulence effect. It is found that the patch with the $$\phi =47\%$$ ϕ = 47 % porosity, has the maximum critical velocity over which the sediment transport initiates. We found the optimum porosity has the minimum sediment erosion and creates negative vorticity sources near the bed that increases the critical velocity. This signifies an optimum porosity for the onset of sediment transport consistent with the porosity of mangroves in nature. The phenomenological model is elucidated based on an analysis of the vorticity evolution equation for viscous incompressible flows. For the optimum porous patch, a sink of vorticity was formed which yielded to lower the near-bed turbulence and vorticity. The minimum velocity fluctuations were sufficient to initiate the boundary layer transition, however, the viscous dissipation dominated the turbulence production to obstruct the sediment transport. This work identified the pivotal role of mangrove root porosity in sediment transport in terms of velocity and its derivatives in wall-bounded flows. Our work also provides insight into the sediment transport and erosion processes that govern the evolution of the shapes of shorelines. 
    more » « less
  4. Abstract

    Dissolved Oxygen (DO) fluxes across the air‐water and sediment‐water interface (AWI and SWI) are two major processes that govern the amount of oxygen available to living organisms in aquatic ecosystems. Aquatic vegetation generates different scales of turbulence that change the flow structure and affect gas transfer mechanisms at AWI and SWI. A series of laboratory experiments with rigid cylinder arrays to mimic vegetation was conducted in a recirculating race‐track flume with a lightweight sediment bed. 2D Planar Particle Image Velocimetry was used to characterize the flow field under different submergence ratios and array densities to access the effect of vegetation‐generated turbulence on gas transfer. Gas transfer rate across AWI was determined by DO re‐aeration curves. The effective diffusion coefficient for gas transfer flux across SWI was estimated by the difference between near‐bed and near‐surface DO concentrations. When sediment begins to mobilize, near‐bed suspended sediment provides a negative buoyancy term that increases the critical Reynolds number for the surface gas transfer process according to a modified Surface Renewal model for vegetated flows. A new Reynolds number dependence model using near‐bed turbulent kinetic energy as an indicator is proposed to provide a universal prediction for the interfacial flux across SWI in flows with aquatic vegetation. This study provides critical information and useful models for future studies on water quality management and ecosystem restoration in natural water environments such as lakes, rivers, and wetlands.

     
    more » « less
  5. An incoming canonical spatially developing turbulent boundary layer (SDTBL) over a 2-D curved hill is numerically investigated via the Reynolds-averaged Navier–Stokes (RANS) equations plus two eddy-viscosity models: the K−ω SST (henceforth SST) and the Spalart–Allmaras (henceforth SA) turbulence models. A spatially evolving thermal boundary layer has also been included, assuming temperature as a passive scalar (Pr = 0.71) and a turbulent Prandtl number, Prt, of 0.90 for wall-normal turbulent heat flux modeling. The complex flow with a combined strong adverse/favorable streamline curvature-driven pressure gradient caused by concave/convex surface curvatures has been replicated from wind-tunnel experiments from the literature, and the measured velocity and pressure fields have been used for validation purposes (the thermal field was not experimentally measured). Furthermore, direct numerical simulation (DNS) databases from the literature were also employed for the incoming turbulent flow assessment. Concerning first-order statistics, the SA model demonstrated a better agreement with experiments where the turbulent boundary layer remained attached, for instance, in Cp, Cf, and Us predictions. Conversely, the SST model has shown a slightly better match with experiments over the flow separation zone (in terms of Cp and Cf) and in Us profiles just upstream of the bubble. The Reynolds analogy, based on the St/(Cf/2) ratio, holds in zero-pressure gradient (ZPG) zones; however, it is significantly deteriorated by the presence of streamline curvature-driven pressure gradient, particularly due to concave wall curvature or adverse-pressure gradient (APG). In terms of second-order statistics, the SST model has better captured the positively correlated characteristics of u′ and v′ or positive Reynolds shear stresses ( > 0) inside the recirculating zone. Very strong APG induced outer secondary peaks in and turbulence production as well as an evident negative slope on the constant shear layer. 
    more » « less