Rate constants and molecular recombination pathways of oxygen from quasi-classical trajectory simulations of the O3 system

Molecular Dynamics Simulation of Thermal Nonequilibrium and Chemical Reaction Processes in Hydrogen Combustion

Using reactive force field (ReaxFF) and molecular dynamics simulation, we investigate the combustion process of hydrogen-oxygen systems in initial thermal nonequilibrium states with different translational and rovibrational temperatures for oxygen. The system studied in this work contains 300 oxygen molecules and 700 hydrogen molecules with a density of 7 times the air density. For this system, the characteristic relaxation times of oxygen and hydrogen vibrational energies are 0.173 and 0.249 ns, respectively. 0.6% of hydrogen undergoes a chemical reaction with oxygen during the thermal nonequilibrium relaxation stage. For the distribution of translational energy and vibrational energy of oxygen in the thermal nonequilibrium state, the maximum mean error of the statistical distribution in the simulation and the Boltzmann distribution at temperature calculated from the average kinetic energy of molecules is about 2.25 × 10-5. At the same time, it was observed in the simulation that many-body interactions play a certain role in the combustion process. Furthermore, we compare the ignition time and temperature rise behavior of different combustion mechanisms and molecular dynamics simulations starting from the thermal equilibrium state. These results will provide meaningful references for the construction of thermal nonequilibrium combustion chemical reaction mechanisms.

Quasi-classical trajectory analysis of three-body collision induced recombination in neutral nitrogen and oxygen

We present theory and a simulation framework to model three-body collisions and gas phase recombination in dilute atom/diatom mixtures of pure oxygen (O/O2) and nitrogen (N/N2) using the Quasi-Classical Trajectory method. We formulate a three-body collision rate constant based on the lifetimes of binary collisions and initialize three-body collisions by sampling the arrival time of a third body within the lifetimes of pre-simulated binary collisions. We use this method to calculate distributions of recombined product energies, probabilities of recombination, and recombination rate constants through different collision pathways. Long-lived binary atom-diatom collisions are observed, but are too rare to play a dominant role in the recombination process for shock-heated air near the equilibrium conditions studied. The resulting recombination rate constants are within an order of magnitude of the predictions of detailed balance. Notably, the recombination simulation framework does not appeal to the principle of detailed balance and could be useful for studying conditions far from equilibrium.