Monte Carlo simulation of excitation and ionization collisions with complexity reduction

Kinetic simulation of plasmas with detailed excitation and ionization collisions presents a significant computational challenge due to the multiscale feature of the collisional rates. In the present work, we propose a complexity reduction method based on atomic level grouping for modeling excitation and ionization collisions. High order of accuracy of the reduction method is realized by allowing an internal distribution within each group. We apply the reduction method to the standard Monte Carlo collision algorithm to model an atomic Hydrogen plasma. Numerical results suggest that the stiffness of the collisional kinetics can be significantly reduced with minimal loss in accuracy.