🔬 condensed matter

A theory of locally impenetrable elastic tubes

This paper presents a reduced-order variational theory for locally impenetrable elastic tubes with circular cross-sections, deriving governing equations and demonstrating that their configurations consist of standard Kirchhoff rod segments connected to regions of constant Frenet curvature.

Original authors: Krishnan Suryanarayanan, Harmeet Singh

Published 2026-07-30

📖 8 min read🧠 Deep dive

CC BY 4.0

Original authors: Krishnan Suryanarayanan, Harmeet Singh

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a world where the rules of physics are a bit more polite than usual. In our everyday reality, if you try to push your hand through a solid wall, it stops. Matter simply refuses to occupy the same space as other matter. This is a fundamental rule of the universe: things can't pass through themselves. Scientists call this "impenetrability." Usually, when engineers design things like bridges, airplane wings, or even the tiny springs in your watch, they assume this rule is always true, but they don't always check it explicitly in their math. They just hope the materials don't get too squished.

However, there is a special branch of physics that studies long, thin, flexible things—like a garden hose, a piece of spaghetti, or a DNA strand. These are called "slender bodies." When these things bend, they can get into trouble. If you bend a hose too sharply, the inside walls might crash into each other, which is impossible in real life. Standard math models for these objects often ignore this crash until it's too late, leading to predictions that look like the hose is folding through itself like a ghost. This paper dives into the messy, tricky middle ground: what happens when a long, bendy tube bends so hard that it almost touches itself, and how we can write new rules to stop that from happening in our calculations.

The authors, Krishnan Suryanarayanan and Harmeet Singh, have built a new mathematical theory for these "locally impenetrable elastic tubes." Think of it as a rulebook for a very picky, very bendy noodle that refuses to let its own skin touch its own skin, even for a split second.

The Problem with the "Ghost Noodle"

To understand their discovery, imagine you have a long, flexible garden hose hanging from a hook. If you pull the two ends closer together, the hose sags in the middle. As you keep pulling the ends in, the sag gets deeper and the curve at the bottom gets tighter and tighter.

In the old, standard way of doing math (called Kirchhoff rod theory), if you pull the ends close enough, the model predicts the curve at the bottom will get infinitely sharp. It's like the hose tries to fold into a perfect, sharp point. In the real world, this is impossible because the hose has thickness; the inside of the bend would hit the outside of the bend before it could get that sharp. The old math ignores this and says, "Sure, it folds into a point!" resulting in a "ghost noodle" that passes through itself.

The authors asked: What if we force the math to respect the rule that the hose cannot touch itself? What does the shape look like then?

The "Active" and "Inactive" Zones

Their big idea is that when the hose gets close to touching itself, it doesn't just bend smoothly anymore. Instead, it splits into two different types of behavior, like a tube with two different personalities.

The Inactive Zone: This is the part of the tube that is bending gently. Here, the tube acts like a normal, standard spring. It follows the usual rules of bending, and the math is simple and linear.
The Active Zone: This is the part where the tube is bending as hard as it possibly can without touching itself. In this zone, the tube hits a "speed limit" for how sharp it can turn. The authors found that once the tube hits this limit, it doesn't just stop bending; it locks into a shape with a constant curve.

Imagine you are trying to bend a stiff wire. If you try to bend it too hard, it resists. But in this new theory, the tube doesn't just resist; it changes its internal rules. In the "active" zone, the tube becomes a perfect arc of a circle. It's as if the tube says, "I will bend this much, and no more. If you push harder, I will just make this curved section longer, but I won't get any sharper."

How They Figured It Out

The authors used a clever mathematical trick called a "variational scheme." You can think of this as a game of finding the lowest energy state. Imagine the tube is trying to find the most comfortable position to hang in, like a cat looking for the perfect spot on a couch.

Usually, the cat just lies down. But in this game, there's a rule: "You cannot curl up so tight that your tail touches your nose." The authors added this rule to the math using a special "slack function." This function acts like a safety valve. If the tube tries to bend too much, the safety valve kicks in, and the math forces the tube to switch from its normal bending mode to this new "constant curve" mode.

They discovered that the tube creates a specific region where this constant curve happens. This region starts at the point of maximum stress (the bottom of the hang) and spreads out as you pull the ends closer. The boundaries between the "normal bending" part and the "constant curve" part are sharp, almost like a step in a staircase, rather than a smooth ramp.

What They Found in Three Scenarios

To prove their theory works, they tested it on three different scenarios, and the results were quite surprising compared to the old models.

1. The Hanging Flexible Tube (The Catenary)
First, they looked at a tube that has no stiffness at all (a "fully flexible" tube), like a heavy rope. In the old math, if you pull the ends of a rope close together, the bottom gets infinitely sharp. In their new theory, the rope hits the "impenetrability limit" and forms a perfect circular arc at the bottom.

The Surprise: In the old model, the rope has no internal bending force (moment) because it's perfectly flexible. But in the new model, because the tube is forced to stay in that constant curve shape, it actually generates internal forces and moments just to hold that shape. The "active" part of the rope suddenly becomes stiff, even though the material itself is floppy.

2. The Soft Elastic Tube
Next, they looked at a tube that does have some stiffness, like a soft rubber hose. They found that if the rubber is soft enough, it will still try to bend too sharply.

The Surprise: Even though the rubber is trying to be flexible, the "no-touching" rule kicks in. They identified a specific threshold (a scaling parameter) where the tube switches from behaving like a normal spring to forming that constant-curve "active" zone. If the tube is stiffer than this threshold, it never hits the limit. But if it's softer, the limit takes over, and the tube forms a circular arc at the bottom, just like the floppy rope.

3. The Twisted 3D Tube
Finally, they took a tube and twisted it like a pretzel. In the old math, as you twist a rod, the bending tends to concentrate in a tiny, sharp spot in the middle, eventually causing the rod to poke through itself.

The Surprise: With the new theory, instead of a sharp, impossible point, the tube nucleates a helix (a spiral) in the middle. As you twist it more, this spiral doesn't get sharper; it just gets longer, spreading out along the rod. The tube essentially turns a dangerous, sharp kink into a safe, spreading spiral.

Why This Matters

This isn't just about garden hoses. The authors suggest this theory could help us understand things like knots, clasps, and even how DNA packs itself inside a cell. When you tie a knot, the rope has to bend around itself. If the rope is thick, it can't bend as sharply as a thin thread. This new theory provides a way to calculate the "ideal shape" of a knot or a clasp without the math breaking down or predicting impossible shapes.

They also mention that this could help explain some weird physics puzzles, like why a falling chain sometimes accelerates faster than gravity (the "chain fountain" effect). The idea is that the links in the chain can't bend past a certain point, and this "impenetrability" creates a force that pushes the chain up.

In short, Suryanarayanan and Singh have given us a better map for navigating the world of bendy things. They showed that when objects get close to touching themselves, they don't just break the rules; they change the rules entirely, creating new zones of behavior that are just as fascinating as the bending itself. The tube doesn't just bend; it adapts, locking into a perfect curve to keep its own skin safe.

Technical Summary: A Theory of Locally Impenetrable Elastic Tubes

Problem Statement
In continuum mechanics, the assumption that matter cannot penetrate itself is fundamental. While global impenetrability (preventing distinct material points from overlapping) is often the primary focus in modeling self-contacting bodies like knots and clasps, local impenetrability (preventing infinitesimally close points from overlapping) is typically treated merely as a post-hoc criterion to discard unphysical solutions rather than an actively enforced constraint. Standard Kirchhoff rod theory (KRT), widely used for slender bodies, reduces a tube to an inextensible centerline with an orthonormal director frame. However, KRT does not inherently prevent the centerline curvature from exceeding a limit dictated by the tube's cross-sectional radius, potentially leading to configurations where the tube's material locally penetrates itself. This paper addresses the need for a reduced-order theory that actively enforces local impenetrability for elastic tubes with uniform circular cross-sections.

Methodology
The authors develop a variational framework to derive a theory for locally impenetrable elastic tubes. The methodology proceeds as follows:

Kinematic Reduction: Starting from the 3D deformation gradient $F$ , the condition of local impenetrability ( $\det F > 0$ ) is reduced to an inequality constraint on the Frenet curvature ( $\kappa$ ) of the tube's centerline. For an inextensible tube of radius $t$ , this yields the constraint $\kappa < 1/t$ .
Variational Formulation: The inequality constraint is converted into an equality constraint using a slack function $g(s)$ , such that $\kappa - \hat{\kappa} + g(s)^2 = 0$ (where $\hat{\kappa} = 1/t$ ). This constraint is incorporated into the action functional via a Lagrange multiplier $\Lambda(s)$ .
Derivation of Governing Equations: By taking variations of the action functional with respect to the position, director frame, and slack function, the authors derive:
- Force and Moment Balance: Standard equilibrium equations valid throughout the tube.
- Constitutive Relations: A piecewise definition of the internal moment. In "inactive" regions (where $\kappa < \hat{\kappa}$ ), the relation is linear (standard KRT). In "active" regions (where $\kappa = \hat{\kappa}$ ), the relation becomes nonlinear, with the Lagrange multiplier $\Lambda$ acting as a variable stiffness term to maintain the curvature bound.
- Jump Conditions: Continuity conditions for position, tangent, internal force, internal moment, and the Hamiltonian function at the interfaces between active and inactive regions.
Conservation Laws: The authors demonstrate that standard first integrals of KRT (conservation of force, moment, twist, and Hamiltonian) persist across both active and inactive regions, despite the modified constitutive behavior.

Key Contributions

Active Enforcement of Local Impenetrability: The paper provides a rigorous mathematical framework where the local impenetrability constraint is not just a filter but an active part of the equilibrium equations, leading to the nucleation and propagation of "active regions."
Two-Phase Constitutive Behavior: The theory reveals that a deformed tube comprises segments of standard Kirchhoff solutions connected to segments of constant Frenet curvature. The transition between these "phases" is governed by jump conditions derived naturally from the variational principle.
Nonlinear Constitutive Response: In active regions, the relationship between internal moments and bending strains becomes nonlinear. Notably, for fully flexible tubes (zero bending stiffness), the active constraint generates non-zero internal moments and shear forces, phenomena absent in standard flexible string theory.
Preservation of Integrals: The work proves that the fundamental conservation laws of Kirchhoff rod theory remain valid even in the presence of the active local impenetrability constraint.

Results
The theory is illustrated through three numerical examples:

Fully Flexible Tube Hanging Under Self-Weight:
- When the ends of a flexible tube are brought close, the curvature at the lowest point reaches the limit $\hat{\kappa}$ .
- An active region of constant curvature nucleates and spreads.
- Unlike standard catenary solutions, the active region supports non-zero bending moments and shear forces.
- The tension in the active region is lower than predicted by standard theory.
- The curvature exhibits a jump discontinuity at the boundary of the active region, consistent with regularity results for self-contacting tubes.
Elastic Tube Hanging Under Self-Weight:
- For extremely soft materials (low bending modulus), the peak curvature can still reach the impenetrability bound even with bending elasticity.
- Once the bound is reached, an active region forms.
- While the differences in shape and internal forces are small compared to standard KRT for soft materials, the theory captures the activation of the constraint where standard theory would predict unphysical penetration.
Highly Twisted 3D Elastic Tube:
- Under increasing twisting moments, standard KRT predicts curvature concentrating into a small region, eventually leading to local penetration.
- The present theory shows that as the curvature bound is reached, a helical active region nucleates at the mid-length and spreads along the rod.
- Significant quantitative and qualitative differences are observed in the deformed configurations, bending strains, and internal force/moment distributions compared to standard KRT, particularly in the active region.

Significance and Claims
The authors claim that this theory offers a necessary refinement for modeling slender bodies where local self-avoidance is critical. The paper positions its work as a step toward computing ideal shapes of knots and clasps that satisfy both global and local impenetrability, noting that previous models often violated local bounds.

The authors suggest the theory's relevance extends to:

Chain Paradoxes: Modeling systems like falling chains or "chain fountains" where discrete links cannot bend beyond a maximum curvature due to local impenetrability, potentially explaining faster-than-gravity accelerations.
Optimal Packing: Addressing the packing of filaments in confined volumes (e.g., DNA), where local impenetrability imposes hard constraints on bending that global constraints alone do not capture.

The paper concludes modestly, stating that while it successfully incorporates local impenetrability, the ultimate extension would involve incorporating global contact to form a complete theory of self-avoiding elastic rods, a task acknowledged as non-trivial.