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: Subgrid-scale pressure field of scale-enriched large eddy simulations using Gabor modes
With the continuing progress in large eddy simulations (LES), and ever increasing computational resources, it is currently possible to numerically solve the time-dependent and anisotropic large scales of turbulence in a wide variety of flows. For some applications this large-scale resolution is satisfactory. However, a wide range of engineering problems involve flows at very large Reynolds numbers where the subgrid-scale dynamics of a practical LES are critically important to design and yet are out of reach given the com- putational demands of solving the Navier Stokes equations; this difficulty is particularly relevant in wall-bounded turbulence where even the large scales are often below the implied filter width of modest cost wall modeled LES. In this paper we briefly introduce a scale enrichment procedure which leverages spatially and spectrally localized Gabor modes. The method provides a physically consistent description of the small-scale velocity field without solving the full nonlinear equations. The enrichment procedure is appraised against its ability to predict small-scale contributions to the pressure field. We find that the method accurately extrapolates the pressure spectrum and recovers pressure variance of the full field remarkably well when compared to a computationally expensive, highly resolved LES. The analysis is conducted both in a priori and a posteriori settings for the case of homogeneous isotropic turbulence.  more » « less
Award ID(s):
1803378
PAR ID:
10334895
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Physical review fluids
Volume:
6
ISSN:
2469-9918
Page Range / eLocation ID:
124607
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    A turbulence enrichment model for subfilter-scale motions in large eddy simulations (LES) is comprehensively evaluated in the context of a posteriori analysis. The paper further develops the Gabor mode enrichment model first introduced in Ghate & Lele (J. Fluid Mech., vol. 819, 2017, pp. 494–539) by analysing three key requisites of LES enrichment using solenoidal small-scale velocity fields: (a) consistent spectral extrapolation and improvement of resolved single- and two-point second-order correlations; (b) ability to accurately capture the flow physics responsible for temporal decorrelation at small scales; and (c) accurate representation of spatially localized and intermittent interscale energy transfer between scales resolved by the coarse-grid LES and subfilter scales. We argue that the spatially and spectrally localized Gabor wavepackets offer an optimal basis to represent small-scale turbulence within quasi-homogeneous regions, although the alignment of fine-scale vorticity with large-scale strain appears to be somewhat overemphasized. Consequently, we interpret the resulting subfilter scales as those induced by a set of spatially dispersed Burgers–Townsend vortices with orientations determined by the larger scale velocity gradients resolved by the coarse-grid LES. Enrichment of coarse-grid simulations of two high Reynolds number flow configurations, homogeneous isotropic turbulence and a rough-wall turbulent boundary layer show promising results. 
    more » « less
  2. Extensive experimental evidence highlights that scalar turbulence exhibits anomalous diffusion and stronger intermittency levels at small scales compared to that in fluid turbulence. This renders the corresponding subgrid-scale dynamics modeling for scalar turbulence a greater challenge to date. We develop a new large eddy simulation (LES) paradigm for efficiently and dynamically nonlocal LES modeling of the scalar turbulence. To this end, we formulate the underlying nonlocal model starting from the filtered Boltzmann kinetic transport equation, where the divergence of subgrid-scale scalar fluxes emerges as a fractional-order Laplacian term in the filtered advection–diffusion model, coding the corresponding superdiffusive nature of scalar turbulence. Subsequently, we develop a robust data-driven algorithm for estimation of the fractional (noninteger) Laplacian exponent, where we, on the fly, calculate the corresponding model coefficient employing a new dynamic procedure. Our a priori tests show that our new dynamically nonlocal LES paradigm provides better agreement with the ground-truth filtered direct numerical simulation data in comparison to the conventional static and dynamic Prandtl–Smagorinsky models. Moreover, in order to analyze the numerical stability and assessing the model's performance, we carry out comprehensive a posteriori tests. They unanimously illustrate that our new model considerably outperforms other existing functional models, correctly predicting the backscattering phenomena and, at the same time, providing higher correlations at small-to-large filter sizes. We conclude that our proposed nonlocal subgrid-scale model for scalar turbulence is amenable for coarse LES and very large eddy simulation frameworks even with strong anisotropies, applicable to environmental applications. 
    more » « less
  3. Langmuir turbulence consists of Langmuir circulation (LC) generated at the surface of rivers, lakes, bays, and oceans by the interaction between the wind-driven shear and surface gravity waves. In homogeneous shallow water, LC can extend to the bottom of the water column and interact with the bottom boundary layer. Large-eddy simulation (LES) of LC in shallow water was performed with the finite volume method and various forms of subgrid-scale (SGS) model characterized by different near-wall treatments of the SGS eddy viscosity. The wave forcing relative to wind forcing in the LES was set following the field measurements of full-depth LC during the presence of LC engulfing a water column 15 m in depth in the coastal ocean, reported in the literature. It is found that the SGS model can greatly impact the structure of LC in the lower half of the water column. Results are evaluated in terms of (1) the Langmuir turbulence velocity statistics and (2) the lateral (crosswind) length scale and overall cell structure of LC. LES with an eddy viscosity with velocity scale in terms of S and Ω (where S is the norm of the strain rate tensor and Ω is the norm of the vorticity tensor) and a Van Driest wall damping function (referred to as the S-Omega model) is found to provide best agreement with pseudo-spectral LES in terms of the lateral length scale and overall cell structure of LC. Two other SGS models, namely the dynamic Smagorinsky model and the wall-adapting local-eddy viscosity model are found to provide less agreement with pseudo- spectral LES, for example, as they lead to less coherent bottom convergence of the cells and weaker associ ated upward transport of slow downwind moving fluid. Finally, LES with the S-Omega SGS model is also found to lead to good agreement with physical measurements of LC in the coastal ocean in terms of Langmuir turbulence decay during periods of surface heating 
    more » « less
  4. For atmospheric turbulence, multiplying an estimate of the convection velocity with the integral time scale is useful for estimating the integral length scale. Velocity scales that have been used to estimate the convection velocity include the local mean velocity, the ratio of $$e$$-folding length and time scales, and the ratio of a prescribed spatial separation and the time lag at which the space-time autocorrelation peaks. A knowledge gap is the lack of evaluation of these velocity scales directly against the convection velocity, especially for canopy flows where previous studies have reported somewhat inconsistent results. The objective of this work is to assess the ability of each candidate velocity scale to estimate the convection velocity in canopy flows. Firstly, large-eddy simulation (LES) results of neutral canopy flows are used to compare these velocity scales to directly quantified convection velocity. When the direction of interest roughly aligns with the mean pressure gradient force (specifically, for an angle of $$7.5^\circ$$ or smaller), all candidate velocity scales other than the local mean wind component approximate the convection velocity fairly well. When the direction of interest departs from the mean pressure gradient force for more than $$15^\circ$$, the ability of each velocity scale to approximate the convection velocity changes substantially. Secondly, data collected during the Canopy Horizontal Array Turbulence Study (CHATS) are used as an example of interpreting estimates of the convection velocity in the field with the guidance from LES findings. Because observational periods are never perfectly neutral, the guidance does not involve direct comparison between observed and simulated velocity scales, but focuses on uncertainties of velocity scale estimates and potential caution needed when using these estimates. 
    more » « less
  5. Accurate simulation of turbulent flows is of crucial importance in many branches of science and engineering. Direct numerical simulation (DNS) provides the highest fidelity means of capturing all intricate physics of turbulent transport. However, the method is computationally expensive because of the wide range of turbulence scales that must be accounted for in such simulations. Large eddy simulation (LES) provides an alternative. In such simulations, the large scales of the flow are resolved, and the effects of small scales are modelled. Reconstruction of the DNS field from the low-resolution LES is needed for a wide variety of applications. Thus the construction of super-resolution methodologies that can provide this reconstruction has become an area of active research. In this work, a new physics-guided neural network is developed for such a reconstruction. The method leverages the partial differential equation that underlies the flow dynamics in the design of spatio-temporal model architecture. A degradation-based refinement method is also developed to enforce physical constraints and to further reduce the accumulated reconstruction errors over long periods. Detailed DNS data on two turbulent flow configurations are used to assess the performance of the model. 
    more » « less