pith. sign in

arxiv: 2509.16800 · v1 · submitted 2025-09-20 · 🧮 math.AP · physics.flu-dyn

The slender body free boundary problem

Laurel Ohm This is my paper

Pith reviewed 2026-05-18 15:27 UTC · model grok-4.3

classification 🧮 math.AP physics.flu-dyn
keywords slender body theoryfree boundary problemStokes fluidelastic filamentinextensibility constraintNeumann-to-Dirichlet mapEuler-Bernoulli beamfluid-structure interaction
0
0 comments X

The pith

Existence of solutions is established for the evolution of an inextensible elastic filament in a Stokes fluid via the slender body Neumann-to-Dirichlet map.

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

The paper sets out to prove the existence of solutions to the evolution equation for a closed elastic filament that cannot stretch as it bends and moves through a surrounding viscous fluid. The filament is modeled as a one-dimensional curve whose elasticity follows Euler-Bernoulli theory, while the fluid interaction is captured by a special map that treats the filament as having a small but fixed thickness. This work matters because it supplies the rigorous backing for many practical simulations of thin biological or synthetic structures in liquids, such as swimming microorganisms or flexible polymers. The proof combines an analysis of the leading term in the fluid response map with a careful determination of the tension needed to keep the filament length fixed.

Core claim

We consider the slender body free boundary problem describing the evolution of an inextensible, closed elastic filament immersed in a Stokes fluid in R^3. The filament elasticity is governed by Euler-Bernoulli beam theory, and the coupling between this 1D elasticity law and the surrounding 3D fluid is governed by the slender body Neumann-to-Dirichlet map. We develop a solution theory for the filament evolution under this coupling. Our analysis relies on two main ingredients: an extraction of the principal symbol of the slender body NtD map and a detailed treatment of the tension determination problem for enforcing the inextensibility constraint.

What carries the argument

The principal symbol of the slender body Neumann-to-Dirichlet map, used to couple the one-dimensional Euler-Bernoulli elasticity to the three-dimensional fluid while the tension determination procedure enforces the inextensibility constraint.

If this is right

  • The model provides a mathematical foundation for computational simulations of slender filaments in fluids.
  • It justifies the use of reduced one-dimensional models for filaments with small but positive cross-sectional radius in dynamic settings.
  • The tension determination procedure ensures the filament length remains fixed throughout the evolution.
  • The solution theory supports analysis of closed inextensible structures evolving under combined bending and fluid forces.

Where Pith is reading between the lines

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

  • The existence result could be leveraged to study the long-time behavior or stability of particular filament configurations such as rings or helices.
  • Numerical methods that discretize the one-dimensional filament equation with the approximated map may inherit convergence guarantees from this theory.
  • The same combination of principal symbol extraction and tension solving might extend to filaments with slowly varying thickness or to interactions with boundaries.

Load-bearing premise

The principal symbol of the slender body Neumann-to-Dirichlet map extracted in prior work remains valid and sufficient for the current free-boundary evolution problem with the inextensibility constraint.

What would settle it

A concrete counterexample consisting of an initial closed filament shape and short time interval on which the reduced one-dimensional system with the extracted symbol produces dynamics that diverge from a full three-dimensional Stokes simulation as the filament radius tends to zero would show the symbol extraction is not sufficient.

Figures

Figures reproduced from arXiv: 2509.16800 by Laurel Ohm.

Figure 1
Figure 1. Figure 1: Snapshot of a filament Σϵ(t) with centerline X(s, t). To analyze (4)-(7), it will be useful for us to reformulate the slender body free boundary problem as the following curve evolution: ∂X ∂t = −Lϵ(X) [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
read the original abstract

We consider the slender body free boundary problem describing the evolution of an inextensible, closed elastic filament immersed in a Stokes fluid in $\mathbb{R}^3$. The filament elasticity is governed by Euler-Bernoulli beam theory, and the coupling between this 1D elasticity law and the surrounding 3D fluid is governed by the slender body Neumann-to-Dirichlet (NtD) map, which treats the filament as a 3D object with constant cross-sectional radius $0<\epsilon\ll1$. This map serves as a mathematical justification for slender body theories wherein such 3D-1D couplings play a central role. We develop a solution theory for the filament evolution under this coupling. Our analysis relies on two main ingredients: (1) an extraction of the principal symbol of the slender body NtD map, from the author's previous work, and (2) a detailed treatment of the tension determination problem for enforcing the inextensibility constraint. Our work provides a mathematical foundation for various computational models in which a slender filament evolves according to a 1D elasticity law in a 3D fluid. This forms a key development in our broader program to place slender body theories on firm theoretical footing.

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 develops a solution theory for the free-boundary evolution of a closed, inextensible elastic filament immersed in a 3D Stokes fluid. The filament is modeled by Euler-Bernoulli beam theory coupled to the fluid via the slender-body Neumann-to-Dirichlet map (with small cross-sectional radius ε), and inextensibility is enforced by determining a tension field. The analysis rests on two ingredients: extraction of the principal symbol of the NtD map from the author's prior work, and a separate treatment of the tension-determination problem.

Significance. If the central claims hold, the work supplies the first rigorous existence theory for an inextensible slender-body free-boundary problem. This directly supports the mathematical justification of slender-body approximations used in computational models of biological filaments and provides a foundation for further analysis of related 3D-1D fluid-structure systems.

major comments (2)
  1. [§2.2] §2.2 and the statement of the evolution equation: the principal symbol of the slender-body NtD map is imported unchanged from the author's previous paper, yet the present work adds the inextensibility constraint and closes the free-boundary coupling. No re-derivation or perturbation argument is supplied showing that the leading-order symbol survives these modifications; this step is load-bearing for all subsequent a-priori estimates.
  2. [§4] §4 (Tension determination problem): the tension is treated as a lower-order correction that enforces the inextensibility constraint, but the manuscript does not provide explicit mapping properties or symbol estimates demonstrating that the tension operator remains a compact perturbation of the principal part once the free-boundary evolution is closed. Without such control the fixed-point or iteration scheme used for local existence cannot be justified.
minor comments (2)
  1. [§3] The dependence of all constants on the small parameter ε is not tracked uniformly through the estimates; a remark clarifying the ε-independent nature of the final existence time would improve readability.
  2. [Notation] Notation for the filament centerline and its tangent vector is introduced in §1 but reused with slight variations in later sections; a single consolidated table of symbols would reduce ambiguity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the constructive comments. The referee correctly identifies two points where additional justification would strengthen the exposition. We address each comment below and will revise the manuscript accordingly to incorporate the requested clarifications.

read point-by-point responses

  1. Referee: [§2.2] §2.2 and the statement of the evolution equation: the principal symbol of the slender-body NtD map is imported unchanged from the author's previous paper, yet the present work adds the inextensibility constraint and closes the free-boundary coupling. No re-derivation or perturbation argument is supplied showing that the leading-order symbol survives these modifications; this step is load-bearing for all subsequent a-priori estimates.

    Authors: We agree that an explicit argument confirming the invariance of the leading-order symbol under the added inextensibility constraint would improve clarity. The tension enters the evolution equation only through a tangential projection that does not affect the highest-order terms in the slender-body NtD symbol; the free-boundary coupling is already encoded in the map derived in our prior work. Nevertheless, to make this transparent, we will insert a short perturbation calculation in the revised §2.2 that isolates the principal symbol and verifies that the inextensibility correction remains of strictly lower order. revision: yes

  2. Referee: [§4] §4 (Tension determination problem): the tension is treated as a lower-order correction that enforces the inextensibility constraint, but the manuscript does not provide explicit mapping properties or symbol estimates demonstrating that the tension operator remains a compact perturbation of the principal part once the free-boundary evolution is closed. Without such control the fixed-point or iteration scheme used for local existence cannot be justified.

    Authors: The tension is recovered by solving a scalar, nonlocal but elliptic equation obtained by enforcing the inextensibility constraint pointwise along the filament. This operator is of order at most zero and therefore compact relative to the principal part generated by bending and the leading term of the NtD map. We acknowledge, however, that the manuscript would benefit from a concise statement of the precise mapping properties and symbol estimates for the closed system. We will add these estimates, together with a brief justification that they suffice for the contraction mapping argument, in the revised §4. revision: yes

Circularity Check

1 steps flagged

Solution theory depends on principal symbol from author's previous work as a main ingredient

specific steps
  1. self citation load bearing [Abstract]
    "Our analysis relies on two main ingredients: (1) an extraction of the principal symbol of the slender body NtD map, from the author's previous work, and (2) a detailed treatment of the tension determination problem for enforcing the inextensibility constraint."

    The solution theory for filament evolution under the slender-body coupling is constructed to rely on the principal symbol extracted in the same author's earlier paper. While the tension treatment adds new analysis, the load-bearing use of the prior symbol extraction means the central result depends on self-cited work without re-derivation or re-justification shown inside the coupled free-boundary system.

full rationale

The paper states that its analysis relies on two main ingredients, one of which is the principal symbol extraction from the author's prior work. This is a self-citation that carries load for the central solution theory. However, the paper also contributes an independent detailed treatment of the tension determination problem to enforce inextensibility, so the overall claim retains substantial independent content and does not reduce entirely to the self-citation. No self-definitional, fitted-prediction, or other circular reductions are present in the given material.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of the slender-body Neumann-to-Dirichlet approximation for small but positive filament radius and on the solvability of the tension problem for the inextensibility constraint; both are treated as domain-standard assumptions rather than newly derived here.

axioms (2)
  • domain assumption The slender body Neumann-to-Dirichlet map approximation holds for 0 < epsilon << 1 and supplies the correct leading-order coupling between 1D elasticity and 3D Stokes flow.
    Invoked to justify treating the filament as a 3D object with constant cross-sectional radius while reducing the fluid problem to a 1D evolution equation.
  • domain assumption The inextensibility constraint can be enforced by solving a tension determination problem that remains well-posed under the extracted principal symbol.
    Central to the second main ingredient listed in the abstract.

pith-pipeline@v0.9.0 · 5731 in / 1506 out tokens · 49100 ms · 2026-05-18T15:27:38.603144+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

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

    math.AP 2026-04 unverdicted novelty 6.0

    Global well-posedness is established for a nonlocal curve evolution of an immersed elastic filament, together with convergence to resistive force theory as the filament thickness approaches zero.

Reference graph

Works this paper leans on

57 extracted references · 57 canonical work pages · cited by 1 Pith paper

  1. [1]

    Alazard, N

    T. Alazard, N. Burq, and C. Zuily. On the cauchy problem for gravity water waves.Inventiones mathematicae, 198(1):71–163, 2014

    work page 2014

  2. [2]

    Alazard and O

    T. Alazard and O. Lazar. Paralinearization of the muskat equation and application to the cauchy problem. Archive for Rational Mechanics and Analysis, 237(2):545–583, 2020

    work page 2020

  3. [3]

    Alazard and G

    T. Alazard and G. M´ etivier. Paralinearization of the dirichlet to neumann operator, and regularity of three- dimensional water waves.Communications in Partial Differential Equations, 34(12):1632–1704, 2009

    work page 2009

  4. [4]

    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

  5. [5]

    H. I. Andersson, E. Celledoni, L. Ohm, B. Owren, and B. K. Tapley. An integral model based on slender body theory, with applications to curved rigid fibers.Physics of Fluids, 33(4):041904, 2021

    work page 2021

  6. [6]

    Bahouri, J.-Y

    H. Bahouri, J.-Y. Chemin, and R. Danchin.Fourier Analysis and Nonlinear Partial Differential Equations. Springer, 2011

    work page 2011

  7. [7]

    Batchelor

    G. Batchelor. Slender-body theory for particles of arbitrary cross-section in Stokes flow.J. Fluid Mech., 44(3):419– 440, 1970

    work page 1970

  8. [8]

    J. T. Beale, T. Y. Hou, and J. S. Lowengrub. Growth rates for the linearized motion of fluid interfaces away from equilibrium.Communications on Pure and Applied Mathematics, 46(9):1269–1301, 1993

    work page 1993

  9. [9]

    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

  10. [10]

    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

  11. [11]

    Cameron and R

    S. Cameron and R. M. Strain. Critical local well-posedness for the fully nonlinear peskin problem.Communica- tions on Pure and Applied Mathematics, 77(2):901–989, 2024

    work page 2024

  12. [12]

    Chen and Q.-H

    K. Chen and Q.-H. Nguyen. The peskin problem with initial data.SIAM Journal on Mathematical Analysis, 55(6):6262–6304, 2023

    work page 2023

  13. [13]

    Cortez and M

    R. Cortez and M. Nicholas. Slender body theory for Stokes flows with regularized forces.Commun. Appl. Math. Comput. Sci., 7(1):33–62, 2012

    work page 2012

  14. [14]

    R. Cox. The motion of long slender bodies in a viscous fluid part 1. general theory.J. Fluid Mech., 44(4):791–810, 1970

    work page 1970

  15. [15]

    L. Euler. Methodus inveniendi curvas lineas maximi minimive proprietate gaudentes sive solution problematis isometrici latissimo sensu accepti. 1744

    work page

  16. [16]

    P. T. Flynn and H. Q. Nguyen. The vanishing surface tension limit of the muskat problem.Communications in Mathematical Physics, 382:1205–1241, 2021

    work page 2021

  17. [17]

    Gancedo, R

    F. Gancedo, R. Granero-Belinch´ on, and S. Scrobogna. Global existence in the lipschitz class for the n-peskin problem.arXiv preprint arXiv:2011.02294, 2020

    work page arXiv 2011

  18. [18]

    Garc´ ıa-Ju´ arez and S

    E. Garc´ ıa-Ju´ arez and S. V. Haziot. Critical well-posedness for the 2d peskin problem with general tension.arXiv preprint arXiv:2311.10157, 2023

    work page arXiv 2023

  19. [19]

    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

  20. [20]

    Garc´ ıa-Ju´ arez, Y

    E. Garc´ ıa-Ju´ arez, Y. Mori, and R. M. Strain. The peskin problem with viscosity contrast.Analysis & PDE, 16(3):785–838, 2023

    work page 2023

  21. [21]

    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

  22. [22]

    Gray and G

    J. Gray and G. Hancock. The propulsion of sea-urchin spermatozoa.Journal of Experimental Biology, 32(4):802– 814, 1955

    work page 1955

  23. [23]

    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

  24. [24]

    T. Y. Hou, J. S. Lowengrub, and M. J. Shelley. Removing the stiffness from interfacial flows with surface tension. Journal of Computational Physics, 114(2):312–338, 1994

    work page 1994

  25. [25]

    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

  26. [26]

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

    work page 1980

  27. [27]

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

    work page 1976

  28. [28]

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

    work page 1976

  29. [29]

    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. THE SLENDER BODY FREE BOUNDARY PROBLEM 35

    work page 1992

  30. [30]

    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

  31. [31]

    D. Lannes. Well-posedness of the water-waves equations.Journal of the American Mathematical Society, 18(3):605–654, 2005

    work page 2005

  32. [32]

    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

  33. [33]

    R. Levien. The elastica: a mathematical history. Technical report, Technical Report No. UCB/EECS-2008-103, 2008

    work page 2008

  34. [34]

    L. Li, H. Manikantan, D. Saintillan, and S. E. Spagnolie. The sedimentation of flexible filaments.J. Fluid Mech., 735:705–736, 2013

    work page 2013

  35. [35]

    Lin and J

    F.-H. Lin and J. Tong. Solvability of the stokes immersed boundary problem in two dimensions.Communications on Pure and Applied Mathematics, 72(1):159–226, 2019

    work page 2019

  36. [36]

    Matsutani

    S. Matsutani. Euler’s elastica and beyond.Journal of Geometry and Symmetry in Physics, 17:45–86, 2010

    work page 2010

  37. [37]

    Maxian, A

    O. Maxian, A. Mogilner, and A. Donev. Integral-based spectral method for inextensible slender fibers in stokes flow.Physical Review Fluids, 6(1):014102, 2021

    work page 2021

  38. [38]

    Maxian, B

    O. Maxian, B. Sprinkle, C. S. Peskin, and A. Donev. Hydrodynamics of a twisting, bending, inextensible fiber in stokes flow.Physical Review Fluids, 7(7):074101, 2022

    work page 2022

  39. [39]

    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

  40. [40]

    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

  41. [41]

    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

  42. [42]

    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

  43. [43]

    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

  44. [44]

    Y. Mori, A. Rodenberg, and D. Spirn. Well-posedness and global behavior of the peskin problem of an immersed elastic filament in stokes flow.Communications on Pure and Applied Mathematics, 72(5):887–980, 2019

    work page 2019

  45. [45]

    H. Q. Nguyen and B. Pausader. A paradifferential approach for well-posedness of the muskat problem.Archive for Rational Mechanics and Analysis, 237(1):35–100, 2020

    work page 2020

  46. [46]

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

    work page 2011

  47. [47]

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

    work page arXiv 2024

  48. [48]

    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

  49. [49]

    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

  50. [50]

    Pozrikidis.Boundary integral and singularity methods for linearized viscous flow

    C. Pozrikidis.Boundary integral and singularity methods for linearized viscous flow. Cambridge University Press, 1992

    work page 1992

  51. [51]

    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

  52. [52]

    S. E. Spagnolie and E. Lauga. Comparative hydrodynamics of bacterial polymorphism.Phys. Rev. Lett., 106(5):058103, 2011

    work page 2011

  53. [53]

    J. Tong. Regularized stokes immersed boundary problems in two dimensions: Well-posedness, singular limit, and error estimates.Communications on Pure and Applied Mathematics, 74(2):366–449, 2021

    work page 2021

  54. [54]

    Tong and D

    J. Tong and D. Wei. Geometric properties of the 2-d peskin problem.arXiv preprint arXiv:2304.09556, 2023

    work page arXiv 2023

  55. [55]

    Tornberg and K

    A.-K. Tornberg and K. Gustavsson. A numerical method for simulations of rigid fiber suspensions.J. Comput. Phys., 215(1):172–196, 2006

    work page 2006

  56. [56]

    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

  57. [57]

    C. H. Wiggins and R. E. Goldstein. Flexive and propulsive dynamics of elastica at low reynolds number.Physical Review Letters, 80(17):3879, 1998. Department of Mathematics, University of Wisconsin - Madison, Madison, WI 53706 Email address:lohm2@wisc.edu

    work page 1998