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.
We developed a novel autonomously dynamic nonlocal turbulence model for the large and very large eddy simulation (LES, VLES) of homogeneous isotropic turbulent flows. The model is based on a generalized (integer-to-noninteger)-order Laplacian of the filtered velocity field, and a novel dynamic model has been formulated to avoid the need for tuning the model constant. Three data-driven approaches were introduced for the determination of the fractional-order to have a model that is totally free of any tuning parameter. Our analysis includes both the a priori and the a posteriori tests. In the former test, using a high-fidelity and well-resolved dataset from direct numerical simulations (DNSs), we computed the correlation coefficients for the stress components of the subgrid-scale (SGS) stress tensor and the one we get directly from the DNS results. Moreover, we compared the probability density function of the ensemble-averaged SGS forces for different filter sizes. In the latter, we employed our new model along with other conventional models including the static and dynamic Smagorinsky models into our pseudo-spectral solver and tested the final predicted quantities. The results of the newly developed model exhibit an expressive agreement with the ground-truth DNS results in all components of the SGS stress and forces. Also, the model exhibits promising results in the VLES region as well as the LES region, which could be remarkably important for cost-efficient nonlocal turbulence modeling, e.g., in meteorological and environmental applications.
more » « less- Award ID(s):
- 1923201
- NSF-PAR ID:
- 10363364
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- Physics of Fluids
- Volume:
- 34
- Issue:
- 3
- ISSN:
- 1070-6631
- Page Range / eLocation ID:
- Article No. 035104
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
The traditional approach of using the Monin–Obukhov similarity theory (MOST) to model near-surface processes in large-eddy simulations (LESs) can lead to significant errors in natural convection. In this study, we propose an alternative approach based on feedforward neural networks (FNNs) trained on output from direct numerical simulation (DNS). To evaluate the performance, we conduct both a priori and a posteriori tests. In the a priori (offline) tests, we compare the statistics of the surface shear stress and heat flux, computed from filtered DNS input variables, to the stress and flux obtained from the filtered DNS. Additionally, we investigate the importance of various input features using the Shapley additive explanations value and the conditional average of the filter grid cells. In the a posteriori (online) tests, we implement the trained models in the System for Atmospheric Modeling (SAM) LES and compare the LES-generated surface shear stress and heat flux with those in the DNS. Our findings reveal that vertical velocity, a traditionally overlooked flow quantity, is one of the most important input features for determining the wall fluxes. Increasing the number of input features improves the a priori test results but does not always improve the model performance in the a posteriori tests because of the differences in input variables between the LES and DNS. Last, we show that physics-aware FNN models trained with logarithmic and scaled parameters can well extrapolate to more intense convection scenarios than in the training dataset, whereas those trained with primitive flow quantities cannot.more » « less
Significance Statement The traditional near-surface turbulence model, based on a shear-dominated boundary layer flow, does not represent near-surface turbulence in natural convection. Using a feedforward neural network (FNN), we can construct a more accurate model that better represents the near-surface turbulence in various flows and reveals previously overlooked controlling factors and process interactions. Our study shows that the FNN-generated models outperform the traditional model and highlight the importance of the near-surface vertical velocity. Furthermore, the physics-aware FNN models exhibit the potential to extrapolate to convective flows of various intensities beyond the range of the training dataset, suggesting their broader applicability for more accurate modeling of near-surface turbulence.
-
The transport and deposition of firebrand particles is an important fire spread mechanism in wildland fires. These particles can be transported by wind over large distances and can ignite secondary fires upon landing. The transport of firebrands by wind is a complex, multiscale process that is largely controlled by interactions between the firebrand particles and the atmospheric wind. To account for the complex temporal evolution of atmospheric turbulence over large scales, the use of large-eddy simulation (LES) techniques is necessary. However, filtering of subgrid-scale (SGS) turbulence in LES hinders the accuracy of particle transport models. In this work, we employ a Lagrangian SGS model in an LES framework to investigate the effects of small-scale turbulence on the transport of mass- and size-changing firebrand particles. The impact of SGS turbulence was analyzed by comparing landing and trajectory statistics for firebrand and regular (fixed size and mass) particles under different Stokes numbers. It was found that the presence of SGS turbulence modifies the particle transport behavior, which is characterized by smaller spanwise dispersions but larger travel distances along the streamwise direction compared with particles under no SGS turbulence. As expected, the enhanced velocity field produced by the SGS model has larger influence on the statistics of firebrand particles compared with regular particles due to the time-evolving reduction in particle mass and size induced by pyrolysis.more » « less
-
We employ numerically implicit subgrid-scale modeling provided by the well-known streamlined upwind/Petrov–Galerkin stabilization for the finite element discretization of advection–diffusion problems in a Large Eddy Simulation (LES) approach. Whereas its original purpose was to provide sufficient algorithmic dissipation for a stable and convergent numerical method, more recently, it has been utilized as a subgrid-scale (SGS) model to account for the effect of small scales, unresolvable by the discretization. The freestream Mach number is 2.5, and direct comparison with a DNS database from our research group, as well as with experiments from the literature of adiabatic supersonic spatially turbulent boundary layers, is performed. Turbulent inflow conditions are generated via our dynamic rescaling–recycling approach, recently extended to high-speed flows. Focus is given to the assessment of the resolved Reynolds stresses. In addition, flow visualization is performed to obtain a much better insight into the physics of the flow. A weak compressibility effect is observed on thermal turbulent structures based on two-point correlations (IC vs. supersonic). The Reynolds analogy (u′ vs. t′) approximately holds for the supersonic regime, but to a lesser extent than previously observed in incompressible (IC) turbulent boundary layers, where temperature was assumed as a passive scalar. A much longer power law behavior of the mean streamwise velocity is computed in the outer region when compared to the log law at Mach 2.5. Implicit LES has shown very good performance in Mach 2.5 adiabatic flat plates in terms of the mean flow (i.e., Cf and UVD+). iLES significantly overpredicts the peak values of u′, and consequently Reynolds shear stress peaks, in the buffer layer. However, excellent agreement between the turbulence intensities and Reynolds shear stresses is accomplished in the outer region by the present iLES with respect to the external DNS database at similar Reynolds numbers.more » « less
-
In this study, we conduct a parametric analysis to evaluate the sensitivities of wall-modeled large-eddy simulation (LES) with respect to subgrid-scale (SGS) models, mesh resolution, wall boundary conditions and mesh anisotropy. While such investigations have been conducted for attached/flat-plate flow configurations, systematic studies specifically targeting turbulent flows with separation are notably sparse. To bridge this gap, our study focuses on the flow over a two-dimensional Gaussian-shaped bump at a moderately high Reynolds number, which involves smooth-body separation of a turbulent boundary layer under pressure-gradient and surface- curvature effects. In the simulations, the no-slip condition at the wall is replaced by three different forms of boundary condition based on the thin boundary layer equations and the mean wall-shear stress from high-fidelity numerical simulation to avoid the additional complexity of modeling the wall-shear stress. Various statistics, including the mean separation bubble size, mean velocity profile, and dissipation from SGS model, are compared and analyzed. The results reveal that capturing the separation bubble strongly depends on the choice of SGS model. While simulations approach grid convergence with resolutions nearing those of wall-resolved LES meshes, above this limit, the LES predictions exhibit intricate sensitivities to mesh resolution. Furthermore, both wall boundary conditions and the anisotropy of mesh cells exert discernible impacts on the turbulent flow predictions, yet the magnitudes of these impacts vary based on the specific SGS model chosen for the simulation.more » « less
