Modeling and Analysis of Chemical Kinetics for Hypersonic Flows in Air

Abstract

Gas-phase chemical kinetics are relevant for hypersonic flows, but they are currently modeled in CFD using empirical assumptions and decades-old experimental data. Recent advances in quantum chemistry have enabled the construction of accurate potential energy surfaces (PESs) for diatom-diatom interactions in air. Using these PESs, a database of simulated interactions is generated and analyzed; N2 + N2, N2 + N, N2 + O2, O2 + O2, and O2 + O reactions are considered. The conditions studied range from 4000 K to 30,000 K and include thermal equilibrium and nonequilibrium test sets. The nitrogen dissociation rate is found to be similar for collision partners N2, N, and O2. The oxygen dissociation rate, in contrast, is moderately dependent on partner species; O2 is approximately 2 to 3 times more effective than partner N2. Oxygen dissociation with partner N2 is therefore found to be substantially overpredicted by current CFD models, which is consistent with the limited experimental data available for this reaction. The presence of N is known from experiments to promote nitrogen dissociation; this augmentation is found to be due to increased vibrational relaxation, rather than an increased dissociation rate as described by current CFD models. Similar observations are made for oxygen dissociation with partner O, due to a combination of vibrational and electronic energy relaxation. Using only the shock tube data that informed popular CFD models, it was impossible to isolate the effect of increased relaxation from increased dissociation. The change in vibrational energy per dissociation, a necessary input to CFD, is found to be very sensitive to the degree of thermal nonequilibrium. This dependence is not well predicted by any existing chemical kinetics models; correctly describing this term fundamentally changes the thermal evolution of a gas in CFD. The mechanics of dissociation are similar for all reactions studied, so a series of aggregate analyses on all dissociation reactions is performed. Vibration is found to have a more pronounced effect on dissociation than rotation, due to rotation increasing the centrifugal barrier. The classic Marrone-Treanor preferential dissociation model is found to accurately describe all data in the nonequilibrium test sets, but it neglects the effect of rotational energy on dissociation. A modified model is proposed that describes rates to within 22% and vibrational energy changes to within 4% of the dissociation energy, for all dissociation reactions and conditions. For this work, we have considered Boltzmann or approximately Boltzmann distributions, but the population of high-energy molecules is known to be depleted in a dissociating gas ensemble. Various kinetics models based on Boltzmann distributions are implemented in US3D, a production CFD solver designed for hypersonic flows. As expected, the dissociation rate is overpredicted compared to the benchmark data. Work remains, therefore, to account for the non-Boltzmann distributions that exist in reality. These data and insights about dissociation can form the basis for next-generation chemical kinetics models for CFD.