Normalized dissipation rate in a moderate Taylor Reynolds number flow

2017 ◽  
Vol 818 ◽  
pp. 184-204 ◽  
Author(s):  
Alejandro J. Puga ◽  
John C. LaRue

Time-resolved velocity measurements are obtained using a hot-wire in a nearly homogeneous and isotropic flow downstream of an active grid for a range of Taylor Reynolds numbers from$191$to$659$. It is found that the dimensionless dissipation rate,$C_{\unicode[STIX]{x1D716}}$, is nearly a constant for sufficiently high values of Taylor Reynolds number,$R_{\unicode[STIX]{x1D706},u_{q}}$, and is approximately equal to$0.87$. This value is approximately$5\,\%$less than the value reported by Boset al.(Phys. Fluids, vol. 19 (4), 2007, 045101), which is obtained using DNS/LES (direct numerical simulation combined with large eddy simulation) for decaying homogeneous and isotropic turbulence, and is in excellent agreement with the active grid experiment of Thormann & Meneveau (Phys. Fluids, vol. 26 (2), 2014, 025112.). The results presented herein show that deviation from isotropy may cause inconsistencies in the computation of$C_{\unicode[STIX]{x1D716}}$. As a result, it is suggested that the velocity scale be the square root of the turbulence kinetic energy. The integral length scale measurements obtained from the longitudinal velocity correlation are in close agreement with the integral length scale measured from the peak of the energy spectrum,$\unicode[STIX]{x1D705}E_{11}(\unicode[STIX]{x1D705})$, where$\unicode[STIX]{x1D705}$is the wavenumber and$E_{11}(\unicode[STIX]{x1D705})$is the one-dimensional power spectrum of the downstream velocity.

2015 ◽  
Vol 772 ◽  
pp. 361-385 ◽  
Author(s):  
Anikesh Pal ◽  
Sutanu Sarkar

Direct numerical simulations are performed to study the evolution of a towed stratified wake subject to external turbulence in the background. A field of isotropic turbulence is combined with an initial turbulent wake field and the combined wake is simulated in a temporally evolving framework similar to that of Rind & Castro (J. Fluid Mech., vol. 710, 2012a, p. 482). Simulations are performed for external turbulence whose initial level varies between zero and a moderate intensity of up to 7 % relative to the free stream and whose initial integral length scale is of the same order as that of the wake turbulence. A series of simulations are carried out at a Reynolds number of 10 000 and Froude number of 3. Background turbulence, especially at a level of 3 % or above, is found to have substantial quantitative effects in the stratified simulations. Turbulence inside the wake increases due to the entrainment of external turbulence, and the energy transfer through turbulent production from mean to fluctuating velocity also increases, leading to reduced mean velocity. The profiles of normalized mean and turbulence quantities in the stratified wake exhibit little change in the vertical direction but the horizontal spread increases in comparison to the case with undisturbed background. The spatial organization of the internal wave field is disrupted even at the 1 % level of external turbulence. However, key characteristics of stratified wakes such as the formation of coherent pancake vortices and the long lifetime of the mean wake are robust to the presence of fluctuations in the background. A corresponding series of simulations for the unstratified situation is carried out at the same Reynolds number of 10 000 and with similar levels of external turbulence. The change of mean and turbulence statistics is found to be weaker in the unstratified cases compared with the corresponding stratified cases and also weaker relative to that found by Rind & Castro (J. Fluid Mech., vol. 710, 2012a, p. 482) at a similar level of external turbulence relative to the free stream and similar integral length scale. Theoretical arguments and additional simulations are provided to show that the level of external turbulence relative to wake turbulence (dissimilar between the present investigation and Rind & Castro (J. Fluid Mech., vol. 710, 2012a, p. 482)) is a key governing parameter in both stratified and unstratified backgrounds.


1998 ◽  
Vol 360 ◽  
pp. 249-271 ◽  
Author(s):  
H. DÜTSCH ◽  
F. DURST ◽  
S. BECKER ◽  
H. LIENHART

Time-averaged LDA measurements and time-resolved numerical flow predictions were performed to investigate the laminar flow induced by the harmonic in-line oscillation of a circular cylinder in water at rest. The key parameters, Reynolds number Re and Keulegan–Carpenter number KC, were varied to study three parameter combinations in detail. Good agreement was observed for Re=100 and KC=5 between measurements and predictions comparing phase-averaged velocity vectors. For Re=200 and KC=10 weakly stable and non-periodic flow patterns occurred, which made repeatable time-averaged measurements impossible. Nevertheless, the experimentally visualized vortex dynamics was reproduced by the two-dimensional computations. For the third combination, Re=210 and KC=6, which refers to a totally different flow regime, the computations again resulted in the correct fluid behaviour. Applying the widely used model of Morison et al. (1950) to the computed in-line force history, the drag and the added-mass coefficients were calculated and compared for different grid levels and time steps. Using these to reproduce the force functions revealed deviations from those originally computed as already noted in previous studies. They were found to be much higher than the deviations for the coarsest computational grid or the largest time step. The comparison of several in-line force coefficients with results obtained experimentally by Kühtz (1996) for β=35 confirmed that force predictions could also be reliably obtained by the computations.


2016 ◽  
Vol 806 ◽  
pp. 356-412 ◽  
Author(s):  
Michael S. Dodd ◽  
Antonino Ferrante

Droplets in turbulent flows behave differently from solid particles, e.g. droplets deform, break up, coalesce and have internal fluid circulation. Our objective is to gain a fundamental understanding of the physical mechanisms of droplet–turbulence interaction. We performed direct numerical simulations (DNS) of 3130 finite-size, non-evaporating droplets of diameter approximately equal to the Taylor length scale and with 5 % droplet volume fraction in decaying isotropic turbulence at initial Taylor-scale Reynolds number $\mathit{Re}_{\unicode[STIX]{x1D706}}=83$. In the droplet-laden cases, we varied one of the following three parameters: the droplet Weber number based on the r.m.s. velocity of turbulence ($0.1\leqslant \mathit{We}_{rms}\leqslant 5$), the droplet- to carrier-fluid density ratio ($1\leqslant \unicode[STIX]{x1D70C}_{d}/\unicode[STIX]{x1D70C}_{c}\leqslant 100$) or the droplet- to carrier-fluid viscosity ratio ($1\leqslant \unicode[STIX]{x1D707}_{d}/\unicode[STIX]{x1D707}_{c}\leqslant 100$). In this work, we derive the turbulence kinetic energy (TKE) equations for the two-fluid, carrier-fluid and droplet-fluid flow. These equations allow us to explain the pathways for TKE exchange between the carrier turbulent flow and the flow inside the droplet. We also explain the role of the interfacial surface energy in the two-fluid TKE equation through the power of the surface tension. Furthermore, we derive the relationship between the power of surface tension and the rate of change of total droplet surface area. This link allows us to explain how droplet deformation, breakup and coalescence play roles in the temporal evolution of TKE. Our DNS results show that increasing $\mathit{We}_{rms}$, $\unicode[STIX]{x1D70C}_{d}/\unicode[STIX]{x1D70C}_{c}$ and $\unicode[STIX]{x1D707}_{d}/\unicode[STIX]{x1D707}_{c}$ increases the decay rate of the two-fluid TKE. The droplets enhance the dissipation rate of TKE by enhancing the local velocity gradients near the droplet interface. The power of the surface tension is a source or sink of the two-fluid TKE depending on the sign of the rate of change of the total droplet surface area. Thus, we show that, through the power of the surface tension, droplet coalescence is a source of TKE and breakup is a sink of TKE. For short times, the power of the surface tension is less than $\pm 5\,\%$ of the dissipation rate. For later times, the power of the surface tension is always a source of TKE, and its magnitude can be up to 50 % of the dissipation rate.


1995 ◽  
Vol 117 (3) ◽  
pp. 597-603 ◽  
Author(s):  
G. J. Van Fossen ◽  
R. J. Simoneau ◽  
C. Y. Ching

This experiment investigated the effects of free-stream turbulence intensity, length scale, Reynolds number, and leading-edge velocity gradient on stagnation-region heat transfer. Heat transfer was measured in the stagnation region of four models with elliptical leading edges downstream of five turbulence-generating grids. Stagnation-region heat transfer augmentation increased with decreasing length scale but ann optimum scale was not found. A correlation was developed that fit heat transfer data for isotropic turbulence to within ±4 percent but did not predict data for anisotropic turbulence. Stagnation heat transfer augmentation caused by turbulence was unaffected by the velocity gradient. The data of other researchers compared well with the correlation. A method of predicting heat transfer downstream of the stagnation point was developed.


1997 ◽  
Vol 345 ◽  
pp. 307-345 ◽  
Author(s):  
SHIGEO KIDA ◽  
SUSUMU GOTO

A set of integro-differential equations in the Lagrangian renormalized approximation (Kaneda 1981) is rederived by applying a perturbation method developed by Kraichnan (1959), which is based upon an extraction of direct interactions among Fourier modes of a velocity field and was applied to the Eulerian velocity correlation and response functions, to the Lagrangian ones for homogeneous isotropic turbulence. The resultant set of integro-differential equations for these functions has no adjustable free parameters. The shape of the energy spectrum function is determined numerically in the universal range for stationary turbulence, and in the whole wavenumber range in a similarly evolving form for the freely decaying case. The energy spectrum in the universal range takes the same shape in both cases, which also agrees excellently with many measurements of various kinds of real turbulence as well as numerical results obtained by Gotoh et al. (1988) for a decaying case as an initial value problem. The skewness factor of the longitudinal velocity derivative is calculated to be −0.66 for stationary turbulence. The wavenumber dependence of the eddy viscosity is also determined.


Author(s):  
Ian Proudman ◽  
G. K. Batchelor

AbstractIn this paper, the theoretical double and triple velocity correlation functions, f(r), g(r) and h(r), which correspond to Heisenberg's spectrum of isotropic turbulence, are obtained numerically for two Reynolds numbers. One set of these correlations is for the limiting case of infinite Reynolds number. In addition, a method is developed for deriving the approximate form of the double correlations for any Reynolds number, which is not too small, from the corresponding correlations for infinite Reynolds number. These theoretical correlations are then compared with the results of experiment.


This paper discusses a type of turbulence in a uniform stream which is next to isotropic turbulence in order of simplicity. Instead of spherical symmetry, or isotropy, axially symmetical turbulence possesses symmetry about an axis which in practice is usually the direction of mean flow. The analysis is developed with the aid of invariant theory, as suggested by a previous paper by Robertson. The form of the fundamental velocity correlation is obtained, and scales of axisymmetric turbulence are defined. The results of greatest practical interest concern the time rates of change of the mean squares of the lateral and longitudinal velocity components. The rates of change involve two terms, the first representing viscous dissipation, and the second representing a transfer of energy from one component to the other due to the finite correlation between the velocity and pressure at neighbouring points. The effect of the velocity-pressure correlation is to bring the two velocity components towards equality, while the effect of the viscous dissipation will only be towards equality if an inequality between the curvatures at the origin of two particular velocity correlation coefficient curves, both of which are measurable, is obeyed. The rates of change of the mean squares of the vorticity components are also obtained.


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.


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.


Sign in / Sign up

Export Citation Format

Share Document