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: Derivation of Canopy Resistance in Turbulent Flow from First-Order Closure Models
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
Award ID(s):
1754893
PAR ID:
10092844
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Water
Volume:
10
Issue:
12
ISSN:
2073-4441
Page Range / eLocation ID:
1782
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Dimensional analysis suggests that the dissipation length scale ( $$\ell _{{\it\epsilon}}=u_{\star }^{3}/{\it\epsilon}$$ ) is the appropriate scale for the shear-production range of the second-order streamwise structure function in neutrally stratified turbulent shear flows near solid boundaries, including smooth- and rough-wall boundary layers and shear layers above canopies (e.g. crops, forests and cities). These flows have two major characteristics in common: (i) a single velocity scale, i.e. the friction velocity ( $$u_{\star }$$ ) and (ii) the presence of large eddies that scale with an external length scale much larger than the local integral length scale. No assumptions are made about the local integral scale, which is shown to be proportional to $$\ell _{{\it\epsilon}}$$ for the scaling analysis to be consistent with Kolmogorov’s result for the inertial subrange. Here $${\it\epsilon}$$ is the rate of dissipation of turbulent kinetic energy (TKE) that represents the rate of energy cascade in the inertial subrange. The scaling yields a log-law dependence of the second-order streamwise structure function on ( $$r/\ell _{{\it\epsilon}}$$ ), where $$r$$ is the streamwise spatial separation. This scaling law is confirmed by large-eddy simulation (LES) results in the roughness sublayer above a model canopy, where the imbalance between local production and dissipation of TKE is much greater than in the inertial layer of wall turbulence and the local integral scale is affected by two external length scales. Parameters estimated for the log-law dependence on ( $$r/\ell _{{\it\epsilon}}$$ ) are in reasonable agreement with those reported for the inertial layer of wall turbulence. This leads to two important conclusions. Firstly, the validity of the $$\ell _{{\it\epsilon}}$$ -scaling is extended to shear flows with a much greater imbalance between production and dissipation, indicating possible universality of the shear-production range in flows near solid boundaries. Secondly, from a modelling perspective, $$\ell _{{\it\epsilon}}$$ is the appropriate scale to characterize turbulence in shear flows with multiple externally imposed length scales. 
    more » « less
  2. Abstract The need for operational models describing the friction factorfin streams remains undisputed given its utility across a plethora of hydrological and hydraulic applications concerned with shallow inertial flows. For small-scale roughness elements uniformly covering the wetted parameter of a wide channel, the Darcy-Weisbachf = 8(u*/Ub)2is widely used at very high Reynolds numbers, whereu*is friction velocity related to the surface kinematic stress,Ub = Q/Ais bulk velocity,Qis flow rate, andAis cross-sectional area orthogonal to the flow direction. In natural streams, the presence of vegetation introduces additional complications to quantifyingf, the subject of the present work. Turbulent flow through vegetation are characterized by a number of coherent vortical structures: (i) von Karman vortex streets in the lower layers of vegetated canopies, (ii) Kelvin-Helmholtz as well as attached eddies near the vegetation top, and (iii) attached eddies well above the vegetated layer. These vortical structures govern the canonical mixing lengths for momentum transfer and their influence onfis to be derived. The main novelty is that the friction factor of vegetated flow can be expressed asfv = 4Cd(Uv/Ub)2whereUvis the spatially averaged velocity within the canopy volume, andCdis a local drag coefficient per unit frontal area derived to include the aforemontioned layer-wise effects of vortical structures within and above the canopy along with key vegetation properties. The proposed expression is compared with a number of empirical relations derived for vegetation under emergent and submerged conditions as well as numerous data sets covering a wide range of canopy morphology, densities, and rigidity. It is envisaged that the proposed formulation be imminently employed in eco-hydraulics where the interaction between flow and vegetation is being sought. 
    more » « less
  3. Abstract In shallow water systems like coral reefs, bottom friction is an important term in the momentum balance. Parameterizations of bottom friction require a representation of canopy geometry, which can be conceptualized as an array of discrete obstacles or a continuous surface. Here, we assess the implications of using obstacle‐ and surface‐based representations to estimate geometric properties needed to parameterize drag. We collected high‐resolution reef topography data using a scanning multibeam sonar that resolved individual coral colonies within a set of 100‐m2reef patches primarily composed of moundingPoritescorals. The topography measurements yielded 1‐cm resolution continuous surfaces consisting of a single elevation value for each position in a regular horizontal grid. These surfaces were analyzed by (1) defining discrete obstacles and quantifying their properties (dimensions, shapes), and (2) computing properties of the elevation field (root mean square (rms) elevations, rms slopes, spectra). We then computed the roughness density (i.e., frontal area per unit plan area) using both analysis approaches. The obstacle and surface‐based estimates of roughness density did not agree, largely because small‐scale topographic variations contributed significantly to total frontal area. These results challenge the common conceptualization of shallow‐water canopies as obstacle arrays, which may not capture significant contributions of high‐wavenumber roughness to total frontal area. In contrast, the full range of roughness length scales present in natural reefs is captured by the continuous surface representation. Parameterizations of bottom friction over reef topography could potentially be improved by representing the contributions of all length scales to total frontal area and drag. 
    more » « less
  4. Stably stratified roughness sublayer flows are ubiquitous yet remain difficult to represent in models and to interpret using field experiments. Here, continuous high‐frequency potential temperature profiles from the forest floor up to 6.5 times the canopy height observed with distributed temperature sensing (DTS) are used to link eddy topology to roughness sublayer stability correction functions and coupling between air layers within and above the canopy. The experiments are conducted at two forest stands classified as hydrodynamically sparse and dense. Near‐continuous profiles of eddy sizes (length scales) and effective mixing lengths for heat are derived from the observed profiles using a novel conditional sampling approach. The approach utilizes potential temperature isoline fluctuations from a statically stable background state. The transport of potential temperature by an observed eddy is assumed to be conserved (adiabatic movement) and we assume that irreversible heat exchange between the eddy and the surrounding background occurs along the (vertical) periphery of the eddy. This assumption is analogous to Prandtl's mixing‐length concept, where momentum is transported rapidly vertically and then equilibrated with the local mean velocity gradient. A distinct dependence of the derived length scales on background stratification, height above ground, and canopy characteristics emerges from the observed profiles. Implications of these findings for (1) the failure of Monin–Obukhov similarity in the roughness sublayer and (2) above‐canopy flow coupling to the forest floor are examined. The findings have practical applications in terms of analysing similar DTS data sets with the proposed approach, modelling roughness sublayer flows, and interpreting nocturnal eddy covariance measurements above tall forested canopies. 
    more » « less
  5. 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