Navier-Stokes incompresible versus comprimible

"Incompresible versus comprimible" se reduce a densidad. ¿Se mantiene constante o cambia?

Pruébelo. Llene una jeringa con agua y presione el émbolo. El agua se mueve, pero en condiciones cotidianas no se comprime notablemente. El agua resiste la compresión con tanta fuerza que tratarla como incompresible es una excelente aproximación. Incompresible. Ahora llena esa jeringa con aire y sella el extremo. Empuja el émbolo hacia adentro y sentirás que el aire cede, la misma masa de aire ahora empaquetada en menos volumen a medida que se comprime bajo el pulgar. Eso es flujo comprimible.

En un flujo incompresible, la densidad $\rho$ es constante en todo el fluido, y cada pequeña porción mantiene su volumen a medida que se mueve a través del espacio. El flujo compresible es diferente. La densidad se convierte en una variable, libre de cambiar de un lugar a otro y de un momento a otro. El aire alrededor de un motor a reacción, el gas en una explosión, la atmósfera a gran escala: todo comprimible, todo impulsado por variaciones de densidad.

¿Por qué preocuparse? Porque esta distinción remodela el aspecto de las ecuaciones de Navier-Stokes, lo que predicen y lo difíciles que son de analizar y resolver.

El supuesto de incompresibilidad fija el campo de densidad $\rho$ en una constante positiva en todo el dominio del flujo. Físicamente, esto significa que cada elemento material de volumen preserva su volumen bajo el mapa de flujo. La configuración compresible eleva $\rho(x,t)$ a incógnita completa, gobernada por su propia ecuación de evolución.

Estas no son dos notaciones para el mismo sistema. Difieren en el número de incógnitas, en la estructura de las ecuaciones de restricción y en el carácter de la presión. El sistema incompresible tiene $d+1$ incógnitas escalares ($u$ y $p$ en $\mathbb{R}^d$); el sistema compresible suele tener $d+2$ campos primarios (velocidad, densidad y energía interna o temperatura), con la presión determinada mediante una ecuación de estado.

La distinción también cambia el tipo de EDP de la presión: elíptica en el modelo incompresible, donde la presión se determina por una resolución espacial instantánea, frente a una estructura mixta hiperbólica-parabólica en el caso viscoso compresible. Los modos acústicos de la parte inviscida o hiperbólica se propagan a velocidad finita del sonido, mientras que la viscosidad y la conducción térmica introducen difusión parabólica desde el punto de vista matemático.

Las ecuaciones incompresibles de Navier-Stokes describen fluidos cuya densidad es constante. Son la versión que aparece en el Problema del Milenio de Clay y la versión en la que se centra este sitio.

El sistema tiene dos partes. La ecuación de momento:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f$$

y la restricción de incompresibilidad:

$$\nabla \cdot u = 0$$

La restricción $\nabla \cdot u = 0$ dice que el campo de velocidades es libre de divergencia: el líquido no se acumula ni se diluye en ninguna parte. Todo lo que fluye hacia una pequeña región debe salir al mismo ritmo. Esta única condición reemplaza toda la ecuación de densidad. La densidad no cambia, por lo que no necesitas una ecuación para rastrearla.

La presión juega un papel especial aquí. No está determinado por una ley termodinámica (como la ley de los gases ideales). En cambio, se ajusta instantáneamente en todas partes para mantener el flujo libre de divergencias. Matemáticamente, $p$ resuelve una ecuación de Poisson derivada de la restricción. Los cambios de presión se propagan infinitamente rápido. No hay "velocidad del sonido" en un flujo incompresible.

El sistema incompresible de Navier-Stokes tiene dos campos desconocidos: velocidad $u$ y presión $p$. Esa simplicidad es engañosa. El término no lineal $(u \cdot \nabla)u$ todavía hace que el sistema sea extremadamente difícil en tres dimensiones.

El sistema incompresible de Navier-Stokes en $\mathbb{R}^3$ con viscosidad cinemática $\nu > 0$:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f, \qquad x \in \mathbb{R}^3,\; t > 0,$$

$$\nabla \cdot u = 0, \qquad u(x,0) = u_0(x).$$

La condición libre de divergencia $\nabla \cdot u = 0$ es una restricción puntual, no una ecuación de evolución. Codifica la preservación del volumen local: el mapa de flujo $\Phi_t$ satisface $\det(D\Phi_t) = 1$ para todo $t$.

Aplicar el operador de divergencia a la ecuación de momento y usar la incompresibilidad produce la ecuación de Poisson de presión:

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

Esta es una ecuación elíptica para $p$ en cada tiempo fijo. La presión no es una variable termodinámica independiente; en la interpretación estándar de EDP, actúa como un multiplicador de Lagrange que impone la restricción libre de divergencia, determinada global e instantáneamente por el campo de velocidades. La información se propaga a velocidad infinita a través de la presión, una diferencia estructural del sistema comprimible que no se puede ocultar.

Las incógnitas son $u : \mathbb{R}^3 \times [0,T) \to \mathbb{R}^3$ y $p : \mathbb{R}^3 \times [0,T) \to \mathbb{R}$. Para la formulación de Clay (Fefferman, 2000), $u_0 \in C^\infty(\mathbb{R}^3)$ está libre de divergencia y la pregunta es si $u$ permanece en $C^\infty(\mathbb{R}^3 \times [0,\infty))$ con energía acotada.

Las ecuaciones compresibles de Navier-Stokes gobiernan flujos en los que la densidad varía. Sistema más grande. Más incógnitas. Más ecuaciones.

Sigues teniendo una ecuación de cantidad de movimiento, pero ahora la densidad $\rho$ aparece explícitamente:

$$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) = -\nabla p + \nabla \cdot \tau + \rho f$$

La restricción $\nabla \cdot u = 0$ desaparece. En su lugar aparece una ecuación de continuidad que sigue la evolución de la densidad:

$$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

Esto dice que la masa se conserva: la densidad cambia porque el flujo comprime o expande parcelas de fluido.

El sistema también necesita una ecuación de energía y una ecuación de estado, una relación termodinámica como $p = \rho R T$ (la ley de los gases ideales) que vincula presión, densidad y temperatura. La presión ya no es un mero ejecutor pasivo de una restricción. Su parte acústica viaja a la velocidad del sonido, mientras que el sistema viscoso completo también contiene efectos difusivos.

El sistema compresible es esencial para la aerodinámica a altas velocidades, la dinámica astrofísica de gases, la combustión y cualquier flujo donde importen los cambios de densidad. Pero es un objeto matemático genuinamente distinto del sistema incompresible: más incógnitas, más ecuaciones y una estructura de EDP enteramente diferente.

El sistema compresible de Navier-Stokes acopla la velocidad $u(x,t)$, la densidad $\rho(x,t)$, la presión $p(x,t)$ y la energía interna específica $e(x,t)$ (o la temperatura $\theta$). En forma de conservación:

Continuidad: $$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

Momento: $$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) + \nabla p = \nabla \cdot \tau + \rho f$$

Energía: $$\partial_t (\rho E) + \nabla \cdot ((\rho E + p)u) = \nabla \cdot (\tau \cdot u) + \nabla \cdot (\kappa \nabla \theta) + \rho f \cdot u$$

donde $E = e + \tfrac{1}{2}|u|^2$ es la energía específica total, $\tau$ es el tensor viscoso de tensiones (para un fluido newtoniano, $\tau = \mu(\nabla u + \nabla u^T) + \lambda (\nabla \cdot u)I$ con viscosidad volumétrica $\lambda$), y $\kappa$ es la conductividad térmica.

El cierre requiere una ecuación de estado, por ejemplo $p = (\gamma - 1)\rho e$ para un gas ideal con índice adiabático $\gamma$.

En la parte inviscida o hiperbólica del sistema compresible, las perturbaciones acústicas se propagan a velocidad finita, con velocidad característica del sonido $c = \sqrt{\partial p / \partial \rho |_s}$. El sistema viscoso completo de Navier-Stokes-Fourier es mixto hiperbólico-parabólico, así que la difusión viscosa y térmica introduce propagación infinita desde el punto de vista matemático. Las ecuaciones compresibles admiten ondas de choque, rarefacciones y discontinuidades de contacto que no tienen análogo en el flujo incompresible.

¿Cuándo importa la compresibilidad? Un número decide: el número de Mach.

$$\text{Ma} = \frac{|u|}{c}$$

$|u|$ es la velocidad del flujo. $c$ es la velocidad del sonido. Su relación indica qué tan rápido se mueve el flujo en comparación con la velocidad a la que las perturbaciones de presión pueden propagarse a través del medio, y esa comparación determina si se pueden ignorar con seguridad los cambios de densidad o si dominarán la física.

Cuando $\text{Ma} < 0.3$, la densidad cambia en menos de aproximadamente un 5 %. Las ecuaciones incompresibles funcionan. Aire en una habitación, agua en una tubería, viento alrededor de un edificio: todos flujos de bajo Mach donde las perturbaciones de presión viajan mucho más rápido que el propio flujo que la densidad apenas se mueve.

Por encima de $\text{Ma} \approx 0.3$, la compresibilidad comienza a afectar, y alrededor de $\text{Ma} \approx 1$ se llega al régimen transónico donde aparecen bolsas supersónicas locales y se forman ondas de choque. Aviones de combate. Boquillas de cohetes. Reingresando a la nave espacial.

No es un interruptor binario. La mayoría de los flujos de fluidos cotidianos, y el Problema del Milenio de Clay, se asientan firmemente en el régimen de bajo Mach donde se aplican las ecuaciones incompresibles.

El número de Mach $\text{Ma} = |u|/c$ parametriza la importancia de la compresibilidad, donde $c = \sqrt{\partial p / \partial \rho |_s}$ es la velocidad isentrópica del sonido. Formalmente, las ecuaciones incompresibles surgen como el límite bajo de Mach del sistema compresible.

La expansión asintótica en potencias de $\text{Ma}^2$ (ver Klainerman & Majda, 1981, 1982; Schochet, 1986) muestra que cuando $\text{Ma} \to 0$ con datos iniciales adecuados, las soluciones comprimibles convergen a la solución incompresible. En no dimensionalizaciones estándar de bajo Mach con datos bien preparados, a menudo se escribe la presión como un fondo termodinámico casi espacialmente uniforme más una corrección dinámica más pequeña que impone la restricción de incompresibilidad en el límite.

Este es un límite singular: la velocidad del sonido $c \to \infty$ y el carácter hiperbólico del sistema compresible degenera a la ecuación de presión elíptica del flujo incompresible. Los modos acústicos se vuelven infinitamente rápidos y se desacoplan de la dinámica del vórtice.

El régimen $\text{Ma} < 0.3$ es una heurística de ingeniería. Refleja la observación empírica de que la variación de la densidad relativa $\delta\rho / \rho \sim \text{Ma}^2 / 2$ se mantiene por debajo de $\sim$5% en este rango. La justificación matemática es el teorema de convergencia para el límite bajo de Mach, que requiere datos iniciales bien preparados y condiciones de contorno compatibles.

El Problema del Milenio de Clay plantea una pregunta precisa: dada una velocidad inicial suave y sin divergencias en $\mathbb{R}^3$, ¿el sistema incompresible de Navier-Stokes siempre produce una solución suave que existe para siempre?

¿Por qué incompresible específicamente? Tres razones.

Primero, ya es bastante difícil. Las ecuaciones 3D incompresibles se han resistido a la prueba de regularidad global desde el trabajo fundacional de Leray en 1934. Agregar densidad variable, termodinámica y ondas de choque haría que el problema fuera mucho más difícil, no más manejable.

En segundo lugar, la dificultad es pura mecánica de fluidos. El sistema incompresible aísla el desafío matemático central, la competencia entre la advección no lineal $(u \cdot \nabla)u$ y la disipación viscosa $\nu \Delta u$, sin complicaciones termodinámicas ni acústicas. Es el ámbito más limpio para plantear la cuestión de la regularidad.

En tercer lugar, la física está limpia. Las ecuaciones incompresibles modelan los flujos cotidianos más comunes. Si pueden producir singularidades a partir de datos suaves es una cuestión fundamental sobre la consistencia matemática de la mecánica de fluidos clásica.

El sistema compresible tiene sus propios problemas profundamente abiertos (existencia de soluciones globales con grandes datos, formación e interacción de choques), pero esos son problemas diferentes con estructuras diferentes. El premio Clay se centra en el caso incompresible porque esa es la pregunta específica sobre regularidad que Fefferman formuló para 3D Navier-Stokes.

La formulación oficial de Clay (Fefferman, 2000) especifica el sistema incompresible en $\mathbb{R}^3$:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0, \qquad u|_{t=0} = u_0,$$

con $u_0 \in C^\infty(\mathbb{R}^3)$ libre de divergencia, y la pregunta es si $u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$ con $\int_{\mathbb{R}^3} |u(x,t)|^2\,dx$ acotado para todo $t \geq 0$.

La elección del sistema incompresible está motivada matemáticamente. El problema abierto clave, la brecha entre las soluciones débiles de Leray-Hopf (que existen globalmente pero pueden no ser únicas o suaves) y las soluciones suaves clásicas (que existen localmente pero pueden explotar), es específica de las ecuaciones 3D incompresibles. En 2D, se conoce la regularidad global para soluciones suaves e incompresibles de Navier-Stokes; el caso 3D permanece abierto.

El sistema compresible introduce dificultades cualitativamente diferentes: formación de choques (que ocurre incluso para ecuaciones de Euler con datos suaves), estados de vacío ($\rho \to 0$) y el acoplamiento entre la vorticidad y los modos acústicos. Estos son problemas abiertos importantes, pero son estructuralmente distintos de la cuestión de la regularidad incompresible.

El problema de la incompresibilidad aísla la competencia entre la no linealidad de energía supercrítica y la disipación viscosa. La estimación de energía natural da $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$, que no alcanza el espacio crítico de escala $L^\infty_t \dot{H}^{1/2}_x$ por media derivada en 3D. Cerrar esta brecha, o demostrar que no se puede cerrar, es el meollo del Problema del Milenio.

Empiece aquí. ¿Quiere separar todos los términos del sistema incompresible y explicar desde cero el significado físico y la función matemática de cada pieza? ¿Qué son las ecuaciones de Navier-Stokes?

¿De dónde viene este sistema? Derivación de las ecuaciones de Navier-Stokes.

Si reduce la viscosidad, obtendrá las ecuaciones de Euler, que son un siglo más antiguas, parecen más simples en la página y, en algunos aspectos, son incluso más difíciles de entender matemáticamente porque se pierde el efecto suavizante del término de difusión. Euler vs. Navier-Stokes.

El premio. El problema de existencia y suavidad de Navier-Stokes.

Rutas sugeridas desde aquí:

• ¿Qué son las ecuaciones de Navier-Stokes?, el sistema incompresible completo con análisis término por término
• Derivación de las ecuaciones de Navier-Stokes, a partir de la mecánica del continuo axiomas del sistema EDP
• Euler vs. Navier-Stokes, el límite invisible $\nu \to 0$ y el papel de la viscosidad en la regularidad
• El problema de existencia y suavidad de Navier-Stokes, la formulación precisa de Clay, la teoría de Leray-Hopf y la brecha entre soluciones débiles y fuertes