Numerical schemes for a multi-species quantum BGK model

TL;DR

This paper develops and analyzes a numerical scheme for a multi-species quantum BGK model, addressing the challenges of computing quantum local equilibrium parameters and demonstrating convergence properties.

Contribution

It introduces an IMEX scheme for the quantum BGK model with a novel Lagrange multiplier method for equilibrium parameter updates.

Findings

The scheme effectively handles stiffness in the collision operator.

Convergence of mean velocity and temperature is established both analytically and numerically.

Quantum component inclusion affects the convergence of physical temperature.

Abstract

This work is devoted to the numerical implementation of the quantum Bhatnagar- Gross-Krook (BGK) model for gas mixtures consisting of classical and quantum particles (fermions, bosons). We consider the model proposed by Bae, Klingenberg, Pirner, and Yun in 2021 and implement an Implicit-Explicit (IMEX) scheme due to the stiffness of the collision operator. A major obstacle is updating the parameters of quantum local equilibrium, which requires computing by inverting the relation between density and energy at every grid point in space and time. We address this difficulty by using the Lagrange multiplier method to minimize a potential function subject to constraints defined by specific moment equalities. Moreover, we analyze the convergence of mean velocity and temperature between the species both analytically and numerically. When a quantum component is included, we observe that the converging quantity is physical temperature, not the kinetic temperature. This differs from the mixture of classical species.