The kinetic modeling of low-pressure ( torr) stationary nitrogen discharges and the corresponding afterglows is reviewed. It is shown that a good description of the overall behavior of nitrogen plasmas requires a deep understanding of the coupling between different kinetics. The central role is played by ground-state vibrationally excited molecules, N2, which have a strong influence on the shape of the electron energy distribution function, on the creation and destruction of electronically excited states, on the gas heating, dissociation and on afterglow emissions. N2 molecules are actually the hinge ensuring a strong link between the various kinetics. The noticeable task done by electronically excited metastable molecules, in particular N2 and N2, is also pointed out. Besides contributing to the same phenomena as vibrationally excited molecules, these electronic metastable states play also a categorical role in ionization. Furthermore, vibrationally excited molecules in high v levels are in the origin of the peaks observed in the flowing afterglow for the concentrations of several species, such as N2, N2, N2 + and electrons, which occur downstream from the discharge after a dark zone as a consequence of the V-V pumping-up mechanism.

This work explains the need for plasma models, introduces arguments for choosing the type of model that better fits the purpose of each study, and presents the basics of the most common nonequilibrium low-temperature plasma models and the information available from each one, along with an extensive list of references for complementary in-depth reading. The paper presents the following models, organised according to the level of multi-dimensional description of the plasma: kinetic models, based on either a statistical particle-in-cell/Monte-Carlo approach or the solution to the Boltzmann equation (in the latter case, special focus is given to the description of the electron kinetics); multi-fluid models, based on the solution to the hydrodynamic equations; global (spatially-average) models, based on the solution to the particle and energy rate-balance equations for the main plasma species, usually including a very complete reaction chemistry; mesoscopic models for plasma–surface interaction, adopting either a deterministic approach or a stochastic dynamical Monte-Carlo approach. For each plasma model, the paper puts forward the physics context, introduces the fundamental equations, presents advantages and limitations, also from a numerical perspective, and illustrates its application with some examples. Whenever pertinent, the interconnection between models is also discussed, in view of multi-scale hybrid approaches.