The Role of Coherent Structures in Subgrid-Scale Energy Transfer in Turbulent Channel Flow

Author(s):  
Vinay Natrajan ◽  
Kenneth Christensen
2012 ◽  
Vol 712 ◽  
pp. 92-128 ◽  
Author(s):  
Jiarong Hong ◽  
Joseph Katz ◽  
Charles Meneveau ◽  
Michael P. Schultz

AbstractThis paper focuses on turbulence structure in a fully developed rough-wall channel flow and its role in subgrid-scale (SGS) energy transfer. Our previous work has shown that eddies of scale comparable to the roughness elements are generated near the wall, and are lifted up rapidly by large-scale coherent structures to flood the flow field well above the roughness sublayer. Utilizing high-resolution and time-resolved particle-image-velocimetry datasets obtained in an optically index-matched facility, we decompose the turbulence into large (${\gt }\lambda $), intermediate ($3\text{{\ndash}} 6k$), roughness ($1\text{{\ndash}} 3k$) and small (${\lt }k$) scales, where $k$ and $\lambda (\lambda / k= 6. 8)$ are roughness height and wavelength, respectively. With decreasing distance from the wall, there is a marked increase in the ‘non-local’ SGS energy flux directly from large to small scales and in the fraction of turbulence dissipated by roughness-scale eddies. Conditional averaging is used to show that a small fraction of the flow volume (e.g. 5 %), which contains the most intense SGS energy transfer events, is responsible for a substantial fraction (50 %) of the energy flux from resolved to subgrid scales. In streamwise wall-normal ($x\text{{\ndash}} y$) planes, the averaged flow structure conditioned on high SGS energy flux exhibits a large inclined shear layer containing negative vorticity, bounded by an ejection below and a sweep above. Near the wall the sweep is dominant, while in the outer layer the ejection is stronger. The peaks of SGS flux and kinetic energy within the inclined layer are spatially displaced from the region of high resolved turbulent kinetic energy. Accordingly, some of the highest correlations occur between spatially displaced resolved velocity gradients and SGS stresses. In wall-parallel $x\text{{\ndash}} z$ planes, the conditional flow field exhibits two pairs of counter-rotating vortices that induce a contracting flow at the peak of SGS flux. Instantaneous realizations in the roughness sublayer show the presence of the counter-rotating vortex pairs at the intersection of two vortex trains, each containing multiple $\lambda $-spaced vortices of the same sign. In the outer layer, the SGS flux peaks within isolated vortex trains that retain the roughness signature, and the distinct pattern of two counter-rotating vortex pairs disappears. To explain the planar signatures, we propose a flow consisting of U-shaped quasi-streamwise vortices that develop as spanwise vorticity is stretched in regions of high streamwise velocity between roughness elements. Flow induced by adjacent legs of the U-shaped structures causes powerful ejections, which lift these vortices away from the wall. As a sweep is transported downstream, its interaction with the roughness generates a series of such events, leading to the formation of inclined vortex trains.


2013 ◽  
Vol 721 ◽  
pp. 541-577 ◽  
Author(s):  
Amin Rasam ◽  
Geert Brethouwer ◽  
Arne V. Johansson

AbstractIn Marstorpet al. (J. Fluid Mech., vol. 639, 2009, pp. 403–432), an explicit algebraic subgrid stress model (EASSM) for large-eddy simulation (LES) was proposed, which was shown to considerably improve LES predictions of rotating and non-rotating turbulent channel flow. In this paper, we extend that work and present a new explicit algebraic subgrid scalar flux model (EASSFM) for LES, based on the modelled transport equation of the subgrid-scale (SGS) scalar flux. The new model is derived using the same kind of methodology that leads to the explicit algebraic scalar flux model of Wikströmet al. (Phys. Fluids, vol. 12, 2000, pp. 688–702). The algebraic form is based on a weak equilibrium assumption and leads to a model that depends on the resolved strain-rate and rotation-rate tensors, the resolved scalar-gradient vector and, importantly, the SGS stress tensor. An accurate prediction of the SGS scalar flux is consequently strongly dependent on an accurate description of the SGS stresses. The new EASSFM is therefore primarily used in connection with the EASSM, since this model can accurately predict SGS stresses. The resulting SGS scalar flux is not necessarily aligned with the resolved scalar gradient, and the inherent dependence on the resolved rotation-rate tensor makes the model suitable for LES of rotating flow applications. The new EASSFM (together with the EASSM) is validated for the case of passive scalar transport in a fully developed turbulent channel flow with and without system rotation. LES results with the new model show good agreement with direct numerical simulation data for both cases. The new model predictions are also compared to those of the dynamic eddy diffusivity model (DEDM) and improvements are observed in the prediction of the resolved and SGS scalar quantities. In the non-rotating case, the model performance is studied at all relevant resolutions, showing that its predictions of the Nusselt number are much less dependent on the grid resolution and are more accurate. In channel flow with wall-normal rotation, where all the SGS stresses and fluxes are non-zero, the new model shows significant improvements over the DEDM predictions of the resolved and SGS quantities.


2019 ◽  
Vol 31 (4) ◽  
pp. 045105 ◽  
Author(s):  
D. Dupuy ◽  
A. Toutant ◽  
F. Bataille

2019 ◽  
Vol 863 ◽  
pp. 1190-1203 ◽  
Author(s):  
Sabarish B. Vadarevu ◽  
Sean Symon ◽  
Simon J. Illingworth ◽  
Ivan Marusic

We study the evolution of velocity fluctuations due to an isolated spatio-temporal impulse using the linearized Navier–Stokes equations. The impulse is introduced as an external body force in incompressible channel flow at $Re_{\unicode[STIX]{x1D70F}}=10\,000$. Velocity fluctuations are defined about the turbulent mean velocity profile. A turbulent eddy viscosity is added to the equations to fix the mean velocity as an exact solution, which also serves to model the dissipative effects of the background turbulence on large-scale fluctuations. An impulsive body force produces flow fields that evolve into coherent structures containing long streamwise velocity streaks that are flanked by quasi-streamwise vortices; some of these impulses produce hairpin vortices. As these vortex–streak structures evolve, they grow in size to be nominally self-similar geometrically with an aspect ratio (streamwise to wall-normal) of approximately 10, while their kinetic energy density decays monotonically. The topology of the vortex–streak structures is not sensitive to the location of the impulse, but is dependent on the direction of the impulsive body force. All of these vortex–streak structures are attached to the wall, and their Reynolds stresses collapse when scaled by distance from the wall, consistent with Townsend’s attached-eddy hypothesis.


Sign in / Sign up

Export Citation Format

Share Document