A Conservative, Entropic Multispecies BGK Model

We derive a conservative multispecies BGK model that follows the spirit of the original, single species BGK model by making the specific choice to conserve species masses, total momentum, and total kinetic energy and to satisfy Boltzmann’s H-Theorem. The derivation emphasizes the connection to the Boltzmann operator which allows for direct inclusion of information from higher-fidelity collision physics models. We also develop a complete hydrodynamic closure via the Chapman-Enskog expansion, including a general procedure to generate symmetric diffusion coefficients based on this model. We numerically investigate velocity and temperature relaxation in dense plasmas and compare the model with previous multispecies BGK models and discuss the trade-offs that are made in defining and using them. In particular, we demonstrate that the BGK model in the NRL plasma formulary does not conserve momentum or energy in general.