Por qué el problema de Navier-Stokes es difícil

Muchas de las ecuaciones que la gente encuentra por primera vez en física son lineales: duplica la entrada y la respuesta se duplica. Navier-Stokes no es así.

¿Navier-Stokes? No lineal. La velocidad del fluido afecta su propia tasa de cambio, lo que significa que el fluido se empuja a sí mismo. Imagínese tratar de predecir adónde irá una multitud cuando el movimiento de cada persona depende de lo que hacen todos los que la rodean, y lo que hacen esas personas depende de todos los que las rodean, en una espiral hacia afuera para siempre. Esa es la situación que estás viendo.

El culpable es el término de autointeracción $(u \cdot \nabla)u$. Crea bucles de retroalimentación donde las pequeñas perturbaciones se amplifican hasta convertirse en grandes, y es por eso que la turbulencia de fluidos es tan tremendamente compleja (consulte subproblemas para obtener más información).

La no linealidad convectiva $(u \cdot \nabla)u$ es el obstáculo fundamental. En la formulación de vorticidad $\omega = \nabla \times u$, la ecuación se convierte en

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

El término de estiramiento del vórtice $(\omega \cdot \nabla)u$ no tiene signo. Puede amplificar la vorticidad sin límites. En 2D este término desaparece (ya que $\omega$ es un escalar perpendicular al flujo), por lo que se conoce la regularidad global 2D (Ladyzhenskaya, 1969). En 3D, el estiramiento del vórtice es el principal mecanismo candidato para la explosión en tiempo finito.

Aquí está el quid de la cuestión: la no linealidad es cuadrática en $u$. La estimación de energía $H^1$ da $\|\nabla u\|_{L^2}$, pero controlar $(u \cdot \nabla)u$ en $L^2$ normalmente necesita información más sólida como $u \in L^\infty$ junto con $\nabla u \in L^2$, o un control de escala crítica equivalente; la clase energética por sí sola no proporciona esto.

Las ecuaciones de Navier-Stokes tienen una simetría de escala. Amplíe una solución, haga todo más pequeño y más rápido en las cantidades correctas y obtendrá otra solución perfectamente válida. Esa simetría es matemáticamente natural, pero analíticamente peligrosa.

¿Por qué? La única cantidad que podemos controlar de manera confiable es la energía total del fluido, y se encuentra en una escala completamente incorrecta, lo que nos informa sobre el panorama general pero no dice absolutamente nada sobre lo que está sucediendo a escalas microscópicas donde realmente se formaría una explosión.

Piense en monitorear el uso total de electricidad de una ciudad para detectar una chispa. ¿Útil? Seguro. ¿Lo suficientemente fino? Ni siquiera cerca. Esa brecha es todo el problema.

Bajo la escala natural $u_\lambda(x,t) = \lambda u(\lambda x, \lambda^2 t)$, el espacio de Sobolev crítico es $\dot{H}^{1/2}(\mathbb{R}^3)$ (equivalente a $L^3$). Una cantidad es:

• Subcrítica si la norma se reduce con el cambio de escala, lo que significa que captura el comportamiento a gran escala pero omite la concentración a pequeña escala, por ejemplo, $\|u\|_{L^2}$
• Crítica si es invariante de escala, por ejemplo, $\|u\|_{L^3}$, $\|u\|_{\dot{H}^{1/2}}$
• Supercrítico si la norma crece con el reescalamiento, lo que significa que la concentración a pequeña escala se vuelve más difícil de controlar

La desigualdad energética da control de $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$. Ambos componentes son subcríticos:

$$\|u_\lambda\|_{L^2} = \lambda^{-1/2} \|u\|_{L^2}, \quad \|u_\lambda\|_{L^2_t \dot{H}^1_x} = \lambda^{-1/2} \|u\|_{L^2_t \dot{H}^1_x}$$

Entonces la estimación de energía da cero control a pequeña escala. La no linealidad puede, en principio, superar la disipación a escalas finas. Pasar de la clase de energía subcrítica a una norma crítica es la dificultad central.

Observa un río. El movimiento fluido se vuelve caótico. Turbulento. Los grandes remolinos se fragmentan en otros más pequeños, que a su vez se fragmentan en otros aún más pequeños, cayendo en cascada hasta escalas microscópicas donde la viscosidad finalmente suaviza las cosas.

Kolmogorov describió esta cascada de energía en 1941, y las ecuaciones de Navier-Stokes la capturan maravillosamente. Pero esto es lo que mantiene a la gente despierta por la noche: ¿qué pasa si la energía se concentra en regiones cada vez más pequeñas más rápido de lo que la viscosidad puede disiparla? Eso es una explosión.

¿Puede suceder realmente? ¿O siempre gana la viscosidad? Esa es la pregunta abierta, punto. Para conocer el puente físico entre el número de Reynolds y esta imagen a pequeña escala, consulte Número de Reynolds, turbulencia y por qué son importantes las pequeñas escalas.

La teoría K41 de Kolmogorov predice un espectro de energía $E(k) \sim \varepsilon^{2/3} k^{-5/3}$ en el rango inercial $k_f \ll k \ll k_\eta$, donde $k_\eta \sim (\varepsilon/\nu^3)^{1/4}$ es el número de onda de disipación de Kolmogorov. Flujo de energía constante a través de escalas. Ésa es la imagen clara.

La pregunta sobre la regularidad plantea algo más oscuro: ¿puede esta cascada degenerar? ¿Puede la escala de disipación $k_\eta^{-1}$ reducirse a cero en un tiempo finito, con $\|\nabla u\|_{L^2} \to \infty$ mientras que la energía total permanece finita?

La conjetura de la anomalía de disipación de Onsager (1949) dice que sí, en cierto sentido: en el límite de viscosidad evanescente $\nu \to 0$, la disipación de energía persiste y las soluciones débiles de Euler pueden disipar energía. Confirmado para exponentes de Hölder inferiores a $1/3$ (Isett, 2018; Buckmaster et al., 2018). Pero lo que esto significa para la regularidad de Navier-Stokes sigue siendo realmente poco claro. Para conocer una intuición a nivel de régimen, consulte El número de Reynolds, la turbulencia y por qué son importantes las pequeñas escalas.

La presión en las ecuaciones de Navier-Stokes es extraña. No es una variable independiente en absoluto; la velocidad la determina completamente mediante una única restricción: el fluido es incompresible, por lo que no se puede exprimir.

Esto hace que la presión sea no local. En el modelo incompresible, cambiar la velocidad en una región afecta el campo de presión globalmente, porque la presión está determinada por todo el campo de velocidad en ese momento.

Para el análisis, eso es devastador. No se puede estudiar lo que sucede en un solo punto sin tener en cuenta todo el fluido a la vez. ¿Razonamiento local? Olvídalo.

La restricción de incompresibilidad $\nabla \cdot u = 0$ determina la presión a través de la ecuación de Poisson

$$-\Delta p = \nabla \cdot ((u \cdot \nabla)u) = \partial_i \partial_j (u_i u_j)$$

so $p = (-\Delta)^{-1} \partial_i \partial_j (u_i u_j)$, que involucra transformadas de Riesz (operadores integrales singulares). La presión es una función no local de la velocidad, y esta no localidad es el obstáculo principal para las estimaciones puntuales o locales en el espacio.

Los argumentos del principio máximo estándar fallan aquí: aunque el término viscoso $\nu \Delta u$ es disipativo, el gradiente de presión $-\nabla p$ puede concentrar energía de regiones distantes. La teoría de Caffarelli-Kohn-Nirenberg maneja esto a través de desigualdades de energía locales en cilindros parabólicos, pero extraer regularidad puntual de estos sigue siendo el paso difícil.

Para las ecuaciones 2D incompresibles de Navier-Stokes en la configuración estándar, se conocen soluciones globales suaves; esto quedó establecido en obras clásicas, incluida la de Ladyzhenskaya (1969). Vea por qué 2D es más fácil.

¿Tres dimensiones? Todo se desmorona y la razón se reduce a un mecanismo: el estiramiento del vórtice. En 2D, los vórtices pueden girar y fusionarse pero no pueden estirarse. En 3D, el fluido puede agarrar tubos de vórtice y hacerlos cada vez más delgados, concentrando potencialmente hasta el último bit de energía en un filamento infinitamente delgado.

¿Puede este estiramiento llegar al infinito en un tiempo finito, o la viscosidad siempre interviene? Ésa es la pregunta del millón. Literalmente.

La dicotomía es marcada.

2D: La vorticidad $\omega$ es un escalar que satisface $\partial_t \omega + u \cdot \nabla \omega = \nu \Delta \omega$. El principio máximo da $\|\omega(t)\|_{L^\infty} \leq \|\omega_0\|_{L^\infty}$, BKM implica regularidad global, y el término de estiramiento del vórtice $(\omega \cdot \nabla)u$ es idénticamente cero.

3D: Vorticidad $\omega \in \mathbb{R}^3$ satisface $\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega$, donde el término de extensión $(\omega \cdot \nabla)u$ puede amplificar $|\omega|$ superlinealmente (formalmente $\sim |\omega|^2$ a través de Biot-Savart). No hay principio máximo disponible.

La enstrofia $\|\omega\|_{L^2}^2$ satisface

$$\frac{d}{dt}\|\omega\|_{L^2}^2 \leq C\|\omega\|_{L^2}^2 \|\nabla u\|_{L^\infty} - 2\nu \|\nabla \omega\|_{L^2}^2$$

Controlar $\|\nabla u\|_{L^\infty}$ requiere $\omega \in L^\infty$, lo que requiere controlar $\|\nabla u\|_{L^\infty}$. Circular. Ninguna técnica existente ha roto este bucle.

Este artículo es parte de El problema.

Estos obstáculos han llevado a los matemáticos a descomponer el problema en subproblemas y desarrollar enfoques especializados para cada uno.

Para conocer el contexto de por qué el término viscoso ayuda pero no es suficiente, consulte Euler vs. Navier-Stokes.

Este artículo es parte de El problema.

El problema de la regularidad se divide en partes manejables; consulte Subproblemas. Las herramientas analíticas creadas para atacar estos obstáculos se tratan en Enfoques. ¿Por qué el término viscoso ayuda pero aún no cierra la brecha? Euler frente a Navier-Stokes.