Saverio E. Spagnolie

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. A general axisymmetric swimmer is described as a linear combination of fundamental solutions to the Stokes equations: a Stokeslet dipole, a source dipole, a Stokeslet quadrupole, and a rotlet dipole. The effects of nearby walls or stress-free surfaces on swimming trajectories are described through the contribution of each singularity, and we address the question of how accurately this multipole approach captures the wall effects observed in full numerical solutions of the Stokes equations. The reduced model is used to provide simple but accurate predictions of the wall-induced attraction and pitching dynamics for model Janus particles, ciliated organisms, and bacteria-like polar swimmers. Transitions in attraction and pitching behaviour as functions of body geometry and propulsive mechanism are described. The reduced model may help to explain a number of recent experimental results.

Surprising behaviors in flapping locomotion with passive pitching

To better understand the role of wing and fin flexibility in flapping locomotion, we study through experiment and numerical simulation a freely moving wing that can “pitch” passively as it is actively heaved in a fluid. We observe a range of flapping frequencies corresponding to large horizontal velocities, a regime of underperformance relative to a clamped (nonpitching) flapping wing, and a surprising, hysteretic regime in which the flapping wing can move horizontally in either direction (despite left/right symmetry being broken by the specific mode of pitching). The horizontal velocity is shown to peak when the flapping frequency is near the immersed system’s resonant frequency. Unlike for the clamped wing, we find that locomotion is achieved by vertically flapped symmetric wings with even the slightest pitching flexibility, and the system exhibits a continuous departure from the Stokesian regime. The phase difference between the vertical heaving motion and consequent pitching changes continuously with ...

Locomotion of Helical Bodies in Viscoelastic Fluids: Enhanced Swimming at Large Helical Amplitudes

The motion of a rotating helical body in a viscoelastic fluid is considered. In the case of force-free swimming, the introduction of viscoelasticity can either enhance or retard the swimming speed and locomotive efficiency, depending on the body geometry, fluid properties, and the body rotation rate. Numerical solutions of the Oldroyd-B equations show how previous theoretical predictions break down with increasing helical radius or with decreasing filament thickness. Helices of large pitch angle show an increase in swimming speed to a local maximum at a Deborah number of order unity. The numerical results show how the small-amplitude theoretical calculations connect smoothly to the large-amplitude experimental measurements.

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 (nonsmooth). Physically, biological organisms (or engineered microswimmers) 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 for a given 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 i...

Comparative hydrodynamics of bacterial polymorphism.

Most bacteria swim through fluids by rotating helical flagella which can take one of 12 distinct polymorphic shapes, the most common of which 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 efficient of the 12 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. Hence, although energetic costs of locomotion are small for bacteria, fluid mechanical forces may have played a significant role in the evolution of the flagellum.

Shape-changing bodies in fluid: Hovering, ratcheting, and bursting

Motivated by recent experiments on the hovering of passive bodies, we demonstrate how a simple shape-changing body can hover or ascend in an oscillating background flow. We study this ratcheting effect through numerical simulations of the two-dimensional Navier–Stokes equations at intermediate Reynolds number. This effect could describe a viable means of locomotion or transport in such environments as a tidal pool with wave-driven sloshing. We also consider the velocity burst achieved by a body through a rapid increase in its aspect ratio, which may contribute to the escape dynamics of such organisms as jellyfish.

The sedimentation of flexible filaments

Author(s): Li, L; Manikantan, H; Saintillan, D; Spagnolie, SE | Abstract: The dynamics of a flexible filament sedimenting in a viscous fluid are explored analytically and numerically. Compared with the well-studied case of sedimenting rigid rods, the introduction of filament compliance is shown to cause a significant alteration in the long-time sedimentation orientation and filament geometry. A model is developed by balancing viscous, elastic and gravitational forces in a slender-body theory for zero-Reynolds-number flows, and the filament dynamics are characterized by a dimensionless elasto-gravitation number. Filaments of both non-uniform and uniform cross-sectional thickness are considered. In the weakly flexible regime, a multiple-scale asymptotic expansion is used to obtain expressions for filament translations, rotations and shapes. These are shown to match excellently with full numerical simulations. Furthermore, we show that trajectories of sedimenting flexible filaments, unlike their rigid counterparts, are restricted to a cloud whose envelope is determined by the elasto-gravitation number. In the highly flexible regime we show that a filament sedimenting along its long axis is susceptible to a buckling instability. A linear stability analysis provides a dispersion relation, illustrating clearly the competing effects of the compressive stress and the restoring elastic force in the buckling process. The instability travels as a wave along the filament opposite the direction of gravity as it grows and the predicted growth rates are shown to compare favourably with numerical simulations. The linear eigenmodes of the governing equation are also studied, which agree well with the finite-amplitude buckled shapes arising in simulations.

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, which corresponds to the potential flow created by a source dipole at the sphere center. Unlike other swimming mechanisms which rely on the presentation of a small cross section ...

A bug on a raft: recoil locomotion in a viscous fluid

The locomotion of a body through an inviscid incompressible fluid, such that the flow remains irrotational everywhere, is known to depend on inertial forces and on both the shape and the mass distribution of the body. In this paper we consider the influence of fluid viscosity on such inertial modes of locomotion. In particular we consider a free body of variable shape and study the centre-of-mass and centre-of-volume variations caused by a shifting mass distribution. We call this recoil locomotion . Numerical solutions of a finite body indicate that the mechanism is ineffective in Stokes flow but that viscosity can significantly increase the swimming speed above the inviscid value once Reynolds numbers are in the intermediate range 50–300. To study the problem analytically, a model which is an analogue of Taylors swimming sheet is introduced. The model admits analysis at fixed, arbitrarily large Reynolds number for deformations of sufficiently small amplitude. The analysis confirms the significant increase of swimming velocity above the inviscid value at intermediate Reynolds numbers.

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.