The importance of O3 excited potential energy surfaces in O2-O high-temperature kinetics

Gaussian Process Regression for State-to-State Integral Cross Sections: The Case of the O + O2 Collision Dissociation Reactions

Research on hypersonic vehicles has become increasingly important worldwide in recent years. However, accurately simulating the dynamics of the nonequilibrium high-temperature reactions that are in the hypersonic flow around the vehicles presents a significant challenge as a large number of states and transitions are accessible even for the smallest atom-diatom reaction systems. It is quite difficult, sometimes even impossible, to exhaustively investigate all relevant combinations or determine high-dimensional analytical representations for the state-to-state reaction probabilities. In this study, we used Gaussian process regression (GPR) to fit a model based on only 807 QCT data for training. The confidence interval of the GPR prediction and the Kullback-Leibler (KL) divergence were used to help minimize the sampling amount of data for fitting the converged GPR model. The model aims to predict the state-to-state integral cross section (ICS) of the O + O2 → 3O dissociation reaction under random initial conditions (Et, v, j). In total, it took almost a month to obtain this converged GPR model, but it took only a few seconds to predict the ICS value for any initial condition. For 330 initial conditions not included in the training set, the mean-square error (MSE) between the QCT-calculated ICSs and the GPR-predicted ones is only 0.08 Å2 and the R2 is 0.9986, indicating that the GPR model can replace the direct expensive QCT calculation with high accuracy. Finally, we calculated the equilibrium dissociation rate coefficients based on the StS ICS values predicted by the GPR model, and the results were in good agreement with available experimental and theoretical results. Thus, this study provides an effective and accurate approach to the extensive direct state-to-state reaction dynamic calculations.

State-to-state study of non-equilibrium recombination of oxygen and nitrogen molecules

Rapidly cooled mixtures are of interest for several applications, including hypersonic flows due to the presence of strong cooling temperature gradients in regions such as hypersonic boundary layers and expanding nozzles. There have been very few studies of rapidly cooled mixtures using the high-fidelity rovibrational databases afforded by ab initio potential energy surfaces. This work makes use of existing rovibrational state-specific databases to study rapidly cooled mixtures. In particular, we seek to understand the importance of thermal non-equilibrium in recombining mixtures using both rovibrational and vibrational state-to-state methods for oxygen and nitrogen molecules. We find that although there is significant non-equilibrium during recombination, it is well captured by the vibrational state-specific approach. Finally, we compare the global recombination rate computed based on the state-specific recombination rate coefficients and the global recombination rate computed based on the time local dissociation rate coefficient, which is reversed using the principle of detailed balance. The local dissociation rate coefficient is computed by weighting the state-specific dissociation rate coefficients with the state-specific distribution of energy states. We find a large difference between these rates, highlighting a potential source of errors in hypersonic flow predictions.

Rovibrationally state-specific collision model for the O2(Σg-3) + O(P3) system in DSMC

A rovibrationally state-specific collision model for the O₂(Σg-3)+O(P3) system is presented for direct simulation Monte Carlo, including rotation-vibration-translation energy transfer, exchange, dissociation, and recombination processes. The two-step binary collision approach is employed to model recombination reactions. Two available cross section databases by Andrienko/Boyd and Esposito/Capitelli are employed for the rovibrationally resolved model (rv-STS) and vibrationally resolved model (v-STS), respectively. The difference between rv-STS and v-STS comes from two contributions: the multisurface factor of dissociation (f_MS) and the rotational averaging process. The dissociation cross section with the constant f_MS is typically larger than with the variable f_MS, especially for the low vibrational energy states. On the other hand, the cross sections resulting from the rotationally averaged database are found to underpredict the dissociation rate coefficient at low temperatures. In the rovibrational heating case, the rv-STS predicts faster relaxation than the v-STS, which also shows a lower quasi-steady-state temperature than v-STS. In the rovibrational cooling case, the rv-STS shows a faster relaxation than v-STS, which also presents a thermal non-equilibrium between rovibrational and translational mode during the cooling process.

Energy exchange rate coefficients from vibrational inelastic O2(Σg-3) + O2(Σg-3) collisions on a new spin-averaged potential energy surface

A new spin-averaged potential energy surface (PES) for non-reactive O₂(Σg-3) + O₂(Σg-3) collisions is presented. The potential is formulated analytically according to the nature of the principal interaction components, with the main van der Waals contribution described through the improved Lennard-Jones model. All the parameters involved in the formulation, having a physical meaning, have been modulated in restricted variation ranges, exploiting a combined analysis of experimental and ab initio reference data. The new PES is shown to be able to reproduce a wealth of different physical properties, ranging from the second virial coefficients to transport properties (shear viscosity and thermal conductivity) and rate coefficients for inelastic scattering collisions. Rate coefficients for the vibrational inelastic processes of O₂, including both vibration-to-vibration (V-V) and vibration-to-translation/rotation (V-T/R) energy exchanges, were then calculated on this PES using a mixed quantum-classical method. The effective formulation of the potential and its combination with an efficient, yet accurate, nuclear dynamics treatment allowed for the determination of a large database of V-V and V-T/R energy transfer rate coefficients in a wide temperature range.

Data-Inspired and Physics-Driven Model Reduction for Dissociation: Application to the O2 + O System

This work presents an in-depth discussion on the nonequilibrium dissociation of O₂ molecules colliding with O atoms, combining quasi-classical trajectory calculations, master equation, and dimensionality reduction. A rovibrationally resolved database for all of the elementary collisional processes is constructed by including all nine adiabatic electronic states of O₃ in the QCT calculations. A detailed analysis of the ab initio data set reveals that for a rovibrational level, the probability of dissociating is mostly dictated by its deficit in internal energy compared to the centrifugal barrier. Because of the assumption of rotational equilibrium, the conventional vibrational-specific calculations fail to characterize such a dependence. Based on this observation, a new physics-based grouping strategy for application to coarse-grained models is proposed. By relying on a hybrid technique made of rovibrationally resolved excitation coupled to coarse-grained dissociation, the new approach is compared to the vibrational-specific model and the direct solution of the rovibrational state-to-state master equation. Simulations are performed in a zero-dimensional isothermal and isochoric chemical reactor for a wide range of temperatures (1500-20,000 K). The study shows that the main contribution to the model inadequacy of vibrational-specific approaches originates from the incapability of characterizing dissociation, rather than the energy transfers. Even when constructed with only twenty groups, the new reduced-order model outperforms the vibrational-specific one in predicting all of the QoIs related to dissociation kinetics. At the highest temperature, the accuracy in the mole fraction is improved by 2000%.