Turbulent kinetic energy production and flow structures in flows past smooth and rough walls

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.

1996 ◽  
Vol 118 (4) ◽  
pp. 857-863 ◽  
Author(s):  
O. Kwon ◽  
F. E. Ames

This paper describes a velocity and length scale approach to low-Reynolds-number k–ε modeling, which formulates the eddy viscosity on the normal component of turbulence and a length scale. The normal component of turbulence is modeled based on the dissipation and distance from the wall and is bounded by the isotropic condition. The model accounts for the anisotropy of the dissipation and the reduced length of mixing in the near wall region. The kinetic energy and dissipation rate were computed from the k and ε transport equations of Durbin (1993). The model was tested for a wide range of turbulent flows and proved to be superior to other k–ε based models.


Author(s):  
T. Kawamura ◽  
T. Nakatani

Direct numerical simulations of homogeneous shear turbulent flows containing deformable bubbles were carried out for clarifying the mechanism of drag reduction by microbubbles. The results show that presence of bubbles can suppress or enhance the development of turbulence depending on condition. The dissipation rate of turbulent kinetic energy is always increased by bubbles, while the production rate can be either increased or decreased depending on the turbulent and shear Reynolds numbers. As a result, the growth rate of turbulent kinetic energy can be either increased or decreased by bubbles depending on conditions. It was shown that the production rate tends to decrease at smaller shear Reynolds number, at larger turbulent Reynolds number, and at larger Weber number. Based on the results, a hypothetical model to explain the dependency on the Reynolds numbers has been proposed.


1988 ◽  
Vol 110 (2) ◽  
pp. 187-194 ◽  
Author(s):  
S. A. Shirazi ◽  
C. R. Truman

An anisotropic form of a low-Reynolds-number two-equation turbulence model has been implemented in a numerical solution for incompressible turbulent flow between corotating parallel disks. Transport equations for turbulent kinetic energy and dissipation rate were solved simultaneously with the governing equations for the mean-flow variables. Comparisons with earlier mixing-length predictions and with measurements are presented. Good agreement between the present predictions and the measurements of velocity components and turbulent kinetic energy was obtained. The low-Reynolds-number two-equation model was found to model adequately the near-wall region as well as the effects of rotation and streamline divergence, which required ad hoc assumptions in the mixing-length model.


2021 ◽  
Vol 2119 (1) ◽  
pp. 012011
Author(s):  
K G Dobroselsky

Abstract Vortex flow structures in a turbulent wake behind a circular Teflon cylinder immersed in an incoming flow with a change in pressure for the Reynolds number Re = 2.2×105 have been experimentally studied using a two-dimensional image (2D-PIV) of particles in a closed-circuit water tunnel. The obtained results are presented in the form of time-averaged velocity fields, Reynolds stresses, and distributions of turbulent kinetic energy. The flow data showed that the size of the wake flow region, Reynolds stresses and turbulent kinetic energy change depending on the pressure in the flow. As a result of a 20% reduction in pressure, the size of the vortex zone in the wake increases by about 20%.


2020 ◽  
Vol 24 (5 Part A) ◽  
pp. 2729-2741
Author(s):  
Zhenchuan Wang ◽  
Guoli Qi ◽  
Meijun Li

The turbulence model fails in supercritical fluid-flow and heat transfer simulation, owing to the drastic change of thermal properties. The inappropriate buoyancy effect model and the improper turbulent Prandtl number model are several of these factors lead to the original low-Reynolds number turbulence model unable to predict the wall temperature for vertically heated tubes under the deteriorate heat transfer conditions. This paper proposed a simplified improved method to modify the turbulence model, using the generalized gradient diffusion hypothesis approximation model for the production term of the turbulent kinetic energy due to the buoyancy effect, using a turbulence Prandtl number model for the turbulent thermal diffusivity instead of the constant number. A better agreement was accomplished by the improved turbulence model compared with the experimental data. The main reason for the over-predicted wall temperature by the original turbulence model is the misuse of the buoyancy effect model. In the improved model, the production term of the turbulent kinetic energy is much higher than the results calculated by the original turbulence model, especially in the boundary-layer. A more accurate model for the production term of the turbulent kinetic energy is the main direction of further modification for the low Reynolds number turbulence model.


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.


2019 ◽  
Vol 1 (2) ◽  
pp. 265-282
Author(s):  
Marta Puzdrowska ◽  
Tomasz Heese

The paper presents an analysis the spatial distribution of turbulent kinetic energy (TKE) for bolt fishways, including the impact of additional spillway slots and fixed channel development. The research was done for two models, each containing a different arrangement of slots. The presented results of research for bolt fishways were obtained as an effect of laboratory tests. The measurements were done for three components of instant flow velocity magnitude (speed). Analysis of the results was done for a 3D flow structure using Matlab software. In the case of bolt fishways, significant differences were noted for the method of velocity and TKE distribution, in reference to research comprising channels with biological development. It was stated that a reason for this is the flexible development of the channel. The occurrence of extreme TKE values in the chamber (pool) is strictly associated with the characteristics of interaction zones between various flow structures. It was also stated that the lower the parapet of the slot’s spillway shelf is in the fishway’s partition, the higher TKE could be expected just downstream of the section. These establishments may be important for the designing process in the case of fish passes of various types of construction.


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.


Sign in / Sign up

Export Citation Format

Share Document