Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ab Initio Kinetics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

16 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Classical trajectory studies of gas phase reaction dynamics and kinetics using ab initio potential energy surfaces

Strategies for constructing global potential energy surfaces from a limited number of accurate ab initio electronic energy calculations are discussed. Generally, these data are concentrated in small regions of configuration space (e.g., in the vicinity of saddle points and energy minima) and difficulties arise in generating a potential function that is globally well-behaved. Efficient computer codes for carrying out classical trajectory calculations on vector and parallel processors are also described. Illustrations are given from recent work on the following chemical systems: Ca + HF yields CaF + H, H + H + H2 yields H2 + H2, N + O2 yields NO + O and O + N2 yields NO + N. The dynamics and kinetics of metathesis, dissociation, recombination, energy transfer and complex formation processes will be discussed.

Jaffe, Richard L.↗
NASA NTRSDate not supplied

Constraints on Triple Oxygen Isotope Kinetics

Triple oxygen isotope methods (e.g., using δ 18 O, and Δ′ 17 O) have been used to classify meteorites, to constrain atmospheric chemistry, and to estimate primary productivity, among other uses. Such utility relies on the mass dependency of most chemical reactions: δ 17 O is proportional to δ 18 O by roughly a factor of 0.5. However, in detail, mass-dependency factors vary. While constraints on mass dependency are well established for equilibrium isotope effects, similar constraints for kinetic isotope effects are lacking. Using an ab initio Monte Carlo approach (see envelopes in Figure 1), we find the range of kinetic isotope effects (i.e., resulting from unidirectional chemical reactions) contains and exceeds the range of accessible equilibrium isotope effects for the 16 O 17 O 18 O system (Hayles and Killingsworth, 2022). Accessible kinetic isotope effects can yield Δ′ 17 O variations as large 0.8 ‰ from a single kinetic fractionation step. Reactions that yield the largest deviations from equilibrium are those with the largest molecular mass difference between decomposition fragments, typically the two products such as a gas and solid generated from the thermal decomposition of an initial solid. Comparison of DFT transition state model results against previously published experiments on the thermal decomposition of calcite (CaCO 3 ) (Miller et al., 2002) and dehydroxylation of brucite (Mg(OH) 2 ) (Clayton and Mayeda, 2009) show they are matched with mass-dependent kinetic isotope effects. This approach can be applied to any system with more than two isotopes.

Oxygen isotopes↗
NASA NTRS1985

Dynamic Creep Buckling: Analysis of Shell Structures Subjected to Time-dependent Mechanical and Thermal Loading

The objective of the present research is to develop a general mathematical model and solution methodologies for analyzing the structural response of thin, metallic shell structures under large transient, cyclic, or static thermomechanical loads. Among the system responses associated with these loads and conditions are thermal buckling, creep buckling, and ratcheting. Thus geometric and material nonlinearities (of high order) can be anticipated and must be considered in developing the mathematical model. A complete, true ab-initio rate theory of kinematics and kinetics for continuum and curved thin structures, without any restriction on the magnitude of the strains or the deformations, was formulated. The time dependence and large strain behavior are incorporated through the introduction of the time rates of metric and curvature in two coordinate systems: fixed (spatial) and convected (material). The relations between the time derivative and the covariant derivative (gradient) were developed for curved space and motion, so the velocity components supply the connection between the equations of motion and the time rates of change of the metric and curvature tensors.

Simitses, G. J.↗
NASA NTRS2011

Kinetic Monte Carlo Simulation of Oxygen and Cation Diffusion in Yttria-Stabilized Zirconia

Yttria-stabilized zirconia (YSZ) is of interest to the aerospace community, notably for its application as a thermal barrier coating for turbine engine components. In such an application, diffusion of both oxygen ions and cations is of concern. Oxygen diffusion can lead to deterioration of a coated part, and often necessitates an environmental barrier coating. Cation diffusion in YSZ is much slower than oxygen diffusion. However, such diffusion is a mechanism by which creep takes place, potentially affecting the mechanical integrity and phase stability of the coating. In other applications, the high oxygen diffusivity of YSZ is useful, and makes the material of interest for use as a solid-state electrolyte in fuel cells. The kinetic Monte Carlo (kMC) method offers a number of advantages compared with the more widely known molecular dynamics simulation method. In particular, kMC is much more efficient for the study of processes, such as diffusion, that involve infrequent events. We describe the results of kinetic Monte Carlo computer simulations of oxygen and cation diffusion in YSZ. Using diffusive energy barriers from ab initio calculations and from the literature, we present results on the temperature dependence of oxygen and cation diffusivity, and on the dependence of the diffusivities on yttria concentration and oxygen sublattice vacancy concentration. We also present results of the effect on diffusivity of oxygen vacancies in the vicinity of the barrier cations that determine the oxygen diffusion energy barriers.

Good, Brian↗
NASA NTRS1978

Transport properties in the Jovian atmosphere

Transport properties in a Jupiter-like atmosphere (89 mol % hydrogen and 11 mol % helium) are obtained by using the method of the kinetic theory of gases. The transport collision integrals are calculated by fitting various two-body semiempirical interaction potentials for which the collision integrals are tabulated to ab initio quantum mechanical calculations of the two-body interactions. The collision integrals are used to calculate the binary diffusion coefficients, viscosity, and 'total' thermal conductivity of the pure gases and the gas mixtures at 1-atm pressure from 1000 K to 25,000 K.

Biolsi, L.↗
NASA NTRS2011

Kinetic Monte Carlo Investigation of the Effects of Vacancy Pairing on Oxygen Diffusivity in Yttria-Stabilized Zirconia

Yttria-stabilized zirconia s high oxygen diffusivity and corresponding high ionic conductivity, and its structural stability over a broad range of temperatures, have made the material of interest for use in a number of applications, for example, as solid electrolytes in fuel cells. At low concentrations, the stabilizing yttria also serves to increase the oxygen diffusivity through the presence of corresponding oxygen vacancies, needed to maintain charge neutrality. At higher yttria concentration, however, diffusivity is impeded by the larger number of relatively high energy migration barriers associated with yttrium cations. In addition, there is evidence that oxygen vacancies preferentially occupy nearest-neighbor sites around either dopant or Zr cations, further affecting vacancy diffusion. We present the results of ab initio calculations that indicate that it is energetically favorable for oxygen vacancies to occupy nearest-neighbor sites adjacent to Y ions, and that the presence of vacancies near either species of cation lowers the migration barriers. Kinetic Monte Carlo results from simulations incorporating this effect are presented and compared with results from simulations in which the effect is not present.

Good, Brian S.↗
NASA NTRSDate not supplied

Overview of NASA Research Activities in Shock Layer Kinetics and Radiation

Shock layer radiation is an important heating mechanism for entry probes to most planetary destinations and re-entry to Earth from beyond low Earth orbit. Our understanding of shock layer radiation phenomena has improved tremendously over the last decade thanks to NASA investment in fundamental radiation research through Entry Systems Modeling’s (ESM) Shock Layer Kinetics and Radiation (SLKR) task. This talk will overview the recent activities within SLKR, including validation of models through ground testing, flight instrumentation and remote observations; development of first principles ab initio calculations for reaction mechanisms and spectroscopic databases and advancing numerical and computational methods for prediction of radiation in reacting hypersonic flows. 1 Sr. Research Scientist, Aerothermodynamics Branch, and AIAA Associate Fellow.

Brett A. Cruden↗
NASA NTRS1997

Glycolonitrile oligomerization: structure of isolated oxazolines, potential heterocycles on the early earth

A study of glycolonitrile polymerization has led to the isolation and characterization of two 2,5-dihydro-4-aminooxazoles, 4 and 5. Previous reports have misassigned these structures as s-triazines or pyrimidines. X-ray diffraction analysis of crystals of 4 and an acetylated oxazole derivative of 5 (6) confirm the proposed structures. Ab initio computations are used to assess the relative thermodynamic stability of three trimer isomers (an s-triazine, an aminohydroxypyrimidine, and an aminooxazoline), and the results indicate that 4 is a novel kinetic product. Mechanistic considerations rationalize kinetic oxazole formation over the more customary triazine or pyrimidine trimers.

Non-NASA Center↗
NASA NTRS1994

Dissociative Recombination without a Curve Crossing

Ab initio calculations show that a curve crossing is not always needed for a high dissociative- recombination cross section. For HeH(+), in which no neutral states cross the ion potential curve, dissociative recombination is driven by the nuclear kinetic-energy operator on adiabatic potential curves. The kinetic-energy derivative operator allows for capture into repulsive curves that are outside of the classical turning points for the nuclear motion. The dominant dissociative route is the C (2)Sigma(+) state leading to H(n = 2) atoms. An analogous mechanism is proposed for the dissociative recombination of H3(+).

Guberman, Steven L.↗
NASA NTRS2015

Integrated Physics-based Modeling and Experiments for Improved Prediction of Combustion Dynamics in Low-Emission Systems

Concurrent simulation and experiment was undertaken to assess the ability of a hybrid RANS-LES model to predict combustion dynamics in a single-element lean direct-inject (LDI) combustor showing self-excited instabilities. High frequency pressure modes produced by Fourier and modal decomposition analysis were compared quantitatively, and trends with equivalence ratio and inlet temperature were compared qualitatively. High frequency OH PLIF and PIV measurements were also taken. Submodels for chemical kinetics and primary and secondary atomization were also tested against the measured behavior. For a point-wise comparison, the amplitudes matched within a factor of two. The dependence on equivalence ratio was matched. Preliminary results from simulation using an 18-reaction kinetics model indicated instability amplitudes closer to measurement. Analysis of the simulations suggested a band of modes around 1400 Hz were due to a vortex bubble breakdown and a band of modes around 6 kHz were due to a precessing vortex core hydrodynamic instability. The primary needs are directly coupled and validated ab initio models of the atomizer free surface flow and the primary atomization processes, and more detailed study of the coupling between the 3D swirling flow and the local thermoacoustics in the diverging venturi section.

combustion↗
NASA NTRS1997

An Accurate Potential Energy Surface for H2O

We have carried out extensive high quality ab initio electronic structure calculations of the ground state potential energy surface (PES) and dipole moment function (DMF) for H2O. A small adjustment is made to the PES to improve the agreement of line positions from theory and experiment. The theoretical line positions are obtained from variational ro-vibrational calculations using the exact kinetic energy operator. For the lines being fitted, the root-mean-square error was reduced from 6.9 to 0.08 /cm. We were then able to match 30,092 of the 30,117 lines from the HITRAN 96 data base to theoretical lines, and 80% of the line positions differed less than 0.1 /cm. About 3% of the line positions in the experimental data base appear to be incorrect. Theory predicts the existence of many additional weak lines with intensities above the cutoff used in the data base. To obtain results of similar accuracy for HDO, a mass dependent correction to the PH is introduced and is parameterized by simultaneously fitting line positions for HDO and D2O. The mass dependent PH has good predictive value for T2O and HTO. Nonadiabatic effects are not explicitly included. Line strengths for vibrational bands summed over rotational levels usually agree well between theory and experiment, but individual line strengths can differ greatly. A high temperature line list containing about 380 million lines has been generated using the present PES and DMF

Schwenke, David W.↗
NASA NTRS1984

Quantum chemical studies of a model for peptide bond formation. 3. Role of magnesium cation in formation of amide and water from ammonia and glycine

The SN2 reaction between glycine and ammonia molecules with magnesium cation Mg2+ as a catalyst has been studied as a model reaction for Mg(2+)-catalyzed peptide bond formation using the ab initio Hartree-Fock molecular orbital method. As in previous studies of the uncatalyzed and amine-catalyzed reactions between glycine and ammonia, two reaction mechanisms have been examined, i.e., a two-step and a concerted reaction. The stationary points of each reaction including intermediate and transition states have been identified and free energies calculated for all geometry-optimized reaction species to determine the thermodynamics and kinetics of each reaction. Substantial decreases in free energies of activation were found for both reaction mechanisms in the Mg(2+)-catalyzed amide bond formation compared with those in the uncatalyzed and amine-catalyzed amide bond formation. The catalytic effect of the Mg2+ cation is to stabilize both the transition states and intermediate, and it is attributed to the neutralization of the developing negative charge on the electrophile and formation of a conformationally flexible nonplanar five-membered chelate ring structure.

Magnesium/chemistry↗
NASA NTRS1977

The molecular composition of dense interstellar clouds

Presented in this paper is an ab initio chemical model for dense interstellar clouds that incorporates 598 grain surface reactions, with small grains providing the reaction area. Gas-phase molecules are depleted through collisions with grains. The abundances of 372 chemical species are calculated as a function of time and are found to be of sufficient magnitude to explain most observations. Peak abundances are achieved on time scales of the order of 100,000 to 1 million years, depending on cloud density and kinetic temperature. The reaction rates for ion-molecule chemistry are approximately the same, indicating that surface and gas-phase chemistry may be coupled in certain regions. The composition of grain mantles is shown to be a function of grain radius. In certain grain-size ranges, large molecules containing two or more heavy atoms are more predominant than lighter 'ices' - H2O, NH3, and CH4. It is possible that absorption due to these large molecules in the mantle may contribute to the observed 3-micron band in astronomical spectra.

Allen, M.↗
NASA NTRS2022

Recommended Direct Simulation Monte Carlo Collision Model Parameters for Reacting Methane Flows

While hydrocarbon combustion and gas-surface interactions have been simulated for a variety of applications in continuum solvers, such processes have been rarely investigated in Direct Simulation Monte Carlo (DSMC). Thus Variable Hard/Soft Sphere (VHS/VSS)collision parameters are not available for most molecules of interest in the decomposition of methane, a significant pyrolysis gas product and potential atmospheric species for Titan re-entries. As these properties are essential for accurate modeling of gas transport and thermochemical effects, a method must be devised to compute them for use in DSMC. In this study, the methodology outlined in Stephani et al is utilized to compute the VHS/VSS parameters by fitting the Ω(1,1)and Ω(2,2) collision integrals [1]. These collision integrals are typically computed from intermolecular potentials and are used in CFD simulations to calculate transport properties such as viscosity, thermal conductivity, and diffusion. This approach ensures consistency between micro-scale collisions and macro-scale transport properties. To produce the necessary collision integrals needed for the fitting routine, we utilize the method described by Laricchiuta et al, where collision integrals are computed based on the constituent particles’ polarizability and effective electron number [2].Preliminary computations of pure-species Ω(1,1)and Ω(2,2)collision integrals produced from the coupled approach of Stephani et al and Laricchiuta et al for methane show good agreement with those computed via available ab-initio Potential Energy Surface Computations, as seen in Figure 1. Above 300K, no more than 12% deviation is observed between the two results over the temperature range of the fitted PES data [3].Current results of this study include full collision-specificVSS transport properties for a complex-species mixture involving major byproducts of pure methane decomposition. Targeted final results will include a comprehensive list of VSS parameters relevant toDSMC simulations of Titan atmospheric entry, based on the reduced kinetic model of Savajano et al [4].

M R Gosma↗
NASA NTRS2002

Theoretical Study of Decomposition Pathways for HArF and HKrF

To provide theoretical insights into the stability and dynamics of the new rare gas compounds HArF and HKrF, reaction paths for decomposition processes HRgF to Rg + HF and HRgF to H + Rg + F (Rg = Ar, Kr) are calculated using ab initio electronic structure methods. The bending channels, HRgF to Rg + HF, are described by single-configurational MP2 and CCSD(T) electronic structure methods, while the linear decomposition paths, HRgF to H + Rg + F, require the use of multi-configurational wave functions that include dynamic correlation and are size extensive. HArF and HKrF molecules are found to be energetically stable with respect to atomic dissociation products (H + Rg + F) and separated by substantial energy barriers from Rg + HF products, which ensure their kinetic stability. The results are compatible with experimental data on these systems.

Chaban, Galina M.↗