Unstable Interface Dynamics for Gravity Stokes Flow

Resumen Técnico: Dinámica de Interfaz Inestable para el Flujo de Stokes por Gravedad

Planteamiento del Problema
Este artículo investiga la dinámica de largo tiempo y los mecanismos de inestabilidad del problema de interfaz nítida para el flujo de Stokes por gravedad en dos dimensiones. El sistema modela dos fluidos incompresibles y homogéneos con densidades constantes distintas ($\rho_+$ y $\rho_-$) que evolucionan bajo la ley de Stokes con una fuerza gravitacional. La interfaz $\Gamma(t)$ separa la región superior $\Omega_+(t)$ de la región inferior $\Omega_-(t)$ en un dominio periódicamente horizontal ($\mathbb{T} \times \mathbb{R}$).

El estudio se centra en dos regímenes primarios:

1. Régimen Inestable (Rayleigh-Taylor): El fluido más denso se encuentra sobre el fluido más ligero ($\rho_+ < \rho_-$). Aunque la existencia global de soluciones para datos iniciales suaves es conocida (p. ej., [12]), el comportamiento a largo tiempo de las cantidades geométricas (perímetro, curvatura) permanece mayormente inexplorado.
2. Régimen Estable: El fluido más denso se encuentra debajo del fluido más ligero ($\rho_+ > \rho_-$). El artículo examina si las configuraciones iniciales estables pueden evolucionar hacia estados inestables, investigando específicamente las "inestabilidades de vuelco" (turning instabilities), donde la interfaz deja de ser una gráfica en un tiempo finito.

El sistema se analiza mediante la Ecuación de Dinámica de Contorno (CDE) derivada en [11]:
$$\partial_t z(\alpha, t) = (\rho_- - \rho_+) \int_{\mathbb{T}} S(z(\alpha, t) - z(\beta, t)) \dot{z}^\perp(\beta, t) z^2(\beta, t) d\beta$$
donde $S$ es el Stokeslet $\text{x}_1$-periódico.

Metodología
Los autores emplean una combinación de técnicas analíticas y simulaciones numéricas:

• Análisis de Energía Potencial: Una herramienta central es un funcional de energía potencial modulada $E(t) = \int_{\mathbb{T}\times\mathbb{R}} x_2(\rho_s(x_2) - \rho(x,t)) dx$, donde $\rho_s$ representa el estado estratificado inestable. Los autores demuestran que $E(t)$ es monótonamente creciente y su derivada temporal se relaciona con la norma $L^2$ de un operador no local aplicado a la densidad.
• Explotación de Simetrías: El análisis depende fuertemente de simetrías específicas de los datos iniciales: simetría central ($z(\alpha) = -z(-\alpha)$) y simetría par con respecto a $x_1 = \pm \pi/2$. Estas simetrías se demuestran como conservadas por el flujo y permiten la reducción del dominio a una franja restringida para estimaciones geométricas.
• Dicotomía Geométrica: Para probar resultados de crecimiento, los autores establecen una dicotomía basada en la cantidad $\delta(t) = \|( -\Delta)^{-1} \partial_1 \rho(t) \|_{L^2}^2$. Demuestran que si el perímetro $L(t)$ no crece lo suficientemente rápido, la cantidad $\delta(t)$ debe volverse pequeña, lo que a su vez obliga a la curvatura máxima $K(t)$ a crecer rápidamente.
• Analiticidad y Fingering (Formación de Dedos): Para el análisis de fingering, los autores utilizan la analiticidad de las soluciones en una franja compleja. Aplican un lema sobre la distribución de puntos críticos de funciones analíticas para estimar el número de "dedos" (extremos locales) que pueden desarrollarse.
• Construcción de la Inestabilidad de Vuelco: Para probar la ruptura en tiempo finito en el régimen estable, los autores construyen datos iniciales específicos que satisfacen condiciones de simetría par/impar. Analizan la derivada temporal de la pendiente en el origen para mostrar que una interfaz estable puede desarrollar una tangente vertical y volcarse.
• Simulaciones Numéricas: Los autores implementan un esquema numérico de tipo espectral utilizando diferencias centrales y la regla de Simpson para las integrales singulares, regularizadas mediante viscosidad artificial, para visualizar el crecimiento de la longitud de la interfaz y la curvatura.

Contribuciones Clave y Resultados

1. Crecimiento en Tiempo Infinito en el Régimen Inestable (Teorema 1.1):
   El artículo demuestra que, para datos iniciales inestables que satisfacen condiciones de simetría específicas, la interfaz exhibe un crecimiento ilimitado en al menos una cantidad geométrica (perímetro $L(t)$ o curvatura máxima $K(t)$) cuando $t \to \infty$. Específicamente, para cualquier $\epsilon > 0$:
   $$\limsup_{t \to \infty} t^{-\frac{1}{17} + \epsilon} (L(t) + K(t)) = \infty$$
   Este resultado establece que las soluciones no se estabilizan hacia un estado estacionario, sino que exhiben un crecimiento de largo tiempo genérico, contrastando con el resultado de existencia global que se mantiene independientemente de la estratificación de densidad.

2. Estimación del Fenómeno de Fingering (Teorema 1.3):
   Los autores proporcionan un límite inferior en el crecimiento del número de dedos (puntos críticos de la pendiente de la interfaz) para una clase de soluciones analíticas. Demuestran que el número de dedos $F(t)$ crece al menos como $C(1 + \sqrt{t})$. Esto conecta el crecimiento analítico de la solución con el fenómeno físico de viscous fingering.

3. Inestabilidad de Vuelco en Tiempo Finito en el Régimen Estable (Teorema 1.4):
   Contrario a la intuición de que la estratificación estable previene la inestabilidad, los autores demuestran la existencia de datos iniciales suaves en el régimen estable ($\rho_+ > \rho_-$) que evolucionan hacia una inestabilidad de vuelco en tiempo finito. La interfaz deja de ser una gráfica y la condición de Rayleigh-Taylor se rompe. El artículo construye una familia de tales curvas iniciales y prueba que el vuelco ocurre en puntos simétricos específicos.

4. Validación Numérica:
   Las simulaciones numéricas confirman las predicciones teóricas, mostrando un crecimiento significativo en la longitud de la interfaz y la curvatura en el corto plazo para configuraciones inestables, así como el desarrollo de dedos en la evolución de la interfaz.

Significancia y Reivindicaciones
El artículo afirma abordar un vacío en la comprensión de la dinámica de largo tiempo del sistema de transporte-Stokes. Mientras que la existencia global fue establecida en trabajos previos [12], el comportamiento de estas soluciones globales era desconocido. Los autores argumentan que sus resultados demuestran que:

• La existencia global no implica estabilidad; las configuraciones inestables conducen a un crecimiento geométrico ilimitado.
• Las configuraciones estables no son inmunes a la inestabilidad; datos iniciales suaves específicos pueden evolucionar hacia singularidades de vuelco, rompiendo la propiedad de gráfica.
• El crecimiento de las cantidades geométricas (perímetro y curvatura) es una característica genérica del régimen inestable, apoyando la conexión con los fenómenos de viscous fingering observados en flujos de medios porosos.

Los autores señalan explícitamente que la tasa de crecimiento derivada en el Teorema 1.1 ($t^{1/17}$) probablemente no es exacta, ya que la evidencia numérica sugiere una tasa de crecimiento lineal, pero la prueba analítica proporciona un límite inferior riguroso para la inestabilidad. El trabajo destaca la distinción en la dinámica del flujo de Stokes respecto al problema de Muskat, que es más singular, particularmente en cuanto a la naturaleza del amortiguamiento lineal y los mecanismos específicos de inestabilidad.