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: Revisiting a Drag Partition Model For Canopy-Like Roughness Elements
Abstract Turbulent flows over a large surface area (S) covered bynobstacles experience an overall drag due to the presence of the ground and the protruding obstacles into the flow. The drag partition between the roughness obstacles and the ground is analyzed using an analytical model proposed by Raupach (Boundary-Layer Meteorol 60:375-395, 1992) and is hereafter referred to as R92. The R92 is based on the premise that the wake behind an isolated roughness element can be described by a shelter areaAand a shelter volumeV. The individual sizes ofAandVwithout any interference from other obstacles can be determined from scaling analysis for the spread of wakes. To upscale from an individual roughness element ton/Selements where wakes may interact, R92 adopted a background stress re-normalizing instead of reducingAorVwith each element addition. This work demonstrates that R92’s approach results in a linear background stress reduction inAandVonly when the ratio ofn/Sis small, due to a low probability of wake interactions. This probabilistic nature suggests that up-scaling from individual to multiple roughness elements can be re-formulated using stochastic averaging methods proposed here. The two approaches are shown to recover R92 under plausible conditions. An alternative scaling for the shelter volume is also proposed here using thermodynamic arguments of work and dissipation though the final outcome remains similar to R92. Comparisons between R92 and available data spanning more than two decades after R92 on blocks and vegetation-like roughness elements confirm the practical utility of R92. The agreement between R92 and this updated databases of experiments and simulations confirm the potential use of R92 in large-scale models provided that the relevant parameters accommodate certain features of the roughness element type (cube versus vegetation-like) and, to a lesser extent, their configuration throughoutS. Last, a comparison between R92 and models based on first-order closure principles with constant mixing length suggests that R92 can outperform such models when evaluated across a wide range of roughness densities.  more » « less
Award ID(s):
2028633
PAR ID:
10550258
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Boundary-Layer Meteorology
Volume:
190
Issue:
11
ISSN:
0006-8314
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We present experiments on oscillating hydrofoils undergoing combined heaving and pitching motions, paying particular attention to connections between propulsive efficiency and coherent wake features extracted using modal analysis. Time-averaged forces and particle image velocimetry measurements of the flow field downstream of the foil are presented for a Reynolds number of Re=11000 and Strouhal numbers in the range St=0.16--0.35. These conditions produce 2S and 2P wake patterns, as well as a near-momentumless wake structure. A triple decomposition using the optimized dynamic mode decomposition method is employed to identify dominant modal components (or coherent structures) in the wake. These structures can be connected to wake instabilities predicted using spatial stability analyses. Examining the modal components of the wake provides insightful explanations into the transition from drag to thrust production, and conditions that lead to peak propulsive efficiency. In particular, we find modes that correspond to the primary vortex development in the wakes. Other modal components capture elements of bluff body shedding at Strouhal numbers below the optimum for peak propulsive efficiency and characteristics of separation for Strouhal numbers higher than the optimum. 
    more » « less
  2. We present experiments on oscillating hydrofoils undergoing combined heaving and pitching motions, paying particular attention to connections between propulsive efficiency and coherent wake features extracted using modal analysis. Time-averaged forces and particle image velocimetry measurements of the flow field downstream of the foil are presented for a Reynolds number of Re=11000 and Strouhal numbers in the range St=0.16–0.35 . These conditions produce 2S and 2P wake patterns, as well as a near-momentumless wake structure. A triple decomposition using the optimized dynamic mode decomposition method is employed to identify dominant modal components (or coherent structures) in the wake. These structures can be connected to wake instabilities predicted using spatial stability analyses. Examining the modal components of the wake provides insightful explanations into the transition from drag to thrust production, and conditions that lead to peak propulsive efficiency. In particular, we find modes that correspond to the primary vortex development in the wakes. Other modal components capture elements of bluff body shedding at Strouhal numbers below the optimum for peak propulsive efficiency and characteristics of separation for Strouhal numbers higher than the optimum. 
    more » « less
  3. 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
  4. Large Eddy Simulations (LES) of neutral flow over regular arrays of cuboids are conducted to explore connections between momentum (z 0m ) and scalar (z 0s ) roughness lengths in urban environments, and how they are influenced by surface geometry. As LES resolves the obstacles but not the micro‐scale boundary layers attached to them, the aforementioned roughness lengths are analyzed at two distinct spatial scales. At the micro‐scale (roughness of individual facets, e.g. roofs), it is assumed that both momentum and scalar transfer are governed by accepted arguments for smooth walls that form the basis for the LES wall model. At the macro‐scale, the roughness lengths are representative of the aggregate effects of momentum and scalar transfer over the resolved roughness elements of the whole surface, and hence they are directly computed from the LES. The results indicate that morphologically‐based parameterizations for macro‐scale z 0m are adequate overall. The relation between the momentum and scalar macro‐roughness values, as conventionally represented by log(z 0m /z 0s ) and assumed to scale with urn:x-wiley:00359009:media:qj3839:qj3839-math-0001 (where Re * is a roughness Reynolds number), is then interpreted using surface renewal theory (SRT). SRT predicts n = 1/4 when only Kolmogorov‐scale eddies dominate the scalar exchange, whereas n = 1/2 is predicted when large eddies limit the renewal dynamics. The latter is found to better capture the LES results. However, both scaling relations indicate that z 0s decreases when z 0m increases for typical urban geometries and scales. This is opposite to how their relation is usually modeled for urban canopies (i.e. z 0s /z 0m is a fixed value smaller than unity). 
    more » « less
  5. Abstract Turbulent mixing of scalars within canopies is investigated using a flume experiment with canopy‐like rods of heighthmounted to the channel bed. The data comprised a time sequence of high‐resolution images of a dye recorded in a plane parallel to the bed atz/h= 0.2. Image processing shows that von Kármán wakes shed by canopy drag and downward turbulent transport from upper canopy layers impose distinct scaling regimes on the scalar spectrum. Measures from information theory are then used to explore the dominant directionality of the interaction between small and large scales underlying these two spectral regimes, showing that the arrival of sweeps from aloft establishes an inertial‐range spectrum with forward “information” cascade. In contrast, wake growth with downstream distance leads to persistent upscale transfer (inverse cascade) of scalar variance, which hints at their nondiffusive character and the significance of the stem diameter as an active length scale in canopy turbulence. 
    more » « less