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

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.

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.


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.


Author(s):  
B. Hu ◽  
S. Banerjee ◽  
K. Liu ◽  
D. Rajamohan ◽  
J. M. Deur ◽  
...  

We performed Large Eddy Simulation (LES) of a turbulent non-reacting n-Heptane spray jet, referred to as Spray H in the Engine Combustion Network (ECN), and executed a data analysis focused on key LES metrics such as fraction of resolved turbulent kinetic energy and similarity index. In the simulation, we used the dynamic structure model for the sub-grid stress, and the Lagrangian-based spray-parcel models coupled with the blob-injection model. The finest mesh-cell size used was characterized by an Adaptive Mesh Refinement (AMR) cell size of 0.0625 mm. To obtain ensemble statistics, we performed 28 numerical realizations of the simulation. Demonstrated by the comparison with experimental data in a previous study [7], this LES has accurately predicted global quantities, such as liquid and vapor penetrations. The analysis in this work shows that 14 realizations of LES are sufficient to provide a reasonable representation of the average flow behavior that is benchmarked against the 28-realization ensemble. With the current mesh, numerical schemes, and sub-grid scale turbulence model, more than 95% of the turbulent kinetic energy is directly resolved in the flow regions of interest. The large-scale flow structures inferred from a statistical analysis reveal a region of disorganized flow around the peripheral region of the spray jet, which appears to be linked to the entrainment process.


Author(s):  
Salar Taghizadeh ◽  
Sumanta Acharya ◽  
Kong Ling ◽  
Yousef Kanani ◽  
Xuan Ge

This study presents a transient three-dimensional numerical study on fluid flow and heat transfer of flat-tube array using large eddy simulation (LES) covering both laminar and turbulent flow regimes. The simulations were performed in a rectangular region containing only one tube with periodic conditions specified on all boundaries. A staggered flat-plate array was first studied, and an existing solution was used for validation purpose. The numerical models were then applied to an in-line array composed of flat tubes with an aspect ratio of 0.25 and fixed tube spacings. By varying the in-flow velocity, the tube array was studied over a wide range of Reynolds number (600–12000). Temperature, velocity, and turbulent kinetic energy distributions as well as the interactions between them are presented and analyzed. Furthermore, the local heat transfer rate was analyzed along the various parts of the tube (leading edge, flat-top and wake or trailing-edge regions). Heat transfer correlation for each region of the tube and the entire tube array is proposed.


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.


2015 ◽  
Vol 143 (3) ◽  
pp. 828-844 ◽  
Author(s):  
Jerôme Schalkwijk ◽  
Harmen J. J. Jonker ◽  
A. Pier Siebesma ◽  
Fred C. Bosveld

Abstract Results are presented of two large-eddy simulation (LES) runs of the entire year 2012 centered at the Cabauw observational supersite in the Netherlands. The LES is coupled to a regional weather model that provides the large-scale information. The simulations provide three-dimensional continuous time series of LES-generated turbulence and clouds, which can be compared in detail to the extensive observational dataset of Cabauw. The LES dataset is available from the authors on request. This type of LES setup has a number of advantages. First, it can provide a more statistical approach to the study of turbulent atmospheric flow than the more common case studies, since a diverse but representative set of conditions is covered, including numerous transitions. This has advantages in the design and evaluation of parameterizations. Second, the setup can provide valuable information on the quality of the LES model when applied to such a wide range of conditions. Last, it also provides the possibility to emulate observation techniques. This might help detect limitations and potential problems of a variety of measurement techniques. The LES runs are validated through a comparison with observations from the observational supersite and with results from the “parent” large-scale model. The long time series that are generated, in combination with information on the spatial structure, provide a novel opportunity to study time scales ranging from seconds to seasons. This facilitates a study of the power spectrum of horizontal and vertical wind speed variance to identify the dominant variance-containing time scales.


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$ .


Sign in / Sign up

Export Citation Format

Share Document