We present direct numerical simulation of heavy inertial particles (dust) immersed in two-dimensional turbulent flow (gas). The dust particles are modeled as monodispersed heavy particles capable of modifying the flow through two-way coupling. By varying the Stokes number (St) and the mass-loading parameter (φm), we study the clustering phenomenon and the gas phase kinetic energy spectra. We find that the dust-dust correlation dimension (d2) also depends on φm. In particular, clustering decreases as mass loading (φm) is increased. In the kinetic energy spectra of gas we show (i) the emergence of a different scaling regime and that (ii) the scaling exponent in this regime is not universal but a function of both St and φm. Using a scale-by-scale enstrophy budget analysis we show that in this emerged scaling regime, which we call the dust-dissipative range, viscous dissipation due to the gas balances the back-reaction from the dust.

We consider the kinematic fluctuation dynamo problem in a flow that is random, white-in-time, with both solenoidal and potential components. This model is a generalization of the well-studied Kazantsev model. If both the solenoidal and potential parts have the same scaling exponent, then, as the compressibility of the flow increases, the growth rate decreases but remains positive. If the scaling exponents for the solenoidal and potential parts differ, in particular if they correspond to typical Kolmogorov and Burgers values, we again find that an increase in compressibility slows down the growth rate but does not turn it off. The slow down is, however, weaker and the critical magnetic Reynolds number is lower than when both the solenoidal and potential components display the Kolmogorov scaling. Intriguingly, we find that there exist cases, when the potential part is smoother than the solenoidal part, for which an increase in compressibility increases the growth rate. We also find that the critical value of the scaling exponent above which a dynamo is seen is unity irrespective of the compressibility. Finally, we realize that the dimension d = 3 is special, as for all other values of d the critical exponent is higher and depends on the compressibility.

We present an extensive numerical study of the time irreversibility of the dynamics of heavy inertial particles in three-dimensional, statistically homogeneous, and isotropic turbulent flows. We show that the probability density function (PDF) of the increment, W(tau), of a particle's energy over a time scale tau is non-Gaussian, and skewed toward negative values. This implies that, on average, particles gain energy over a period of time that is longer than the duration over which they lose energy. We call this slow gain and fast loss. We find that the third moment of W(tau) scales as tau(3) for small values of tau. We show that the PDF of power-input p is negatively skewed too; we use this skewness Ir as a measure of the time irreversibility and we demonstrate that it increases sharply with the Stokes number St for small St; this increase slows down at St similar or equal to 1. Furthermore, we obtain the PDFs of t(+) and t(-), the times over which p has, respectively, positive or negative signs, i.e., the particle gains or loses energy. We obtain from these PDFs a direct and natural quantification of the slow gain and fast loss of the energy of the particles, because these PDFs possess exponential tails from which we infer the characteristic loss and gain times t(loss) and t(gain), respectively, and we obtain t(loss) < t(gain) for all the cases we have considered. Finally, we show that the fast loss of energy occurs with greater probability in the strain-dominated region than in the vortical one; in contrast, the slow gain in the energy of the particles is equally likely in vortical or strain-dominated regions of the flow.

We perform direct numerical simulations of a bidisperse suspension of heavy spherical particles in forced, homogeneous, and isotropic three-dimensional turbulence. We compute the joint distribution of relative particle distances and longitudinal relative velocities between particles of different inertia. For a pair of particles with small difference in their inertias we compare our results with recent theoretical predictions [Meibohm et al., Phys. Rev. E 96, 061102 (2017)] for the shape of this distribution. We also compute the moments of relative velocities as a function of particle separation and compare with the theoretical predictions. We observe good agreement. For a pair of particles that are very different from each other-one is heavy and the other one has negligible inertia-we give a theory to calculate their root-mean-square relative velocity. This theory also agrees well with the results of our simulations.

We study suspensions of deformable (viscoelastic) spheres in a Newtonian solvent in planeCouette geometry, by means of direct numerical simulations. We find that in the limit of vanishing inertia, the effective viscosity mu of the suspension increases as the volume fraction occupied by the spheres Phi increases and decreases as the elastic modulus of the spheres G decreases; the function mu(Phi,G) collapses to a universal function mu(Phi(e)) with a reduced effective volume fraction Phi(e)(Phi,G). Remarkably, the function mu(Phi(e)) is the well- known Eilers fit that describes the rheology of suspension of rigid spheres at all Phi. Our results suggest different ways to interpret the macrorheology of blood.

We use direct numerical simulations to calculate the joint probability density function of the relative distance R and relative radial velocity component V-R for a pair of heavy inertial particles suspended in homogeneous and isotropic turbulent flows. At small scales the distribution is scale invariant, with a scaling exponent that is related to the particle-particle correlation dimension in phase space, D-2. It was argued [K. Gustavsson and B. Mehlig, Phys. Rev. E 84, 045304 (2011); J. Turbul. 15, 34 (2014)] that the scale invariant part of the distribution has two asymptotic regimes: (1) vertical bar V-R vertical bar << R, where the distribution depends solely on R, and (2) vertical bar V-R vertical bar >> R, where the distribution is a function of vertical bar V-R vertical bar alone. The probability distributions in these two regimes are matched along a straight line: vertical bar V-R vertical bar = z*R. Our simulations confirm that this is indeed correct. We further obtain D-2 and z* as a function of the Stokes number, St. The former depends nonmonotonically on St with aminimum at about St approximate to 0.7 and the latter has only a weak dependence on St.

The effect of turbulence on the mass transfer between a fluid and embedded small heavy inertial particles that experience surface reactions is studied. For simplicity, the surface reaction, which takes place when a gas phase reactant is converted to a gas phase product at the external surface of the particles, is unimolar and isothermal. Two effects are identified. The first effect is due to the relative velocity between the fluid and the particles, and a model for the relative velocity is presented. The second effect is due to the clustering of particles, where the mass transfer rate is inhibited due to the rapid depletion of the consumed species inside the dense particle clusters. This last effect is relevant for large Damkohler numbers, where the Damkohler number is defined as the ratio of the turbulent and chemical time scales, and it may totally control the mass transfer rate for Damkohler numbers larger than unity. A model that describes how this effect should be incorporated into existing simulation tools that utilize the Reynolds averaged Navier-Stokes approach is presented.

We perform direct numerical simulations (DNS) of passive heavy inertial particles (dust) in homogeneous and isotropic two-dimensional turbulent flows (gas) for a range of Stokes number, St < 1. We solve for the particles using both a Lagrangian and an Eulerian approach (with a shock-capturing scheme). In the latter, the particles are described by a dust-density field and a dust-velocity field. We find the following: the dust-density field in our Eulerian simulations has the same correlation dimension d(2) as obtained from the clustering of particles in the Lagrangian simulations for St < 1; the cumulative probability distribution function of the dust density coarse grained over a scale r, in the inertial range, has a left tail with a power-law falloff indicating the presence of voids; the energy spectrum of the dust velocity has a power-law range with an exponent that is the same as the gas-velocity spectrum except at very high Fourier modes; the compressibility of the dust-velocity field is proportional to St(2). We quantify the topological properties of the dust velocity and the gas velocity through their gradient matrices, called A and B, respectively. Our DNS confirms that the statistics of topological properties of B are the same in Eulerian and Lagrangian frames only if the Eulerian data are weighed by the dust density. We use this correspondence to study the statistics of topological properties of A in the Lagrangian frame from our Eulerian simulations by calculating density-weighted probability distribution functions. We further find that in the Lagrangian frame, the mean value of the trace of A is negative and its magnitude increases with St approximately as exp(-C/St) with a constant C approximate to 0.1. The statistical distribution of different topological structures that appear in the dust flow is different in Eulerian and Lagrangian (density-weighted Eulerian) cases, particularly for St close to unity. In both of these cases, for small St the topological structures have close to zero divergence and are either vortical (elliptic) or strain dominated (hyperbolic, saddle). As St increases, the contribution to negative divergence comes mostly from saddles and the contribution to positive divergence comes from both vortices and saddles. Compared to the Eulerian case, the Lagrangian (density-weighted Eulerian) case has less outward spirals and more converging saddles. Inward spirals are the least probable topological structures in both cases.

We show, by direct numerical simulations, that heavy inertial particles (characterized by Stokes number St) in inhomogeneously forced statistically stationary isothermal turbulent flows cluster at the minima of mean-square turbulent velocity. Two turbulent transport processes, turbophoresis and turbulent diffusion together determine the spatial distribution of the particles. If the turbulent diffusivity is assumed to scale with turbulent root-mean-square velocity, as is the case for homogeneous turbulence, the turbophoretic coefficient can be calculated. Indeed, for the above assumption, the non-dimensional product of the turbophoretic coefficient and the rms velocity is shown to increase with St for small St, reach a maxima for St approximate to 10 and decrease as similar to St(-0.33) for large St.

The effect of turbulence on the heterogeneous (solid-fluid) reactions of solid particles is studied numerically with Direct Numerical Simulations (DNS). A simplified reaction system is used, where the solid-fluid reaction is represented by a single isothermal reaction step. It is found that, due to the clustering of particles by the isotropic turbulence, the overall reaction rate is entirely controlled by the turbulence for large Damkohler numbers. The particle clustering significantly slows down the reaction rate for increasing Damkohler numbers which reaches an asymptotic limit that can be analytically derived. This implies that the effect of turbulence on heterogeneously reacting particles should be included in models that are used in CFD simulations of e. g. char burnout in combustors or gasifiers. Such a model, based on the chemical and turbulent time scales, is here proposed for the heterogeneous reaction rate in the presence of turbulence.