🔢 mathematics

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

This paper presents a thermodynamically consistent framework for modelling polyatomic gases by integrating the GENERIC formalism with an 11-moment description to unify rational and extended irreversible thermodynamics, ensuring second-law compliance and validating the model against DSMC simulations and experiments for shock wave structures in nitrogen.

Original authors: Masrakain Ahmad, Anil Kumar, Anirudh Singh Rana

Published 2026-08-26

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Masrakain Ahmad, Anil Kumar, Anirudh Singh Rana

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Gases are not just invisible clouds of uniform particles; they are complex collections of molecules that can spin, vibrate, and store energy in ways that simple atoms cannot. When a gas is calm and stable, these internal movements settle into a predictable balance with the forward motion of the molecules. However, when a gas is pushed violently—such as when a spacecraft re-enters the atmosphere or when a shock wave ripples through the air—this balance breaks. The energy that usually flows smoothly between the molecules' forward speed and their internal spinning gets out of sync. One part of the gas might heat up rapidly while another part lags behind, creating a chaotic, non-equilibrium state. Understanding exactly how this happens is critical for engineers designing high-speed vehicles and for scientists trying to predict how heat moves in extreme environments. For decades, the standard tools used to model these gases have been accurate enough for gentle flows but fail when the gas is pushed to its limits, often missing the subtle delays and energy exchanges that define these violent moments.

A team of researchers at the Birla Institute of Technology and Science in India has developed a new way to describe these chaotic moments, creating a mathematical framework that respects the fundamental laws of thermodynamics while capturing the messy reality of polyatomic gases. Their work focuses on a specific type of gas molecule that has internal parts capable of rotating and vibrating, like nitrogen, which makes up most of our atmosphere. The researchers built a model that tracks eleven different properties of the gas simultaneously, rather than just the few basic ones like pressure and temperature used in older methods. This new approach treats the gas as a system where reversible movements, like the smooth flow of a river, and irreversible movements, like the friction that generates heat, are clearly separated but mathematically linked. By doing this, they ensured that their model never violates the second law of thermodynamics, which dictates that disorder, or entropy, must always increase in a closed system. This is a significant achievement because many previous attempts to improve gas models either broke this fundamental law or became so complicated that they were impossible to use reliably.

The core of this new framework is a structure that acts like a set of rules for how the gas evolves over time. It uses two distinct mathematical tools to govern the gas's behavior: one tool handles the reversible, energy-conserving motions, while the other handles the irreversible, energy-dissipating processes that create heat. This separation ensures that the model remains physically realistic at every step. The researchers included specific variables to track the temperature of the molecules' forward motion separately from the temperature of their internal spinning and vibration. They also introduced a "temperature tensor," which is a way of describing how the temperature might be different in different directions within the gas, a phenomenon that occurs when the gas is under extreme stress. By weaving these elements together, the team created a set of equations that naturally produces the correct amount of entropy, guaranteeing that the model behaves like a real gas would.

To test whether their theory actually worked in the real world, the researchers applied their model to a classic problem in physics: the structure of a shock wave in nitrogen gas. A shock wave is a sharp boundary where gas properties change almost instantly, creating a steep gradient that is difficult to model. The team simulated how this wave would look in nitrogen gas moving at speeds ranging from 1.7 to 3.8 times the speed of sound. They compared their results against two other sources of truth: direct computer simulations that track individual gas molecules, known as Direct Simulation Monte Carlo, and actual experimental measurements taken in wind tunnels. The results were striking. The new model matched the experimental data and the detailed computer simulations with high precision across the entire range of speeds tested. It successfully predicted how the density of the gas changed and how the translational and internal temperatures evolved as the shock wave passed.

In contrast, older models that treat the gas as having a single temperature or that use simpler approximations showed noticeable errors, particularly in the region where the gas is far from equilibrium. The new model captured the subtle lag between the heating of the molecules' forward motion and the heating of their internal rotation, a delay that is crucial for accurate predictions. The researchers found that their approach remained accurate up to Mach 3.8, a speed where many other models begin to struggle. While the model did show slight deviations at even higher speeds, likely because the assumption that the gas molecules have fixed internal properties becomes less valid under extreme heat, the performance at moderate speeds was robust. This suggests that the framework provides a reliable and mathematically sound way to predict how complex gases behave when they are pushed to the edge of stability.

The significance of this work lies not just in the numbers it produces, but in the structure it provides. By grounding the model in a rigorous thermodynamic framework, the researchers have created a tool that is both physically consistent and mathematically stable. This means that engineers and scientists can use it to design better thermal protection systems for spacecraft or to understand heat transfer in micro-scale devices without worrying that the underlying math will produce impossible results. The study demonstrates that it is possible to capture the rich, complex behavior of polyatomic gases without sacrificing the fundamental laws that govern the universe. As the researchers look toward the future, they plan to refine the model further to handle even more extreme conditions where the internal properties of the molecules themselves might change, but for now, they have established a solid foundation for understanding the chaotic dance of gas molecules in the most violent environments.

Technical Summary: Non-Equilibrium Modelling of Polyatomic Gases: Generic Framework and 11-Moment Model

Problem Statement
Accurate modeling of polyatomic gases under non-equilibrium conditions is critical for applications ranging from aerospace re-entry to microscale heat transfer. Unlike monatomic gases, polyatomic gases possess internal degrees of freedom (rotational and vibrational) that exchange energy with translational motion, creating complex relaxation dynamics. Classical Navier–Stokes–Fourier (NSF) equations, limited to first-order accuracy in the Knudsen number, fail to capture strong non-equilibrium effects. While extended models (e.g., Regularized Gaussian 11-moment or RG11) offer higher accuracy, they often lack a closed-form entropy equation, making rigorous thermodynamic consistency difficult to verify. Furthermore, existing moment closures sometimes struggle to maintain hyperbolicity and realizability in nonlinear regimes. The challenge lies in developing a model that captures internal energy exchanges, ensures mathematical well-posedness, and strictly adheres to the laws of thermodynamics, particularly the second law.

Methodology
The authors propose a thermodynamically consistent framework by integrating the General Equation for Non-Equilibrium Reversible-Irreversible Coupling (GENERIC) with an 11-moment description. The methodology proceeds as follows:

Kinetic Theory Foundation: Starting from the Boltzmann equation for polyatomic gases, the authors define the distribution function $f$ considering translational velocity and internal energy. They introduce distinct fields for translational temperature ( $\theta_{tr}$ ), internal temperature ( $\theta_{in}$ ), and a temperature tensor ( $\Theta_{ij}$ ).
Moment Closure: An 11-moment system is derived by maximizing entropy under constraints for mass, momentum, total energy, and the temperature tensor. This yields a specific distribution function ( $f_{11}$ ) combining a Gaussian for translational modes and a Gamma distribution for internal modes.
GENERIC Formulation: The evolution equations are structured within the GENERIC framework:
- State Variables: Mass density ( $\rho$ ), momentum density ( $\rho v_i$ ), temperature tensor ( $\Theta_{ij}$ ), and internal temperature ( $\theta_{in}$ ).
- Reversible Dynamics: Encoded via a Poisson matrix ( $L$ ) that satisfies the Jacobi identity and degeneracy conditions, ensuring energy and entropy conservation in the absence of dissipation.
- Irreversible Dynamics: Encoded via a friction matrix ( $M$ ) that is symmetric and positive semi-definite. This matrix couples thermodynamic fluxes (heat fluxes, stress rates) to thermodynamic forces (gradients of inverse temperatures and tensor inverses).
Constitutive Relations: The closure relations for heat fluxes and production terms are derived explicitly to ensure non-negative entropy production. The model incorporates coupling between translational and internal degrees of freedom through phenomenological coefficients ( $\zeta_{11}, \zeta_{12}, \zeta_{22}$ ) and a coupling parameter ( $\bar{\omega}$ ), derived from a hydrodynamic framework for spherical molecules.
Numerical Validation: The model is applied to a one-dimensional steady shock wave structure in nitrogen gas. The governing equations are reduced to a system of ordinary differential equations (ODEs) and solved using a finite difference midpoint method. Results are compared against Direct Simulation Monte Carlo (DSMC) simulations and experimental data.

Key Contributions

Unified Framework: The paper unifies Rational Irreversible Thermodynamics (RIT) and Extended Irreversible Thermodynamics (EIT) into a single GENERIC-compatible structure. This ensures that the second law of thermodynamics is satisfied by construction through the explicit definition of Poisson and friction matrices.
Thermodynamic Consistency: Unlike previous moment models (e.g., RG11), this formulation provides a closed-form entropy equation and an explicit expression for entropy production, guaranteeing non-negative entropy generation.
11-Moment Polyatomic Model: The authors present the first demonstration of a polyatomic gas model that combines GENERIC and RET-inspired moment closures. The model explicitly tracks translational and internal energy exchanges, as well as the temperature tensor, offering a physically meaningful representation of non-equilibrium behavior.
Mathematical Structure: The system is shown to be hyperbolic-parabolic (specifically, hyperbolic-incompletely parabolic), with real eigenvalues for the inviscid subsystem, ensuring local well-posedness in the realizable state space.

Results
The model was tested on shock wave structures in nitrogen gas across Mach numbers up to 3.8.

Density Profiles: The 11-moment model demonstrates strong agreement with both experimental data and DSMC simulations for normalized density profiles across the tested Mach range.
Temperature Profiles: The model accurately captures the separation between translational ( $\theta_{tr}$ ) and internal ( $\theta_{in}$ ) temperatures. It shows significantly better agreement with DSMC data compared to the two-temperature (2T) model, particularly in the non-equilibrium region of the shock.
Limitations: At higher Mach numbers (beyond 4), slight deviations appear in the downstream translational temperature profiles. The authors attribute this to the assumption of a constant relaxation parameter $\delta$ (implying temperature-independent specific heat capacities), which becomes restrictive under strong temperature gradients where energy redistribution becomes temperature-dependent.

Significance and Claims
The authors modestly claim that the primary contribution of this work is not necessarily superior predictive accuracy over all existing models, but rather the demonstration of a transparent, thermodynamically consistent structure for polyatomic gas dynamics.

The GENERIC framework explicitly identifies reversible and irreversible contributions, entropy flux, and translational–internal energy coupling.
The model automatically satisfies conservation laws and the second law of thermodynamics due to its built-in structure.
The agreement with DSMC and experimental data validates the specific 11-moment closure developed within this framework, rather than the GENERIC structure itself.
The work establishes a foundation for future extensions, such as incorporating temperature-dependent relaxation parameters, while maintaining the rigorous thermodynamic structure.

More like this