Lyapunov Unstable Motion Bifurcating from a Circular Vortex Filament

1. Planteamiento del Problema

El artículo aborda la dinámica y estabilidad de los filamentos de vórtice cerrados en $\mathbb{R}^3$, gobernados por la Ecuación de Inducción Localizada (Localized Induction Equation - LIE):
$$x_t = x_s \times x_{ss}$$
donde $x(s,t)$ es el vector de posición del filamento parametrizado por la longitud de arco $s$.

El foco específico es la estabilidad del filamento circular $x^R$, que representa un anillo de vórtice. El problema central surge de una distinción sutil pero crítica en la teoría de estabilidad:

• Estabilidad Orbital: La solución perturbada permanece cerca de la órbita del filamento circular (es decir, cerca del conjunto de todas las rotaciones y traslaciones del círculo).
• Estabilidad de Lyapunov: La solución perturbada permanece cerca de la trayectoria específica del filamento circular en el tiempo, sin permitir el "desfase" por simetrías.

En un trabajo previo [2], los autores demostraron que el filamento circular es orbitalmente estable bajo perturbaciones asimétricas, pero inestable de Lyapunov debido al crecimiento lineal de los modos de traslación. Sin embargo, no se habían construido explícitamente soluciones que ilustraran esta brecha: soluciones que se mantengan uniformemente cerca de la órbita (satisfaciendo estimaciones de estabilidad orbital) pero que se alejen secularmente de la trayectoria de referencia (violando la estabilidad de Lyapunov).

2. Metodología

Para abordar este problema, los autores emplean una combinación de análisis geométrico, teoría de bifurcación y análisis funcional:

1. Ansatz de Movimiento de Tornillo Axial (Axial Screw Motion):
   Se define una clase específica de soluciones donde la forma del filamento permanece invariante en el tiempo, pero el filamento rota alrededor del eje de simetría ($e_3$), se traslada a lo largo de este eje y experimenta un deslizamiento tangencial constante. La forma general es:
   $$x(s, t) = Q_{\Omega t}^z y(s + ct) + V t e_3$$
   donde $y(s)$ es la forma del perfil, $\Omega$ es la velocidad angular, $c$ es la velocidad de deslizamiento tangencial y $V$ es la velocidad de traslación axial.

2. Reducción del Sistema:
   Sustituyendo el ansatz en la LIE, el problema de evolución parcial se reduce a una ecuación diferencial ordinaria (ODE) no lineal para el perfil $y(s)$:
   $$y_s \times y_{ss} = c y_s + \Omega e_3 \times y + V e_3$$
   Los autores reformulan este sistema en términos de perturbaciones $(z, \delta V)$ alrededor del perfil circular $x_0^R$, descomponiendo la perturbación en componentes tangencial ($u$), normal ($v$) y binormal ($w$) respecto al marco de Frenet-Serret del círculo.

3. Teoría de Bifurcación de Crandall-Rabinowitz:
   El núcleo de la demostración es aplicar el teorema de bifurcación de Crandall-Rabinowitz.

   • Se identifica un valor crítico de velocidad angular $\Omega_k = \frac{\sqrt{k^2-1}}{R^2}$ para $k \geq 2$.
   • Se analiza el operador linealizado alrededor de la solución trivial (el círculo).
   • Se demuestra que en $\Omega_k$, el operador linealizado tiene un núcleo unidimensional (bifurcación simple), lo que permite la existencia de una rama de soluciones no triviales que emergen del círculo.

4. Análisis de Estabilidad y Deriva Secular:
   Una vez construidas las soluciones bifurcadas, se analiza su comportamiento asintótico en el tiempo comparando la velocidad de traslación $V_\lambda$ de la nueva solución con la velocidad $1/R$ del círculo original.

3. Contribuciones Clave

• Construcción Explícita de Soluciones: Los autores prueban la existencia de una familia continua de soluciones no triviales (movimientos de tornillo axial) que bifurcan del filamento circular.
• Cuantificación de la Inestabilidad de Lyapunov: Demuestran que estas nuevas soluciones, aunque permanecen cerca de la variedad de simetrías del círculo (órbita), se alejan linealmente en el tiempo de la trayectoria de referencia específica debido a una diferencia en la velocidad de traslación axial.
• Resolución de la Brecha Teórica: El trabajo proporciona el primer ejemplo explícito que satisface las estimaciones de estabilidad orbital pero falla en la estabilidad de Lyapunov para este sistema físico, clarificando la relación entre ambos conceptos en el contexto de la LIE.
• Diferenciación con Trabajos Previos: Se distingue claramente de las soluciones de Kida [29] y García-Vega [19]. Mientras que estos últimos construyen soluciones no triviales, su enfoque en el vector tangente pierde información sobre la velocidad de traslación global, impidiéndoles distinguir entre estabilidad orbital y de Lyapunov. Este trabajo, al trabajar con la posición completa, captura dicha deriva.

4. Resultados Principales

El Teorema 1.5 es el resultado central del artículo:

• Existencia de Bifurcación: Para cada entero $k \geq 2$, existe un valor crítico de velocidad angular $\Omega_k = \frac{\sqrt{k^2-1}}{R^2}$ y un $\varepsilon_k > 0$ tal que existe un mapeo suave $\lambda \mapsto (\Omega_\lambda, y_\lambda, c_\lambda, V_\lambda)$ que define una familia de movimientos de tornillo axial $x(\lambda)$.
• Estabilidad Orbital (Propiedad 1): Para cualquier $\lambda$ pequeño, la distancia entre la solución $x(\lambda)(t)$ y la variedad de simetrías $\Sigma$ (círculos rotados y trasladados) está acotada uniformemente en el tiempo:
  $$\text{dist}(x(\lambda)(t), \Sigma) \leq C|\lambda|$$
  Esto confirma que la solución no "explota" ni se aleja de la forma circular en términos de geometría relativa.
• Inestabilidad de Lyapunov (Propiedad 2): Para $\lambda \neq 0$, existe una deriva secular. La distancia entre la solución perturbada y la trayectoria de referencia $x^R(t)$ crece linealmente con el tiempo:
  $$|x(\lambda)(s, t) - x^R(s, t)| \geq \gamma t \quad \text{para } t \to \infty$$
• Mecanismo de Deriva: La causa de esta inestabilidad es la diferencia en la velocidad de traslación a lo largo del eje de simetría. La velocidad $V_\lambda$ de la solución bifurcada se expande como:
  $$V_\lambda = \frac{1}{R} - \left( \frac{k^2(k^2-1)}{2R^3} \right)\lambda^2 + o(\lambda^2)$$
  Dado que $V_\lambda < 1/R$ para $\lambda \neq 0$, la solución se mueve más lentamente que el círculo de referencia, acumulando una separación lineal en el tiempo.

5. Significado e Impacto

Este trabajo es fundamental para la comprensión matemática de la dinámica de fluidos, específicamente para la ecuación de Euler y sus aproximaciones (como la LIE):

1. Rigor en la Estabilidad: Establece rigurosamente que la estabilidad orbital no implica estabilidad de Lyapunov para filamentos de vórtice, incluso cuando las perturbaciones son pequeñas y la forma geométrica se mantiene casi inalterada.
2. Implicaciones Físicas: Sugiere que en flujos reales donde se modelan anillos de vórtice, pequeñas asimetrías pueden llevar a que el anillo se desplace a una velocidad diferente a la predicha por el modelo de anillo perfecto, acumulando errores de posición significativos a largo plazo, aunque su forma siga siendo un anillo.
3. Herramientas Analíticas: La aplicación exitosa de la teoría de bifurcación de Crandall-Rabinowitz a un sistema de ecuaciones diferenciales parciales no lineales con restricciones geométricas (conservación de la longitud de arco) ofrece un marco metodológico para estudiar otros problemas de estabilidad en dinámica de fluidos geométrica.

En resumen, el artículo cierra una brecha teórica importante al construir explícitamente soluciones que demuestran cómo un sistema puede ser estable en su forma (orbitalmente) pero inestable en su posición absoluta (Lyapunov), proporcionando una comprensión más profunda de la dinámica de los filamentos de vórtice.