Analysis of Moment Closures Using φ-Divergences for Rarefied Dynamics with Binary Collisions and Their Galerkin Discretizations
Abstract
This work introduces a robust deterministic framework for approximating solutions of the Boltzmann equation with binary collisions by discretizing their dependence on time, position, and velocity using Galerkin methods. By employing a family of parametric Galerkin closures based on φ-divergences in velocity space, we derive rigorous hierarchies of moment equations that govern fluid dynamic variables. Addressing the limitation that these closures alone do not guarantee dissipation of a φ-divergence entropy for the true binary collision operator, we restore this property by formulating a compatible approximate collision operator tailored to each closure. This constructed operator intrinsically retains fundamental physical properties essential for high-fidelity flow simulations, including Galilean invariance, exact conservation of mass, momentum, and energy, and strict dissipation of a φ-divergence entropy. Furthermore, we show that the resulting closed moment systems are symmetric-dissipative, yielding Cauchy problems that are well-posed locally in time. To translate this mathematical foundation into an efficient computational tool, we discretize the position and time variables with an entropy-stable discontinuous Galerkin (DG) finite element method. The fully implicit, entropy-stable space-time approach enables time steps far beyond typical CFL-limited step sizes and the direct computation of steady states. The robustness and accuracy of the methodology are verified and validated through numerical simulations on the supersonic nozzle flow of argon, mass flow through a channel, and heat transfer between parallel walls, demonstrating agreement with analytical benchmarks, experimental measurements, and stochastic particle simulations.
Related
- A general synthetic iterative scheme with high-order spatial discretization for rarefied gas mixture flows
- An asymptotic-preserving reduced-order method for parametrised rarefied gas flow by proper generalised decomposition
- A nonoverlapping spectral additive Schwarz method for interior penalty discontinuous Galerkin discretizations of Anisotropic Elliptic Problems
- Sparse-grid discretizations of BGG sequences
- Divergences in Nonpolynomial Quantum Field Theories
BibTeX
@misc{abdelmalik2026analysis,
title = {{Analysis of Moment Closures Using $φ$-Divergences for Rarefied Dynamics with Binary Collisions and Their Galerkin Discretizations}},
author = {Michael R. A. Abdelmalik and Irene M. Gamba and Torsten Keßler and Sergej Rjasanow},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2608.03640},
eprint = {2608.03640},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2608.03640},
}
RIS
TY - UNPB
TI - Analysis of Moment Closures Using φ-Divergences for Rarefied Dynamics with Binary Collisions and Their Galerkin Discretizations
AU - Michael R. A. Abdelmalik
AU - Irene M. Gamba
AU - Torsten Keßler
AU - Sergej Rjasanow
PY - 2026
DA - 2026/09/03
PB - arXiv
DO - 10.48550/arxiv.2608.03640
UR - https://arxiv.org/abs/2608.03640
AB - This work introduces a robust deterministic framework for approximating solutions of the Boltzmann equation with binary collisions by discretizing their dependence on time, position, and velocity using Galerkin methods. By employing a family of parametric Galerkin closures based on φ-divergences in velocity space, we derive rigorous hierarchies of moment equations that govern fluid dynamic variables. Addressing the limitation that these closures alone do not guarantee dissipation of a φ-divergence entropy for the true binary collision operator, we restore this property by formulating a compatible approximate collision operator tailored to each closure. This constructed operator intrinsically retains fundamental physical properties essential for high-fidelity flow simulations, including Galilean invariance, exact conservation of mass, momentum, and energy, and strict dissipation of a φ-divergence entropy. Furthermore, we show that the resulting closed moment systems are symmetric-dissipative, yielding Cauchy problems that are well-posed locally in time. To translate this mathematical foundation into an efficient computational tool, we discretize the position and time variables with an entropy-stable discontinuous Galerkin (DG) finite element method. The fully implicit, entropy-stable space-time approach enables time steps far beyond typical CFL-limited step sizes and the direct computation of steady states. The robustness and accuracy of the methodology are verified and validated through numerical simulations on the supersonic nozzle flow of argon, mass flow through a channel, and heat transfer between parallel walls, demonstrating agreement with analytical benchmarks, experimental measurements, and stochastic particle simulations.
ER -