Source author record

Julia M. Yeomans

Julia M. Yeomans appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

46works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

46 published item(s)

preprint2026arXiv

Bridging Elastic and Active Turbulence

Remarkably, even under negligible inertia, the addition of microstructural agents can generate chaotic flow fields. Such behavior can arise in polymer solutions, leading to elastic turbulence, or from active, self-driven particles, which generate active turbulence. Here, we demonstrate a close and hitherto unrecognized connection between these two classes of turbulence. Specifically, we reveal that their continuum descriptions are analogous at the macroscopic level, such that polymeric fluids can be interpreted as a deformable analogue of contractile active matter. Moreover, our numerical results for Kolmogorov flow demonstrate that the transition into the well-known traveling arrowhead structures in elastic turbulence is marked by the emergence of $\pm 1/2$ topological defects, long recognized as a defining feature of active turbulence, in the polymer director field. Importantly, these coherent structures originate from a transverse instability driven by activity-like gradients generated by anisotropically stretched, contractile polymers. At sufficiently strong activity, the system undergoes a transition into a flow-suppressed state characterized by weak polymer stretching and ordering, a behavior that can be explained by analogy with the spontaneous-flow transition observed in channel-confined active nematics.

preprint2022arXiv

Geometrical control of interface patterning underlies active matter invasion

Interaction between active materials and the boundaries of geometrical confinement is key to many emergent phenomena in active systems. For living active matter consisting of animal cells or motile bacteria, the confinement boundary is often a deformable interface, and it has been unclear how activity-induced interface dynamics might lead to morphogenesis and pattern formation. Here we studied the evolution of bacterial active matter confined by a deformable boundary. We discovered that an ordered morphological pattern emerged at the interface characterized by periodically-spaced interfacial protrusions; behind the interfacial protrusions, bacterial swimmers self-organized into multicellular clusters displaying +1/2 nematic defects. Subsequently, a hierarchical sequence of transitions from interfacial protrusions to creeping branches allowed the bacterial active drop to rapidly invade surrounding space with a striking self-similar branch pattern. We found that this interface patterning is controlled by the local curvature of the interface, a phenomenon we denote as collective curvature sensing. Using a continuum active model, we revealed that the collective curvature sensing arises from enhanced active stresses near high-curvature regions, with the active length-scale setting the characteristic distance between the interfacial protrusions. Our findings reveal a protrusion-to-branch transition as a novel mode of active matter invasion and suggest a new strategy to engineer pattern formation of active materials.

preprint2022arXiv

Self-sustained oscillations of active viscoelastic matter

Models of active nematics in biological systems normally require complexity arising from the hydrodynamics involved at the microscopic level as well as the viscoelastic nature of the system. Here we show that a minimal, space-independent, model based on the temporal alignment of active and polymeric particles provides an avenue to predict and study their coupled dynamics within the framework of dynamical systems. In particular, we examine, using analytical and numerical methods, how such a simple model can display self-sustained oscillations in an activity-driven viscoelastic shear flow.

preprint2022arXiv

Steering self-organisation through confinement

Self-organisation is the spontaneous emergence of spatio-temporal structures and patterns from the interaction of smaller individual units. Examples are found across many scales in very different systems and scientific disciplines, from physics, materials science and robotics to biology, geophysics and astronomy. Recent research has highlighted how self-organisation can be both mediated and controlled by confinement. Confinement occurs through interactions with boundaries, and can function as either a catalyst or inhibitor of self-organisation. It can then become a means to actively steer the emergence or suppression of collective phenomena in space and time. Here, to provide a common framework for future research, we examine the role of confinement in self-organisation and identify overarching scientific challenges across disciplines that need to be addressed to harness its full scientific and technological potential. This framework will not only accelerate the generation of a common deeper understanding of self-organisation but also trigger the development of innovative strategies to steer it through confinement, with impact, e.g., on the design of smarter materials, tissue engineering for biomedicine and crowd management.

preprint2020arXiv

Active Inter-cellular Forces in Collective Cell Motility

The collective behaviour of confluent cell sheets is strongly influenced both by polar forces, arising through cytoskeletal propulsion and by active inter-cellular forces, which are mediated by interactions across cell-cell junctions. We use a phase-field model to explore the interplay between these two contributions and compare the dynamics of a cell sheet when the polarity of the cells aligns to (i) their main axis of elongation, (ii) their velocity, and (iii) when the polarity direction executes a persistent random walk.In all three cases, we observe a sharp transition from a jammed state (where cell rearrangements are strongly suppressed) to a liquid state (where the cells can move freely relative to each other) when either the polar or the inter-cellular forces are increased. In addition, for case (ii) only, we observe an additional dynamical state, flocking (solid or liquid), where the majority of the cells move in the same direction. The flocking state is seen for strong polar forces, but is destroyed as the strength of the inter-cellular activity is increased.

preprint2020arXiv

Active nematics with anisotropic friction: the decisive role of the flow aligning parameter

We use continuum simulations to study the impact of anisotropic hydrodynamic friction on the emergent flows of active nematics. We show that, depending on whether the active particles align with or tumble in their collectively self-induced flows, anisotropic friction can result in markedly different patterns of motion. In a flow-aligning regime and at high anisotropic friction, the otherwise chaotic flows are streamlined into flow lanes with alternating directions, reproducing the experimental laning state that has been obtained by interfacing microtubule-motor protein mixtures with smectic liquid crystals. Within a flow-tumbling regime, however, we find that no such laning state is possible. Instead, the synergistic effects of friction anisotropy and flow tumbling can lead to the emergence of bound pairs of topological defects that align at an angle to the easy flow direction and navigate together throughout the domain. In addition to confirming the mechanism behind the laning states observed in experiments, our findings emphasise the role of the flow aligning parameter in the dynamics of active nematics.

preprint2020arXiv

Degenerate states, emergent dynamics and fluid mixing by magnetic rotors

We investigate the collective motion of magnetic rotors suspended in a viscous fluid under an uniform rotating magnetic field. The rotors are positioned on a square lattice, and low Reynolds hydrodynamics is assumed. For a $3 \times 3$ array of magnets, we observe three characteristic dynamical patterns as the external field strength is varied: a synchronized pattern, an oscillating pattern, and a chessboard pattern. The relative stability of these depends on the competition between the energy due to the external magnetic field and the energy of the magnetic dipole-dipole interactions among the rotors. We argue that the chessboard pattern can be understood as an alternation in the stability of two degenerate states, characterized by striped and spin-ice configurations, as the applied magnetic field rotates. For larger arrays, we observe propagation of slip waves that are similar to metachronal waves. The rotor arrays have potential as microfluidic devices that can mix fluids and create vortices of different sizes.

preprint2020arXiv

Polar jets of swimming bacteria condensed by a patterned liquid crystal

Active matter exhibits remarkable collective behavior in which flows, continuously generated by active particles, are intertwined with the orientational order of these particles. The relationship remains poorly understood as the activity and order are difficult to control independently. Here we demonstrate important facets of this interplay by exploring dynamics of swimming bacteria in a liquid crystalline environment with pre-designed periodic splay and bend in molecular orientation. The bacteria are expelled from the bend regions and condense into polar jets that propagate and transport cargo unidirectionally along the splay regions. The bacterial jets remain stable even when the local concentration exceeds the threshold of bending instability in a non-patterned system. Collective polar propulsion and different role of bend and splay are explained by an advection-diffusion model and by numerical simulations that treat the system as a two-phase active nematic. The ability of prepatterned liquid crystalline medium to streamline the chaotic movements of swimming bacteria into polar jets that can carry cargo along a predesigned trajectory opens the door for potential applications in cell sorting, microscale delivery and soft microrobotics.

preprint2019arXiv

Activity induced nematic order in isotropic liquid crystals

We use linear stability analysis to show that an isotropic phase of elongated particles with dipolar flow fields can develop nematic order as a result of their activity. We argue that ordering is favoured if the particles are flow-aligning and is strongest if the wavevector of the order perturbation is neither parallel nor perpendicular to the nematic director. Numerical solutions of the hydrodynamic equations of motion of an active nematic confirm the results. The instability is contrasted to the well-known hydrodynamic instability of an ordered active nematic.

preprint2019arXiv

Mesoscale modelling of polymer aggregate digestion

We use mesoscale simulations to gain insight into the digestion of biopolymers by studying the break-up dynamics of polymer aggregates (boluses) bound by physical cross-links. We investigate aggregate evolution, establishing that the linking bead fraction and the interaction energy are the main parameters controlling stability with respect to diffusion. We show $\textit{via}$ a simplified model that chemical breakdown of the constituent molecules causes aggregates that would otherwise be stable to disperse. We further investigate breakdown of biopolymer aggregates in the presence of fluid flow. Shear flow in the absence of chemical breakdown induces three different regimes depending on the flow Weissenberg number ($Wi$). i) At $Wi \ll 1$, shear flow has a negligible effect on the aggregates. ii) At $Wi \sim 1$, the aggregates behave approximately as solid bodies and move and rotate with the flow. iii) At $Wi \gg 1$, the energy input due to shear overcomes the attractive cross-linking interactions and the boluses are broken up. Finally, we study bolus evolution under the combined action of shear flow and chemical breakdown, demonstrating a synergistic effect between the two at high reaction rates.

preprint2016arXiv

Active micromachines: Microfluidics powered by mesoscale turbulence

Dense active matter, from bacterial suspensions and microtubule bundles driven by motor proteins to cellular monolayers and synthetic Janus particles, is characterised by mesoscale turbulence, the emergence of chaotic flow structures. By immersing an ordered array of symmetric rotors in an active fluid, we introduce a microfluidic system that exploits spontaneous symmetry breaking in mesoscale turbulence to generate work. The lattice of rotors self-organises into a spin-state where neighbouring discs continuously rotate in permanent alternating directions due to combined hydrodynamic and elastic effects. Our virtual prototype demonstrates a new research direction for the design of micromachines powered by the nematohydrodynamic properties of active turbulence.

preprint2016arXiv

Defect-mediated morphologies in growing cell colonies

Morphological trends in growing colonies of living cells are at the core of physiological and evolutionary processes. Using active gel equations, which include cell division, we show that shape changes during the growth can be regulated by the dynamics of topological defects in the orientation of cells. The friction between the dividing cells and underlying substrate drives anisotropic colony shapes toward more isotropic morphologies, by mediating the number density and velocity of topological defects. We show that the defects interact with the interface at a specific interaction range, set by the vorticity length scale of flows within the colony, and that the cells predominantly reorient parallel to the interface due to division-induced active stresses.

preprint2016arXiv

Hydrodynamics of Micro-swimmers in Films

One of the principal mechanisms by which surfaces and interfaces affect microbial life is by perturbing the hydrodynamic flows generated by swimming. By summing a recursive series of image systems we derive a numerically tractable approximation to the three-dimensional flow fields of a Stokeslet (point force) within a viscous film between a parallel no-slip surface and no-shear interface and, from this Green's function, we compute the flows produced by a force- and torque-free micro-swimmer. We also extend the exact solution of Liron & Mochon (1976) to the film geometry, which demonstrates that the image series gives a satisfactory approximation to the swimmer flow fields if the film is sufficiently thick compared to the swimmer size, and we derive the swimmer flows in the thin-film limit. Concentrating on the thick film case, we find that the dipole moment induces a bias towards swimmer accumulation at the no-slip wall rather than the water-air interface, but that higher-order multipole moments can oppose this. Based on the analytic predictions we propose an experimental method to find the multipole coefficient that induces circular swimming trajectories, allowing one to analytically determine the swimmer's three-dimensional position under a microscope.

preprint2016arXiv

Pore emptying transition during nucleation in hydrophobic nanopores

Using the 2D Ising model we study the generic properties of nucleation in hydrophobic nanopores. To explore the pathways to nucleation of a spin-up phase from a metastable spin-down phase we perform umbrella sampling and transition path sampling simulations. We find that for narrow pores the nucleation occurs on the surface outside the pore. For wide pores the nucleation starts in the pore, and continues outside the filled pore. Intriguingly, we observe a pore emptying transition for a range of intermediate pore widths: a pre-critical nucleus fills the pore, continues to expand outside of the filled pore, but then suddenly gets expelled from the pore before reaching its critical size.

preprint2016arXiv

Stabilization of active matter by flow-vortex lattices and defect ordering

Active systems, from bacterial suspensions to cellular monolayers, are continuously driven out of equilibrium by local injection of energy from their constituent elements and exhibit turbulent-like and chaotic patterns. Here we demonstrate both theoretically and through numerical simulations, that the crossover between wet active systems, whose behaviour is dominated by hydrodynamics, and dry active matter where any flow is screened, can be achieved by using friction as a control parameter. Moreover we discover unexpected vortex ordering at this wet-dry crossover. We show that the self organisation of vortices into lattices is accompanied by the spatial ordering of topological defects leading to active crystal-like structures. The emergence of vortex lattices, which leads to the positional ordering of topological defects, suggests potential applications in the design and control of active materials.

preprint2016arXiv

The macroscopic pancake bounce

We demonstrate that the so-called pancake bounce of millimetric water droplets on surfaces patterned with hydrophobic posts [Nat. Phys. 10, 515 (2014)] can be reproduced on larger scales. In our experiment, a bed of nails plays the role of the structured surface and a water balloon models the water droplet. The macroscopic version largely reproduces the features of the microscopic experiment, including the Weber number dependence and the reduced contact time for pancake bouncing. The scalability of the experiment confirms the mechanisms of pancake bouncing, and allows us to measure the force exerted on the surface during the bounce. The experiment is simple and inexpensive and is an example where front-line research is accessible to student projects.

preprint2015arXiv

Cell division: a source of active stress in cellular monolayers

We introduce the notion of cell division-induced activity and show that the cell division generates extensile forces and drives dynamical patterns in cell assemblies. Extending the hydrodynamic models of lyotropic active nematics we describe turbulent-like velocity fields that are generated by the cell division in a confluent monolayer of cells. We show that the experimentally measured flow field of dividing Madin-Darby Canine Kidney (MDCK) cells is reproduced by our modeling approach. Division-induced activity acts together with intrinsic activity of the cells in extensile and contractile cell assemblies to change the flow and director patterns and the density of topological defects. Finally we model the evolution of the boundary of a cellular colony and compare the fingering instabilities induced by cell division to experimental observations on the expansion of MDCK cell cultures.

preprint2015arXiv

Driven active and passive nematics

We investigate similarities in the micro-structural dynamics between externally driven and actively driven nematics. Walls, lines of strong deformations in the director field, and topological defects are characteristic features of an active nematic. Similar structures form in driven passive nematics when there are inhomogeneities in imposed velocity gradients due to non-linear flow fields or geometrical constraints. Specifically, pressure driven flow of a tumbling passive nematic in an expanding-contracting channel produces walls and defects similar to those seen in active nematics. We also study the response of active nematics to external driving, confirming that imposed shear suppresses the hydrodynamic instabilities. We show that shear fields can lead to wall alignments and the localisation of active turbulence.

preprint2015arXiv

Intrinsic Free Energy in Active Nematics

Basing our arguments on the theory of active liquid crystals, we demonstrate both analytically and numerically, that the activity can induce an effective free energy which enhances ordering in extensile systems of active rods and in contractile suspensions of active discs. We argue that this occurs because any ordering fluctuation is enhanced by the flow field it produces. A phase-diagram in the temperature-activity plane compares ordering due to a thermodynamic free energy to that resulting from the activity. We also demonstrate that activity can drive variations in concentration, but for a different physical reason that relies on the separation of hydrodynamic and diffusive timescales.

preprint2015arXiv

Multi-Particle Collision Dynamics Algorithm for Nematic Fluids

Research on transport, self-assembly and defect dynamics within confined, flowing liquid crystals requires versatile and computationally efficient mesoscopic algorithms to account for fluctuating nematohydrodynamic interactions. We present a multi-particle collision dynamics (MPCD) based algorithm to simulate liquid-crystal hydrodynamic and director fields in two and three dimensions. The nematic-MPCD method is shown to successfully reproduce the features of a nematic liquid crystal, including a nematic-isotropic phase transition with hysteresis in 3D, defect dynamics, isotropic Frank elastic coefficients, tumbling and shear alignment regimes and boundary condition dependent order parameter fields.

preprint2015arXiv

Thermal analogue of gimbal lock in a colloidal ferromagnetic Janus rod

We report an entropy-driven orientational hopping transition in a magnetically confined colloidal Janus rod. In a magnetic field, the sedimented rod randomly hops between horizontal and vertical states: the latter state comes at a substantial gravitational cost at no reduction of magnetic potential energy. The probability distribution over the angles of the rod shows that the presence of an external magnetic field leads to the emergence of a metastable vertical state separated from the ground state by an effective barrier. This barrier does not come from the potential energy but rather from the vast gain in phase space available to the rod as it approaches the vertical state. The loss of rotational degree of freedom that gives rise to this effect is a statistical mechanical analogue of the phenomenon of gimbal lock from classical mechanics.

preprint2015arXiv

Upstream swimming in microbiological flows

Interactions between microorganisms and their complex flowing environments are essential in many biological systems. We develop a model for microswimmer dynamics in non-Newtonian Poiseuille flows. We predict that swimmers in shear-thickening (-thinning) fluids migrate upstream more (less) quickly than in Newtonian fluids and demonstrate that viscoelastic normal stress differences reorient swimmers causing them to migrate upstream at the centreline, in contrast to well-known boundary accumulation in quiescent Newtonian fluids. Based on these observations, we suggest a sorting mechanism to select microbes by swimming speed.

preprint2014arXiv

Active nematic materials with substrate friction

Active turbulence in dense active systems is characterized by high vorticity on a length scale that is large compared to that of individual entities. We describe the properties of active turbulence as momentum propagation is screened by frictional damping. As friction is increased, the spacing between the walls in the nematic director field decreases as a consequence of the more rapid velocity decays. This leads to, first, a regime with more walls and an increased number of topological defects, and then to a jammed state in which the walls deliminate bands of opposing flow, analogous to the shear bands observed in passive complex fluids.

preprint2014arXiv

Biphasic, Lyotropic, Active Nematics

We perform dynamical simulations of a two-dimensional active nematic fluid in coexistence with an isotropic fluid. Drops of active nematic become elongated, and an effective anchoring develops at the nematic-isotropic interface. The activity also causes an undulatory instability of the interface. This results in defects of positive topological charge being ejected into the nematic, leaving the interface with a diffuse negative charge. Quenching the active lyotropic fluid results in a steady state in which phase-separating domains are elongated and then torn apart by active stirring.

preprint2014arXiv

Stirring by swimmers in confined microenvironments

We consider the tracer diffusion $D_{rr}$ that arises from the run-and-tumble motion of low Reynolds number swimmers, such as bacteria. In unbounded dilute suspensions, where the dipole swimmers move in uncorrelated runs of length $λ$, an exact solution showed that $D_{rr}$ is independent of $λ$. Here we verify this result in numerical simulations for a particular model swimmer, the spherical squirmer. We also note that in confined microenvironments, such as microscopic droplets, microfluidic devices and bacterial microzones in marine ecosystems, the size of the system can be comparable to $λ$. We show that this effect alone reduces the value of $D_{rr}$ in comparison to its bulk value, and predict a scaling form for its relative decrease.

preprint2014arXiv

Vorticity, Defects and Correlations in Active Turbulence

We describe a numerical investigation of a continuum model of an active nematic, concentrating on the regime of active turbulence. Results are presented for the effect of three parameters, activity, elastic constant and rotational diffusion constant, on the order parameter and flow fields. Defects and distortions in the director field act as sources of vorticity, and thus vorticity is strongly correlated to the director field. In particular the characteristic length of decay of vorticity and order parameter correlations is controlled by the defect density. By contrast the decay of velocity correlations is determined by a balance between activity and dissipation. We highlight the role of microscopic flow generation mechanisms in determining the flow patterns and characteristic scales of active turbulence and contrast the behaviour of extensile and contractile active nematics.

preprint2013arXiv

Enhanced motility of a microswimmer in rigid and elastic confinement

We analyse the effect of confining rigid and elastic boundaries on the motility of a model dipolar microswimmer. Flexible boundaries are deformed by the velocity field of the swimmer in such a way that the motility of both extensile and contractile swimmers is enhanced. The magnitude of the increase in swimming velocity is controlled by the ratio of the swimmer-advection and elastic timescales, and the dipole moment of the swimmer. We explain our results by considering swimming between inclined rigid boundaries.

preprint2013arXiv

Fluid mixing by curved trajectories of microswimmers

We consider the tracer diffusion $D_{rr}$ that arises from the run-and-tumble motion of low Reynolds number swimmers, such as bacteria. Assuming a dilute suspension, where the bacteria move in uncorrelated runs of length $λ$, we obtain an exact expression for $D_{rr}$ for dipolar swimmers in three dimensions, hence explaining the surprising result that this is independent of $λ$. We compare $D_{rr}$ to the contribution to tracer diffusion from entrainment.

preprint2013arXiv

Instabilities and Topological Defects in Active Nematics

We study a continuum model of an extensile active nematic to show that mesoscale turbulence develops in two stages: (i) ordered regions undergo an intrinsic hydrodynamic instability generating walls, lines of stong bend deformations, (ii) the walls relax by forming oppositely charged pairs of defects. Both creation and annihilation of defect pairs reinstate nematic regions which undergo further instabilities, leading to a dynamic steady state. We compare this with the development of active turbulence in a contractile active nematic.

preprint2013arXiv

Modelling unidirectional liquid spreading on slanted microposts

A lattice Boltzmann algorithm is used to simulate the slow spreading of drops on a surface patterned with slanted micro-posts. Gibb's pinning of the interface on the sides or top of the posts leads to unidirectional spreading over a wide range of contact angles and inclination angles of the posts. Regimes for spreading in no, one or two directions are identified, and shown to agree well with a two-dimensional theory proposed in Chu, Xiao and Wang (Nature Materials, 9, 413). A more detailed numerical analysis of the contact line shapes allows us to understand deviations from the two dimensional model, and to identify the shapes of the pinned interfaces.

preprint2013arXiv

Velocity Correlations in an Active Nematic

The flow properties of a continuum model for an active nematic is studied and compared with recent experiments on suspensions of microtubule bundles and molecular motors. The velocity correlation length is found to be independent of the strength of the activity while the characteristic velocity scale increases monotonically as the activity is increased, both in agreement with the experimental observations. We interpret our results in terms of the creation and annihilation dynamics of a gas of topological defects.

preprint2013arXiv

Viscous fingering at ultralow interfacial tension

We experimentally study the viscous fingering instability in a fluid-fluid phase separated colloid-polymer mixture by means of laser scanning confocal microscopy and microfluidics. We focus on three aspects of the instability. (i) The interface between the two demixed phases has an ultralow surface tension, such that we can address the role of thermal interface fluctuations. (ii) We image the interface in three dimensions allowing us to study the interplay between interface curvature and flow. (iii) The displacing fluid wets all walls completely, in contrast to traditional viscous fingering experiments, in which the displaced fluid wets the walls. We also perform lattice Boltzmann simulations, which help to interpret the experimental observations.

preprint2012arXiv

A circle swimmer at low Reynolds number

Swimming in circles occurs in a variety of situations at low Reynolds number. Here we propose a simple model for a swimmer that undergoes circular motion, generalising the model of a linear swimmer proposed by Najafi and Golestanian (Phys. Rev. E 69, 062901 (2004)). Our model consists of three solid spheres arranged in a triangular configuration, joined by two links of time-dependent length. For small strokes, we discuss the motion of the swimmer as a function of the separation angle between its links. We find that swimmers describe either clockwise or anticlockwise circular motion depending on the tilting angle in a non-trivial manner. The symmetry of the swimmer leads to a quadrupolar decay of the far flow field. We discuss the potential extensions and experimental realisation of our model.

preprint2012arXiv

Anisotropic wetting and de-wetting of drops on substrates patterned with polygonal posts

We present results showing how water drops, produced by ink-jet printing, spread on surfaces patterned with lattices of diamond or triangular posts. Considering post widths typically ~7 m and lattice spacings between 15-40 m, we observe drop shapes with 3,4 and 6-fold symmetry, depending on both the symmetry of the lattice and the shape of the posts. This is a result of the different mechanisms of interface pinning and depinning which depend on the direction of the contact line motion with respect to the post shape. Lattice Boltzmann simulations are used to describe these mechanisms in detail for triangular posts. We also follow the motion of the contact line as the drops evaporate showing that they tend to return to their original shape. To explain this we show that the easy direction for movement is the same for spreading and drying drops. We compare the behaviour of small drops with that of larger drops created by jetting several drops at the same position. We find that the contact line motion is unexpectedly insensitive to drop volume, even when a spherical cap of fluid forms above the posts. The findings are relevant to microfluidic applications and to the control of drop shapes in ink-jet printing.

preprint2012arXiv

Confined active nematic flow in cylindrical capillaries

We use numerical modelling to study the flow patterns of an active nematic confined in a cylindrical capillary, considering both planar and homeotropic boundary conditions. We find that active flow emerges not only along the capillary axis but also within the plane of the capillary, where radial vortices are formed. If topological defects are imposed by the boundary conditions, they act as local pumps driving the flow. At higher activity we demonstrate escape of the active defects and flow into the third dimension, indicating the importance of dimensionality in active materials. We argue that measuring the magnitude of the active flow as a function of the capillary radius allows determination of a value for the activity coefficient.

preprint2012arXiv

Easier sieving through narrower pores: fluctuations and barrier crossing in flow-driven polymer translocation

We show that the injection of polymer chains into nanochannels becomes easier as the channel becomes narrower. This counter intuitive result arises because of a decrease in the diffusive time scale of the chains with increasing confinement. The results are obtained by extending the de Gennes blob model of confined polymers, and confirmed by hybrid molecular dynamics - lattice-Boltzmann simulations.

preprint2012arXiv

Fluid transport by individual microswimmers

We discuss the path of a tracer particle as a microswimmer moves past on an infinite straight trajectory. If the tracer is sufficiently far from the path of the swimmer it moves in a closed loop. As the initial distance between the tracer and the path of the swimmer $ρ$ decreases, the tracer is displaced a small distance backwards (relative to the direction of the swimmer velocity). For much smaller tracer-swimmer separations, however, the tracer displacement becomes positive and diverges as $ρ\to 0$. To quantify this behaviour we calculate the Darwin drift, the total volume swept out by a material sheet of tracers, initially perpendicular to the swimmer path, during the swimmer motion. We find that the drift can be written as the sum of a {\em universal} term which depends on the quadrupolar flow field of the swimmer, together with a non-universal contribution given by the sum of the volumes of the swimmer and its wake. The formula is compared to exact results for the squirmer model and to numerical calculations for a more realistic model swimmer.

preprint2012arXiv

Meso-scale turbulence in living fluids

Turbulence is ubiquitous, from oceanic currents to small-scale biological and quantum systems. Self-sustained turbulent motion in microbial suspensions presents an intriguing example of collective dynamical behavior amongst the simplest forms of life, and is important for fluid mixing and molecular transport on the microscale. The mathematical characterization of turbulence phenomena in active non-equilibrium fluids proves even more difficult than for conventional liquids or gases. It is not known which features of turbulent phases in living matter are universal or system-specific, or which generalizations of the Navier-Stokes equations are able to describe them adequately. Here, we combine experiments, particle simulations, and continuum theory to identify the statistical properties of self-sustained meso-scale turbulence in active systems. To study how dimensionality and boundary conditions affect collective bacterial dynamics, we measured energy spectra and structure functions in dense Bacillus subtilis suspensions in quasi-2D and 3D geometries. Our experimental results for the bacterial flow statistics agree well with predictions from a minimal model for self-propelled rods, suggesting that at high concentrations the collective motion of the bacteria is dominated by short-range interactions. To provide a basis for future theoretical studies, we propose a minimal continuum model for incompressible bacterial flow. A detailed numerical analysis of the 2D case shows that this theory can reproduce many of the experimentally observed features of self-sustained active turbulence.

preprint2012arXiv

Surface evolver simulations of drops on microposts

An important feature in the design of superhydrophobic surfaces is their robustness against collapse from the Cassie-Baxter configuration to the Wenzel state. Upon such a transition a surface loses its properties of low adhesion and friction. We describe how to adapt the Surface Evolver algorithm to predict the parameters and mechanism of the collapse transition on posts of arbitrary shape. In particular, contributions to the free energy evaluated over the solid-liquid surface are reduced to line integrals to give good convergence. The algorithm is validated for straight, vertical and inclined, posts. Numerical results for curved posts with a horizontal section at their ends show that these are more efficient in stabilising the Cassie state than straight posts, and identify whether the interface first depins from the post sides or the post tips.

preprint2011arXiv

Length-dependent translocation of polymers through nanochannels

We consider the flow-driven translocation of single polymer chains through nanochannels. Using analytical calculations based on the de Gennes blob model and mesoscopic numerical simulations, we estimate the threshold flux for the translocation of chains of different number of monomers. The translocation of the chains is controlled by the competition between entropic and hydrodynamic effects, which set a critical penetration length for the chain before it can translocate through the channel. We demonstrate that the polymers show two different translocation regimes depending on how their length under confinement compares to the critical penetration length. For polymer chains longer than the threshold, the translocation process is insensitive to the number of monomers in the chain as predicted in Sakaue {\it et al.}, {\it Euro. Phys. Lett.}, {\bf 72} 83 (2005). However, for chains shorter than the critical length we show that the translocation process is strongly dependent on the length of the chain. We discuss the possible relevance of our results to biological transport.

preprint2010arXiv

CUDA simulations of active dumbbell suspensions

We describe and analyze CUDA simulations of hydrodynamic interactions in active dumbbell suspensions. GPU-based parallel computing enables us not only to study the time-resolved collective dynamics of up to a several hundred active dumbbell swimmers but also to test the accuracy of effective time-averaged models. Our numerical results suggest that the stroke-averaged model yields a relatively accurate description down to distances of only a few times the dumbbell's length. This is remarkable in view of the fact that the stroke-averaged model is based on a far-field expansion. Thus, our analysis confirms that stroke-averaged far-field equations of motion may provide a useful starting point for the derivation of hydrodynamic field equations.

preprint2010arXiv

Lévy Fluctuations and Tracer Diffusion in Dilute Suspensions of Algae and Bacteria

Swimming microorganisms rely on effective mixing strategies to achieve efficient nutrient influx. Recent experiments, probing the mixing capability of unicellular biflagellates, revealed that passive tracer particles exhibit anomalous non-Gaussian diffusion when immersed in a dilute suspension of self-motile Chlamydomonas reinhardtii algae. Qualitatively, this observation can be explained by the fact that the algae induce a fluid flow that may occasionally accelerate the colloidal tracers to relatively large velocities. A satisfactory quantitative theory of enhanced mixing in dilute active suspensions, however, is lacking at present. In particular, it is unclear how non-Gaussian signatures in the tracers' position distribution are linked to the self-propulsion mechanism of a microorganism. Here, we develop a systematic theoretical description of anomalous tracer diffusion in active suspensions, based on a simplified tracer-swimmer interaction model that captures the typical distance scaling of a microswimmer's flow field. We show that the experimentally observed non-Gaussian tails are generic and arise due to a combination of truncated Lévy statistics for the velocity field and algebraically decaying time correlations in the fluid. Our analytical considerations are illustrated through extensive simulations, implemented on graphics processing units to achieve the large sample sizes required for analyzing the tails of the tracer distributions.

preprint2010arXiv

Swimmer-tracer scattering at low Reynolds number

Understanding the stochastic dynamics of tracer particles in active fluids is important for identifying the physical properties of flow generating objects such as colloids, bacteria or algae. Here, we study both analytically and numerically the scattering of a tracer particle in different types of time-dependent, hydrodynamic flow fields. Specifically, we compare the tracer motion induced by an externally driven colloid with the one generated by various self-motile, multi-sphere swimmers. Our results suggest that force-free swimmers generically induce loop-shaped tracer trajectories. The specific topological structure of these loops is determined by the hydrodynamic properties of the microswimmer. Quantitative estimates for typical experimental conditions imply that the loops survive on average even if Brownian motion effects are taken into account.

preprint2009arXiv

Hydrodynamic Synchronisation of Model Microswimmers

We define a model microswimmer with a variable cycle time, thus allowing the possibility of phase locking driven by hydrodynamic interactions between swimmers. We find that, for extensile or contractile swimmers, phase locking does occur, with the relative phase of the two swimmers being, in general, close to 0 or pi, depending on their relative position and orientation. We show that, as expected on grounds of symmetry, self T-dual swimmers, which are time-reversal covariant, do not phase-lock. We also discuss the phase behaviour of a line of tethered swimmers, or pumps. These show oscillations in their relative phases reminiscent of the metachronal waves of cilia.