scholarly journals Higher-order dissipation in the theory of homogeneous isotropic turbulence

2016 ◽  
Vol 803 ◽  
pp. 250-274 ◽  
Author(s):  
Norbert Peters ◽  
Jonas Boschung ◽  
Michael Gauding ◽  
Jens Henrik Goebbert ◽  
Reginald J. Hill ◽  
...  

The two-point theory of homogeneous isotropic turbulence is extended to source terms appearing in the equations for higher-order structure functions. For this, transport equations for these source terms are derived. We focus on the trace of the resulting equations, which is of particular interest because it is invariant and therefore independent of the coordinate system. In the trace of the even-order source term equation, we discover the higher-order moments of the dissipation distribution, and the individual even-order source term equations contain the higher-order moments of the longitudinal, transverse and mixed dissipation distribution functions. This shows for the first time that dissipation fluctuations, on which most of the phenomenological intermittency models are based, are contained in the Navier–Stokes equations. Noticeably, we also find the volume-averaged dissipation $\unicode[STIX]{x1D700}_{r}$ used by Kolmogorov (J. Fluid Mech., vol. 13, 1962, pp. 82–85) in the resulting system of equations, because it is related to dissipation correlations.

2019 ◽  
Vol 864 ◽  
pp. 244-272 ◽  
Author(s):  
L. Djenidi ◽  
R. A. Antonia ◽  
S. L. Tang

The problem of homogeneous isotropic turbulence (HIT) is revisited within the analytical framework of the Navier–Stokes equations, with a view to assessing rigorously the consequences of the scale invariance (an exact property of the Navier–Stokes equations) for any Reynolds number. The analytical development, which is independent of the 1941 (K41) and 1962 (K62) theories of Kolmogorov for HIT for infinitely large Reynolds number, is applied to the transport equations for the second- and third-order moments of the longitudinal velocity increment, $(\unicode[STIX]{x1D6FF}u)$. Once the normalised equations and the constraints required for complying with the scale-invariance property of the equations are presented, results derived from these equations and constraints are discussed and compared with measurements. It is found that the fluid viscosity, $\unicode[STIX]{x1D708}$, and the mean kinetic energy dissipation rate, $\overline{\unicode[STIX]{x1D716}}$ (the overbar denotes spatial and/or temporal averaging), are the only scaling parameters that make the equations scale-invariant. The analysis further leads to expressions for the distributions of the skewness and the flatness factor of $(\unicode[STIX]{x1D6FF}u)$ and shows that these distributions must exhibit plateaus (of different magnitudes) in the dissipative and inertial ranges, as the Taylor microscale Reynolds number $Re_{\unicode[STIX]{x1D706}}$ increases indefinitely. Also, the skewness and flatness factor of the longitudinal velocity derivative become constant as $Re_{\unicode[STIX]{x1D706}}$ increases; this is supported by experimental data. Further, the analysis, backed up by experimental evidence, shows that, beyond the dissipative range, the behaviour of $\overline{(\unicode[STIX]{x1D6FF}u)^{n}}$ with $n=2$, 3 and 4 cannot be represented by a power law of the form $r^{\unicode[STIX]{x1D701}_{n}}$ when the Reynolds number is finite. It is shown that only when $Re_{\unicode[STIX]{x1D706}}\rightarrow \infty$ can an $n$-thirds law (i.e. $\overline{(\unicode[STIX]{x1D6FF}u)^{n}}\sim r^{\unicode[STIX]{x1D701}_{n}}$, with $\unicode[STIX]{x1D701}_{n}=n/3$) emerge, which is consistent with the onset of a scaling range.


2012 ◽  
Vol 700 ◽  
pp. 1-4 ◽  
Author(s):  
R. M. Kerr

AbstractSince the advent of cluster computing over 10 years ago there has been a steady output of new and better direct numerical simulation of homogeneous, isotropic turbulence with spectra and lower-order statistics converging to experiments and many phenomenological models. The next step is to directly compare these simulations to new models and new mathematics, employing the simulated data sets in novel ways, especially when experimental results do not exist or are poorly converged. For example, many of the higher-order moments predicted by the models converge slowly in experiments. The solution with a simulation is to do what an experiment cannot. The calculation and analysis of Yeung, Donzis & Sreenivasan (J. Fluid Mech., this issue, vol. 700, 2012, pp. 5–15) represents the vanguard of new simulations and new numerical analysis that will fill this gap. Where individual higher-order moments of the vorticity squared (the enstrophy) and kinetic energy dissipation might be converging slowly, they have focused upon ratios between different moments that have better convergence properties. This allows them to more fully explore the statistical distributions that eventually must be modelled. This approach is consistent with recent mathematics that focuses upon temporal intermittency rather than spatial intermittency. The principle is that when the flow is nearly singular, during ‘bad’ phases, when global properties can go up and down by many orders of magnitude, if appropriate ratios are taken, convergence rates should improve. Furthermore, in future analysis it might be possible to use these ratios to gain new insights into the intermittency and regularity properties of the underlying equations.


2021 ◽  
Vol 931 ◽  
Author(s):  
Fujihiro Hamba

The energy spectrum is commonly used to describe the scale dependence of turbulent fluctuations in homogeneous isotropic turbulence. In contrast, one-point statistical quantities, such as the turbulent kinetic energy, are mainly employed for inhomogeneous turbulence models. Attempts have been made to describe the scale dependence of inhomogeneous turbulence using the second-order structure function and two-point velocity correlation. However, unlike the energy spectrum, expressions for the energy density in the scale space fail to satisfy the requirement of being non-negative. In this study, a new expression for the scale-space energy density based on filtered velocities is proposed to clarify the reasons behind the negative values of the energy density and to obtain a better understanding of inhomogeneous turbulence. The new expression consists of homogeneous and inhomogeneous parts; the former is always non-negative, while the latter can be negative because of the turbulence inhomogeneity. Direct numerical simulation data of homogeneous isotropic turbulence and a turbulent channel flow are used to evaluate the two parts of the energy density and turbulent energy. It was found that the inhomogeneous part of the turbulent energy shows non-zero values near the wall and at the centre of a channel flow. In particular, the inhomogeneous part of the energy density changes its sign depending on the scale. A concave profile of the filtered-velocity variance at the wall accounts for the negative value of the energy density in the region very close to the wall.


2017 ◽  
Vol 833 ◽  
pp. 717-744 ◽  
Author(s):  
Shaowu Pan ◽  
Eric Johnsen

While Stokes’ hypothesis of neglecting bulk viscous effects is exact for monatomic gases and unlikely to strongly affect the dynamics of fluids whose bulk-to-shear viscosity ratio is small and/or of weakly compressible turbulence, it is unclear to what extent this assumption holds for compressible, turbulent flows of gases whose bulk viscosity is orders of magnitude larger than their shear viscosities (e.g. $\text{CO}_{2}$). Our objective is to understand the mechanisms by which bulk viscosity and the associated phenomena affect moderately compressible turbulence, in particular energy transfer and dissipation. Using direct numerical simulations of the compressible Navier–Stokes equations, we study the decay of compressible, homogeneous, isotropic turbulence for ratios of bulk-to-shear viscosity ranging from 0 to 1000. Our simulations demonstrate that bulk viscosity increases the decay rate of turbulent kinetic energy; whereas enstrophy exhibits little sensitivity to bulk viscosity, dilatation is reduced by over two orders of magnitude within the first two eddy-turnover times. Via a Helmholtz decomposition of the flow, we determine that the primary action of bulk viscosity is to damp the dilatational velocity fluctuations and reduce dilatational–solenoidal exchanges, as well as pressure–dilatation coupling. In short, bulk viscosity renders compressible turbulence incompressible by reducing energy transfer between translational and internal degrees of freedom. Our results indicate that for gases whose bulk viscosity is of the order of their shear viscosity (e.g. hydrogen) the turbulence is not significantly affected by bulk viscous dissipation, in which case neglecting bulk viscosity is acceptable in practice. However, in problems involving compressible, turbulent flows of gases like $\text{CO}_{2}$ whose bulk viscosities are thousands of times greater than their shear viscosities, bulk viscosity cannot be ignored.


2012 ◽  
Vol 707 ◽  
pp. 74-110 ◽  
Author(s):  
Krishnendu Sinha

AbstractInteraction of turbulent fluctuations with a shock wave plays an important role in many high-speed flow applications. This paper studies the amplification of enstrophy, defined as mean-square fluctuating vorticity, in homogeneous isotropic turbulence passing through a normal shock. Linearized Navier–Stokes equations written in a frame of reference attached to the unsteady shock wave are used to derive transport equations for the vorticity components. These are combined to obtain an equation that describes the evolution of enstrophy across a time-averaged shock wave. A budget of the enstrophy equation computed using results from linear interaction analysis and data from direct numerical simulations identifies the dominant physical mechanisms in the flow. Production due to mean flow compression and baroclinic torques are found to be the major contributors to the enstrophy amplification. Closure approximations are proposed for the unclosed correlations in the production and baroclinic source terms. The resulting model equation is integrated to obtain the enstrophy jump across a shock for a range of upstream Mach numbers. The model predictions are compared with linear theory results for varying levels of vortical and entropic fluctuations in the upstream flow. The enstrophy model is then cast in the form of$k$–$\epsilon $equations and used to compute the interaction of homogeneous isotropic turbulence with normal shocks. The results are compared with available data from direct numerical simulations. The equations are further used to propose a model for the amplification of turbulent viscosity across a shock, which is then applied to a canonical shock–boundary layer interaction. It is shown that the current model is a significant improvement over existing models, both for homogeneous isotropic turbulence and in the case of complex high-speed flows with shock waves.


2014 ◽  
Vol 31 (10) ◽  
pp. 2131-2144 ◽  
Author(s):  
Robert C. Jackson ◽  
Greg M. McFarquhar

Abstract A recent study showed that the ratio of the number of distribution functions derived from 2D cloud probes (2DCs) with standard tips to those with antishatter tips used during the 2008 Indirect and Semidirect Aerosol Campaign (ISDAC) and Instrumentation Development and Education in Airborne Science 2011 (IDEAS-2011) was greater than 1 for ice crystals with maximum dimension D < 500 μm. To assess the applicability of 2DC data obtained without antishatter tips previously used in parameterization schemes for numerical models and remote sensing retrievals, the impacts of artifacts on bulk microphysical and scattering properties were examined by quantifying differences between such properties derived from 2DCs with standard and antishatter tips, and with and without the use of shatter detection algorithms using the ISDAC and IDEAS-2011 data. Using either modified tips or algorithms changed the quantities dominated by higher-order moments, such as ice water content, bulk extinction, effective radius, mass-weighted terminal velocity, median mass diameter, asymmetry parameter, and single-scatter albedo, at wavenumbers from 5 to 100 cm−1 and wavelengths of 0.5–5 μm by less than 20%. This is significantly less than the fractional changes quantities dominated by lower-order moments, such as number concentration. The results suggest that model parameterizations and remote sensing techniques based on higher-order moments of ice particle size distributions obtained in conditions similar to those sampled during IDEAS-2011 and ISDAC derived from 2DCs are not substantially biased by shattered remnants.


Sign in / Sign up

Export Citation Format

Share Document