Governing equations
The starting point for a mathematical description of almost any type of fluid flow is the classical set of Navier–Stokes equations. To describe particle-laden flows, we must modify these equations to account for the effect of the particles on the carrier, or vice versa, or both - a suitable choice of such added complications depend on a variety of the parameters, for instance, how dense the particles are, how concentrated they are, or whether or not they are chemically reactive. In most real world cases, the particles are very small and occur in low concentrations, hence the dynamics are governed primarily by the continuous phase. A possible way to represent the dynamics of the carrier phase is by the following modified Navier-Stokes momentum equation: : where is a momentum source or sink term, arising from the presence of the particle phase. The above equation is an Eulerian equation, that is, the dynamics are understood from the viewpoint of a fixed point in space. The dispersed phase is typically (though not always) treated in a Lagrangian framework, that is, the dynamics are understood from the viewpoint of fixed particles as they move through space. A usual choice of momentum equation for a particle is: : where represents the carrier phase velocity and represents the particle velocity. is the particle relaxation time, and represents a typical timescale of the particle's reaction to changes in the carrier phase velocity - loosely speaking, this can be thought of as the particle's inertia with respect to the fluid with contains it. The interpretation of the above equation is that particle motion is hindered by a drag force. In reality, there are a variety of other forces which act on the particle motion (such as gravity, Basset history and added mass) – as described through for instance the Basset–Boussinesq–Oseen equation. However, for many physical examples, in which the density of the particle far exceeds the density of the medium, the above equation is sufficient. A typical assumption is that the particles are spherical, in which case the drag is modeled usingCoupling
If the mass fraction of the dispersed phase is small, then ''one-way coupling'' between the phases is a reasonable assumption; that is, the dynamics of the particle phase are affected by the carrier phase, but the reverse is not the case. However, if the mass fraction of the dispersed phase is large, the interaction of the dynamics between the two phases must be considered - this is ''two-way coupling''. A problem with the Lagrangian treatment of the dispersed phase is that once the number of particles becomes large, it may require a prohibitive amount of computational power to track a sufficiently large sample of particles required for statistical convergence. In addition, if the particles are sufficiently light, they behave essentially like a second fluid. In this case, an Eulerian treatment of the dispersed phase is sensible.Modeling
Like all fluid dynamics-related disciplines, the modelling of particle-laden flows is an enormous challenge for researchers - this is because most flows of practical interest arePreferential migration
An interesting aspect of particle-laden flows is preferential migration of particles to certain regions within the fluid flow. This is often characterized by the Stokes number (St) of the particles. At low St, particles tend to act as tracers and are uniformly distributed. At high St, particles are heavy and are influenced less by the fluid and more by its inertia. At intermediate St, particles are affected by both the fluid motion and its inertia, which gives rise to several interesting behaviors. This is especially true in wall-bounded flows where there is a velocity gradient near the wall. One of the earliest works describing preferential migration is the experimental work of Segre and Silberberg. They showed that a neutrally buoyant particle in a laminar pipe flow comes to an equilibrium position between the wall and the axis. This is referred to as the Segré–Silberberg effect. Saffman explained this in terms of the force acting on the particle when it experiences a velocity gradient across it. Feng et al. have studied this through detailed direct numerical simulations and have elaborated on the physical mechanism of this migration. Recently it was found that even for non-neutrally buoyant particles similar preferential migration occurs . At low St, the particles tend to settle at an equilibrium position while for high St, the particles begin to oscillate about the center of the channel. The behavior becomes even interesting in turbulent flows. Here, the turbophoretic force (transport of particles down gradients of turbulent kinetic energy) causes a high concentration of particles near the walls. Experimental and particle-resolved DNS studies have explained the mechanism of this migration in terms of the Saffman lift and the turbophoretic force . These preferential migration are of significant importance to several applications where wall-bounded particle-laden flows are encountered and is an active area of research.Further reading
* Mashayek, F. and Pandya, R. V. R. (1921), ''Progress in Energy and Combustion Science'' 20, 196, 196–212. * Jebakumar, A. S., Premnath K. N. and Abraham J. (2016), ''Computers and Fluids'' 124, 208, 208–219.References
{{Reflist Fluid mechanics Fluid dynamics