Search

We present an asymptotic theory for solving the dynamics of slender autophoretic loops and knots. Our formulation is valid for non-intersecting 3D centrelines, with arbitrary chemical patterning and varying (circular) cross-sectional radius, allowing a broad class of slender active loops and knots to be studied. The theory is amenable to closed-form solutions in simpler cases, allowing us to analytically derive the swimming speed of chemically patterned tori, and the pumping strength (stresslet) of a uniformly active slender torus. Using simple numerical solutions of our asymptotic equations, we then elucidate the behaviour of many exotic active particle geometries, such as a bumpy uniformly active torus that spins and a Janus trefoil knot, which rotates as it swims forwards., Comment: 24 pages with 7 figures, accepted in Physical Review Fluids

Filaments are ubiquitous within the microscopic world. They occur frequently in both biological and industrial environments and display varied and rich dynamics. Their wide range of applications has spurred the development of a special branch of asymptotics focused on the behaviour of filaments, called slender-body theory. Slender-body theories are typically computationally efficient and focus on the mechanics of an isolated fibre that is not too curved. However, slender-body theories that work beyond these standard limits are needed to explore more complex systems. Recently, we developed tubular-body theory for slow viscous flows, an approach similar to slender-body theory that allows the hydrodynamic traction on any isolated cable-like body in a highly viscous fluid to be determined exactly. In this paper, we extend tubular-body theory to model filaments near plane interfaces by performing an similar expansion on the single-layer boundary integral equations for bodies by a plane interface. In the derivation of the new theory, called tubular-body theory for interfaces, we established a criteria for the convergence of the tubular-body theory series representation, before comparing the result to boundary integral simulations for a prolate spheroid by a wall. The tubular-body theory for interfaces simulations are found to capture the lubrication effects when close to the plane wall. Finally we simulate the hydrodynamics of a helix beneath a free interface and a plane wall to demonstrate the broad applicability of the technique.

The slow viscous flow through a doubly-periodic array of cylinders does not have an analytical solution. However, as a reduced model for the flow within fibrous porous media, this solution is important for many real-world systems. We asymptotically determine the flow around a doubly-periodic array of infinite slender cylinders, by placing doubly-periodic two-dimensional singularity solutions within the cylinder and expanding the no-slip condition on the cylinder's surface in powers of the cylinder radius. The asymptotic solution provides a closed-form estimate for the flow and forces as a function of the radius and the dimensions of the cell. The force is compared to results from lattice-Boltzmann simulations of low-Reynolds-number flows in the same geometry, and the accuracy of the no-slip condition on the surface of the cylinder, predicted by the asymptotic theory, is checked. Finally, the behaviour of the flow, flux, force and effective permeability of the cell is investigated as a function of the geometric parameters. The structure of the asymptotic permeability is consistent with other models for the flow parallel to an array of rods. These models could be used to help understand the flows within porous systems composed of fibres and systems involving periodic arrays such as deterministic lateral displacement.

Cable-like bodies play a key role in many interdisciplinary systems but are hard to simulate. Asymptotic theories, called slender-body theories, are effective but apply in specific regimes and can be hard to extend beyond leading order. In this letter we develop an exact slender-body-like theory for the surface traction of cable-like bodies in viscous flow. This theory expresses the traction as a series of solutions to a well-behaved one-dimensional Fredholm integral equation of the second kind. This process can be simply generalised to other systems., Comment: 7 pages 3 figures

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.

Slender-body approximations have been successfully used to explain many phenomena in low-Reynolds number fluid mechanics. These approximations typically use a line of singularity solutions to represent the flow. These singularities can be difficult to implement numerically because they diverge at their origin. Hence people have regularized these singularities to overcome this issue. This regularization blurs the force over a small blob therefore removing the divergent behaviour. However it is unclear how best to regularize the singularities to minimize errors. In this paper we investigate if a line of regularized Stokeslets can describe the flow around a slender body. This is achieved by comparing the asymptotic behaviour of the flow from the line of regularized Stokeslets with the results from slender-body theory. We find that the flow far from the body can be captured if the regularization parameter is proportional to the radius of the slender body. This is consistent with what is assumed in numerical simulations and provides a choice for the proportionality constant. However more stringent requirements must be placed on the regularization blob to capture the near field flow outside a slender body. This inability to replicate the local behaviour indicates that many regularizations cannot satisfy the non-slip boundary conditions on the bodies surface to leading order, with one of the most commonly used regularizations showing an angular dependency of velocity along any cross section. This problem can be overcome with compactly supported blobs { and we construct one such example blob which could be effectively used to simulate the flow around a slender body

The viscous drag on a slender rod by a wall is important to many biological and industrial systems. This drag critically depends on the separation between the rod and the wall and can be approximated asymptotically in specific regimes, namely far from, or very close to, the wall, but is typically determined numerically for general separations. In this note we determine an asymptotic representation of the local drag for a slender rod parallel to a wall which is valid for all separations. This is possible through matching the behaviour of a rod close to the wall and a rod far from the wall. We show that the leading order drag in both these regimes has been known since 1981 and that they can used to produce a composite representation of the drag which is valid for all separations. This is in contrast to a sphere above a wall, where no simple uniformly valid representation exists. We estimate the error on this composite representation as the separation increases, discuss how the results could be used as resistive-force theory and demonstrate their use on a two-hinged swimmer above a wall., Comment: 25 pages 6 figures

The geometric phase techniques for swimming in viscous flows express the net displacement of a swimmer as a path integral of a field in configuration space. This representation can be transformed into an area integral for simple swimmers using Stokes theorem. Since this transformation applies for any loop, the integrand of this area integral can be used to help design these swimmers. However, the extension of this Stokes theorem technique to more complicated swimmers is hampered by problems with variables that do not commute and by how to visualise and understand the higher dimensional spaces. In this paper, we develop a treatment for each of these problems, thereby allowing the displacement of general swimmers in any environment to be designed and understood similarly to simple swimmers. The net displacement arising from non-commuting variables is tackled by embedding the integral into a higher dimensional space, which can then be visualised through a suitability constructed surface. These methods are developed for general swimmers and demonstrated on {three} benchmark examples: Purcell's two-hinged swimmer, an axisymmetric squirmer in free space {and an axisymmetric squirmer approaching a free interface}. We show in particular that for swimmers with more than two modes of deformation, there exists an infinite set of strokes that generate each net displacement. Hence, in the absence of additional restrictions, general microscopic swimmers do not have a single stroke that maximises their displacement., Comment: 42 pages, 11 figures

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., Comment: 37 pages, 5 main figures

Since their development in 2001, regularised stokeslets have become a popular numerical tool for low-Reynolds number flows since the replacement of a point force by a smoothed blob overcomes many computational difficulties associated with flow singularities (Cortez, 2001, \textit{SIAM J. Sci. Comput.} \textbf{23}, 1204). The physical changes to the flow resulting from this process are, however, unclear. In this paper, we analyse the flow induced by general regularised stokeslets. An explicit formula for the flow from any regularised stokeslet is first derived, which is shown to simplify for spherically symmetric blobs. Far from the centre of any regularised stokeslet we show that the flow can be written in terms of an infinite number of singularity solutions provided the blob decays sufficiently rapidly. This infinite number of singularities reduces to a point force and source dipole for spherically symmetric blobs. Slowly-decaying blobs induce additional flow resulting from the non-zero body forces acting on the fluid. We also show that near the centre of spherically symmetric regularised stokeslets the flow becomes isotropic, which contrasts with the flow anisotropy fundamental to viscous systems. The concepts developed are used to { identify blobs that reduce regularisation errors. These blobs contain regions of negative force in order to counter the flows produced in the regularisation process, but still retain a form convenient for computations., Comment: 32 pages, 7 figures

The current understanding of motility through body shape deformation of microorganisms and the knowledge of fluid flows at the microscale provides ample examples for mimicry and design of soft microrobots. In this work, a two-dimensional spiral is presented that is capable of rotating by non-reciprocal curling deformations. The body of the microswimmer is a ribbon consisting of a thermo-responsive hydrogel bilayer with embedded plasmonic gold nanorods. Such a system allows fast local photothermal heating and non-reciprocal bending deformation of the hydrogel bilayer under non-equilibrium conditions. We show that the spiral acts as a spring capable of large deformations thanks to its low stiffness, which is tunable by the swelling degree of the hydrogel and the temperature. Tethering the ribbon to a freely rotating microsphere enables rotational motion of the spiral by stroboscopic irradiation. The efficiency of the rotor is estimated using resistive force theory for Stokes flow. The present research demonstrates microscopic locomotion by the shape change of a spiral and may find applications in the field of microfluidics, or soft micro-robotics., Comment: 5 Figures, Version 1.0

Control on microscopic scales depends critically on our ability to manipulate interactions with different physical fields. The creation of micro-machines therefore requires us to understand how multiple fields, such as surface capillary or electro-magnetic, can be used to produce predictable behaviour. Recently, a spinning micro-raft system was developed that exhibited both static and dynamic self-assembly [Wang et al. (2017) Sci. Adv. 3, e1602522]. These rafts employed both capillary and magnetic interactions and, at a critical driving frequency, would suddenly change from stable orbital patterns to static assembled structures. In this paper, we explain the dynamics of two interacting micro-rafts through a combination of theoretical models and experiments. This is first achieved by identifying the governing physics of the orbital patterns, the assembled structures, and the collapse separately. We find that the orbital patterns are determined by the short range capillary interactions between the disks, while the explanations of the other two behaviours only require the capillary far field. Finally we combine the three models to explain the dynamics of a new micro-raft experiment., Comment: 10 pages 7 figures

Many microorganisms and artificial microswimmers use helical appendages in order to generate locomotion. Though often rotated so as to produce thrust, some species of bacteria such Spiroplasma, Rhodobacter sphaeroides and Spirochetes induce movement by deforming a helical-shaped body. Recently, artificial devices have been created which also generate motion by deforming their helical body in a non-reciprocal way (Mourran et al., Adv. Mater., 29, 1604825, 2017). Inspired by these systems, we investigate the transport of a deforming helix within a viscous fluid. Specifically, we consider a swimmer that maintains a helical centreline and a single handedness while changing its helix radius, pitch and wavelength uniformly across the body. We first discuss how a deforming helix can create a non-reciprocal translational and rotational swimming stroke and identify its principle direction of motion. We then determine the leading-order physics for helices with small helix radius before considering the general behaviour for different configuration parameters and how these swimmers can be optimised. Finally, we explore how the presence of walls, gravity, and defects in the centreline allow the helical device to break symmetries, increase its speed, and generate transport in directions not available to helices in bulk fluids.

The incompressible Stokes equations can classically be recast in a boundary integral (BI) representation, which provides a general method to solve low-Reynolds number problems analytically and computationally. Alternatively, one can solve the Stokes equations by using an appropriate distribution of flow singularities of the right strength within the boundary, a method particularly useful to describe the dynamics of long slender objects for which the numerical implementation of the BI representation becomes cumbersome. While the BI approach is a mathematical consequence of the Stokes equations, the singularity method involves making judicious guesses that can only be justified a posteriori. In this paper we use matched asymptotic expansions to derive an algebraically accurate slender-body theory directly from the BI representation able to handle arbitrary surface velocities and surface tractions. This expansion procedure leads to sets of uncoupled linear equations and to a single one-dimensional integral equation identical to that derived by Keller and Rubinow (1976) and Johnson (1979) using the singularity method. Hence we show that it is a mathematical consequence of the BI approach that the leading-order flow around a slender body can be represented using a distribution of singularities along its centreline. Furthermore when derived from either the single-layer or double-layer modified BI representation, general slender solutions are only possible in certain types of flow, in accordance with the limitations of these representations.

The low-Reynolds number hydrodynamics of slender ribbons is accurately captured by slender-ribbon theory, an asymptotic solution to the Stokes equation which assumes that the three length scales characterising the ribbons are well separated. We show in this paper that the force distribution across the width of an isolated ribbon located in a infinite fluid can be determined analytically, irrespective of the ribbon's shape. This, in turn, reduces the surface integrals in the slender-ribbon theory equations to a line integral analogous to the one arising in slender-body theory to determine the dynamics of filaments. This result is then used to derive analytical solutions to the motion of a rigid plate ellipsoid and a ribbon torus and to propose a ribbon resistive-force theory, thereby extending the resistive-force theory for slender filaments.

Experiments involving the two-dimensional passive diffusion of colloidal boomerangs tracked off their centre of mobility have shown striking non-Gaussian tails in their probability distribution function [Chakrabarty et al., Soft Matter 12, 4318 (2016)]. This in turn can lead to anomalous diffusion characteristics, including mean drift. In this paper, we develop a general theoretical explanation for these measurements. The idea relies on calculating the two-dimensional probability densities at the centre of mobility of the particle, where all distributions are Gaussian, and then transforming them to a different reference point. Our model clearly captures the experimental results, without any fitting parameters, and demonstrates that the one-dimensional probability distributions may also exhibit strongly non-Gaussian tops. These results indicate that the choice of tracking point can cause a considerable departure from Gaussian statistics, potentially causing some common modelling techniques to fail., Comment: 7 pages, 4 figures, submitted. LK and ML contributed equally

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"., Comment: 8 pages, 6 figures

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.

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.

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\"obius 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.

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.

We derive a governing equation for a Kelvin wave supported on a vortex line in a Bose-Einstein condensate, in a rotating cylindrically symmetric parabolic trap. From this solution the Kelvin wave dispersion relation is determined. In the limit of an oblate trap and in the absence of longitudinal trapping our results are consistent with previous work. We show that the derived Kelvin wave dispersion in the general case is in quantitative agreement with numerical calculations of the Bogoliubov spectrum and offer a significant improvement upon previous analytical work., Comment: 5 pages with 1 figure

We derive a set of equations that describe the shape and behaviour of a single perturbed vortex line in a Bose-Einstein condensate. Through the use of a matched asymptotic expansion and a unique coordinate transform a relation for a vortex's velocity, anywhere along the line, is found in terms of the trapping, rotation, and distortion of the line at that location. This relation is then used to find a set of differential equations that give the line's specific shape and motion. This work corrects a previous similar derivation by Anatoly A. Svidzinsky and Alexander L. Fetter [Phys. Rev. A \textbf{62}, 063617 (2000)], and enables a comparison with recent numerical results., Comment: 12 pages with 3 figures