Published online by Cambridge University Press: 19 March 2010
This study investigates the two-way coupling effects of finite-size solid spherical particles on decaying isotropic turbulence using direct numerical simulation with an immersed boundary method. We fully resolve all the relevant scales of turbulence around freely moving particles of the Taylor length-scale size, 1.2≤d/λ≤2.6. The particle diameter and Stokes number in terms of Kolmogorov length- and time scales are 16≤d/η≤35 and 38≤τp/τk≤178, respectively, at the time the particles are released in the flow. The particles mass fraction range is 0.026≤φm≤1.0, corresponding to a volume fraction of 0.01≤φv≤0.1 and density ratio of 2.56≤ρp/ρf≤10. The maximum number of dispersed particles is 6400 for φv=0.1. The typical particle Reynolds number is of O(10). The effects of the particles on the temporal development of turbulence kinetic energy E(t), its dissipation rate (t), its two-way coupling rate of change Ψp(t) and frequency spectra E(ω) are discussed.
In contrast to particles with d < η, the effect of the particles in this study, with d > η, is that E(t) is always smaller than that of the single-phase flow. In addition, Ψp(t) is always positive for particles with d > η, whereas it can be positive or negative for particles with d < η.
Movie 1. Instantaneous contours of sijsij in xz-plane (1.5 < t < 3) for single-phase flow (case A, left) and particle-laden flow (case F, right). In case F, the spherical uniform-size particles appear to be of different sizes because the xz-plane intersects with them while they are moving in the third direction (perpendicular to the plane). Red regions of high sijsij are formed downstream of the particles, thus increasing the dissipation rate of TKE. The speed of these animations is 1000 times slower than the actual speed for the case of liquid water carrier fluid.
Movie 2. TKE in xz-plane (1.5 < t < 3) for single-phase flow (case A, left) and particle-laden flow (case F, right). The high TKE regions (red) in case A are reduced by the particles in case F.
Movie 3. Instantaneous contours of ωy in xz-plane (1.5 < t < 3) for single-phase flow (case A, left) and particle-laden flow (case F, right). Counter-rotating (blue and red) vortices are created adjacent to the particle surface, thus increasing the ensemble-average enstrophy <ω2> (see Table 5).