Kelvin wave and soliton propagation in classical viscous vortex filaments

This paper uses numerical simulations of the three-dimensional incompressible Navier-Stokes equations to demonstrate that Kelvin waves in classical viscous vortex filaments align with Lord Kelvin's predictions, while also revealing the existence and collision dynamics of solitons and proposing a feasible experimental method for their generation.

The Gist

Imagine a fluid, like water or air, not as a smooth, continuous sheet, but as a collection of invisible, twisting ropes. These ropes are called vortex filaments. You've seen them in action: the swirling funnel of a tornado, the smoke ring from a chimney, or the tiny spirals left behind by an airplane wing.

This paper is a scientific detective story about what happens when you poke, prod, or wiggle these invisible ropes. The researchers, Elio Sterkers and Giorgio Krstulovic, used powerful computer simulations to see how these fluid ropes behave, specifically looking for two special types of "waves" that travel along them.

Here is the story of their discovery, broken down into simple concepts.

1. The Invisible Ropes and the "Local" Rule

First, imagine a long, thin rope made of spinning water. In the past, scientists used a simplified rule called the Local Induced Approximation (LIA) to predict how these ropes move.

Think of the LIA like this: Imagine the rope is made of tiny beads. The LIA rule says, "To know where a bead moves next, just look at the curve right next to it." It ignores the rest of the rope. It's like driving a car and only looking at the bumper in front of you, ignoring the road ahead.

While this rule is a bit "crude" (it misses some long-distance effects), it predicts two very cool things:

1. Kelvin Waves: Ripples that travel along the rope like a snake slithering.
2. Solitons: A special, self-contained "hump" or pulse that travels down the rope without changing shape, like a perfect wave in a canal that never breaks.

2. The Big Question: Do These Exist in "Real" Fluids?

The tricky part is that the LIA rule is a mathematical fantasy. Real fluids (like water in a tank or air in the sky) have viscosity (thickness/stickiness). In the real world, things usually slow down and fade away due to friction.

The big question the authors asked was: "Do these perfect mathematical waves actually exist in real, sticky fluids, or are they just math ghosts?"

To find out, they didn't use a bathtub; they used a super-computer to simulate the full, complex laws of fluid motion (the Navier-Stokes equations), which account for all the stickiness and friction.

3. The Findings: The Waves Are Real!

The Snake Ripples (Kelvin Waves)

First, they tested the "snake" ripples. They shook the virtual rope and watched the ripples travel.

• The Result: The ripples moved exactly as the 19th-century physicist Lord Kelvin predicted over 100 years ago.
• The Analogy: It's like plucking a guitar string. Even though the air is "sticky," the note still rings out with the exact pitch the mathematician said it should. This confirmed that even in real fluids, these waves are a fundamental part of how energy moves.

The Self-Healing Hump (Solitons)

Next, they tried to create the "soliton"—that perfect, shape-shifting pulse. They shaped the virtual rope into a specific hump and let it go.

• The Result: The hump traveled down the rope, keeping its shape for a long time, even though the fluid was sticky. It was like a surfer riding a wave that never seemed to lose speed or flatten out.
• The Collision: They sent two of these humps toward each other. In a perfect mathematical world, they would pass right through each other like ghosts. In their simulation, they collided, twisted the rope, and briefly created a tiny ring of fluid that popped off (like a bubble popping). But after the chaos, two smaller humps survived and kept traveling!
• The Analogy: Imagine two people running toward each other on a trampoline. They jump, collide, and for a split second, the trampoline gets messy. But then, they both bounce back up and keep running in opposite directions, still intact.

4. How to Make One in a Lab (The "Ring Throw")

The most exciting part is their proposal for how to actually see this in a real laboratory.

They realized that a soliton carries momentum (push). So, if you want to create one, you need to give the rope a sudden push.

• The Experiment: Imagine you have a long, vertical rope of spinning water. You take a small, horizontal ring of spinning water (like a smoke ring) and throw it at the rope.
• The Magic: When the ring hits the rope, they reconnect (like two rubber bands snapping together). The ring disappears, but its "push" is transferred to the rope.
• The Result: That sudden push creates a soliton that zooms away down the rope, carrying the energy of the thrown ring.

Why Does This Matter?

This isn't just about cool physics tricks. It helps us understand:

• Tornadoes and Weather: Since tornadoes are giant vortex filaments, understanding how energy moves along them (via these waves) could help us predict their behavior.
• Airplane Safety: Understanding how these "ropes" interact helps engineers design planes that don't get caught in dangerous wake turbulence.
• Energy Transfer: It shows us how energy moves from big swirls to tiny, invisible swirls, which is a key mystery in how fluids become turbulent.

The Bottom Line

The paper proves that even in our messy, sticky, real-world fluids, nature still follows some beautiful, mathematical rules. The "perfect" waves and pulses predicted by mathematicians a century ago aren't just theory; they are real, observable phenomena that we can create, study, and perhaps one day use to predict the weather or design better aircraft.

Technical Summary

1. Problem Statement

Vortex filaments are highly rotating, localized fluid structures found in both classical fluids (e.g., tornadoes, aircraft wakes) and quantum fluids (e.g., superfluid helium, Bose-Einstein condensates). While the dynamics of these filaments in superfluids are well-studied due to their topological nature and quantized circulation, their behavior in classical viscous fluids remains less understood, particularly regarding nonlinear excitations.

The study addresses three primary questions:

1. Kelvin Waves (KWs): Can the dispersion relation of helical excitations (Kelvin waves) propagating along classical viscous vortex filaments be measured and does it align with Lord Kelvin's theoretical predictions derived for inviscid fluids?
2. Solitons: Do Hasimoto solitons (nonlinear solitary waves predicted by the Local Induced Approximation, LIA) exist and remain stable in classical viscous fluids governed by the Navier-Stokes equations, despite the breakdown of the LIA's assumptions (e.g., ignoring non-local interactions and viscosity)?
3. Experimental Feasibility: Can such solitons be generated in a controlled laboratory setting, and what is the mechanism for their creation?

2. Methodology

The authors employed Direct Numerical Simulations (DNS) of the three-dimensional incompressible Navier-Stokes (NS) equations to model vortex filaments in classical fluids.

• Governing Equations: The simulations solved the NS equations for velocity $\mathbf{v}$ and pressure $p$ with kinematic viscosity $\nu$:
  $$\partial_t \mathbf{v} + \mathbf{v} \cdot \nabla \mathbf{v} = -\nabla p + \nu \nabla^2 \mathbf{v}, \quad \nabla \cdot \mathbf{v} = 0$$
• Initial Conditions:
  • Vorticity fields were initialized using a line integral representation of a vortex filament with a regularized core model (solid-body rotation) of size $a_0$.
  • Kelvin Wave Study: Perturbations were introduced along an almost straight filament with random phases to analyze wave propagation.
  • Soliton Study: Initial conditions were set using the analytical Hasimoto soliton solution derived from the LIA equation.
  • Collision Study: Two conjugate solitons were initialized to study their interaction.
  • Generation Study: A small vortex ring was thrown toward a straight filament to induce reconnection and trigger a soliton.
• Numerical Setup:
  • Code: GHOST (pseudo-spectral code).
  • Resolution: $512^2 \times 1024$ grid points.
  • Domain: Periodic box with four alternating-sign vortices to minimize boundary effects.
  • Reynolds Numbers: Simulations were conducted at $Re_\nu = \Gamma/\nu \approx 10^3 - 2.5 \times 10^3$.
• Analysis Techniques:
  • Vortex position tracking via vorticity-weighted averaging.
  • Fourier analysis in space and time to extract dispersion relations.
  • Comparison of numerical results against analytical predictions from the LIA and Lord Kelvin's theory.

3. Key Contributions and Results

A. Validation of Kelvin Wave Dispersion

• The authors successfully measured the dispersion relation of Kelvin waves in a viscous classical fluid.
• Result: The numerical data showed excellent agreement with Lord Kelvin's theoretical prediction for an inviscid solid-body rotation core:
  $$\omega_k = \frac{\Gamma}{4\pi} k^2 \left( b - \log(ka_0) \right)$$
• Observation: While viscous effects caused the vortex core size $a_0$ to increase over time (blurring the curves), the dispersion relation remained robust for moderate Reynolds numbers, confirming that KWs can propagate with minimal damping in classical fluids.

B. Existence and Stability of Hasimoto Solitons

• The study numerically demonstrated the existence of Hasimoto solitons in the Navier-Stokes regime.
• Propagation: The soliton propagated along the filament, maintaining its shape and oscillating coherently.
• Viscous Effects: Over time, the soliton's amplitude decayed slowly while its width increased, consistent with viscous dissipation. However, it traveled more than ten times its original width before significant distortion, proving its stability in a viscous environment.
• Parameter Fitting: By fitting the soliton speed, the authors determined a phenomenological LIA constant $\Lambda_{fit} \approx 1.857$, validating the LIA model's predictive power even in viscous flows.

C. Soliton Collision and Reconnection

• Elastic vs. Inelastic: Unlike the ideal 1D Nonlinear Schrödinger (1dNLS) equation where solitons collide elastically, the classical viscous collision was inelastic.
• Mechanism: When two solitons collided, the opposing filament segments approached closely, triggering a vortex reconnection.
• Outcome: A vortex ring was ejected from the filament, and the original solitons survived as smaller structures (roughly 5 times smaller in amplitude) traveling away from the collision site. This highlights the role of non-local interactions and reconnection in breaking the integrability of the system in real fluids.

D. Experimental Generation Proposal

• The authors proposed and numerically validated a method to generate solitons in the lab: colliding a vortex ring with a straight vortex filament.
• Mechanism: The ring transfers momentum to the filament. Upon reconnection, the excess momentum is absorbed by the filament, triggering the formation of a soliton.
• Quantitative Agreement: The generated soliton's amplitude ($A$) and width ($\lambda$) matched theoretical estimates derived from momentum and energy conservation laws (with fitted loss coefficients $\gamma_1 \approx 0.41, \gamma_2 \approx 0.32$). The soliton speed also matched the theoretical prediction using the fitted $\Lambda$.

4. Significance

• Bridging Classical and Quantum Hydrodynamics: The work establishes that phenomena previously thought to be exclusive to superfluids (quantized vortices) or idealized mathematical models (LIA) are observable in classical viscous fluids.
• Wave Turbulence: By confirming that Kelvin waves propagate with low damping in classical fluids, the study opens new avenues for applying wave turbulence theory to classical vortex filaments. This is crucial for understanding energy cascades toward small scales in turbulence.
• Experimental Roadmap: The paper provides a concrete, feasible protocol for generating and observing vortex solitons in laboratory experiments (e.g., water tanks), moving the field from purely theoretical/superfluid contexts to accessible classical fluid dynamics.
• Relevance to Meteorology and Safety: Since vortex filaments are central to tornado dynamics and aircraft safety, understanding these nonlinear excitations and their generation mechanisms (via reconnection) could improve forecasting and safety protocols.

In conclusion, Sterkers and Krstulovic demonstrate that despite the strong approximations of the LIA model, its predictions regarding Kelvin waves and solitons remain phenomenologically valid in classical viscous fluids, provided the Reynolds number is sufficiently high. This validates the use of integrable system theory as a robust tool for analyzing complex fluid dynamics.