🔢 mathematics

Non-equilibrium modelling of polyatomic gases: generic framework and 11-moment model

Este artículo presenta un marco termodinámicamente consistente para modelar gases poliatómicos integrando el formalismo GENERIC con una descripción de 11 momentos para unificar la termodinámica racional y la termodinámica irreversible extendida, asegurando el cumplimiento de la segunda ley y validando el modelo frente a simulaciones DSMC y experimentos para estructuras de ondas de choque en nitrógeno.

Autores originales: Masrakain Ahmad, Anil Kumar, Anirudh Singh Rana

Publicado 2026-08-26

📖 6 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Masrakain Ahmad, Anil Kumar, Anirudh Singh Rana

Artículo original bajo licencia CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

Los gases no son solo nubes invisibles de partículas uniformes; son colecciones complejas de moléculas que pueden girar, vibrar y almacenar energía de formas que los átomos simples no pueden. Cuando un gas está tranquilo y estable, estos movimientos internos se asientan en un equilibrio predecible con el movimiento hacia adelante de las moléculas. Sin embargo, cuando un gas es empujado violentamente —como cuando una nave espacial reentra en la atmósfera o cuando una onda de choque recorre el aire— este equilibrio se rompe. La energía que normalmente fluye suavemente entre la velocidad de avance de las moléculas y su giro interno se desincroniza. Una parte del gas puede calentarse rápidamente mientras otra parte se queda rezagada, creando un estado caótico de no equilibrio. Comprender exactamente cómo sucede esto es crítico para los ingenieros que diseñan vehículos de alta velocidad y para los científicos que intentan predecir cómo se mueve el calor en entornos extremos. Durante décadas, las herramientas estándar utilizadas para modelar estos gases han sido lo suficientemente precisas para flujos suaves, pero fallan cuando el gas es llevado a sus límites, perdiendo a menudo los sutiles retrasos e intercambios de energía que definen estos momentos violentos.

Un equipo de investigadores del Instituto de Tecnología y Ciencia de Birla en la India ha desarrollado una nueva forma de describir estos momentos caóticos, creando un marco matemático que respeta las leyes fundamentales de la termodinámica mientras captura la realidad desordenada de los gases poliatómicos. Su trabajo se centra en un tipo específico de molécula de gas que tiene partes internas capaces de rotar y vibrar, como el nitrógeno, que constituye la mayor parte de nuestra atmósfera. Los investigadores construyeron un modelo que rastrea once propiedades diferentes del gas simultáneamente, en lugar de solo las pocas propiedades básicas como la presión y la temperatura utilizadas en métodos más antiguos. Este nuevo enfoque trata al gas como un sistema donde los movimientos reversibles, como el flujo suave de un río, y los movimientos irreversibles, como la fricción que genera calor, están claramente separados pero matemáticamente vinculados. Al hacer esto, se aseguraron de que su modelo nunca viole la segunda ley de la termodinámica, que dicta que el desorden, o la entropía, debe siempre aumentar en un sistema cerrado. Este es un logro significativo porque muchos intentos previos de mejorar los modelos de gases o bien rompían esta ley fundamental o se volvían tan complicados que eran imposibles de usar de manera fiable.

El núcleo de este nuevo marco es una estructura que actúa como un conjunto de reglas para cómo el gas evoluciona a lo largo del tiempo. Utiliza dos herramientas matemáticas distintas para gobernar el comportamiento del gas: una herramienta maneja los movimientos reversibles que conservan la energía, mientras que la otra maneja los procesos irreversibles que disipan la energía y crean calor. Esta separación asegura que el modelo sea físicamente realista en cada paso. Los investigadores incluyeron variables específicas para rastrear la temperatura del movimiento hacia adelante de las moléculas por separado de la temperatura de su giro y vibración interna. También introdujeron un "tensor de temperatura", que es una forma de describir cómo la temperatura podría ser diferente en distintas direcciones dentro del gas, un fenómeno que ocurre cuando el gas está bajo un estrés extremo. Al entrelazar estos elementos, el equipo creó un conjunto de ecuaciones que produce naturalmente la cantidad correcta de entropía, garantizando que el modelo se comporte como lo haría un gas real.

Para probar si su teoría realmente funcionaba en el mundo real, los investigadores aplicaron su modelo a un problema clásico de la física: la estructura de una onda de choque en el nitrógeno gaseoso. Una onda de choque es un límite agudo donde las propiedades del gas cambian casi instantáneamente, creando un gradiente pronunciado que es difícil de modelar. El equipo simuló cómo se vería esta onda en nitrógeno gaseoso moviéndose a velocidades que oscilan entre 1.7 y 3.8 veces la velocidad del sonido. Compararon sus resultados con otras dos fuentes de verdad: simulaciones computacionales directas que rastrean moléculas de gas individuales, conocidas como Método de Monte Carlo de Simulación Directa, y mediciones experimentales reales tomadas en túneles de viento. Los resultados fueron sorprendentes. El nuevo modelo coincidió con los datos experimentales y las detalladas simulaciones computacionales con alta precisión en todo el rango de velocidades probadas. Predijo con éxito cómo cambiaba la densidad del gas y cómo evolucionaban las temperaturas traslacional y de rotación interna a medida que pasaba la onda de choque.

En contraste, los modelos más antiguos que tratan al gas como si tuviera una sola temperatura o que utilizan aproximaciones más simples mostraron errores notables, particularmente en la región donde el gas está lejos del equilibrio. El nuevo modelo capturó el sutil retraso entre el calentamiento del movimiento hacia adelante de las moléculas y el calentamiento de su rotación interna, un retraso que es crucial para predicciones precisas. Los investigadores encontraron que su enfoque se mantenía preciso hasta Mach 3.8, una velocidad donde muchos otros modelos comienzan a tener dificultades. Aunque el modelo mostró ligeras desviaciones a velocidades aún mayores, probablemente debido a que la suposición de que las moléculas de gas tienen propiedades internas fijas se vuelve menos válida bajo calor extremo, el rendimiento a velocidades moderadas fue robusto. Esto sugiere que el marco proporciona una forma fiable y matemáticamente sólida para predecir cómo se comportan los gases complejos cuando son empujados al límite de la estabilidad.

La importancia de este trabajo no reside solo en los números que produce, sino en la estructura que proporciona. Al fundamentar el modelo en un marco termodinámico riguroso, los investigadores han creado una herramienta que es tanto físicamente consistente como matemáticamente estable. Esto significa que los ingenieros y científicos pueden usarlo para diseñar mejores sistemas de protección térmica para naves espaciales o para entender la transferencia de calor en dispositivos a microescala sin preocuparse de que la matemática subyacente produzca resultados imposibles. El estudio demuestra que es posible capturar el rico y complejo comportamiento de los gases poliatómicos sin sacrificar las leyes fundamentales que gobiernan el universo. Mientras los investigadores miran hacia el futuro, planean refinar el modelo aún más para manejar condiciones todavía más extremas donde las propiedades internas de las propias moléculas podrían cambiar, pero por ahora, han establecido una base sólida para comprender la danza caótica de las moléculas de gas en los entornos más violentos.

Resumen Técnico: Modelado de No Equilibrio de Gases Poliatómicos: Marco Genérico y Modelo de 11 Momentos

Planteamiento del Problema
El modelado preciso de gases poliatómicos bajo condiciones de no equilibrio es crítico para aplicaciones que van desde la reentrada aeroespacial hasta la transferencia de calor a microescala. A diferencia de los gases monoatómicos, los gases poliatómicos poseen grados de libertad internos (rotacionales y vibracionales) que intercambian energía con el movimiento traslacional, creando dinámicas de relajación complejas. Las ecuaciones clásicas de Navier–Stokes–Fourier (NSF), limitadas a una precisión de primer orden en el número de Knudsen, no logran capturar efectos de no equilibrio fuertes. Aunque los modelos extendidos (por ejemplo, el de Gauss Regularizado de 11 momentos o RG11) ofrecen una mayor precisión, a menudo carecen de una ecuación de entropía de forma cerrada, lo que dificulta la verificación rigurosa de la consistencia termodinámica. Además, los cierres de momentos existentes a veces tienen dificultades para mantener la hiperbolicidad y la realizabilidad en regímenes no lineales. El desafío radica en desarrollar un modelo que capture los intercambios de energía interna, asegure la buena definición matemática (well-posedness) y se adhiera estrictamente a las leyes de la termodinámica, particularmente a la segunda ley.

Metodología
Los autores proponen un marco termodinámicamente consistente integrando la Ecuación General para el Acoplamiento Reversible-Irreversible de No Equilibrio (GENERIC) con una descripción de 11 momentos. La metodología procede de la siguiente manera:

Fundamento de la Teoría Cinética: Partiendo de la ecuación de Boltzmann para gases poliatómicos, los autores definen la función de distribución $f$ considerando la velocidad traslacional y la energía interna. Introducen campos distintos para la temperatura traslacional ( $\theta_{tr}$ ), la temperatura interna ( $\theta_{in}$ ) y un tensor de temperatura ( $\Theta_{ij}$ ).
Cierre de Momentos: Se deriva un sistema de 11 momentos maximizando la entropía bajo restricciones de masa, momento, energía total y el tensor de temperatura. Esto produce una función de distribución específica ( $f_{11}$ ) que combina una Gaussiana para los modos traslacionales y una distribución Gamma para los modos internos.
Formulación GENERIC: Las ecuaciones de evolución se estructuran dentro del marco GENERIC:
- Variables de Estado: Densidad de masa ( $\rho$ ), densidad de momento ( $\rho v_i$ ), tensor de temperatura ( $\Theta_{ij}$ ) y temperatura interna ( $\theta_{in}$ ).
- Dinámica Reversible: Codificada mediante una matriz de Poisson ( $L$ ) que satisface la identidad de Jacobi y las condiciones de degeneración, asegurando la conservación de la energía y la entropía en ausencia de disipación.
- Dinámica Irreversible: Codificada mediante una matriz de fricción ( $M$ ), que es simétrica y semidefinida positiva. Esta matriz acopla los flujos termodinámicos (flujos de calor, tasas de esfuerzo) con las fuerzas termodinámicas (gradientes de temperaturas inversas y de las inversas de los tensores).
Relaciones Constitutivas: Las relaciones de cierre para los flujos de calor y los términos de producción se derivan explícitamente para asegurar una producción de entropía no negativa. El modelo incorpora el acoplamiento entre los grados de libertad traslacionales e internos a través de coeficientes fenomenológicos ( $\zeta_{11}, \zeta_{12}, \zeta_{22}$ ) y un parámetro de acoplamiento ( $\bar{\omega}$ ), derivados de un marco hidrodinámico para moléculas esféricas.
Validación Numérica: El modelo se aplica a la estructura de una onda de choque estacionaria unidimensional en gas nitrógeno. Las ecuaciones gobernantes se reducen a un sistema de ecuaciones diferenciales ordinarias (ODEs) y se resuelven utilizando un método de punto medio de diferencias finitas. Los resultados se comparan con simulaciones de Monte Carlo de Método Directo de Simulación (DSMC) y datos experimentales.

Contribuciones Clave

Marco Unificado: El artículo unifica la Termodinámica de Procesos Irreversibles Racional (RIT) y la Termodinámica Irreversible Extendida (EIT) en una estructura única compatible con GENERIC. Esto asegura que la segunda ley de la termodinámica se satisfaga por construcción mediante la definición explícita de las matrices de Poisson y de fricción.
Consistencia Termodinámica: A diferencia de otros modelos de momentos previos (por ejemplo, RG11), esta formulación proporciona una ecuación de entropía de forma cerrada y una expresión explícita de la producción de entropía, garantizando una generación de entropía no negativa.
Modelo de 11 Momentos Poliatómico: Los autores presentan la primera demostración de un modelo de gas poliatómico que combina GENERIC y cierres de momentos inspirados en RET. El modelo rastrea explícitamente los intercambios de energía traslacional e interna, así como el tensor de temperatura, ofreciendo una representación físicamente significativa del comportamiento de no equilibrio.
Estructura Matemática: Se demuestra que el sistema es hiperbólico-parabólico (específicamente, hiperbólico-incompletamente parabólico), con autovalores reales para el subsistema de fluido sin viscosidad (inviscid), lo que asegura la buena definición local en el espacio de estados realizables.

Resultados
El modelo fue probado en estructuras de ondas de choque en nitrógeno a números de Mach de hasta 3.8.

Perfiles de Densidad: El modelo de 11 momentos demuestra una fuerte concordancia tanto con los datos experimentales como con las simulaciones DSMC para los perfiles de densidad normalizados en todo el rango de Mach probado.
Perfiles de Temperatura: El modelo captura con precisión la separación entre las temperaturas traslacional ( $\theta_{tr}$ ) e interna ( $\theta_{in}$ ). Muestra una concordancia significativamente mejor con los datos de DSMC en comparación con el modelo de dos temperaturas (2T), particularmente en la región de no equilibrio del choque.
Limitaciones: A números de Mach más altos (superiores a 4), aparecen ligeras desviaciones en los perfiles de temperatura traslacional aguas abajo. Los autores atribuyen esto al supuesto de un parámetro de relajación $\delta$ constante (que implica capacidades caloríficas independientes de la temperatura), lo cual resulta restrictivo bajo gradientes de temperatura fuertes donde la redistribución de energía depende de la temperatura.

Significancia y Reivindicaciones
Los autores afirman modestamente que la contribución principal de este trabajo no es necesariamente una precisión predictiva superior sobre todos los modelos existentes, sino la demostración de una estructura transparente y termodinámicamente consistente para la dinámica de gases poliatómicos.

El marco GENERIC identifica explícitamente las contribuciones reversibles e irreversibles, el flujo de entropía y el acoplamiento de la energía traslacional e interna.
El modelo satisface automáticamente las leyes de conservación y la segunda ley de la termodinámica debido a su estructura integrada.
El acuerdo con DSMC y los datos experimentales valida el cierre de 11 momentos específico desarrollado dentro de este marco, más que la estructura GENERIC en sí misma.
El trabajo establece una base para futuras extensiones, como la incorporación de parámetros de relajación dependientes de la temperatura, manteniendo al mismo tiempo la rigurosa estructura termodinámica.

Más como este