Blaise Delmotte

Unstable fronts and motile structures formed by microrollers

Collections of rolling colloids are shown to pinch off into motile clusters resembling droplets sliding down a windshield. These stable dynamic structures are formed through a fingering instability that relies on hydrodynamic interactions alone. Condensation of objects into stable clusters occurs naturally in equilibrium1 and driven systems2,3,4,5. It is commonly held that potential interactions6, depletion forces7, or sensing8 are the only mechanisms which can create long-lived compact structures. Here we show that persistent motile structures can form spontaneously from hydrodynamic interactions alone, with no sensing or potential interactions. We study this structure formation in a system of colloidal rollers suspended and translating above a floor, using both experiments and large-scale three-dimensional simulations. In this system, clusters originate from a previously unreported fingering instability, where fingers pinch off from an unstable front to form autonomous ‘critters’, whose size is selected by the height of the particles above the floor. These critters are a stable state of the system, move much faster than individual particles, and quickly respond to a changing drive. With speed and direction set by a rotating magnetic field, these active structures offer interesting possibilities for guided transport, flow generation, and mixing at the microscale.

Hydrodynamics of suspensions of passive and active rigid particles: a rigid multiblob approach

Author(s): Usabiaga, FB; Kallemov, B; Delmotte, B; Bhalla, APS; Griffith, BE; Donev, A | Abstract:

Large-scale simulation of steady and time-dependent active suspensions with the force-coupling method

We present a new development of the force-coupling method (FCM) to address the accurate simulation of a large number of interacting micro-swimmers. Our approach is based on the squirmer model, which we adapt to the FCM framework, resulting in a method that is suitable for simulating semi-dilute squirmer suspensions. Other effects, such as steric interactions, are considered with our model. We test our method by comparing the velocity field around a single squirmer and the pairwise interactions between two squirmers with exact solutions to the Stokes equations and results given by other numerical methods. We also illustrate our methods ability to describe spheroidal swimmer shapes and biologically-relevant time-dependent swimming gaits. We detail the numerical algorithm used to compute the hydrodynamic coupling between a large collection (104-105) of micro-swimmers. Using this methodology, we investigate the emergence of polar order in a suspension of squirmers and show that for large domains, both the steady-state polar order parameter and the growth rate of instability are independent of system size. These results demonstrate the effectiveness of our approach to achieve near continuum-level results, allowing for better comparison with experimental measurements while complementing and informing continuum models.

A general formulation of Bead Models applied to flexible fibers and active filaments at low Reynolds number

This contribution provides a general framework to use Lagrange multipliers for the simulation of low Reynolds number fiber dynamics based on Bead Models (BM). This formalism provides an efficient method to account for kinematic constraints. We illustrate, with several examples, to which extent the proposed formulation offers a flexible and versatile framework for the quantitative modeling of flexible fibers deformation and rotation in shear flow, the dynamics of actuated filaments and the propulsion of active swimmers. Furthermore, a new contact model called Gears Model is proposed and successfully tested. It avoids the use of numerical artifices such as repulsive forces between adjacent beads, a source of numerical difficulties in the temporal integration of previous Bead Models.

Simulating Brownian suspensions with fluctuating hydrodynamics.

Fluctuating hydrodynamics has been successfully combined with several computational methods to rapidly compute the correlated random velocities of Brownian particles. In the overdamped limit where both particle and fluid inertia are ignored, one must also account for a Brownian drift term in order to successfully update the particle positions. In this paper, we present an efficient computational method for the dynamic simulation of Brownian suspensions with fluctuating hydrodynamics that handles both computations and provides a similar approximation as Stokesian Dynamics for dilute and semidilute suspensions. This advancement relies on combining the fluctuating force-coupling method (FCM) with a new midpoint time-integration scheme we refer to as the drifter-corrector (DC). The DC resolves the drift term for fluctuating hydrodynamics-based methods at a minimal computational cost when constraints are imposed on the fluid flow to obtain the stresslet corrections to the particle hydrodynamic interactions. With the DC, this constraint needs only to be imposed once per time step, reducing the simulation cost to nearly that of a completely deterministic simulation. By performing a series of simulations, we show that the DC with fluctuating FCM is an effective and versatile approach as it reproduces both the equilibrium distribution and the evolution of particulate suspensions in periodic as well as bounded domains. In addition, we demonstrate that fluctuating FCM coupled with the DC provides an efficient and accurate method for large-scale dynamic simulation of colloidal dispersions and the study of processes such as colloidal gelation.

Brownian dynamics of confined suspensions of active microrollers

We develop efficient numerical methods for performing many-body Brownian dynamics simulations of a recently observed fingering instability in an active suspension of colloidal rollers sedimented above a wall [M. Driscoll, B. Delmotte, M. Youssef, S. Sacanna, A. Donev, and P. Chaikin, Nat. Phys. (2016), preprint arXiv:1609.08673. We present a stochastic Adams-Bashforth integrator for the equations of Brownian dynamics, which has the same cost but is more accurate than the widely used Euler-Maruyama scheme, and use a random finite difference to capture the stochastic drift proportional to the divergence of the configuration-dependent mobility matrix. We generate the Brownian increments using a Krylov method and show that for particles confined to remain in the vicinity of a no-slip wall by gravity or active flows, the number of iterations is independent of the number of particles. Our numerical experiments with active rollers show that the thermal fluctuations set the characteristic height of the colloids above the wall, both in the initial condition and the subsequent evolution dominated by active flows. The characteristic height in turn controls the time scale and wavelength for the development of the fingering instability.

Minimal model for a hydrodynamic fingering instability in microroller suspensions

We derive a minimal continuum model to investigate the hydrodynamic mechanism behind the fingering instability recently discovered in a suspension of microrollers near a floor [Driscoll et al. Nature Physics, 2016]. Our model, consisting of two continuous lines of rotlets, exhibits a linear instability driven only by hydrodynamics interactions, and reproduces the lengthscale selection observed in large scale particle simulations and in experiments. By adjusting only one parameter, the distance between the two lines, our dispersion relation exhibits quantitative agreement with the simulations and qualitative agreement with experimental measurements. Our linear stability analysis indicate that this instability is caused by the combination of the advective and transverse flows generated by the microrollers near a no-slip surface. Our simple model offers an interesting formalism to characterize other hydrodynamic instabilities that have not been yet well understood, such as size scale selection in suspensions of particles sedimenting adjacent to a wall, or the recently observed formations of traveling phonons in systems of confined driven particles.

Hydrodynamic interactions among large populations of swimming micro-organisms

Bio-convection, biofilm forming or mechanics of reproduction are connected to the motility and collective behaviour of micro-organisms. For instance, spermatozoa suspensions exhibit coherent motion whose frequency and lifespan are strongly correlated to semen fertility (Moore et al. 2002). As interplays in many-bodied systems result in intricate patterns, the understanding of these requires an in-depth knowledge of the suspension microstructure and statistics. Representative and reliable statistics require a large number Np of interactive swimmers that many simulation methods can hardly afford [e.g. Np ¼ 40 in Mehandia and Nott (2008) or Np # 216 in Ishikawa et al. (2008)]. In the following, a spherical swimmer model is derived from the classical low Reynolds number framework and implemented in the force-coupling method (FCM) for large populations. Resulting statistics reveals non-trivial spatial arrangements of swimmers depending on their swimming gait.

Leveraging Collective Effects in Externally Driven Colloidal Suspensions: Experiments and Simulations

In this review article, we focus on collective motion in externally driven colloidal suspensions, as well as how these collective effects can be harnessed for use in microfluidic applications. We highlight the leading role of hydrodynamic interactions in the self-assembly, emergent behavior, transport, and mixing properties of colloidal suspensions. A special emphasis is given to recent numerical methods to simulate driven colloidal suspensions at large scales. In combination with experiments, they help us to understand emergent dynamics and to identify control parameters for both individual and collective motion in colloidal suspensions.