Modelo de Boltzmann en red para dinámica de fluidos compresibles no ideales

Información Básica

• ID del artículo: 2510.14712
• Título: Lattice Boltzmann model for non-ideal compressible fluid dynamics
• Autores: S. A. Hosseini, M. Feinberg, I. V. Karlin (ETH Zurich)
• Clasificación: physics.flu-dyn
• Revista de publicación: Journal of Fluid Mechanics
• Enlace del artículo: https://arxiv.org/abs/2510.14712

Resumen

En este artículo se propone un método de Boltzmann en red para simular flujos de fluidos no ideales compresibles. El método utiliza una red de primeros vecinos, introduciendo un conjunto coherente de términos correctivos a través de un atractor cuasiequilibrio, garantizando una tasa de disipación de Navier-Stokes definida positiva e invariante de Galileo. Esta construcción evita la necesidad de plantillas extendidas o regularización temporal, manteniendo simultáneamente estabilidad numérica y consistencia termodinámica en un amplio rango de estados de flujo. El modelo resultante reproduce con precisión la dinámica de fluidos en los niveles de Euler y Navier-Stokes. Como validación rigurosa, los autores demuestran por primera vez dentro del marco de Boltzmann en red la simulación cuantitativamente exacta de interacciones gota-onda de choque con números de Mach de hasta 1.47.

Antecedentes de investigación y motivación

Definición del problema

La dinámica de fluidos compresibles no ideales es una rama emergente y de rápido desarrollo de la mecánica de fluidos, impulsada principalmente por la aparición de métodos y tecnologías que operan en estados cercanos al crítico, transcríticos y supercríticos. Esta clase de fluidos incluye:

1. Fluidos supercríticos: Fluidos bajo condiciones de presión y temperatura por encima del punto crítico
2. Vapores densos: Gases cuyo factor de compresibilidad se desvía del comportamiento de gas ideal
3. Fluidos BZT: Fluidos de complejidad molecular elevada con derivadas negativas de dinámica de gases fundamentales

Importancia

1. Demanda de aplicaciones industriales: Tecnologías energéticas como ciclos Rankine orgánicos y turbinas de CO₂ supercrítico
2. Escasez de datos experimentales: Configuraciones experimentales complejas y dificultad en la obtención de datos
3. Necesidad de herramientas numéricas: Requerimiento de herramientas numéricas consistentes y eficientes para comprender fenómenos físicos complejos

Limitaciones de métodos existentes

1. Métodos CFD tradicionales: En regiones de espinodal termodinámicamente inestables, el cuadrado de la velocidad del sonido adiabático se vuelve negativo, causando que las ecuaciones de evolución cambien de hiperbólicas a parabólicas
2. Modelos LBM existentes: Principalmente orientados a aplicaciones de ebullición, careciendo de modelos sistemáticos para flujos compresibles no ideales
3. Problemas de estabilidad numérica: Desafíos de estabilidad bajo números de Mach elevados y efectos no ideales fuertes

Contribuciones principales

1. Propuesta de nuevo modelo LBM de doble función de distribución: Implementación de simulación de flujos no ideales compresibles usando la red clásica de primeros vecinos (D3Q27)
2. Introducción del mecanismo de atractor cuasiequilibrio: Garantía de tasa de disipación definida positiva e invariante de Galileo mediante términos correctivos fᵢ y gᵢ
3. Implementación de control independiente de viscosidad volumétrica: Resolución del problema de viscosidad volumétrica potencialmente negativa en la estructura BGK
4. Primera simulación de interacción gota-onda de choque a alto número de Mach: Simulación cuantitativamente exacta con Ma=1.47 dentro del marco LBM
5. Mantenimiento de consistencia termodinámica: Estabilidad numérica sostenida en un amplio rango de estados de flujo

Detalles del método

Definición de la tarea

Simulación de fluidos compresibles no ideales que satisfacen el conjunto de ecuaciones de Navier-Stokes-Korteweg (NSK):

• Conservación de masa: ∂ₜρ + ∇·(ρu) = 0
• Conservación de momento: ∂ₜ(ρu) + ∇·(ρu⊗u) + ∇P + ∇·T^NS + ∇·T^K = 0
• Conservación de energía: ∂ₜ(ρE) + ∇·((ρE+P)u) + u·(∇·T^K) + ∇·(u·T^NS) + ∇·q = 0

Donde T^K es el tensor de tensión superficial de Korteweg y T^NS es el tensor de esfuerzo viscoso.

Arquitectura del modelo

1. Diseño de doble función de distribución

Utilización de dos conjuntos de funciones de distribución {fᵢ, gᵢ}, i=1,...,27:

• fᵢ: Manejo de conservación de masa y momento
• gᵢ: Manejo de conservación de energía

2. Ecuación de Boltzmann en red

fᵢ(x+cᵢδt, t+δt) = fᵢ + 2β(fᵢᵉᵠ - fᵢ) + (1-β)(fᵢ* - fᵢᵉᵠ)
gᵢ(x+cᵢδt, t+δt) = gᵢ + 2β(gᵢᵉᵠ - gᵢ) + (1-β)(gᵢ* - gᵢᵉᵠ)

Donde β = Pδt/(2μ + Pδt) es el parámetro de relajación.

3. Funciones de distribución de equilibrio

Adopción de forma de producto:

fᵢᵉᵠ = ρ ∏_α Ψᵢα(uα, P/ρ + uα²)

Donde Ψᵢα es la función base:

Ψᵢα(ξα, ζαα) = 1 - cᵢα² + ½[(3cᵢα² - 2)ζαα + cᵢαξα]

4. Términos correctivos de cuasiequilibrio

Corrección de función de distribución f:

• Término de fuerza: Fα = -κρ∂α∇²ρ (fuerza de Korteweg)
• Corrección de esfuerzo: Φαα contiene corrección de tensor de esfuerzo diagonal y control independiente de viscosidad volumétrica

Corrección de función de distribución g:

• Velocidad desplazada: u* = u + δt F/ρ
• Temperatura desplazada: T* = T - δt²F·F/(2ρ²cᵥ)
• Corrección de flujo de calor: Consideración de efectos de gradiente de entalpía en gases no ideales

Puntos de innovación técnica

1. Mecanismo de atractor cuasiequilibrio

Mediante la introducción de fᵢ y gᵢ, el sistema no se relaja simplemente hacia el equilibrio local, sino hacia un estado cuasiequilibrio modificado, lo que garantiza:

• Dinámica macroscópica correcta
• Estabilidad numérica
• Invariancia de Galileo

2. Control independiente de viscosidad volumétrica

En el modelo BGK tradicional, la viscosidad volumétrica η' = ((D+2)/D - ρcs²/P)μ puede ser negativa. Mediante el término correctivo Φ' se implementa control independiente:

Φ' = P((D+2)/D - ρcs²/P - η/μ)(∇·u)

3. Consistencia termodinámica

Mediante diseño cuidadoso de temperatura desplazada T* y correcciones de energía, se garantiza consistencia termodinámica bajo ecuaciones de estado no lineales.

Configuración experimental

Casos de validación

1. Validación básica:
   • Medición de velocidad del sonido
   • Verificación de viscosidad de corte
   • Flujo de Couette térmico
   • Disipación de modo normal
2. Validación de flujos multifásicos:
   • Densidades de coexistencia líquido-vapor
   • Consistencia de interfaz y convergencia
3. Validación de flujos compresibles:
   • Tubo de choque no ideal unidimensional
   • Interacción de onda de choque-columna de líquido

Parámetros de sustancia

Validación principal utilizando nitrógeno (N₂):

• Presión crítica: Pc = 3.4×10⁶ Pa
• Densidad crítica: ρc = 241.96 kg/m³
• Temperatura crítica: Tc = 126.2 K
• R/cv = 0.4

Ecuación de estado

Adopción de ecuación de estado de van der Waals:

P(ρ,T) = ρRT/(1-bρ) - aρ²

Resultados experimentales

Resultados principales

1. Verificación de cantidades físicas básicas

• Velocidad del sonido: Coincidencia perfecta con soluciones analíticas en ramas de líquido y vapor saturados en rango Tr∈0.7,1
• Viscosidad de corte: Mantenimiento de invariancia de Galileo bajo diferentes números de Mach, error <1%
• Flujo de Couette térmico: Consistencia con soluciones analíticas en Pr∈{0.6,1.2,4.9} y Ma∈{0.8,1.2,1.6}

2. Propiedades de flujos multifásicos

• Densidades de coexistencia: Coincidencia completa con construcción de igual área de Maxwell en rango Tr∈0.3,0.99
• Convergencia de interfaz: Exhibición de buenas características de convergencia bajo diferentes resoluciones de malla

3. Tubo de choque no ideal

Simulación exitosa de tres configuraciones diferentes de tubo de choque no ideal:

• Caso I: Onda de rarefacción en región Γ<0
• Caso II: Comportamiento de onda de compresión clásica
• Caso III: Onda de rarefacción mixta cruzando dinámicamente la línea Γ=0

Todos los casos muestran excelente concordancia con datos de referencia.

4. Interacción de onda de choque-gota

• Primera implementación en marco LBM de simulación a alto número de Mach Ma=1.47
• Captura precisa de estructuras de onda complejas: onda incidente, onda transmitida, onda reflejada, vara de Mach
• Evolución de deformación de gota altamente consistente con datos experimentales

Experimentos de ablación

Verificación de convergencia del modelo mediante pruebas de interfaz con diferentes resoluciones de malla:

• Cuatro niveles de resolución con δx desde 5μm hasta 0.1μm
• Exhibición de buenas características de convergencia de malla

Hallazgos clave

1. Estabilidad numérica: Mantenimiento de estabilidad bajo números de Mach elevados y efectos no ideales fuertes
2. Eficiencia computacional: Uso de red estándar de primeros vecinos con costo computacional controlable
3. Fidelidad física: Reproducción precisa de fenómenos físicos complejos en fluidos no ideales

Trabajo relacionado

Modelos LBM multifásicos tradicionales

1. Modelo de gradiente de color (Gunstensen et al. 1991)
2. Modelo de pseudopotencial (Shan & Chen 1993)
3. Modelo de energía libre (Swift et al. 1996)

Desarrollo de LBM compresible

• He & Doolen (2002): Primer intento de flujos no ideales compresibles
• Modelos existentes principalmente orientados a aplicaciones de ebullición
• Carencia de modelos sistemáticos para flujos compresibles no ideales

Ventajas relativas de este trabajo

1. Enfoque sistemático: Marco completo de doble función de distribución
2. Fundamento teórico: Análisis multiescala riguroso
3. Completitud de validación: Validación integral desde cantidades físicas básicas hasta fenómenos de flujo complejos

Conclusiones y discusión

Conclusiones principales

1. Desarrollo exitoso de modelo LBM aplicable a fluidos compresibles no ideales
2. Implementación de estabilidad numérica y consistencia física mediante mecanismo de atractor cuasiequilibrio
3. Primera implementación en marco LBM de simulación de interacción gota-onda de choque a alto número de Mach
4. Mantenimiento de consistencia termodinámica del modelo en amplio rango de estados de flujo

Limitaciones

1. Modelo de tiempo de relajación único: Uso actual de modelo de colisión BGK simple, posiblemente insuficiente para simulación de turbulencia
2. Restricciones de ecuación de estado: Validación principal de ecuación de van der Waals, requiriéndose verificación adicional de otras ecuaciones complejas
3. Validación tridimensional limitada: Concentración de validación principalmente en configuraciones unidimensionales y bidimensionales

Direcciones futuras

1. Extensión a modelos de colisión más robustos: Para simulación de configuraciones turbulentas
2. Ecuaciones de estado más complejas: Como Peng-Robinson, entre otras
3. Aplicaciones tridimensionales a gran escala: Aplicación a problemas de ingeniería práctica

Evaluación profunda

Fortalezas

1. Rigor teórico: Análisis multiescala completo y derivaciones matemáticas
2. Innovación metodológica: El mecanismo de atractor cuasiequilibrio constituye una innovación importante
3. Completitud de validación: Validación sistemática desde física básica hasta aplicaciones complejas
4. Valor práctico: Resolución de problemas importantes de aplicación ingenieril
5. Eficiencia numérica: Uso de red estándar evitando complejidad computacional

Insuficiencias

1. Modelo de colisión simple: Aplicabilidad limitada del modelo BGK a flujos complejos
2. Dependencia de parámetros: Selección de ciertos parámetros correctivos carece de orientación sistemática
3. Validación tridimensional insuficiente: Necesidad de mayor validación en geometrías complejas tridimensionales

Impacto

1. Contribución académica: Apertura de nuevas direcciones para LBM en campo de flujos compresibles no ideales
2. Valor ingenieril: Provisión de herramienta numérica importante para tecnología de fluidos supercríticos
3. Significado metodológico: El concepto de atractor cuasiequilibrio puede influir en otros métodos numéricos

Escenarios de aplicación

1. Tecnología de fluidos supercríticos: Turbinas de CO₂, ciclos Rankine orgánicos
2. Flujos multifásicos: Ruptura de gotas, procesos de cambio de fase
3. Dinámica de ondas de choque: Interacciones complejas de ondas en flujos de alta velocidad
4. Simulación numérica directa: Investigación física fundamental requiriendo alta precisión

Referencias

1. Guardone, A. et al. (2024). Nonideal compressible fluid dynamics of dense vapors and supercritical fluids. Annual Review of Fluid Mechanics, 56, 241-269.
2. Hosseini, S.A. & Karlin, I.V. (2023). Lattice Boltzmann for non-ideal fluids: Fundamentals and practice. Physics Reports, 1030, 1-137.
3. Shan, X. & Chen, H. (1993). Lattice Boltzmann model for simulating flows with multiple phases and components. Physical Review E, 47(3), 1815.