Conceitos essenciais
An efficient iterative solver is proposed to implicitly update the moment equations for a multi-species BGK model with collision frequencies that depend on individual species temperatures. The method is proven to be convergent under mild time step restrictions that are independent of the stiffness of the collision operator.
Resumo
The paper presents an efficient iterative solver for the implicit update of the moment equations in a multi-species Bhatnagar-Gross-Krook (M-BGK) model. The key features of the model are:
- The M-BGK model uses a sum of relaxation operators, each with a Maxwellian-like target, to capture binary interactions between species.
- The collision frequencies in the model depend on the individual species temperatures, which are not conserved by the collision dynamics.
- This temperature dependence makes the implicit update of the moment equations challenging, as the average velocities and temperatures need to be solved for implicitly.
The proposed iterative solver, based on a Gauss-Seidel-type (GST) approach, is shown to be convergent under mild time step restrictions that are independent of the stiffness of the collision operator. The analysis proves that the differences in the velocity and temperature iterates satisfy contraction-type bounds, ensuring convergence of the method.
The key steps in the analysis are:
- Decomposing the velocity and temperature iterates into null space and range space components, and bounding the differences in these components separately.
- Deriving bounds on the source terms in the temperature equation that depend on the velocity iterates.
- Establishing a time step restriction that ensures the contraction of the iterates, independent of the stiffness of the collision operator.
The proposed solver allows the implicit treatment of the collision term in the M-BGK model, while maintaining stability and efficiency, even in the fluid regime where the time step is determined by particle advection rather than collisions.