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.
Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, May 14 until 2:00 AM ET on Friday, May 15 due to maintenance. We apologize for the inconvenience.


Title: Lagrangian investigation of the interface dynamics in single-mode Rayleigh–Taylor instability
We apply Lagrangian particle tracking to the two-dimensional single-mode Rayleigh–Taylor (RT) instability to study the dynamical evolution of fluid interface. At the onset of the nonlinear RT stage, we select three ensembles of tracer particles located at the bubble tip, at the spike tip, and inside the spiral of the mushroom structure, which cover most of the interfacial region as the instability develops. Conditional statistics performed on the three sets of particles and over different RT evolution stages, such as the trajectory curvature, velocity, and acceleration, reveals the temporal and spatial flow patterns characterizing the single-mode RT growth. The probability density functions of tracer particle velocity and trajectory curvature exhibit scalings compatible with local flow topology, such as the swirling motion of the spiral particles. Large-scale anisotropy of RT interfacial flows, measured by the ratio of horizontal to vertical kinetic energy, also varies for different particle ensembles arising from the differing evolution patterns of the particle acceleration. In addition, we provide direct evidence to connect the RT bubble re-acceleration to its interaction with the transported fluid from the spike side, due to the shear driven Kelvin–Helmholtz instability. Furthermore, we reveal that the secondary RT instability inside the spiral, which destabilizes the spiraling motion and induces complex flow structures, is generated by the centrifugal acceleration.  more » « less
Award ID(s):
2020249
PAR ID:
10519055
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
aip.org
Date Published:
Journal Name:
Physics of Fluids
Volume:
35
Issue:
10
ISSN:
1070-6631
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Sharp edge structures have been demonstrated as an efficient way of generating acoustic streaming in microfluidic devices, which finds numerous applications in fluid mixing, pumping, particle actuation, and cell lysis. A sharp tip capillary is widely available means of generating sharp structures without the need of microfabrication, which has been used for studying enzyme kinetics, droplet digital PCR, and mass spectrometry analysis. In this work, we studied the influence of liquid inside the vibrating glass capillary on its efficiency of generating acoustic streaming. Using fluorescence microscopy and fluorescent particles, we observed that adding liquid to the inside of the vibrating glass capillary changed the streaming patterns as well as led to increased streaming velocity. Based on the observed streaming patterns, we hypothesized the liquid present in the capillary changed vibration mode of the capillary, which is matched with COMSOL simulations. Finally, the utility of the liquid filled vibrating capillary was demonstrated for higher energy efficiency for fluid mixing and mass spectrometry experiments. This study will provide useful guidance when optimizing the efficiency of vibrating sharp tip capillary systems. 
    more » « less
  2. The encapsulation of active particles, such as bacteria or active colloids, inside a droplet gives rise to a non-trivial shape dynamics and droplet displacement. To understand this behaviour, we derive an asymptotic solution for the fluid flow about a deformable droplet containing an active particle, modelled as a Stokes-flow singularity, in the case of small shape distortions. We develop a general solution for any Stokes singularity and apply it to compute the flows and resulting droplet velocity due to common singularity representations of active particles, such as Stokeslets, rotlets and stresslets. The results show that offsetting of the active particle from the centre of the drop breaks symmetry and excites a large number of generally non-axisymmetric shape modes as well as particle and droplet motion. In the case of a swimming stresslet singularity, a run-and-tumble locomotion results in superdiffusive droplet displacement. The effect of interfacial properties is also investigated. Surfactants adsorbed at the droplet interface counteract the internal flow and arrest the droplet motion for all Stokes singularities except the Stokeslet. Our results highlight strategies to steer the flows of active particles and create autonomously navigating containers. 
    more » « less
  3. Abstract The Rayleigh–Taylor (RT) instability is ubiquitously observed, yet has traditionally been studied using ideal fluid models. Collisionality can vary strongly across the fluid interface, and previous work demonstrates the necessity of kinetic models to completely capture dynamics in certain collisional regimes. Where previous kinetic simulations used spatially and temporally constant collision frequency, this work presents five-dimensional (two spatial, three velocity dimensions) continuum-kinetic simulations of the RT instability using a more realistic spatially varying collision frequency. Three cases of collisional variation are explored for two Atwood numbers: low to intermediate, intermediate to high, and low to high. The low-to-intermediate case exhibits no RT instability growth, while the intermediate-to-high case is similar to a fluid-limit kinetic case with interface widening biased toward the lower-collisionality region. A novel contribution of this work is the low-to-high collisionality case that shows significantly altered instability growth through an upward movement of the interface and damped spike growth due to increased free-streaming particle diffusion in the lower region. Contributions to the energy flux from the non-Maxwellian portions of the distribution function are not accessible to fluid models and are greatest in magnitude in the spike and regions of low collisionality. Increasing the Atwood number results in greater RT instability growth and reduced upward interface movement. Deviation of the distribution function from Maxwellian is inversely proportional to collision frequency and concentrated around the fluid interface. The linear phase of RT instability growth is well described by theoretical linear growth rates accounting for viscosity and diffusion. 
    more » « less
  4. Particle-laden flows of sedimenting solid particles or droplets in a carrier gas have strong inter-phase coupling. Even at low particle volume fractions, the two-way coupling can be significant due to the large particle to gas density ratio. In this semi-dilute regime, the slip velocity between phases leads to sustained clustering that strongly modulates the overall flow. The analysis of perturbations in homogeneous shear reveals the process by which clusters form: (i) the preferential concentration of inertial particles in the stretching regions of the flow leads to the formation of highly concentrated particle sheets, (ii) the thickness of the latter is controlled by particle-trajectory crossing, which causes a local dispersion of particles, (iii) a transverse Rayleigh–Taylor instability, aided by the shear-induced rotation of the particle sheets towards the gravity normal direction, breaks the planar structure into smaller clusters. Simulations in the Euler–Lagrange formalism are compared to Euler–Euler simulations with the two-fluid and anisotropic-Gaussian methods. It is found that the two-fluid method is unable to capture the particle dispersion due to particle-trajectory crossing and leads instead to the formation of discontinuities. These are removed with the anisotropic-Gaussian method which derives from a kinetic approach with particle-trajectory crossing in mind. 
    more » « less
  5. The Lagrangian particle method is widely used to understand scalar tracer concentration fields in models of the atmosphere and oceans. Simulating virtual particles provides an alternative description of advection to the Eulerian representation in models and aids in identifying pathways, timescales, and connectivity. Atmospheric and oceanic models solve advection‐diffusion‐reaction equations to simulate tracers, in which only the advective component is captured by traditional Lagrangian approaches. In this work, we report a novel method that closes tracer budgets on Lagrangian trajectories in a manner consistent with Eulerian budgets in finite‐volume models. The scalar tracer concentrations on grid cell walls are derived from the model advection scheme and then interpolated inside grid boxes along streamlines. The divergence of the diffusive flux and reaction terms are interpolated based on velocity and tracer concentration, ensuring the tracer budget closes in terms of both trajectory and volume integrals. Compared to the Eulerian budget analysis, which considers a fixed volume, our method quantifies the tracer evolution within a volume that moves along with the flow. We demonstrate the method using a case study of Southern Ocean biogeochemistry. Another case study involves analyzing the heat budget of the 2011 Western Australian marine heat wave. The method bridges the gap between Eulerian budget and Lagrangian particle analyses by representing the advective processes with particle movements and interpolating the diffusive and reactive processes onto trajectories in a way consistent with the finite‐volume description. 
    more » « less