skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Reciprocal swimming at intermediate Reynolds number
In Stokes flow, Purcell's scallop theorem forbids objects with time-reversible (reciprocal) swimming strokes from moving. In the presence of inertia, this restriction is eased and reciprocally deforming bodies can swim. A number of recent works have investigated dimer models that swim reciprocally at intermediate Reynolds numbers $${\textit Re} \approx 1$$ –1000. These show interesting results (e.g. switches of the swim direction as a function of inertia) but the results vary and seem to be case specific. Here, we introduce a general model and investigate the behaviour of an asymmetric spherical dimer of oscillating length for small-amplitude motion at intermediate $${\textit {Re}}$$ . In our analysis we make the important distinction between particle and fluid inertia, both of which need to be considered separately. We asymptotically expand the Navier–Stokes equations in the small-amplitude limit to obtain a system of linear partial differential equations. Using a combination of numerical (finite element) and analytical (reciprocal theorem, method of reflections) methods we solve the system to obtain the dimer's swim speed and show that there are two mechanisms that give rise to motion: boundary conditions (an effective slip velocity) and Reynolds stresses. Each mechanism is driven by two classes of sphere–sphere interactions, between one sphere's motion and (1) the oscillating background flow induced by the other's motion, and (2) a geometric asymmetry induced by the other's presence. We can thus unify and explain behaviours observed in other works. Our results show how sensitive, counterintuitive and rich motility is in the parameter space of finite inertia of particles and fluid.  more » « less
Award ID(s):
1753148
PAR ID:
10412701
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
952
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The hydrodynamics of swimming at the millimeter-to-centimeter scale often present the challenge of having both viscous and inertial effects playing nontrivial roles. Inertial forces arise from the momentum of a moving fluid, while viscous forces come from friction within the flow. The non-dimensional Reynolds number (Re) compares the magnitudes of the inertial and viscous forces within a flow. At low Re (≪ 1), viscous forces dominate; at higher Re (≫ 1), inertial forces are more important. Efforts to understand the hydrodynamics of swimming have mainly focused on the extremes of fully viscous-dominated (Re ≪ 1) or inertia-dominated flow (Re ≫ 1). However, many animals swim in an intermediate regime, where inertia and viscosity are both significant. As an impactful and generalizable case study, we focus on ctenophores (comb jellies), a type of marine zooplankton. Ctenophores swim via the coordinated rowing of numerous highly flexible appendages (ctenes), with Reynolds numbers on the order of 10-100. Their locomotory dynamics present a unique opportunity to study the scaling of rowing (drag-based propulsion) across the low to intermediate Reynolds number range. With a combination of animal experiments, reduced-order analytical modeling, and physical-robotic modeling, we investigate how the kinematic and geometric variables of beating ctenes vary across Re, and how they affect swimming (including force production, speed, and maneuverability). Using animal experiments, we quantify the spatiotemporal asymmetry of beating ctenes across a wide range of animal sizes and Re. With our reduced-order model—the first to incorporate adequate formulations for the viscous-inertial nature of this regime—we explore the maneuverability and agility displayed by ctenophores, and show that by controlling the kinematics of their distributed appendages, ctenophores are capable of nearly omnidirectional swimming. Finally, we use a compliant robotic model that mimics ctenophore rowing kinematics to study rowing performance with direct calculation of thrust and lift forces distributed along the propulsor. These experiments shed new light on the relationship between motion asymmetries and thrust and lift production. This combination of animal experiments, analytical modeling, and physical modeling is the most detailed study of low to intermediate Re rowing to date, and provides a foundation for future applications in bio-inspired design. 
    more » « less
  2. Abstract Many organisms use flexible appendages for locomotion, feeding, and other functional behaviors. The efficacy of these behaviors is determined in large part by the fluid dynamics of the appendage interacting with its environment. For oscillating appendages at low Reynolds numbers, viscosity dominates over inertia, and appendage motion must be spatially asymmetric to generate net flow. At high Reynolds numbers, viscous forces are negligible and appendage motion is often also temporally asymmetric, with a fast power stroke and a slow recovery stroke; such temporal asymmetry does not affect the produced flow at low Reynolds numbers. At intermediate Reynolds numbers, both viscous and inertial forces play non-trivial roles—correspondingly, both spatial and temporal asymmetry can strongly affect overall propulsion. Here we perform experiments on three robotic paddles with different material flexibilities and geometries, allowing us to explore the effects of motion asymmetry (both spatial and temporal) on force production. We show how a flexible paddle’s time-varying shape throughout the beat cycle can reorient the direction of the produced force, generating both thrust and lift. We also evaluate the propulsive performance of the paddle by introducing a new quantity, which we term ‘integrated efficiency’. This new definition of propulsive efficiency can be used to directly evaluate an appendage’s performance independently from full-body swimming dynamics. Use of the integrated efficiency allows for accurate performance assessment, generalization, and comparison of oscillating appendages in both robotic devices and behaving organisms. Finally, we show that a curved flexible paddle generates thrust more efficiently than a straight paddle, and produces spatially asymmetric motion—thereby improving performance—without the need for complex actuation and controls, opening new avenues for bioinspired technology development. 
    more » « less
  3. This paper addresses peristaltic flow induced in a non-axisymmetric annular tube by a periodic small-amplitude wave of arbitrary shape propagating axially along its inner surface, assumed to be a circular cylinder. The study is motivated by recent in vivo experimental observations pertaining to the flow of cerebrospinal fluid along the perivascular spaces of cerebral arteries. The analysis employs the lubrication approximation, describing low-Reynolds-number peristaltic flow in the long-wavelength approximation. Closed-form analytic expressions are derived for the average pumping rate in infinitely long tubes and also in tubes of finite length. Consideration is also given to the transverse motion arising in non-axisymmetric tubes. For small-amplitude waves, the solution is reduced to the integration of a parameter-free Stokes-flow problem, which is solved for relevant cross-sectional shapes, with closed-form analytical results derived for thin canals. 
    more » « less
  4. Abstract Assemblies of self-rotating particles are gaining interest as a novel realization of active matter with unique collective behaviors such as edge currents and non-trivial dynamic states. Here, we develop a continuum model for a system of fluid-embedded spinners by coarse-graining the equations of motion of the discrete particles. We apply the model to explore mixtures of clockwise and counterclockwise rotating spinners. We find that the dynamics is sensitive to fluid inertia; in the inertialess system, after transient turbulent-like motion the spinners segregate and form steady traffic lanes. At small but finite Reynolds number instead, the turbulent-like motion persists and the system exhibits a chirality breaking transition leading to a single rotation sense state. Our results shed light on the dynamic behavior of non-equilibrium materials exemplified by active spinners. 
    more » « less
  5. The motion of a disk in a Langmuir film bounding a liquid substrate is a classical hydrodynamic problem, dating back to Saffman (J. Fluid Mech., vol. 73, 1976, p. 593) who focused upon the singular problem of translation at large Boussinesq number,$${\textit {Bq}}\gg 1$$. A semianalytic solution of the dual integral equations governing the flow at arbitrary$${\textit {Bq}}$$was devised by Hugheset al.(J. Fluid Mech., vol. 110, 1981, p. 349). When degenerated to the inviscid-film limit$${\textit {Bq}}\to 0$$, it produces the value$$8$$for the dimensionless translational drag, which is$$50\,\%$$larger than the classical$$16/3$$-value corresponding to a free surface. While that enhancement has been attributed to surface incompressibility, the mathematical reasoning underlying the anomaly has never been fully elucidated. Here we address the inviscid limit$${\textit {Bq}}\to 0$$from the outset, revealing a singular mechanism where half of the drag is contributed by the surface pressure. We proceed beyond that limit, considering a nearly inviscid film. A naïve attempt to calculate the drag correction using the reciprocal theorem fails due to an edge singularity of the leading-order flow. We identify the formation of a boundary layer about the edge of the disk, where the flow is primarily in the azimuthal direction with surface and substrate stresses being asymptotically comparable. Utilising the reciprocal theorem in a fluid domain tailored to the asymptotic topology of the problem produces the drag correction$$(8\,{\textit {Bq}}/{\rm \pi} ) [ \ln (2/{\textit {Bq}}) + \gamma _E+1]$$,$$\gamma _E$$being the Euler–Mascheroni constant. 
    more » « less