skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Boundary layer dynamics and bottom friction in combined wave–current flows over large roughness elements
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
Award ID(s):
2123709 1435530 2123707 2123708
PAR ID:
10344679
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
931
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Coral reef roughness produces turbulent boundary layers and bottom stresses that are important for reef metabolism monitoring and reef circulation modeling. However, there is some uncertainty as to whether field methods for estimating bottom stress are applicable in shallow canopy environments as found on coral reefs. Friction velocities () and drag coefficients () were estimated using five independent methods and compared across 14 sites on a shallow forereef (2–9 m deep) in Palau with large and spatially variable coral roughness elements (0.4–1 m tall). The methods included the following: (a) momentum balance closure, (b) log‐fitting to velocity profiles, (c) Reynolds stresses, (d) turbulence dissipation, and (e) roughness characterization from digital elevation models (DEMs). Both velocity profiles and point turbulence measurements indicated good agreement with log‐layer scaling, suggesting that measurements were taken within a well‐developed turbulent boundary layer and that canopy effects were minimal. However, estimated from the DEMs, momentum budget and log‐profile fitting were consistently larger than those estimated from direct turbulence measurements. Near‐bed Reynolds stresses only contributed about 1/3 of the total bottom stress and drag produced by the reef. Thus, effects of topographical heterogeneity that induce mean velocity fluxes, dispersive stresses, and form drag are expected to be important. This decoupling of total drag and local turbulence implies that both rates of mass transfer as well as values of fluxes inferred from concentration measurements may be proportional to smaller, turbulence‐derived values of rather than to those based on larger‐scale flow structure. 
    more » « less
  2. We present direct numerical simulation results of a wave-current boundary layer in a current-dominated flow regime (wave driven to steady current ratio of 0.34) over bumpy walls for hydraulically smooth flow conditions (wave orbital excursion to roughness ratio of 10). The turbulent, wave-current channel flow has a friction Reynolds number of $350$ and a wave Reynolds number of $351$ . At the lower boundary, a bumpy wall is introduced with a direct forcing immersed boundary method, while the top wall has a free-slip boundary condition. Despite the hydraulically smooth nature of the wave-driven flow, the phase variations of the turbulent statistics for the bumpy wall case were found to vary substantially when compared with the flat wall case. Results show that the addition of weak waves to a steady current over flat walls has a negligible effect on the turbulence or bottom drag. However, the addition of weak waves to a steady current over bumpy walls has a significant effect through enhancement of the Reynolds stress (RS) accompanied by a drag coefficient increase of $$11\,\%$$ relative to the steady current case. This enhancement occurs just below the top of the roughness elements during the acceleration portion of the wave cycle: Turbulent kinetic energy (TKE) is subsequently transported above the roughness elements to a maximum height of roughly twice the turbulent Stokes length. We analyse the TKE and RS budgets to understand the mechanisms behind the alterations in the turbulence properties due to the bumpy wall. The results provide a mechanistic picture of the differences between bumpy and flat walls in wave-current turbulent boundary layers and illustrate the importance of bumpy features even in weakly energetic wave conditions. 
    more » « less
  3. Aiming to study the rough-wall turbulent boundary layer structure over differently arranged roughness elements, an experimental study was conducted on flows with regular and random roughness. Varying planform densities of truncated cone roughness elements in a square staggered pattern were investigated. The same planform densities were also investigated in random arrangements. Velocity statistics were measured via two-component laser Doppler velocimetry and stereoscopic particle image velocimetry. Friction velocity, thickness, roughness length and zero-plane displacement, determined from spatially averaged flow statistics, showed only minor differences between the regular and random arrangements at the same density. Recent a priori morphometric and statistical drag prediction methods were evaluated against experimentally determined roughness length. Observed differences between regular and random surface flow parameters were due to the presence of secondary flows which manifest as high-momentum pathways and low-momentum pathways in the streamwise velocity. Contrary to expectation, these secondary flows were present over the random surfaces and not discernible over the regular surfaces. Previously identified streamwise-coherent spanwise roughness heterogeneity does not seem to be present, suggesting that such roughness heterogeneity is not necessary to sustain secondary flows. Evidence suggests that the observed secondary flows were initiated at the front edge of the roughness and sustained over irregular roughness. Due to the secondary flows, local turbulent boundary layer profiles do not scale with local wall shear stress but appear to scale with local turbulent shear stress above the roughness canopy. Additionally, quadrant analysis shows distinct changes in the populations of ejection and sweep events. 
    more » « less
  4. Quantification of roughness effects on free surface flows is unquestionably necessary when describing water and material transport within ecosystems. The conventional hydrodynamic resistance formula empirically shows that the Darcy–Weisbach friction factor f~(r/hw)1/3 describes the energy loss of flowing water caused by small-scale roughness elements characterized by size r (< 
    more » « less
  5. We introduce an analytical model that describes the vertical structure of Ekman boundary layer flows coupled to the Monin-Obukhov Similarity Theory (MOST) surface layer repre- sentation, which is valid for conventionally neutral (CNBL) and stable (SBL) atmospheric conditions. The model is based on a self-similar profile of horizontal stress for both CNBL and SBL flows that merges the classic 3/2 power law profile with a MOST-consistent stress profile in the surface layer. The velocity profiles are then obtained from the Ekman momentum balance equation. The same stress model is used to derive a new self-consistent Geostrophic Drag Law (GDL). We determine the ABL height (h) using an equilibrium boundary layer height model and parameterize the surface heat flux for quasi-steady SBL flows as a function of a prescribed surface temperature cooling rate. The ABL height and GDL equations can then be solved together to obtain the friction velocity (u∗) and the cross-isobaric angle (α0) as a function of known input parameters such as the Geostrophic wind speed and surface roughness (z0). We show that the model predictions agree well with simulation data from the literature and newly generated Large Eddy Simulations (LES). These results indicate that the proposed model provides an efficient and relatively accurate self-consistent approach for predicting the mean wind velocity distribution in CNBL and SBL flows. 
    more » « less