pith. sign in

arxiv: 2606.10661 · v1 · pith:FRCPCLUNnew · submitted 2026-06-09 · ⚛️ physics.flu-dyn

Far-field approximations for multi-timescale microswimmers near a boundary

Pith reviewed 2026-06-27 11:58 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords microswimmersfar-field approximationhydrodynamic interactionsboundariesmultiscale analysisforce dipolehovering statesswimmer trajectories
0
0 comments X

The pith

Incorporating time dependence and higher-order singularities into far-field models of microswimmers near boundaries reveals hovering states absent from simpler force-dipole models.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper extends minimal force-dipole models by incorporating higher-order flow singularities and using multiscale analysis to account for rapid oscillations in swimmer shape. This allows the model to handle dynamics closer to boundaries where leading-order terms fail. The extended model increases the dimensionality of the parameter space and predicts three behaviors: crashing, escaping, and hovering. Hovering is a new behavior that only emerges when time-dependence is included in the higher-order models. Sympathetic readers would care because averaged models may overlook stable positions that real microswimmers could exploit near surfaces.

Core claim

We demonstrate that the inclusion of time-dependence into these higher order models significantly expands the reachable parameter space, in particular by increasing its dimensionality. In these extended dynamics, we observe three distinct behaviours: crashing, escaping and hovering. Notably, hovering states are absent from the dynamics predicted by the simplest models, but are observed in more complex models.

What carries the argument

Multiscale analysis incorporating time-dependent higher-order flow singularities to extend the force-dipole far-field model.

If this is right

  • The parameter space for swimmer dynamics gains an extra dimension from time dependence.
  • Three behaviors become possible: crashing into the boundary, escaping away from it, and hovering at constant distance.
  • Hovering requires the inclusion of both higher-order singularities and time dependence.
  • The model remains applicable closer to the boundary than the basic dipole version.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Real microswimmers with periodic propulsion might achieve stable hovering near walls in ways missed by averaged models.
  • The multiscale method could extend to other problems involving oscillating swimmers in confined fluid flows.
  • Observations of hovering in experiments would support using higher-order terms over simple dipoles.

Load-bearing premise

The far-field approximation stays valid near the boundary even when higher-order singularities are added, and the separation between fast oscillation and slow swimming timescales holds.

What would settle it

Numerical simulation or physical experiment of an oscillating microswimmer near a wall showing hovering at the distance predicted by the higher-order time-dependent model but not by the time-independent dipole model.

Figures

Figures reproduced from arXiv: 2606.10661 by Benjamin J. Walker, Mohit P. Dalwadi, Sara Drummond-Curtis.

Figure 1
Figure 1. Figure 1: FIG. 1: Schematic of a microswimmer above a fixed boundary. We denote the distance between the centre of the [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Hovering dynamics of a single swimmer modelled as a force dipole and source dipole near a stress-free [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Escape dynamics of a single swimmer modelled as a force dipole and source dipole near a stress-free surface. [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: The proportion of the main behaviours (hovering, escaping or crashing) for a single swimmer modelled as a [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Phase portraits in ( [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Hovering dynamics of a single swimmer modelled as a force dipole and quadrupole near a stress-free surface. [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Escape dynamics of a single swimmer modelled as a force dipole and quadrupole near a stress-free surface. [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: The proportion of the main behaviours (hovering, escaping or crashing) for a single swimmer modelled as a [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Phase portraits in ( [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: The proportion of the main behaviours (hovering, escaping or crashing) for a single swimmer modelled as a [PITH_FULL_IMAGE:figures/full_fig_p015_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11: The proportion of the main behaviours (hovering, escaping or crashing) for a single swimmer modelled as a [PITH_FULL_IMAGE:figures/full_fig_p015_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12: The effect of a changing effective swimming speed on the prevalent behaviours for a single swimmer [PITH_FULL_IMAGE:figures/full_fig_p017_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13: The effect of a changing effective swimming speed on the prevalent behaviours for a single swimmer [PITH_FULL_IMAGE:figures/full_fig_p019_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14: Dynamics of a single swimmer modelled as a combination of a force dipole, a source dipole, and a [PITH_FULL_IMAGE:figures/full_fig_p020_14.png] view at source ↗
read the original abstract

Hydrodynamic interactions with boundaries can significantly affect the trajectories of microscale swimmers. In simple swimmer models, a common assumption is that swimmer shape remains constant, essentially averaging over the rapid oscillations in geometry and associated fluid flows that often are the source of propulsion. Previous work in minimal force-dipole models has shown how the inclusion of time-dependent swimmer changes can lead to a fundamentally wider class of behaviours than for their classic (implicitly averaged) counterparts. However, since force dipole models correspond to the leading-order term in the far-field description of the swimmer-induced flow, they break down as the swimmer approaches a boundary and predictions can become qualitatively inaccurate. Here, we extend the minimal force-dipole model by incorporating higher order flow singularities, systematically accounting for rapid oscillations in shape and singularity strength through a multiscale analysis. We demonstrate that the inclusion of time-dependence into these higher order models significantly expands the reachable parameter space, in particular by increasing its dimensionality. In these extended dynamics, we observe three distinct behaviours: crashing, escaping and hovering. Notably, hovering states are absent from the dynamics predicted by the simplest models, but are observed in more complex models.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript extends minimal force-dipole models of microswimmers near a no-slip boundary by systematically incorporating higher-order singularities (quadrupole and beyond) whose strengths vary with the swimmer's rapid shape oscillations. A multiscale asymptotic analysis separates fast and slow timescales to derive reduced ODEs governing the slow evolution of height and orientation. The resulting dynamics exhibit three regimes—crashing, escaping, and hovering—with hovering states appearing only when the higher-order terms are retained, thereby expanding the reachable parameter space relative to leading-order models.

Significance. If the far-field construction remains quantitatively reliable, the work supplies a concrete, higher-dimensional reduced model that predicts qualitatively new stable behaviors (hovering) absent from dipole-only approximations. This could serve as a practical bridge between minimal models and full boundary-integral simulations for microswimmer navigation near surfaces.

major comments (2)
  1. [§2] §2 and the asymptotic construction: the reduced ODEs for height and orientation are derived under the assumption that the swimmer-boundary distance remains ≫ swimmer radius even after inclusion of higher singularities; no explicit error bound, validity criterion, or matched near-field asymptotics is supplied to justify that the hovering fixed points lie inside the regime where the far-field truncation is accurate.
  2. [Results section (hovering states)] The central claim that hovering states are absent from simplest models but present in higher-order ones rests on the quantitative accuracy of the multiscale reduction near the boundary; without a direct comparison to full Stokes solutions or an a-posteriori check of the neglected near-field contributions, it is unclear whether the new states survive in the true hydrodynamics.
minor comments (2)
  1. [§2] Notation for the time-dependent singularity strengths is introduced without a clear table or appendix listing the explicit expressions for each multipole coefficient as functions of the fast phase.
  2. [Figures] Figure captions for the phase portraits should explicitly state the parameter values used and whether the trajectories are obtained from the reduced ODEs or from direct simulation.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The comments highlight important points about the regime of validity of the far-field model. We respond to each major comment below and indicate planned revisions.

read point-by-point responses
  1. Referee: [§2] §2 and the asymptotic construction: the reduced ODEs for height and orientation are derived under the assumption that the swimmer-boundary distance remains ≫ swimmer radius even after inclusion of higher singularities; no explicit error bound, validity criterion, or matched near-field asymptotics is supplied to justify that the hovering fixed points lie inside the regime where the far-field truncation is accurate.

    Authors: We agree that the manuscript does not supply explicit a-posteriori error bounds or matched near-field asymptotics. The derivation in §2 relies on the standard far-field ordering h ≫ a, and the reported hovering equilibria satisfy this ordering for the parameter values examined. We will revise §2 to include a brief discussion of the expected truncation error (based on the relative magnitude of the retained quadrupole and higher terms) and a simple numerical check confirming that the hovering fixed points remain outside the near-field region for the plotted trajectories. revision: yes

  2. Referee: [Results section (hovering states)] The central claim that hovering states are absent from simplest models but present in higher-order ones rests on the quantitative accuracy of the multiscale reduction near the boundary; without a direct comparison to full Stokes solutions or an a-posteriori check of the neglected near-field contributions, it is unclear whether the new states survive in the true hydrodynamics.

    Authors: The manuscript presents a systematic far-field reduction; the central claim is therefore that, within this reduced model, inclusion of time-dependent higher-order singularities qualitatively enlarges the phase space to include stable hovering states absent from the dipole truncation. We do not assert that the hovering states are quantitatively accurate for the full Stokes problem at small h. We will add an explicit caveat in the Results section noting the far-field limitation and that direct validation against boundary-integral simulations lies outside the present scope. revision: partial

Circularity Check

0 steps flagged

No circularity: derivation uses standard asymptotic expansions independent of outputs

full rationale

The paper constructs its multiscale far-field model via systematic inclusion of higher-order singularities and time-dependent strengths, starting from the known breakdown of the leading dipole near boundaries. This yields extended ODEs whose solutions exhibit new behaviors such as hovering; these outcomes are not presupposed in the construction or obtained by fitting parameters to the target dynamics. Prior force-dipole literature is cited for context and differentiation rather than as a load-bearing uniqueness theorem or ansatz that forces the present results. No self-definitional steps, renamed empirical patterns, or predictions that reduce to fitted inputs appear in the derivation chain.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Based solely on the abstract, no specific free parameters, axioms, or invented entities can be identified; the approach relies on standard far-field singularity expansions and timescale separation common to the field.

pith-pipeline@v0.9.1-grok · 5741 in / 1238 out tokens · 34713 ms · 2026-06-27T11:58:06.300543+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

53 extracted references · 2 canonical work pages

  1. [1]

    (19a), independent of the value ofhor the system parameters

    Existence of trivial steady states We note dθ/dt= 0 forθ=nπ,n∈Zin Eq. (19a), independent of the value ofhor the system parameters. We therefore refer to the steady states withθ= 0, πastrivial, corresponding to the swimmer facing away from and towards the boundary, respectively. Without loss of generality, we only consider s <0 as we can map to s >0 by shi...

  2. [2]

    (19b) yields the steady height as a function ofθ h= 2scosθ p(1 + 3 cos 2θ).(22) 8 By substituting Eq

    Existence of non-trivial steady states To find the remaining steady states, which we refer to asnon-trivialsteady states, we note that seeking dh/dt= 0 in Eq. (19b) yields the steady height as a function ofθ h= 2scosθ p(1 + 3 cos 2θ).(22) 8 By substituting Eq. (22) into Eq. (19a) and factoring out 3 ssinθ/32h 3, we find that the non-trivialθsteady states ...

  3. [3]

    (27a) for all values ofhand system parameters, and continue to refer toθ= 0, πas trivial

    Existence of trivial steady states We again find that dθ/dt= 0 forθ=nπ,n∈Zin Eq. (27a) for all values ofhand system parameters, and continue to refer toθ= 0, πas trivial. We assume q >0 as we can map to q <0 by shifting byθbyπ. Substitutingθ= 0,±π into Eq. (27b) the trivial steady states are θ= 0, h= q p andθ=π, h=− q p .(28) We find that the eigenvalues ...

  4. [4]

    Existence of non-trivial steady states To find the remaining, non-trivial, steady states we proceed as before and seek solutions to dh/dt= 0 in Eq. (27b). The steady height as a function ofθat fixed points is h=− 2qcosθ p(1 + 3 cosθ) (1−3 cos 2θ).(30) We substitute Eq. (30) into Eq. (27a) and factor out 3qsinθ/32h 4 (the roots of which correspond to the t...

  5. [5]

    (19a) does not depend on the swimming speed, we have the same trivial steady statesθ= 0, θ=πas Section III B

    Trivial steady states As theθequation Eq. (19a) does not depend on the swimming speed, we have the same trivial steady statesθ= 0, θ=πas Section III B. Now, however, the corresponding heights are solutions to cubic polynomials: 4uh3 + 2ph− s= 0,(40a) −4uh 3 + 2ph+ s= 0,(40b) forθ= 0 andθ=π, respectively. Forh >0 in the domain, we find three cases: no solu...

  6. [6]

    However, theh equation Eq

    Non-trivial steady states For the remaining, non-trivial, steady states, we seek the steady values ofhas a function ofθ. However, theh equation Eq. (39) corresponds to solving the cubic equation 8ucosθh 3 + p(1 + 3 cos 2θ)h−2 scosθ= 0.(41) Instead of seeking solutions directly, we factor out 3 sinθ/32h 4 from theθequation Eq. (19a) and write h= Bss 2pcosθ...

  7. [7]

    Klapper and J

    I. Klapper and J. Dockery, Mathematical description of microbial biofilms, SIAM Rev.52, 221 (2010)

  8. [8]

    S. B. Guttenplan and D. B. Kearns, regulation of flagellar motility during biofilm formation, FEMS Microbiol. Rev.37, 849 (2013)

  9. [9]

    Suarez and A

    S. Suarez and A. A. Pacey, Sperm transport in the female reproductive tract, Hum. Reprod. Update12, 23 (2006)

  10. [10]

    Elgeti, U

    J. Elgeti, U. B. Kaupp, and G. Gompper, Hydrodynamics of sperm cells near surfaces, Biophys. J.99, 1018 (2010)

  11. [11]

    Rothschild, Non-random distribution of bull spermatozoa in a drop of sperm suspension, Nature198, 1221 (1963)

  12. [12]

    P. D. Frymier, R. M. Ford, H. C. Berg, and P. T. Cummings, Three-dimensional tracking of motile bacteria near a solid planar surface., PNAS92, 6195 (1995)

  13. [13]

    Bianchi, F

    S. Bianchi, F. Saglimbeni, G. Frangipane, D. Dell’Arciprete, and R. Di Leonardo, 3D dynamics of bacteria wall entrapment at a water–air interface, Soft Matter15, 3397 (2019)

  14. [14]

    Berg and L

    H. Berg and L. Turner, Chemotaxis of bacteria in glass capillary arrays. Escherichia coli, motility, microchannel plate, and light scattering, Biophys. J.58, 919 (1990)

  15. [15]

    A. P. Berke, L. Turner, H. C. Berg, and E. Lauga, Hydrodynamic attraction of swimming microorganisms by surfaces, Phys. Rev. Lett.101, 038102 (2008)

  16. [16]

    Di Leonardo, D

    R. Di Leonardo, D. Dell’Arciprete, L. Angelani, and V. Iebba, Swimming with an image, Phys. Rev. Lett.106, 038101 (2011)

  17. [17]

    Kantsler, J

    V. Kantsler, J. Dunkel, M. Polin, and R. E. Goldstein, Ciliary contact interactions dominate surface scattering of swimming eukaryotes, PNAS110, 1187 (2013)

  18. [18]

    Bianchi, F

    S. Bianchi, F. Saglimbeni, and R. Di Leonardo, Holographic imaging reveals the mechanism of wall entrapment in swimming bacteria, Phys. Rev. X7, 10.1103/physrevx.7.011010 (2017)

  19. [19]

    S. E. Spagnolie and E. Lauga, Hydrodynamics of self-propulsion near a boundary: predictions and accuracy of far-field approximations, J. Fluid Mech.700, 105 (2012)

  20. [20]

    Elgeti and G

    J. Elgeti and G. Gompper, Microswimmers near surfaces, Eur. Phys. J. Spec. Top.225, 2333 (2016)

  21. [21]

    E. F. Yeo, B. J. Walker, P. Pearce, and M. P. Dalwadi, A shear-induced limit on bacterial surface adhesion in fluid flow, PNAS123, e2516069123 (2026)

  22. [22]

    Lauga, W

    E. Lauga, W. R. DiLuzio, G. M. Whitesides, and H. A. Stone, Swimming in circles: motion of bacteria near solid boundaries, Biophys. J.90, 400 (2006)

  23. [23]

    Or and R

    Y. Or and R. M. Murray, Dynamics and stability of a class of low reynolds number swimmers near a wall, Phys. Rev. E: Stat. Nonlinear Soft Matter Phys.79, 045302 (2009)

  24. [24]

    H. Shum, E. A. Gaffney, and D. J. Smith, Modelling bacterial behaviour close to a no-slip plane boundary: the influence of bacterial geometry, Proc. R. Soc. A466, 1725 (2010)

  25. [25]

    Giacch´ e, T

    D. Giacch´ e, T. Ishikawa, and T. Yamaguchi, Hydrodynamic entrapment of bacteria swimming near a solid surface, Phys. Rev. E82, 056309 (2010)

  26. [26]

    Manabe, T

    J. Manabe, T. Omori, and T. Ishikawa, Shape matters: entrapment of a model ciliate at interfaces, J. Fluid Mech.892, A15 (2020)

  27. [27]

    Lauga,The fluid dynamics of cell motility, Cambridge Texts in Applied Mathematics (Cambridge University Press, 2020)

    E. Lauga,The fluid dynamics of cell motility, Cambridge Texts in Applied Mathematics (Cambridge University Press, 2020)

  28. [28]

    Drescher, R

    K. Drescher, R. E. Goldstein, N. Michel, M. Polin, and I. Tuval, Direct Measurement of the Flow Field around Swimming Microorganisms, Phys. Rev. Lett.105, 168101 (2010). 23

  29. [29]

    Drescher, J

    K. Drescher, J. Dunkel, L. H. Cisneros, S. Ganguly, and R. E. Goldstein, Fluid dynamics and noise in bacterial cell–cell and cell–surface scattering, PNAS108, 10940 (2011)

  30. [30]

    Ishimoto and E

    K. Ishimoto and E. A. Gaffney, Squirmer dynamics near a boundary, Phys. Rev. E88, 062702 (2013)

  31. [31]

    D. J. Smith, E. A. Gaffney, J. R. Blake, and J. C. Kirkman-Brown, Human sperm accumulation near surfaces: a simulation study, J. Fluid Mech.621, 289 (2009)

  32. [32]

    P. H. Htet, D. Das, and E. Lauga, Hydrodynamic hovering of swimming bacteria above surfaces, Phys. Rev. Res.6, L032070 (2024)

  33. [33]

    B. J. Walker, K. Ishimoto, C. Moreau, E. A. Gaffney, and M. P. Dalwadi, Emergent rheotaxis of shape-changing swimmers in Poiseuille flow, J. Fluid Mech.944, R2 (2022)

  34. [34]

    Omori, K

    T. Omori, K. Kikuchi, M. Schmitz, M. Pavlovic, C.-H. Chuang, and T. Ishikawa, Rheotaxis and migration of an unsteady microswimmer, J. Fluid Mech.930, A30 (2022)

  35. [35]

    B. J. Walker, K. Ishimoto, E. A. Gaffney, C. Moreau, and M. P. Dalwadi, Effects of rapid yawing on simple swimmer models and planar Jeffery’s orbits, Phys. Rev. Fluids7, 023101 (2022)

  36. [36]

    E. A. Gaffney, M. P. Dalwadi, C. Moreau, K. Ishimoto, and B. J. Walker, Canonical orbits for rapidly deforming planar microswimmers in shear flow, Phys. Rev. Fluids7, L022101 (2022)

  37. [37]

    M. P. Dalwadi, C. Moreau, E. A. Gaffney, K. Ishimoto, and B. J. Walker, Generalised Jeffery’s equations for rapidly spinning particles. Part 1. Spheroids, J. Fluid Mech.979, A1 (2024)

  38. [38]

    M. P. Dalwadi, C. Moreau, E. A. Gaffney, B. J. Walker, and K. Ishimoto, Generalised Jeffery’s equations for rapidly spinning particles. Part 2. Helicoidal objects with chirality, J. Fluid Mech.979, A2 (2024)

  39. [39]

    M. P. Dalwadi, Rapidly yawing spheroids in viscous shear flow: emergent loss of symmetry, J. Fluid Mech.1009, A27 (2025)

  40. [40]

    B. J. Walker, K. Ishimoto, and E. A. Gaffney, Systematic parameterizations of minimal models of microswimming, Phys. Rev. Fluids8, 034102 (2023)

  41. [41]

    M. H. Holmes,Introduction to perturbation methods, Vol. 20 (Springer Science & Business Media, 2012)

  42. [42]

    C. M. Bender and S. A. Orszag,Advanced mathematical methods for scientists and engineers I(New York, NY, 1999)

  43. [43]

    A. T. Chwang and T. Y.-T. Wu, Hydromechanics of low-Reynolds-number flow. Part 2. Singularity method for Stokes flows, J. Fluid Mech.67, 787 (1975)

  44. [44]

    Kim and S

    S. Kim and S. J. Karrila,Microhydrodynamics: principles and selected applications(Butterworth-Heinemann, 2013)

  45. [45]

    F. P. Bretherton, The motion of rigid particles in a shear flow at low Reynolds number, J. Fluid Mech.14, 284 (1962)

  46. [46]

    G. B. Jeffery, The motion of ellipsoidal particles immersed in a viscous fluid, Proc. R. Soc. Lond. A102, 161 (1922)

  47. [47]

    J. S. Guasto, K. A. Johnson, and J. P. Gollub, Oscillatory flows induced by microorganisms swimming in two dimensions, Phys. Rev. Lett.105, 168102 (2010)

  48. [48]

    Ishimoto, H

    K. Ishimoto, H. Gadˆ elha, E. A. Gaffney, D. J. Smith, and J. Kirkman-Brown, Coarse-graining the fluid flow around a human sperm, Phys. Rev. Lett.118, 124501 (2017)

  49. [49]

    S. F. Schoeller and E. E. Keaveny, From flagellar undulations to collective motion: predicting the dynamics of sperm suspensions, Journal of The Royal Society Interface15, 20170834 (2018)

  50. [50]

    L. F. Shampine and M. W. Reichelt, The matlab ode suite (1997), https://doi.org/10.1137/S1064827594276424

  51. [51]

    Golestanian, J

    R. Golestanian, J. M. Yeomans, and N. Uchida, Hydrodynamic synchronization at low Reynolds number, Soft Matter7, 3074 (2011)

  52. [52]

    Elgeti and G

    J. Elgeti and G. Gompper, Emergence of metachronal waves in cilia arrays, PNAS110, 4470 (2013)

  53. [53]

    Friedrich, Hydrodynamic synchronization of flagellar oscillators, Eur

    B. Friedrich, Hydrodynamic synchronization of flagellar oscillators, Eur. Phys. J. Spec. Top.225, 2353 (2016)