Modal and non-modal linear stability analysis of channel flow with a dilute particle suspension is presented where particles are assumed to be solid, spherical, and heavy. The two-way coupling between particle and fluid flow is therefore modeled by the Stokes drag only. The results are presented as function of the particle relaxation time and mass fraction. First, we consider exponentially growing perturbations and extend previous findings showing the potential for a significant increase of the critical Reynolds number. The largest stabilization is observed when the ratio between the particle relaxation time and the oscillation period of the wave is of order one. By examining the energy budget, we show that this stabilization is due to the increase of the dissipation caused by the Stokes drag. The observed stabilization has led to the hypothesis that dusty flows can be more stable. However, transition to turbulence is most often subcritical in canonical shear flows where non-modal growth mechanisms are responsible for the initial growth of external disturbances. The non-modal analysis of the particle-laden flow, presented here for the first time, reveals that the transient energy growth is, surprisingly, increased by the presence of particles, in proportion to the particle mass fraction. The generation of streamwise streaks via the lift-up mechanism is still the dominant disturbance-growth mechanism in the particle laden flow; the length scales of the most dangerous disturbances are unaffected, while the initial disturbance growth can be delayed. These results are explained in terms of a dimensionless parameter relating the particle relaxation time to the time scale of the instability. The presence of a dilute solid phase therefore may not always work as a flow-control strategy for maintaining the flow as laminar. Despite the stabilizing effect on modal instabilities, non-modal mechanisms are still strong in internal flows seeded with heavy particles. Our results indicate that the initial stages of transition in dilute suspensions of small particles are similar to the stages in a single phase flow.

We present direct numerical simulations of subcritical transition to turbulence in a particle-laden channel flow, with particles assumed rigid, spherical, and heavier than the fluid. The equations describing the fluid flow are solved with an Eulerian mesh, whereas those describing the particle dynamics are solved by Lagrangian tracking. Two-way coupling between fluid and particles is modeled with Stokes drag. The numerical code is first validated against previous results from linear stability: the nonmodal growth of streamwise vortices resulting in streamwise streaks is still the most efficient mechanism for linear disturbance amplification at subcritical conditions as for the case of a single phase fluid. To analyze the full nonlinear transition, we examine two scenarios well studied in the literature: (1) transition initiated by streamwise independent counter-rotating streamwise vortices and one three-dimensional mode and (2) oblique transition, initiated by the nonlinear interaction of two symmetric oblique waves. The threshold energy for transition is computed, and it is demonstrated that for both scenarios the transition may be facilitated by the presence of particles at low number density. This is due to the fact that particles may introduce in the system detrimental disturbances of length scales not initially present. At higher concentrations, conversely, we note an increase of the disturbance energy needed for transition. The threshold energy for the oblique scenario shows a more significant increase in the presence of particles, by a factor about four. Interestingly, for the streamwise-vortex scenario the time at which transition occurs increases with the particle volume fraction when considering disturbances of equal initial energy. These results are explained by considering the reduced amplification of oblique modes in the two-phase flow. The results from these two classical scenarios indicate that, although linear stability analysis shows hardly any effect on optimal growth, particles do influence secondary instabilities and streak breakdown. These effects can be responsible of the reduced drag observed in turbulent channel flow laden with heavy particles.

No abstract