A reduced-order NLTE kinetic model for radiating plasmas of outer envelopes of stellar atmospheres

TL;DR

This paper introduces a simplified, reduced-order NLTE kinetic model for radiating plasmas in stellar atmospheres, effectively balancing accuracy and computational efficiency by grouping energy states.

Contribution

It develops a novel grouping strategy for energy states in NLTE models, enabling accurate simulations with fewer unknowns in stellar plasma applications.

Findings

Maxwell-Boltzmann grouping maintains accuracy with fewer unknowns.

Two groups suffice for accurate modeling without line radiation.

Adding one group captures line radiation effects effectively.

Abstract

The present work proposes a self-consistent reduced-order NLTE kinetic model for radiating plasmas such as are found in the outer layers of stellar atmospheres. Starting from the most up-to-date set of ab-initio and experimental data, the highly complex collisional-radiative kinetic mechanism is simplified by lumping the bound energy states in groups. Different grouping strategies are investigated, such as uniform and Maxwell-Boltzmann. The reduced set of governing equations for the material gas and the radiation field is obtained based on a moment method. Applications consider the steady flow across a shock wave in partially ionized hydrogen. The results show that adopting a Maxwell-Boltzmann grouping allows, on the one hand, for a substantial reduction of the number of unknowns and, on the other, to maintain accuracy for both gas and radiation quantities. It is observed that, when neglecting line radiation, the use of two groups already leads to a very accurate resolution of the photo-ionization precursor, internal relaxation and radiative cooling regions. The inclusion of line radiation requires adopting just one additional group to account for optically thin losses in the alpha, beta and gamma lines of the Balmer and Paschen series. This trend has been observed for a wide range of shock wave velocities.