Source author record

Eric Lauga

Eric Lauga 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

90works
10topics
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

90 published item(s)

preprint2022arXiv

Elastohydrodynamic synchronization of rotating bacterial flagella

To rotate continuously without jamming, the flagellar filaments of bacteria need to be locked in phase. While several models have been proposed for eukaryotic flagella, the synchronization of bacterial flagella is less well understood. Starting from a reduced model of flexible and hydrodynamically-coupled bacterial flagella, we rigorously coarse-grain the equations of motion using the method of multiple scales, and hence show that bacterial flagella generically synchronize to zero phase difference via an elastohydrodynamic mechanism. Remarkably, the far-field rate of synchronization is maximized at an intermediate value of elastic compliance, with surprising implications for bacteria.

preprint2022arXiv

Jet-driven viscous locomotion of confined thermoresponsive microgels

We consider the dynamics of micro-sized, asymmetrically-coated thermoresponsive hydrogel ribbons (microgels) under periodic heating and cooling in the confined space between two planar surfaces. As the result of the temperature changes, the volume and thus the shape of the slender microgel change, which lead to repeated cycles of bending and elastic relaxation, and to net locomotion. Small devices designed for biomimetic locomotion need to exploit flows that are not symmetric in time (non-reciprocal) to escape the constraints of the scallop theorem and undergo net motion. Unlike other biological slender swimmers, the non-reciprocal bending of the gel centreline is not sufficient here to explain for the overall swimming motion. We show instead that the swimming of the gel results from the flux of water periodically emanating from (or entering) the gel itself due to its shrinking (or swelling). The associated flows induce viscous stresses that lead to a net propulsive force on the gel. We derive a theoretical model for this hypothesis of jet-driven propulsion, which leads to excellent agreement with our experiments.

preprint2022arXiv

Order and Information in the Patterns of Spinning Magnetic Micro-disks at the Air-water Interface

The application of the Shannon entropy to study the relationship between information and structures has yielded insights into molecular and material systems. However, the difficulty in directly observing and manipulating atoms and molecules hampers the ability of these systems to serve as model systems for further exploring the links between information and structures. Here, we use, as a model experimental system, hundreds of spinning magnetic micro-disks self-organizing at the air-water interface to generate various spatiotemporal patterns with varying degrees of orders. Using the neighbor distance as the information-bearing variable, we demonstrate the links among information, structure, and interactions. Most importantly, we establish a direct link between information and structure without using explicit knowledge of interactions. Finally, we show that the Shannon entropy by neighbor distances is a powerful observable in characterizing structural changes. Our findings are relevant for analyzing natural self-organizing systems and for designing collective robots.

preprint2021arXiv

Asymptotic theory of hydrodynamic interactions between slender filaments

Hydrodynamic interactions (HIs) are important in biophysics research because they influence both the collective and the individual behaviour of microorganisms and self-propelled particles. For instance, HIs at the micro-swimmer level determine the attraction or repulsion between individuals, and hence their collective behaviour. Meanwhile, HIs between swimming appendages (e.g. cilia and flagella) influence the emergence of swimming gaits, synchronised bundles and metachronal waves. In this study, we address the issue of HIs between slender filaments separated by a distance larger than their contour length (d>L) by means of asymptotic calculations and numerical simulations. We first derive analytical expressions for the extended resistance matrix of two arbitrarily-shaped rigid filaments as a series expansion in inverse powers of d/L>1. The coefficients in our asymptotic series expansion are then evaluated using two well-established methods for slender filaments, resistive-force theory (RFT) and slender-body theory (SBT), and our asymptotic theory is verified using numerical simulations based on SBT for the case of two parallel helices. The theory captures the qualitative features of the interactions in the regime d/L>1, which opens the path to a deeper physical understanding of hydrodynamically governed phenomena such as inter-filament synchronisation and multiflagellar propulsion. To demonstrate the usefulness of our results, we next apply our theory to the case of two helices rotating side-by-side, where we quantify the dependence of all forces and torques on the distance and phase difference between them. Using our understanding of pairwise HIs, we then provide physical intuition for the case of a circular array of rotating helices. Our theoretical results will be useful for the study of HIs between bacterial flagella, nodal cilia, and slender microswimmers.

preprint2020arXiv

3-D Printed Swimming Microtori for Cargo Transport and Flow Manipulation

Through billions of years of evolution, microorganisms mastered unique swimming behaviors to thrive in complex fluid environments. Limitations in nanofabrication have thus far hindered the ability to design and program synthetic swimmers with the same abilities. Here we encode multi-behavioral responses in artificial swimmers such as microscopic, self-propelled tori using nanoscale 3D printing. We show experimentally and theoretically that the tori continuously transition between two primary swimming modes in response to a magnetic field. The tori also manipulate and transport other artificial swimmers, bimetallic nanorods, as well as passive colloidal particles. In the first behavioral mode, the tori accumulate and transport nanorods; in the second mode, nanorods align along the tori's self-generated streamlines. Our results indicate that such shape-programmed microswimmers have the potential to manipulate biological active matter, e.g. bacteria or cells.

preprint2020arXiv

Collective stiffening of soft hair assemblies

Many living systems use assemblies of soft and slender structures whose deflections allow them to mechanically probe their immediate environment. In this work, we study the collective response of artificial soft hair assemblies to a shear flow by imaging their deflections. At all hair densities, the deflection is found to be proportional to the local shear stress with a proportionality factor that decreases with density. The measured collective stiffening of hairs is modeled both with a microscopic elastohydrodynamic model that takes into account long range hydrodynamic hair-hair interactions and a phenomenological model that treats the hair assemblies as an effective porous medium. While the microscopic model is in reasonable agreement with the experiments at low hair density, the phenomenological model is found to be predictive across the entire density range.

preprint2020arXiv

Geometrical constraints on the tangling of bacterial flagellar filaments

Many species of bacteria swim through viscous environments by rotating multiple helical flagella. The filaments gather behind the cell body and form a close helical bundle, which propels the cell forward during a "run". The filaments inside the bundle cannot be continuously actuated, nor can they easily unbundle, if they are tangled around one another. The fact that bacteria can passively form coherent bundles, i.e. bundles which do not contain tangled pairs of filaments, may appear surprising given that flagella are actuated by uncoordinated motors. In this article, we establish the theoretical conditions under which a pair of rigid helical filaments can form a tangled bundle, and we compare these constraints with experimental data collected from the literature. Our results suggest that bacterial flagella are too straight and too far apart to form tangled bundles based on their intrinsic, undeformed geometry alone. This makes the formation of coherent bundles more robust against the passive nature of the bundling process, where the position of individual filaments cannot be controlled.

preprint2020arXiv

Irreversible hydrodynamic trapping by surface rollers

A colloidal particle driven by externally actuated rotation can self-propel parallel to a rigid boundary by exploiting the hydrodynamic coupling that surfaces induce between translation and rotation. As such a roller moves along the boundary it generates local vortical flows, which can be used to trap and transport passive cargo particles. However, the details and conditions for this trapping mechanism have not yet been fully understood. Here, we show that the trapping of cargo is accomplished through time-irreversible interactions between the cargo and the boundary, leading to its migration across streamlines into a steady flow vortex next to the roller. The trapping mechanism is explained analytically with a two dimensional model, investigated numerically in three dimensions for a wide range of parameters and is shown to be analogous to the deterministic lateral displacement (DLD) technique used in microfluidics for the separation of differently sized particles. The several geometrical parameters of the problem are analysed and we predict that thin, disc-like rollers offer the most favourable trapping conditions.

preprint2020arXiv

Selectively Controlled Magnetic Microrobots with Opposing Helices

Magnetic microrobots that swim through liquid media are of interest for minimally invasive medical procedures, bioengineering, and manufacturing. Many of the envisaged applications, such as micromanipulation and targeted cargo delivery, necessitate the use and adequate control of multiple microrobots, which will increase the velocity, robustness, and efficacy of a procedure. While various methods involving heterogeneous geometries, magnetic properties, and surface chemistries have been proposed to enhance independent control, the main challenge has been that the motion between all microwsimmers remains coupled through the global control signal of the magnetic field. Katsamba and Lauga proposed transchiral microrobots, a theoretical design with magnetized spirals of opposite handedness. The competition between the spirals can be tuned to give an intrinsic nonlinearity that each device can function only within a given band of frequencies. This allows individual microrobots to be selectively controlled by varying the frequency of the rotating magnetic field. Here we present the experimental realization and characterization of transchiral micromotors composed of independently driven magnetic helices. We show a swimming micromotor that yields negligible net motion until a critical frequency is reached and a micromotor that changes its translation direction as a function of the frequency of the rotating magnetic field. This work demonstrates a crucial step towards completely decoupled and addressable swimming magnetic microrobots.

preprint2020arXiv

Stochastic dynamics of dissolving active particles

The design of artificial microswimmers has generated significant research interest in recent years, for promise in applications such as nanomotors and targeted drug-delivery. However, many current designs suffer from a common problem, namely the swimmers remain in the fluid indefinitely, posing risks of clogging and damage. Inspired by recently proposed experimental designs, we investigate mathematically the dynamics of degradable active particles. We develop and compare two distinct chemical models for the decay of a swimmer, taking into account the material composition and nature of the chemical or enzymatic reaction at its surface. These include a model for dissolution without a reaction, as well as models for a reacting swimmer studied in the limit of large and small Damköhler number. A new dimensionless parameter emerges that allows the classification of colloids into ballistic and diffusive type. Using this parameter, we perform an asymptotic analysis to derive expressions for colloid lifetimes and their total mean-squared displacement from release and validate these by numerical Monte Carlo simulations of the associated Langevin dynamics. Supported by general scaling relationships, our theoretical results provide new insight into the experimental applicability of a wide range of designs for degradable active colloids.

preprint2020arXiv

Stokes flow due to point torques and sources in a spherical geometry

Solutions to the Stokes equations written in terms of a small number of hydrodynamic image singularities have been a useful tool in theoretical and numerical computations for nearly fifty years. In this article, we extend the catalogue of known solutions by deriving the flow expressions due to a general point torque and point source in the presence of a stationary sphere with either a no-slip or a stress-free (no shear) boundary condition. For an axisymmetric point torque and a no-slip sphere the image system simplifies to a single image point torque, reminiscent of the solution for a point charge outside an equipotential sphere in electrostatics. By symmetry, this also gives a simple representation of the solution due to an axisymmetric point torque inside a rigid spherical shell. In all remaining cases, the solution can be described by a collection of physically intuitive point and line singularities. Our results will be useful for the theoretical modelling of the propulsion of microswimmers and efficient numerical implementation of far-field hydrodynamic interactions in this geometry.

preprint2020arXiv

Swirling Instability of the Microtubule Cytoskeleton

In the cellular phenomena of cytoplasmic streaming, molecular motors carrying cargo along a network of microtubules entrain the surrounding fluid. The piconewton forces produced by individual motors are sufficient to deform long microtubules, as are the collective fluid flows generated by many moving motors. Studies of streaming during oocyte development in the fruit fly $D.~melanogaster$ have shown a transition from a spatially-disordered cytoskeleton, supporting flows with only short-ranged correlations, to an ordered state with a cell-spanning vortical flow. To test the hypothesis that this transition is driven by fluid-structure interactions we study a discrete-filament model and a coarse-grained continuum theory for motors moving on a deformable cytoskeleton, both of which are shown to exhibit a $swirling~instability$ to spontaneous large-scale rotational motion, as observed.

preprint2016arXiv

Elastohydrodynamic Synchronization of Adjacent Beating Flagella

It is now well established that nearby beating pairs of eukaryotic flagella or cilia typically synchronize in phase. A substantial body of evidence supports the hypothesis that hydrodynamic coupling between the active filaments, combined with waveform compliance, provides a robust mechanism for synchrony. This elastohydrodynamic mechanism has been incorporated into `bead-spring' models in which flagella are represented by microspheres tethered by radial springs as internal forces drive them about orbits. While these low-dimensional models reproduce the phenomenon of synchrony, their parameters are not readily relatable to those of flagella. More realistic models which reflect the elasticity of the axonemes and active force generation take the form of fourth-order nonlinear PDEs. While computational studies have shown synchrony, the effects of hydrodynamic coupling between nearby filaments governed by such models have been theoretically examined only in the regime of interflagellar distances $d$ large compared to flagellar length $L$. Yet, in many biological situations $d/L \ll 1$. Here, we first present an asymptotic analysis of the hydrodynamic coupling between two filaments in the regime $d/L \ll 1$, and find that the form of the coupling is independent of the details of the internal forces that govern the motion of the filaments. The analysis is like the localized induction approximation for vortex filament motion, extended to the case of mutual induction. To understand how the coupling mechanism leads to synchrony of extended objects, we introduce a heuristic model of flagellar beating, a single fourth-order nonlinear PDE whose form is derived from symmetry considerations, the physics of elasticity, and the overdamped nature of the dynamics. Analytical and numerical studies of this model illustrate how synchrony between two filaments is achieved through the asymptotic coupling.

preprint2016arXiv

Flagellar flows around bacterial swarms

Flagellated bacteria on nutrient-rich substrates can differentiate into a swarming state and move in dense swarms across surfaces. A recent experiment measured the flow in the fluid around an Escherichia coli swarm (Wu, Hosu and Berg, 2011 Proc. Natl. Acad. Sci. USA 108 4147). A systematic chiral flow was observed in the clockwise direction (when viewed from above) ahead of the swarm with flow speeds of about $10~μ$m/s, about 3 times greater than the radial velocity at the edge of the swarm. The working hypothesis is that this flow is due to the action of cells stalled at the edge of a colony that extend their flagellar filaments outwards, moving fluid over the virgin agar. In this work we quantitatively test his hypothesis. We first build an analytical model of the flow induced by a single flagellum in a thin film and then use the model, and its extension to multiple flagella, to compare with experimental measurements. The results we obtain are in agreement with the flagellar hypothesis. The model provides further quantitative insight into the flagella orientations and their spatial distributions as well as the tangential speed profile. In particular, the model suggests that flagella are on average pointing radially out of the swarm and are not wrapped tangentially.

preprint2016arXiv

Flow analysis of the low-Reynolds number swimmer C. elegans

Swimming cells and microorganisms are a critical component of many biological processes. In order to better interpret experimental studies of low Reynolds number swimming, we combine experimental and numerical methods to perform an analysis of the flow-field around the swimming nematode Caenorhabditis elegans. We first use image processing and particle tracking velocimetry to extract the body shape, kinematics, and flow-fields around the nematode. We then construct a three-dimensional model using the experimental swimming kinematics and employ a boundary element method to simulate flow-fields, obtaining very good quantitative agreement with experiment. We use this numerical model to show that calculation of flow shear rates using purely planar data results in significant underestimates of the true three-dimensional value. Applying symmetry arguments, validated against numerics, we demonstrate that the out-of-plane contribution can be accounted for via incompressibility and therefore simply calculated from particle tracking velocimetry. Our results show how fundamental fluid mechanics considerations may be used to improve the accuracy of measurements in biofluiddynamics.

preprint2016arXiv

Hydrodynamic interactions between nearby slender filaments

Cellular biology abound with filaments interacting through fluids, from intracellular microtubules, to rotating flagella and beating cilia. While previous work has demonstrated the complexity of capturing nonlocal hydrodynamic interactions between moving filaments, the problem remains difficult theoretically. We show here that when filaments are closer to each other than their relevant length scale, the integration of hydrodynamic interactions can be approximately carried out analytically. This leads to a set of simplified local equations, illustrated on a simple model of two interacting filaments, which can be used to tackle theoretically a range of problems in biology and physics.

preprint2016arXiv

Microscale Flow Dynamics of Ribbons and Sheets

Numerical study of the hydrodynamics of thin sheets and ribbons presents difficulties associated with resolving multiple length scales. To circumvent these difficulties, asymptotic methods have been developed to describe the dynamics of slender fibres and ribbons. However, such theories entail restrictions on the shapes that can be studied, and often break down in regions where standard boundary element methods are still impractical. In this paper we develop a regularised stokeslet method for ribbons and sheets in order to bridge the gap between asymptotic and boundary element methods. The method is validated against the analytical solution for plate ellipsoids, as well as the dynamics of ribbon helices and an experimental microswimmer. We then demonstrate the versatility of this method by calculating the flow around a double helix, and the swimming dynamics of a microscale "magic carpet".

preprint2016arXiv

Phoretic flow induced by asymmetric confinement

Internal phoretic flows due to the interactions of solid boundaries with local chemical gradients may be created using chemical patterning. Alternatively, we demonstrate here that internal flows might also be induced by geometric asymmetries of chemically-homogeneous surfaces. We characterise the circulatory flow created in a cavity enclosed between two eccentric cylindrical walls of uniform chemical activity. Local gradients of the diffusing solute induce a slip flow along the surface of the cylinders, leading to a circulatory bulk flow pattern which can be solved analytically in the diffusive limit. The flow strength can be controlled by adjusting the relative positions of the cylinders and an optimal configuration is identified. These results provide a model system for tunable phoretic pumps.

preprint2016arXiv

Rotation of slender swimmers in isotropic-drag media

The drag anisotropy of slender filaments is a critical physical property allowing swimming in low-Reynolds number flows, and without it linear translation is impossible. Here we show that, in contrast, net rotation can occur under isotropic drag. We first demonstrate this result formally by considering the consequences of the force- and torque-free conditions on swimming bodies and we then illustrate it with two examples (a simple swimmers made of three rods and a model bacterium with two helical flagellar filaments). Our results highlight the different role of hydrodynamic forces in generating translational vs.~rotational propulsion.

preprint2016arXiv

Slender-ribbon theory

Ribbons are long narrow strips possessing three distinct material length scales (thickness, width, and length) which allow them to produce unique shapes unobtainable by wires or filaments. For example when a ribbon has half a twist and is bent into a circle it produces a Möbius strip. Significant effort has gone into determining the structural shapes of ribbons but less is know about their behavior in viscous fluids. In this paper we determine, asymptotically, the leading-order hydrodynamic behavior of a slender ribbon in Stokes flows. The derivation, reminiscent of slender-body theory for filaments, assumes that the length of the ribbon is much larger than its width, which itself is much larger than its thickness. The final result is an integral equation for the force density on a mathematical ruled surface, termed the ribbon plane, located inside the ribbon. A numerical implementation of our derivation shows good agreement with the known hydrodynamics of long flat ellipsoids, and successfully captures the swimming behavior of artificial microscopic swimmers recently explored experimentally. We also study the asymptotic behavior of a ribbon bent into a helix, that of a twisted ellipsoid, and we investigate how accurately the hydrodynamics of a ribbon can be effectively captured by that of a slender filament. Our asymptotic results provide the fundamental framework necessary to predict the behavior of slender ribbons at low Reynolds numbers in a variety of biological and engineering problems.

preprint2015arXiv

A reciprocal theorem for boundary-driven channel flows

In a variety of physical situations, a bulk viscous flow is induced by a distribution of surface velocities, for example in diffusiophoresis (as a result of chemical gradients) and above carpets of cilia (as a result of biological activity). When such boundary-driven flows are used to pump fluids, the primary quantity of interest is the induced flow rate. In this letter we propose a method, based on the reciprocal theorem of Stokes flows, to compute the net flow rate for arbitrary flow distribution and periodic pump geometry using solely stress information from a dual Poiseuille-like problem. After deriving the general result we apply it to straight channels of triangular, elliptic and rectangular geometries and quantify the relationship between bulk motion and surface forcing.

preprint2015arXiv

A regularised singularity approach to phoretic problems

An efficient, accurate, and flexible numerical method is proposed for the solution of the swimming problem of one or more autophoretic particles in the purely-diffusive limit. The method relies on successive boundary element solutions of the Laplacian and the Stokes flow equations using regularised Green's functions for swift, simple implementations, an extension of the well-known method of "regularised stokeslets" for Stokes flow problems. The boundary element method is particularly suitable for phoretic problems, since no quantities in the domain bulk are required to compute the swimming velocity. For time-dependent problems, the method requires no re-meshing and simple boundaries such as a plane wall may be added at no increase to the size of the linear system through the method of images. The method is validated against two classical examples for which an analytical or semi-analytical solution is known, a two-sphere system and a Janus particle, and provides a rigorous computational pipeline to address further problems with complex geometry and multiple bodies.

preprint2015arXiv

Arbitrary axisymmetric steady streaming: Flow, force and propulsion

A well-developed method to induce mixing on microscopic scales is to exploit flows generated by steady streaming. Steady streaming is a classical fluid dynamics phenomenon whereby a time-periodic forcing in the bulk or along a boundary is enhanced by inertia to induce a non-zero net flow. Building on classical work for simple geometrical forcing and motivated by the complex shape oscillations of elastic capsules and bubbles, we develop the mathematical framework to quantify the steady streaming of a spherical body with arbitrary axisymmetric time-periodic boundary conditions. We compute the flow asymptotically for small-amplitude oscillations of the boundary in the limit where the viscous penetration length scale is much smaller than the body. In that case, the flow has a boundary layer structure and the fluid motion is solved by asymptotic matching. Our results, presented in the case of no-slip boundary conditions and extended to include the motion of vibrating free surfaces, recovers classical work as particular cases. We illustrate the flow structure given by our solution and propose one application of our results for small-scale force-generation and synthetic locomotion.

preprint2015arXiv

Bacterial hydrodynamics

Bacteria predate plants and animals by billions of years. Today, they are the world's smallest cells yet they represent the bulk of the world's biomass, and the main reservoir of nutrients for higher organisms. Most bacteria can move on their own, and the majority of motile bacteria are able to swim in viscous fluids using slender helical appendages called flagella. Low-Reynolds-number hydrodynamics is at the heart of the ability of flagella to generate propulsion at the micron scale. In fact, fluid dynamic forces impact many aspects of bacteriology, ranging from the ability of cells to reorient and search their surroundings to their interactions within mechanically and chemically-complex environments. Using hydrodynamics as an organizing framework, we review the biomechanics of bacterial motility and look ahead to future challenges.

preprint2015arXiv

Geometric capture and escape of a microswimmer colliding with an obstacle

Motivated by recent experiments, we consider the hydrodynamic capture of a microswimmer near a stationary spherical obstacle. Simulations of model equations show that a swimmer approaching a small spherical colloid is simply scattered. In contrast, when the colloid is larger than a critical size it acts as a passive trap: the swimmer is hydrodynamically captured along closed trajectories and endlessly orbits around the colloidal sphere. In order to gain physical insight into this hydrodynamic scattering problem, we address it analytically. We provide expressions for the critical trapping radius, the depth of the "basin of attraction," and the scattering angle, which show excellent agreement with our numerical findings. We also demonstrate and rationalize the strong impact of swimming-flow symmetries on the trapping efficiency. Finally, we give the swimmer an opportunity to escape the colloidal traps by considering the effects of Brownian, or active, diffusion. We show that in some cases the trapping time is governed by an Ornstein-Uhlenbeck process, which results in a trapping time distribution that is well-approximated as inverse-Gaussian. The predictions again compare very favorably with the numerical simulations. We envision applications of the theory to bioremediation, microorganism sorting techniques, and the study of bacterial populations in heterogeneous or porous environments.

preprint2015arXiv

Geometric pumping in autophoretic channels

Many microfluidic devices use macroscopic pressure differentials to overcome viscous friction and generate flows in microchannels. In this work, we investigate how the chemical and geometric properties of the channel walls can drive a net flow by exploiting the autophoretic slip flows induced along active walls by local concentration gradients of a solute species. We show that chemical patterning of the wall is not required to generate and control a net flux within the channel, rather channel geometry alone is sufficient. Using numerical simulations, we determine how geometric characteristics of the wall influence channel flow rate, and confirm our results analytically in the asymptotic limit of lubrication theory.

preprint2015arXiv

Non-decaying hydrodynamic interactions along narrow channels

Particle-particle interactions are of paramount importance in every multi-body system as they determine the collective behaviour and coupling strength. Many well-known interactions like electro-static, van der Waals or screened Coulomb, decay exponentially or with negative powers of the particle spacing r. Similarly, hydrodynamic interactions between particles undergoing Brownian motion decay as 1/r in bulk, and are assumed to decay in small channels. Such interactions are ubiquitous in biological and technological systems. Here we confine two particles undergoing Brownian motion in narrow, microfluidic channels and study their coupling through hydrodynamic interactions. Our experiments show that the hydrodynamic particle-particle interactions are distance-independent in these channels. This finding is of fundamental importance for the interpretation of experiments where dense mixtures of particles or molecules diffuse through finite length, water-filled channels or pore networks.

preprint2015arXiv

Phase-separation models for swimming enhancement in complex fluids

Swimming cells often have to self-propel through fluids displaying non-Newtonian rheology. While past theoretical work seems to indicate that stresses arising from complex fluids should systematically hinder low-Reynolds number locomotion, experimental observations suggest that locomotion enhancement is possible. In this paper we propose a physical mechanism for locomotion enhancement of microscopic swimmers in a complex fluid. It is based on the fact that micro-structured fluids will generically phase-separate near surfaces, leading to the presence of low-viscosity layers which promote slip and decrease viscous friction near the surface of the swimmer. We use two models to address the consequence of this phase separation: a nonzero apparent slip length for the fluid and then an explicit modeling of the change of viscosity in a thin layer near the swimmer. Considering two canonical setups for low-Reynolds number locomotion, namely the waving locomotion of a two-dimensional sheet and that of a three-dimensional filament, we show that phase-separation systematically increases the locomotion speeds, possibly by orders of magnitude. We close by confronting our predictions with recent experimental results.

preprint2015arXiv

Small-amplitude swimmers can self-propel faster in viscoelastic fluids

Many small organisms self-propel in viscous fluids using travelling wave-like deformation of their bodies or appendages. Examples include small nematodes moving through soil using whole-body undulations or spermatozoa swimming through mucus using flagellar waves. When self-propulsion occurs in a non-Newtonian fluid, one fundamental question is whether locomotion will occur faster or slower than in a Newtonian environment. Here we consider the general problem of swimming using small-amplitude periodic waves in a viscoelastic fluid described by the classical Oldroyd-B constitutive relationship. Using Taylor's swimming sheet model, we show that if all travelling waves move in the same direction, the locomotion speed of the organism is systematically decreased. However, if we allow waves to travel in two opposite directions, we show that this can lead to enhancement of the swimming speed, which is physically interpreted as due to asymmetric viscoelastic damping of waves with different frequencies. A change of the swimming direction is also possible. By analysing in detail the cases of swimming using two or three travelling waves, we demonstrate that swimming can be enhanced in a viscoelastic fluid for all Deborah numbers below a critical value or, for three waves or more, only for a finite, non-zero range of Deborah numbers, in which case a finite amount of elasticity in the fluid is required to increase the swimming speed.

preprint2015arXiv

Theory of locomotion through complex fluids

Microorganisms such as bacteria often swim in fluid environments that cannot be classified as Newtonian. Many biological fluids contain polymers or other heterogeneities which may yield complex rheology. For a given set of boundary conditions on a moving organism, flows can be substantially different in complex fluids, while non-Newtonian stresses can alter the gait of the microorganisms themselves. Heterogeneities in the fluid may also be characterized by length scales on the order of the organism itself leading to additional dynamic complexity. In this chapter we present a theoretical overview of small-scale locomotion in complex fluids with a focus on recent efforts quantifying the impact of non-Newtonian rheology on swimming microorganisms.

preprint2014arXiv

Asymmetric steady streaming as a mechanism for acoustic propulsion of rigid bodies

Recent experiments showed that standing acoustic waves could be exploited to induce self-propulsion of rigid metallic particles in the direction perpendicular to the acoustic wave. We propose in this paper a physical mechanism for these observations based on the interplay between inertial forces in the fluid and the geometrical asymmetry of the particle shape. We consider an axisymmetric rigid near-sphere oscillating in a quiescent fluid along a direction perpendicular to its symmetry axis. The kinematics of oscillations can be either prescribed or can result dynamically from the presence of an external oscillating velocity field. Steady streaming in the fluid, the inertial rectification of the time-periodic oscillating flow, generates steady stresses on the particle which, in general, do not average to zero, resulting in a finite propulsion speed along the axis of the symmetry of the particle and perpendicular to the oscillation direction. Our derivation of the propulsion speed is obtained at leading order in the Reynolds number and the deviation of the shape from that of a sphere. The results of our model are consistent with the experimental measurements, and more generally explains how time periodic forcing from an acoustic field can be harnessed to generate autonomous motion.

preprint2014arXiv

Dynamics of swimming bacteria at complex interfaces

Flagellated bacteria exploiting helical propulsion are known to swim along circular trajectories near surfaces. Fluid dynamics predicts this circular motion to be clockwise (CW) above a rigid surface (when viewed from inside the fluid) and counter-clockwise (CCW) below a free surface. Recent experimental investigations showed that complex physicochemical processes at the nearby surface could lead to a change in the direction of rotation, both at solid surfaces absorbing slip-inducing polymers and interfaces covered with surfactants. Motivated by these results, we use a far-field hydrodynamic model to predict the kinematics of swimming near three types of interfaces: clean fluid-fluid interface, slipping rigid wall, and a fluid interface covered by incompressible surfactants. Representing the helical swimmer by a superposition of hydrodynamic singularities, we first show that in all cases the surfaces reorient the swimmer parallel to the surface and attract it, both of which are a consequence of the Stokes dipole component of the swimmer flow field. We then show that circular motion is induced by a higher-order singularity, namely a rotlet dipole, and that its rotation direction (CW vs. CCW) is strongly affected by the boundary conditions at the interface and the bacteria shape. Our results suggest thus that the hydrodynamics of complex interfaces provide a mechanism to selectively stir bacteria.

preprint2014arXiv

Enhanced active swimming in viscoelastic fluids

Swimming microorganisms often self propel in fluids with complex rheology. While past theoretical work indicates that fluid viscoelasticity should hinder their locomotion, recent experiments on waving swimmers suggest a possible non-Newtonian enhancement of locomotion. We suggest a physical mechanism, based on fluid-structure interaction, leading to swimming in a viscoelastic fluid at a higher speed than in a Newtonian one. Using Taylor's two-dimensional swimming sheet model, we solve for the shape of an active swimmer as a balance between the external fluid stresses, the internal driving moments, and the passive elastic resistance. We show that this dynamic balance leads to a generic transition from hindered rigid swimming to enhanced flexible locomotion. The results are physically interpreted as due to a viscoelastic suction increasing the swimming amplitude in a non-Newtonian fluid and overcoming viscoelastic damping.

preprint2014arXiv

Generalized squirming motion of a sphere

A number of swimming microorganisms such as ciliates ($\textit{Opalina}$) and multicellular colonies of flagellates ($\textit{Volvox}$) are approximately spherical in shape and swim using beating arrays of cilia or short flagella covering their surfaces. Their physical actuation on the fluid may be mathematically modeled as the generation of surface velocities on a continuous spherical surface -- a model known in the literature as squirming, which has been used to address various aspects of the biological physics of locomotion. Previous analyses of squirming assumed axisymmetric fluid motion and hence restricted all swimming kinematics to take place along a line. In this paper we generalize squirming to three spatial dimensions. We derive analytically the flow field surrounding a spherical squirmer with arbitrary surface motion, and use it to derive its three-dimensional translational and rotational swimming kinematics. We then use our results to physically interpret the flow field induced by the swimmer in terms of fundamental flow singularities up to terms decaying spatially as $\sim 1/r^3$. Our results will enable to develop new models in biological physics, in particular in the area of hydrodynamic interactions and collective locomotion.

preprint2014arXiv

Locomotion in complex fluids: Integral theorems

The biological fluids encountered by self-propelled cells display complex microstructures and rheology. We consider here the general problem of low-Reynolds number locomotion in a complex fluid. {Building on classical work on the transport of particles in viscoelastic fluids,} we demonstrate how to mathematically derive three integral theorems relating the arbitrary motion of an isolated organism to its swimming kinematics {in a non-Newtonian fluid}. These theorems correspond to three situations of interest, namely (1) squirming motion in a linear viscoelastic fluid, (2) arbitrary surface deformation in a weakly non-Newtonian fluid, and (3) small-amplitude deformation in an arbitrarily non-Newtonian fluid. Our final results, valid for a wide-class of {swimmer geometry,} surface kinematics and constitutive models, at most require mathematical knowledge of a series of Newtonian flow problems, and will be useful to quantity the locomotion of biological and synthetic swimmers in complex environments.

preprint2014arXiv

Mixing by microorganisms in stratified fluids

We examine the vertical mixing induced by the swimming of microorganisms at low Reynolds and Péclet numbers in a stably stratified ocean, and show that the global contribution of oceanic microswimmers to vertical mixing is negligible. We propose two approaches to estimating the mixing efficiency, $η$, or the ratio of the rate of potential energy creation to the total rate-of-working on the ocean by microswimmers. The first is based on scaling arguments and estimates $η$ in terms of the ratio between the typical organism size, $a$, and an intrinsic length scale for the stratified flow, $\ell = \left ( νκ/ N^2 \right )^{1/4}$, where $ν$ is the kinematic viscosity, $κ$ the diffusivity, and $N$ the buoyancy frequency. In particular, for small organisms in the relevant oceanic limit, $a / \ell \ll 1$, we predict the scaling $η\sim (a / \ell)^3$. The second estimate of $η$ is formed by solving the full coupled flow-stratification problem by modeling the swimmer as a regularized force dipole, and computing the efficiency numerically. Our computational results, which are examined for all ratios $a/\ell$, validate the scaling arguments in the limit $a / \ell \ll 1$ and further predict $η\approx 1.2 \left ( a / \ell \right )^3$ for vertical swimming and $η\approx 0.15 \left ( a / \ell \right )^3$ for horizontal swimming. These results, relevant for any stratified fluid rich in biological activity, imply that the mixing efficiency of swimming microorganisms in the ocean is at very most 8\% and is likely smaller by at least two orders of magnitude.

preprint2014arXiv

Optimal propulsive flapping in Stokes flows

Swimming fish and flying insects use the flapping of fins and wings to generate thrust. In contrast, microscopic organisms typically deform their appendages in a wavelike fashion. Since a flapping motion with two degrees of freedom is able, in theory, to produce net forces from a time-periodic actuation at all Reynolds number, we compute in this paper the optimal flapping kinematics of a rigid spheroid in a Stokes flow. The hydrodynamics for the force generation and energetics of the flapping motion is solved exactly. We then compute analytically the gradient of a flapping efficiency in the space of all flapping gaits and employ it to derive numerically the optimal flapping kinematics as a function of the shape of the flapper and the amplitude of the motion. The kinematics of optimal flapping are observed to depend weakly on the flapper shape and are very similar to the figure-eight motion observed in the motion of insect wings. Our results suggest that flapping could be a exploited experimentally as a propulsion mechanism valid across the whole range of Reynolds numbers.

preprint2014arXiv

Phoretic self-propulsion at finite Péclet numbers

Phoretic self-propulsion is a unique example of force- and torque-free motion on small scales. The classical framework describing the flow field around a particle swimming by self-diffusiophoresis neglects the advection of the solute field by the flow and assumes that the chemical interaction layer is thin compared to the particle size. In this paper we quantify and characterize the effect of solute advection on the phoretic swimming of a sphere. We first rigorously derive the regime of validity of the thin-interaction layer assumption at finite values of the Péclet number (Pe). Within this assumption, we solve computationally the flow around Janus phoretic particles and examine the impact of solute advection on propulsion and the flow created by the particle. We demonstrate that although advection always leads to a decrease of the swimming speed and flow stresslet at high values of the Péclet number, an increase can be obtained at intermediate values of Pe. This possible enhancement of swimming depends critically on the nature of the chemical interactions between the solute and the surface. We then derive an asymptotic analysis of the problem at small Pe allowing to rationalize our computational results. Our computational and theoretical analysis is accompanied by a parallel study of the role of reactive effects at the surface of the particle on swimming (Damköhler number).

preprint2014arXiv

Stochastic dynamics of active swimmers in linear flows

Most classical work on the hydrodynamics of low-Reynolds-number swimming addresses deterministic locomotion in quiescent environments. Thermal fluctuations in fluids are known to lead to a Brownian loss of the swimming direction. As most cells or synthetic swimmers are immersed in external flows, we consider theoretically in this paper the stochastic dynamics of a model active particle (a self-propelled sphere) in a steady general linear flow. The stochasticity arises both from translational diffusion in physical space, and from a combination of rotary diffusion and run-and-tumble dynamics in orientation space. We begin by deriving a general formulation for all components of the long-time mean square displacement tensor for a swimmer with a time-dependent swimming velocity and whose orientation decorrelates due to rotary diffusion alone. This general framework is applied to obtain the convectively enhanced mean-squared displacements of a steadily-swimming particle in three canonical linear flows (extension, simple shear, and solid-body rotation). We then show how to extend our results to the case where the swimmer orientation also decorrelates on account of run-and-tumble dynamics. Self-propulsion in general leads to the same long-time temporal scalings as for passive particles in linear flows but with increased coefficients. In the particular case of solid-body rotation, the effective long-time diffusion is the same as that in a quiescent fluid, and we clarify the lack of flow-dependence by briefly examining the dynamics in elliptic linear flows. By comparing the new active terms with those obtained for passive particles we see that swimming can lead to an enhancement of the mean-square displacements by orders of magnitude, and could be relevant for biological organisms or synthetic swimming devices in fluctuating environmental or biological flows.

preprint2014arXiv

The geometry and wetting of capillary folding

Capillary forces are involved in a variety of natural phenomena, ranging from droplet breakup to the physics of clouds. The forces from surface tension can also be exploited in industrial application provided the length scales involved are small enough. Recent experimental investigations showed how to take advantage of capillarity to fold planar structures into three-dimensional configurations by selectively melting polymeric hinges joining otherwise rigid shapes. In this paper we use theoretical calculations to quantify the role of geometry and fluid wetting on the final folded state. Considering folding in two and three dimensions, studying both hydrophilic and hydrophobic situations with possible contact angle hysteresis, and addressing the shapes to be folded to be successively infinite, finite, curved, kinked, elastic, we are able to derive an overview of the geometrical parameter space available for capillary folding.

preprint2014arXiv

The Optimal Swimming Sheet

Propulsion at microscopic scales is often achieved through propagating traveling waves along hair-like organelles called flagella. Taylor's two-dimensional swimming sheet model is frequently used to provide insight into problems of flagellar propulsion. We derive numerically the large-amplitude waveform of the two-dimensional swimming sheet that yields optimum hydrodynamic efficiency; the ratio of the squared swimming speed to the rate-of-working of the sheet against the fluid. Using the boundary element method, we show the optimal waveform is a front-back symmetric regularized cusp that is 25% more efficient than the optimal sine-wave. This optimal two-dimensional shape is smooth, qualitatively different from the kinked form of Lighthill's optimal three-dimensional flagellum, not predicted by small-amplitude theory, and different from the smooth circular-arc-like shape of active elastic filaments.

preprint2014arXiv

The other optimal Stokes drag profile

The lowest drag shape of fixed volume in Stokes flow has been known for some 40 years. It is front-back symmetric and similar to an American football with ends tangent to a cone of 60 degrees. The analogous convex axisymmetric shape of fixed surface area, which may be of interest for particle design in chemistry and colloidal science, is characterized in this paper. This "other" optimal shape has a surface vorticity proportional to the mean surface curvature, which is used with a local analysis of the flow near the tip to show that the front and rear ends are tangent to a cone of angle 30.8 degrees. Using the boundary element method, we numerically represent the shape by expanding its tangent angle in terms decaying odd Legendre modes, and show that it has 11.3% lower drag than a sphere of equal surface area, significantly more pronounced than for the fixed-volume optimal.

preprint2014arXiv

The passive diffusion of Leptospira interrogans

Motivated by recent experimental measurements, the passive diffusion of the bacterium \textit{Leptospira interrogans} is investigated theoretically. By approximating the cell shape as a straight helix and using the slender-body-theory approximation of Stokesian hydrodynamics, the resistance matrix of \textit{Leptospira} is first determined numerically. The passive diffusion of the helical cell is then obtained computationally using a Langevin formulation which is sampled in time in a manner consistent with the experimental procedure. Our results are in excellent quantitative agreement with the experimental results with no adjustable parameters.

preprint2014arXiv

The wobbling-to-swimming transition of rotated helices

A growing body of work aims at designing and testing micron-scale synthetic swimmers. One method, inspired by the locomotion of flagellated bacteria, consists of applying a rotating magnetic field to a rigid, helically-shaped, propeller attached to a magnetic head. When the resulting device, termed an artificial bacteria flagellum, is aligned perpendicularly to the applied field, the helix rotates and the swimmer moves forward. Experimental investigation of artificial bacteria flagella shows that at low frequency of the applied field, the axis of the helix does not align perpendicularly to the field but wobbles around the helix, with an angle increasing as the inverse of the field frequency. By numerical computations and asymptotic analysis, we provide a theoretical explanation for this wobbling behavior. We numerically demonstrate the wobbling-to-swimming transition as a function of the helix geometry and the dimensionless Mason number which quantifies the ratio of viscous to magnetic torques. We then employ an asymptotic expansion for near-straight helices to derive an analytical estimate for the wobbling angle allowing to rationalize our computations and past experimental results. These results can help guide future design of artificial helical swimmers.

preprint2014arXiv

Viscous pumping inspired by flexible propulsion

Fluid-suspended microorganisms have evolved different swimming and feeding strategies in order to cope with an environment dominated by viscous effects. For instance ciliated organisms rely on the collective motion of flexible appendices to move and feed. By performing a non-reciprocal motion, flexible filaments can produce a net propulsive force, or pump fluid, in the absence of inertia. Inspired by such fundamental concept, we propose a strategy to produce macroscopic pumping and mixing in creeping flow. We measure experimentally the net motion of a Newtonian viscous fluid induced by the reciprocal motion of a flapper. When the flapper is rigid no net motion is induced. In contrast, when the flapper is made of a flexible material, a net fluid pumping is measured. We quantify the effectiveness of this pumping strategy and show that optimal pumping is achieved when the length of the flapper is on the same order as the elasto-hydrodynamic penetration length. We finally discuss the possible applications of flexible impellers in mixing operations at low Reynolds numbers.

preprint2014arXiv

Waving transport and propulsion in a generalized Newtonian fluid

Cilia and flagella are hair-like appendages that protrude from the surface of a variety of eukaryotic cells and deform in a wavelike fashion to transport fluids and propel cells. Motivated by the ubiquity of non-Newtonian fluids in biology, we address mathematically the role of shear-dependent viscosities on both the waving flagellar locomotion and ciliary transport by metachronal waves. Using a two-dimensional waving sheet as model for the kinematics of a flagellum or an array of cilia, and allowing for both normal and tangential deformation of the sheet, we calculate the flow field induced by a small-amplitude deformation of the sheet in a generalized Newtonian Carreau fluid up to order four in the dimensionless waving amplitude. Shear-thinning and shear-thickening fluids are seen to always induce opposite effects. When the fluid is shear-thinning, the rate of working of the sheet against the fluid is always smaller than in the Newtonian fluid, and the largest gain is obtained for antiplectic metachronal waves. Considering a variety of deformation kinematics for the sheet, we further show that in all cases transport by the sheet is more efficiency in a shear-thinning fluid, and in most cases the transport speed in the fluid is also increased. Comparing the order of magnitude of the shear-thinning contributions with past work on elastic effects as well as the magnitude of the Newtonian contributions, our theoretical results suggest that the impact of shear-dependent viscosities on transport could play a major biological role.

preprint2013arXiv

Crawling scallop: Friction-based locomotion with one degree of freedom

Fluid-based locomotion at low Reynolds number is subject to the constraints of the scallop theorem, which dictate that body kinematics identical under a time-reversal symmetry (in particular, those with a single degree of freedom) cannot display locomotion on average. The implications of the theorem naturally compel one to ask whether similar symmetry constraints exist for locomotion in different environments. In this work we consider locomotion along a surface where forces are described by isotropic Coulomb friction. To address whether motions with a single degree of freedom can lead to transport, we analyze a model system consisting of two bodies whose separation distance undergoes periodic time variations. The behavior of the two-body system is entirely determined by the kinematic specification of their separation, the friction forces, and the mass of each body. We show that the constraints of the scallop theorem can be escaped in frictional media if two asymmetry conditions are met at the same time: the frictional forces of each body against the surface must be distinct and the time-variation of the body-body separation must vary asymmetrically in time (so quick-slow or slow-quick in the extension-contraction phases). Our results are demonstrated numerically and interpreted using asymptotic expansions.

preprint2013arXiv

Fluid elasticity increases the locomotion of flexible swimmers

We conduct experiments with flexible swimmers to address the impact of fluid viscoelasticity on their locomotion. The swimmers are composed of a magnetic head actuated in rotation by a frequency-controlled magnetic field and a flexible tail whose deformation leads to forward propulsion. We consider both viscous Newtonian and glucose-based Boger fluids with similar viscosities. We find that the elasticity of the fluid systematically enhances the locomotion speed of the swimmer, and that this enhancement increases with Deborah number. Using Particle Image Velocimetry to visualize the flow field, we find a significant difference in the amount of shear between the rear and leading parts of the swimmer head. We conjecture that viscoelastic normal stresses lead to a net elastic forces in the swimming direction and thus a faster swimming speed.

preprint2013arXiv

Fluid transport by active elastic membranes

A flexible membrane deforming its shape in time can self-propel in a viscous fluid. Alternatively, if the membrane is anchored, its deformation will lead to fluid transport. Past work in this area focused on situations where the deformation kinematics of the membrane were prescribed. Here we consider models where the deformation of the membrane is not prescribed, but instead the membrane is internally forced. Both the time-varying membrane shape, and the resulting fluid motion, result then from a balance between prescribed internal active stresses, internal passive resistance, and external viscous stresses. We introduce two specific models for such active internal forcing: one where a distribution of active bending moments is prescribed, and one where active inclusions exert normal stresses on the membrane by pumping fluid through it. In each case, we asymptotically calculate the membrane shape and the fluid transport velocities for small forcing amplitudes, and recover our results using scaling analysis.

preprint2013arXiv

Hydrodynamic fluctuations in confined emulsions

When an ensemble of particles interact hydrodynamically, they generically display large-scale transient structures such as swirls in sedimenting particles [1], or colloidal strings in sheared suspensions [2]. Understanding these nonequilibrium fluctuations is a very difficult problem, yet they are of great importance for a wide range of processes including pigment deposition in cosmetic or paint films, the transport of microfluidic droplets, ... All these samples concern rigidly confined fluids, which we consider in this paper. We address the collective dynamics of non-Brownian droplets cruising in a shallow microchannel. We provide a comprehensive characterization of their spatiotemporal density fluctuations. We show that density excitations freely propagate at all scales, and in all directions even though the particles are neither affected by potential forces nor by inertia. We introduce a theory which quantitatively accounts for our experimental findings. By doing so we demonstrate that the fluctuation spectrum of this nonequilibrium system is shaped by the combination of truly long-range hydrodynamic interactions and local collisions.

preprint2013arXiv

Low-Reynolds number swimming in a capillary tube

We use the boundary element method to study the low-Reynolds number locomotion of a spherical model microorganism in a circular tube. The swimmer propels itself by tangen- tial or normal surface motion in a tube whose radius is on the order of the swimmer size. Hydrodynamic interactions with the tube walls significantly affect the average swimming speed and power consumption of the model microorganism. In the case of swimming parallel to the tube axis, the locomotion speed is always reduced (resp. increased) for swimmers with tangential (resp. normal) deformation. In all cases, the rate of work nec- essary for swimming is increased by confinement. Swimmers with no force-dipoles in the far field generally follow helical trajectories, solely induced by hydrodynamic interactions with the tube walls, and in qualitative agreement with recent experimental observations for Paramecium. Swimmers of the puller type always display stable locomotion at a lo- cation which depends on the strength of their force dipoles: swimmers with weak dipoles (small α) swim in the centre of the tube while those with strong dipoles (large α) swim near the walls. In contrast, pusher swimmers and those employing normal deformation are unstable and end up crashing into the walls of the tube. Similar dynamics is observed for swimming into a curved tube. These results could be relevant for the future design of artificial microswimmers in confined geometries.

preprint2013arXiv

Micropropulsion and microrheology in complex fluids via symmetry breaking

Many biological fluids have polymeric microstructures and display non-Newtonian rheology. We take advantage of such nonlinear fluid behavior and combine it with geometrical symmetry-breaking to design a novel small-scale propeller able to move only in complex fluids. Its propulsion characteristics are explored numerically in an Oldroyd-B fluid for finite Deborah numbers while the small Deborah number limit is investigated analytically using a second-order fluid model. We then derive expressions relating the propulsion speed to the rheological properties of the complex fluid, allowing thus to infer the normal stress coefficients in the fluid from the locomotion of the propeller. Our simple mechanism can therefore be used either as a non-Newtonian micro-propeller or as a micro-rheometer.

preprint2013arXiv

Orientational order in concentrated suspensions of spherical microswimmers

We use numerical simulations to probe the dynamics of concentrated suspensions of spherical microswimmers interacting hydrodynamically. Previous work in the dilute limit predicted orientational instabilities of aligned suspensions for both pusher and puller swimmers, which we confirm computationally. Unlike previous work, we show that isotropic suspensions of spherical swimmers are also always unstable. Both types of initial conditions develop long-time polar order, of a nature which depends on the hydrodynamic signature of the swimmer but very weakly on the volume fraction up to very high volume fractions.

preprint2013arXiv

Shape of optimal active flagella

Many eukaryotic cells use the active waving motion of flexible flagella to self-propel in viscous fluids. However, the criteria governing the selection of particular flagellar waveforms among all possible shapes has proved elusive so far. To address this question, we derive computationally the optimal shape of an internally-forced periodic planar flagellum deforming as a travelling wave. The optimum is here defined as the shape leading to a given swimming speed with minimum energetic cost. To calculate the energetic cost though, we consider the irreversible internal power expanded by the molecular motors forcing the flagellum, only a portion of which ending up dissipated in the fluid. This optimisation approach allows us to derive a family of shapes depending on a single dimensionless number quantifying the relative importance of elastic to viscous effects: the Sperm number. The computed optimal shapes are found to agree with the waveforms observed on spermatozoon of marine organisms, thus suggesting that these eukaryotic flagella might have evolved to be mechanically optimal.

preprint2013arXiv

Spontaneous autophoretic motion of isotropic particles

Suspended colloidal particles interacting chemically with a solute are able to self-propel by autophoretic motion when they are asymmetrically patterned (Janus colloids). Here we demonstrate that the chemical anisotropy is not a necessary condition to achieve locomotion. The non linear interplay between surface osmotic flows and solute advection can produce spontaneous, and self-sustained motion of isotropic particles. We solve, for a spherical particle, the classical nonlinear autophoretic theoretical framework at arbitrary Péclet number. For a given set of material parameters, there exists a critical particle size, or equivalently a critical Peclet number, above which spontaneous autophoretic motion occurs. The flow induced by the particle further displays a hierarchy of instabilities associated with quantized critical Peclet numbers. Using numerical solutions of the full (unsteady) diffusiophoretic problem we confirm our analytical predictions and show that, above the instability threshold, the isotropic particles reach a steady swimming state with broken front-back symmetry in the concentration field and the hydrodynamic signature of a "pusher" swimmer. This instability to propulsion could be relevant to the high-throughput production of self-propelled particles.

preprint2013arXiv

Taylor's swimming sheet: Analysis and improvement of the perturbation series

In G.I. Taylor's historic paper on swimming microorganisms, a two dimensional sheet was proposed as a model for flagellated cells passing traveling waves as a means of locomotion. Using a perturbation series, Taylor computed swimming speeds up to fourth order in amplitude. Here we systematize that expansion so that it can be carried out formally to arbitrarily high order. The resultant series diverges for an order one value of the wave amplitude, but may be transformed into series with much improved convergence properties and which yield results comparing favorably to those obtained numerically via a boundary integral method for moderate and large values of the wave amplitudes.

preprint2012arXiv

Active and driven hydrodynamic crystals

Motivated by the experimental ability to produce monodisperse particles in microfluidic devices, we study theoretically the hydrodynamic stability of driven and active crystals. We first recall the theoretical tools allowing to quantify the dynamics of elongated particles in a confined fluid. In this regime hydrodynamic interactions between particles arise from a superposition of potential dipolar singularities. We exploit this feature to derive the equations of motion for the particle positions and orientations. After showing that all five planar Bravais lattices are stationary solutions of the equations of motion, we consider separately the case where the particles are passively driven by an external force, and the situation where they are self-propelling. We first demonstrate that phonon modes propagate in driven crystals, which are always marginally stable. The spatial structure of the eigenmodes depend solely on the symmetries of the lattices, and on the orientation of the driving force. For active crystals, the stability of the particle positions and orientations depends not only on the symmetry of the crystals but also on the perturbation wavelengths and on the crystal density. Unlike unconfined fluids, the stability of active crystals is independent of the nature of the propulsion mechanism at the single particle level. The square and rectangular lattices are found to be linearly unstable at short wavelengths provided the volume fraction of the crystals is high enough. Differently, hexagonal, oblique, and face-centered crystals are always unstable. Our work provides a theoretical basis for future experimental work on flowing microfluidic crystals.

preprint2012arXiv

Elastocapillary self-folding: buckling, wrinkling and collapse of floating filaments

When a flexible filament is confined to a fluid interface, the balance between capillary attraction, bending resistance, and tension from an external source can lead to a self-buckling instability. We perform an analysis of this instability and provide analytical formulae that compare favorably with the results of detailed numerical computations. The stability and long-time dynamics of the filament are governed by a single dimensionless elastocapillary number quantifying the ratio between capillary to bending stresses. Complex, folded filament configurations such as loops, needles, and racquet shapes may be reached at longer times, and long filaments can undergo a cascade of self-folding events.

preprint2012arXiv

Hydrodynamics of confined active fluids

We theoretically describe the dynamics of swimmer populations confined in thin liquid films. We first demonstrate that hydrodynamic interactions between confined swimmers only depend on their shape and are independent of their specific swimming mechanism. We also show that due to friction with the walls, confined swimmers do not reorient due to flow gradients but the flow field itself. We then quantify the consequences of these microscopic interaction rules on the large-scale hydrodynamics of isotropic populations. We investigate in details their stability and the resulting phase behavior, highlighting the differences with conventional active, three-dimensional suspensions. Two classes of polar swimmers are distinguished depending on their geometrical polarity. The first class gives rise to coherent directed motion at all scales whereas for the second class we predict the spontaneous formation of coherent clusters (swarms).

preprint2012arXiv

Hydrodynamics of self-propulsion near a boundary: predictions and accuracy of far-field approximations

The swimming trajectories of self-propelled organisms or synthetic devices in a viscous fluid can be altered by hydrodynamic interactions with nearby boundaries. We explore a multipole description of swimming bodies and provide a general framework for studying the fluid-mediated modifications to swimming trajectories. The validity of the far-field description is probed for a selection of model swimmers of varying geometry and propulsive activity by comparison with full numerical simulations. The reduced model is then used to deliver simple but accurate predictions of hydrodynamically generated wall attraction and pitching dynamics, and may help to explain a number of experimental results.

preprint2012arXiv

Hydrodynamics of the double-wave structure of insect spermatozoa flagella

In addition to conventional planar and helical flagellar waves, insect sperm flagella have also been observed to display a double-wave structure characterized by the presence of two superimposed helical waves. In this paper, we present a hydrodynamic investigation of the locomotion of insect spermatozoa exhibiting the double-wave structure, idealized here as superhelical waves. Resolving the hydrodynamic interactions with a non-local slender body theory, we predict the swimming kinematics of these superhelical swimmers based on experimentally collected geometric and kinematic data. Our consideration provides insight into the relative contributions of the major and minor helical waves to swimming; namely, propulsion is due primarily to the minor wave, with negligible contribution from the major wave. We also explore the dependence of the propulsion speed on geometric and kinematic parameters, revealing counter-intuitive results, particularly for the case when the minor and major helical structures are of opposite chirality.

preprint2012arXiv

Kinematics of the most efficient cilium

In a variety of biological processes, eukaryotic cells use cilia to transport flow. Although cilia have a remarkably conserved internal molecular structure, experimental observations report very diverse kinematics. To address this diversity, we determine numerically the kinematics and energetics of the most efficient cilium. Specifically, we compute the time-periodic deformation of a wall-bound elastic filament leading to transport of a surrounding fluid at minimum energetic cost, where the cost is taken to be the positive work done by all internal molecular motors. The optimal kinematics are found to strongly depend on the cilium bending rigidity through a single dimensionless number, the Sperm number, and closely resemble the two-stroke ciliary beating pattern observed experimentally.

preprint2012arXiv

Locomotion by tangential deformation in a polymeric fluid

In several biologically relevant situations, cell locomotion occurs in polymeric fluids with Weissenberg {number} larger than one. Here we present results of three-dimensional numerical simulations for the steady locomotion of a self-propelled body in a model polymeric (Giesekus) fluid at low Reynolds number. Locomotion is driven by steady tangential deformation at the surface of the body (so-called squirming motion). In the case of a spherical squirmer, we show that the swimming velocity is systematically less than that in a Newtonian fluid, with a minimum occurring for Weissenberg numbers of order one. The rate of work done by the swimmer always goes up compared to that occurring in the Newtonian solvent alone, but is always lower than the power necessary to swim in a Newtonian fluid with the same viscosity. The swimming efficiency, defined as the ratio between the rate of work necessary to pull the body at the swimming speed in the same fluid and the rate of work done by swimming, is found to always be increased in a polymeric fluid. Further analysis reveals that polymeric stresses break the Newtonian front-back symmetry in the flow profile around the body. In particular, a strong negative elastic wake is present behind the swimmer, which correlates with strong polymer stretching, and its intensity increases with Weissenberg number and viscosity contrasts. {The velocity induced by the squirmer is found to decay in space faster than} in a Newtonian flow, with a strong {dependence} on the polymer relaxation time and viscosity. Our computational results are also extended to prolate spheroidal swimmers and smaller polymer stretching are obtained for slender shapes compared to bluff swimmers. The swimmer with an aspect ratio of two is found to be the most hydrodynamically efficient.

preprint2012arXiv

Unsteady feeding and optimal strokes of model ciliates

The flow field created by swimming microorganisms not only enables their locomotion but also leads to advective transport of nutrients. In this paper we address analytically and computationally the link between unsteady feeding and unsteady swimming on a model microorganism, the spherical squirmer, actuating the fluid in a time-periodic manner. We start by performing asymptotic calculations at low Péclet number (Pe) on the advection-diffusion problem for the nutrients. We show that the mean rate of feeding as well as its fluctuations in time depend only on the swimming modes of the squirmer up to order Pe^(3/2), even when no swimming occurs on average, while the influence of non-swimming modes comes in only at order Pe^2. We also show that generically we expect a phase delay between feeding and swimming of 1/8th of a period. Numerical computations for illustrative strokes at finite Pe confirm quantitatively our analytical results linking swimming and feeding. We finally derive, and use, an adjoint-based optimization algorithm to determine the optimal unsteady strokes maximizing feeding rate for a fixed energy budget. The overall optimal feeder is always the optimal steady swimmer. Within the set of time-periodic strokes, the optimal feeding strokes are found to be equivalent to those optimizing periodic swimming for all values of the Péclet number, and correspond to a regularization of the overall steady optimal.

preprint2012arXiv

Viscous Marangoni propulsion

Marangoni propulsion is a form of locomotion wherein an asymmetric release of surfactant by a body located at the surface of a liquid leads to its directed motion. We present in this paper a mathematical model for Marangoni propulsion in the viscous regime. We consider the case of a thin rigid circular disk placed at the surface of a viscous fluid and whose perimeter has a prescribed concentration of an insoluble surfactant, to which the rest of its surface is impenetrable. Assuming a linearized equation of state between surface tension and surfactant concentration, we derive analytically the surfactant, velocity and pressure fields in the asymptotic limit of low Capillary, Peclet and Reynolds numbers. We then exploit these results to calculate the Marangoni propulsion speed of the disk. Neglecting the stress contribution from Marangoni flows is seen to over-predict the propulsion speed by 50%.

preprint2011arXiv

Comparative hydrodynamics of bacterial polymorphism

Most bacteria swim through fluids by rotating helical flagella which can take one of twelve distinct polymorphic shapes. The most common helical waveform is the "normal" form, used during forward swimming runs. To shed light on the prevalence of the normal form in locomotion, we gather all available experimental measurements of the various polymorphic forms and compute their intrinsic hydrodynamic efficiencies. The normal helical form is found to be the most hydrodynamically efficient of the twelve polymorphic forms by a significant margin - a conclusion valid for both the peritrichous and polar flagellar families, and robust to a change in the effective flagellum diameter or length. The hydrodynamic optimality of the normal polymorph suggests that, although energetic costs of locomotion are small for bacteria, fluid mechanical forces may have played a significant role in the evolution of the flagellum.

preprint2011arXiv

Enhanced diffusion by reciprocal swimming

Purcell's scallop theorem states that swimmers deforming their shapes in a time-reversible manner ("reciprocal" motion) cannot swim. Using numerical simulations and theoretical calculations we show here that in a fluctuating environment, reciprocal swimmers undergo, on time scales larger than that of their rotational diffusion, diffusive dynamics with enhanced diffusivities, possibly by orders of magnitude, above normal translational diffusion. Reciprocal actuation does therefore lead to a significant advantage over non-motile behavior for small organisms such as marine bacteria.

preprint2011arXiv

Extensibility enables locomotion under isotropic drag

Anisotropic viscous drag is usually believed to be a requirement for the low Reynolds number locomotion of slender bodies such as flagella and cilia. Here we show that locomotion under isotropic drag is possible for extensible slender bodies. After general considerations, a two-ring swimmer and a model dinoflagellate flagellum are studied analytically to illustrate how extensibility can be exploited for self-propulsion without drag anisotropy. This new degree of freedom could be useful for some complex swimmer geometries and locomotion in complex fluid environments where drag anisotropy is weak or even absent.

preprint2011arXiv

High-Speed Propulsion of Flexible Nanowire Motors: Theory and Experiments

Micro/nano-scale propulsion has attracted considerable recent attention due to its promise for biomedical applications such as targeted drug delivery. In this paper, we report on a new experimental design and theoretical modelling of high-speed fuel-free magnetically-driven propellers which exploit the flexibility of nanowires for propulsion. These readily prepared nanomotors display both high dimensional propulsion velocities (up to ~ 21 micrometer per second) and dimensionless speeds (in body lengths per revolution) when compared with natural microorganisms and other artificial propellers. Their propulsion characteristics are studied theoretically using an elastohydrodynamic model which takes into account the elasticity of the nanowire and its hydrodynamic interaction with the fluid medium. The critical role of flexibility in this mode of propulsion is illustrated by simple physical arguments, and is quantitatively investigated with the help of an asymptotic analysis for small-amplitude swimming. The theoretical predictions are then compared with experimental measurements and we obtain good agreement. Finally, we demonstrate the operation of these nanomotors in a real biological environment (human serum), emphasizing the robustness of their propulsion performance and their promise for biomedical applications.

preprint2011arXiv

Optimal feeding is optimal swimming for all Péclet numbers

Cells swimming in viscous fluids create flow fields which influence the transport of relevant nutrients, and therefore their feeding rate. We propose a modeling approach to the problem of optimal feeding at zero Reynolds number. We consider a simplified spherical swimmer deforming its shape tangentially in a steady fashion (so-called squirmer). Assuming that the nutrient is a passive scalar obeying an advection-diffusion equation, the optimal use of flow fields by the swimmer for feeding is determined by maximizing the diffusive flux at the organism surface for a fixed rate of energy dissipation in the fluid. The results are obtained through the use of an adjoint-based numerical optimization implemented by a Legendre polynomial spectral method. We show that, to within a negligible amount, the optimal feeding mechanism consists in putting all the energy expended by surface distortion into swimming - so-called treadmill motion - which is also the solution maximizing the swimming efficiency. Surprisingly, although the rate of feeding depends strongly on the value of the Péclet number, the optimal feeding stroke is shown to be essentially independent of it, which is confirmed by asymptotic analysis. Within the context of steady actuation, optimal feeding is therefore found to be equivalent to optimal swimming for all Péclet numbers.

preprint2011arXiv

Synchronization of flexible sheets

When swimming in close proximity, some microorganisms such as spermatozoa synchronize their flagella. Previous work on swimming sheets showed that such synchronization requires a geometrical asymmetry in the flagellar waveforms. Here we inquire about a physical mechanism responsible for such symmetry-breaking in nature. Using a two-dimensional model, we demonstrate that flexible sheets with symmetric internal forcing, deform when interacting with each other via a thin fluid layer in such a way as to systematically break the overall waveform symmetry, thereby always evolving to an in-phase conformation where energy dissipation is minimized. This dynamics is shown to be mathematically equivalent to that obtained for prescribed waveforms in viscoelastic fluids, emphasizing the crucial role of elasticity in symmetry-breaking and synchronization.

preprint2010arXiv

A two-dimensional model of low-Reynolds number swimming beneath a free surface

Biological organisms swimming at low Reynolds number are often influenced by the presence of rigid boundaries and soft interfaces. In this paper we present an analysis of locomotion near a free surface with surface tension. Using a simplified two-dimensional singularity model, and combining a complex variable approach with conformal mapping techniques, we demonstrate that the deformation of a free surface can be harnessed to produce steady locomotion parallel to the interface. The crucial physical ingredient lies in the nonlinear hydrodynamic coupling between the disturbance flow created by the swimmer and the free boundary problem at the fluid surface.

preprint2010arXiv

Efficiency optimization and symmetry-breaking in a model of ciliary locomotion

A variety of swimming microorganisms, called ciliates, exploit the bending of a large number of small and densely-packed organelles, termed cilia, in order to propel themselves in a viscous fluid. We consider a spherical envelope model for such ciliary locomotion where the dynamics of the individual cilia are replaced by that of a continuous overlaying surface allowed to deform tangentially to itself. Employing a variational approach, we determine numerically the time-periodic deformation of such surface which leads to low-Reynolds locomotion with minimum rate of energy dissipation (maximum efficiency). Employing both Lagrangian and Eulerian points of views, we show that in the optimal swimming stroke, individual cilia display weak asymmetric beating, but that a significant symmetry-breaking occurs at the organism level, with the whole surface deforming in a wave-like fashion reminiscent of metachronal waves of biological cilia. This wave motion is analyzed using a formal modal decomposition, is found to occur in the same direction as the swimming direction, and is interpreted as due to a spatial distribution of phase-differences in the kinematics of individual cilia. Using additional constrained optimizations, as well as a constructed analytical ansatz, we derive a complete optimization diagram where all swimming efficiencies, swimming speeds, and amplitude of surface deformation can be reached, with the mathematically optimal swimmer, of efficiency one half, being a singular limit. Biologically, our work suggests therefore that metachronal waves may allow cilia to propel cells forward while reducing the energy dissipated in the surrounding fluid.

preprint2010arXiv

Hydrodynamic friction of fakir-like super-hydrophobic surfaces

A fluid droplet located on a super-hydrophobic surface makes contact with the surface only at small isolated regions, and is mostly in contact with the surrounding air. As a result, a fluid in motion near such a surface experiences very low friction, and super-hydrophobic surfaces display strong drag-reduction in the laminar regime. Here we consider theoretically a super-hydrophobic surface composed of circular posts (so called fakir geometry) located on a planar rectangular lattice. Using a superposition of point forces with suitably spatially-dependent strength, we derive the effective surface slip length for a planar shear flow on such a fakir surface as the solution to an infinite series of linear equations. In the asymptotic limit of small surface coverage by the posts, the series can be interpreted as Riemann sums, and the slip length can be obtained analytically. For posts on a square lattice, our analytical results are in excellent quantitative agreement with previous numerical computations.

preprint2010arXiv

Jet propulsion without inertia

A body immersed in a highly viscous fluid can locomote by drawing in and expelling fluid through pores at its surface. We consider this mechanism of jet propulsion without inertia in the case of spheroidal bodies, and derive both the swimming velocity and the hydrodynamic efficiency. Elementary examples are presented, and exact axisymmetric solutions for spherical, prolate spheroidal, and oblate spheroidal body shapes are provided. In each case, entirely and partially porous (i.e. jetting) surfaces are considered, and the optimal jetting flow profiles at the surface for maximizing the hydrodynamic efficiency are determined computationally. The maximal efficiency which may be achieved by a sphere using such jet propulsion is 12.5%, a significant improvement upon traditional flagella-based means of locomotion at zero Reynolds number. Unlike other swimming mechanisms which rely on the presentation of a small cross section in the direction of motion, the efficiency of a jetting body at low Reynolds number increases as the body becomes more oblate, and limits to approximately 162% in the case of a flat plate swimming along its axis of symmetry. Our results are discussed in the light of slime extrusion mechanisms occurring in many cyanobacteria.

preprint2010arXiv

Passive hydrodynamic synchronization of two-dimensional swimming cells

Spermatozoa flagella are known to synchronize when swimming in close proximity. We use a model consisting of two-dimensional sheets propagating transverse waves of displacement to demonstrate that fluid forces lead to such synchronization passively. Using two distinct asymptotic descriptions (small amplitude and long wavelength), we derive the synchronizing dynamics analytically for ar- bitrarily shaped waveforms in Newtonian fluids, and show that phase locking will always occur for sufficiently asymmetric shapes. We characterize the effect of the geometry of the waveforms and the separation between the swimmers on the synchronizing dynamics, the final stable conformations, and the energy dissipated by the cells. For two closely-swimming cells, synchronization always oc- curs at the in-phase or opposite-phase conformation, depending solely on the geometry of the cells. In contrast, the work done by the swimmers is always minimized at the in-phase conformation. As the swimmers get further apart, additional fixed points arise at intermediate values of the relative phase. In addition, computations for large-amplitude waves using the boundary integral method reveal that the two asymptotic limits capture all the relevant physics of the problem. Our results provide a theoretical framework to address other hydrodynamic interactions phenomena relevant to populations of self-propelled organisms.

preprint2010arXiv

Propulsion by passive filaments and active flagella near boundaries

Confinement and wall effects are known to affect the kinematics and propulsive characteristics of swimming microorganisms. When a solid body is dragged through a viscous fluid at constant velocity, the presence of a wall increases fluid drag, and thus the net force required to maintain speed has to increase. In contrast, recent optical trapping experiments have revealed that the propulsive force generated by human spermatozoa is decreased by the presence of boundaries. Here, we use a series of simple models to analytically elucidate the propulsive effects of a solid boundary on passively actuated filaments and model flagella. For passive flexible filaments actuated periodically at one end, the presence of the wall is shown to increase the propulsive forces generated by the filaments in the case of displacement-driven actuation, while it decreases the force in the case of force-driven actuation. In the case of active filaments as models for eukaryotic flagella, we demonstrate that the manner in which a solid wall affects propulsion cannot be known a priori, but is instead a nontrivial function of the flagellum frequency, wavelength, its material characteristics, the manner in which the molecular motors self-organize to produce oscillations (prescribed activity model or self-organized axonemal beating model), and the boundary conditions applied experimentally to the tethered flagellum. In particular, we show that in some cases, the increase in fluid friction induced by the wall can lead to a change in the waveform expressed by the flagella which results in a decrease in their propulsive force.

preprint2010arXiv

Pumping by flapping in a viscoelastic fluid

In a world without inertia, Purcell's scallop theorem states that in a Newtonian fluid a time-reversible motion cannot produce any net force or net flow. Here we consider the extent to which the nonlinear rheological behavior of viscoelastic fluids can be exploited to break the constraints of the scallop theorem in the context of fluid pumping. By building on previous work focusing on force generation, we consider a simple, biologically-inspired geometrical example of a flapper in a polymeric (Oldroyd-B) fluid, and calculate asymptotically the time-average net fluid flow produced by the reciprocal flapping motion. The net flow occurs at fourth order in the flapping amplitude, and suggests the possibility of transporting polymeric fluids using reciprocal motion in simple geometries even in the absence of inertia. The induced flow field and pumping performance are characterized and optimized analytically. Our results may be useful in the design of micro-pumps handling complex fluids.

preprint2010arXiv

Stokesian jellyfish: Viscous locomotion of bilayer vesicles

Motivated by recent advances in vesicle engineering, we consider theoretically the locomotion of shape-changing bilayer vesicles at low Reynolds number. By modulating their volume and membrane composition, the vesicles can be made to change shape quasi-statically in thermal equilibrium. When the control parameters are tuned appropriately to yield periodic shape changes which are not time-reversible, the result is a net swimming motion over one cycle of shape deformation. For two classical vesicle models (spontaneous curvature and bilayer coupling), we determine numerically the sequence of vesicle shapes through an enthalpy minimization, as well as the fluid-body interactions by solving a boundary integral formulation of the Stokes equations. For both models, net locomotion can be obtained either by continuously modulating fore-aft asymmetric vesicle shapes, or by crossing a continuous shape-transition region and alternating between fore-aft asymmetric and fore-aft symmetric shapes. The obtained hydrodynamic efficiencies are similar to that of other low Reynolds number biological swimmers, and suggest that shape-changing vesicles might provide an alternative to flagella-based synthetic microswimmers.

preprint2009arXiv

Shaking-induced motility in suspensions of soft active particles

We investigate theoretically the collective dynamics of soft active particles living in a viscous fluid. We focus on a minimal model for active but non-motile particles consisting of $N>1$ elastic dimers deformed by active stresses and interacting hydrodynamically. We first derive a set of effective equations of motion for the positions of the particles. We then exploit these equations in two experimentally-relevant cases: uncorrelated random internal stresses, and uniform monochromatic external shaking. In both cases, we show that small groups of intrinsically non-motile particles can display non-trivial modes of locomotion resulting from the hydrodynamic correlations between the particle-conformation fluctuations. In addition, we demonstrate that a coherent shaking yields spatial ordering in suspension of soft particles interacting solely through the fluid.

preprint2009arXiv

The friction of a mesh-like super-hydrophobic surface

When a liquid droplet is located above a super-hydrophobic surface, it only barely touches the solid portion of the surface, and therefore slides very easily on it. More generally, super-hydrophobic surfaces have been shown to lead to significant reduction of viscous friction in the laminar regime, so it is of interest to quantify their effective slipping properties as a function of their geometric characteristics. Most previous studies have considered flows bounded by arrays of either long grooves, or isolated solid pillars on an otherwise flat solid substrate, and for which therefore the surrounding air constitutes the continuous phase. Here we consider instead the case where the super-hydrophobic surface is made of isolated holes in an otherwise continuous no-slip surface, and specifically focus on the mesh-like geometry recently achieved experimentally. We present an analytical method to calculate the friction of such a surface in the case where the mesh is thin. The results for the effective slip length of the surface are computed, compared to simple estimates, and a practical fit is proposed displaying a logarithmic dependance on the area fraction of the solid surface.

preprint2009arXiv

The long-time dynamics of two hydrodynamically-coupled swimming cells

Swimming micro-organisms such as bacteria or spermatozoa are typically found in dense suspensions, and exhibit collective modes of locomotion qualitatively different from that displayed by isolated cells. In the dilute limit where fluid-mediated interactions can be treated rigorously, the long-time hydrodynamics of a collection of cells result from interactions with many other cells, and as such typically eludes an analytical approach. Here we consider the only case where such problem can be treated rigorously analytically, namely when the cells have spatially confined trajectories, such as the spermatozoa of some marine invertebrates. We consider two spherical cells swimming, when isolated, with arbitrary circular trajectories, and derive the long-time kinematics of their relative locomotion. We show that in the dilute limit where the cells are much further away than their size, and the size of their circular motion, a separation of time scale occurs between a fast (intrinsic) swimming time, and a slow time where hydrodynamic interactions lead to change in the relative position and orientation of the swimmers. We perform a multiple-scale analysis and derive the effective dynamical system - of dimension two - describing the long-time behavior of the pair of cells. We show that the system displays one type of equilibrium, and two types of rotational equilibrium, all of which are found to be unstable. A detailed mathematical analysis of the dynamical systems further allows us to show that only two cell-cell behaviors are possible in the limit of $t\to\infty$, either the cells are attracted to each other (possibly monotonically), or they are repelled (possibly monotonically as well), which we confirm with numerical computations.

preprint2009arXiv

The optimal elastic flagellum

Motile eukaryotic cells propel themselves in viscous fluids by passing waves of bending deformation down their flagella. An infinitely long flagellum achieves a hydrodynamically optimal low-Reynolds number locomotion when the angle between its local tangent and the swimming direction remains constant along its length. Optimal flagella therefore adopt the shape of a helix in three dimensions (smooth) and that of a sawtooth in two dimensions (non-smooth). Physically, biological organisms (or engineered micro-swimmers) must expend internal energy in order to produce the waves of deformation responsible for the motion. Here we propose a physically-motivated derivation of the optimal flagellum shape. We determine analytically and numerically the shape of the flagellar wave which leads to the fastest swimming while minimizing an appropriately-defined energetic expenditure. Our novel approach is to define an energy which includes not only the work against the surrounding fluid, but also (1) the energy stored elastically in the bending of the flagellum, (2) the energy stored elastically in the internal sliding of the polymeric filaments which are responsible for the generation of the bending waves (microtubules), and (3) the viscous dissipation due to the presence of an internal fluid. This approach regularizes the optimal sawtooth shape for two-dimensional deformation at the expense of a small loss in hydrodynamic efficiency. The optimal waveforms of finite-size flagella are shown to depend upon a competition between rotational motions and bending costs, and we observe a surprising bias towards half-integer wave-numbers. Their final hydrodynamic efficiencies are above 6%, significantly larger than those of swimming cells, therefore indicating available room for further biological tuning.

preprint2009arXiv

The transient swimming of a waving sheet

Small-scale locomotion plays an important role in biology. Different modelling approaches have been proposed in the past. The simplest model is an infinite inextensible two-dimensional waving sheet, {originally introduced by Taylor}, which serves as an idealized geometrical model for both spermatozoa locomotion and ciliary transport in Stokes flow. Here we complement classic steady-state calculations by deriving the transient low-Reynolds number swimming speed of such a waving sheet when starting from rest (small-amplitude initial-value problem). We also determine the transient fluid flow in the `pumping' setup where the sheet is not free to move but instead generates a net fluid flow around it. The time scales for these two problems, which in general govern transient effects in transport and locomotion at low Reynolds numbers, are also derived using physical arguments.

preprint2009arXiv

Two-dimensional flagellar synchronization in viscoelastic fluids

Experimental studies have demonstrated that spermatozoa synchronize their flagella when swimming in close proximity. In a Newtonian fluid, it was shown theoretically that such synchronization arises passively due to hydrodynamic forces between the two swimmers if their waveforms exhibit a front-back geometrical asymmetry. Motivated by the fact that most biological fluids possess a polymeric microstructure, we address here synchronization in a viscoelastic fluid analytically. Using a two-dimensional infinite sheet model we show that the presence of polymeric stresses removes the geometrical asymmetry constraint, and therefore even symmetric swimmers synchronize. Such synchronization occurs on asymptotically faster time scales than in a Newtonian fluid, and the swimmers are seen to be driven into a stable in-phase conformation minimizing the energy dissipated in the surrounding fluid.

preprint2005arXiv

Swimming in circles: Motion of bacteria near solid boundaries

Near a solid boundary, E. coli swims in clockwise circular motion. We provide a hydrodynamic model for this behavior. We show that circular trajectories are natural consequences of force-free and torque-free swimming, and the hydrodynamic interactions with the boundary, which also leads to a hydrodynamic trapping of the cells close to the surface. We compare the results of the model with experimental data and obtain reasonable agreement. In particular, we show that the radius of curvature of the trajectory increases with the length of the bacterium body.

preprint2004arXiv

Shear-dependent apparent slip on hydrophobic surfaces: The Mattress Model

Recent experiments (Zhu & Granick (2001) Phys. Rev. Lett. 87 096105) have measured a large shear dependent fluid slip at partially wetting fluid-solid surfaces. We present a simple model for such slip, motivated by the recent observations of nanobubbles on hydrophobic surfaces. The model considers the dynamic response of bubbles to change in hydrodynamic pressure due to the oscillation of a solid surface. Both the compression and diffusion of gas in the bubbles decrease the force on the oscillating surface by a ``leaking mattress'' effect, thereby creating an apparent shear-dependent slip. With bubbles similar to those observed by atomic force microscopy to date, the model is found to lead to force decreases consistent with the experimental measurements of Zhu & Granick.