Two common definitions of the spatially local rate of kinetic energy cascade at some scale
This content will become publicly available on November 10, 2024
Drag for wallbounded flows is directly related to the spatial flux of spanwise vorticity outward from the wall. In turbulent flows a key contribution to this wallnormal flux arises from nonlinear advection and stretching of vorticity, interpretable as a cascade. We study this process using numerical simulation data of turbulent channel flow at friction Reynolds number
 Award ID(s):
 2103874
 NSFPAR ID:
 10511081
 Publisher / Repository:
 CUP
 Date Published:
 Journal Name:
 Journal of Fluid Mechanics
 Volume:
 974
 ISSN:
 00221120
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

in turbulent flows are (i) the cubic velocity difference term appearing in the ‘scaleintegrated local Kolmogorov–Hill’ equation (structurefunction approach), and (ii) the subfilterscale energy flux term in the transport equation for subgridscale kinetic energy (filtering approach). We perform a comparative study of both quantities based on direct numerical simulation data of isotropic turbulence at Taylorscale Reynolds number 1250. While in the past observations of negative subfilterscale energy flux (backscatter) have led to debates regarding interpretation and relevance of such observations, we argue that the interpretation of the local structurefunctionbased cascade rate definition is unambiguous since it arises from a divergence term in scale space. Conditional averaging is used to explore the relationship between the local cascade rate and the local filtered viscous dissipation rate as well as filtered velocity gradient tensor properties such as its invariants. We find statistically robust evidence of inverse cascade when both the largescale rotation rate is strong and the largescale strain rate is weak. Even stronger net inverse cascading is observed in the ‘vortex compression’$\ell$ ,$R>0$ quadrant, where$Q>0$ and$R$ are velocity gradient invariants. Qualitatively similar but quantitatively much weaker trends are observed for the conditionally averaged subfilterscale energy flux. Flow visualizations show consistent trends, namely that spatially, the inverse cascade events appear to be located within largescale vortices, specifically in subregions when$Q$ is large.$R$ 
Researchers have hypothesized that the poststall lift benefit of bird’s alular feathers, or alula, stems from the maintenance of an attached leadingedge vortex (LEV) over their thinprofiled, outer hand wing. Here, we investigate the connection between the alula and LEV attachment via flow measurements in a wind tunnel. We show that a model alula, whose wetted area is 1 % that of the wing, stabilizes a recirculatory afttilted LEV on a steadily translating unswept wing at poststall angles of attack. The attached vortex is the result of the alula’s ability to smoothly merge otherwise separate leading and sideedge vortical flows. We identify two key processes that facilitate this merging: (i) the steering of spanwise vorticity generated at the wing’s leading edge back to the wing plane and (ii) an aftlocated wall jet of highmagnitude roottotip spanwise flow ( ${>}80\,\%$ that of the freestream velocity). The former feature induces LEV rollup while the latter tilts LEV vorticity aft and evacuates this flow toward the wing tip via an outboard vorticity flux. We identify the alula’s streamwise position (relative to the leading edge of the thin wing) as important for vortex steering and the alula’s cant angle as important for highmagnitude spanwise flow generation. These findings advance our understanding of the likely ways birds leverage LEVs to augment slow flight.more » « less

null (Ed.)We use an online database of a turbulent channelflow simulation at $Re_\tau =1000$ (Graham et al. J. Turbul. , vol. 17, issue 2, 2016, pp. 181–215) to determine the origin of vorticity in the nearwall buffer layer. Following an experimental study of Sheng et al. ( J. Fluid Mech. , vol. 633, 2009, pp.17–60), we identify typical ‘ejection’ and ‘sweep’ events in the buffer layer by local minima/maxima of the wall stress. In contrast to their conjecture, however, we find that vortex lifting from the wall is not a discrete event requiring $\sim$ 1 viscous time and $\sim$ 10 wall units, but is instead a distributed process over a space–time region at least $1\sim 2$ orders of magnitude larger in extent. To reach this conclusion, we exploit a rigorous mathematical theory of vorticity dynamics for Navier–Stokes solutions, in terms of stochastic Lagrangian flows and stochastic Cauchy invariants, conserved on average backward in time. This theory yields exact expressions for vorticity inside the flow domain in terms of vorticity at the wall, as transported by viscous diffusion and by nonlinear advection, stretching and rotation. We show that Lagrangian chaos observed in the buffer layer can be reconciled with saturated vorticity magnitude by ‘virtual reconnection’: although the Eulerian vorticity field in the viscous sublayer has a single sign of spanwise component, opposite signs of Lagrangian vorticity evolve by rotation and cancel by viscous destruction. Our analysis reveals many unifying features of classical fluids and quantum superfluids. We argue that ‘bundles’ of quantized vortices in superfluid turbulence will also exhibit stochastic Lagrangian dynamics and satisfy stochastic conservation laws resulting from particle relabelling symmetry.more » « less

Largeeddy simulation was used to model turbulent atmospheric surface layer (ASL) flow over canopies composed of streamwisealigned rows of synthetic trees of height,
, and systematically arranged to quantify the response to variable streamwise spacing,$h$ , and spanwise spacing,$\delta _1$ , 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 Reynoldsaveraged counterrotating secondary cells, the latter associated with the$\delta _2$  and$k$ type response. No study has addressed the confluence of both, and results herein show secondary flow polarity reversal across ‘critical’ values of$d$ and$\delta _1$ . For$\delta _2$ and$\delta _2/\delta \lesssim 1$ , where$\gtrsim 2$ is the flow depth, the counterrotating secondary cells are aligned such that upwelling and downwelling, respectively, occurs above the elements. The streamwise spacing$\delta$ regulates this transition, with secondary cell reversal occurring first for the largest$\delta _1$ 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,$k$ , where$z_{0,{Eff.}} = \alpha \sigma _h$ is a proportionality constant and$\alpha$ is height root mean square. We report$\sigma _h$ , the value reported over many decades for a broad range of rough surfaces, for$\alpha \approx 10^{1}$ type cases at small$k$ , whereas the transition to$\delta _2$ type arrangements necessitates larger$d$ . Though preliminary, results highlight the nontrivial response to variation of streamwise and spanwise spacing.$\delta _2$ 
The statistical properties of uniform momentum zones (UMZs) are extracted from laboratory and field measurements in rough wall turbulent boundary layers to formulate a set of stochastic models for the simulation of instantaneous velocity profiles. A spatiotemporally resolved velocity dataset, covering a field of view of
, was obtained in the atmospheric surface layer using superlargescale particle image velocimetry (SLPIV), as part of the Grandscale Atmospheric Imaging Apparatus (GAIA). Wind tunnel data from a previous study are included for comparison (Heisel$8 \times 9\,{\rm m}^2$ et al. ,J. Fluid Mech. , vol. 887, 2020, R1). The probability density function of UMZ attributes such as their thickness, modal velocity and averaged vertical velocity are built at varying elevations and modelled using lognormal and Gaussian distributions. Inverse transform sampling of the distributions is used to generate synthetic steplike velocity profiles that are spatially and temporally uncorrelated. Results show that in the wide range of wallnormal distances and up to$Re_\tau$ investigated here, shear velocity scaling is manifested in the velocity jump across shear interfaces between adjacent UMZs, and attached eddy behaviour is observed in the linear proportionality between UMZ thickness and their wall normal location. These very same characteristics are recovered in the generated instantaneous profiles, using both fully stochastic and datadriven hybrid stochastic (DHS) models, which address, in different ways, the coupling between modal velocities and UMZ thickness. Our method provides a stochastic approach for generating an ensemble of instantaneous velocity profiles, consistent with the structural organisation of UMZs, where the ensemble reproduces the logarithmic mean velocity profile and recovers significant portions of the Reynolds stresses and, thus, of the streamwise and vertical velocity variability.$\sim O(10^6)$