Vibrational Modeling with an Anharmonic Oscillator Model in Direct Simulation Monte Carlo

Vehicles undergoing hypersonic speed experience extreme aerothermodynamic conditions. Real gas effects cannot be neglected, and thus internal degrees of freedom of molecules being partially/fully excited must be carefully predicted in order to accurately capture the physics of the flowfield. Within direct simulation Monte Carlo solvers, a harmonic oscillator (HO) model, where the quantum levels are evenly spaced, is typically used for vibrational energy. A more realistic model is an anharmonic oscillator (aHO), in which the energy between quantum levels is not evenly spaced. In this work, the Morse-aHO model is compared against HO. The Morse-aHO model is implemented in the dsmcFoam+ solver, and the numerical results are in excellent agreement with analytical and potential energy surface solutions for the partition function, mean vibrational energy, and degrees of freedom. A method for measuring the vibrational temperature of the gas when using the anharmonic model in a direct simulation Monte Carlo solver is presented, which is essential for returning macroscopic fields. For important thermophysical properties of molecular oxygen, such as the specific heat capacity, it is shown that the aHO and HO models begin to diverge at temperatures above 1000 K, making the use of HO questionable for all but low-enthalpy flows. For the same gas, including the electronic energy mode significantly improves the accuracy of the specific heat prediction, compared to experimental data, for temperatures above 2000 K. For relaxation from a state of thermal nonequilibrium, it is shown that the aHO model results in a slightly lower equilibrium temperature. When applied to hypersonic flow over a cylinder, the aHO model results in a smaller shock standoff distance and lower peak temperatures.

[1]  Luc Mieussens,et al.  IXV post-flight reconstruction and analysis of the aerothermodynamic measurements along the rarefied portion of the reentry trajectory , 2021 .

[2]  T. Schwartzentruber,et al.  Non-Boltzmann vibrational energy distributions and coupling to dissociation rate. , 2019, The Journal of chemical physics.

[3]  C. Munz,et al.  Combining particle-in-cell and direct simulation Monte Carlo for the simulation of reactive plasma flows , 2019, Physics of Fluids.

[4]  Liwu Liu,et al.  High-temperature partition functions, specific heats and spectral radiative properties of diatomic molecules with an improved calculation of energy levels , 2018 .

[5]  Stephen M. Longshaw,et al.  dsmcFoam+: An OpenFOAM based direct simulation Monte Carlo solver , 2017, Comput. Phys. Commun..

[6]  Manuel Torrilhon,et al.  The 35-Moment System with the Maximum-Entropy Closure for Rarefied Gas Flows , 2017 .

[7]  Thomas E. Schwartzentruber,et al.  Nonequilibrium Gas Dynamics and Molecular Simulation , 2017 .

[8]  D. Griffiths.  Introduction to Quantum Mechanics , 2016 .

[9]  Alina A. Alexeenko,et al.  Ab initio-informed maximum entropy modeling of rovibrational relaxation and state-specific dissociation with application to the O2 + O system. , 2016, The Journal of chemical physics.

[10]  D. Levin,et al.  A study of internal energy relaxation in shocks using molecular dynamics based models. , 2015, The Journal of chemical physics.

[11]  Alina A. Alexeenko,et al.  Effect of O2 + O ab initio and Morse additive pairwise potentials on dissociation and relaxation rates for nonequilibrium flow calculations , 2015 .

[12]  G. Candler,et al.  An improved potential energy surface and multi-temperature quasiclassical trajectory calculations of N2 + N2 dissociation reactions. , 2015, The Journal of chemical physics.

[13]  Marco Panesi,et al.  General multi-group macroscopic modeling for thermo-chemical non-equilibrium gas mixtures. , 2015, The Journal of chemical physics.

[14]  F. Esposito,et al.  N2, O2, NO state-to-state vibrational kinetics in hypersonic boundary layers: The problem of rescaling rate coefficients to uniform vibrational ladders , 2015 .

[15]  T. Schwartzentruber,et al.  Rovibrational coupling in molecular nitrogen at high temperature: An atomic-level study , 2014 .

[16]  Thomas E. Schwartzentruber,et al.  Progress and Future Prospects for Particle-Based Simulation of Hypersonic Flow , 2013 .

[17]  A. Soufiani,et al.  High-Temperature and Nonequilibrium Partition Function and Thermodynamic Data of Diatomic Molecules , 2009 .

[18]  M. Capitelli,et al.  Recombination-Assisted Nitrogen Dissociation Rates Under Nonequilibrium Conditions , 2008 .

[19]  J. Loureiro,et al.  Vibrational distributions in N2 with an improved calculation of energy levels using the RKR method , 2008 .

[20]  J. L. Stollery,et al.  Hypersonic and High-Temperature Gas Dynamics – Second edition J.D. Anderson American Institute of Aeronautics and Astronautics, 1801 Alexander Bell Drive, Suite 500, Reston, VA, USA. 20191-4344. 2006. 811pp. Illustrated. $64.95 (AIAA members), $90.95 (non-members). ISBN 1-56347-780-7. , 2007, The Aeronautical Journal (1968).

[21]  L. Marraffa,et al.  High-Temperature Thermodynamic Properties of Mars-Atmosphere Components , 2004 .

[22]  Katsuhisa Koura,et al.  Direct simulation Monte Carlo study of rotational nonequilibrium in shock wave and spherical expansion of nitrogen using classical trajectory calculations , 2002 .

[23]  I. Adamovich,et al.  Semiclassical modeling of state-specific dissociation rates in diatomic gases , 2000 .

[24]  Hrvoje Jasak,et al.  A tensorial approach to computational continuum mechanics using object-oriented techniques , 1998 .

[25]  Katsuhisa Koura,et al.  Monte Carlo direct simulation of rotational relaxation of nitrogen through high total temperature shock waves using classical trajectory calculations , 1998 .

[26]  Katsuhisa Koura,et al.  4 Carlo direct simulation of rotational relaxation of diatomic molecules using classical trajectory calculations: Nitrogen shock wave , 1997 .

[27]  Stefan Dietrich,et al.  Scalar and Parallel Optimized Implementation of the Direct Simulation Monte Carlo Method , 1996 .

[28]  Stephen M. Ruffin,et al.  Prediction of Vibrational Relaxation in Hypersonic Expanding Flows Part 2: Results , 1995 .

[29]  S. Ruffin,et al.  Prediction of vibrational relaxation in hypersonic expanding flows. I: Model development , 1995 .

[30]  I. Adamovich,et al.  Vibrational relaxation and dissociation behind shock waves. Part 1 - Kinetic rate models. , 1995 .

[31]  K. Koura.  A set of model cross sections for the Monte Carlo simulation of rarefied real gases: Atom–diatom collisions , 1994 .

[32]  James B. Anderson,et al.  Direct Monte Carlo simulation of chemical reaction systems: Internal energy transfer and an energy‐dependent unimolecular reaction , 1993 .

[33]  Iain D. Boyd,et al.  Models for direct Monte Carlo simulation of coupled vibration-dissociation , 1993 .

[34]  Graeme A. Bird,et al.  Definition of mean free path for real gases , 1983 .

[35]  Claus Borgnakke,et al.  Statistical collision model for Monte Carlo simulation of polyatomic gas mixture , 1975 .

[36]  P. V. Marrone,et al.  Chemical Relaxation with Preferential Dissociation from Excited Vibrational Levels , 1963 .

[37]  P. V. Marrone,et al.  Effect of Dissociation on the Rate of Vibrational Relaxation , 1962 .

[38]  J. D. Teare,et al.  Theory of Radiation from Luminous Shock Waves in Nitrogen , 1959 .

[39]  R. L. Roy,et al.  LEVEL: A computer program for solving the radial Schrödinger equation for bound and quasibound levels , 2017 .

[40]  Hugh M. Hulburt,et al.  Potential Energy Functions for Diatomic Molecules , 1941 .