pith. sign in

arxiv: 2604.09496 · v1 · submitted 2026-04-10 · 🧮 math.AP · physics.flu-dyn

A nonlocal curve evolution for an immersed elastic filament: global existence and convergence to resistive force theory

Laurel Ohm This is my paper

Pith reviewed 2026-05-10 17:18 UTC · model grok-4.3

classification 🧮 math.AP physics.flu-dyn
keywords nonlocal curve evolutionelastic filamentStokes flowglobal well-posednessresistive force theorypseudodifferential operatorinextensible curve
0
0 comments X

The pith

Global well-posedness holds for a nonlocal curve evolution of an elastic filament in Stokes flow, with convergence to resistive force theory as radius shrinks to zero.

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

The paper examines a nonlocal model for the motion of an inextensible elastic filament immersed in a three-dimensional viscous fluid, where the fluid interaction is encoded by a pseudodifferential force-to-velocity operator. It establishes that solutions to this curve evolution exist for all time in the natural energy space. The analysis further demonstrates that these solutions converge to the simpler resistive force theory dynamics in the limit as the filament cross-sectional radius approaches zero. A sympathetic reader would care because the result supplies a rigorous foundation for resistive force theory, a widely used approximation in studies of microscopic biological filaments, and indicates that more complete three-dimensional free-boundary problems may inherit global existence when large-scale geometry is suitably controlled.

Core claim

We show global well-posedness for the curve evolution in the natural energy space. This loosely suggests that the full evolution may be globally well-posed if the large-scale geometry is controlled. Furthermore, we prove convergence to resistive force theory dynamics as ε→0, which illustrates how resistive force theory emerges from more detailed models.

What carries the argument

The pseudodifferential force-to-velocity operator that encodes fluid effects on the filament centerline and interpolates between resistive force theory at low wavenumbers and a Stokes boundary-value problem at high wavenumbers.

If this is right

  • Solutions to the nonlocal curve evolution exist globally in time when posed in the natural energy space.
  • The filament dynamics converge to those of resistive force theory in the limit as the cross-sectional radius ε tends to zero.
  • The convergence illustrates the emergence of resistive force theory from a more detailed fluid-structure interaction model.
  • The global well-posedness result loosely indicates that the corresponding full three-dimensional free-boundary problem remains globally well-posed provided large-scale geometry is controlled.

Where Pith is reading between the lines

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

  • Quantitative error estimates between the nonlocal model and resistive force theory could be derived for any fixed positive ε using the established convergence.
  • Numerical methods built on the nonlocal operator may inherit long-time stability from the global existence theorem.
  • The operator's interpolation between low- and high-wavenumber regimes may extend to other immersed slender-body problems such as flexible fibers or cilia.
  • Establishing control on large-scale geometry would allow the same energy methods to treat filament configurations with nontrivial topology.

Load-bearing premise

The model with the given pseudodifferential operator exactly captures the principal part of the full free-boundary problem for a filament of small but positive constant radius in Stokes flow.

What would settle it

A concrete smooth initial curve for which the solution develops a singularity in finite time, or a direct numerical comparison in which the ε→0 limit fails to recover resistive-force-theory velocities for a chosen test filament, would falsify the global existence or convergence statements.

Figures

Figures reproduced from arXiv: 2604.09496 by Laurel Ohm.

Figure 1
Figure 1. Figure 1: We consider X(s, t) as the centerline of an immersed 3D filament with constant cross sectional radius 0 < ϵ ≪ 1. The simplest choice of force-to-velocity map which still incorporates meaningful hydrodynamic effects is the resistive force theory (local slender body theory) [12, 16, 18, 37] approximation Lϵ,RFT(X) = |log ϵ| 4π (I + Xs ⊗ Xs). (1.4) The map Lϵ,RFT includes only the leading order O(|log ϵ|) eff… view at source ↗
Figure 2
Figure 2. Figure 2: Plots of the inverse tangential (left) and normal (right) direction mul￾tipliers mt ϵ (k) −1 and mn ϵ (k) −1 versus 2πϵ|k|. Here the inverses are used to more clearly display the linear (∼ ϵ|k|) behavior at high wavenumbers |k| ≳ 1 ϵ . For k ̸= 0, the multipliers mt ϵ (k) and mn ϵ (k) are precisely the tangential and normal eigenvalues, respectively, of the slender body Neumann-to-Dirichlet map in the spec… view at source ↗
read the original abstract

We consider a nonlocal curve evolution belonging to a hierarchy of models for the dynamics of an inextensible elastic filament in a 3D Stokes fluid. This model captures the principal part of a full free boundary problem for an elastic filament in Stokes flow. The fluid effects on the filament evolution are encoded in a pseudodifferential force-to-velocity operator which may be regarded as an interpolation between resistive force theory at low wavenumbers and a Stokes boundary value problem at high wavenumbers. Here the curve is considered to be the centerline of a 3D filament with constant cross sectional radius $\epsilon>0$. We show global well-posedness for the curve evolution in the natural energy space. This loosely suggests that the full evolution may be globally well-posed if the large-scale geometry is controlled. Furthermore, we prove convergence to resistive force theory dynamics as $\epsilon\to 0$, which illustrates how resistive force theory emerges from more detailed 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 paper studies a nonlocal inextensible curve evolution equation modeling the centerline of a thin elastic filament of radius ε immersed in 3D Stokes flow. Fluid effects are incorporated via a pseudodifferential force-to-velocity operator that interpolates between resistive-force theory at low wavenumbers and a full Stokes boundary-value problem at high wavenumbers. The central results are global well-posedness of the evolution in the natural energy space and convergence of the dynamics to resistive-force theory as ε→0.

Significance. If the proofs hold, the global existence result supplies a rigorous existence theory for a nonlocal regularization of filament dynamics that remains controlled at high frequencies, while the convergence theorem provides a precise mathematical link showing how resistive-force theory arises as the leading-order approximation to a more detailed nonlocal model. These are load-bearing contributions for the analysis of free-boundary fluid-structure problems involving slender bodies.

major comments (2)
  1. [Main existence theorem / energy estimates section] The global well-posedness statement (presumably Theorem 1.1 or the main existence theorem) relies on a priori energy estimates that close in the natural Sobolev space; the precise control of the nonlocal operator at high wavenumbers and the role of the inextensibility constraint in preventing filament self-intersection must be verified in detail, as any gap here would undermine the claim that the result “loosely suggests” global well-posedness for the full free-boundary problem.
  2. [Convergence theorem / asymptotic analysis section] The convergence proof as ε→0 (presumably the second main theorem) must specify the precise function-space topology and the rate (if any) at which the nonlocal velocity converges to the resistive-force-theory velocity; without an explicit error estimate or compactness argument, it is unclear whether the limit passage is strong enough to pass to the limit inside the nonlinear inextensibility constraint.
minor comments (2)
  1. [Abstract / Introduction] The abstract’s phrase “loosely suggests” is imprecise; replace it in the introduction with a concrete statement of what geometric control on the large-scale shape would be needed to extend the result to the full free-boundary problem.
  2. [Model formulation section] Notation for the pseudodifferential operator and the precise symbol of the force-to-velocity map should be introduced once and used consistently; several minor inconsistencies in the definition of the interpolation between low- and high-wavenumber regimes appear in the preliminary sections.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major comment in turn below, providing clarifications on the structure of the proofs and indicating the revisions we will make.

read point-by-point responses

  1. Referee: [Main existence theorem / energy estimates section] The global well-posedness statement (presumably Theorem 1.1 or the main existence theorem) relies on a priori energy estimates that close in the natural Sobolev space; the precise control of the nonlocal operator at high wavenumbers and the role of the inextensibility constraint in preventing filament self-intersection must be verified in detail, as any gap here would undermine the claim that the result “loosely suggests” global well-posedness for the full free-boundary problem.

    Authors: The a priori estimates are derived in Section 3 by multiplying the evolution equation by the velocity field and integrating by parts, exploiting the fact that the pseudodifferential force-to-velocity operator is self-adjoint and positive semi-definite. Its symbol, analyzed in Appendix A, coincides with the resistive-force-theory symbol at low frequencies and transitions to the full Stokes symbol (with the appropriate decay) at high frequencies; this guarantees the high-wavenumber dissipation needed to close the H^1 estimate without additional regularity. The inextensibility constraint is enforced pointwise by a Lagrange multiplier that projects the velocity onto the tangent bundle; because the evolution preserves arc-length parametrization exactly (see Lemma 2.3), the curve remains embedded in the energy space for all time. The phrase “loosely suggests” in the abstract is deliberately cautious and refers only to the fact that the same energy structure would be available in a full free-boundary setting once the large-scale geometry is controlled; we will add a short paragraph in the introduction making this distinction explicit. No changes to the proofs themselves are required, but we will expand the discussion of the symbol estimates and the preservation of embedding. revision: partial

  2. Referee: [Convergence theorem / asymptotic analysis section] The convergence proof as ε→0 (presumably the second main theorem) must specify the precise function-space topology and the rate (if any) at which the nonlocal velocity converges to the resistive-force-theory velocity; without an explicit error estimate or compactness argument, it is unclear whether the limit passage is strong enough to pass to the limit inside the nonlinear inextensibility constraint.

    Authors: Theorem 4.1 states that the solutions converge strongly in C([0,T]; H^1) as ε→0. The proof proceeds by obtaining uniform bounds from the energy estimates, applying the Aubin-Lions lemma to extract a strongly convergent subsequence in L^2, and then identifying the limit by testing against smooth test functions. The difference between the nonlocal velocity operator and its resistive-force-theory counterpart is controlled by an explicit O(ε) estimate in L^2 that follows from the symbol expansion (Lemma A.2). Because the inextensibility constraint is a linear pointwise condition on the velocity (vanishing normal component), strong L^2 convergence of the velocities is sufficient to pass to the limit inside the constraint. We will insert a dedicated paragraph after the statement of Theorem 4.1 that records the precise topology, the compactness argument, and the O(ε) rate, together with a short proof that the limit satisfies the inextensible resistive-force-theory equation. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper presents a mathematical analysis of a nonlocal curve evolution PDE model for an elastic filament in Stokes flow. The central results are a global well-posedness theorem in the natural energy space and a convergence result to resistive force theory as the filament radius ε tends to zero. These are established via standard techniques for nonlocal PDEs, including energy estimates, a priori bounds, and compactness arguments that derive from the structure of the pseudodifferential operator and the inextensibility constraint. No parameters are fitted to data, no quantities are defined in terms of the quantities they are used to predict, and no load-bearing steps reduce to self-citations or ansatzes imported from prior work by the same author. The model assumptions (curve as centerline, interpolation between RFT and Stokes) are stated explicitly as modeling choices rather than derived results, leaving the existence and convergence proofs self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard PDE assumptions for curve evolutions in Stokes flow; no free parameters, invented entities, or ad-hoc constants are visible from the abstract.

axioms (2)
  • domain assumption The curve represents the centerline of a filament with constant radius ε > 0
    Explicitly stated as the geometric setting for the model.
  • domain assumption Fluid effects are captured by a pseudodifferential force-to-velocity operator interpolating resistive force theory and Stokes flow
    This operator defines the nonlocal evolution equation.

pith-pipeline@v0.9.0 · 5457 in / 1375 out tokens · 99476 ms · 2026-05-10T17:18:17.364629+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

43 extracted references · 43 canonical work pages · 1 internal anchor

  1. [1]

    Albritton and L

    D. Albritton and L. Ohm. Rods in flows: the PDE theory of immersed elastic filaments.arXiv preprint arXiv:2503.14440, 2025

    work page arXiv 2025

  2. [2]

    Brezis and P

    H. Brezis and P. Mironescu. Gagliardo–nirenberg inequalities and non-inequalities: The full story.Annales de l’Institut Henri Poincar´ e C, Analyse non lin´ eaire, 35(5):1355–1376, 2018

    work page 2018

  3. [3]

    Camalet and F

    S. Camalet and F. J¨ ulicher. Generic aspects of axonemal beating.New Journal of Physics, 2(1):24, 2000

    work page 2000

  4. [4]

    Camalet, F

    S. Camalet, F. J¨ ulicher, and J. Prost. Self-organized beating and swimming of internally driven filaments.Physical review letters, 82(7):1590, 1999

    work page 1999

  5. [5]

    El Alaoui-Faris, J.-B

    Y. El Alaoui-Faris, J.-B. Pomet, S. R´ egnier, and L. Giraldi. Optimal actuation of flagellar magnetic microswim- mers.Physical Review E, 101(4):042604, 2020

    work page 2020

  6. [6]

    W. Fenchel. On the differential geometry of closed space curves.Bulletin of the AMS, 57, 1951

    work page 1951

  7. [7]

    B. M. Friedrich, I. H. Riedel-Kruse, J. Howard, and F. J¨ ulicher. High-precision tracking of sperm swimming fine structure provides strong test of resistive force theory.Journal of Experimental Biology, 213(8):1226–1234, 2010

    work page 2010

  8. [8]

    Gadˆ elha, E

    H. Gadˆ elha, E. Gaffney, D. Smith, and J. Kirkman-Brown. Nonlinear instability in flagellar dynamics: a novel modulation mechanism in sperm migration?Journal of The Royal Society Interface, 7(53):1689–1697, 2010

    work page 2010

  9. [9]

    Gadˆ elha and E

    H. Gadˆ elha and E. A. Gaffney. Flagellar ultrastructure suppresses buckling instabilities and enables mammalian sperm navigation in high-viscosity media.Journal of The Royal Society Interface, 16(152):20180668, 2019

    work page 2019

  10. [10]

    Garc´ ıa-Ju´ arez, P.-C

    E. Garc´ ıa-Ju´ arez, P.-C. Kuo, and Y. Mori. The immersed inextensible interface problem in 2d stokes flow.SIAM Journal on Mathematical Analysis, 57(4):3454–3487, 2025

    work page 2025

  11. [11]

    G¨ otz.Interactions of fibers and flow: asymptotics, theory and numerics

    T. G¨ otz.Interactions of fibers and flow: asymptotics, theory and numerics. Doctoral dissertation, University of Kaiserslautern, 2000

    work page 2000

  12. [12]

    Gray and G

    J. Gray and G. Hancock. The propulsion of sea-urchin spermatozoa.J. Exp. Biol., 32(4):802–814, 1955. 36 LAUREL OHM

    work page 1955

  13. [13]

    A. L. Hall-McNair, T. D. Montenegro-Johnson, H. Gadˆ elha, D. J. Smith, and M. T. Gallagher. Efficient imple- mentation of elastohydrodynamics via integral operators.Physical Review Fluids, 4(11):113101, 2019

    work page 2019

  14. [14]

    Hines and J

    M. Hines and J. Blum. Bend propagation in flagella. i. derivation of equations of motion and their simulation. Biophysical Journal, 23(1):41–57, 1978

    work page 1978

  15. [15]

    S. Hu, J. Zhang, and M. J. Shelley. Enhanced clamshell swimming with asymmetric beating at low reynolds number.Soft Matter, 2022

    work page 2022

  16. [16]

    Johnson and C

    R. Johnson and C. Brokaw. Flagellar hydrodynamics. a comparison between resistive-force theory and slender- body theory.Biophysical journal, 25(1):113–127, 1979

    work page 1979

  17. [17]

    R. E. Johnson. An improved slender-body theory for Stokes flow.J. Fluid Mech., 99(02):411–431, 1980

    work page 1980

  18. [18]

    J. B. Keller and S. Rubinow. Swimming of flagellated microorganisms.Biophysical Journal, 16(2):151–170, 1976

    work page 1976

  19. [19]

    J. B. Keller and S. I. Rubinow. Slender-body theory for slow viscous flow.J. Fluid Mech., 75(4):705–714, 1976

    work page 1976

  20. [20]

    N. Koiso. On the motion of a curve towards elastica.Actes de la Table Ronde de G´ eom´ etrie Diff´ erentielle (Luminy, 1992), 1:403–436, 1996

    work page 1992

  21. [21]

    Kuo, M.-C

    P.-C. Kuo, M.-C. Lai, Y. Mori, and A. Rodenberg. The tension determination problem for an inextensible interface in 2d stokes flow.Research in the Mathematical Sciences, 10(4):46, 2023

    work page 2023

  22. [22]

    E. Lauga. Floppy swimming: Viscous locomotion of actuated elastica.Physical Review E, 75(4):041916, 2007

    work page 2007

  23. [23]

    Lauga and C

    E. Lauga and C. Eloy. Shape of optimal active flagella.Journal of Fluid Mechanics, 730, 2013

    work page 2013

  24. [24]

    Lauga and T

    E. Lauga and T. R. Powers. The hydrodynamics of swimming microorganisms.Reports on Progress in Physics, 72(9):096601, 2009

    work page 2009

  25. [25]

    Montenegro-Johnson, H

    T. Montenegro-Johnson, H. Gadelha, and D. J. Smith. Spermatozoa scattering by a microchannel feature: an elastohydrodynamic model.Royal Society Open Science, 2(3):140475, 2015

    work page 2015

  26. [26]

    Moreau, F

    C. Moreau, F. Alouges, A. Lefebvre-Lepot, and J. Levillain. The n-link model for slender rods in a viscous fluid: well-posedness and convergence to classical elastohydrodynamics equations.HAL Id: hal-04944051, 2025

    work page 2025

  27. [27]

    Moreau, L

    C. Moreau, L. Giraldi, and H. Gadˆ elha. The asymptotic coarse-graining formulation of slender-rods, bio-filaments and flagella.Journal of the Royal Society Interface, 15(144):20180235, 2018

    work page 2018

  28. [28]

    Mori and L

    Y. Mori and L. Ohm. Accuracy of slender body theory in approximating force exerted by thin fiber on viscous fluid.Studies in Applied Mathematics, 2021

    work page 2021

  29. [29]

    Mori and L

    Y. Mori and L. Ohm. Well-posedness and applications of classical elastohydrodynamics for a swimming filament. Nonlinearity, 36(3):1799, 2023

    work page 2023

  30. [30]

    Y. Mori, L. Ohm, and D. Spirn. Theoretical justification and error analysis for slender body theory.Communi- cations on Pure and Applied Mathematics, 73(6):1245–1314, 2020

    work page 2020

  31. [31]

    Y. Mori, L. Ohm, and D. Spirn. Theoretical justification and error analysis for slender body theory with free ends.Archive for Rational Mechanics and Analysis, 235(3):1905–1978, 2020

    work page 1905

  32. [32]

    D. ¨Oelz. On the curve straightening flow of inextensible, open, planar curves.SeMA Journal, 54(1):5–24, 2011

    work page 2011

  33. [33]

    L. Ohm. A free boundary problem for an immersed filament in 3D Stokes flow.arXiv preprint arXiv:2404.04737, 2024

    work page arXiv 2024

  34. [34]

    L. Ohm. Well-posedness of a viscoelastic resistive force theory and applications to swimming.Journal of Nonlinear Science, 34(5):82, 2024

    work page 2024

  35. [35]

    L. Ohm. On an angle-averaged Neumann-to-Dirichlet map for thin filaments.Archive for Rational Mechanics and Analysis, 249(1):8, 2025

    work page 2025

  36. [36]

    L. Ohm. The slender body free boundary problem.arXiv preprint arXiv:2509.16800, 2025

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  37. [37]

    Pironneau and D

    O. Pironneau and D. Katz. Optimal swimming of flagellated micro-organisms.Journal of Fluid Mechanics, 66(2):391–415, 1974

    work page 1974

  38. [38]

    M. J. Shelley and T. Ueda. The Stokesian hydrodynamics of flexing, stretching filaments.Physica D: Nonlinear Phenomena, 146(1-4):221–245, 2000

    work page 2000

  39. [39]

    S. E. Spagnolie and E. Lauga. The optimal elastic flagellum.Physics of Fluids, 22(3):455, 2010

    work page 2010

  40. [40]

    Tornberg and M

    A.-K. Tornberg and M. J. Shelley. Simulating the dynamics and interactions of flexible fibers in Stokes flows. Journal of Computational Physics, 196(1):8–40, 2004

    work page 2004

  41. [41]

    B. J. Walker, K. Ishimoto, and E. A. Gaffney. Efficient simulation of filament elastohydrodynamics in three dimensions.Physical Review Fluids, 5(12):123103, 2020

    work page 2020

  42. [42]

    C. H. Wiggins and R. E. Goldstein. Flexive and propulsive dynamics of elastica at low reynolds number.Physical Review Letters, 80(17):3879, 1998

    work page 1998

  43. [43]

    C. H. Wiggins, D. Riveline, A. Ott, and R. E. Goldstein. Trapping and wiggling: elastohydrodynamics of driven microfilaments.Biophysical Journal, 74(2):1043–1060, 1998. Department of Mathematics, University of Wisconsin - Madison, Madison, WI 53706 Email address:lohm2@wisc.edu

    work page 1998