Laboratory experiments explored the impact of vegetation patchiness on channel‐averaged turbulence and sediment transport. Stems were clustered into 16 randomly distributed circular patches of decreasing diameter. For the same channel velocity, the sediment transport increased with total stem number but decreased as stems were clustered into smaller patch diameters, occupying a smaller fraction of the bed area. The channel‐averaged turbulence, which also declined with increased clustering, was shown to be a good predictor for sediment transport at the channel scale. Previous models for uniform vegetation were adapted to predict both the channel‐averaged turbulence and sediment transport as a function of the total number of stems and degree of clustering, represented by the fraction of bed covered by patches. This provides a way for numerical modelers to represent the impact of subgrid‐scale vegetation patchiness on sediment transport.
- Award ID(s):
- 1945685
- NSF-PAR ID:
- 10314488
- Date Published:
- Journal Name:
- Advances in water resources
- Volume:
- 161
- ISSN:
- 0309-1708
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
Abstract The impacts of aquatic vegetation on bed load transport rate and bedform characteristics were quantified using flume measurements with model emergent vegetation. First, a model for predicting the turbulent kinetic energy,
k t , in vegetated channels from channel average velocityU and vegetation volume fractionϕ was validated for mobile sediment beds. Second, using data from several studies, the predictedk t was shown to be a good predictor of bed load transport rate,Q s , allowingQ s to be predicted fromU andϕ for vegetated channels. The control ofQ s byk t was explained by statistics of individual grain motion recorded by a camera, which showed that the number of sediment grains in motion per bed area was correlated withk t . Third, ripples were observed and characterized in channels with and without model vegetation. For low vegetation solid volume fraction (ϕ ≤ 0.012), the ripple wavelength was constrained by stem spacing. However, at higher vegetation solid volume fraction (ϕ =0.025), distinct ripples were not observed, suggesting a transition to sheet flow, which is sediment transport over a plane bed without the formation of bedforms. The fraction of the bed load flux carried by migrating ripples decreased with increasingϕ , again suggesting that vegetation facilitated the formation of sheet flow. -
A three-dimensional Eulerian two-phase flow solver, SedFoam, has been developed for various sediment transport applications. The solver has demonstrated success in modeling sheet flow and bedforms driven by oscillatory flows using a Reynolds-averaged Navier–Stokes (RANS) formulation. However, the accuracy of the RANS formulation for more complex flows, such as scour around structures, requires further evaluation. SedFoam has recently been enhanced to incorporate two-phase large-eddy simulation (LES) capability. In this study, RANS and LES approaches are tested via a three-dimensional case of wave-induced local scour around a single vertical circular pile. Two laboratory experiments, one with an erodible bed and the other with a rigid bed, were chosen for simulation, with both experiments having a Keulegan-Carpenter (KC) number of 10. The k-ω turbulence closure was selected for the RANS simulation, and the dynamic Lagrangian subgrid closure was chosen for the LES simulation. Numerical results reveal that both RANS and LES simulations can resolve lee-wake vortices, although the vortices are significantly weaker in the RANS simulation. In comparison with the LES results, the RANS approach fails to predict horseshoe vortex with sufficient intensity, leading to an underestimation of scour hole depth development. Although the scour depths develop at a very similar rate in the early stage, the scour depth predicted by the RANS simulation quickly reaches equilibrium, while the LES simulation follows the measured trend. These findings indicate that a turbulence-resolving methodology, i.e. LES, is necessary for accurate scour simulations.more » « less
-
Abstract Laboratory experiments examined the impact of model vegetation on turbulence and resuspension. The turbulent kinetic energy increased with increasing velocity and increasing solid volume fraction, but did not depend on stem diameter. The vegetation‐generated turbulence dominated the total turbulence inside canopies. For the same sediment size, the critical turbulent kinetic energy at which resuspension was initiated was the same for both vegetated and bare beds, which resulted in a critical velocity that decreased with increasing solid volume fraction. Both the critical turbulence and critical velocity for resuspension had no dependence on stem diameter. However, for denser canopies and/or a canopy of smaller stem size, a greater energy slope is required to initiate resuspension. This study provides a way to predict the onset of resuspension in regions with vegetation, an important threshold for sediment transport and landscape evolution.
-
In the coastal ocean, interactions of waves and currents with large roughness elements, similar in size to wave orbital excursions, generate drag and dissipate energy. These boundary layer dynamics differ significantly from well-studied small-scale roughness. To address this problem, we derived spatially and phase-averaged momentum equations for combined wave–current flows over rough bottoms, including the canopy layer containing obstacles. These equations were decomposed into steady and oscillatory parts to investigate the effects of waves on currents, and currents on waves. We applied this framework to analyse large-eddy simulations of combined oscillatory and steady flows over hemisphere arrays (diameter $D$ ), in which current ( $U_c$ ), wave velocity ( $U_w$ ) and period ( $T$ ) were varied. In the steady momentum budget, waves increase drag on the current, and this is balanced by the total stress at the canopy top. Dispersive stresses from oscillatory flow around obstacles are increasingly important as $U_w/U_c$ increases. In the oscillatory momentum budget, acceleration in the canopy is balanced by pressure gradient, added-mass and form drag forces; stress gradients are small compared to other terms. Form drag is increasingly important as the Keulegan–Carpenter number $KC=U_wT/D$ and $U_c/U_w$ increase. Decomposing the drag term illustrates that a quadratic relationship predicts the observed dependences of steady and oscillatory drag on $U_c/U_w$ and $KC$ . For large roughness elements, bottom friction is well represented by a friction factor ( $f_w$ ) defined using combined wave and current velocities in the canopy layer, which is proportional to drag coefficient and frontal area per unit plan area, and increases with $KC$ and $U_c/U_w$ .more » « less