Implicit Large Eddy Simulation of Two-Dimensional Homogeneous Turbulence Using Weighted Compact Nonlinear Scheme

2009 ◽  
Vol 131 (6) ◽  
Author(s):  
Keiichi Ishiko ◽  
Naofumi Ohnishi ◽  
Kazuyuki Ueno ◽  
Keisuke Sawada
Keyword(s):  
Energy Spectrum ◽  
Inertial Range ◽  
Homogeneous Turbulence ◽  
Two Dimensional ◽  
Third Order ◽  
The Third ◽  
Eddy Simulation ◽  
Numerical Viscosity ◽  
Large Eddy ◽  
Implicit Large Eddy Simulation

For the aim of computing compressible turbulent flowfield involving shock waves, an implicit large eddy simulation (LES) code has been developed based on the idea of monotonically integrated LES. We employ the weighted compact nonlinear scheme (WCNS) not only for capturing possible shock waves but also for attaining highly accurate resolution required for implicit LES. In order to show that WCNS is a proper choice for implicit LES, a two-dimensional homogeneous turbulence is first obtained by solving the Navier–Stokes equations for incompressible flow. We compare the inertial range in the computed energy spectrum with that obtained by the direct numerical simulation (DNS) and also those given by the different LES approaches. We then obtain the same homogeneous turbulence by solving the equations for compressible flow. It is shown that the present implicit LES can reproduce the inertial range in the energy spectrum given by DNS fairly well. A truncation of energy spectrum occurs naturally at high wavenumber limit indicating that dissipative effect is included properly in the present approach. A linear stability analysis for WCNS indicates that the third order interpolation determined in the upwind stencil introduces a large amount of numerical viscosity to stabilize the scheme, but the same interpolation makes the scheme weakly unstable for waves satisfying kΔx≈1. This weak instability results in a slight increase in the energy spectrum at high wavenumber limit. In the computed result of homogeneous turbulence, a fair correlation is shown to exist between the locations where the magnitude of ∇×ω becomes large and where the weighted combination of the third order interpolations in WCNS deviates from the optimum ratio to increase the amount of numerical viscosity. Therefore, the numerical viscosity involved in WCNS becomes large only at the locations where the subgrid-scale viscosity can arise in ordinary LES. This suggests the reason why the present implicit LES code using WCNS can resolve turbulent flowfield reasonably well.

Can the atmospheric kinetic energy spectrum be explained by two-dimensional turbulence?

1999 ◽  
Vol 388 ◽  
pp. 259-288 ◽  
Author(s):  
ERIK LINDBORG
Keyword(s):  
Energy Spectrum ◽  
Structure Function ◽  
Structure Functions ◽  
Separation Distance ◽  
Inertial Range ◽  
Order Structure ◽  
Two Dimensional ◽  
Third Order ◽  
Kinetic Energy Spectrum ◽  
The Third

The statistical features of turbulence can be studied either through spectral quantities, such as the kinetic energy spectrum, or through structure functions, which are statistical moments of the difference between velocities at two points separated by a variable distance. In this paper structure function relations for two-dimensional turbulence are derived and compared with calculations based on wind data from 5754 airplane flights, reported in the MOZAIC data set. For the third-order structure function two relations are derived, showing that this function is generally positive in the two-dimensional case, contrary to the three-dimensional case. In the energy inertial range the third-order structure function grows linearly with separation distance and in the enstrophy inertial range it grows cubically with separation distance. A Fourier analysis shows that the linear growth is a reflection of a constant negative spectral energy flux, and the cubic growth is a reflection of a constant positive spectral enstrophy flux. Various relations between second-order structure functions and spectral quantities are also derived. The measured second-order structure functions can be divided into two different types of terms, one of the form r2/3, giving a k−5/3-range and another, including a logarithmic dependence, giving a k−3-range in the energy spectrum. The structure functions agree better with the two-dimensional isotropic relation for larger separations than for smaller separations. The flatness factor is found to grow very fast for separations of the order of some kilometres. The third-order structure function is accurately measured in the interval [30, 300] km and is found to be positive. The average enstrophy flux is measured as Πω≈1.8×10−13 s−3 and the constant in the k−3-law is measured as [Kscr ]≈0.19. It is argued that the k−3-range can be explained by two-dimensional turbulence and can be interpreted as an enstrophy inertial range, while the k−5/3-range can probably not be explained by two-dimensional turbulence and should not be interpreted as a two-dimensional energy inertial range.

A note on Kolmogorov's third-order structure-function law, the local isotropy hypothesis and the pressure–velocity correlation

1996 ◽  
Vol 326 ◽  
pp. 343-356 ◽  
Author(s):  
Erik Lindborg
Keyword(s):  
Structure Function ◽  
Strain Tensor ◽  
Inertial Range ◽  
Order Structure ◽  
Two Dimensional ◽  
Mean Flow ◽  
Velocity Correlation ◽  
Third Order ◽  
The Third ◽  
Point Pressure

We show that Kolmogorov's (1941b) inertial-range law for the third-order structure function can be derived from a dynamical equation including pressure terms and mean flow gradient terms. A new inertial-range law, relating the two-point pressure–velocity correlation to the single-point pressure–strain tensor, is also derived. This law shows that the two-point pressure–velocity correlation, just like the third-order structure function, grows linearly with the separation distance in the inertial range. The physical meaning of both this law and Kolmogorov's law is illustrated by a Fourier analysis. An inertial-range law is also derived for the third-order velocity–enstrophy structure function of two-dimensional turbulence. It is suggested that the second-order vorticity structure function of two-dimensional turbulence is constant and scales with$\epsilon ^{2/3}_\omega$in the enstrophy inertial range, εωbeing the enstrophy dissipation. Owing to the constancy of this law, it does not imply a Fourier-space inertial-range law, and therefore it is not equivalent to thek−1law for the enstrophy spectrum, suggested by Kraichnan (1967) and Batchelor (1969).

Applicability of Kolmogorov's and Monin's equations of turbulence

1997 ◽  
Vol 353 ◽  
pp. 67-81 ◽  
Author(s):  
REGINALD J. HILL
Keyword(s):  
Structure Function ◽  
Structure Functions ◽  
Inertial Range ◽  
Homogeneous Turbulence ◽  
Order Structure ◽  
Third Order ◽  
Local Isotropy ◽  
The Third ◽  
Moderate Reynolds Number ◽  
Wide Range

The equation relating second- and third-order velocity structure functions was presented by Kolmogorov; Monin attempted to derive that equation on the basis of local isotropy. Recently, concerns have been raised to the effect that Kolmogorov's equation and an ancillary incompressibility condition governing the third-order structure function were proven only on the restrictive basis of isotropy and that the statistic involving pressure that appears in the derivation of Kolmogorov's equation might not vanish on the basis of local isotropy. These concerns are resolved. In so doing, results are obtained for the second- and third-order statistics on the basis of local homogeneity without use of local isotropy. These results are applicable to future studies of the approach toward local isotropy. Accuracy of Kolmogorov's equation is shown to be more sensitive to anisotropy of the third-order structure function than to anisotropy of the second-order structure function. Kolmogorov's 4/5 law for the inertial range of the third-order structure function is obtained without use of the incompressibility conditions on the second- and third-order structure functions. A generalization of Kolmogorov's 4/5 law, which applies to the inertial range of locally homogeneous turbulence at very large Reynolds numbers, is shown to also apply to the energy-containing range for the more restrictive case of stationary, homogeneous turbulence. The variety of derivations of Kolmogorov's and Monin's equations leads to a wide range of applicability to experimental conditions, including, in some cases, turbulence of moderate Reynolds number.

Large eddy simulation of aircraft wake vortices within homogeneous turbulence - Crow instability

AIAA Journal ◽  
2000 ◽  
Vol 38 ◽  
pp. 292-300 ◽  
Author(s):  
Jongil Han ◽  
Yuh-Lang Lin ◽  
David G. Schowalter ◽  
S. P. Arya ◽  
Fred H. Proctor
Keyword(s):  
Large Eddy Simulation ◽  
Homogeneous Turbulence ◽  
Wake Vortices ◽  
Eddy Simulation ◽  
Large Eddy ◽  
Aircraft Wake ◽  
Crow Instability

Large eddy simulation investigation of the canonical shock–turbulence interaction

2018 ◽  
Vol 858 ◽  
pp. 500-535 ◽  
Author(s):  
N. O. Braun ◽  
D. I. Pullin ◽  
D. I. Meiron
Keyword(s):  
Nonlinear Effects ◽  
Streamwise Velocity ◽  
Inertial Range ◽  
Eddy Simulation ◽  
Weakly Compressible ◽  
Turbulent Statistics ◽  
Large Eddy ◽  
Reynolds Number Effects ◽  
Fully Developed Turbulent Flow ◽  
Mach Numbers

High resolution large eddy simulations (LES) are performed to study the interaction of a stationary shock with fully developed turbulent flow. Turbulent statistics downstream of the interaction are provided for a range of weakly compressible upstream turbulent Mach numbers $M_{t}=0.03{-}0.18$, shock Mach numbers $M_{s}=1.2{-}3.0$ and Taylor-based Reynolds numbers $Re_{\unicode[STIX]{x1D706}}=20{-}2500$. The LES displays minimal Reynolds number effects once an inertial range has developed for $Re_{\unicode[STIX]{x1D706}}>100$. The inertial range scales of the turbulence are shown to quickly return to isotropy, and downstream of sufficiently strong shocks this process generates a net transfer of energy from transverse into streamwise velocity fluctuations. The streamwise shock displacements are shown to approximately follow a $k^{-11/3}$ decay with wavenumber as predicted by linear analysis. In conjunction with other statistics this suggests that the instantaneous interaction of the shock with the upstream turbulence proceeds in an approximately linear manner, but nonlinear effects immediately downstream of the shock significantly modify the flow even at the lowest considered turbulent Mach numbers.

Wall-Modelled Implicit Large-Eddy Simulation of the RA16SC1 Highlift Configuration

2013 ◽  
Author(s):  
Michael Meyer ◽  
Stefan Hickel ◽  
Christian Breitsamter ◽  
Nikolaus Adams
Keyword(s):  
Large Eddy Simulation ◽  
Eddy Simulation ◽  
Large Eddy ◽  
Implicit Large Eddy Simulation
Sign in / Sign up

Export Citation Format

Share Document