Transport equations for the mean energy and temperature dissipation rates in grid turbulence

2000 ◽  
Vol 28 (2) ◽  
pp. 143-151 ◽  
Author(s):  
T. Zhou ◽  
R. A. Antonia ◽  
L. Danaila ◽  
F. Anselmet
2021 ◽  
Vol 930 ◽  
Author(s):  
S.L. Tang ◽  
R.A. Antonia ◽  
L. Djenidi

Transport equations for the normalized moments of the longitudinal velocity derivative ${F_{n + 1}}$ (here, $n$ is $1, 2, 3\ldots$ ) are derived from the Navier–Stokes (N–S) equations for shearless grid turbulence. The effect of the (large-scale) streamwise advection of ${F_{n + 1}}$ by the mean velocity on the normalized moments of the velocity derivatives can be expressed as $C_1 {F_{n + 1}}/Re_\lambda$ , where $C_1$ is a constant and $Re_\lambda$ is the Taylor microscale Reynolds number. Transport equations for the normalized odd moments of the transverse velocity derivatives ${F_{y,n + 1}}$ (here, $n$ is 2, 4, 6), which should be zero if local isotropy is satisfied, are also derived and discussed in sheared and shearless grid turbulence. The effect of the (large-scale) streamwise advection term on the normalized moments of the velocity derivatives can also be expressed in the form $C_2 {F_{y,n + 1}}/Re_\lambda$ , where $C_2$ is a constant. Finally, the contribution of the mean shear in the transport equation for ${F_{n + 1}}$ can be modelled as $15 B/Re_\lambda$ , where $B$ ( $=S^*{S_{s,n + 1}}$ ) is the product of the non-dimensional shear parameter $S^*$ and the normalized mixed longitudinal-transverse velocity derivatives ${{S_{s,n + 1}}}$ ; if local isotropy is satisfied, $S_{s,n + 1}$ should be zero. These results indicate that if ${F_{n + 1}}$ , ${F_{y,n + 1}}$ and $B$ do not increase as rapidly as $Re_\lambda$ , then the effect of the large-scale structures on small-scale turbulence will disappear when $Re_\lambda$ becomes sufficiently large.


2013 ◽  
Vol 730 ◽  
pp. 593-606 ◽  
Author(s):  
L. Djenidi ◽  
S. F. Tardu ◽  
R. A. Antonia

AbstractA long-time direct numerical simulation (DNS) based on the lattice Boltzmann method is carried out for grid turbulence with the view to compare spatially averaged statistical properties in planes perpendicular to the mean flow with their temporal counterparts. The results show that the two averages become equal a short distance downstream of the grid. This equality indicates that the flow has become homogeneous in a plane perpendicular to the mean flow. This is an important result, since it confirms that hot-wire measurements are appropriate for testing theoretical results based on spatially averaged statistics. It is equally important in the context of DNS of grid turbulence, since it justifies the use of spatial averaging along a lateral direction and over several realizations for determining various statistical properties. Finally, the very good agreement between temporal and spatial averages validates the comparison between temporal (experiments) and spatial (DNS) statistical properties. The results are also interesting because, since the flow is stationary in time and spatially homogeneous along lateral directions, the equality between the two types of averaging provides strong support for the ergodic hypothesis in grid turbulence in planes perpendicular to the mean flow.


Author(s):  
Hossam Elasrag ◽  
Shaoping Li

Simulations for the Cambridge swirl bluff-body spray burner are performed near blow-out conditions. A hybrid stress blended eddy simulation (SBES) model is utilized for sub-grid turbulence closure. SBES blends the RANS-SST model at the boundary layer with large eddy simulation dynamic Smagorinsky model outside the boundary layer. The injected N-heptane spray droplets are tracked using a typical Eulerian-Lagrangian approach. Heat transfer coupling between the bluff-body walls and the near-walls fluid is accounted for by coupling the solid and fluid energy equations at the bluff-body surface. Mixing and chemistry are modeled using the Flamelet Generated Manifold (FGM) model. The study investigates how successful the FGM model is in predicting finite rate effects like local extinction and flame lift-off height. To this end, two near blow-out spray flames, the H1S1 (75% to blow-out) and H1S2 (88% to blow-out) are simulated. Good results are shown matching the spray Sauter mean diameter (SMD) and axial velocity mean and rms experimental data. The results also show that the FGM model captured reasonably well the flame structure and lift-off height as well as the spray pattern. Overall the spray droplets mean D32 and mean axial velocity were under-predicted, while the rms distribution matched reasonably well for the H1S1 flame. The mean flame brush lift-off height is estimated based on the statistically stationary mean flame brush and is estimated to be around 6 mm from the bluff-body base. Instantaneous local flame extinction is observed. The H1S2 flame, however, showed similar but slightly better match with the measurements for the mean spray data compared to the H1S1 flame, with slight under-prediction for D32 at Z = 10 mm and Z = 20 mm. Future work will investigate the sensitivity of the simulation to the spray boundary conditions and grid resolution.


1982 ◽  
Vol 104 (4) ◽  
pp. 523-527 ◽  
Author(s):  
J. Courchesne ◽  
A. Laneville

This paper describes an experimental evaluation of the effects of the intensity and scale of turbulence on the drag coefficient of two-dimensional rectangular cylinders exposed to grid turbulence. It is observed that the mean drag coefficient is principally influenced, for a given cylinder, by the intensity of turbulence and that the scale of turbulence plays a secondary role.


1990 ◽  
Vol 216 ◽  
pp. 35-70 ◽  
Author(s):  
S. Veeravalli ◽  
Z. Warhaft

We experimentally investigate dispersion from a heated line source placed in the central region of a turbulence mixing layer. Recently described by Veeravalli & Warhaft (1989) the mixing layer has no mean shear and consists of gradients in the velocity variance and scale; it is formed from a composite grid of constant solidity from which two distinct velocity scales are formed, one on either side of the stream. Mixing is effected by intermittent turbulent penetration and diffusion. The dispersion measurements were carried out in the convective regime where both plume flapping and fine-scale turbulent mixing play a role, the latter becoming more dominant as the plume evolves. The mean and variance temperature profiles are strongly skewed (with larger tails on the low turbulence side of the flow) in the earlier stages of the plume development. Here, in the convective range, the median and peak of the mean plume are deflected toward the large-scale region. As the flow evolves the profiles become more symmetrical but as the plume enters the turbulent diffusive stage there is evidence that the profiles again became asymmetric but now with longer tails in the high turbulence side of the flow (owing to the higher diffusivity). The temperature variance and heat flux budgets are highly asymmetric but tend to exhibit many of the characteristics of the budget of a line source in decaying homogeneous grid turbulence which is also presented here. However, a distinct region of negative production (counter-gradient heat flux) is found in the temperature variance budget and this is shown to be a consequence of the asymmetry of the transverse velocity probability density function in the mixing layer. Temperature spectra, both of the time series and of the intermittency function, across the plume are described. They are shown to peak at high wavenumbers in the centre and edge of the plume and at lower wavenumbers in the intermediate region. Their form is shown to change as the plume develops fine-scale structure and flapping becomes less important.


2015 ◽  
Vol 777 ◽  
pp. 151-177 ◽  
Author(s):  
S. L. Tang ◽  
R. A. Antonia ◽  
L. Djenidi ◽  
H. Abe ◽  
T. Zhou ◽  
...  

The transport equation for the mean turbulent energy dissipation rate $\overline{{\it\epsilon}}$ along the centreline of a fully developed channel flow is derived by applying the limit at small separations to the two-point budget equation. Since the ratio of the isotropic energy dissipation rate to the mean turbulent energy dissipation rate $\overline{{\it\epsilon}}_{iso}/\overline{{\it\epsilon}}$ is sufficiently close to 1 on the centreline, our main focus is on the isotropic form of the transport equation. It is found that the imbalance between the production of $\overline{{\it\epsilon}}$ due to vortex stretching and the destruction of $\overline{{\it\epsilon}}$ caused by the action of viscosity is governed by the diffusion of $\overline{{\it\epsilon}}$ by the wall-normal velocity fluctuation. This imbalance is intrinsically different from the advection-driven imbalance in decaying-type flows, such as grid turbulence, jets and wakes. In effect, the different types of imbalance represent different constraints on the relation between the skewness of the longitudinal velocity derivative $S_{1,1}$ and the destruction coefficient $G$ of enstrophy in different flows, thus resulting in non-universal approaches of $S_{1,1}$ towards a constant value as the Taylor microscale Reynolds number, $R_{{\it\lambda}}$, increases. For example, the approach is slower for the measured values of $S_{1,1}$ along either the channel or pipe centreline than along the axis in the self-preserving region of a round jet. The data for $S_{1,1}$ collected in different flows strongly suggest that, in each flow, the magnitude of $S_{1,1}$ is bounded, the value being slightly larger than 0.5.


2000 ◽  
Vol 416 ◽  
pp. 117-149 ◽  
Author(s):  
D. LIVESCU ◽  
F. A. JABERI ◽  
C. K. MADNIA

The structure and development of the scalar wake produced by a single line source are studied in decaying isotropic turbulence. The incompressible Navier–Stokes and the passive-scalar transport equations are solved via direct numerical simulations (DNS). The velocity and the scalar fields are generated by simulating Warhaft's (1984) experiment. The results for mean and r.m.s. scalar statistics are in good agreement with those obtained from the experiment. The structure of the scalar wake is examined first. At initial times, most of the contribution to the scalar variance is due to the flapping of the wake around the centreline. Near the end of the turbulent convective regime, the wake develops internal structure and the contribution of the flapping component to the scalar variance becomes negligible. The influence of the source size on the development of the scalar wake has been examined for source sizes ranging from the Kolmogorov microscale to the integral scale. After an initial development time, the half-widths of mean and scalar r.m.s. wakes grow at rates independent of the source size. The mixing in the scalar wake is studied by analysing the evolution of the terms in the transport equations for mean, scalar flux, variance, and scalar dissipation. The DNS results are used to test two types of closures for the mean and the scalar variance equations. For the time range simulated, the gradient diffusion model for the scalar flux and the commonly used scalar dissipation model are not supported by the DNS data. On the other hand, the model based on the unconditional probability density function (PDF) method predicts the scalar flux reasonably well near the end of the turbulent convective regime for the highest Reynolds number examined. The scalar source size does not significantly influence the models' predictions, although it appears that the time-scale ratio of mechanical dissipation to scalar dissipation approaches an asymptotic value earlier for larger source sizes.


2013 ◽  
Vol 43 (4) ◽  
pp. 824-835 ◽  
Author(s):  
Zhao Jing ◽  
Lixin Wu

Abstract Profiles of potential density obtained from CTD measurements during the Hawaii Ocean Time series (HOT) program in the vicinity of the island of Oahu, Hawaii, are used to evaluate low-frequency variability of turbulent kinetic dissipation rates based on a finescale parameterization method. A distinct seasonal cycle, as well as an increasing trend of dissipation rates, is found in the upper 300–600 m. The trend is mainly due to the much weaker diapycnal mixing in the first four years of the record, that is, 1988–92. In the upper 300–600 m, enhanced diapycnal mixing is found under anticyclonic eddies with the mean dissipation rate about 53% larger than that under eddy-free conditions. The modulation of dissipation rates by anticyclonic eddies becomes more evident with increasing eddy strength. The role of cyclonic eddies in modulating diapycnal mixing is almost negligible compared with that of anticyclonic eddies. The mean dissipation rate under cyclonic eddies is comparable to that under eddy-free conditions with a difference of less than 10%. Seasonality of the dissipation rates is partly modulated by the seasonal variation of anticyclonic eddies.


1998 ◽  
Vol 374 ◽  
pp. 29-57 ◽  
Author(s):  
R. A. ANTONIA ◽  
T. ZHOU ◽  
Y. ZHU

All components of the fluctuating vorticity vector have been measured in decaying grid turbulence using a vorticity probe of relatively simple geometry (four X-probes, i.e. a total of eight hot wires). The data indicate that local isotropy is more closely satisfied than global isotropy, the r.m.s. vorticities being more nearly equal than the r.m.s. velocities. Two checks indicate that the performance of the probe is satisfactory. Firstly, the fully measured mean energy dissipation rate 〈ε〉 is in good agreement with the value inferred from the rate of decay of the mean turbulent energy 〈q2〉 in the quasi-homogeneous region; the isotropic mean energy dissipation rate 〈εiso〉 agrees closely with this value even though individual elements of 〈ε〉 indicate departures from isotropy. Secondly, the measured decay rate of the mean-square vorticity 〈ω2〉 is consistent with that of 〈q2〉 and in reasonable agreement with the isotropic form of the transport equation for 〈ω2〉. Although 〈ε〉≃〈εiso〉, there are discernible differences between the statistics of ε and εiso; in particular, εiso is poorly correlated with either ε or ω2. The behaviour of velocity increments has been examined over a narrow range of separations for which the third-order longitudinal velocity structure function is approximately linear. In this range, transverse velocity increments show larger departures than longitudinal increments from predictions of Kolmogorov (1941). The data indicate that this discrepancy is only partly associated with differences between statistics of locally averaged ε and ω2, the latter remaining more intermittent than the former across this range. It is more likely caused by a departure from isotropy due to the small value of Rλ, the Taylor microscale Reynolds number, in this experiment.


Sign in / Sign up

Export Citation Format

Share Document