Por qué Navier-Stokes 2D es más fácil que 3D

Las ecuaciones de Navier-Stokes describen cómo se mueven los fluidos. Funcionan en 2D (plano, como el agua que se esparce sobre una mesa) y en 3D (la vida real, como las corrientes oceánicas que giran alrededor de un submarino o el viento que pasa por un rascacielos). Mismas ecuaciones. Casi idéntico.

Aquí está el giro. En 2D, los matemáticos pueden probar que las ecuaciones siempre se comportan bien, que las matemáticas nunca fallan y que las soluciones se mantienen suaves por toda la eternidad. ¿En 3D? Nadie lo sabe. Ni una sola persona en la Tierra. El fluido podría hacer algo tan violento y repentino que las matemáticas dejen de funcionar por completo, y demostrar si eso puede suceder es el Problema del Premio del Milenio de Clay, valorado en un millón de dólares.

Esto no es sólo "3D es más difícil porque hay más cosas". Un mecanismo específico en 3D no existe en 2D. Lo cambia todo.

Las ecuaciones incompresibles de Navier-Stokes en $\mathbb{R}^n$ (o un dominio periódico $\mathbb{T}^n$) con $\nu > 0$ son

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0.$$

Para $n = 2$, la existencia global y la unicidad de soluciones clásicas para datos iniciales suaves y libres de divergencia es un teorema. Las referencias clave son Ladyzhenskaya (1959), basada en trabajos anteriores de Leray (1934). El resultado se extiende a soluciones suaves en $\mathbb{R}^2$, $\mathbb{T}^2$ y dominios acotados con condiciones de contorno estándar.

Para $n = 3$, la existencia global de soluciones clásicas a partir de datos suaves arbitrarios es abierta. Este es el contenido del Problema del Milenio de Clay formulado por Fefferman (2000). Leray (1934) estableció la existencia global de soluciones débiles en $L^2(\mathbb{R}^3)$, pero la unicidad y regularidad de estas soluciones siguen sin resolverse.

La brecha entre $n = 2$ y $n = 3$ no es contabilidad. Refleja una diferencia estructural en la ecuación de vorticidad, las propiedades de escala de la no linealidad y las estimaciones a priori disponibles. Cada uno de estos se examina en las secciones siguientes.

La pregunta del millón sólo se refiere al 3D. ¿Por qué? Porque el 2D ya está hecho. Finalizado. Los matemáticos demostraron hace décadas que las soluciones bidimensionales de Navier-Stokes siempre permanecen suaves, sin importar las condiciones iniciales que se les presenten, sin importar cuánto tiempo se espere. No se necesita premio por un problema resuelto.

Así que la verdadera pregunta no es "¿por qué es difícil el 3D?" Es "¿por qué el 2D es fácil y el 3D difícil?" ¿Qué se rompe exactamente cuando agregas esa tercera dimensión?

El problema de Clay (Fefferman 2000) considera el problema de Cauchy para Navier-Stokes incompresible en $\mathbb{R}^3$ con viscosidad $\nu > 0$ y datos iniciales suaves y libres de divergencia $u_0 \in C^\infty(\mathbb{R}^3)$ que satisfacen las condiciones de decaimiento adecuados. La pregunta: ¿existe una solución suave $u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$ con energía acotada?

La afirmación análoga para $\mathbb{R}^2$ es un teorema. Ladyzhenskaya demostró la existencia global y la unicidad de soluciones fuertes para 2D Navier-Stokes con $u_0 \in H^1(\mathbb{R}^2)$, y el arranque proporciona $C^\infty$ regularidad para obtener datos suaves. La prueba se basa en estimaciones a priori específicas de dos dimensiones que no se extienden a tres.

El Problema del Milenio es, por tanto, puramente tridimensional. Los resultados parciales en 3D (las soluciones débiles de Leray, 1934; el teorema de regularidad parcial de Caffarelli-Kohn-Nirenberg, 1982; varios criterios de regularidad condicional) no llegan a resolver la cuestión de la regularidad global completa. Cada resultado parcial resalta una brecha específica en nuestro control de las soluciones 3D.

2D tiene un arma secreta. Se llama vorticidad: cuánto gira el fluido en cada punto.

En 2D, la vorticidad es solo un número. Eso es todo. En sentido horario o antihorario, rápido o lento. Y esto es lo que hace que dos dimensiones sean tan notablemente diferentes de tres: estos pequeños remolinos pueden flotar a través del fluido y desvanecerse gradualmente debido a la fricción, pero nunca, bajo ninguna circunstancia, pueden volverse más fuertes de lo que eran al principio. ¿Giro máximo en el tiempo cero? Ese es el giro máximo que jamás verás.

¿Por qué es importante? Todo se deriva de ello. La velocidad se mantiene suave. La presión se mantiene suave. La solución sigue funcionando para siempre, sin importar lo absurdamente lejos que vayamos en el futuro, porque esa única restricción de vorticidad actúa como la primera ficha de dominó de una cadena que derriba todas las demás fichas a la vista.

En 2D, la vorticidad $\omega = \partial_1 u_2 - \partial_2 u_1$ es un escalar que satisface

$$\partial_t \omega + (u \cdot \nabla)\omega = \nu \Delta \omega.$$

Esta es una ecuación de advección-difusión para la vorticidad escalar $\omega$. El sistema completo sigue siendo no lineal porque $u$ se recupera de $\omega$ a través de la ley de Biot-Savart $u = \nabla^\perp (-\Delta)^{-1} \omega$, pero la ecuación no tiene un término fuente que se estire en vórtice: el lado derecho contiene solo el término laplaciano, no el término $(\omega \cdot \nabla)u$ que aparece en 3D.

El máximo escalar El principio se aplica directamente: $\|\omega(t)\|_{L^\infty} \leq \|\omega_0\|_{L^\infty}$ para todo $t \geq 0$. Simultáneamente, las normas $L^p$ de $\omega$ no son crecientes para todo $1 \leq p \leq \infty$.

A partir del límite $L^\infty$ en $\omega$, se recupera $u \in L^\infty(0,T; W^{1,p})$ para todo $p < \infty$ mediante estimaciones de Calderón-Zygmund sobre el Núcleo de Biot-Savart. Se obtiene una mayor regularidad al diferenciar la ecuación de vorticidad y aplicar bootstrap parabólico: cada derivada espacial de $\omega$ satisface una ecuación parabólica con coeficientes controlables, por lo que los límites se propagan a todos los órdenes.

El principio de máxima vorticidad 2D tiene un pedigrí más antiguo para el caso no viscoso. Wolibner (1933) demostró la existencia global de Euler 2D en espacios de Hölder, y Yudovich (1963) estableció la unicidad de las soluciones de Euler 2D de vorticidad limitada. Con viscosidad ($\nu > 0$), el suavizado parabólico solo fortalece estas estimaciones. La prueba de Ladyzhenskaya de la regularidad global de Navier-Stokes 2D se basa en esta estructura, combinada con la desigualdad de interpolación de Ladyzhenskaya $\|f\|_{L^4}^2 \leq C \|f\|_{L^2} \|\nabla f\|_{L^2}$ (válida en 2D, con una estructura de exponente diferente a la de su contraparte 3D).

En 3D, la vorticidad no es un número. Es un vector que lleva tanto una dirección como una fuerza, y deberías imaginártelo como pequeños tubos de tornado atravesando el fluido.

Esto es lo que arruina todo. Esos tubos se pueden estirar. Tira de uno como si fuera un caramelo: se adelgaza y gira más rápido. Mucho, mucho más rápido. Esto es estiramiento de vórtice, y es el villano de toda la historia porque significa que el fluido puede amplificar su propia rotación, alimentando energía en escalas cada vez más pequeñas hasta que, posiblemente, la rotación en un solo punto se vuelve infinitamente intensa.

Eso es una explosión. Las matemáticas se rompen.

¿Puede la viscosidad (la fricción interna del fluido) siempre frenar bruscamente el estiramiento antes de que llegue al infinito, o el estiramiento a veces supera la fricción y gana? Nadie lo sabe. Ésa es, literalmente, la pregunta del millón. Este tira y afloja entre el estiramiento y la fricción es la razón por la que el problema es tan difícil.

En 3D, la vorticidad $\omega = \nabla \times u$ satisface

$$\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega.$$

El término $(\omega \cdot \nabla)u$ es el término de estiramiento de vórtice, ausente en 2D. Es cuadrático en el sentido de que $u$ se recupera a partir de $\omega$ mediante la ley 3D de Biot-Savart, así que $(\omega \cdot \nabla)u$ escala aproximadamente como $|\omega|^2$ en el peor caso. Este término permite el crecimiento de $\|\omega\|_{L^\infty}$ y destruye el principio máximo escalar disponible en 2D.

Para Euler 3D, el criterio de Beale-Kato-Majda (1984) establece que una solución suave en $[0, T)$ se prolonga más allá del tiempo $T$ si y solo si

$$\int_0^T \|\omega(t)\|_{L^\infty} \, dt < \infty.$$

Criterios de continuación análogos valen para Navier-Stokes 3D. La explosión, si ocurre, exige que $\|\omega\|_{L^\infty}$ se vuelva no integrable en el tiempo. El término de estiramiento es la fuente en la ecuación de vorticidad que puede amplificar la vorticidad y bloquear un argumento de principio máximo. En 2D, $\|\omega\|_{L^\infty}$ queda acotada por los datos iniciales para todo tiempo; en 3D, controlar esa norma es el problema abierto central.

Progreso parcial: Caffarelli-Kohn-Nirenberg (1982) demostraron que, para cualquier solución débil adecuada de Navier-Stokes 3D, el conjunto de puntos singulares en el espacio-tiempo tiene medida de Hausdorff parabólica unidimensional nula. Las singularidades, si existen, son extremadamente escasas. Pero el teorema no descarta su existencia.

Para el caso no viscoso, Elgindi (2021) demostró formación de singularidades en tiempo finito para Euler 3D con datos iniciales $C^{1,\alpha}$ ($\alpha$ pequeño), usando un mecanismo impulsado por el estiramiento de vórtice a lo largo de un eje de simetría. Esto no implica directamente explosión para Navier-Stokes (la viscosidad todavía podría regularizar), pero sí muestra que el mecanismo de estiramiento es lo bastante fuerte como para producir singularidades sin amortiguamiento viscoso.

El estiramiento de vórtice no es el único problema. Hay una razón estructural más profunda por la que 3D se resiste a la prueba, y aparece cuando haces "zoom" sobre el fluido.

Las ecuaciones de Navier-Stokes tienen un truco de acercamiento. Toma cualquier solución, haz zoom en una región más pequeña, acelera el tiempo en la cantidad correcta y obtendrás otra solución perfectamente válida. Entonces, ¿qué le pasa a la estimación de energía cuando haces zoom?

• En 2D, el zoom mantiene la energía en la misma escala. Los matemáticos llaman a esto escalamiento crítico. Tus estimaciones de energía funcionan en todas las escalas. Grande o pequeño, no pierdes el control.
• En 3D, la cota de energía se vuelve más débil como control de pequeña escala. Eso es escalamiento supercrítico, y es devastador: las estimaciones conocidas pierden agarre precisamente en las escalas donde una explosión se concentraría.

Una analogía: en 2D tu linterna siempre da suficiente luz. En 3D, cuanto más te acercas, más tenue se vuelve, y el fluido se hace más difícil de resolver. Te quedas a oscuras.

Esto no es un pequeño inconveniente técnico que un truco ingenioso vaya a arreglar. Es una pared. Las herramientas matemáticas estándar no controlan Navier-Stokes 3D a pequeña escala. Hace falta algo fundamentalmente nuevo.

Las ecuaciones de Navier-Stokes son invariantes bajo el escalamiento

$$u(x,t) \mapsto \lambda u(\lambda x, \lambda^2 t), \qquad p(x,t) \mapsto \lambda^2 p(\lambda x, \lambda^2 t),$$

para cualquier $\lambda > 0$. La norma $L^2$ se transforma como $\|u_\lambda\|_{L^2(\mathbb{R}^n)} = \lambda^{1-n/2} \|u\|_{L^2}$.

• Para $n = 2$: $\|u_\lambda\|_{L^2} = \|u\|_{L^2}$. Es invariante de escala. La ecuación es crítica en energía.
• Para $n = 3$: $\|u_\lambda\|_{L^2} = \lambda^{-1/2} \|u\|_{L^2}$, que disminuye cuando $\lambda \to \infty$ (al acercar la escala). La ecuación es supercrítica en energía: una norma de energía acotada se vuelve más débil como control de pequeña escala.

El espacio crítico de Sobolev para Navier-Stokes 3D es $\dot{H}^{1/2}(\mathbb{R}^3)$, invariante bajo el escalamiento natural. Pero la identidad de energía solo controla $u$ en $L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$. Eso está media derivada por debajo del nivel crítico. Esta es la brecha de supercriticidad.

En 2D, la identidad de energía $\frac{d}{dt}\|u\|_{L^2}^2 = -2\nu\|\nabla u\|_{L^2}^2$ da exactamente el control al nivel crítico que hace falta; combinada con el principio máximo para la vorticidad, proporciona suficiente regularidad para arrancar hasta $C^\infty$. En 3D, la misma identidad produce una cota más débil que la crítica. No se conoce ninguna estimación adicional a priori que cierre la brecha.

Tao subrayó esta barrera. No se espera que argumentos basados solo en la desigualdad de energía y en el escalamiento resuelvan la regularidad global 3D, así que cualquier prueba exitosa probablemente tendrá que explotar estructura adicional, algo que el análisis de escala por sí solo no ve. Los métodos "ciegos a la supercriticidad" (que tratan la ecuación solo a través de su escalamiento y su estructura energética) no pueden funcionar. La estructura algebraica específica de la ecuación, en particular la condición libre de divergencia y la estructura antisimétrica de la no linealidad, tendría que entrar en juego. Consulte Por qué Navier-Stokes es difícil para un tratamiento más profundo.

La prueba 2D funciona porque la vorticidad permanece limitada y el escalado es fundamental. El 3D no tiene ninguna de las dos cosas. Entonces, ¿qué necesitaría una prueba?

Nadie lo sabe. Pero esto es lo que los investigadores están persiguiendo:

• Encontrar una nueva \"perilla de control\". La vorticidad es la perilla de control del 2D: permanece limitada y todo lo demás se deriva de ese único hecho. En 3D, necesitamos una cantidad diferente, algo que se mantenga dócil independientemente de lo que haga el fluido y que sea lo suficientemente potente como para obligar a toda la solución a permanecer suave para siempre. Nadie lo ha encontrado. Los investigadores han estado buscando durante décadas y todavía no se encuentra.
• Explota la estructura oculta. Los fluidos son incompresibles. No se pueden apretar. Esa restricción limita lo que puede hacer el estiramiento del vórtice, y puede haber patrones geométricos más profundos enterrados en las ecuaciones que nadie ha explotado completamente todavía.
• Demuestre que realmente se rompe. Tal vez las soluciones 3D puedan explotar. Eso sería igualmente enorme. Necesitaría construir una condición inicial específica en la que el estiramiento del vórtice supere la viscosidad y lleve la solución al infinito en un tiempo finito, y para las ecuaciones de Euler más simples (Navier-Stokes sin fricción), la formación de singularidades se ha demostrado en entornos relacionados, pero el caso viscoso permanece completamente abierto.

Para obtener más información sobre lo que se ha probado, consulte Subproblemas de Navier-Stokes.

Una prueba de regularidad global 3D requeriría cerrar la brecha de supercriticidad. En concreto, se necesita una estimación a priori de la forma $\|u(t)\|_X \leq C(\|u_0\|_Y, t)$ para alguna norma $X$ en la escala crítica o por encima de ella, donde $C$ permanezca finita para todo $t$. La cota de energía conocida $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$ está media derivada por debajo de la crítica y no basta.

Varios programas de investigación apuntan a esta brecha:

• Descomposición de perfiles y concentración-compacidad. Adaptados del éxito en ecuaciones dispersivas críticas (Kenig-Merle 2006), estos métodos buscan clasificar perfiles de explosión. Para Navier-Stokes hay resultados parciales (por ejemplo, Gallagher-Koch-Planchon 2016), pero el carácter supercrítico de la energía vuelve más difícil ejecutar el programa completo que en los contextos de ondas o Schrödinger de energía crítica.
• Extensiones de soluciones leves. El marco de Fujita-Kato (1964) da buen planteamiento local en $\dot{H}^{1/2}(\mathbb{R}^3)$ y buen planteamiento global para datos pequeños en espacios críticos ($L^3$, $\dot{H}^{1/2}$, $BMO^{-1}$). La pregunta es si las soluciones de datos grandes pueden continuarse globalmente, lo que exige controlar la norma crítica.
• Criterios de regularidad. Más allá de Beale-Kato-Majda ($\int_0^T \|\omega\|_{\infty} < \infty$), están las condiciones de Prodi-Serrin ($u \in L^p_t L^q_x$ con $2/p + 3/q = 1$, $q > 3$), Escauriaza-Seregin-Šverák ($u \in L^\infty_t L^3_x$, 2003) y otros criterios de punto extremo. Cada uno reduce la regularidad global a una sola estimación a priori, pero demostrar esa estimación sigue abierto.
• Construir explosión. Tao (2016) construyó una solución con explosión para un sistema promediado de Navier-Stokes que respeta la identidad de energía y el escalamiento, pero no la estructura completa libre de divergencia. Esto nos dice que cualquier prueba de regularidad debe usar la estructura geométrica específica de la no linealidad, no solo sus propiedades de escala. Sigue abierto si el verdadero Navier-Stokes admite explosión.

Para el problema no viscoso, la explosión $C^{1,\alpha}$ de Elgindi para Euler 3D (2021) muestra que el estiramiento de vórtice puede producir singularidades por debajo de la regularidad $C^\infty$. La cuestión de la explosión suave de Euler ($C^\infty$) sigue abierta, al igual que la de si la viscosidad puede frenar tales mecanismos en el marco de Navier-Stokes.

Todo lo anterior, en una tabla:

Esto no es un tecnicismo. La brecha entre 2D y 3D es un abismo. La estrategia de prueba que funciona en dos dimensiones no solo "necesita un poco más de trabajo" para pasar a tres; fundamentalmente no puede funcionar porque la estructura matemática de la que depende, el principio máximo para la vorticidad y la criticidad energética que hacen tan manejable al 2D, simplemente no existen en 3D.

Para ver las ecuaciones completas, consulte ¿Qué son las ecuaciones de Navier-Stokes?. Para el problema abierto preciso, consulte Existencia y suavidad de Navier-Stokes. Para saber por qué es tan difícil, consulte Por qué Navier-Stokes es difícil.

Los siguientes contrastes resumen la división matemática:

El resultado 2D no es simplemente un calentamiento de dimensiones inferiores. Es un teorema completo cuyo mecanismo de prueba (el principio de máxima vorticidad combinado con la criticidad energética) no tiene una contraparte tridimensional conocida. Cualquier resolución del problema 3D, ya sea de regularidad o de blow-up, requerirá ideas fundamentalmente nuevas. Para conocer el estado actual de los resultados parciales y los programas de investigación, consulte Subproblemas de Navier-Stokes y Por qué Navier-Stokes es difícil.