A Multispecies ESBGK Model for Gas Mixtures with Variable Hard Sphere Transport: Theory and Verification

Abstract

A multi-species Bhatnagar-Gross-Krook (BGK) model for gas mixtures is presented that achieves the correct species-wise relaxation of velocities, temperatures, and pressure tensors according to the Boltzmann collision integral, as well as the correct mixture Prandtl number, while retaining a single relaxation term per species. The model extends the ellipsoidal statistical BGK (ESBGK) model by introducing relative relaxation targets for each species, derived from the Variable Hard Sphere (VHS) production rates of the Grad 13 approximation. Three approaches for the species relaxation frequency are proposed and analyzed: a Grad 13-based per-species frequency, a mixture-averaged frequency, and an empirical harmonic mean of the two. The model is implemented in the particle-based code PICLas and verified against Direct Simulation Monte Carlo (DSMC) results for a range of test cases, including 0D reservoir relaxation, mass diffusion, supersonic Couette flow, and hypersonic flow around a 70° blunted cone for binary and ternary gas mixtures. Across all test cases, the proposed model reproduces the correct Prandtl number, species temperature, velocity relaxation rates and pressure tensor relaxation, with the empirical relaxation frequency consistently yielding the best agreement with DSMC.