Ecuaciones de Navier–Stokes

1. Por qué importa (problema del milenio)

El agua que fluye alrededor de un barco y el aire que rodea un ala obedecen el mismo balance local: la masa se conserva y las fuerzas cambian la cantidad de movimiento. Las ecuaciones de Navier–Stokes convierten esa idea en un modelo continuo. Escribirlas resulta mucho más sencillo que demostrar qué puede hacer cada una de sus soluciones.

Estado de la investigación comprobado: 13 de octubre de 2026. El 8 de septiembre, OpenAI anunció una construcción de una singularidad en tiempo finito para Navier–Stokes incompresible en tres dimensiones con una fuerza externa suave, junto con un artículo y código Lean. El resultado aborda las alternativas C y D del problema del milenio del Instituto Clay. Se trata de un avance importante anunciado, pero no demuestra la formación de singularidades para Navier–Stokes sin fuerza externa. El resultado separado sobre Euler sin fuerza tiene viscosidad nula y no debe confundirse con este. Anuncio de OpenAI, actualizado el 10 de septiembre.

Aquí importa distinguir entre un resultado anunciado, un enunciado formal comprobado, una demostración legible y la aceptación de la comunidad. La sección 14 presenta el estado de la evaluación con fechas, incluido el examen crítico de octubre. La física que sigue explica por qué una energía acotada puede coexistir con un movimiento local cada vez más extremo.


2. Campos continuos y leyes de conservación

Sean la densidad de masa, la velocidad, la presión y el tensor de tensiones de Cauchy. La fuerza de volumen por unidad de masa es ; un ejemplo es la gravedad.

Un campo asigna un valor a cada posición y tiempo, promediando sobre muchas moléculas. El modelo supone que las escalas de longitud relevantes son lo bastante grandes para que este promedio y una ley local de tensiones tengan sentido. Una singularidad matemática se refiere a las ecuaciones del continuo; no predice que una molécula real alcance una velocidad infinita.

2.1 Conservación de la masa

El primer término mide cuánta masa se acumula localmente; la divergencia mide el flujo neto de masa hacia el exterior. Con densidad constante, la conservación de la masa se reduce a que la divergencia de la velocidad sea nula. Una porción de fluido puede estirarse y cambiar de forma conservando su volumen.

2.2 Conservación de la cantidad de movimiento (Cauchy)

Descomponemos la tensión en presión isótropa y parte viscosa:

La derivada material sigue a una porción de fluido en movimiento. Incluso un campo de velocidad estacionario puede acelerarla: el agua que entra en un paso más estrecho llega a una región con otra velocidad. Eso describe el término de transporte no lineal .


3. Ley constitutiva (fluido newtoniano)

Para un fluido newtoniano isótropo,

donde es la viscosidad dinámica y el segundo coeficiente de viscosidad; la viscosidad volumétrica es .

Entre una placa inferior fija y una placa superior móvil, la velocidad laminar estacionaria crece linealmente con la altura. Esto ilustra el esfuerzo cortante viscoso y la condición de no deslizamiento.

Si los coeficientes de viscosidad son constantes en el espacio, al sustituir en el balance de cantidad de movimiento obtenemos Navier–Stokes compresible:

Definimos la viscosidad cinemática cuando es constante.

La viscosidad responde a la deformación, no a la traslación uniforme. Sus dimensiones son y . En el diagrama de flujo entre placas, una separación y una velocidad de la placa superior dan y un esfuerzo cortante . Esta sencilla solución exacta es suave; afirmar la existencia de una singularidad se refiere a ciertos datos particulares, no a todos los flujos.


4. Forma incompresible

Para densidad constante y

obtenemos

La ecuación expresa que la aceleración es igual a la aceleración debida a la presión, más la difusión viscosa y la aceleración aplicada:

TérminoSignificado
Cambio medido en una posición fija
Cambio debido al movimiento de una porción de fluido a través del campo de velocidad
Las fuerzas de presión redistribuyen la cantidad de movimiento
La viscosidad suaviza las variaciones espaciales
Fuerza externa aplicada por unidad de masa

La incompresibilidad no exige velocidad ni aceleración nulas. La presión actúa como una restricción que mantiene nula la divergencia de la velocidad; en este modelo de densidad constante no constituye una ecuación de estado independiente.

Tomando la divergencia de la ecuación de cantidad de movimiento obtenemos una ecuación de Poisson para la presión (con condiciones de contorno adecuadas):


5. Vorticidad y helicidad

Definimos la vorticidad .
Para un flujo incompresible ():

En componentes:

El término es el estiramiento de vórtices (solo en 3D).
La helicidad permite caracterizar la topología de las líneas de vórtice.

Un tubo de vórtice estirado a lo largo de su eje puede hacerse más delgado e intenso. En la ecuación de vorticidad, el estiramiento puede amplificar la rotación mientras la viscosidad la difunde. Ninguno de los dos efectos decide por sí solo si la dinámica acoplada completa permanece suave.

En un flujo verdaderamente plano , la vorticidad apunta perpendicularmente al plano y los campos no varían en esa dirección. Por tanto, . Esto elimina un mecanismo importante de amplificación en dos dimensiones; no significa que todo vórtice tridimensional deba desarrollar una singularidad.

La circulación antihoraria en el plano tiene una vorticidad que apunta hacia fuera de la página, representada mediante un punto dentro de un círculo.


6. Balance de energía (caso incompresible)

Para una solución suave, multiplicamos la ecuación de cantidad de movimiento por e integramos sobre un dominio con una frontera fija de no deslizamiento. La misma identidad se cumple en una caja periódica o en todo el espacio con decaimiento suficiente:

Para , la energía cinética decrece monótonamente; la viscosidad la disipa a una tasa .

Esta energía está normalizada por la densidad constante: la energía cinética física es multiplicado por la integral del cuadrado de la velocidad. El transporte y la presión solo aportan flujos de frontera, que se anulan bajo estas hipótesis; la integración por partes da el término viscoso no negativo. Una pared móvil o una frontera abierta añade términos de trabajo o de flujo.

6.1 Por qué la energía finita no garantiza suavidad

Una cota integral controla el movimiento total, no su mayor valor local. Como ilustración de escala, elegimos un campo fijo , suave, de soporte compacto y de divergencia nula, y definimos

Entonces crece cuando , mientras que disminuye. Son instantáneas reescaladas, no una solución ni una demostración de formación de singularidades. Muestran por qué velocidades cada vez mayores en un volumen cada vez menor no tienen por qué violar una cota de energía. Controlar las derivadas y la concentración es el desafío adicional de la regularidad.


7. Adimensionalización y parámetros de control

Elegimos escalas características (longitud) y (velocidad), y definimos , , , .

Las ecuaciones incompresibles se convierten en

con el número de Reynolds

Aquí . El tiempo de difusión es y el de advección es ; su cociente es . Un número de Reynolds grande indica que la advección es rápida frente al alisamiento viscoso en la escala elegida. No es un umbral universal de turbulencia ni un criterio de singularidad: también importan la geometría, las perturbaciones y la escala medida. Por ejemplo, , y dan .

Otros grupos adimensionales (según el caso):

  • Mach (compresibilidad),
  • Froude (gravedad),
  • Prandtl (difusión de cantidad de movimiento frente a difusión térmica).

8. Navier–Stokes–Fourier compresible (esquema)

Incluimos la energía interna o la temperatura , con el flujo de calor de Fourier (conductividad térmica , distinta de la difusividad térmica anterior):

Una ecuación de estado (por ejemplo, la de un gas ideal) cierra el sistema: .

El trabajo de compresión y la conversión de energía cinética en energía interna aparecen ahora explícitamente. El problema del milenio incompresible no resuelve la formación de ondas de choque, el transporte de calor ni la teoría de regularidad de este sistema diferente.


9. Teoría matemática (existencia y unicidad)

La existencia local significa que hay una solución durante algún intervalo posterior al instante inicial; la existencia global, que la hay para todos los tiempos futuros. Una solución clásica posee suficientes derivadas para satisfacer las ecuaciones punto por punto; aquí, suave significa diferenciable a cualquier orden. Una solución débil satisface una versión integrada, trasladando las derivadas a funciones de prueba suaves. Son criterios de existencia diferentes.

9.1 Soluciones débiles (Leray–Hopf)

Para datos iniciales de divergencia nula y energía finita en el marco habitual de todo el espacio o de condiciones periódicas, y con fuerzas admisibles, la teoría de Leray–Hopf proporciona soluciones débiles globales que satisfacen una desigualdad de energía. Una solución suave es única dentro de la clase de energía correspondiente mientras siga siendo suave (unicidad débil–fuerte). Por tanto, una posible pérdida de suavidad no significa que dejen de existir todas las soluciones generalizadas. Descripción oficial del problema por Fefferman.

La unicidad depende de la clase de soluciones, de los datos iniciales y de la fuerza. Albritton, Brué y Colombo establecen no unicidad para ciertas fuerzas. Hou, Wang y Yang presentan un resultado de no unicidad asistido por ordenador, sin fuerza externa y con datos iniciales singulares. Estos resultados no deben presentarse como formación de singularidades a partir de datos iniciales suaves sin fuerza externa.

9.2 Criterios de regularidad (3D)

  • Criterio de continuación por la vorticidad (tipo Beale–Kato–Majda): para una solución clásica en el marco habitual suficientemente regular de todo el espacio o de condiciones periódicas, con una fuerza suficientemente regular hasta y a través de ese instante, la pérdida de regularidad en requiere

La explicación de Luo de los criterios de continuación de Navier–Stokes enuncia el criterio sin fuerza externa; la condición sobre la fuerza indicada arriba excluye una pérdida de regularidad impuesta por una entrada externa singular.

  • Prodi–Serrin: una solución débil de la clase de energía, con una fuerza adecuada, es regular si con

en cuyo caso la solución es suave en .

  • Caffarelli–Kohn–Nirenberg: las soluciones débiles adecuadas, bajo las hipótesis del teorema sobre la fuerza, tienen un conjunto singular de medida de Hausdorff parabólica unidimensional nula. Es regularidad parcial; no demuestra que el conjunto singular sea vacío. Contexto y referencias en la descripción oficial del problema.

Son afirmaciones condicionales. Especifican qué control impediría la pérdida de regularidad, pero no proporcionan ese control para todo flujo tridimensional. Los exponentes anteriores reflejan el escalamiento natural ; en tres dimensiones, la norma de energía por sí sola no tiene el escalamiento crítico.

9.3 Regularidad global en 2D

Para datos suaves y compatibles en el marco habitual de todo el espacio o de condiciones periódicas, y con una fuerza suficientemente regular, el flujo incompresible bidimensional permanece suave durante todo tiempo finito. Sin fuerza externa, la enstrofía no aumenta; la fuerza añade un término fuente. La ausencia de estiramiento de vórtices permite un control mucho más fuerte que en 3D.

9.4 Flujos compresibles

Existen soluciones débiles globales de energía finita (para gases con ley , bajo ciertas condiciones), pero la regularidad y la unicidad siguen estando en gran medida abiertas.

9.5 Qué piden realmente las cuatro alternativas de Clay

La formulación oficial acepta una demostración de cualquiera de cuatro enunciados, con viscosidad positiva:

AlternativaDominioFuerza y cuestión planteada
AFuerza nula; soluciones suaves globales de energía finita para todos los datos iniciales suaves permitidos
BPeriódico Fuerza nula; soluciones suaves globales para todos los datos iniciales suaves periódicos permitidos
CAlgunos datos suaves permitidos y una fuerza suave de decaimiento rápido no admiten una solución suave global con energía uniformemente acotada
DPeriódico Algunos datos y fuerzas suaves periódicos permitidos no admiten una solución suave global

En todo el espacio, los datos iniciales y la fuerza cumplen las condiciones de decaimiento indicadas; la fuerza periódica también satisface condiciones de decaimiento temporal. C y D permiten una fuerza no nula, de modo que demostrar una de ellas no refutaría lógicamente los enunciados sin fuerza A o B. «Resolver una alternativa de Clay» y «resolver todas las versiones de la regularidad» son afirmaciones diferentes.


10. Turbulencia y espectros

En la turbulencia homogénea e isótropa con alto (K41), el espectro de energía del intervalo inercial es

con tasa media de disipación y constante de Kolmogórov .

Microescala de Kolmogórov:

La estimación K41 da . Resolver esta escala en tres dimensiones sugiere un número de puntos de malla del orden de , antes de contar el coste de avanzar en el tiempo. Esta estimación heurística motiva la simulación numérica directa (DNS), la simulación de grandes remolinos (LES) y las ecuaciones de Navier–Stokes promediadas de Reynolds (RANS), con diferentes decisiones de resolución y modelización. LES y RANS necesitan hipótesis de cierre para las correlaciones no resueltas.

El espectro describe un modelo estadístico del intervalo inercial, no todos los flujos ni todas las escalas. La intermitencia y las fronteras complican esa descripción. Una solución singular construida específicamente no proporcionaría un modelo universal de turbulencia ni eliminaría la necesidad de estimar los errores numéricos.


11. Condiciones de contorno y geometría

  • Pared sin deslizamiento:
  • Deslizamiento libre: ,
  • Caja periódica (toro) para teoría y DNS
  • Entrada y salida con condiciones apropiadas de tensión o presión

La elección afecta a la buena formulación del problema, al balance de energía y a la estabilidad numérica.

Una caja periódica identifica las caras opuestas; no tiene paredes físicas. En todo el espacio, el decaimiento sustituye a una condición de pared en el infinito. Un resultado en estas geometrías no se transfiere automáticamente a la superficie de un avión o a una tubería: las condiciones de contorno forman parte del problema matemático.


12. Métodos numéricos (guía breve)

  • Métodos de proyección: avanzar , resolver Poisson para y proyectar para obtener .
  • Esquemas estables en energía: la forma antisimétrica del término convectivo conserva la energía discreta.
  • Restricción CFL para esquemas explícitos de advección y difusión: .
  • Control de la divergencia: bases espectrales de divergencia nula o elementos finitos mixtos con pares velocidad–presión estables según la condición inf–sup. La divergencia exactamente nula requiere espacios adecuados; la estabilidad inf–sup por sí sola no la garantiza.

Una simulación puede revelar un pico creciente, pero una malla finita no puede establecer por sí sola un límite infinito. Para distinguir una concentración real de una inestabilidad numérica hacen falta estudios de refinamiento, control de errores de discretización y un argumento analítico. Las demostraciones asistidas por ordenador proporcionan cotas rigurosas; los asistentes de demostración formal comprueban deducciones a partir de definiciones e hipótesis codificadas. Ninguno de estos procedimientos consiste simplemente en una animación de alta resolución.


13. Conexiones y extensiones

  • Ecuaciones de Euler: poner elimina la difusión viscosa. El paso al límite puede ser delicado, especialmente cerca de paredes. El resultado sin fuerza sobre Euler anunciado por separado por OpenAI se refiere a esta ecuación sin viscosidad, no al problema sin fuerza y con viscosidad positiva. Anuncio y artículo separado sobre Euler.
  • Hidrodinámica relativista: la formulación ingenua de Navier–Stokes es acausal; las teorías de segundo orden de Israel–Stewart pueden restaurar la causalidad bajo condiciones adecuadas sobre los coeficientes de transporte y los estados del fluido (un ejemplo riguroso con viscosidad volumétrica).
  • Fluidos cuánticos: análogos superfluidos de Navier–Stokes con modelos de dos fluidos y vórtices cuantizados.

14. Estado actual

A 13 de octubre de 2026: se ha anunciado la construcción tridimensional con una fuerza suave y está siendo examinada; la regularidad de Navier–Stokes tridimensional sin fuerza externa sigue abierta. Esto no invalida la regularidad bidimensional ni la utilidad cotidiana de las soluciones de flujo suave.

14.1 El enunciado preciso y su mecanismo

El teorema 1.1 del artículo de OpenAI afirma: para todo , existe una fuerza suave de soporte compacto en espacio y tiempo y un flujo incompresible suave sobre , que parte del reposo, tales que

La velocidad y la presión tienen un soporte espacial compacto fijo antes del tiempo 1; la fuerza se prolonga suavemente a través de ese instante. El artículo afirma resolver la alternativa C y, mediante la construcción periódica del corolario 10.6, la D. Su convención para la presión incorpora la densidad constante.

El mecanismo propuesto concentra un vórtice que gira en espiral hacia dentro y se estira axialmente. Las correcciones oscilatorias y su transporte de cantidad de movimiento cancelan fuerzas residuales que, de otro modo, serían singulares. No basta con elegir cualquier velocidad que diverja y definir una fuerza a partir de la ecuación: la fuerza y todas sus derivadas deben permanecer suaves a través del instante singular. Esta es una descripción conceptual de la construcción anunciada, no una demostración independiente.

«Para toda viscosidad existe una fuerza» permite que la fuerza dependa de la viscosidad. No significa que una misma fuerza supere cualquier viscosidad, que baste una fuerza arbitrariamente pequeña o que un fluido inicialmente en reposo desarrolle una singularidad sin aporte de energía.

14.2 Lean, examen matemático y premio

El repositorio de Lean publicado contiene enunciados objetivo para todo el espacio y para el caso periódico, además de instrucciones para la comprobación independiente con Comparator. Un verificador formal comprueba deducciones para el enunciado y las hipótesis exactos que se han codificado. Si esas definiciones expresan el problema pretendido y cómo el argumento explica la demostración escrita son cuestiones adicionales. Estas notas no afirman haber recompilado ni verificado de forma independiente el proyecto Lean.

El 6 de octubre, Bastounis, Circelli y Hansen publicaron una crítica sobre la correspondencia entre la exposición en lenguaje natural y Lean: los ejemplos incluyen distintos requisitos de derivadas en el lema 8.6 y un argumento diferente para la estimación de la presión en el lema 10.5. Esto cuestiona que se trate la formalización como una certificación automática de cada paso de la demostración escrita. Por sí solo, no refuta el teorema final de formación de singularidades con fuerza ni muestra que el enunciado formal final sea incorrecto.

El 11 de septiembre, la respuesta del Instituto Clay acogió favorablemente el aparente avance, subrayando que su proceso de evaluación es deliberadamente pausado. Las reglas del premio exigen una publicación que cumpla los requisitos, al menos dos años desde esa publicación y aceptación general por la comunidad matemática antes de considerar la concesión. La fecha del anuncio no inicia por sí sola un plazo establecido para el premio. Aquí no se afirma que se haya concedido ningún premio; OpenAI dice que no tiene intención de reclamarlo.

14.3 Qué establece el trabajo posterior y qué no

Los trabajos de septiembre también deben leerse respetando su alcance. Petrillo y Glimm separan explícitamente el resultado con fuerza del problema sin fuerza aún abierto. La exposición de Lei y Ren estudia la construcción del perfil y deja para más adelante la cancelación oscilatoria adicional; no constituye una verificación independiente completa. Constantin, Ignatova y Vicol obtienen un resultado de regularidad que impide la singularidad bajo hipótesis específicas de geometría, crecimiento y fuerza. Es relevante para los intentos de eliminar la fuerza de este mecanismo, pero no es un teorema de regularidad global para todos los flujos con fuerza analítica.

Las preguntas útiles son ahora precisas: ¿se verifican las hipótesis y conclusiones anunciadas?, ¿puede explicarse el mecanismo con claridad?, ¿qué sobrevive al eliminar la fuerza? La predicción de la turbulencia, la unicidad de soluciones débiles en todas las clases y el flujo compresible siguen siendo problemas distintos.

15. Comprueba lo que has aprendido

  1. Traslación uniforme: con velocidad constante, presión constante y fuerza nula, ¿qué términos de la ecuación incompresible se anulan? Todas las derivadas espaciales y temporales se anulan, de modo que cualquier viscosidad positiva es compatible con este movimiento.
  2. Energía frente a velocidad máxima: en la sección 6.1, haz el cambio de variables . ¿Por qué la energía escala como ? El cuadrado de la amplitud aporta y el volumen tridimensional aporta .
  3. Cambio de escala: manteniendo fijos y , ¿qué le ocurre al número de Reynolds si se reduce a la mitad? También se reduce a la mitad; la importancia relativa de la viscosidad aumenta en la escala menor.
  4. Lectura de un teorema: ¿por qué puede cumplirse la alternativa C mientras A sigue sin decidirse? C permite elegir una fuerza suave no nula, mientras que A solo trata de fuerza nula.
  5. Un verificador de demostraciones: ¿qué más hay que comprobar después de que Lean acepte una demostración? Que las definiciones, hipótesis y conclusión final codificadas coincidan con la pregunta matemática pretendida; también hay que revisar la correspondencia con una exposición separada.