The half-space problem of evaporation and condensation for polyatomic gases and entropy inequalities

This study investigates the steady Boltzmann equation in one spatial variable for a polyatomic single-component gas in a half-space. Inflow boundary conditions are assumed at the half-space boundary, where particles entering the half-space are distributed as a Maxwellian, an equilibrium distribution characterized by macroscopic parameters of the boundary. At the far end, the gas tends to an equilibrium distribution, which is also Maxwellian. Using conservation laws and an entropy inequality, we derive relations between the macroscopic parameters of the boundary and at infinity required for the existence of solutions. The relations vary depending on the sign of the Mach number at infinity, which dictates whether evaporation or condensation takes place at the interface between the gas and the condensed phase. We explore the obtained relations numerically. This investigation reveals that, although the relations are qualitatively similar for various internal degrees of freedom of the gas, clear quantitative differences exist.

💡 Research Summary

The paper addresses the steady Boltzmann equation in a half‑space for a single‑component polyatomic gas, extending classical monatomic analyses to gases with internal degrees of freedom. The authors consider a one‑dimensional spatial domain x>0, with the distribution function f(x,ξ,I) depending on three velocity components ξ and a continuous internal energy variable I≥0. The internal structure is modeled by a power‑law density of states φ(I)=I^{δ/2‑1}, where δ>0 denotes the number of internal degrees of freedom (δ=0 for monatomic, δ=2 for linear diatomic, δ=3 for spherical tops, etc.).

Boundary conditions are of inflow type at the wall (x=0): particles entering the domain follow a Maxwellian M₀⁺ with temperature T₀, density n₀ and zero bulk velocity. At infinity (x→∞) the solution is required to approach another Maxwellian M∞ characterized by temperature T∞, density n∞ and a bulk velocity u∞ directed along the x‑axis. The Mach number at infinity, M∞₁ = u∞/c (c being the sound speed computed from T∞ and δ), determines whether the interface experiences evaporation (M∞₁>0) or condensation (M∞₁<0).

The core of the analysis relies on two fundamental properties of the Boltzmann collision operator: (i) the five conservation laws (mass, three components of momentum, total energy) and (ii) the H‑theorem, which guarantees non‑negative entropy production and identifies Maxwellians as the only zero‑entropy‑production states. By integrating the entropy inequality across the half‑space and using the conservation laws, the authors obtain algebraic relations linking the macroscopic parameters (n₀, T₀, p₀) at the wall to those at infinity (n∞, T∞, p∞).

To make the problem tractable, they first linearize around the far‑field Maxwellian: f = M∞ + √M∞ F, discard quadratic terms, and introduce a possible source term S that is orthogonal to the collision invariants. The resulting linear kinetic equation (ξ₁+u∞)∂ₓF + L F = S, together with the inflow condition F(0,ξ,I)=F₀ for ξ₁+u∞>0 and the decay condition F→0 as x→∞, is analyzed using the spectral theory of the linearized operator L. L is self‑adjoint, non‑negative, and its null space is exactly the five‑dimensional space spanned by the collision invariants. Existing results for the linear half‑space problem (Bernhoff 2022) provide existence, uniqueness, and a classification of solutions in terms of the number of free parameters.

The entropy‑based constraints lead to four distinct regimes, determined solely by the sign and magnitude of M∞₁:

1. Supersonic evaporation (M∞₁>1) – three independent algebraic relations must hold, leaving no free macroscopic parameter.
2. Subsonic evaporation (0≤M∞₁<1) – two relations, one free parameter (e.g., pressure ratio).
3. Subsonic condensation (‑1<M∞₁<0) – a single relation, two free parameters.
4. Supersonic condensation (M∞₁≤‑1) – no algebraic restriction; three parameters can be chosen independently, subject only to physical admissibility (positivity, etc.).

These regimes mirror those found for monatomic gases, confirming that the qualitative structure is robust against the inclusion of internal energy. However, the quantitative bounds differ with δ. The authors compute the admissible domains for the pressure ratio p₀/p∞ and temperature ratio T₀/T∞ for several values of δ (0, 2, 3) by numerically minimizing the entropy production functional. The “maximal entropy production” surfaces—curves in the (p₀/p∞, T₀/T∞) plane for a fixed Mach number—shift noticeably as δ increases: higher internal degrees of freedom allow larger pressure differences for the same Mach number, and the condensation region expands. This reflects the fact that internal modes can store kinetic energy, thereby relaxing the constraints imposed by the external flow.

Importantly, the derivations do not depend on a specific collision kernel; only the micro‑reversibility and positivity assumptions are required. Thus the results hold for a broad class of intermolecular potentials, including hard‑sphere‑like models and more general power‑law kernels.

The paper concludes that entropy inequalities combined with conservation laws provide a powerful, kernel‑independent framework for characterizing half‑space problems of polyatomic gases. The numerical illustrations demonstrate that while the overall topology of admissible parameter sets remains unchanged, internal degrees of freedom produce measurable quantitative effects that must be accounted for in applications such as evaporation cooling, high‑altitude aerodynamics, and plasma–wall interactions. Future work is suggested on extending the analysis to time‑dependent (non‑steady) settings, partially reflective boundary conditions, and rigorous existence proofs for the full nonlinear problem.