State-to-state modeling of a recombining nitrogen plasma experiment

Plasma Modeling and Prebiotic Chemistry: A Review of the State-of-the-Art and Perspectives

We review the recent progress in the modeling of plasmas or ionized gases, with compositions compatible with that of primordial atmospheres. The plasma kinetics involves elementary processes by which free electrons ultimately activate weakly reactive molecules, such as carbon dioxide or methane, thereby potentially starting prebiotic reaction chains. These processes include electron–molecule reactions and energy exchanges between molecules. They are basic processes, for example, in the famous Miller-Urey experiment, and become relevant in any prebiotic scenario where the primordial atmosphere is significantly ionized by electrical activity, photoionization or meteor phenomena. The kinetics of plasma displays remarkable complexity due to the non-equilibrium features of the energy distributions involved. In particular, we argue that two concepts developed by the plasma modeling community, the electron velocity distribution function and the vibrational distribution function, may unlock much new information and provide insight into prebiotic processes initiated by electron–molecule collisions.

Adaptive coarse graining method for energy transfer and dissociation kinetics of polyatomic species

A novel reduced-order method is presented for modeling reacting flows characterized by strong non-equilibrium of the internal energy level distribution of chemical species in the gas. The approach seeks for a reduced-order representation of the distribution function by grouping individual energy states into macroscopic bins, and then reconstructing state population using the maximum entropy principle. This work introduces an adaptive grouping methodology to identify and lump together groups of states that are likely to equilibrate faster with respect to each other. To this aim, two algorithms have been considered: the modified island algorithm and the spectral clustering method. Both methods require a measure of dissimilarity between internal energy states. This is achieved by defining "metrics" based on the strength of the elementary rate coefficients included in the state-specific kinetic mechanism. Penalty terms are used to avoid grouping together states characterized by distinctively different energies. The two methods are used to investigate excitation and dissociation of N₂ (Σg+1) molecules due to interaction with N(Su4) atoms in an ideal chemical reactor. The results are compared with a direct numerical simulation of the state-specific kinetics obtained by solving the master equations for the complete set of energy levels. It is found that adaptive grouping techniques outperform the more conventional uniform energy grouping algorithm by providing a more accurate description of the distribution function, mole fraction and energy profiles during non-equilibrium relaxation.

Energy transfer models in nitrogen plasmas: analysis of N₂(X¹Σg⁺)-N(⁴S(u))-e⁻ interaction

The relaxation of N₂(X¹Σg⁺) molecules in a background gas composed of N((4)S(u)) atoms and free electrons is studied by using an ideal isochoric and isothermic chemical reactor. A rovibrational state-to-state model is developed to study energy transfer process induced by free electron and atomic collisions. The required cross sections and the corresponding rate coefficients are taken from two well-known kinetic databases: NASA Ames kinetic mechanism for the description of the N₂(X¹Σg⁺)-N((4)S(u)) processes and the Phys4Entry database for the electron driven processes, N₂(X¹Σg⁺)-e(-). The evolution of the population densities of each individual rovibrational level is explicitly determined via the numerical solution of the master equation for temperatures ranging from 10000 to 30,000 K. It was found that the distribution of the rovibrational energy levels of N₂(X¹Σg⁺) is strongly influenced by the electron driven collisional processes, which promote the excitation of the low lying vibrational levels. The macroscopic vibrational energy relaxation is governed by the molecule-atom collisions, when free electrons, initially cold are relaxing to the final heat-bath temperature. Thus, the main role of the free electrons is to ensure the equilibration of vibrational and free electron excitation, thus validating the existence of the local equilibrium T(V)-T(e). However, if electrons and heavy particles are assumed to be in equilibrium at the heat bath temperature, electron driven processes dominate the vibrational relaxation. Finally, we have assessed the validity of the Landau-Teller model for the description of the inelastic energy transfer between molecules and free electrons. In the case of free-electron temperatures lower than 10,000 K, Landau-Teller relaxation model gives an accurate description of the vibrational relaxation, while at higher temperatures the error in the predictions can be significant and the model should not be used.

Nonequilibrium shock-heated nitrogen flows using a rovibrational state-to-state method

A rovibrational collisional model is developed to study the internal energy excitation and dissociation processes behind a strong shock wave in a nitrogen flow. The reaction rate coefficients are obtained from the ab initio database of the NASA Ames Research Center. The master equation is coupled with a one-dimensional flow solver to study the nonequilibrium phenomena encountered in the gas during a hyperbolic reentry into Earth's atmosphere. The analysis of the populations of the rovibrational levels demonstrates how rotational and vibrational relaxation proceed at the same rate. This contrasts with the common misconception that translational and rotational relaxation occur concurrently. A significant part of the relaxation process occurs in non-quasi-steady-state conditions. Exchange processes are found to have a significant impact on the relaxation of the gas, while predissociation has a negligible effect. The results obtained by means of the full rovibrational collisional model are used to assess the validity of reduced order models (vibrational collisional and multitemperature) which are based on the same kinetic database. It is found that thermalization and dissociation are drastically overestimated by the reduced order models. The reasons of the failure differ in the two cases. In the vibrational collisional model the overestimation of the dissociation is a consequence of the assumption of equilibrium between the rotational energy and the translational energy. The multitemperature model fails to predict the correct thermochemical relaxation due to the failure of the quasi-steady-state assumption, used to derive the phenomenological rate coefficient for dissociation.