scholarly journals Dissipation rate of turbulent kinetic energy in stably stratified sheared flows

2019 ◽  
Vol 19 (4) ◽  
pp. 2489-2496 ◽  
Author(s):  
Sergej Zilitinkevich ◽  
Oleg Druzhinin ◽  
Andrey Glazunov ◽  
Evgeny Kadantsev ◽  
Evgeny Mortikov ◽  
...  

Abstract. Over the years, the problem of dissipation rate of turbulent kinetic energy (TKE) in stable stratification remained unclear because of the practical impossibility to directly measure the process of dissipation that takes place at the smallest scales of turbulent motion. Poor representation of dissipation causes intolerable uncertainties in turbulence-closure theory and thus in modelling stably stratified turbulent flows. We obtain a theoretical solution to this problem for the whole range of stratifications from neutral to limiting stable; and validate it via (i) direct numerical simulation (DNS) immediately detecting the dissipation rate and (ii) indirect estimates of dissipation rate retrieved via the TKE budget equation from atmospheric measurements of other components of the TKE budget. The proposed formulation of dissipation rate will be of use in any turbulence-closure models employing the TKE budget equation and in problems requiring precise knowledge of the high-frequency part of turbulence spectra in atmospheric chemistry, aerosol science, and microphysics of clouds.

2018 ◽  
Author(s):  
Sergej Zilitinkevich ◽  
Oleg Druzhinin ◽  
Andrey Glazunov ◽  
Evgeny Kadantsev ◽  
Evgeny Mortikov ◽  
...  

Abstract. Over the years, the problem of dissipation rate of turbulent kinetic energy (TKE) in stable stratification remained unclear because of the practical impossibility to directly measure the process of dissipation that takes place at the smallest scales of turbulent motions. Poor representation of dissipation causes intolerable uncertainties in turbulence-closure theory and, thus, in modelling stably stratified turbulent flows. We obtain theoretical solution to this problem for the whole range of stratifications from neutral to limiting stable; and validate it via (i) direct numerical simulation (DNS) immediately detecting the dissipation rate and (ii) indirect estimates of dissipation rate retrieved via the TKE-budget equation from atmospheric measurements of other components of the TKE-budget. The proposed formulation of dissipation rate will be of use in any turbulence-closure models employing the TKE budget equation and in problems requiring precise knowledge of the high-frequency part of turbulence spectra in atmospheric chemistry, aerosol science and microphysics of clouds.


2006 ◽  
Vol 23 (7) ◽  
pp. 964-976 ◽  
Author(s):  
Ren-Chieh Lien ◽  
Eric A. D'Asaro

Abstract This study tests the ability of a neutrally buoyant float to estimate the dissipation rate of turbulent kinetic energy ɛ from its vertical acceleration spectrum using an inertial subrange method. A Lagrangian float was equipped with a SonTek acoustic Doppler velocimeter (ADV), which measured the vector velocity 1 m below the float's center, and a pressure sensor, which measured the float's depth. Measurements were taken in flows where estimates of ɛ varied from 10−8 to 10−3 W kg−1. Previous observational and theoretical studies conclude that the Lagrangian acceleration spectrum is white within the inertial subrange with a level proportional to ɛ. The size of the Lagrangian float introduces a highly reproducible spectral attenuation at high frequencies. Estimates of the dissipation rate of turbulent kinetic energy using float measurements ɛfloat were obtained by fitting the observed spectra to a model spectrum that included the attenuation effect. The ADV velocity measurements were converted to a wavenumber spectrum using a variant of Taylor's hypothesis. The spectrum exhibited the expected −5/3 slope within an inertial subrange. The turbulent kinetic energy dissipation rate ɛADV was computed from the level of this spectrum. These two independent estimates, ɛADV and ɛfloat, were highly correlated. The ratio ɛfloat/ɛADV deviated from one by less than a factor of 2 over the five decades of ɛ measured. This analysis confirms that ɛ can be estimated reliably from Lagrangian float acceleration spectra in turbulent flows. For the meter-sized floats used here, the size of the float and the noise level of the pressure measurements sets a lower limit of ɛfloat > 10−8 W kg−1.


2019 ◽  
Vol 866 ◽  
pp. 897-928 ◽  
Author(s):  
P. Orlandi

Data available in the literature from direct numerical simulations of two-dimensional turbulent channels by Lee & Moser (J. Fluid Mech., vol. 774, 2015, pp. 395–415), Bernardini et al. (J. Fluid Mech., 742, 2014, pp. 171–191), Yamamoto & Tsuji (Phys. Rev. Fluids, vol. 3, 2018, 012062) and Orlandi et al. (J. Fluid Mech., 770, 2015, pp. 424–441) in a large range of Reynolds number have been used to find that $S^{\ast }$ the ratio between the eddy turnover time ($q^{2}/\unicode[STIX]{x1D716}$, with $q^{2}$ being twice the turbulent kinetic energy and $\unicode[STIX]{x1D716}$ the isotropic rate of dissipation) and the time scale of the mean deformation ($1/S$), scales very well with the Reynolds number in the wall region. The good scaling is due to the eddy turnover time, although the turbulent kinetic energy and the rate of isotropic dissipation show a Reynolds dependence near the wall; $S^{\ast }$, as well as $-\langle Q\rangle =\langle s_{ij}s_{ji}\rangle -\langle \unicode[STIX]{x1D714}_{i}\unicode[STIX]{x1D714}_{i}/2\rangle$ are linked to the flow structures, and also the latter quantity presents a good scaling near the wall. It has been found that the maximum of turbulent kinetic energy production $P_{k}$ occurs in the layer with $-\langle Q\rangle \approx 0$, that is, where the unstable sheet-like structures roll-up to become rods. The decomposition of $P_{k}$ in the contribution of elongational and compressive strain demonstrates that the two contributions present a good scaling. However, the good scaling holds when the wall and the outer structures are separated. The same statistics have been evaluated by direct simulations of turbulent flows in the presence of different types of corrugations on both walls. The flow physics in the layer near the plane of the crests is strongly linked to the shape of the surface and it has been demonstrated that the $u_{2}$ (normal to the wall) fluctuations are responsible for the modification of the flow structures, for the increase of the resistance and of the turbulent kinetic energy production.


2012 ◽  
Vol 7 (1) ◽  
pp. 53-69
Author(s):  
Vladimir Dulin ◽  
Yuriy Kozorezov ◽  
Dmitriy Markovich

The present paper reports PIV (Particle Image Velocimetry) measurements of turbulent velocity fluctuations statistics in development region of an axisymmetric free jet (Re = 28 000). To minimize measurement uncertainty, adaptive calibration, image processing and data post-processing algorithms were utilized. On the basis of theoretical analysis and direct measurements, the paper discusses effect of PIV spatial resolution on measured statistical characteristics of turbulent fluctuations. Underestimation of the second-order moments of velocity derivatives and of the turbulent kinetic energy dissipation rate due to a finite size of PIV interrogation area and finite thickness of laser sheet was analyzed from model spectra of turbulent velocity fluctuations. The results are in a good agreement with the measured experimental data. The paper also describes performance of possible ways to account for unresolved small-scale velocity fluctuations in PIV measurements of the dissipation rate. In particular, a turbulent viscosity model can be efficiently used to account for the unresolved pulsations in a free turbulent flow


2021 ◽  
pp. 079
Author(s):  
Ivan Bašták Ďurán ◽  
Pascal Marquet

Le schéma de turbulence Toucans est utilisé dans la configuration opérationnelle Alaro du modèle Aladin depuis début 2015. Son développement a été initié, guidé et en grande partie conçu par Jean-François Geleyn. Ce développement a commencé avec le prédécesseur du schéma Toucans, le schéma « pseudo-pronostique » en énergie cinétique turbulente, lui-même basé sur l'ancien schéma de turbulence de Louis, mais étendu dans Toucans à un schéma pronostique. Le schéma Toucans a pour objectif de traiter de manière cohérente les fonctions qui dépendent de la stabilité verticale de l'atmosphère, de l'influence de l'humidité et des échelles de longueur de la turbulence (de mélange et de dissipation). De plus, de nouvelles caractéristiques ont été ajoutées : une représentation améliorée pour les stratifications très stables (absence de nombre de Richardson critique), une meilleure représentation de l'anisotropie, un paramétrage unifié de la turbulence et des nuages par l'ajout d'une deuxième énergie turbulente pronostique et la paramétrisation des moments du troisième ordre. The Toucans turbulence scheme is a turbulence scheme that is used in the operational Alaro configuration of the Aladin model since early 2015. Its development was initiated, guided and to a large extend authored by Jean-François Geleyn. The development started with the predecessor of the Toucans scheme, the "pseudo-prognostic" turbulent kinetic energy scheme which itself was built on the "Louis" turbulence scheme, but extended to a prognostic scheme. The Toucans scheme aims for a consistent treatment of stability dependency functions, influence of moisture, and turbulence length scales. Additionally, new features were added to the turbulence scheme: improved representation of turbulence in very stable stratification (absence of critical gradient Richardson number), better representation of anisotropy, unified parameterization of turbulence and clouds via addition of second prognostic turbulence energy, and parameterization of third order moments.


Author(s):  
J. D. Schwarzkopf ◽  
C. T. Crowe ◽  
P. Dutta

A dissipation transport equation for the carrier phase of particle-laden turbulent flows was recently developed. This equation shows a new production of dissipation term due to the presence of particles that is related to the velocity difference between the particle and the surrounding fluid. In the development, it was assumed that each coefficient was the sum of the coefficient for single phase flow and a coefficient quantifying the contribution of the particulate phase. The coefficient for the new production term (due to the presence of particles) was found from homogeneous turbulence generation by particles and the coefficient for the dissipation of dissipation term was analyzed using DNS. A numerical model was developed and applied to particles falling in a channel of downward turbulent air flow. Boundary conditions were also developed to ensure that the production of turbulent kinetic energy due to mean velocity gradients and particle surfaces balanced with the turbulent dissipation near the wall. The turbulent kinetic energy is compared with experimental data. The results show attenuation of turbulent kinetic energy with increased particle loading; however the model does under predict the turbulent kinetic energy near the center of the channel. To understand the effect of this additional production of dissipation term (due to particles), the coefficients associated with the production of dissipation due to mean velocity gradients and particle surfaces are varied to assess the effects of the dispersed phase on the carrier phase turbulent kinetic energy across the channel. The results show that this additional term plays a significant role in predicting the turbulent kinetic energy and a reason for under predicting the turbulent kinetic energy near the center of the channel is discussed. It is concluded that the dissipation coefficients play a critical role in predicting the turbulent kinetic energy in particle-laden turbulent flows.


Sign in / Sign up

Export Citation Format

Share Document