Nonuniqueness of Leray-Hopf solutions to the unforced incompressible 3D Navier-Stokes Equation

Aquí presento un resumen técnico detallado del artículo "Nonuniqueness of Leray–Hopf solutions to the unforced incompressible 3D Navier–Stokes Equation" (No unicidad de soluciones de Leray-Hopf para la ecuación de Navier-Stokes incompresible no forzada en 3D), escrito por Thomas Hou, Yixuan Wang y Changhe Yang.

1. El Problema

El artículo aborda uno de los problemas abiertos más centrales en la dinámica de fluidos matemática: la unicidad de las soluciones de Leray-Hopf para las ecuaciones de Navier-Stokes incompresibles en tres dimensiones (3D) sin fuerzas externas.

• Contexto: Jean Leray (1934) y Eberhard Hopf (1951) demostraron la existencia de soluciones débiles globales (soluciones de Leray-Hopf) que satisfacen desigualdades de energía. Sin embargo, la unicidad de estas soluciones sigue siendo un problema abierto (parte del problema del Milenio).
• Objetivo: Demostrar rigurosamente que, bajo ciertas condiciones iniciales, existen infinitas soluciones de Leray-Hopf distintas que evolucionan desde la misma condición inicial.
• Marco de trabajo: Los autores se centran en soluciones auto-similares (invariantes de escala) que surgen de datos iniciales singulares, específicamente aquellos con una singularidad de tipo $O(|x|^{-1})$ en el origen.

2. Metodología: Prueba Asistida por Computadora Rigurosa

La contribución principal del trabajo es el desarrollo de un marco de prueba asistida por computadora rigurosa para verificar la existencia de un perfil auto-similar inestable. A diferencia de trabajos anteriores que dependían de suposiciones espectrales no verificadas o forzaron términos externos, este trabajo resuelve la ecuación exacta sin forzamiento.

La estrategia se divide en tres etapas principales:

A. Construcción de un Perfil Aproximado de Alta Precisión

1. Método Numérico: Utilizan un método de elementos finitos en un dominio acotado para obtener una solución aproximada de la ecuación de Navier-Stokes auto-similar.
2. Optimización: Emplean un descenso de gradiente proyectado para minimizar el valor propio más negativo del operador linealizado, buscando un perfil que tenga un valor propio inestable (parte real positiva en el operador linealizado, o negativo en la formulación de autovalores utilizada en el papel).
3. Interpolación Espectral: Para garantizar la divergencia nula exacta y la suavidad, interpolan la solución de elementos finitos en una base espectral trigonométrica en coordenadas esféricas transformadas ($r = \tan \beta$). Esto permite manejar la condición de frontera en el infinito y la singularidad en el origen de manera analítica.

B. Verificación Rigurosa de la Existencia (Análisis de Perturbación)

Una vez obtenido un candidato numérico $(U, v, \lambda)$ con residuos pequeños, el objetivo es demostrar que existe una solución exacta $(\tilde{U}, \tilde{v}, \tilde{\lambda})$ en su vecindad.

1. Descomposición del Operador: Descomponen el operador linealizado $L$ en una parte coerciva ($L_1$) y una perturbación compacta ($L_2$). La viscosidad proporciona la coercividad necesaria en la norma $L^2/H^1$.
2. Aproximación de Rango Finito: La parte compacta $L_2$ se aproxima mediante un operador de rango finito ($L_{ap}$) construido sobre un subespacio de dimensión pequeña (basado en los modos de mayor energía de $L_2$).
3. Argumento de Punto Fijo: Utilizan el Teorema del Punto Fijo de Schauder. Demuestran que, si el residuo de la solución aproximada es suficientemente pequeño y el operador linealizado restringido al subespacio de rango finito es invertible (verificado numéricamente con aritmética de intervalos), entonces existe una solución exacta en una bola pequeña alrededor de la aproximación.
4. Cálculos Rigurosos: Todos los errores (residuos, errores de cuadratura, cotas $L^\infty$ de derivadas) se acotan rigurosamente utilizando aritmética de intervalos y desigualdades de Bernstein para polinomios trigonométricos.

C. Demostración de la No Unicidad

Una vez establecida la existencia de un perfil auto-similar $\tilde{U}$ con un valor propio inestable $\tilde{\lambda}$:

1. Localización: Modifican el dato inicial singular para hacerlo de soporte compacto (usando una descomposición de Bogovskii y un corte suave).
2. Descomposición Estable-Inestable: Analizan la evolución de la perturbación alrededor del perfil auto-similar. Utilizan una proyección de Riesz para separar la dinámica en una parte estable y una parte inestable.
3. Bifurcación: Demuestran que la existencia del modo inestable permite construir infinitas trayectorias distintas que coinciden en el tiempo inicial (tras la localización) pero divergen posteriormente. Esto genera infinitas soluciones de Leray-Hopf distintas para el mismo dato inicial.

3. Contribuciones Clave

1. Primera Prueba Rigurosa: Es, según los autores, la primera prueba rigurosa asistida por computadora que demuestra la no unicidad de soluciones de Leray-Hopf para Navier-Stokes no forzada en 3D.
2. Eliminación de Suposiciones Espectrales: A diferencia de trabajos previos (como los de Jia-Sverak o Albritton-Brue-Colombo) que asumían la existencia de perfiles inestables o añadían fuerzas externas, este trabajo calcula y verifica la existencia del perfil inestable exacto que satisface la ecuación sin forzamiento.
3. Marco Computacional Innovador: Desarrollan un método que combina elementos finitos, bases espectales divergencia-libres, y verificación de rango finito con aritmética de intervalos. La viscosidad juega un papel crucial al proporcionar coercividad, lo que simplifica el análisis en comparación con el caso de Euler (sin viscosidad).
4. Código Abierto: El código utilizado para la verificación es público, permitiendo la reproducibilidad total de los resultados.

4. Resultados Principales

• Teorema 1 (No Unicidad): Se demuestra la existencia de infinitas soluciones de Leray-Hopf distintas y adecuadas para la ecuación de Navier-Stokes en $\mathbb{R}^3 \times [0, 1]$ con la misma condición inicial de soporte compacto y divergencia nula. Estas soluciones son suaves para tiempos positivos pero fallan en cumplir la condición Prodi-Serrin que garantiza la unicidad.
• Teorema 2 (Existencia del Perfil y Autovalor): Se prueba la existencia de una solución débil $(\tilde{U}, \tilde{v}, \tilde{\lambda})$ al sistema de ecuaciones auto-similares donde $\tilde{U}$ es el perfil y $\tilde{\lambda} < 0$ (en la formulación del papel) corresponde a un valor propio inestable.
• Constantes Verificadas: Los autores reportan constantes numéricas verificadas rigurosamente, como un valor propio $\tilde{\lambda} \approx -0.113$ y residuos del orden de $10^{-7}$, suficientes para cerrar el argumento de punto fijo.

5. Significado e Impacto

Este trabajo representa un hito significativo en la teoría de ecuaciones en derivadas parciales y la dinámica de fluidos:

• Resolución de un Problema Abierto: Proporciona evidencia contundente de que la unicidad de las soluciones de Leray-Hopf falla en 3D, incluso para datos iniciales suaves (tras localización), desafiando la intuición clásica de que la viscosidad debería regularizar y asegurar la unicidad.
• Nueva Metodología: Establece un nuevo estándar para las pruebas asistidas por computadora en problemas no lineales complejos, demostrando que es posible verificar la existencia de soluciones inestables y estructuras de bifurcación con rigor matemático completo.
• Implicaciones Físicas: Sugiere que la turbulencia y la falta de unicidad pueden ser fenómenos intrínsecos a las ecuaciones de Navier-Stokes, no solo artefactos de soluciones débiles mal definidas o forzamientos externos.

En resumen, el artículo combina análisis funcional avanzado, teoría espectral y computación de alta precisión para demostrar que las ecuaciones de Navier-Stokes en 3D admiten múltiples evoluciones temporales a partir de un mismo estado inicial, resolviendo una cuestión fundamental que ha permanecido abierta durante décadas.