skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Instability of a dusty vortex
We investigate the effect of inertial particles dispersed in a circular patch of finite radius on the stability of a two-dimensional Rankine vortex in semi-dilute dusty flows. Unlike the particle-free case where no unstable modes exist, we show that the feedback force from the particles triggers a novel instability. The mechanisms driving the instability are characterized using linear stability analysis for weakly inertial particles and further validated against Eulerian–Lagrangian simulations. We show that the particle-laden vortex is always unstable if the mass loading $M>0$ . Surprisingly, even non-inertial particles destabilize the vortex by a mechanism analogous to the centrifugal Rayleigh–Taylor instability in radially stratified vortex with density jump. We identify a critical mass loading above which an eigenmode $$m$$ becomes unstable. This critical mass loading drops to zero as $$m$$ increases. When particles are inertial, modes that fall below the critical mass loading become unstable, whereas modes above it remain unstable but with lower growth rates compared with the non-inertial case. Comparison with Eulerian–Lagrangian simulations shows that growth rates computed from simulations match well the theoretical predictions. Past the linear stage, we observe the emergence of high-wavenumber modes that turn into spiralling arms of concentrated particles emanating out of the core, while regions of particle-free flow are sucked inward. The vorticity field displays a similar pattern which leads to the breakdown of the initial Rankine structure. This novel instability for a dusty vortex highlights how the feedback force from the disperse phase can induce the breakdown of an otherwise resilient vortical structure.  more » « less
Award ID(s):
2148710
PAR ID:
10435744
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
948
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We investigate the effect of particle inertia on the merger of co-rotating dusty vortex pairs at semi-dilute concentrations. In a particle-free flow, the merger is triggered once the ratio of vortex core size to vortex separation reaches a critical value. The vortex pair separation then decreases monotonically until the two cores merge together. Using Eulerian–Lagrangian simulations of co-rotating particle-laden vortices, we show substantial departure from the vortex dynamics previously established in particle-free flows. Most strikingly, we find that disperse particles with moderate inertia cause the vortex pair to push apart to a separation nearly twice as large as the initial separation. During this stage, the drag force exerted by particles ejected out of the vortex cores on the fluid results in a net repulsive force that pushes the two cores apart. Eventually, the two dusty vortices merge into a single vortex with most particles accumulating outside the core, similar to the dusty Lamb–Oseen vortex described in Shuai & Kasbaoui (J. Fluid Mech., vol 936, 2022, p. A8). For weakly inertial particles, we find that the merger dynamics follows the same mechanics as that of a single-phase flow, albeit with a density that must be adjusted to match the mixture density. For highly inertial particles, the feedback force exerted by the particles on the fluid may stretch the two cores during the merger to a point where each core splits into two, resulting in inner and outer vortex pairs. In this case, the merger occurs in two stages where the inner vortices merge first, followed by the outer ones. 
    more » « less
  2. We study the instability of a dusty simple shear flow where the dust particles are distributed non-uniformly. A simple shear flow is modally stable to infinitesimal perturbations. Also, a band of particles remains unaffected in the absence of any background flow. However, we demonstrate that the combined scenario – comprising a simple shear flow with a localized band of particles – can exhibit destabilization due to their two-way interaction. The instability originates solely from the momentum feedback from the particle phase to the fluid phase. Eulerian–Lagrangian simulations are employed to illustrate the existence of this instability. Furthermore, the results are compared with a linear stability analysis of the system using an Eulerian–Eulerian model. Our findings indicate that the instability has an inviscid origin and is characterized by a critical wavelength below which it is not persistent. We have observed that increasing particle inertia dampens the unstable modes, whereas the strength of the instability increases with the strength of the coupling between the fluid and particle phases. 
    more » « less
  3. We investigate the effect of inertial particles on the stability and decay of a prototypical vortex tube, represented by a two-dimensional Lamb–Oseen vortex. In the absence of particles, the strong stability of this flow makes it resilient to perturbations, whereby vorticity and enstrophy decay at a slow rate controlled by viscosity. Using Eulerian–Lagrangian simulations, we show that the dispersion of semidilute inertial particles accelerates the decay of the vortex tube by orders of magnitude. In this work, mass loading is unity, ensuring that the fluid and particle phases are tightly coupled. Particle inertia and vortex strength are varied to yield Stokes numbers 0.1–0.4 and circulation Reynolds numbers 800–5000. Preferential concentration causes these inertial particles to be ejected from the vortex core forming a ring-shaped cluster and a void fraction bubble that expand outwards. The outward migration of the particles causes a flattening of the vorticity distribution, which enhances the decay of the vortex. The latter is further accelerated by small-scale clustering that causes enstrophy to grow, in contrast with the monotonic decay of enstrophy in single-phase two-dimensional vortices. These dynamics unfold on a time scale that is set by preferential concentration and is two orders of magnitude lower than the viscous time scale. Increasing particle inertia causes a faster decay of the vortex. This work shows that the injection of inertial particles could provide an effective strategy for the control and suppression of resilient vortex tubes. 
    more » « less
  4. We present the Moving Eulerian-Lagrangian Particles (MELP), a novel mesh-free method for simulating incompressible fluid on thin films and foams. Employing a bi-layer particle structure, MELP jointly simulates detailed, vigorous flow and large surface deformation at high stability and efficiency. In addition, we design multi-MELP: a mechanism that facilitates the physically-based interaction between multiple MELP systems, to simulate bubble clusters and foams with non-manifold topological evolution. We showcase the efficacy of our method with a broad range of challenging thin film phenomena, including the Rayleigh-Taylor instability across double-bubbles, foam fragmentation with rim surface tension, recovery of the Plateau borders, Newton black films, as well as cyclones on bubble clusters. 
    more » « less
  5. Summary In this paper, a three‐dimensional numerical solver is developed for suspensions of rigid and soft particles and droplets in viscoelastic and elastoviscoplastic (EVP) fluids. The presented algorithm is designed to allow for the first time three‐dimensional simulations of inertial and turbulent EVP fluids with a large number particles and droplets. This is achieved by combining fast and highly scalable methods such as an FFT‐based pressure solver, with the evolution equation for non‐Newtonian (including EVP) stresses. In this flexible computational framework, the fluid can be modeled by either Oldroyd‐B, neo‐Hookean, FENE‐P, or Saramito EVP models, and the additional equations for the non‐Newtonian stresses are fully coupled with the flow. The rigid particles are discretized on a moving Lagrangian grid, whereas the flow equations are solved on a fixed Eulerian grid. The solid particles are represented by an immersed boundary method with a computationally efficient direct forcing method, allowing simulations of a large numbers of particles. The immersed boundary force is computed at the particle surface and then included in the momentum equations as a body force. The droplets and soft particles on the other hand are simulated in a fully Eulerian framework, the former with a level‐set method to capture the moving interface and the latter with an indicator function. The solver is first validated for various benchmark single‐phase and two‐phase EVP flow problems through comparison with data from the literature. Finally, we present new results on the dynamics of a buoyancy‐driven drop in an EVP fluid. 
    more » « less