SEARCH · Search NASA

Results for “ab initio modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗

NASA NTRSDate not supplied

Elucidation of Electrolyte Decomposition Pathways in Li--O2 Batteries using Ab Initio Computations

Li-O2 batteries offer exciting possibilities for energy-storage applications requiring higher energy content than current Li-ion battery technologies. However, development of commercial Li-O2 batteries faces many challenges including electrolyte decomposition in order to realize the full potential of this emerging technology. In this work, we present ab initio calculations that probe the role of reactive oxygen species including singlet oxygen, superoxide, and peroxides on electrolyte decomposition for potential organic electrolytes for use in Li-O2 batteries. We investigate the influence that the harsh electrochemical environment near the cathode plays in promoting electrolyte decomposition by using increasingly detailed computational models. Our results highlight existing problems with organic electrolytes for Li-O2 battery development and suggest potential design strategies for reducing electrolyte decomposition.

Li-O2 Batteries↗

NASA NTRS1996

Why HOC+ is detectable in interstellar clouds: the rate of the reaction between HOC+ and H2

The recent confirmation by Ziurys and Apponi of the detection of HOC+ toward Sgr B2 (OH), and their identification of the ion in Orion-KL and several other sources show that HOC+ is far more abundant than predicted by previous ion-molecule models. In these models, the reaction HOC(+) + H2 --> HCO(+) + H2 is assumed to rapidly destroy HOC+, based on the results of a prior calculation. We have recalculated the rate of this reaction as a function of temperature using a new ab initio potential surface and a phase space approach to the dynamics which includes tunneling. The newly calculated rate is small (< or = 1 x 10(-10) cm3 s-1) at temperatures under 100 K.

Non-NASA Center↗

NASA NTRS1988

A molecular orbital study of a model of the Mg2+ coordination complex of the self splicing reaction of ribosomal RNA

Recent discoveries have established the fact that RNA is capable of acting as an enzyme. In this study two different types of molecular orbital calculations, INDO and ab initio, were used in an attempt to assess the structural/functional role of the Mg2+ hydrated complex in ribozyme reactions. Preliminary studies indicate that the reaction is multistep and that the Mg2+ complex exerts a stabilizing effect on the intermediate or midpoint of the reaction.

Non-NASA Center↗

NASA NTRS2025

Line List for the A 2 Π−X 2 Σ + and the B 2 Σ + −X 2 Σ + Band Systems of CaF

The spectral analysis of molecule-rich asymptotic giant branch (AGB) stars is challenging. Although other calcium and fluorine bearing molecules are observed in the microwave and in the visible spectrum of AGB stars, CaF has never been detected, despite favorable chemical equilibrium predictions. Yet, measuring the CaF abundance could give more insight into the fluorine budget of AGB stars and allow better simulation of stellar spectra. In this work, we present an analysis of the visible spectrum of CaF obtained with a Fourier transform spectrometer. CaF A 2 Π−X 2 Σ + and the B 2 Σ + −X 2 Σ + band systems were excited with a hollow cathode discharge. Using previous fluorescence spectroscopy measurements, highly accurate ground state constants, and reasonable extrapolation schemes, the strongest features of CaF in the visible spectrum can be accurately modeled for both band systems. Spectroscopic constants are determined for A 2 Π with v ≤ 16 and for B 2 Σ + with v ≤ 20. Ab initio transition dipole moment curves of both transitions were calculated and scaled, and we provide a line list with Einstein A coefficients and oscillator strengths. This line list can be used to simulate spectra of CaF at temperatures and pressures relevant to astrophysical environments.

Léo Lavy↗

NASA NTRS2016

Ames S-32 O-16 O-18 Line List for High-Resolution Experimental IR Analysis

By comparing to the most recent experimental data and spectra of the SO2 628 ν1/ν3 bands (see Ulenikov et al., JQSRT 168 (2016) 29-39), this study illustrates the reliability and accuracy of the Ames-296K SO2 line list, which is accurate enough to facilitate such high-resolution spectroscopic analysis. The SO2 628 IR line list is computed on a recently improved potential energy surface (PES) refinement, denoted Ames-Pre2, and the published purely ab initio CCSD(T)/aug-cc-pVQZ dipole moment surface. Progress has been made in both energy level convergence and rovibrational quantum number assignments agreeing with laboratory analysis models. The accuracy of the computed 628 energy levels and line list is similar to what has been achieved and reported for SO2 626 and 646, i.e. 0.01-0.03 cm(exp −1) for bands up to 5500 cm(exp −1). During the comparison, we found some discrepancies in addition to overall good agreements. The three-IR-list based feature-by-feature analysis in a 0.25 cm(exp −1) spectral window clearly demonstrates the power of the current Ames line lists with new assignments, correction of some errors, and intensity contributions from varied sources including other isotopologues. We are inclined to attribute part of detected discrepancies to an incomplete experimental analysis and missing intensity in the model. With complete line position, intensity, and rovibrational quantum numbers determined at 296 K, spectroscopic analysis is significantly facilitated especially for a spectral range exhibiting such an unusually high density of lines. The computed 628 rovibrational levels and line list are accurate enough to provide alternatives for the missing bands or suspicious assignments, as well as helpful to identify these isotopologues in various celestial environments. The next step will be to revisit the SO2 828 and 646 spectral analyses.

Sulfur dioxide↗

NASA NTRS2016

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Comparisons are made between potential energy surfaces (PES) for N2 + N and N2 + N2 collisions and between rate coefficients for N2 dissociation that were computed using the quasiclassical trajectory method (QCT) on these PESs. For N2 + N we compare the Laganas empirical LEPS surface with one from NASA Ames Research Center based on ab initio quantum chemistry calculations. For N2 + N2 we compare two ab initio PESs (from NASA Ames and from the University of Minnesota). These use different methods for computing the ground state electronic energy for N4, but give similar results. Thermal N2 dissociation rate coefficients, for the 10,000K-30,000K temperature range, have been computed using each PES and the results are in excellent agreement. Quasi-stationary state (QSS) rate coefficients using both PESs have been computed at these temperatures using the Direct Molecular Simulation of Schwartentruber and coworkers. The QSS rate coefficients are up to a factor of 5 lower than the thermal ones and the thermal and QSS values bracket the results of shock-tube experiments. We conclude that the combination of ab initio quantum chemistry PESs and QCT calculations provides an attractive approach for the determination of accurate high-temperature rate coefficients for use in aerothermodynamics modeling.

quasiclassical trajectory method↗

NASA NTRS1997

The rate of the reaction between C2H and C2H2 at interstellar temperatures

The reaction between the radical C2H and the stable hydrocarbon C2H2 is one of the simplest neutral-neutral hydrocarbon reactions in chemical models of dense interstellar clouds and carbon-rich circumstellar shells. Although known to be rapid at temperatures > or = 300 K, the reaction has yet to be studied at lower temperatures. We present here ab initio calculations of the potential surface for this reaction and dynamical calculations to determine its rate at low temperature. Despite a small potential barrier in the exit channel, the calculated rate is large, showing that this reaction and, most probably, more complex analogs contribute to the formation of complex organic molecules in low-temperature sources.

Non-NASA Center↗

NASA NTRSDate not supplied

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗

NASA NTRS1987

Potential surfaces for O atom-polymer reactions

Ab initio quantum chemistry methods are used to study the energetics of interactions of O atoms with organic compounds. Polyethylene (CH2)n has been chosen as the model system to study the interactions of O(3P) and O(1D) atoms with polymers. In particular, H abstraction is investigated and polyethylene is represented by a C3 (propane) oligomeric model. The gradient method, as implemented in the GRADSCF package of programs, is used to determine the geometries and energies of products and reactants. The saddle point, barrier geometry is determined by minimizing the squares of the gradients of the potential with respect to the internal coordinates. To correctly describe the change in bonding during the reaction at least a two configuration MCSCF (multiconfiguration self consistent field) or GVB (generalized valence bond) wave function has to be used. Basis sets include standard Pople and Dunning sets, however, increased with polarization functions and diffuse p functions on both the C and O atoms. The latter is important due to the O(-) character of the wave function at the saddle point and products. Normal modes and vibrational energy levels are given for the reactants, saddle points and products. Finally, quantitative energetics are obtained by implementing a small CAS (complete active space) approach followed by limited configuration interaction (CI) calculations. Comparisons are made with available experimental data.

Laskowski, B. C.↗

NASA NTRS2004

Cross-Linked Nanotube Materials with Variable Stiffness Tethers

The constitutive properties of a cross-linked single-walled carbon nanotube material are predicted with a multi-scale model. The material is modeled as a transversely isotropic solid using concepts from equivalent-continuum modeling. The elastic constants are determined using molecular dynamics simulation. Some parameters of the molecular force field are determined specifically for the cross-linker from ab initio calculations. A demonstration of how the cross-linked nanotubes may affect the properties of a nanotube/polyimide composite is included using a micromechanical analysis.

Frankland, Sarah-Jane V.↗

NASA NTRS2020

Toward an In-Depth Material Model for Cermet Nuclear Thermal Rocket Fuel Elements

The development and qualification of nuclear thermal propulsion (NTP) fuel element technologies would be aided by an in-depth model of material response and failure modes at operating conditions. Integrated computational materials engineering techniques have the potential to provide such a model, as demonstrated here through three case studies focused on a tungsten-uranium mononitride cermet fuel. The first case focuses on the erosion of tungsten (W, also named wolfram), a nominal coating/cladding material, in hot hydrogen. Ab initio techniques are used to calculate erosion rates and thermal expansion at NTP operating conditions. The second focuses on the stability of uranium mononitride (UN) fuels at high temperature and in the presence of hydrogen. Phase diagram techniques reveal potential instabilities and decomposition pathways at high hydrogen concentrations. The third focuses on using microstructure information to predict high temperature mechanical response and failure of tungsten, used in refractory cermet materials. Combined finite element and discrete dislocation dynamics techniques provide mechanical properties in agreement with experimental methods. The integration of these techniques for an all-encompassing material model is discussed.

cermet↗

NASA NTRS1993

Statistical modeling of capture, association, and exit-channel dynamics in the CH3(+)/CH3CN system

The ion-molecule CH3(+) + CH3CN reaction is presently modeled by means of a master equation treatment incorporating weak collisions. The parameter required for the Rice-Ramsberger-Kassel-Markus (RRKM) treatment is derived from an ab initio investigation of the energy surface in question, and a means of including the capture rate coefficients in the RRKM approach is developed in which only the hindered dipole rotation is coupled in the reaction coordinate at large separations. Unimolecular rate coefficients for the (CH3CNCH3+)-activated complex are calculated for all product channels in the 300-600 K temperature range.

Smith, S. C.↗

NASA NTRS1997

Carbon in The Universe: PAHs and Clusters

Following the initial demonstration of this new technique (Science 265 1686 (1994)) and its application to a series of neutral PAHs which have been proposed as condidates for the UIRs (Nature 380, 227 (1996)), we have concentrated on two major aspects of this project. 1. Developing a detailed model for infrared emission spectra of a collection of highly excited PAH molecules, in which experimental bandshapes and temperature-dependent redshifts are used in conjunction with ab initio vibrational frequencies and intensities to simulate the UIR bands. This shows that a collection of nine different cations (as large as ovalene) reproduce the UIR features better than do a collection of the corresponding neutrals, but a detailed match with the UIRs is not obtained. 2. Construction of SPIRES apparatus for the study of PAH ion emission spectra. The design of this experiment is shown and described. Unfortunately a disasterous accident occurred just as we were preparing to start the testing of the ion apparatus. A vacuum implosion occurred, destroying the liquid He cooled monochromator. It has taken us nearly one full year to reconstruct this, and we arc only now in the final testing of the new system. We expect to try the ion experiments by the end of summer.

Saykally, Richard J.↗

NASA NTRS1997

Rotational Energy Transfer of N2 Gas Determined Using a New Ab Initio Potential Energy Surface

Rotational energy transfer between two N2 molecules is a fundamental process of some importance. Exchange is expected to play a role, but its importance is somewhat uncertain. Rotational energy transfer cross sections of N2 also have applications in many other fields including modeling of aerodynamic flows, laser operations, and linewidth analysis in nonintrusive laser diagnostics. A number of N2-N2 rigid rotor potential energy surface (PES) has been reported in the literature.

Huo, Winifred M.↗

NASA NTRS1991

Ab initio study of the electrostatic multipole nature of torsional potentials in CH3SSCH3, CH3SSH, and HOOH

The origin of torsional potentials in H3CSSCH3, H3CSSH, and HOOH and the anisotropy of the local charge distribution has been analyzed in terms of atomic multipoles calculated from the ab initio LCAO-MO-SCF wave function in the 6-31G* basis set. The results indicate that for longer -S-S-bonds the major contribution to these torsional barriers are electrostatic interactions of the atomic multipoles located on two atoms forming the rotated bond. This finding demonstrates the important role of electrostatic 1-2 interatomic interactions, usually neglected in conformational studies. It also opens the possibility to derive directly from accurate ab initio wave functions a simple nonempirical torsional potential involving atomic multipoles of two bonded atoms defining the torsional angle. For shorter -O-O- bonds, use of more precise models and inclusion of 1-3 interactions seems to be necessary.

NASA Program Exobiology↗

NASA NTRS2012

Multiscale Modeling of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2: Application to Lattice Thermal Conductivity

We are developing a multiscale framework in computational modeling for the ultra high temperature ceramics (UHTC) ZrB2 and HfB2. These materials are characterized by high melting point, good strength, and reasonable oxidation resistance. They are candidate materials for a number of applications in extreme environments including sharp leading edges of hypersonic aircraft. In particular, we used a combination of ab initio methods, atomistic simulations and continuum computations to obtain insights into fundamental properties of these materials. Ab initio methods were used to compute basic structural, mechanical and thermal properties. From these results, a database was constructed to fit a Tersoff style interatomic potential suitable for atomistic simulations. These potentials were used to evaluate the lattice thermal conductivity of single crystals and the thermal resistance of simple grain boundaries. Finite element method (FEM) computations using atomistic results as inputs were performed with meshes constructed on SEM images thereby modeling the realistic microstructure. These continuum computations showed the reduction in thermal conductivity due to the grain boundary network.

Lawson, John W.↗

NASA NTRS2012

Multiscale Modeling of UHTC: Thermal Conductivity

Lawson, John W.↗