Non-uniqueness of smooth solutions of the Navier-Stokes equations from almost the same initial conditions

Mediante simulaciones numéricas limpias, este trabajo presenta evidencias de que las ecuaciones de Navier-Stokes pueden admitir soluciones globales distintas a partir de condiciones iniciales casi idénticas, con diferencias tan pequeñas como $10^{-40}$, lo cual aporta perspectivas sobre el problema de la unicidad y existencia asociado al Premio del Milenio.

Lo esencial

¡Claro que sí! Imagina que este artículo es como una historia sobre dos gemelos idénticos que, por un detalle minúsculo, terminan viviendo vidas completamente diferentes.

Aquí tienes la explicación de la investigación de Liao y Qin, traducida a un lenguaje sencillo y con analogías para que cualquiera pueda entenderla:

🌪️ El Gran Misterio: ¿El caos es predecible?

Imagina que el agua que sale de tu grifo o el viento que mueve las nubes sigue unas reglas matemáticas muy estrictas llamadas Ecuaciones de Navier-Stokes. Durante décadas, los matemáticos han creído que si conoces el estado exacto del agua al principio (la "condición inicial"), puedes predecir exactamente cómo se moverá después. Es como decir: "Si suelto esta gota de tinta aquí, sé exactamente dónde estará en un minuto".

Pero, ¿qué pasa si el sistema es caótico? Aquí es donde entra la idea del "efecto mariposa": un aleteo de mariposa hoy podría causar un tornado mañana.

🦋 El Experimento: Dos mundos casi idénticos

Los autores de este papel hicieron un experimento mental (y numérico) muy curioso:

1. Prepararon dos escenarios: Imagina dos tazas de café humeante.

   • Taza A: Tiene un remolino perfecto y simétrico.
   • Taza B: Es exactamente igual a la Taza A, pero con una diferencia tan pequeña que es casi invisible. Imagina que en la Taza B pusimos una gota de agua que es 100 billones de veces más pequeña que un grano de arena (¡una diferencia de $10^{-40}$!).

2. La pregunta: Si las dos tazas son casi idénticas, ¿se comportarán igual mientras se enfrían?

3. El resultado sorprendente:

   • Al principio, las dos tazas parecen iguales.
   • Pero, después de un tiempo, se separan por completo.
   • La Taza A mantiene su simetría perfecta (gira igual en ambos lados).
   • La Taza B se vuelve un caos total, pierde su simetría y su comportamiento estadístico es totalmente diferente.

🛠️ ¿Cómo lo hicieron? (El "Microscopio" Matemático)

Aquí está la parte genial. Normalmente, cuando los científicos usan ordenadores para simular esto, el ordenador comete pequeños errores de redondeo (como si alguien susurrara un ruido de fondo). En sistemas caóticos, esos ruidos pequeños crecen y arruinan la simulación, haciendo que parezca que el sistema es impredecible por culpa del ordenador, no de la física.

Para evitar esto, los autores usaron una técnica llamada Simulación Numérica Limpia (CNS).

• La analogía: Imagina que quieres escuchar un susurro en una habitación ruidosa. En lugar de gritar para que te escuchen (como hacen los métodos normales), ellos construyeron una habitación insonorizada y usaron un micrófono de ultra-alta precisión.
• Gracias a esto, pudieron eliminar el "ruido" del ordenador. Lo que vieron no fue un error de cálculo, sino el comportamiento real de las ecuaciones.

💡 ¿Qué significa esto para el mundo?

El descubrimiento es profundo y un poco inquietante:

• No unicidad: Demuestran que, incluso con condiciones iniciales casi idénticas (casi indistinguibles), las ecuaciones de Navier-Stokes pueden dar dos soluciones globales totalmente diferentes.
• El problema del Milenio: Esto toca uno de los problemas más famosos de matemáticas (el Premio del Milenio del Instituto Clay). Si las ecuaciones pueden tener dos respuestas diferentes para el mismo punto de partida, eso desafía nuestra comprensión de la existencia y unicidad de las soluciones suaves.
• La metáfora final: Es como si tuvieras dos semillas genéticamente idénticas plantadas en el mismo suelo. Una crece como un roble perfecto y la otra, debido a una diferencia de un átomo en su ADN, crece como un helecho retorcido.

En resumen

Este papel nos dice que el caos no es solo "ruido" o "errores de cálculo". Es una propiedad fundamental. Un cambio tan pequeño que ni siquiera podrías medirlo con el instrumento más preciso del universo puede, con el tiempo, destruir la simetría y cambiar por completo el destino de un sistema fluido.

Los autores no dicen que "todo es imposible de predecir", sino que nos dan un mapa para entender que la realidad fluida es mucho más sensible y compleja de lo que pensábamos, y que quizás, las ecuaciones que rigen el universo admiten múltiples caminos para un mismo punto de partida.

Resumen técnico

A continuación presento un resumen técnico detallado del artículo en español, estructurado según los puntos solicitados:

Título: No unicidad de soluciones suaves de las ecuaciones de Navier-Stokes a partir de condiciones iniciales casi idénticas

Autores: Shijun Liao y Shijie Qin
Institución: Universidad de Jiaotong de Shanghái, China.

---

1. El Problema

El artículo aborda uno de los Problemas del Milenio del Instituto Clay: la existencia y unicidad de soluciones suaves para las ecuaciones de Navier-Stokes (NS) en tres dimensiones (y su análogo en 2D en este contexto).

• Contexto: Aunque se acepta que los flujos turbulentos se describen mediante las ecuaciones NS, la unicidad de sus soluciones suaves con energía finita sigue sin estar demostrada.
• Limitación de métodos anteriores: Recientemente, Coiculescu y Palasek demostraron la no unicidad en el espacio crítico $BMO^{-1}$, pero su construcción utiliza datos con energía infinita, lo que limita su aplicabilidad física directa.
• Desafío numérico: Tradicionalmente, es imposible demostrar la no unicidad mediante simulaciones numéricas estándar (como la Simulación Numérica Directa o DNS). Esto se debe al "efecto mariposa" en sistemas caóticos: el ruido numérico (errores de redondeo y truncamiento), aunque minúsculo, crece exponencialmente, contaminando rápidamente la trayectoria exacta y haciendo que los resultados no sean confiables para distinguir entre soluciones matemáticas reales y artefactos numéricos.

2. Metodología: Simulación Numérica Limpia (CNS)

Para superar las limitaciones del ruido numérico, los autores utilizan la Simulación Numérica Limpia (Clean Numerical Simulation - CNS).

• Principio de CNS: A diferencia de la DNS, la CNS reduce tanto el error de truncamiento como el error de redondeo a niveles extremadamente bajos mediante el uso de aritmética de alta precisión (múltiples dígitos significativos) y expansiones de Taylor de alto orden.
• Configuración del modelo:
  • Ecuación: Ecuación de Navier-Stokes incompresible 2D en un dominio cuadrado $[0, 2\pi] \times [0, 2\pi]$ con condiciones de frontera periódicas, bajo forzamiento de Kolmogorov.
  • Parámetros: Número de Reynolds $Re = 2000$ y escala de forzamiento $n_K = 16$.
  • Método numérico: Método pseudo-espectral de Fourier para el espacio y expansión de Taylor de orden $M$ para el tiempo, con paso de tiempo $\Delta t = 10^{-3}$.
• Condiciones Iniciales: Se comparan cuatro casos con condiciones iniciales casi idénticas:
  1. Flow CNS: $\psi_1(x, y, 0) = -\frac{1}{2}[\cos(x+y) + \cos(x-y)]$. (Posee simetría de rotación y traslación).
  2. Flow CNS'$_1$, CNS'$_2$, CNS'$_3$: $\psi_2(x, y, 0) = \psi_1(x, y, 0) + \delta \sin(x+y)$.
     • Se prueban tres valores de perturbación $\delta$: $10^{-10}$, $10^{-20}$ y $10^{-40}$.
     • La diferencia entre las condiciones iniciales es del orden de magnitud de $10^{-40}$.
• Validación: Se utiliza una estrategia de doble verificación donde se ejecuta una segunda simulación con parámetros aún más estrictos (mayor orden $M$ y más dígitos significativos $N_s$) para confirmar que el ruido numérico es despreciable en el intervalo de tiempo estudiado ($t \in [0, 300]$).

3. Contribuciones Clave

• Evidencia Numérica de No Unicidad: Proporcionan los primeros ejemplos numéricos explícitos donde las ecuaciones de Navier-Stokes admiten soluciones globales distintas a partir de condiciones iniciales que difieren en magnitudes tan pequeñas como $10^{-40}$, manteniendo una energía finita.
• Superación del Ruido Numérico: Demuestran que, con la metodología CNS, es posible aislar el comportamiento dinámico real del sistema caótico del ruido artificial, permitiendo observar la divergencia de trayectorias que serían indistinguibles con métodos tradicionales.
• Conjetura Formal: Proponen una conjetura basada en sus hallazgos: Las ecuaciones de Navier-Stokes admiten soluciones globales suaves distintas para todo $t > 0$ a partir de condiciones iniciales $U_1$ y $U_2$ tales que $|U_1 - U_2| \leq \delta$ para cualquier $\delta$ arbitrariamente pequeño (o cuando $\delta \to 0$).

4. Resultados Principales

Los resultados se presentan en el intervalo de tiempo $t \in [0, 300]$:

• Divergencia de Trayectorias: Inicialmente, las soluciones de los casos perturbados (CNS'$_1$, CNS'$_2$, CNS'$_3$) son indistinguibles de la solución base (CNS). Sin embargo, tras un periodo de tiempo, las trayectorias se separan drásticamente debido a la inestabilidad caótica (cascada de expansión de ruido).
• Diferencia en Disipación de Energía:
  • Para $t > 100$, la tasa de disipación de energía cinética espacialmente promediada $\langle D \rangle_A$ del flujo base (CNS) es casi el doble que la de los flujos perturbados.
  • Las soluciones perturbadas convergen a un estado estadístico diferente, aunque entre sí (CNS'$_1$, CNS'$_2$, CNS'$_3$) son estadísticamente idénticas, lo que sugiere que la magnitud de $\delta$ no altera el estado final, solo la sensibilidad inicial.
• Pérdida de Simetría Espacial:
  • La condición inicial $\psi_1$ posee simetría de rotación y traslación.
  • La condición inicial $\psi_2$ (perturbada) solo posee simetría de traslación.
  • Aunque la perturbación es minúscula, el caos amplifica esta diferencia, destruyendo la simetría de rotación en las soluciones perturbadas. En $t=200$, los campos de vorticidad muestran patrones espaciales cualitativamente diferentes entre el caso base y los perturbados.
• Estadísticas Distintas: Las funciones de densidad de probabilidad (PDF) de la disipación de energía y los promedios espacio-temporales difieren significativamente entre la solución base y las soluciones perturbadas.

5. Significado e Impacto

• Implicaciones Matemáticas: Estos resultados actúan como contraejemplos potenciales para la unicidad de soluciones suaves de las ecuaciones de Navier-Stokes con energía finita. Sugieren que la unicidad podría no sostenerse en el sentido estricto para condiciones iniciales arbitrariamente cercanas.
• Guía para Matemáticos: El artículo ofrece una "hoja de ruta" o inspiración para que los matemáticos intenten probar rigurosamente la no unicidad, utilizando la idea de que pequeñas perturbaciones pueden destruir simetrías y llevar a estados globales distintos.
• Validación de Simulaciones: Refuerza la necesidad de métodos de alta precisión (como CNS) en el estudio de la turbulencia, demostrando que las simulaciones tradicionales pueden estar "contaminadas" por ruido numérico que enmascara o altera la física real del sistema.

En resumen, el trabajo utiliza una metodología computacional de ultra-alta precisión para demostrar que, en el régimen turbulento, la sensibilidad extrema a las condiciones iniciales puede llevar a soluciones globalmente distintas a partir de diferencias casi imperceptibles, desafiando la noción convencional de unicidad en las ecuaciones de Navier-Stokes.