skip to main content


Title: Preferential concentration driven instability of sheared gas–solid suspensions
We examine the linear stability of a homogeneous gas–solid suspension of small Stokes number particles, with a moderate mass loading, subject to a simple shear flow. The modulation of the gravitational force exerted on the suspension, due to preferential concentration of particles in regions of low vorticity, in response to an imposed velocity perturbation, can lead to an algebraic instability. Since the fastest growing modes have wavelengths small compared with the characteristic length scale ( $U_{g}/{\it\Gamma}$ ) and oscillate with frequencies large compared with ${\it\Gamma}$ , $U_{g}$ being the settling velocity and ${\it\Gamma}$ the shear rate, we apply the WKB method, a multiple scale technique. This analysis reveals the existence of a number density mode which travels due to the settling of the particles and a momentum mode which travels due to the cross-streamline momentum transport caused by settling. These modes are coupled at a turning point which occurs when the wavevector is nearly horizontal and the most amplified perturbations are those in which a momentum wave upstream of the turning point creates a downstream number density wave. The particle number density perturbations reach a finite, but large amplitude that persists after the wave becomes aligned with the velocity gradient. The growth of the amplitude of particle concentration and fluid velocity disturbances is characterised as a function of the wavenumber and Reynolds number ( $\mathit{Re}=U_{g}^{2}/{\it\Gamma}{\it\nu}$ ) using both asymptotic theory and a numerical solution of the linearised equations.  more » « less
Award ID(s):
1233793
NSF-PAR ID:
10025838
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
770
ISSN:
0022-1120
Page Range / eLocation ID:
85 to 123
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The standard theory of pulsations deals with the frequencies and growth rates of infinitesimal perturbations in a stellar model. Modes which are calculated to be linearly driven should increase their amplitudes exponentially with time; the fact that nearly constant amplitudes are usually observed is evidence that nonlinear mechanisms inhibit the growth of finite amplitude pulsations. Models predict that the mass of DAV convection zones is very sensitive to temperature (i.e., MCZ∝T−90eff) leading to the possibility that even "small amplitude" pulsators may experience significant nonlinear effects. In particular, the outer turning point of finite-amplitude g-mode pulsations can vary with the local surface temperature, producing a reflected wave that is slightly out of phase with that required for a standing wave. This can lead to a lack of coherence of the mode and a reduction in its global amplitude. We compute the size of this effect for specific examples and discuss the results in the context of Kepler and K2 observations. 
    more » « less
  2. We investigate wind wave growth by direct numerical simulations solving for the two-phase Navier–Stokes equations. We consider the ratio of the wave speed $c$ to the wind friction velocity $u_*$ from $c/u_*= 2$ to 8, i.e. in the slow to intermediate wave regime; and initial wave steepness $ak$ from 0.1 to 0.3; the two being varied independently. The turbulent wind and the travelling, nearly monochromatic waves are fully coupled without any subgrid-scale models. The wall friction Reynolds number is 720. The novel fully coupled approach captures the simultaneous evolution of the wave amplitude and shape, together with the underwater boundary layer (drift current), up to wave breaking. The wave energy growth computed from the time-dependent surface elevation is in quantitative agreement with that computed from the surface pressure distribution, which confirms the leading role of the pressure forcing for finite amplitude gravity waves. The phase shift and the amplitude of the principal mode of surface pressure distribution are systematically reported, to provide direct evidence for possible wind wave growth theories. Intermittent and localised airflow separation is observed for steep waves with small wave age, but its effect on setting the phase-averaged pressure distribution is not drastically different from that of non-separated sheltering. We find that the wave form drag force is not a strong function of wave age but closely related to wave steepness. In addition, the history of wind wave coupling can affect the wave form drag, due to the wave crest shape and other complex coupling effects. The normalised wave growth rate we obtain agrees with previous studies. We make an effort to clarify various commonly adopted underlying assumptions, and to reconcile the scattering of the data between different previous theoretical, numerical and experimental results, as we revisit this longstanding problem with new numerical evidence. 
    more » « less
  3. A theoretical study is made of the stability of propagating internal gravity wave modes along a horizontal stratified fluid layer bounded by rigid walls. The analysis is based on the Floquet eigenvalue problem for infinitesimal perturbations to a wave mode of small amplitude. The appropriate instability mechanism hinges on how the perturbation spatial scale relative to the basic-state wavelength, controlled by a parameter $\mu$ , compares to the basic-state amplitude parameter, $\epsilon \ll 1$ . For $\mu ={O}(1)$ , the onset of instability arises due to perturbations that form resonant triads with the underlying wave mode. For short-scale perturbations such that $\mu \ll 1$ but $\alpha =\mu /\epsilon \gg 1$ , this triad resonance instability reduces to the familiar parametric subharmonic instability (PSI), where triads comprise fine-scale perturbations with half the basic-wave frequency. However, as $\mu$ is further decreased holding $\epsilon$ fixed, higher-frequency perturbations than these two subharmonics come into play, and when $\alpha ={O}(1)$ Floquet modes feature broadband spectrum. This broadening phenomenon is a manifestation of the advection of small-scale perturbations by the basic-wave velocity field. By working with a set of ‘streamline coordinates’ in the frame of the basic wave, this advection can be ‘factored out’. Importantly, when $\alpha ={O}(1)$ PSI is replaced by a novel, multi-mode resonance mechanism which has a stabilising effect that provides an inviscid short-scale cut-off to PSI. The theoretical predictions are supported by numerical results from solving the Floquet eigenvalue problem for a mode-1 basic state. 
    more » « less
  4. Particulate matter in the environment, such as sediment, marine debris and plankton, is transported by surface waves. The transport of these inertial particles is different from that of fluid parcels described by Stokes drift. In this study, we consider the transport of negatively buoyant particles that settle in flow induced by surface waves as described by linear wave theory in arbitrary depth. We consider particles that fall under both a linear drag regime in the low Reynolds number limit and in a nonlinear drag regime in the transitional Reynolds number range. Based on an analysis of typical applications, we find that the nonlinear regime is the most widely applicable. From an expansion in the particle Stokes number, we find kinematic expressions for inertial particle motion in waves, and from a multiscale expansion in the dimensionless wave amplitude, we find expressions for the wave-averaged drift velocities. These drift velocities are analogous to Stokes drift and can be used in large-scale models that do not resolve surface waves. We find that the horizontal drift velocity is reduced relative to the Stokes drift of fluid parcels and that the vertical drift velocity is enhanced relative to the particle terminal settling velocity. We also demonstrate that a cloud of settling particles released simultaneously will disperse in the horizontal direction. Finally, we discuss the accuracy of our expressions by comparing against numerical simulations, which show excellent agreement, and against experimental data, which show the same trends. 
    more » « less
  5. Abstract

    The dynamics of geophysical dilute turbulent gas‐particles mixtures depends to a large extent on particle concentration, which in turn depends predominantly on the particle settling velocity. We experimentally investigate air‐particle mixtures contained in a vertical pipe in which the velocity of an ascending air flux matches the settling velocity of glass particles. To obtain local particle concentrations in these mixtures, we use acoustic probing and air pressure measurements and show that these independent techniques yield similar results for a range of particle sizes and particle concentrations. Moreover, we find that in suspensions of small particles (78 μm) the settling velocity increases with the local particle concentration due to the formation of particle clusters. These clusters settle with a velocity that is four times faster than the terminal settling velocity of single particles, and they double settling speeds of the suspensions. In contrast, in suspensions of larger particles (467 μm) the settling velocity decreases with increasing particle concentration. Although particle clusters are still present in this case, the settling velocity is decreased by 30%, which is captured by a hindered settling model. These results suggest an interplay between hindered settling and cluster‐induced enhanced settling, which in our experiments occur respectively at Stokes number O(100) and O(1). We discuss implications for volcanic plumes and pyroclastic currents. Our study suggests that clustering and related enhanced or hindered particle settling velocities should be considered in models of volcanic phenomena and that drag law corrections are needed for reliable predictions and hazard assessment.

     
    more » « less