Lattice Boltzmann model for non-ideal compressible fluid dynamics

Este trabajo presenta un modelo novedoso de Boltzmann en retículo que utiliza retículos de primeros vecinos y términos de corrección cuasi-equilibrio para simular flujos de fluidos compresibles no ideales con consistencia termodinámica y estabilidad numérica, validando con éxito su precisión mediante simulaciones cuantitativas de interacciones entre ondas de choque y gotas a números de Mach de hasta 1.47.

Lo esencial

Imagina que estás intentando simular cómo se comporta un fluido en una computadora. Durante mucho tiempo, las computadoras han sido excelentes simulando fluidos "ideales", como el agua que fluye suavemente en un río o el aire que se mueve lentamente alrededor de un ala. Estos fluidos siguen reglas simples y predecibles.

Pero, ¿qué sucede cuando el fluido está bajo presión y calor extremos, comportándose como un gas denso que casi es un líquido, o un líquido que casi es un gas? Este es el mundo de los fluidos compresibles no ideales. Piénsalo como un fluido que está "estresado" y actúa de manera extraña, negándose a seguir las reglas simples del mundo ideal. Esto ocurre en tecnologías avanzadas como turbinas de CO2 supercrítico y ciclos de energía orgánicos.

El problema es que las herramientas informáticas existentes luchan con estos fluidos estresados. O bien se bloquean, dan respuestas incorrectas o requieren tanta potencia de cálculo que se vuelven inútiles.

Este artículo presenta una nueva y más inteligente forma de simular estos fluidos complicados utilizando un método llamado Método de Boltzmann en Red (LBM). Así es como funciona el nuevo enfoque de los autores, explicado mediante analogías simples:

1. El sistema de "dos carriles"

La mayoría de los antiguos métodos de simulación intentan rastrear todo (masa, velocidad, energía) con un único y complicado conjunto de reglas. Los autores se dieron cuenta de que esto era como intentar conducir un coche mientras simultáneamente haces malabares con una docena de pelotas: se vuelve desordenado e inestable.

En su lugar, construyeron un sistema de dos carriles:

• Carril A (La multitud): Un conjunto de reglas rastrea la densidad y la velocidad del fluido (cuántas partículas hay y hacia dónde van).
• Carril B (La energía): Un segundo conjunto de reglas, separado, rastrea la energía total.
  Al separar estos elementos, la computadora no se confunde. Es como tener un controlador de tráfico dedicado para los coches y otro separado para el combustible, asegurando que ningún sistema haga colapsar al otro.

2. El atractor de "cuasi-equilibrio"

En física, los fluidos naturalmente quieren asentarse en un estado tranquilo llamado "equilibrio". Sin embargo, en estas condiciones extremas, el fluido está siendo constantemente empujado y tirado, por lo que nunca llega a asentarse del todo.

Los autores inventaron un truco inteligente llamado atractor de "cuasi-equilibrio".

• La analogía: Imagina a un perro persiguiendo una pelota. La pelota representa el "estado de calma perfecto". El perro representa el fluido. En una situación normal, el perro corre directamente hacia la pelota.
• El problema: En este fluido extremo, la pelota sigue alejándose o cambiando de forma. Si el perro persigue la pelota a ciegas, podría correr hacia un precipicio (la simulación se vuelve inestable).
• La solución: Los autores le dieron al perro un "GPS" que predice dónde estará la pelota dentro de un instante, basándose en cómo cambian el viento (presión) y el terreno (densidad). Este objetivo "desplazado" permite que el perro corra suavemente sin caer por el precipicio. Esto asegura que la simulación se mantenga estable incluso cuando el fluido se mueve muy rápido o cambia de densidad rápidamente.

3. Arreglando el calor "espurio"

Cuando los fluidos se mueven rápido, generan calor. En los modelos informáticos estándar, el calor a veces fluye en la dirección incorrecta o crea calor "fantasma" falso que no existe en la realidad.

• La analogía: Es como un termostato que piensa que la habitación está congelada porque está midiendo la corriente de aire de una ventana, no la temperatura real de la habitación.
• La solución: Los autores añadieron un "término de corrección" específico a sus ecuaciones. Esto actúa como un filtro que elimina las corrientes falsas, asegurando que el calor fluya exactamente como dicta la física (Ley de Fourier), incluso en estas condiciones extremas y no ideales.

4. Las pruebas de "tubo de choque" y "columna líquida"

Para demostrar que su nuevo método funciona, no solo hicieron matemáticas; realizaron pruebas extremas:

• El tubo de choque: Simularon una explosión repentina de presión (una onda de choque) moviéndose a través de un gas denso. En los gases normales, estas ondas se comportan de una manera. En estos gases "no ideales", las ondas pueden hacer algo extraño: un "choque de rarefacción" (una onda que se expande pero aún actúa como un choque agudo). Su modelo predijo con éxito este comportamiento extraño, que los modelos más antiguos pasaron por alto.
• La columna líquida: Simularon una onda de choque de alta velocidad golpeando una gota de líquido. Esta es una prueba muy difícil porque el choque rebota, se refleja y rompe la gota en pedazos. Su modelo manejó la colisión perfectamente, coincidiendo con experimentos del mundo real donde la gota de líquido se aplana y se expande exactamente como debería.

Por qué esto es importante

Los autores afirman que su método es rápido, estable y preciso. Utiliza una cuadrícula simple (como un tablero de ajedrez estándar) en lugar de necesitar una cuadrícula supercompleja y estirada. Esto significa que los científicos ahora pueden simular estos flujos de fluidos no ideales, extremos y de alta velocidad en computadoras estándar con alta precisión.

En resumen: El artículo presenta un nuevo "manual de conducción" para las simulaciones por computadora que les permite manejar fluidos bajo estrés extremo sin bloquearse. Al utilizar un sistema de dos carriles y un "GPS" inteligente para la energía del fluido, pueden predecir con precisión cómo se comportan estos fluidos complejos en escenarios de alta velocidad, abriendo la puerta a mejores diseños para sistemas de energía avanzados.

Resumen técnico

Resumen Técnico: Modelo de Boltzmann en Red para la Dinámica de Fluidos Compresibles No Ideales

Enunciado del Problema
La dinámica de fluidos compresibles no ideales (NICFD) aborda los comportamientos de los fluidos en regímenes cercanos, transcríticos y supercríticos, donde los fluidos se desvían significativamente de las leyes de los gases ideales. Estos regímenes son críticos para las tecnologías energéticas modernas, como los ciclos Rankine orgánicos y las turbinas de CO₂ supercrítico. Sin embargo, los datos experimentales en estos estados son escasos y difíciles de obtener. Aunque la Simulación Numérica Directa (DNS) es necesaria para comprender la física compleja (por ejemplo, estructuras de choque, turbulencia en vapores densos) y generar bases de datos de ingeniería, las herramientas numéricas existentes a menudo luchan con la consistencia termodinámica, la estabilidad y la recuperación precisa de los límites hidrodinámicos para ecuaciones de estado (EOS) no ideales. Específicamente, los métodos convencionales de Dinámica de Fluidos Computacional (CFD) enfrentan desafíos en la región espínodal termodinámicamente inestable, donde el cuadrado de la velocidad del sonido adiabática se vuelve negativo, requiriendo tratamientos especiales como la regla de áreas iguales de Maxwell. El Método de Boltzmann en Red (LBM), aunque exitoso para flujos incompresibles y multifásicos, ha ignorado históricamente la termodinámica compresible no ideal en sus formulaciones cinéticas, dependiendo a menudo de regularizaciones ad hoc o plantillas extendidas que comprometen la eficiencia o la invariancia galileana.

Metodología
Los autores proponen un modelo cinético novedoso y su realización en Lattice Boltzmann diseñada para recuperar las ecuaciones de Navier–Stokes–Korteweg (NSK) compresibles para fluidos no ideales. La metodología se basa en los siguientes pilares:

1. Enfoque de Función de Doble Distribución (DDF): El modelo emplea dos funciones de distribución: $f$ para masa y momento, y $g$ para energía volumétrica ($\rho E$). Este desacoplamiento permite el uso de redes clásicas de primeros vecinos (específicamente D3Q27) sin requerir plantillas extendidas.
2. Atractores de Cuasi-Equilibrio: El operador de colisión utiliza una estructura tipo BGK de doble relajación que involucra un atractor de equilibrio local ($f^{eq}, g^{eq}$) y un atractor de cuasi-equilibrio desplazado ($f^\star_\lambda, g^\star_\lambda$).
   • El equilibrio local se define utilizando un parámetro de temperatura de referencia $\theta = P/\rho$ (trabajo termodinámico del flujo) en lugar de la temperatura termodinámica $T$, asegurando el ancho correcto de la distribución de Maxwell independiente de la energía interna.
   • El equilibrio desplazado introduce correcciones para la fuerza de Korteweg (tensión superficial), la viscosidad volumétrica y el flujo de calor. Emplea valores desplazados para la velocidad del flujo ($u^\star$), la temperatura de referencia ($\theta^\star$) y la temperatura termodinámica ($T^\star$).
3. Consistencia Termodinámica:
   • La temperatura de referencia $\theta$ se establece en $P/\rho$, una formulación que asegura que el modelo permanezca válido para EOS no ideales (por ejemplo, van der Waals).
   • Se introduce un término de corrección en el flujo de calor ($q^c_\lambda$) para compensar los componentes espurios no de Fourier que surgen de los gradientes de densidad en fluidos no ideales, asegurando la recuperación de la ley de Fourier.
   • El desplazamiento en la temperatura de referencia permite que la viscosidad volumétrica sea un parámetro independiente, ajustable y definido positivo, evitando el problema de la viscosidad volumétrica negativa que puede surgir en los modelos BGK estándar para presiones no ideales.
4. Discretización: El modelo utiliza un esquema explícito de segundo orden derivado de la integración trapezoidal a lo largo de las características. Emplea una red estándar D3Q27 y maneja los términos fuente (fuerza de Korteweg) a través del equilibrio desplazado, evitando errores de discretización explícita de los términos fuente.

Contribuciones Clave

• Nuevo Marco Cinético: El artículo introduce un modelo cinético que recupera las ecuaciones hidrodinámicas objetivo (balance de masa, momento y energía volumétrica) para fluidos compresibles no ideales utilizando una plantilla cinética mínima (primeros vecinos).
• Disipación Definida Positivamente: Al construir el modelo mediante atractores de cuasi-equilibrio, los autores aseguran que las tasas de disipación de Navier–Stokes (viscosidad de corte y volumétrica) sean definidas positivamente e invariantes galileanas, un punto de falla común en intentos anteriores de LBM para flujos compresibles no ideales.
• Robustez Termodinámica: El modelo maneja naturalmente las regiones de inestabilidad termodinámica (espínodal) donde la velocidad del sonido se vuelve imaginaria en el límite inviscido, ya que el LBM hereda la propagación a lo largo de características discretas fijas en lugar de depender de solucionadores de EDP hiperbólicas que requieren tratamientos especiales en estas regiones.
• Validación de Fenómenos No Clásicos: El estudio proporciona la primera validación cuantitativa del LBM para interacciones choque-gota a números de Mach de hasta 1.47 en un contexto de fluido no ideal, demostrando la capacidad del modelo para capturar estructuras de ondas complejas.

Resultados
El modelo fue validado a través de una jerarquía de configuraciones cada vez más complejas:

1. Estabilidad Lineal y Dispersión: Las simulaciones de la velocidad del sonido, la disipación de ondas de corte y la conductividad térmica mostraron un excelente acuerdo con soluciones analíticas para un fluido de van der Waals ajustado a las propiedades del nitrógeno en un rango de números de Mach y temperaturas reducidas.
2. Consistencia Multifásica: El modelo reprodujo con precisión las densidades de coexistencia líquido-vapor y los perfiles de interfaz, demostrando convergencia de segundo orden con el refinamiento de la cuadrícula y acuerdo con la regla de áreas iguales de Maxwell.
3. Tubos de Choque No Clásicos: Se simularon tres casos de tubos de choque 1-D para sondear la dinámica de ondas no clásicas caracterizada por el derivado fundamental de la dinámica de gases ($\Gamma$).
   • Los casos con $\Gamma < 0$ capturaron correctamente la formación de choques de rarefacción (ondas de expansión afiladas) y abanicos de compresión (ondas de compresión dispersas), fenómenos distintos de la dinámica de gases clásica.
   • El modelo transitó exitosamente entre regímenes clásicos ($\Gamma > 0$) y no clásicos ($\Gamma < 0$) dentro de una sola simulación.
4. Interacción Choque-Coluna Líquida: Se realizó una simulación 2D de una onda de choque plana interactuando con una columna de líquido saturado cerca del punto crítico ($M_a = 1.47$).
   • El modelo capturó estructuras de ondas complejas que incluyen choques transmitidos, choques reflejados, tallos de Mach y puntos triples.
   • La comparación cuantitativa de la evolución del ancho de la columna líquida frente a datos experimentales y resultados numéricos previos mostró un buen acuerdo.
   • El análisis termodinámico reveló que la interacción del choque impulsa el vapor circundante a un estado supercrítico mientras la columna líquida permanece marginalmente subcrítica, capturando con precisión el fuerte acoplamiento termodinámico.

Significado y Afirmaciones
Los autores afirman que este trabajo extiende la aplicabilidad de los métodos de Boltzmann en Red a flujos compresibles no ideales de alta velocidad con una plantilla cinética mínima. El significado principal radica en proporcionar un marco termodinámicamente consistente y numéricamente estable que no depende de plantillas extendidas o regularizaciones ad hoc. El modelo se presenta como una teoría cinética reducida que apunta al límite hidrodinámico NSK, capaz de manejar el rango completo de comportamientos no ideales, incluida la región espínodal, sin las inestabilidades numéricas a menudo asociadas con la CFD convencional en estos regímenes. La simulación exitosa de interacciones choque-gota a altos números de Mach sirve como una validación rigurosa, demostrando el potencial del método para investigar flujos complejos en tecnologías supercríticas y procesos impulsados por transiciones de fase.