Compressible and immiscible fluids with arbitrary density ratio

Este artículo presenta una nueva teoría basada en la termodinámica para modelar fluidos compresibles e inmiscibles con relaciones de densidad arbitrariamente altas, abordando las limitaciones de las ecuaciones tradicionales de Navier-Stokes y Euler mediante la derivación de ecuaciones de evolución de la densidad que consideran adecuadamente la verdadera evolución del momento.

Lo esencial

Imagina que estás intentando describir una danza entre dos compañeros muy diferentes: un elefante pesado y de movimiento lento (agua) y una pluma ligera y de movimiento rápido (aire). En el mundo de la física de fluidos, estos son "inmiscibles", lo que significa que no se mezclan como el aceite y el agua, pero sí tienen un límite difuso donde se encuentran.

Durante mucho tiempo, los científicos han luchado por escribir las "reglas de la danza" (las ecuaciones) para sistemas donde la diferencia de densidad es enorme, como cuando uno es 1,000 veces más pesado que el otro, como el elefante frente a la pluma. Las reglas antiguas tenían un fallo importante: trataban el peso del elefante como si fuera un número constante e inalterable, incluso justo en el límite difuso donde el elefante se convierte en pluma. Esto es como intentar describir a una persona convirtiéndose en un fantasma diciendo: "Sigue siendo un humano sólido todo el tiempo", lo cual no tiene sentido.

Aquí hay un desglose sencillo de lo que este artículo propone solucionar.

1. El problema con las reglas antiguas

La forma tradicional de calcular cómo se mueven los fluidos (usando las ecuaciones de Navier-Stokes y Euler) se basa en un atajo llamado aproximación de Boussinesq.

• La analogía: Imagina que estás empujando un carrito de la compra. Si el carrito está lleno de ladrillos, es pesado. Si está vacío, es ligero. Las reglas antiguas asumían que, si estás empujando un carrito que está parcialmente lleno de ladrillos y parcialmente lleno de aire, el peso del carrito nunca cambia mientras lo empujas. Simplemente asume que el peso es un promedio fijo.
• El fallo: En la realidad, a medida que el carrito se mueve y los ladrillos se desplazan (difusión), el peso cambia. Las reglas antiguas ignoraban el hecho de que el "momento" (masa $\times$ velocidad) cambia porque la masa misma está cambiando en el límite. También asumieron que el aire y el agua ocupan una cantidad de espacio perfectamente predecible, ignorando que cuando se mezclan en el límite, el espacio que ocupan puede en realidad oscilar y cambiar debido a la presión.

2. El nuevo enfoque: Minimización de la energía

En lugar de adivinar cómo cambia la densidad, el autor parte de un principio fundamental: la naturaleza siempre intenta usar la menor cantidad de energía posible.

• La analogía: Piensa en una pelota rodando por una colina. No le importa los pasos específicos que da; solo quiere llegar al fondo (menor energía). El autor utiliza este concepto de la "colina de energía" para derivar nuevas reglas sobre cómo interactúan el agua y el aire.
• La innovación clave: El autor introduce un concepto llamado "Volumen Excedente".
  • Imagina que tienes un cubo de agua y un cubo de aire. Si los viertes juntos, podrías esperar que el volumen total sea exactamente la suma de los dos cubos. Pero a nivel microscópico, cuando se encuentran, las moléculas pueden empaquetarse más apretadas o más sueltas, creando un poco de espacio "extra" o "faltante".
  • Las reglas antiguas asumían que este espacio extra era cero en todas partes. Este artículo dice: "No, ese espacio extra existe, cambia de un lugar a otro y afecta a la densidad".

3. Las nuevas "Reglas de la Danza" (Los resultados)

Al contabilizar este "espacio extra" cambiante y utilizar la minimización de la energía, el autor deriva un nuevo conjunto de ecuaciones que hacen tres cosas principales:

A. Una nueva forma de escuchar el sonido (Velocidad del sonido)
El artículo muestra que la velocidad del sonido no es solo un número aleatorio; proviene directamente de cómo cambia la energía del fluido cuando se comprime.

• La metáfora: Piensa en el sonido como una onda en una multitud. La velocidad de esa onda depende de qué tan apretadas estén las personas (moléculas) y cuánta energía tengan. La nueva fórmula calcula esta velocidad de forma natural, sin necesidad de que se le indique de antemano. Incluso sugiere que, en un gas, la velocidad del sonido es aproximadamente la misma que la velocidad promedio a la que rebotan las moléculas del gas.

B. Una nueva regla para la presión y la velocidad (Ley de Bernoulli)
Probablemente hayas oído hablar del principio de Bernoulli: "Cuando un fluido se mueve más rápido, su presión cae".

• El giro: La regla antigua funciona muy bien para el agua fluyendo en una tubería, pero falla cuando hay un salto enorme en la densidad (como el agua chocando con el aire). El autor crea una Ley de Bernoulli Generalizada.
• La metáfora: Imagina un río que fluye hacia una cascada. La regla antigua dice que la energía se mantiene igual. La nueva regla dice: "Espera, a medida que el agua se convierte en niebla (aire), parte de la energía se pierde o se transforma porque la 'densidad de ocupación' del agua está cambiando". La nueva ecuación tiene en cuenta este cambio de energía, lo que la hace precisa incluso cuando el fluido está cambiando su naturaleza de pesado a ligero.

C. El "bulto" en la densidad
Este es quizás el resultado más visual.

• La visión antigua: Si miras el límite entre el agua y el aire, los modelos antiguos decían que la densidad simplemente bajaría suavemente de "agua pesada" a "aire ligero", como una rampa descendente.
• La nueva visión: La matemática del autor predice un bulto. Al cruzar el límite, la densidad en realidad sube ligeramente antes de caer al nivel del aire.
• La metáfora: Imagina a una multitud de personas (agua) intentando pasar por una puerta hacia un pasillo vacío (aire). Mientras se amontonan en el umbral, podrían apretarse un poco más durante un breve instante antes de dispersarse. La nueva teoría predice este "bulto de empaquetamiento", lo cual coincide con lo que las simulaciones computacionales avanzadas (llamadas Teoría del Funcional de la Densidad) ya han visto, pero los modelos simples antiguos pasaron por alto.

Resumen

Este artículo propone una nueva forma de escribir las leyes de la física para fluidos que son muy diferentes en peso (como el agua y el aire).

1. Deja de pretender que el peso es constante en el límite.
2. Admite que el espacio que ocupan las moléculas cambia (volumen excedente).
3. Utiliza el principio de "menor energía" para derivar nuevas reglas que explican cómo viaja el sonido, cómo cambia la presión y por qué la densidad en el límite agua-aire tiene un pequeño "bulto" en ella.

El autor afirma que este nuevo marco de trabajo funciona para cualquier mezcla de fluidos, sin importar cuán diferentes sean sus pesos, abriendo la puerta a simulaciones computacionales más precisas de cosas como la caída de la lluvia, el choque de las olas o el ascenso de las burbujas.

Resumen técnico

Resumen Técnico: Fluidos Compresibles e Inmiscibles con Razones de Densidad Arbitrarias: Teoría de la Interfaz Difusa

Planteamiento del Problema
Los modelos actuales de dinámica de fluidos computacional (CFD) tienen dificultades para simular con precisión sistemas de fluidos inmiscibles con altas razones de densidad, como el agua-aire (razón $\approx 1000$). Los enfoques existentes dependen de dos métodos principales, ambos de los cuales poseen limitaciones significativas:

1. La Aproximación de Boussinesq: Asume una densidad constante ($\rho$) dentro de la interfaz, utilizando la ecuación $\rho(du/dt) = -\nabla p + \dots$. Esta formulación no logra contabilizar la "verdadera" evolución del momento $d(\rho u)/dt$, ignorando los cambios de densidad y la resultante transferencia de momento. Además, no puede generar naturalmente fuerzas de gravedad y flotabilidad sin adiciones ad hoc, ya que estas surgen de las diferencias de densidad.
2. El Método de Interpolación: Define la densidad como una función lineal de la concentración volumétrica ($\rho = \rho_w \phi + \rho_a(1-\phi)$). Este enfoque crea un conflicto fundamental entre la incompresibilidad ($d\rho/dt = 0$) y la difusión ($d\phi/dt \neq 0$). Asimismo, asume un exceso de volumen uniforme (a menudo cero) en todo el dominio, lo cual contradice la realidad física donde el exceso de volumen local depende de la presión local. Por consiguiente, las ecuaciones de continuidad estándar y la interpolación lineal no logran capturar los perfiles de densidad no monotónicos observados en la interfaz agua-aire.

Metodología
Los autores proponen un nuevo marco teórico basado en el principio termodinámico de minimización de la energía, distinto de la derivación tradicional de las ecuaciones de Navier-Stokes y Euler. La metodología implica:

• Reformulación del Volumen y la Densidad: En lugar de asumir un exceso de volumen de cero, el modelo introduce un exceso de volumen espacialmente variable ($v_e$) dentro de la interfaz. Este volumen se redistribuye entre los componentes del fluido, lo que conduce a densidades parciales ($\rho'_w, \rho'_a$) que ya no son constantes, sino que dependen de la presión local y la mezcla. Esto requiere dos ecuaciones de evolución acopladas: una para la concentración de volumen ($\phi$) y otra para la densidad ($\rho$).
• Principio de Disipación-Conservación: Los autores establecen un teorema que vincula la primera ley de la termodinámica (conservación) con la segunda ley (disipación). Al definir la energía total del sistema ($E$) como la suma de la energía libre ($F$), la energía de presión ($P$), la energía cinética ($K$) y la energía libre de pared ($W$), derivan las ecuaciones cinéticas minimizando la disipación de energía.
• Relación de Gibbs-Duhem Modificada: Para resolver las incompatibilidades matemáticas entre los gradientes espaciales y las derivadas funcionales en los modelos existentes, se deriva una relación de Gibbs-Duhem isotérmica modificada. Esta explica explícitamente los cambios de entropía dependientes de la velocidad y el trabajo inercial, lo que conduce a una definición rigurosa de la presión que incluye contribuciones de la energía libre y la entropía de la velocidad.
• Definición de la Energía Cinética: La energía cinética se formula como $K = \int \rho u \cdot u \, d\Omega$ (densidad de momento lineal $\zeta = \rho u$) en lugar del convencional $\frac{1}{2}\rho u^2$. Esta distinción permite la separación del balance de momento (cambiando tanto $u$ como $\rho$) de la disipación de velocidad (cambiando solo $u$).

Contribuciones Clave
La derivación produce una ecuación de evolución de densidad generalizada y tres avances teóricos principales:

1. Ley de Bernoulli Generalizada: Se demuestra que el principio de Bernoulli estándar ($p + \frac{1}{2}\rho u^2 = \text{constante}$) es un caso especial. La forma generalizada tiene en cuenta los cambios de la entropía de la velocidad y la evolución de la densidad en la interfaz: $p + \rho u^2 = \text{constante}$ (bajo condiciones específicas de estado estacionario e invíscido), reconociendo que ocurren saltos de energía a través de las interfaces debido a las diferencias de densidad.
2. Formulación Universal de la Velocidad del Sonido: Se deriva una expresión general para la velocidad del sonido basada en la derivada de la energía libre con respecto a la densidad: $c = \sqrt{d\Psi/d\rho}$. Para gases ideales, esto aproxima la velocidad térmica media de las moléculas.
3. Ecuación de Estado (EOS) Generalizada: La teoría demuestra que la ecuación de evolución de la densidad genera naturalmente una EOS para los fluidos (por ejemplo, $p \sim \rho \ln \rho$ para gases ideales) sin requerir una EOS prescrita, contrastando con los métodos de interfaz difusa estándar.

Resultados y Validación
El artículo valida la teoría a través de tres casos de estudio:

• Propagación del Sonido: La ecuación de la velocidad del sonido derivada es consistente con el análisis de perturbación y predice correctamente la relación de dispersión para las ondas sonoras en el aire.
• EOS de Gas Ideal: En el estado estacionario, el modelo reproduce la ley de los gases ideales, confirmando que la ecuación de evolución de la densidad captura intrínsecamente las relaciones de estado termodinámicas.
• Dinámica de la Interfaz Agua-Aire: Una solución analítica para una interfaz plana de agua-aire en movimiento en un tubo demuestra que el modelo predice un perfil de densidad no monotónico. Específicamente, la densidad aumenta desde el valor del agua en masa hasta un pico máximo en el centro de la interfaz antes de disminuir al valor del gas en masa. Este resultado se alinea con los hallazgos de la Teoría del Funcional de la Densidad (DFT) [11, 12] y contradice la disminución monotónica predicha por los métodos de interpolación lineal. El modelo también reproduce correctamente la ecuación de Young-Laplace para interfaces curvas.

Significado
El artículo afirma proporcionar un marco general para fluidos inmiscibles con razones de densidad arbitrariamente altas, aplicable tanto a regímenes compresibles como incompresibles. Al tratar la evolución de la densidad como una variable dinámica acoplada con el exceso de volumen y la disipación viscosa, el modelo resuelve las contradicciones inherentes a la aproximación de Boussinesq y a los métodos de interpolación lineal. El trabajo sugiere que la disipación viscosa altera intrínsecamente la densidad, desafiando la estricta suposición de incompresibilidad en flujos disipativos. Los autores postulan que esta teoría abre una nueva ventana para la CFD, ofreciendo una descripción físicamente consistente de la transferencia de momento y energía en las interfaces fluido-fluido donde las razones de densidad son extremas.