Collective waves in dense and confined microfluidic droplet arrays

The Stability of a Hydrodynamic Bravais Lattice

We present the results of a theoretical investigation of the stability and collective vibrations of a two-dimensional hydrodynamic lattice comprised of millimetric droplets bouncing on the surface of a vibrating liquid bath. We derive the linearized equations of motion describing the dynamics of a generic Bravais lattice, as encompasses all possible tilings of parallelograms in an infinite plane-filling array. Focusing on square and triangular lattice geometries, we demonstrate that for relatively low driving accelerations of the bath, only a subset of inter-drop spacings exist for which stable lattices may be achieved. The range of stable spacings is prescribed by the structure of the underlying wavefield. As the driving acceleration is increased progressively, the initially stationary lattices destabilize into coherent oscillatory motion. Our analysis yields both the instability threshold and the wavevector and polarization of the most unstable vibrational mode. The non-Markovian nature of the droplet dynamics renders the stability analysis of the hydrodynamic lattice more rich and subtle than that of its solid state counterpart.

Microfluidic formation of crystal-like structures

Crystal-like structures find application in several fields ranging from biomedical engineering to material science. For instance, droplet crystals are critical for high throughput assays and material synthesis, while particle crystals are important for particles and cell encapsulation, Drop-seq technologies, and single-cell analysis. Formation of crystal-like structures relies entirely on the possibility of manipulating with great accuracy the micrometer-size objects forming the crystal. In this context, microfluidic devices offer versatile tools for the precise manipulation of droplets and particles, thus enabling fabrication of crystal-like structures that form due to hydrodynamic interactions among droplets or particles. In this review, we aim at providing an holistic representation of crystal-like structure formation mediated by hydrodynamic interactions in microfluidic devices. We also discuss the physical origin of these hydrodynamic interactions and their relation to parameters such as device geometry, fluid properties, and flow conditions.

Universal motion of mirror-symmetric microparticles in confined Stokes flow

Comprehensive understanding of particle motion in microfluidic devices is essential to unlock additional technologies for shape-based separation and sorting of microparticles like microplastics, cells, and crystal polymorphs. Such particles interact hydrodynamically with confining surfaces, thus altering their trajectories. These hydrodynamic interactions are shape dependent and can be tuned to guide a particle along a specific path. We produce strongly confined particles with various shapes in a shallow microfluidic channel via stop flow lithography. Regardless of their exact shape, particles with a single mirror plane have identical modes of motion: in-plane rotation and cross-stream translation along a bell-shaped path. Each mode has a characteristic time, determined by particle geometry. Furthermore, each particle trajectory can be scaled by its respective characteristic times onto two master curves. We propose minimalistic relations linking these timescales to particle shape. Together these master curves yield a trajectory universal to particles with a single mirror plane.

Particle pairs and trains in inertial microfluidics

Staggered and linear multi-particle trains constitute characteristic structures in inertial microfluidics. Using lattice-Boltzmann simulations, we investigate their properties and stability, when flowing through microfluidic channels. We confirm the stability of cross-streamline pairs by showing how they contract or expand to their equilibrium axial distance. In contrast, same-streamline pairs quickly expand to a characteristic separation but even at long times slowly drift apart. We reproduce the distribution of particle distances with its characteristic peak as measured in experiments. Staggered multi-particle trains initialized with an axial particle spacing larger than the equilibrium distance contract non-uniformly due to collective drag reduction. Linear particle trains, similar to pairs, rapidly expand toward a value about twice the equilibrium distance of staggered trains and then very slowly drift apart non-uniformly. Again, we reproduce the statistics of particle distances and the characteristic peak observed in experiments. Finally, we thoroughly analyze the damped displacement pulse traveling as a microfluidic phonon through a staggered train and show how a defect strongly damps its propagation.

Flowing droplet interface bilayers: A microfluidic tool to control droplet trajectories and to study mechanical properties of unsupported lipid bilayers

We introduce the concept of Flowing Droplet Interface Bilayers (FDIBs) that are made of two droplets maintained in contact due to the presence of an adhesive lipidic surfactant. This system is similar to a flowing dumbbell made of two droplets interconnected by a lipid bilayer and driven by an external flow. Interestingly, such a dumbbell does not show a straight flow trajectory, but it oscillates between the sidewalls while moving along the microchannel. The origin of this unusual motion is hydrodynamic interactions, as demonstrated by analytical calculations and micro particle image velocimentry (µPiV) measurements. The hydrodynamic motion appears to be highly sensitive to the mechanical properties of the lipid bilayer connecting the two droplets (FDIB). Thus, droplet trajectories can be controlled by tuning the lipid bilayer composition, which enables in turn investigating mechanical properties of free-standing lipid bilayers.

Collective vibrations of a hydrodynamic active lattice

Recent experiments show that quasi-one-dimensional lattices of self-propelled droplets exhibit collective instabilities in the form of out-of-phase oscillations and solitary-like waves. This hydrodynamic lattice is driven by the external forcing of a vertically vibrating fluid bath, which invokes a field of subcritical Faraday waves on the bath surface, mediating the spatio-temporal droplet coupling. By modelling the droplet lattice as a memory-endowed system with spatially non-local coupling, we herein rationalize the form and onset of instability in this new class of dynamical oscillator. We identify the memory-driven instability of the lattice as a function of the number of droplets, and determine equispaced lattice configurations precluded by geometrical constraints. Each memory-driven instability is then classified as either a super- or subcritical Hopf bifurcation via a systematic weakly nonlinear analysis, rationalizing experimental observations. We further discover a previously unreported symmetry-breaking instability, manifest as an oscillatory-rotary motion of the lattice. Numerical simulations support our findings and prompt further investigations of this nonlinear dynamical system.

Simultaneous measurement of surface and bilayer tension in a microfluidic chip

Freestanding lipid bilayers are one of the most used model systems to mimic biological cell membranes. To form an unsupported bilayer, we employ two aqueous fingers in a microfluidic chip surrounded by an oily phase that contains lipids. Upon pushing two aqueous fingers forward, their interface becomes decorated with a lipid monolayer and eventually zip to form a bilayer when the monolayers have nanoscopic contact with each other. Using this straightforward approach, the quick and easy bilayer formation is facilitated by oil draining into the microfluidic device material consisting of polydimethylsiloxane. However, the oil drainage limits the lifetime of a bilayer to about 1 h. We demonstrate that this drainage can be managed, resulting in superior bilayer stability and an increased lifetime of several hours when using a pressure-controlled system. Applying different pressures to the aqueous fingers in the microfluidic chip, the formed bilayer can even be bent to a desired curvature. Extracting the contact angle and the resulting curvature of the bilayer region, for a given applied pressure difference, both the bilayer tension and the surface tension of each lipid monolayer can be derived from a single experiment using the Young Laplace pressure equation.

A flowing pair of particles in inertial microfluidics

A flowing pair of particles in inertial microfluidics gives important insights into understanding and controlling the collective dynamics of particles like cells or droplets in microfluidic devices. They are applied in medical cell analysis and engineering. We study the dynamics of a pair of solid particles flowing through a rectangular microchannel using lattice Boltzmann simulations. We determine the inertial lift force profiles as a function of the two particle positions, their axial distance, and the Reynolds number. Generally, the profiles strongly differ between particles leading and lagging in flow and the lift forces are enhanced due to the presence of a second particle. At small axial distances, they are determined by viscous forces, while inertial forces dominate at large separations. We identify cross-streamline pairs as stable fixed points in the lift force profiles and argue that same-streamline configurations are only one-sided stable. Depending on the initial conditions, the two-particle lift forces in combination with the Poiseuille flow give rise to three types of unbound particle trajectories, called moving-apart, passing, and swapping, and one type of bound trajectory, where the particles perform damped oscillations towards the cross-stream line configuration. The damping rate scales with Reynolds number squared, since inertial forces are responsible for driving the particles to their steady-state positions.

Small asymmetric Brownian objects self-align in nanofluidic channels

Although the self-alignment of asymmetric macro-sized objects of a few tens of microns in size have been studied extensively in experiments and theory, access to much smaller length scales is still hindered by technical challenges. We combine molecular dynamics and stochastic rotation dynamics techniques to investigate the self-orientation phenomenon at different length scales, ranging from the micron to the nano scale by progressively increasing the relative strength of diffusion over convection. To this end, we model an asymmetric dumbbell particle in Hele-Shaw flow and explore a wide range of Péclet numbers (Pe) and different particle shapes, as characterized by the size ratio of the two dumbbell spheres (R[combining tilde]). By independently varying these two parameters we analyse the process of self-orientation and characterize the alignment of the dumbbell with the direction of the fluid flow. We identify three different regimes of strong, weak and no alignment and we map out a state diagram in Pe versus R[combining tilde] plane. Based on these results, we estimate dimensional length scales and flow rates for which these findings would be applicable in experiments. Finally, we find that the characteristic reorientation time of the dumbbell is a monotonically decreasing function of the dumbbell anisotropy.

Activity-induced instability of phonons in 1D microfluidic crystals

One-dimensional crystals of passively-driven particles in microfluidic channels exhibit collective vibrational modes reminiscent of acoustic 'phonons'. These phonons are induced by the long-range hydrodynamic interactions among the particles and are neutrally stable at the linear level. Here, we analyze the effect of particle activity - self-propulsion - on the emergence and stability of these phonons. We show that the direction of wave propagation in active crystals is sensitive to the intensity of the background flow. We also show that activity couples, at the linear level, transverse waves to the particles' rotational motion, inducing a new mode of instability that persists in the limit of large background flow, or, equivalently, vanishingly small activity. We then report a new phenomenon of phonons switching back and forth between two adjacent crystals in both passively-driven and active systems, similar in nature to the wave switching observed in quantum mechanics, optical communication, and density stratified fluids. These findings could have implications for the design of commercial microfluidic systems and the self-assembly of passive and active micro-particles into one-dimensional structures.

Mesoscopic modelling and simulation of soft matter

The deformability of soft condensed matter often requires modelling of hydrodynamical aspects to gain quantitative understanding. This, however, requires specialised methods that can resolve the multiscale nature of soft matter systems. We review a number of the most popular simulation methods that have emerged, such as Langevin dynamics, dissipative particle dynamics, multi-particle collision dynamics, sometimes also referred to as stochastic rotation dynamics, and the lattice-Boltzmann method. We conclude this review with a short glance at current compute architectures for high-performance computing and community codes for soft matter simulation.

Two-dimensional flow of driven particles: a microfluidic pathway to the non-equilibrium frontier

We discuss the basic physics of the flow of micron-scale droplets in 2D geometry.