5 results
Measurements of particle dispersion obtained from direct numerical simulations of isotropic turbulence
- Kyle D. Squires, John K. Eaton
-
- Journal:
- Journal of Fluid Mechanics / Volume 226 / May 1991
- Published online by Cambridge University Press:
- 26 April 2006, pp. 1-35
-
- Article
- Export citation
-
Measurements of heavy particle dispersion have been made using direct numerical simulations of isotropic turbulence. The parameters affecting the dispersion of solid particles, namely particle inertia and drift due to body forces were investigated separately. In agreement with the theoretical studies of Reeks, and Pismen & Nir, the effect of particle inertia is to increase the eddy diffusivity over that of the fluid (in the absence of particle drift). The increase in the eddy diffusivity of particles over that of the fluid was between 2 and 16%, in reasonable agreement with the increases reported in Reeks, and Pismen & Nir. The effect of a deterministic particle drift is shown to decrease unequally the dispersion in directions normal and parallel to the particle drift direction. Eddy diffusivities normal and parallel to particle drift are shown to be in good agreement with the predictions of Csanady and the experimental measurements of Wells & Stock.
Partitioning of particle velocities in gas–solid turbulent flows into a continuous field and a spatially uncorrelated random distribution: theoretical formalism and numerical study
- PIERRE FÉVRIER, OLIVIER SIMONIN, KYLE D. SQUIRES
-
- Journal:
- Journal of Fluid Mechanics / Volume 533 / 25 June 2005
- Published online by Cambridge University Press:
- 15 June 2005, pp. 1-46
-
- Article
- Export citation
-
The velocity distribution of dilute suspensions of heavy particles in gas–solid turbulent flows is investigated. A statistical approach – the mesoscopic Eulerian formalism (MEF) – is developed in which an average conditioned on a realization of the turbulent carrier flow is introduced and enables a decomposition of the instantaneous particle velocity into two contributions. The first is a contribution from an underlying continuous turbulent velocity field shared by all the particles – the mesoscopic Eulerian particle velocity field (MEPVF) – that accounts for all particle–particle and fluid–particle two-point correlations. The second contribution corresponds to a distribution – the quasi-Brownian velocity distribution (QBVD) – that represents a random velocity component satisfying the molecular chaos assumption that is not spatially correlated and identified with each particle of the system. The MEF is used to investigate properties of statistically stationary particle-laden isotropic turbulence. The carrier flow is computed using direct numerical simulation (DNS) or large-eddy simulation (LES) with discrete particle tracking employed for the dispersed phase. Particle material densities are much larger than that of the fluid and the force of the fluid on the particle is assumed to reduce to the drag contribution. Computations are performed in the dilute regime for which the influences of inter-particle collisions and fluid-turbulence modulation are neglected. The simulations show that increases in particle inertia increase the contribution of the quasi-Brownian component to the particle velocity. The particle velocity field is correlated at larger length scales than the fluid, with the integral length scales of the MEPVF also increasing with particle inertia. Consistent with the previous work of Abrahamson (1975), the MEF shows that in the limiting case of large inertia, particle motion becomes stochastically equivalent to a Brownian motion with a random spatial distribution of positions and velocities. For the current system of statistically stationary isotropic turbulence, both the DNS and LES show that the fraction of the kinetic energy residing in the mesoscopic field decreases with particle inertia as the square root of the ratio of the total particulate-phase kinetic energy to that of the fluid.
Prediction and investigation of the turbulent flow over a rotating disk
- XIAOHUA WU, KYLE D. SQUIRES
-
- Journal:
- Journal of Fluid Mechanics / Volume 418 / 10 September 2000
- Published online by Cambridge University Press:
- 10 September 2000, pp. 231-264
-
- Article
- Export citation
-
Large-eddy simulation (LES) has been used to predict the statistically three-dimensional turbulent boundary layer (3DTBL) over a rotating disk. LES predictions for six parameter cases were compared to the experimental measurements of Littell & Eaton (1994), obtained at a momentum thickness Reynolds number of 2660. A signal-decomposition scheme was developed by modifying the method of Spalart (1988) to prescribe time-dependent boundary conditions along the radial direction, entrainment towards the disk surface was prescribed by satisfying global mass conservation. Predictions of the mean velocities and r.m.s. fluctuations are in good agreement with data, with the largest discrepancy occurring in the prediction of the wall-normal intensities. The primary and two secondary shear stresses are also in good agreement with the measurements and one-dimensional energy spectra of the velocity fluctuations agree well with established laws, i.e. a −1 slope in the buffer region and −5/3 slope near the edge of the boundary layer.
Conditionally averaged velocities provide new evidence in support of the structural model of Littell & Eaton (1994) concerning the interaction of mean-flow three-dimensionality and shear-stress producing structures. Inside the buffer region under strong ejections, the conditionally averaged crossflow (radial) velocity is larger than the unconditioned mean, and the profile conditioned on strong sweeps is smaller than the mean. This is consistent with the notion that streamwise vortices having the same sign as the mean streamwise vorticity, and beneath the peak crossflow location, are mostly responsible for strong sweep events; streamwise vortices with opposite sign as the mean streamwise vorticity promote strong ejections. Comparison of two-point spatial correlations with previous measurements in two-dimensional turbulent boundary layers (2DTBLs) indicates interesting structural similarities, e.g. the correlation of wall pressure and surface-normal velocity fluctuations is an odd function of streamwise separation, being positive downstream and negative upstream. These similarities offer quantitative indirect support to the hypothesis advanced by Littell & Eaton (1994) and Johnston & Flack (1996) that structural models describing 2DTBLs may be employed as a baseline in (equilibrium) 3DTBL structural studies.
Direct numerical simulation of turbulence modulation by particles in isotropic turbulence
- MARC BOIVIN, OLIVIER SIMONIN, KYLE D. SQUIRES
-
- Journal:
- Journal of Fluid Mechanics / Volume 375 / 25 November 1998
- Published online by Cambridge University Press:
- 25 November 1998, pp. 235-263
-
- Article
- Export citation
-
The modulation of isotropic turbulence by particles has been investigated using direct numerical simulation (DNS). The particular focus of the present work is on the class of dilute flows in which particle volume fractions and inter-particle collisions are negligible. Gravitational settling is also neglected and particle motion is assumed to be governed by drag with particle relaxation times ranging from the Kolmogorov scale to the Eulerian time scale of the turbulence and particle mass loadings up to 1. The velocity field was made statistically stationary by forcing the low wavenumbers of the flow. The calculations were performed using 963 collocation points and the Taylor-scale Reynolds number for the stationary flow was 62. The effect of particles on the turbulence was included in the Navier–Stokes equations using the point-force approximation in which 963 particles were used in the calculations. DNS results show that particles increasingly dissipate fluid kinetic energy with increased loading, with the reduction in kinetic energy being relatively independent of the particle relaxation time. Viscous dissipation in the fluid decreases with increased loading and is larger for particles with smaller relaxation times. Fluid energy spectra show that there is a non-uniform distortion of the turbulence with a relative increase in small-scale energy. The non-uniform distortion significantly affects the transport of the dissipation rate, with the production and destruction of dissipation exhibiting completely different behaviours. The spectrum of the fluid–particle energy exchange rate shows that the fluid drags particles at low wavenumbers while the converse is true at high wavenumbers for small particles. A spectral analysis shows that the increase of the high-wavenumber portion of the fluid energy spectrum can be attributed to transfer of the fluid–particle covariance by the fluid turbulence. This in turn explains the relative increase of small-scale energy caused by small particles observed in the present simulations as well as those of Squires & Eaton (1990) and Elghobashi & Truesdell (1993).
Numerical investigation of the turbulent boundary layer over a bump
- XIAOHUA WU, KYLE D. SQUIRES
-
- Journal:
- Journal of Fluid Mechanics / Volume 362 / 10 May 1998
- Published online by Cambridge University Press:
- 10 May 1998, pp. 229-271
-
- Article
- Export citation
-
Large-eddy simulation (LES) has been used to calculate the flow of a statistically two-dimensional turbulent boundary layer over a bump. Subgrid-scale stresses in the filtered Navier–Stokes equations were closed using the dynamic eddy viscosity model. LES predictions for a range of grid resolutions were compared to the experimental measurements of Webster et al. (1996). Predictions of the mean flow and turbulence intensities are in good agreement with measurements. The resolved turbulent shear stress is in reasonable agreement with data, though the peak is over-predicted near the trailing edge of the bump. Analysis of the flow confirms the existence of internal layers over the bump surface upstream of the summit and along the downstream trailing at plate, and demonstrates that the quasi-step increases in skin friction due to perturbations in pressure gradient and surface curvature selectively enhance near-wall shear production of turbulent stresses and are responsible for the formation of the internal layers. Though the flow experiences a strong adverse pressure gradient along the rear surface, the boundary layer is unique in that intermittent detachment occurring near the wall is not followed by mean-flow separation. Certain turbulence characteristics in this region are similar to those previously reported in instantaneously separating boundary layers. The present investigation also explains the driving mechanism for the surprisingly rapid return to equilibrium over the trailing flat plate found in the measurements of Webster et al. (1996), i.e. the simultaneous uninterrupted development of an inner energy-equilibrium region and the monotonic decay of elevated turbulence shear production away from the wall. LES results were also used to examine response of the dynamic eddy viscosity model. While subgrid-scale dissipation decreases/increases as the turbulence is attenuated/enhanced, the ratio of the (averaged) forward to reverse energy transfers predicted by the model is roughly constant over a significant part of the layer. Model predictions of backscatter, calculated as the percentage of points where the model coefficient is negative, show a rapid recovery downstream similar to the mean-flow and turbulence quantities.