Search

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.

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 order. Using the neighbor distance as the information-bearing variable, we demonstrate the links among information, structure, and interactions. We establish a direct link between information and structure without using explicit knowledge of interactions. Last, 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., Science Advances, 8 (2), ISSN:2375-2548

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., 32 pages, 7 figures

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

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 fields, 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., Sci. Adv., 2017, 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., European Union (European Union); European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme; Max Planck Society; Alexander von Humboldt Foundation

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 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 and thereby removing 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 no-slip boundary conditions on the body’s surface to leading order, with one of the most commonly used blobs 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 can be effectively used to simulate the flow around a slender body.

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 that is 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 & Rubinow (J. Fluid Mech., vol. 75, 1976, p. 705) and Johnson (J. Fluid Mech., vol. 99, 1979, p. 411) 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 the 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., 7 pages, 4 figures, submitted. LK and ML contributed equally

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 versus rotational propulsion.

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., This work was funded in part by the European Union through a CIG grant and a ERC Consolidator grant to EL and by the Cambridge Trust.

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 Mobius 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 as 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.