skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Friday, September 29 until 11:59 PM ET on Saturday, September 30 due to maintenance. We apologize for the inconvenience.

Title: A volume penalization immersed boundary method for flow interactions with aquatic vegetation
Abstract: A volume-penalization immersed boundary (VPIB) method was developed to study flow interactions with aquatic vegetation. The model has been validated with data from laboratory experiments and previous high-fidelity models with satisfactory results. Sensitivity analyzes on both penalty parameter and thickness parameter were conducted, and optimal values for these parameters are recommended. The validated model has been applied to study the effects of swaying motion of vegetation stems on the flow dynamics at both vegetate-stem scale and patch scale. The swaying motion of the vegetation stem is prescribed following a cubic law that peaks at the top and decreases to zero at the bottom. At stem-scale, the hydrodynamics depend on the Keulegan Carpenter number (KC), which is defined as the maximum excursion of the vegetation stem to the diameter of the stem. Simulations with three KC values were carried out. For KC≥1, the flow turbulence is significantly enhanced by the swaying motion of the stem, and turbulence becomes more isotropic in the wake. The swaying motion of vegetation stems caused a 5% increase of the bottom shear stress at the shoulders of the stem, and the effect is negligible in the wake. At patch-scale, the hydrodynamics depend on the effective Keulegan Carpenter number based on the patch size of the vegetation patch, and the solid volume fraction for dense vegetation canopy. Solid volume fraction was varied while maintaining the same effective Keulegan Carpenter in the simulations. When the effective Keulgen Carpenter number is small (KC<1), effects of the swaying motion of vegetation stems on the large patch-scale dynamics are not significant, including both the turbulence statistics and the bottom stress.  more » « less
Award ID(s):
Author(s) / Creator(s):
Date Published:
Journal Name:
Advances in water resources
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    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.

    more » « less
  2. 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,kt, in vegetated channels from channel average velocityUand vegetation volume fractionϕwas validated for mobile sediment beds. Second, using data from several studies, the predictedktwas shown to be a good predictor of bed load transport rate,Qs, allowingQsto be predicted fromUandϕfor vegetated channels. The control ofQsbyktwas 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 withkt. 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.

    more » « less
  3. 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
  4. 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.

    more » « less
  5. 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