scholarly journals Statistical Analysis of Dynamic Subgrid Modeling Approaches in Large Eddy Simulation

Aerospace ◽  
2021 ◽  
Vol 8 (12) ◽  
pp. 375
Author(s):  
Mohammad Khalid Hossen ◽  
Asokan Mulayath Variyath ◽  
Jahrul M. Alam

In large eddy simulation (LES) of turbulent flows, dynamic subgrid models would account for an average cascade of kinetic energy from the largest to the smallest scales of the flow. Yet, it is unclear which of the most critical dynamical processes can ensure the criterion mentioned above. Furthermore, evidence of vortex stretching being the primary mechanism of the cascade is not out of the question. In this article, we study essential statistical characteristics of vortex stretching. Our numerical results demonstrate that vortex stretching rate provides the energy dissipation rate necessary for modeling subgrid-scale turbulence. We have compared the interaction of subgrid stresses with the filtered quantities among four models using invariants of the velocity gradient tensor. The individual and the joint probability of vortex stretching and strain amplification show that vortex stretching rate is highly correlated with the energy cascade rate. Sheet-like flow structures are correlated with viscous dissipation, and vortex tubes are more stretched than compressed. The overall results indicate that the stretching mechanism extracts energy from the large-scale straining motion and passes it onto small-scale stretched vortices.

Author(s):  
Mohammad Khalid Hossen ◽  
Asokan Mulayath Variyath ◽  
Jahrul M Alam

In large eddy simulation (LES) of turbulent flows, the most critical dynamical processes to be considered by dynamic subgrid models to account for an average cascade of kinetic energy from the largest to the smallest scales of the flow is not fully clear. Furthermore, evidence of vortex stretching being the primary mechanism of the cascade is not out of the question. In this article, we study some essential statistical characteristics of vortex stretching and its role in dynamic approaches of modeling subgrid-scale turbulence. We have compared the interaction of subgrid stresses with the filtered quantities among four models using invariants of the velocity gradient tensor. This technique is a single unified approach to studying a wide range of length scales in the turbulent flow. In addition, it also provides a rational basis for the statistical characteristics a subgrid model must serve in physical space to ensure an appropriate cascade of kinetic energy. Results indicate that the stretching mechanism extracts energy from the large-scale straining motion and passes it onto small-scale stretched vortices.


1992 ◽  
Vol 238 ◽  
pp. 325-336 ◽  
Author(s):  
M. Germano

Explicit or implicit filtered representations of chaotic fields like spectral cut-offs or numerical discretizations are commonly used in the study of turbulence and particularly in the so-called large-eddy simulations. Peculiar to these representations is that they are produced by different filtering operators at different levels of resolution, and they can be hierarchically organized in terms of a characteristic parameter like a grid length or a spectral truncation mode. Unfortunately, in the case of a general implicit or explicit filtering operator the Reynolds rules of the mean are no longer valid, and the classical analysis of the turbulence in terms of mean values and fluctuations is not so simple.In this paper a new operatorial approach to the study of turbulence based on the general algebraic properties of the filtered representations of a turbulence field at different levels is presented. The main results of this analysis are the averaging invariance of the filtered Navier—Stokes equations in terms of the generalized central moments, and an algebraic identity that relates the turbulent stresses at different levels. The statistical approach uses the idea of a decomposition in mean values and fluctuations, and the original turbulent field is seen as the sum of different contributions. On the other hand this operatorial approach is based on the comparison of different representations of the turbulent field at different levels, and, in the opinion of the author, it is particularly fitted to study the similarity between the turbulence at different filtering levels. The best field of application of this approach is the numerical large-eddy simulation of turbulent flows where the large scale of the turbulent field is captured and the residual small scale is modelled. It is natural to define and to extract from the resolved field the resolved turbulence and to use the information that it contains to adapt the subgrid model to the real turbulent field. Following these ideas the application of this approach to the large-eddy simulation of the turbulent flow has been produced (Germano et al. 1991). It consists in a dynamic subgrid-scale eddy viscosity model that samples the resolved scale and uses this information to adjust locally the Smagorinsky constant to the local turbulence.


2001 ◽  
Vol 436 ◽  
pp. 353-391 ◽  
Author(s):  
J. C. R. HUNT ◽  
N. D. SANDHAM ◽  
J. C. VASSILICOS ◽  
B. E. LAUNDER ◽  
P. A. MONKEWITZ ◽  
...  

Recent research is making progress in framing more precisely the basic dynamical and statistical questions about turbulence and in answering them. It is helping both to define the likely limits to current methods for modelling industrial and environmental turbulent flows, and to suggest new approaches to overcome these limitations. Our selective review is based on the themes and new results that emerged from more than 300 presentations during the Programme held in 1999 at the Isaac Newton Institute, Cambridge, UK, and on research reported elsewhere. A general conclusion is that, although turbulence is not a universal state of nature, there are certain statistical measures and kinematic features of the small-scale flow field that occur in most turbulent flows, while the large-scale eddy motions have qualitative similarities within particular types of turbulence defined by the mean flow, initial or boundary conditions, and in some cases, the range of Reynolds numbers involved. The forced transition to turbulence of laminar flows caused by strong external disturbances was shown to be highly dependent on their amplitude, location, and the type of flow. Global and elliptical instabilities explain much of the three-dimensional and sudden nature of the transition phenomena. A review of experimental results shows how the structure of turbulence, especially in shear flows, continues to change as the Reynolds number of the turbulence increases well above about 104 in ways that current numerical simulations cannot reproduce. Studies of the dynamics of small eddy structures and their mutual interactions indicate that there is a set of characteristic mechanisms in which vortices develop (vortex stretching, roll-up of instability sheets, formation of vortex tubes) and another set in which they break up (through instabilities and self- destructive interactions). Numerical simulations and theoretical arguments suggest that these often occur sequentially in randomly occurring cycles. The factors that determine the overall spectrum of turbulence were reviewed. For a narrow distribution of eddy scales, the form of the spectrum can be defined by characteristic forms of individual eddies. However, if the distribution covers a wide range of scales (as in elongated eddies in the ‘wall’ layer of turbulent boundary layers), they collectively determine the spectra (as assumed in classical theory). Mathematical analyses of the Navier–Stokes and Euler equations applied to eddy structures lead to certain limits being defined regarding the tendencies of the vorticity field to become infinitely large locally. Approximate solutions for eigen modes and Fourier components reveal striking features of the temporal, near-wall structure such as bursting, and of the very elongated, spatial spectra of sheared inhomogeneous turbulence; but other kinds of eddy concepts are needed in less structured parts of the turbulence. Renormalized perturbation methods can now calculate consistently, and in good agreement with experiment, the evolution of second- and third-order spectra of homogeneous and isotropic turbulence. The fact that these calculations do not explicitly include high-order moments and extreme events, suggests that they may play a minor role in the basic dynamics. New methods of approximate numerical simulations of the larger scales of turbulence or ‘very large eddy simulation’ (VLES) based on using statistical models for the smaller scales (as is common in meteorological modelling) enable some turbulent flows with a non-local and non-equilibrium structure, such as impinging or convective flows, to be calculated more efficiently than by using large eddy simulation (LES), and more accurately than by using ‘engineering’ models for statistics at a single point. Generally it is shown that where the turbulence in a fluid volume is changing rapidly and is very inhomogeneous there are flows where even the most complex ‘engineering’ Reynolds stress transport models are only satisfactory with some special adaptation; this may entail the use of transport equations for the third moments or non-universal modelling methods designed explicitly for particular types of flow. LES methods may also need flow-specific corrections for accurate modelling of different types of very high Reynolds number turbulent flow including those near rigid surfaces.This paper is dedicated to the memory of George Batchelor who was the inspiration of so much research in turbulence and who died on 30th March 2000. These results were presented at the last fluid mechanics seminar in DAMTP Cambridge that he attended in November 1999.


2010 ◽  
Vol 132 (2) ◽  
Author(s):  
F. Abbasian ◽  
S. D. Yu ◽  
J. Cao

Large eddy simulation (LES) is employed in this paper to model the axial flow along a circular array of rods with a focus on anisotropic large-scale turbulence. The circular array consists of four whole rods and eight half rods, with a pitch-to-diameter ratio of 1.08. A dynamic Smagorinsky model with SIMPLE coupling method and a bounded central difference scheme are used to reduce numerical errors. The high demands for computations of the three-dimensional turbulent flows are afforded through parallel processing and utilization of 20 processors. The numerical results obtained using LES are compared with independent experimental data available in the literature; good agreement is achieved. The LES model was developed to accurately predict (i) the dependence of turbulence intensity and dominant frequency on the gap size and (ii) the turbulence structure in different directions.


Author(s):  
A. RINOSHIKA ◽  
Y. ZHENG ◽  
E. SHISHIDO

The three-dimensional orthogonal wavelet multi-resolution technique was applied to analyze flow structures of various scales around an externally mounted vehicle mirror. Firstly, the three-dimensional flow of mirror wake was numerically analyzed at a Reynolds number of 105 by using the large-eddy simulation (LES). Then the instantaneous velocity and vorticity were decomposed into the large-, intermediate- and relatively small-scale components by the wavelet multi-resolution technique. It was found that a three-dimensional large-scale vertical vortex dominates the mirror wake flow and makes a main contribution to vorticity concentration. Some intermediate- and relatively small-scale vortices were extracted from the LES and were clearly identifiable.


2018 ◽  
Vol 855 ◽  
pp. 371-407 ◽  
Author(s):  
W. Cheng ◽  
D. I. Pullin ◽  
R. Samtaney

We present wall-resolved large-eddy simulation (LES) of flow with free-stream velocity $\boldsymbol{U}_{\infty }$ over a cylinder of diameter $D$ rotating at constant angular velocity $\unicode[STIX]{x1D6FA}$ , with the focus on the lift crisis, which takes place at relatively high Reynolds number $Re_{D}=U_{\infty }D/\unicode[STIX]{x1D708}$ , where $\unicode[STIX]{x1D708}$ is the kinematic viscosity of the fluid. Two sets of LES are performed within the ( $Re_{D}$ , $\unicode[STIX]{x1D6FC}$ )-plane with $\unicode[STIX]{x1D6FC}=\unicode[STIX]{x1D6FA}D/(2U_{\infty })$ the dimensionless cylinder rotation speed. One set, at $Re_{D}=5000$ , is used as a reference flow and does not exhibit a lift crisis. Our main LES varies $\unicode[STIX]{x1D6FC}$ in $0\leqslant \unicode[STIX]{x1D6FC}\leqslant 2.0$ at fixed $Re_{D}=6\times 10^{4}$ . For $\unicode[STIX]{x1D6FC}$ in the range $\unicode[STIX]{x1D6FC}=0.48{-}0.6$ we find a lift crisis. This range is in agreement with experiment although the LES shows a deeper local minimum in the lift coefficient than the measured value. Diagnostics that include instantaneous surface portraits of the surface skin-friction vector field $\boldsymbol{C}_{\boldsymbol{f}}$ , spanwise-averaged flow-streamline plots, and a statistical analysis of local, near-surface flow reversal show that, on the leeward-bottom cylinder surface, the flow experiences large-scale reorganization as $\unicode[STIX]{x1D6FC}$ increases through the lift crisis. At $\unicode[STIX]{x1D6FC}=0.48$ the primary-flow features comprise a shear layer separating from that side of the cylinder that moves with the free stream and a pattern of oscillatory but largely attached flow zones surrounded by scattered patches of local flow separation/reattachment on the lee and underside of the cylinder surface. Large-scale, unsteady vortex shedding is observed. At $\unicode[STIX]{x1D6FC}=0.6$ the flow has transitioned to a more ordered state where the small-scale separation/reattachment cells concentrate into a relatively narrow zone with largely attached flow elsewhere. This induces a low-pressure region which produces a sudden decrease in lift and hence the lift crisis. Through this process, the boundary layer does not show classical turbulence behaviour. As $\unicode[STIX]{x1D6FC}$ is further increased at constant $Re_{D}$ , the localized separation zone dissipates with corresponding attached flow on most of the cylinder surface. The lift coefficient then resumes its increasing trend. A logarithmic region is found within the boundary layer at $\unicode[STIX]{x1D6FC}=1.0$ .


2002 ◽  
Vol 124 (4) ◽  
pp. 892-903 ◽  
Author(s):  
L. Persson ◽  
C. Fureby ◽  
N. Svanstedt

The ability to predict complex engineering flows is limited by the available turbulence models and the present-day computer capacity. In Reynolds averaged numerical simulations (RANS), which is the most prevalent approach today, equations for the mean flow are solved in conjunction with a model for the statistical properties of the turbulence. Considering the limitations of RANS and the desire to study more complex flows, more sophisticated methods are called for. An approach that fulfills these requirements is large-eddy simulation (LES) which attempts to resolve the dynamics of the large-scale flow, while modeling only the effects of the small-scale fluctuations. The limitations of LES are, however, closely tied to the subgrid model, which invariably relies on the use of eddy-viscosity models. Turbulent flows of practical importance involve inherently three-dimensional unsteady features, often subjected to strong inhomogeneous effects and rapid deformation that cannot be captured by isotropic models. As an alternative to the filtering approach fundamental to LES, we here consider the homogenization method, which consists of finding a so-called homogenized problem, i.e. finding a homogeneous “material” whose overall response is close to that of the heterogeneous “material” when the size of the inhomogeneity is small. Here, we develop a homogenization-based LES-model using a multiple-scales expansion technique and taking advantage of the scaling properties of the Navier-Stokes equations. To study the model simulations of forced homogeneous isotropic turbulence and channel flow are carried out, and comparisons are made with LES, direct numerical simulation and experimental data.


2008 ◽  
Vol 614 ◽  
pp. 207-252 ◽  
Author(s):  
M. BINI ◽  
W. P. JONES

A large-eddy-based methodology for the simulation of turbulent sprays is discussed. The transport equations for the spatially filtered gas phase variables, in which source terms accounting for the droplet effects are added, are solved together with a probabilistic description of the liquid phase. The probabilistic approach for the liquid phase is based on the transport equation for the spatially filtered joint probability density function of the variables required in order to describe the state of the liquid phase. In this equation, unclosed terms representing the filtered Lagrangian rates of change of the variables describing the spray are present. General modelling ideas for subgrid-scale (SGS) effects are proposed. The capabilities of the approach and the validity of the closure models, with particular with respect to the SGS dispersion, are investigated through application to a dilute particle-laden turbulent mixing layer. It is demonstrated that the formulation is able to reproduce very closely the measured properties of both the continuous and dispersed phases. The large-eddy simulation (LES) results are also found to be entirely consistent with the experimentally observed characteristics of droplet–gas turbulence interactions. Consistent with direct numerical simulation (DNS) studies of isotropic turbulence laden with particles where the entire turbulence spectrum is found to be modulated by the presence of particles, the present investigation, which comprises the effects of particle transport upon the large-scale vortical structures of a turbulent shear flow, highlights what appears to be a selective behaviour; few large-scale frequencies gain energy whereas the remaining modes are damped.


2007 ◽  
Vol 64 (12) ◽  
pp. 4445-4457 ◽  
Author(s):  
M. Antonelli ◽  
R. Rotunno

Abstract This paper describes results from a large-eddy simulation (LES) model used in an idealized setting to simulate the onset of the sea breeze. As the LES is capable of simulating boundary layer–scale, three-dimensional turbulence along with the mesoscale sea-breeze circulation, a parameterization of the planetary boundary layer was unnecessary. The basic experimental design considers a rotating, uniformly stratified, resting atmosphere that is suddenly heated at the surface over the “land” half of the domain. To focus on the simplest nontrivial problem, the diurnal cycle, effects of moisture, interactions with large-scale winds, and coastline curvature were all neglected in this study. The assumption of a straight coastline allows the use of a rectangular computational domain that extends to 50 km on either side of the coast, but only 5 km along the coast, with 100-m grid intervals so that the small-scale turbulent convective eddies together with the mesoscale sea breeze may be accurately computed. Through dimensional analysis of the simulation results, the length and velocity scales characterizing the simulated sea breeze as functions of the externally specified parameters are identified.


Sign in / Sign up

Export Citation Format

Share Document