skip to main content


Title: Impact of gaps on the flow statistics in an emergent rigid canopy
High-resolution large eddy simulations and complementary laboratory experiments using particle image velocimetry were performed to provide a detailed quantitative assessment of flow response to gaps in cylinder arrays. The base canopy consists of a dense array of emergent rigid cylinders placed in a regular staggered pattern. The gaps varied in length from [Formula: see text] to 24, in intervals of 4 d, where d is the diameter of the cylinders. The analysis was performed under subcritical conditions with Froude numbers [Formula: see text] and bulk Reynolds numbers [Formula: see text]. Results show that the gaps affect the flow statistics at the upstream and downstream proximity of the canopy. The affected zone was [Formula: see text] for the mean flow and [Formula: see text] for the second-order statistics. Dimensionless time-averaged streamwise velocity within the gap exhibited minor variability with gap spacing; however, in-plane turbulent kinetic energy, k, showed a consistent decay rate when normalized with that at [Formula: see text] from the beginning of the gap. The emergent canopy acts as a passive turbulence generator for the gap flow for practical purposes. The streamwise dependence of k follows an exponential trend within [Formula: see text] and transitions to a power-law at [Formula: see text]. The substantially lower maximum values of k within the gap compared to k within the canopy evidence a limitation of gap measurements representative of canopy flow statistics. We present a base framework for estimating representative in-canopy statistics from measurements in the gap.  more » « less
Award ID(s):
1753200
NSF-PAR ID:
10407280
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Physics of Fluids
Volume:
34
Issue:
6
ISSN:
1070-6631
Page Range / eLocation ID:
066601
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Large-eddy simulation was used to model turbulent atmospheric surface layer (ASL) flow over canopies composed of streamwise-aligned rows of synthetic trees of height,$h$, and systematically arranged to quantify the response to variable streamwise spacing,$\delta _1$, and spanwise spacing,$\delta _2$, between adjacent trees. The response to spanwise and streamwise heterogeneity has, indeed, been the topic of a sustained research effort: the former resulting in formation of Reynolds-averaged counter-rotating secondary cells, the latter associated with the$k$- and$d$-type response. No study has addressed the confluence of both, and results herein show secondary flow polarity reversal across ‘critical’ values of$\delta _1$and$\delta _2$. For$\delta _2/\delta \lesssim 1$and$\gtrsim 2$, where$\delta$is the flow depth, the counter-rotating secondary cells are aligned such that upwelling and downwelling, respectively, occurs above the elements. The streamwise spacing$\delta _1$regulates this transition, with secondary cell reversal occurring first for the largest$k$-type cases, as elevated turbulence production within the canopy necessitates entrainment of fluid from aloft. The results are interpreted through the lens of a benchmark prognostic closure for effective aerodynamic roughness,$z_{0,{Eff.}} = \alpha \sigma _h$, where$\alpha$is a proportionality constant and$\sigma _h$is height root mean square. We report$\alpha \approx 10^{-1}$, the value reported over many decades for a broad range of rough surfaces, for$k$-type cases at small$\delta _2$, whereas the transition to$d$-type arrangements necessitates larger$\delta _2$. Though preliminary, results highlight the non-trivial response to variation of streamwise and spanwise spacing.

     
    more » « less
  2. This work experimentally investigates the pressure-driven flow of a pure Newtonian fluid through three-dimensional (3D) porous media models. The porous media model consists of square arrays of rods that also could be interpreted as a periodic tandem rod arrangement. We employed a time-resolved three-dimensional particle tracking velocimetry (3D Shake-the-Box) technique for a range of Reynolds numbers [Formula: see text] to observe flow structures and vortex formation between the rods in porous media structures with different porosities of [Formula: see text] which corresponds to the spacing ratio of [Formula: see text], where L is the distance between the centers of the rods, and D is the diameter of the rods. For all the examined cases, we further analyzed the effect of the Reynolds number and the spacing ratio on the instantaneous and averaged patterns of velocity, vorticity, and the other flow parameters after obtaining the two-dimensional velocity fields using the bin-averaging method. We observed both symmetrical and asymmetrical patterns of structure and recirculation regions between the rods depending on the Reynolds number and spacing ratio. Increasing the Reynolds number reduced the symmetrical patterns of flow structures with respect to the centerline of the gap region, while the spacing ratio was randomly affecting the symmetry degree. Vortex shedding was considerable for the two examined high Reynolds numbers of Re = 444 and Re = 890 behind the upstream rod as the porosity increased. The backward movement of the reattachment point has been observed by increasing the Reynolds number. 
    more » « less
  3. The effect of nozzle surface features on the overall atomization behavior of a liquid jet is analyzed in the present computational work by adopting three representative geometries, namely a single X-ray tomography scan of the Engine Combustion Network’s Spray A nozzle (Unprocessed), a spline reconstruction of multiple scans (Educated), and a purely external flow configuration. The latter configuration is often used in fundamental jet atomization studies. Numerically, the two-phase flow is solved based on algebraic volume-of-fluid methodology utilizing the OpenFoam solver, interFoam. Quantitative characterization of the surface features concerning the first two geometries reveals that while both of them have similar levels of cylindrical asymmetries, the nozzle configuration pertaining to the Unprocessed geometry has much larger surface features along the streamwise direction than the Educated geometry. This produces for the Unprocessed configuration a much larger degree of non-axial velocity components in the flow exiting the orifice and also a more pronounced disturbance of the liquid surface in the first few diameters downstream of the nozzle orifice. Furthermore, this heightened level of surface destabilization generates a much shorter intact liquid core length, that is, it produces faster primary atomization. The surprising aspect of this finding is that the differences between the Unprocessed and Educated geometries are of [Formula: see text](1) μm, and they are able to produce [Formula: see text](1) mm effects in the intact liquid core length. In spite of more pronounced atomization for the Unprocessed geometry, the magnitude of the turbulent liquid kinetic energy is roughly the same as the Educated geometry. This highlights the important role of mean field quantities, in particular, non-axial velocity components, in precipitating primary atomization. At the other end of the spectrum, the external-only configuration has the mildest level of surface disturbances in the near field resulting in the longest intact liquid core length. 
    more » « less
  4. Dominant flow features in the near and intermediate wake of a horizontal-axis wind turbine are studied at near field-scale Reynolds numbers. Measurements of the axial velocity component were performed using a nano-scale hot-wire anemometer and analyzed using spectral methods to reveal the extent and evolution of the flow features. Experiments were conducted at a range of Reynolds numbers, of [Formula: see text], based on the rotor diameter and freestream velocity. Five different downstream locations were surveyed, between [Formula: see text], including the near wake, transition to the intermediate wake, and the intermediate wake. Three dominant wake features are identified and studied: the tip vortices, an annular shear layer in the wake core, and wake meandering. The tip vortices are shown to have a broadband influence in the flow in their vicinity, which locally alters the turbulence in that area. It is shown that shedding in the wake core and wake meandering are two distinct and independent low frequency features, and the wake meandering persists into the intermediate wake, whereas the signatures of the core shedding vanish early in the near wake. 
    more » « less
  5. Abstract

    The effect of hydraulic resistance on the downstream evolution of the water surface profilehin a sloping channel covered by a uniform dense rod canopy following the instantaneous collapse of a dam was examined using flume experiments. Near the head of the advancing wavefront, wherehmeets the rods, the conventional picture of a turbulent boundary layer was contrasted to a distributed drag force representation. The details of the boundary layer around the rod and any interferences between rods were lumped into a drag coefficientCd. The study demonstrated the following: In the absence of a canopy, the Ritter solution agreed well with the measurements. When the canopy was represented by an equivalent wall friction as common when employing Manning's formula with constant roughness, it was possible to match the measured wavefront speed but not the precise shape of the water surface profile. However, upon adopting a distributed drag force with a constantCd, the agreement between measured and modeledhwas quite satisfactory at all positions and times. The measurements and model calculations suggested that the shape ofhnear the wavefront was quasilinear with longitudinal distance for a constantCd. The computed constantCd(≈0.4) was surprisingly much smaller than theCd(≈1) reported in uniform flow experiments with staggered cylinders for the same element Reynolds number. This finding suggested that drag reduction mechanisms associated with unsteadiness, nonuniformity, transient waves, and other flow disturbances were more likely to play a role when compared to conventional sheltering effects.

     
    more » « less