Sheets of Spectral Data of Stokes Waves in Weakly Nonlinear Models

Resumen Técnico

1. Planteamiento del Problema

El estudio se centra en la estabilidad espectral de ondas de Stokes de pequeña amplitud ($\varepsilon$) en una familia general de modelos unidireccionales débilmente no lineales de la forma:
$$u_t + Lu + (u^2)_x = 0$$
donde $L$ es un multiplicador de Fourier con símbolo puramente imaginario $i\omega(k)$, y $\omega(k)$ es la relación de dispersión lineal del modelo.

Los autores investigan cómo las perturbaciones espectrales inestables (autovalores con parte real positiva) surgen a partir de colisiones de autovalores en el estado plano ($\varepsilon = 0$). Se identifican dos tipos principales de inestabilidades:

1. Inestabilidades de Alta Frecuencia: Surgen de colisiones de autovalores no nulos en un exponente de Floquet $p_0 \neq 0$. Estas generan "isolas" de inestabilidad (curvas cerradas en el plano complejo) conocidas como isolas.
2. Inestabilidad de Benjamin-Feir (Modulacional): Surge de la colisión de autovalores en el origen ($\lambda = 0, p_0 = 0$), generando una estructura en forma de "ocho" (lemniscata).

El objetivo es organizar estos datos espectrales en "hojas" (sheets) bidimensionales en el espacio $(\text{Re}\lambda, \text{Im}\lambda, \varepsilon)$ y proporcionar aproximaciones analíticas que abarquen modelos con relaciones de dispersión tanto suaves como discontinuas (como la ecuación de Akers-Milewski).

2. Metodología

Los autores desarrollan un método de perturbación unificado que expande los datos espectrales (autovalores, autofunciones y exponente de Floquet) en series de potencias de la amplitud de la onda $\varepsilon$.

• Parametrización Innovadora: A diferencia de enfoques anteriores donde el exponente de Floquet se expande y sus correcciones son parámetros libres, los autores eligen el ratio $\beta$ de las amplitudes de las autofunciones del estado plano que colisionan como el parámetro libre.

  • Esto simplifica las condiciones de solvabilidad a sistemas lineales en cada orden de la perturbación.
  • Variar $\beta$ y $\varepsilon$ genera una hoja de datos espectrales 2D.
  • Las rebanadas de esta hoja a amplitud fija revelan las curvas de inestabilidad (isolas o lemniscatas).

• Tratamiento de Discontinuidades: Se adapta el marco perturbativo para manejar relaciones de dispersión $\omega(k)$ que no son diferenciables (específicamente con saltos en $k=0$), lo cual es crucial para la ecuación de Akers-Milewski.

• Validación Numérica: Las predicciones asintóticas se comparan rigurosamente con espectros numéricos calculados mediante dos métodos independientes:

  1. Método de Hill de Fourier-Floquet.
  2. Método de continuación cuasi-Newton para pares propios individuales.

3. Contribuciones Clave

1. Método de Perturbación Directo y Unificado: Se presenta un método que trata simultáneamente las inestabilidades de alta frecuencia y la de Benjamin-Feir para la familia general de modelos (1). El método no solo proporciona el comportamiento líder de los autovalores inestables, sino también sus autofunciones asociadas y las tasas de crecimiento más inestables.
2. Estructura de Resonancias y "Hojas" Espectrales:
   • Se demuestra que la estructura de las hojas espectrales está organizada por la interacción resonante entre la onda portadora y los modos colisionantes.
   • Resonancias Triadas ($|k_1 - k_2| = 1$): Generan isolas elípticas centradas en el autovalor del estado plano, con ejes mayores que escalan como $O(\varepsilon)$.
   • Resonancias Cuartetas ($|k_1 - k_2| = 2$): Generan isolas elípticas cuyos centros se desplazan $O(\varepsilon^2)$ del autovalor del estado plano, con ejes mayores que también escalan como $O(\varepsilon^2)$.
   • Se establece que la condición de firma de Krein es necesaria y suficiente para el inicio de estas inestabilidades de alta frecuencia.
3. Comparación Directa de Tasas de Crecimiento: Al aplicar el mismo marco a múltiples modelos (Whitham, Capillary-Whitham, Kawahara, Akers-Milewski), los autores cuantifican cómo compiten las inestabilidades de alta frecuencia con la de Benjamin-Feir. Se determina cuándo una domina a la otra según la amplitud y el número de onda.
4. Análisis de Modelos con Dispersión Discontinua: Se logra, por primera vez, una aproximación analítica del espectro de Benjamin-Feir para la ecuación de Akers-Milewski, cuya relación de dispersión tiene una discontinuidad en $k=0$. Esto demuestra la robustez del método más allá de los modelos de agua clásicos con dispersión suave.

4. Resultados Principales

• Geometría de las Inestabilidades:
  • Las inestabilidades de alta frecuencia (triadas) aparecen como elipses centradas en el eje imaginario con tamaño $O(\varepsilon)$.
  • Las inestabilidades de cuartetas aparecen como elipses desplazadas $O(\varepsilon^2)$ con tamaño $O(\varepsilon^2)$.
  • La inestabilidad de Benjamin-Feir se aproxima mediante una lemniscata (forma de ocho) centrada en el origen.
• Dominancia de Inestabilidades:
  • Si existen inestabilidades de triadas, estas dominan a la de Benjamin-Feir a amplitudes pequeñas porque crecen como $O(\varepsilon)$ frente a $O(\varepsilon^2)$.
  • En ausencia de triadas, la competencia es entre cuartetas ($O(\varepsilon^2)$) y Benjamin-Feir ($O(\varepsilon^2)$). La dominancia depende de coeficientes específicos del modelo.
  • Si no hay ni triadas ni cuartetas, Benjamin-Feir es la fuente dominante de inestabilidad.
• Validación: Existe un acuerdo asintótico correcto entre las predicciones analíticas y los resultados numéricos para modelos representativos (Whitham, Kawahara, Akers-Milewski). Se observan desviaciones cuantitativas significativas en ciertos casos (como el modelo Whitham para cuartetas) donde los términos de orden superior juegan un papel importante, aunque el exponente de error confirma la validez del orden líder.
• Caso Akers-Milewski: Se derivó exitosamente la lemniscata de Benjamin-Feir para este modelo, requiriendo un análisis hasta el orden $O(\varepsilon^3)$ debido a la cancelación de términos en los órdenes inferiores causada por la discontinuidad de la dispersión.

5. Significado e Impacto

Este trabajo es fundamental porque:

1. Unificación: Proporciona un marco teórico único que conecta fenómenos de inestabilidad previamente estudiados por separado (alta frecuencia y modulacional) bajo una misma metodología perturbativa.
2. Generalidad: Extiende la teoría de estabilidad de ondas de Stokes a modelos con relaciones de dispersión no suaves, abriendo la puerta al análisis riguroso de ecuaciones más complejas y físicas (como la de Akers-Milewski).
3. Herramienta Predictiva: La organización de los datos en "hojas espectrales" y la parametrización mediante $\beta$ ofrecen una herramienta poderosa para predecir y visualizar la estructura de inestabilidad en diversos sistemas de ondas dispersivas.
4. Base para Futuras Investigaciones: La estructura de bifurcación revelada sirve como base para futuros resultados de existencia rigurosa de espectros inestables y para estudios de dinámica no lineal cerca del umbral de inestabilidad.

En resumen, el artículo ofrece una descripción analítica profunda y validada numéricamente de cómo surgen las inestabilidades en ondas de Stokes débiles, resolviendo problemas abiertos sobre la comparación de tasas de crecimiento y la aplicabilidad a modelos con discontinuidades.