Non-Isothermal Maxwell-Stefan Equations

• Non-isothermal Maxwell-Stefan equations extend classical diffusion by coupling mass transport with temperature dynamics, incorporating thermal gradients into the driving forces.
• They are derived from kinetic models such as the Boltzmann and BGK frameworks, leading to complex closure relations and energy balances.
• The systems display rich analytical properties, including entropy-based symmetrization and rigorous well-posedness under diverse thermodynamic conditions.

Non-isothermal Maxwell–Stefan equations are multicomponent diffusion systems in which the species balances are coupled to a temperature or energy balance, so that temperature is an unknown of the transport problem rather than a fixed parameter. In contrast with the classical isothermal theory, the diffusive driving forces involve temperature-weighted composition gradients such as $\nabla(c_i T)$, $\nabla(\rho_i\theta)$, or $\nabla(n_i\theta)$, and the thermal equation may contain Fourier conduction, enthalpy transport by species diffusion, or purely advective transport, depending on the constitutive setting. The subject connects kinetic theory, continuum thermodynamics, and cross-diffusion PDE analysis: formal and rigorous derivations from Boltzmann and BGK models are available, as are local and global well-posedness results for several non-isothermal variants (Hutridurga et al., 2017, Georgiadis et al., 2023, Chen et al., 5 Aug 2025).

1. Defining structure and principal formulations

Across the literature, the unknowns are species concentrations, mass densities, or number densities—typically denoted by $c_i$, $\rho_i$, or $n_i$—together with diffusive fluxes $J_i$ or species velocities $u_i$, and a common temperature $T$ or $\theta$. The common structural core is a family of species continuity equations,

$\nabla(\rho_i\theta)$0

closed by Maxwell–Stefan friction relations expressing balance between thermodynamic driving forces and pairwise interspecies drag. In non-isothermal models, those driving forces are no longer purely compositional.

A representative algebraic closure is

$\nabla(\rho_i\theta)$1

with $\nabla(\rho_i\theta)$2 after thermodynamic specialization and with a separate energy balance

$\nabla(\rho_i\theta)$3

in the vanishing-barycentric-velocity formulation (Georgiadis et al., 2023). In other derivations, the macroscopic flux law is written directly in Maxwell–Stefan form,

$\nabla(\rho_i\theta)$4

or, at the number-density level,

$\nabla(\rho_i\theta)$5

again coupled to an energy equation (Hutridurga et al., 2017, Jüngel et al., 17 Sep 2025).

Formulation	Diffusive driving force	Thermal equation or flux
Monatomic Boltzmann limit	$\nabla(\rho_i\theta)$6	$\nabla(\rho_i\theta)$7 (Hutridurga et al., 2017)
Zero-barycentric-velocity model	$\nabla(\rho_i\theta)$8	$\nabla(\rho_i\theta)$9 (Georgiadis et al., 2023)
BGK diffusion limit	$\nabla(n_i\theta)$0	convective plus diffusive heat flux with $\nabla(n_i\theta)$1 (Jüngel et al., 17 Sep 2025)
Maxwell–Stefan–Fourier simplification	$\nabla(n_i\theta)$2	$\nabla(n_i\theta)$3 after normalization (Georgiadis, 11 Jun 2026)

This variety is substantive rather than notational. The literature contains several structurally distinct non-isothermal Maxwell–Stefan systems, mainly differing in the closure imposed on the total flux, the pressure or barycentric constraint, and the constitution of the heat flux.

2. Kinetic origins and asymptotic derivations

A central line of development derives non-isothermal Maxwell–Stefan equations from multispecies kinetic models under diffusive scaling. For monatomic gaseous mixtures, one starts from a system of Boltzmann equations

$\nabla(n_i\theta)$4

assumes that the solution remains close to a specieswise local Maxwellian with small velocities of order $\nabla(n_i\theta)$5, and passes formally to the limit $\nabla(n_i\theta)$6. In the derivation for Maxwellian molecules, pairwise elastic collisions conserve mass, total momentum, and kinetic energy, and the asymptotic analysis yields both species balances and a temperature equation. A key consequence is the synchronization of species temperatures: $\nabla(n_i\theta)$7 so that all species converge to a common limit temperature $\nabla(n_i\theta)$8 (Hutridurga et al., 2017).

A related derivation replaces the local-Maxwellian ansatz by a maximum entropy principle. In that setting, the scaled Boltzmann system for a mixture of monatomic, inert, ideal gases is combined with constrained entropy maximization, yielding an exact scaled Maxwellian

$\nabla(n_i\theta)$9

The resulting momentum balance contains the non-isothermal driving term $c_i$0, and the diffusion limit produces the system

$c_i$1

together with

$c_i$2

The same framework yields an explicit kinetic entropy production quadratic in the relative velocities $c_i$3 (Anwasia et al., 2021).

The first rigorous Boltzmann-to-non-isothermal-Maxwell–Stefan asymptotic result appears in the small-data setting of the multi-species Boltzmann equations under diffusive scaling,

$c_i$4

There, the Maxwellian built from the non-isothermal Maxwell–Stefan solution is not an exact local equilibrium for the mixture collision operator, because cross-species interactions produce a nonzero defect $c_i$5. The analysis therefore combines perturbation around global equilibrium, a local coercivity estimate for the operator linearized around the non-equilibrium local Maxwellian, and uniform-in-$c_i$6 hypocoercive estimates (Chen et al., 5 Aug 2025).

Formal BGK diffusion limits yield further variants. For the multispecies Gross–Krook model, the leading-order temperatures again synchronize,

$c_i$7

and the limit system contains the Maxwell–Stefan-type balance

$c_i$8

together with an energy equation whose flux contains both a convective term with $c_i$9 and a diffusive term

$\rho_i$0

This differs from non-isothermal derivations whose energy balance has only convective transport by the species velocities (Jüngel et al., 17 Sep 2025).

3. Closure relations, pressure constraints, and reduced temperature dynamics

Non-isothermal Maxwell–Stefan systems are typically underdetermined before an additional closure is imposed. In the monatomic formal Boltzmann limit, the limiting variables are $\rho_i$1, giving $\rho_i$2 unknowns but only $\rho_i$3 independent equations. In the isothermal theory, one usually imposes $\rho_i$4, but in the non-isothermal derivation based on $\rho_i$5 the authors propose instead

$\rho_i$6

calling it the decoupling closure relation. When $\rho_i$7 is constant, this reduces to the classical isothermal constraint $\rho_i$8 (Hutridurga et al., 2017).

Under that closure, the total concentration solves

$\rho_i$9

while the temperature obeys an advection-type equation. In the PDE analysis of the associated cross-diffusion system, the temperature equation is written as

$n_i$0

and the temperature is represented along characteristics of the velocity field

$n_i$1

This decoupling is the technical basis for the reduction of the original Maxwell–Stefan system to a quasi-linear parabolic system for $n_i$2 concentrations (Hutridurga et al., 2017).

A different closure appears in thermodynamically consistent zero-barycentric-velocity models. There, the mixture velocity is fixed to $n_i$3, the compatibility constraint is

$n_i$4

and thermodynamic consistency requires the pressure constraint

$n_i$5

Because the barycentric velocity vanishes, the total density is then time-independent: $n_i$6 The algebraic Maxwell–Stefan system is singular, and inversion is performed on the subspace orthogonal to the density vector rather than by an additional Fick-type law for the total density (Georgiadis et al., 2023).

In the rigorous Boltzmann asymptotics for small perturbations, the closure

$n_i$7

again produces a heat equation for $n_i$8,

$n_i$9

while the relation $J_i$0 allows the temperature equation to be rewritten as

$J_i$1

Different closures therefore lead to materially different thermal equations: transport-type in some settings, nonlinear parabolic in others (Chen et al., 5 Aug 2025).

4. Thermodynamic structure, entropy variables, and dissipation

The mathematically most developed non-isothermal Maxwell–Stefan theories are built from a free-energy closure. In the vanishing-barycentric-velocity model, the Helmholtz free energy density of species $J_i$2 is

$J_i$3

from which one obtains

$J_i$4

With these constitutive laws the driving force simplifies to $J_i$5, and the heat flux becomes

$J_i$6

(Georgiadis et al., 2023).

The singular Maxwell–Stefan algebraic system is then rewritten as $J_i$7, where the matrix $J_i$8 has kernel $J_i$9. The Bott–Duffin inverse $u_i$0 exists and is symmetric positive definite on the constraint subspace

$u_i$1

which yields $u_i$2. This inversion underlies the entropy-variable formulation.

The mathematical entropy density is

$u_i$3

and the entropy variables are

$u_i$4

In these variables the fluxes can be written as

$u_i$5

with full Onsager matrix

$u_i$6

positive semidefinite. This positive semidefiniteness is the structural reason the PDE behaves as a nonlinear parabolic system in entropy variables (Georgiadis et al., 2023).

Entropy production is also explicit in kinetic derivations. For the maximum-entropy Boltzmann closure, the total entropy production has the form

$u_i$7

which is nonnegative and identifies Maxwell–Stefan friction as a kinetic dissipation mechanism (Anwasia et al., 2021). In the BGK diffusion limit, under the normalization $u_i$8, the macroscopic entropy production splits into a diffusion/friction term involving $u_i$9 and a thermal term involving $T$0 (Jüngel et al., 17 Sep 2025).

5. Analytical results and well-posedness theory

The first PDE well-posedness results for a non-isothermal Maxwell–Stefan system were local in time. Starting from the system

$T$1

one eliminates one species, rewrites the flux law as $T$2, analyzes the singular Maxwell–Stefan matrix by Perron–Frobenius arguments, and reduces the problem to a quasi-linear parabolic system

$T$3

Because the reduced diffusion matrix has spectrum bounded away from zero and $T$4 remains positive, the system is normally elliptic, and $T$5-maximal regularity yields local-in-time existence and uniqueness of strong solutions (Hutridurga et al., 2017).

Global weak existence and weak-strong uniqueness have been proved for a different non-isothermal Maxwell–Stefan system with heat conduction and vanishing barycentric velocity. Under $T$6, positive initial data, fixed total density bounded away from vacuum, and conductivity growth

$T$7

there exists a global weak solution with $T$8 and $T$9 almost everywhere. The proof uses the boundedness-by-entropy method: time discretization by implicit Euler, higher-order regularization, a discrete entropy inequality, and compactness arguments based on Aubin–Lions. If a strong solution exists with bounded positive densities and bounded temperature, then any weak solution with the same initial data coincides with it; the proof proceeds through a relative entropy inequality and Grönwall’s lemma (Georgiadis et al., 2023).

For the small-perturbation closure-based system arising in the rigorous Boltzmann asymptotics, global classical solutions are available in Sobolev spaces for $\theta$0, with positivity preserved for all time. The proof combines coercivity of the Maxwell–Stefan matrix on $\theta$1, the heat equation for $\theta$2, the relation $\theta$3, Sobolev product estimates, and a nonlinear iteration scheme. On the kinetic side, the same work establishes a unique global classical solution to the perturbed Boltzmann equation, uniform in Knudsen number $\theta$4, thereby rigorously justifying the non-isothermal Maxwell–Stefan asymptotic limit (Chen et al., 5 Aug 2025).

A further analytical direction treats degenerate heat conduction in a Maxwell–Stefan–Fourier system with

$\theta$5

In that case the entropy inequality controls only $\theta$6, not $\theta$7 or $\theta$8, so the identification of the heat flux cannot rely on standard entropy methods. The main new tool is a renormalized energy estimate that yields

$\theta$9

which provides compactness of the temperature and permits passage to the limit in

$\nabla(\rho_i\theta)$00

This gives global weak solutions and, in addition, an entropy production defect measure accounting for possible loss of dissipation in the degenerate limit (Georgiadis, 11 Jun 2026).

6. Extensions, variants, and relation to the isothermal theory

Non-isothermal Maxwell–Stefan equations admit several physically distinct extensions. For reactive mixtures of polyatomic gases with continuous internal energy, the diffusive limit of a reactive Boltzmann system yields

$\nabla(\rho_i\theta)$01

together with the non-isothermal Maxwell–Stefan law

$\nabla(\rho_i\theta)$02

and a temperature equation containing both transport and the reaction source $\nabla(\rho_i\theta)$03. Here temperature enters the flux law explicitly through the term $\nabla(\rho_i\theta)$04, and the coefficients $\nabla(\rho_i\theta)$05 as well as the chemical production term $\nabla(\rho_i\theta)$06 are derived from the collision kernels (Anwasia et al., 2019).

A broader thermodynamic generalization arises from a GENERIC/Hamiltonian framework for binary mixtures with multiple partial velocities and partial temperatures. After reduction to a single-temperature, single-momentum description, the generalized diffusion-velocity law contains chemical-potential gradients, temperature gradients, and electric forcing; after reduction to mechanical equilibrium it becomes a non-isothermal variant of the dusty gas model. In this framework, temperature differences between constituents are identified as a possible new source of the Soret coefficient (Vágner et al., 2022).

The relation to the isothermal Maxwell–Stefan theory is precise in many models. Constant temperature reduces $\nabla(\rho_i\theta)$07, $\nabla(\rho_i\theta)$08, or $\nabla(\rho_i\theta)$09 to the classical composition-driven forces, and the closure $\nabla(\rho_i\theta)$10 collapses to $\nabla(\rho_i\theta)$11 when $\nabla(\rho_i\theta)$12 is constant (Hutridurga et al., 2017). The isothermal theory supplies important operator-theoretic and entropy-based baselines: normal ellipticity from Perron–Frobenius analysis for the reduced diffusion operator (Bothe, 2010), and entropy-variable symmetrization for singular Maxwell–Stefan matrices coupled to incompressible flow (Chen et al., 2013).

Two recurrent simplifications from the isothermal literature do not transfer automatically. First, the equimolar or zero-total-flux condition is not universal in non-isothermal settings: some models retain $\nabla(\rho_i\theta)$13 through the zero-barycentric-velocity constraint, whereas others replace it by a Fick-type law for the total concentration (Georgiadis et al., 2023, Hutridurga et al., 2017). Second, the specieswise local Maxwellian used in asymptotic derivations need not be a true local equilibrium of the full mixture collision operator; the rigorous Boltzmann analysis shows that cross-interactions obstruct that property and generate a nontrivial source term in the perturbation equation (Chen et al., 5 Aug 2025).

The modern theory therefore presents non-isothermal Maxwell–Stefan equations not as a single canonical PDE, but as a family of thermodynamically and kinetically motivated cross-diffusion systems. Their unifying feature is the replacement of purely isothermal diffusion by a coupled mass–heat dynamics in which temperature directly enters the driving forces, the entropy production, and the analytic structure of the equations.