Pore-resolved direct numerical simulations are performed to investigate the interactions between streamflow turbulence and groundwater flow through a randomly packed porous sediment bed for three permeability Reynolds numbers,$$Re_K=2.56$$, 5.17 and 8.94, representative of natural stream or river systems. Time–space averaging is used to quantify the Reynolds stress, form-induced stress, mean flow and shear penetration depths, and mixing length at the sediment–water interface (SWI). The mean flow and shear penetration depths increase with$$Re_K$$and are found to be nonlinear functions of non-dimensional permeability. The peaks and significant values of the Reynolds stresses, form-induced stresses, and pressure variations are shown to occur in the top layer of the bed, which is also confirmed by conducting simulations of just the top layer as roughness elements over an impermeable wall. The probability distribution functions (p.d.f.s) of normalized local bed stress are found to collapse for all Reynolds numbers, and their root-mean-square fluctuations are assumed to follow logarithmic correlations. The fluctuations in local bed stress and resultant drag and lift forces on sediment grains are mainly a result of the top layer; their p.d.f.s are symmetric with heavy tails, and can be well represented by a non-Gaussian model fit. The bed stress statistics and the pressure data at the SWI potentially can be used in providing better boundary conditions in modelling of incipient motion and reach-scale transport in the hyporheic zone.
more »
« less
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
- PAR ID:
- 10412944
- 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
-
-
A large eddy simulation (LES) is performed for a turbulent open channel flow over a porous sediment bed at permeability Reynolds number of ReK∼2.56 (Reτ = 270) representative of aquatic systems. A continuum approach based on the upscaled, volume-averaged Navier−Stokes (VaNS) equations is used by defining smoothly varying porosity across the sediment water interface (SWI) and modeling the drag force in the porous bed using a modified Ergun equation with Forchheimer corrections for inertial terms. The results from the continuum approach are compared with a pore-resolved direct numerical simulation (PR-DNS) in which turbulent flow over a randomly packed sediment bed of monodispersed particles is investigated [Karra et al.,J. Fluid Mech. 971, A23 (2023)] A spatially varying porosity profile generated from the pore-resolved DNS is used in the continuum approach. Mean flow, Reynolds stress statistics, and net momentum exchange between the freestream and the porous bed are compared between the two studies, showing reasonably good agreement. Small deviations within the transitional region between the sediment bed and the freestream as compared to the PR-DNS results are attributed to the local protrusions of particles in a randomly packed bed that are absent in the continuum approach but are present in the PR-DNS. A better representation of the effective permeability in the top transition layer that accounts for roughness effect of exposed particles is necessary. The continuum approach significantly reduces the computational cost, thereby making it suitable to study hyporheic exchange of mass and momentum in large scale aquatic domains with combined influence of bedform and bed roughness.more » « less
-
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
-
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
-
Abstract We took field observations on the shallow shoals of South San Francisco Bay to examine how sediment‐induced stratification affects the mean flow and mixing of momentum and sediment throughout the water column. A Vectrino Profiler measured near‐bed velocity and suspended sediment concentration profiles, which we used to calculate profiles of turbulent sediment and momentum fluxes. Additional turbulence statistics were calculated using data from acoustic Doppler velocimeters placed throughout the water column. Results showed that sediment‐induced stratification, which was set up by strong near‐bed wave shear, can reduce the frictional bottom drag felt by the mean flow. Measured turbulence statistics suggest that this drag reduction is caused by stratification suppressing near‐bed turbulent fluxes and reducing turbulent kinetic energy dissipation. Turbulent sediment fluxes, however, were not shown to be limited by sediment‐induced stratification. Finally, we compared our results to a common model parameterization which characterizes stratification through a stability parameter modification to the turbulent eddy viscosity and suggest a new nondimensional parameter that may be better suited to represent stratification when modeling oscillatory boundary layer flows.more » « less
An official website of the United States government

