Enfoques del problema de Navier-Stokes

El enfoque más antiguo comienza con la energía. Un fluido en movimiento transporta energía cinética y la viscosidad la devora, como la fricción que detiene las cosas. La energía total solo puede disminuir con el tiempo, asumiendo que nada bombea energía desde el exterior.

Leray vio esto en 1934 e hizo un movimiento clave: usar la energía destinada a demostrar que tiene que existir una solución global débil con energía cinética finita. Construir soluciones aproximadas, suavizadas artificialmente. Demuestre que todos obedecen el límite de energía. Toma un límite. Algo debe sobrevivir en ese límite, y lo hace.

Pero aquí está el truco. Los límites de energía son instrumentos contundentes. Garantizan que el fluido tiene una energía total finita, claro, pero no pueden decirte que la velocidad permanece finita en cada punto del espacio y el tiempo. Esa brecha entre "energía finita" y "suave en todas partes" es exactamente el problema de regularidad, y ha estado abierta durante noventa años.

Enlaces de artículos: Leray, Sur le mouvement d'un liquide visqueux emplissant El espacio (1934); Hopf, Über die Anfangswertaufgabe für die hydrodynamischen Grundgleichungen (1951).

La construcción de Leray-Hopf se desarrolla en dos variantes principales: la aproximación de Galerkin (proyectando sobre subespacios de dimensión finita) o la regularización de Friedrichs (suavizando la no linealidad mediante convolución). Ambos comparten el mismo esqueleto de cuatro pasos. La estrategia estándar tiene cuatro pasos:

1. Aproximado: Resolver el sistema suavizado $\partial_t u_\varepsilon + (J_\varepsilon u_\varepsilon \cdot \nabla) u_\varepsilon = \nu \Delta u_\varepsilon - \nabla p_\varepsilon$ en subespacios de dimensión finita.
2. Cota de energía: La estimación a priori $\|u_\varepsilon(t)\|_{L^2}^2 + 2\nu \int_0^t \|\nabla u_\varepsilon\|_{L^2}^2 \leq \|u_0\|_{L^2}^2$ se mantiene uniformemente en $\varepsilon$.
3. Compacidad: Extrae una subsecuencia débilmente convergente $u_\varepsilon \rightharpoonup u$ en $L^2_t \dot{H}^1_x$ usando el lema de Aubin-Lions.
4. Pasa al límite: El término no lineal converge por una fuerte convergencia $L^2_{\text{loc}}$ de $u_\varepsilon$.

La solución débil resultante satisface la desigualdad de energía, no la igualdad. La energía se puede perder en momentos irregulares. Y ese es el problema: la brecha entre la clase energética $L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$ y la suavidad real es precisamente lo que no podemos cerrar.

Enlaces de artículos: Leray, Sur le mouvement d'un liquide visqueux emplissant l'espace (1934); Hopf, Über die Anfangswertaufgabe für die hydrodynamischen Grundgleichungen (1951).

El enfoque de Caffarelli-Kohn-Nirenberg (1982) no intenta demostrar una suavidad total. Pregunta algo completamente distinto: ¿qué tan malas pueden ser realmente las singularidades?

Apenas malas en absoluto. Su teorema de regularidad $\varepsilon$ dice que si ciertas cantidades locales invariantes de escala son lo suficientemente pequeñas en una región espacio-temporal pequeña, la solución allí es automáticamente suave. Y dado que la energía total es finita, simplemente no hay suficiente "presupuesto" para que coexistan muchos puntos singulares.

Piénselo de esta manera. Una pared puede tener grietas. Pero la longitud total de todas esas grietas combinadas es cero, lo que significa que el conjunto singular es extremadamente pequeño en el sentido teórico de la medida parabólica (medida parabólica unidimensional de Hausdorff cero).

Enlaces de artículos: Caffarelli-Kohn-Nirenberg, Regularidad parcial de soluciones débiles adecuadas de las ecuaciones de Navier-Stokes (1982); Albritton-Barker-Prange, regularidad Epsilon para las ecuaciones de Navier-Stokes mediante unicidad débil-fuerte.

El teorema de regularidad $\varepsilon$ de CKN: para soluciones débiles adecuadas, la pequeñez de ciertas cantidades locales invariantes de escala (que involucran tanto $u$ como $p$ en un cilindro parabólico $Q_r(z_0) = B_r(x_0) \times (t_0 - r^2, t_0)$) implica regularidad (continuidad de Hölder) cerca del punto $z_0 = (x_0, t_0)$.

La prueba combina la desigualdad de energía local con una iteración de tipo Campanato: si la energía invariante de escala es pequeña, un argumento de arranque muestra que $u$ está acotado, luego es Hölder y luego suave por la teoría clásica de Schauder.

La estimación dimensional $\mathcal{P}^1(\Sigma) = 0$ se obtiene mediante un argumento de recubrimiento de Vitali. Los puntos singulares fuerzan cotas inferiores sobre cantidades localizadas de disipación/energía invariantes de escala; el recubrimiento controla entonces la medida parabólica unidimensional del conjunto singular mediante una medida local de energía, que se anula a escalas pequeñas.

Enlaces de artículos: Caffarelli-Kohn-Nirenberg, Regularidad parcial de soluciones débiles adecuadas de las ecuaciones de Navier-Stokes (1982); Albritton-Barker-Prange, regularidad Epsilon para las ecuaciones de Navier-Stokes mediante unicidad débil-fuerte.

Aquí hay una reducción drástica de todo el problema. Beale, Kato y Majda demostraron en 1984 que, para las ecuaciones 3D de Euler, la explosión sólo puede ocurrir si se pierde el control de la vorticidad. Posteriormente se establecieron criterios análogos para Navier-Stokes. Eso es todo. Una condición.

La vorticidad mide el giro local. El criterio BKM dice: mantenga el giro máximo limitado en la norma correcta y la solución se mantendrá suave. Todo lo demás se alinea automáticamente.

Una familia de cantidades para controlar. Desafortunadamente, controlarlos ha resultado ser exactamente tan difícil como el problema original. La reducción es limpia. La ejecución sigue fuera de alcance.

Enlaces del artículo: Beale-Kato-Majda, Comentarios sobre la ruptura de soluciones suaves para las ecuaciones tridimensionales de Euler (1984); Kozono-Taniuchi, Estimaciones bilineales en BMO y las ecuaciones de Navier-Stokes (2000).

El teorema original de Beale-Kato-Majda (1984) es para las ecuaciones tridimensionales de Euler. Para Navier-Stokes, criterios de continuación análogos implican que una solución suave $u$ en $[0, T^*)$ se extiende más allá de $T^*$ siempre que

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

donde $\omega = \nabla \times u$ es la vorticidad. Los refinamientos incluyen:

• Kozono-Taniuchi (2000): $\|\omega\|_{L^\infty}$ puede ser reemplazado por $\|\omega\|_{\mathrm{BMO}}$
• Variantes del espacio Besov: el control Besov crítico o límite también puede servir como criterio de continuación
• Criterios de dirección restringida: Las condiciones de tipo Serrin en los componentes de $\nabla u$ también pueden servir como criterios de continuación (ver, por ejemplo, Beirão da Veiga, 1995)

Estos criterios se conectan con la imagen de estiramiento del vórtice: cualquier singularidad de tiempo finito debe forzar que la vorticidad se acumule demasiado rápido para que la integral de tiempo anterior permanezca finita.

Enlaces de artículos: Beale-Kato-Majda, Comentarios sobre la ruptura de soluciones suaves para las ecuaciones tridimensionales de Euler (1984); Kozono-Taniuchi, Estimaciones bilineales en BMO y las ecuaciones de Navier-Stokes (2000); Chemin-Planchon, Límites automejorantes para las ecuaciones de Navier-Stokes (2012).

Un ángulo más moderno trabaja con espacios funcionales (como $L^3$ o $\dot{H}^{1/2}$) que se sitúan justo en el borde de lo que permite la simetría de escala. Esos son los espacios críticos, y ahí viven los resultados de regularidad más afilados.

La lógica es limpia: si puedes demostrar que una solución se mantiene dentro de ciertas cotas críticas, la suavidad sigue automáticamente. Varios equipos lo han demostrado, construyendo todo un menú de criterios de regularidad.

El problema es la brecha. Los métodos energéticos solo dan control energético supercrítico, por debajo de la escala crítica. Lo que necesitamos son cotas críticas. Esa brecha es estrecha, a veces de una sola derivada de regularidad, pero ha resistido todos los intentos de cerrarla.

Enlaces de artículos: Koch-Tataru, Buen planteamiento para las ecuaciones de Navier-Stokes (2001); Kenig-Koch, un enfoque alternativo a la regularidad de las ecuaciones de Navier-Stokes en espacios críticos; Gallagher-Koch-Planchon, un enfoque de descomposición de perfiles para el criterio de regularidad $L^\infty_t(L^3_x)$ Navier-Stokes. Para una comparación detallada de por qué la criticidad energética funciona en 2D pero falla en 3D, consulte Por qué 2D es más fácil que 3D.

Programas principales en regularidad de espacios críticos:

• Koch-Tataru (2001): Buen planteamiento global para datos pequeños en $\text{BMO}^{-1}$. La estimación clave es una cota de punto fijo para la aplicación bilineal de Duhamel en el espacio de soluciones de flujo de calor de Koch-Tataru construido sobre datos en $\text{BMO}^{-1}$, no una mera estimación estática de producto en $\text{BMO}^{-1}$. Este es uno de los resultados más afilados conocidos para métodos perturbativos.
• Gallagher-Koch-Planchon (2013): Enfoque de descomposición de perfiles para el criterio de regularidad $L^\infty_t L^3_x$ de Navier-Stokes. Toda sucesión de soluciones con norma crítica acotada tiene una subsucesión que se descompone en perfiles asintóticamente desacoplados.

La obstrucción central es que no se conoce ningún funcional coercitivo que esté a la vez controlado por la evolución y sea crítico con respecto al escalamiento de Navier-Stokes.

Enlaces de artículos: Koch-Tataru, Buen planteamiento para las ecuaciones de Navier-Stokes (2001); Kenig-Koch, un enfoque alternativo a la regularidad de las ecuaciones de Navier-Stokes en espacios críticos; Gallagher-Koch-Planchon, Un enfoque de descomposición de perfiles para el criterio de regularidad $L^\infty_t(L^3_x)$ Navier-Stokes.

La teoría moderna de la EDP se basa en gran medida en el análisis armónico. La idea central: dividir una función en ondas de diferentes frecuencias, de la misma manera que dividirías un acorde musical en notas individuales. Excepto que aquí las "notas" son oscilaciones espaciales de la velocidad del fluido a escalas tremendamente diferentes.

La descomposición de Littlewood-Paley hace exactamente esto. Divida el campo de velocidad en componentes escala por escala. Realice un seguimiento de cómo fluye la energía entre ellos. De repente, la intuición física informal de la "cascada de energía" se convierte en algo sobre lo que realmente se pueden demostrar teoremas, y los teoremas son precisos. Estos métodos han producido muchos de los mejores resultados en criterios de regularidad y tasas de explosión.

Enlaces de artículos: Cannone-Meyer, descomposición de Littlewood-Paley y ecuaciones de Navier-Stokes (1995); Gallagher-Koch-Planchon, Un enfoque de descomposición de perfiles para el criterio de regularidad $L^\infty_t(L^3_x)$ Navier-Stokes.

La teoría de Littlewood-Paley descompone $u = \sum_j \Delta_j u$ donde $\Delta_j$ localiza las frecuencias $|\xi| \sim 2^j$. Aplicado a la no linealidad del transporte $(u \cdot \nabla)u$:

La descomposición en paraproducto de la no linealidad $(u \cdot \nabla)u$ se divide en interacciones de frecuencia baja-alta, alta-baja y alta-alta:

$$(u \cdot \nabla)u = T_u \nabla u + T_{\nabla u} u + R(u, \nabla u)$$

donde $T$ es el paraproducto y $R$ el resto. Cada pieza tiene diferentes propiedades de regularidad en los espacios de Besov $\dot{B}^s_{p,q}$.

Resultados clave al utilizar esta maquinaria:

• Espacios de Chemin-Lerner: $\widetilde{L}^\rho_T \dot{B}^s_{p,q}$ proporcionan el marco natural para una buen planteamiento crítica: los mapas de forma bilineal de Navier-Stokes $\widetilde{L}^\infty_T \dot{B}^{-1+3/p}_{p,q} \times \widetilde{L}^1_T \dot{B}^{1+3/p}_{p,q} \to \widetilde{L}^1_T \dot{B}^{-1+3/p}_{p,q}$.
• Cannone-Meyer: Los métodos de Littlewood-Paley brindan una formulación limpia de ondas/Besov de la teoría de los datos pequeños.

Enlaces de artículos: Cannone-Meyer, descomposición de Littlewood-Paley y ecuaciones de Navier-Stokes (1995); Gallagher-Koch-Planchon, Un enfoque de descomposición de perfiles para el criterio de regularidad $L^\infty_t(L^3_x)$ Navier-Stokes.

Aquí hay un instinto completamente diferente. En lugar de rastrear números (normas, energías), estos métodos estudian la forma de la solución: cómo se doblan los tubos de vórtice, cómo se organizan en el espacio las regiones de intensa rotación.

La idea clave es que la explosión no se trata sólo de que algo crezca. Se trata de que el fluido se organice en una configuración geométrica muy específica. Si puede demostrar que la configuración es imposible (porque contradice la estructura de disipación de energía, la incompresibilidad o ambas), habrá descartado la explosión sin siquiera calcular una norma.

Este punto de vista geométrico se ha convertido en un punto de vista que ha inspirado varios criterios de regularidad rigurosos junto con métodos puramente analíticos. Y se siente diferente. Pregunta ¿qué forma adopta el desastre? en lugar de ¿qué tan grande puede llegar a ser este número?

Enlaces de artículos: Constantin-Fefferman, Restricciones geométricas sobre soluciones potencialmente singulares para Euler tridimensional ecuaciones (1993); Albritton-Barker-Prange, Suavizado localizado y concentración para las ecuaciones de Navier-Stokes en el medio espacio.

Los enfoques geométrico-topológicos explotan restricciones estructurales que son invisibles para los métodos puramente analíticos:

• Geometría de línea de vórtice: Constantin y Fefferman (1993) demostraron que si el campo de dirección de vorticidad $\hat{\omega} = \omega/|\omega|$ es Lipschitz en regiones de alta vorticidad, la solución es regular. La explosión requiere que la dirección de la vorticidad desarrolle una singularidad simultáneamente con la magnitud.
• Argumentos de incompatibilidad: si una configuración de explosión está geométricamente restringida (por ejemplo, a través de límites de empaquetamiento en el número de regiones de concentración independientes que se ajustan a los presupuestos de energía y disipación), se puede derivar una contradicción sin estimar directamente las normas críticas.
• Partición de casos (especulativa/programática): una estrategia propuesta clasifique cada región espacial como perteneciente a uno de un número finito de escenarios (por ejemplo, localmente regular, tipo I, tipo II, densamente empaquetado) e intente mostrar que cada escenario da regularidad o transfiere el problema a un argumento de conteo acotado. Esto sigue siendo un programa de investigación más que un resultado establecido.

La página de prueba de este sitio explora argumentos de este tipo; nada aquí se presenta como una prueba formal completa del problema del Milenio.

Enlaces de artículos: Constantin-Fefferman, Restricciones geométricas en soluciones potencialmente singulares para las ecuaciones tridimensionales de Euler (1993); Albritton-Barker-Prange, Suavizado localizado y concentración para las ecuaciones de Navier-Stokes en el medio espacio.

Esto tomó a mucha gente por sorpresa. Las soluciones débiles del método de Leray (sección 1) resultan ser no únicas, al menos cuando hay forzamiento externo.

El arma es la integración convexa, una técnica nacida en problemas de geometría y adaptada a ecuaciones de fluidos por De Lellis y Székelyhidi a partir de 2009. La idea es construir soluciones "salvajes" apilando iterativamente correcciones de alta frecuencia que, en conjunto, satisfacen la ecuación pero se comportan de forma errática.

Para 3D Navier-Stokes, Buckmaster y Vicol (2019) demostraron no unicidad de soluciones débiles de energía finita por debajo de la clase de Leray-Hopf. La no unicidad para Euler viene del programa de integración convexa de De Lellis-Székelyhidi, Isett y Buckmaster. Luego, en 2022, Albritton, Brué y Colombo demostraron que incluso las soluciones de Leray-Hopf de 3D Navier-Stokes no son únicas cuando hay fuerza externa. Sigue abierto si la no unicidad persiste para las ecuaciones no forzadas de Navier-Stokes.

¿Por qué importa esto? Porque "existe una solución débil" ha sido el gran titular desde 1934. Ahora sabemos que eso no fija una respuesta única. La pregunta se afila: ¿qué solución, si alguna, es la físicamente correcta?

Enlaces de artículos: De Lellis-Székelyhidi, flujos continuos disipativos de Euler (2013); Buckmaster-Vicol, No unicidad de las soluciones débiles de la ecuación de Navier-Stokes (2019); Albritton-Brué-Colombo, No unicidad de las soluciones de Leray de las ecuaciones forzadas de Navier-Stokes (2022).

La integración convexa para ecuaciones de fluidos se origina en el programa De Lellis-Székelyhidi (2009-2013), que adapta la técnica de incrustación isométrica $C^1$ de Nash-Kuiper para construir soluciones débiles de Euler que disipan energía. Etapas clave:

• De Lellis-Székelyhidi (2013): existencia de flujos de Euler disipativos continuos ($C^0$) en $\mathbb{T}^3$ (regularidad de Hölder por debajo de $1/5$ lograda en el trabajo posterior de Buckmaster-De Lellis-Isett-Székelyhidi de 2015; luego mejorada a $< 1/3$ por Isett, 2018, resolviendo el lado flexible de la conjetura de Onsager).
• Buckmaster-Vicol (2019): no unicidad de soluciones débiles de energía finita para Navier-Stokes 3D, con soluciones en $C^0_t H^\beta_x$ para algún $\beta > 0$. La construcción usa flujos de Beltrami intermitentes como bloques de construcción, agregando correcciones oscilatorias en cada paso de la iteración mientras mantiene el control del esfuerzo de Reynolds. Esto queda por debajo de la clase energética de Leray-Hopf, así que no contradice directamente la unicidad de Leray.
• Albritton-Brué-Colombo (2022): no unicidad de soluciones de Leray-Hopf para las ecuaciones 3D forzadas de Navier-Stokes. La prueba construye una solución autosimilar inestable de fondo y usa un mecanismo de inestabilidad para ramificarse en soluciones distintas de Leray-Hopf a partir de los mismos datos iniciales. Esto muestra que la desigualdad de energía por sí sola no selecciona una solución única cuando hay forzamiento.

La pregunta central abierta es si la no unicidad persiste para las ecuaciones de Navier-Stokes no forzadas en la clase de Leray-Hopf. El resultado con forzamiento muestra que la desigualdad energética no es un principio de selección suficiente, pero no resuelve si la ecuación no forzada tiene una estructura adicional que restaure la unicidad.

Enlaces del artículo: De Lellis-Székelyhidi, Flujos de Euler continuos disipativos (2013); Buckmaster-Vicol, No unicidad de las soluciones débiles de la ecuación de Navier-Stokes (2019); Albritton-Brué-Colombo, No unicidad de las soluciones de Leray de las ecuaciones forzadas de Navier-Stokes (2022).

¿Podemos al menos descartar ciertas estrategias de prueba? Terence Tao demostró en 2016 que sí, podemos. Y el resultado es aleccionador.

Tao construyó una versión modificada de las ecuaciones de Navier-Stokes, un sistema "promediado", que mantiene muchas características estructurales clave de las ecuaciones reales: la identidad energética, la forma en que crece la enstrofia (una medida de la intensidad de la vorticidad), la simetría de escala. Pero en este sistema modificado, las soluciones explotan en un tiempo finito.

La implicación descarta familias amplias de estrategias de prueba. Cualquier prueba de que la suavidad global sea válida para las ecuaciones reales debe utilizar alguna propiedad estructural específica de la no linealidad verdadera que el sistema promediado no tiene. No se puede demostrar la regularidad utilizando únicamente límites de energía, escalamiento y crecimiento de enstrofia. Esas herramientas por sí solas son consistentes con la explosión.

Esto no significa que las ecuaciones reales exploten. Dice que familias enteras de estrategias de prueba son callejones sin salida. La prueba final (si se mantiene la regularidad) debe ser más contundente que un argumento energético genérico. Mucho más nítido.

Enlaces de artículos: Tao, blow-up en tiempo finito para una ecuación tridimensional promediada de Navier-Stokes (2016).

Tao (2016) considera un sistema de la forma

$$\partial_t u + \tilde{B}(u, u) = \nu \Delta u - \nabla p, \quad \nabla \cdot u = 0,$$

donde $\tilde{B}$ es un operador bilineal que concuerda con la verdadera no linealidad de Navier-Stokes $(u \cdot \nabla)u$ en los siguientes sentidos:

• es un multiplicador de Fourier de orden 1, preservando el escalamiento $u(x,t) \mapsto \lambda u(\lambda x, \lambda^2 t)$
• satisface la misma identidad energética: $\langle \tilde{B}(u,u), u \rangle = 0$
• se reproduce la estructura de crecimiento de la enstrofia

Para este sistema promediado, Tao construye datos iniciales suaves cuya solución explota en un tiempo finito. El mecanismo de explosión programa una secuencia de concentraciones de enstrofia cada vez más agudas a escalas cada vez más pequeñas, con cada etapa duplicando la enstrofia en una cascada controlada.

La obstrucción que esto crea: cualquier cantidad supercrítica que esté (a) controlada por la evolución en la clase energética y (b) invariante bajo la escala de Navier-Stokes no puede por sí sola descartar la explosión, porque también estaría controlada en el sistema promediado, que sí explota. Una prueba de regularidad debe explotar la estructura de cancelación algebraica específica del término de transporte verdadero $(u \cdot \nabla)u$ que el operador promediado $\tilde{B}$ no comparte.

Enlaces en papel: Tao, explosión de tiempo finito para una ecuación tridimensional promediada de Navier-Stokes (2016).

Este artículo es parte de Progreso.

Desde la prueba de existencia de Leray de 1934 a través de la integración convexa y las barreras de prueba de Tao, estas son las principales estrategias que la gente ha lanzado al problema 3D de Navier-Stokes. Ninguno ha resuelto por completo el problema de la regularidad 3D. Para conocer el contexto sobre cómo la viscosidad da forma a las matemáticas en comparación con las ecuaciones invisibles de Euler, consulte Euler vs. Navier-Stokes. Para conocer el estado actual, consulte ¿Se ha resuelto el problema de Navier-Stokes? Para conocer la declaración formal exacta, regrese a El problema del Milenio.

Este artículo es parte de Progreso.

Los enfoques anteriores representan los principales hilos rigurosos en la literatura sobre regularidad y unicidad hasta 2022. Ninguna combinación ha resuelto el problema 3D completo. El campo sigue moviéndose. Los métodos de prueba asistidos por computadora son un área activa que se trata por separado en este sitio.

Para comparar los sistemas viscosos y no viscosos (y por qué la viscosidad ayuda pero no es suficiente), consulte Euler vs. Navier-Stokes. Para conocer los subproblemas a los que se dirigen estos enfoques, consulte Subproblemas. Para conocer los obstáculos para escalar, consulte Por qué es difícil.