SEARCH · Search NASA
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.
Recently, a semiempirical method for alloys based on equivalent crystal theory was introduced. The method successfully predicts the concentration dependence of the heat of formation and lattice parameter of binary alloys. A study of the parameters of the method is presented, along with new results for (gamma)Fe-Pd and (gamma)Fe-Ni alloys.
A computer program is described which (1) calculates thermodynamic functions (heat capacity, enthalpy, entropy, and free energy) for several optional forms of the partition function, (2) fits these functions to empirical equations by means of a least-squares fit, and (3) calculates, as a function of temperture, heats of formation and equilibrium constants. The program provides several methods for calculating ideal gas properties. For monatomic gases, three methods are given which differ in the technique used for truncating the partition function. For diatomic and polyatomic molecules, five methods are given which differ in the corrections to the rigid-rotator harmonic-oscillator approximation. A method for estimating thermodynamic functions for some species is also given.
Quantum mechanical computational methods are employed for an ab initio investigation of: (1) the molecular properties of the lowest isomers of the ClO dimer; and (2) predicted molecular and thermochemical properties. Techniques employed include electron correlation and particularly singles and doubles coupled-cluster (CCSD) theory with or without perturbational estimates of the effects of connected triple excitations. The isomers ClOClO and ClClO2 are found to have higher energies than the ClOOCl isomer, and the theoretical vibrational frequencies of the isomers are well correlated with experimental data. Experimental values of the heat of formation for the isomers are also compared with calculations based on an isodesmic reaction with Cl2O, H2O, and HOOH.
The structure of ClO2 has been calculated for the X 2A-double prime ground state using unrestricted Hartree-Fock (UHF), unrestricted second-order Moller-Plesset perturbation (UMP2), configuration interaction employing single and double excitation (CISD), and quadratic configuration interaction (QCI) ab initio molecular orbital methods. Calculations using UMP2 and CISD wave functions predict a ClO bond length of 1.728 +/- 0.01 A. The single-configuration-based QCI in the singles and doubles space with perturbation inclusion of triple substitutions, denoted QCISD(T), yield a ClO bond length of 2.205 A. The QCI results are consistent with results of Jensen (1990) who showed that the ClO bond length is 2.181 A using annihilated self-consistent methods (AUMP2). The thermochemistry of ClO2 radical has been calculated using MP2 and QCI methods employing an isodesmic scheme. Our scheme predicts the heat of formation for ClO2 at 0 K to be 24.6 +/- 2 kcal/mol.
The potential energy surface for NH2+O has been characterized using complete active space self-consistent field (CASSCF)/derivative calculations to determine stationary point geometries and frequencies followed by internally contracted configuration interaction (ICCI) calculations to determine the energetics. The calculations predict a NO bond strength of 85.8 kcal/mol for NH2O. The barrier for isomerization of NH2O to trans-HNOH is predicted to be 48.0 kcal/mol and the barriers for H+HNO forming NH2O and NHOH are predicted to be 2.1 and 8.3 kcal/mol, respectively (all corrected for zero-point energy). The computed heats of formation for NH2O and cis- and trans-HNOH are in good agreement with the present results. The barrier for H + HNO yields H2 + NO is computed to be about 0.3 kcal/mol.
Libraries of thermodynamic data and transport properties are given for individual species in the form of least-squares coefficients. Values of C(sup 0)(sub p)(T), H(sup 0)(T), and S(sup 0)(T) are available for 1130 solid, liquid, and gaseous species. Viscosity and thermal conductivity data are given for 155 gases. The original C(sup 0)(sub p)(T) values were fit to a fourth-order polynomial with integration constants for H(sup 0)(T) and S(sup 0)(T). For each species the integration constant for H(sup 0)(T) includes the heat of formation. Transport properties have a different functional form. The temperature range for most of the data is 300 to 5000 K, although some of the newer thermodynamic data have a range of 200 to 6000 K. Because the species are mainly possible products of reaction, the data are useful for chemical equilibrium and kinetics computer codes. Much of the data has been distributed for several years with the NASA Lewis equilibrium program CET89. The thermodynamic properties of the reference elements were updated along with about 175 species that involve the elements carbon, hydrogen, oxygen, and nitrogen. These sets of data will be distributed with the NASA Lewis personal computer program for calculating chemical equilibria, CETPC.
The Bozzolo, Ferrante and Smith (BFS) method for alloys is applied to the study of the defect structure of B2 FeAI alloys. First-principles Linear Muffin Tin Orbital calculations are used to determine the input parameters to the BFS method used in this work. The calculations successfully determine the phase field of the B2 structure, as well as the dependence with composition of the lattice parameter. Finally, the method is used to perform 'static' simulations where instead of determining the ground state configuration of the alloy with a certain concentration of vacancies, a large number of candidate ordered structures are studied and compared, in order to determine not only the lowest energy configurations but other possible metastable states as well. The results provide a description of the defect structure consistent with available experimental data. The simplicity of the BFS method also allows for a simple explanation of some of the essential features found in the concentration dependence of the heat of formation, lattice parameter and the defect structure.
The BFS method for alloys is applied to the study of 200 alloys obtained from adding Cu and Au impurities to a Ni3Al matrix. We analyze the trends in the bulk properties of these alloys (heat of formation, lattice parameter, and bulk modulus) and detect specific alloy compositions for which these quantities have particular values. A detailed analysis of the atomic interactions that lead to the preferred ordering patterns is presented.
We report calculations of the minimum energy pathways connecting CH2 + N2 to diazomethane and diazirine, for the rearrangement of diazirine to diazomethane, for the dissociation of diazirine to HCN2+H, and of diazomethane to CH2N+N. The calculations use Complete Active Space Self-Consistent Field (CASSCF) derivative methods to characterize the stationary points and Internally Contracted Configuration Interaction (ICCI) to determine the energetics. The calculations suggest a potential new source of prompt NO from the reaction CH2 with N2 to give diazirine, and subsequent reaction of diazirine with hydrogen abstracters to form doublet HCN2, which leads to HCN+N(S-4) on the previously studied CH+N2 surface. The calculations also predict accurate 0 K heats of formation of 77.7 kcal/mol and 68.0 kcal/mol for diazirine and diazomethane, respectively.
We report calculations of the minimum energy pathways connecting (1)CH2+N2 to diazomethane and diazirine, for the rearrangement of diazirine to diazomethane, for the dissociation of diazirine to HCN2+H, and of diazomethane to CH2N+N. The calculations use complete active space self-consistent field (CASSCF) derivative methods to characterize the stationary points and internally contracted configuration interaction (ICCI) to determine the energetics. The calculations suggest a potential new source of prompt NO from the reaction of (1)CH2 with N2 to give diazirine, and subsequent reaction of diazirine with hydrogen abstracters to form doublet HCN2, which leads to HCN+N(S-4) on the previously studied CH+N2 Surface. The calculations also predict accurate 0 K heats of formation of 77.7 kcal/mol and 68.0 kcal/mol for diazirine and diazomethane, respectively.
A one-formula surrogate fuel formulation and its quasi-global combustion kinetics model are developed to support the design of injectors and thrust chambers of kerosene-fueled rocket engines. This surrogate fuel model depicts a fuel blend that properly represents the general physical and chemical properties of kerosene. The accompanying gaseous-phase thermodynamics of the surrogate fuel is anchored with the heat of formation of kerosene and verified by comparing a series of one-dimensional rocket thrust chamber calculations. The quasi-global combustion kinetics model consists of several global steps for parent fuel decomposition, soot formation, and soot oxidation, and a detailed wet-CO mechanism. The final thermophysics formulations are incorporated with a computational fluid dynamics model for prediction of the combustor efficiency of an uni-element, tri-propellant combustor and the radiation of a kerosene-fueled thruster plume. The model predictions agreed reasonably well with those of the tests.
Due to formal and computational advances in coupled-cluster theory over the past few years, it is now possible to obtain very accurate molecular geometries, vibrational frequencies, heats of formation, binding energies, and vertical electronic excitation energies. The present monograph reviews the advances that have lead to the current state-of-the art, and also summarizes selected examples from the published literature.
The first spectra of the A (2)Pi(sub 3/2) from X (2)Pi(sub 3/2) electronic transition of BrO using Fourier transform ultraviolet spectroscopy are obtained. Broadband vibrational spectra acquired at 298 +/- 2 K and 228 +/- 5 K, as well as high-resolution rotational spectra of the A from X 7,0 and 12,0 vibrational bands are presented. Wavenumber positions for the spectra are obtained with high accuracy, and cross section assignments are made, incorporating the existing literature. With 35 cm(exp -1) (0.40 nm) resolution the absolute cross section at the peak of the 7,0 band is determined to be (1.58 +/- 0.12) x 10(exp -17) sq cm/molecule at 298 +/- 2 K and (1.97 +/- 0.15) x 10(exp -17) sq cm/molecule at 228 +/- 5 K. BrO dissociation energies are determined with a graphical Birge-Sponer technique, using Le Roy-Bernstein theory to place an upper limit on the extrapolation. From the ground-state dissociation energy, D(sub o)" = 231.0 +/- 1.7 kJ/mol, the heat of formation of BrO(g) is calculated, del(sub f)H(0 K) = 133.7 +/- 1.7 kJ/mol and del(sub f)H(298.15 K) = 126.2 +/- 1.7 kJ/mol. Cross sections for the high-resolution 7,0 and 12,0 rotational peaks are the first to be reported. The band structures are modeled, and improved band origins, rotational constants, centrifugal distortion constants, and linewidths are determined. In particular, J-dependent linewidths and lifetimes are observed for the both the 7,0 and 12,0 bands.
A one-formula surrogate fuel formulation and its quasi-global combustion kinetics model are developed to support the design of injectors and thrust chambers of kerosene-fueled rocket engines. This surrogate fuel model depicts a fuel blend that properly represents the general physical and chemical properties of kerosene. The accompanying gaseous-phase thermodynamics of the surrogate fuel is anchored with the heat of formation of kerosene and verified by comparing a series of one-dimensional rocket thrust chamber calculations. The quasi-global combustion kinetics model consists of several global steps for parent fuel decomposition, soot formation, and soot oxidation and a detailed wet-CO mechanism to complete the combustion process. The final thermophysics formulations are incorporated with a computational fluid dynamics model for prediction of the combustion efficiency of an unielement, tripropellant combustor and the radiation of a kerosene-fueled thruster plume. The model predictions agreed reasonably well with those of the tests.
The ability of modern state-of-the art ab initio quantum chemical techniques to characterize reliably the gas-phase molecular structure, vibrational spectrum, electronic spectrum, and thermal stability of chlorine oxide and nitrogen oxide species will be demonstrated by presentation of some example studies. In particular the geometrical structures, vibrational spectra, and heats of formation Of ClNO2, CisClONO, and trans-ClONO are shown to be in excellent agreement with the available experimental data, and where the experimental data are either not known or are inconclusive, the ab initio results are shown to fill in the gaps and to resolve the experimental controversy. In addition, ab initio studies in which the electronic spectra and the characterization of excited electronic states of ClONO2, HONO2, ClOOC17 ClOOH, and HOOH will also be presented. Again where available, the ab initio results are compared to experimental observations, and are used to aid in the interpretation of the experimental studies.
In this work we demonstrate that density functional theory (DFT) methods make an important contribution to understanding chemical systems and are an important additional method for the computational chemist. We report calibration calculations obtained with different functionals for the 55 G2 molecules to justify our selection of the B3LYP functional. We show that accurate geometries and vibrational frequencies obtained at the B3LYP level can be combined with traditional methods to simplify the calculation of accurate heats of formation. We illustrate the application of the B3LYP approach to a variety of chemical problems from the vibrational frequencies of polycyclic aromatic hydrocarbons to transition metal systems. We show that the B3LYP method typically performs better than the MP2 method at a significantly lower computational cost. Thus the B3LYP method allows us to extend our studies to much larger systems while maintaining a high degree of accuracy. We show that for transition metal systems, the B3LYP bond energies are typically of sufficient accuracy that they can be used to explain experimental trends and even differentiate between different experimental values. We show that for boron clusters the B3LYP energetics are not as good as for many of the other systems presented, but even in this case the B3LYP approach is able to help understand the experimental trends.
InXn atomization energies are computed for n = 1-3 and X = H, Cl, and CH3. The geometries and frequencies are determined using density functional theory. The atomization energies are computed at the coupled cluster level of theory. The complete basis set limit is obtained by extrapolation. The scalar relativistic effect is computed using the Douglas-Kroll approach. While the heats of formation for InH, InCl and InCl3 are in good agreement with experiment, the current results show that the experimental value for In(CH3)3 must be wrong.
The chondrule cooling rates used in most chondrule-formation models appear to be too low. Recent petrographic evidence indicates that the amount of crystal (especially olivine) growth that occurred after the last melting event was about 30 smaller than the grain sizes simulated in order to estimate cooling rates. The smaller amount of growth leads to an upwards revision of cooling rates by about a factor of 1000. Most chondrules are porphyritic. They consist of large and small crystals of olivine and, less commonly, pyroxene immersed in a mesostasis having a plagioclase-rich composition. In the most primitive chondrites the mesostasis is often vitreous. Because the large majority of chondrules contain FeS, it is clear that the nebula had cooled below the FeS condensation temperature (ca. 650 K) before chondrule formation occurred. The high FeO/(FeO+MgO) ratios of some chondrules require still lower nebular temperatures (less than 500 K). The traditional view has been that porphyritic chondrules formed in a single heating/cooling event and many laboratory experiments have been carried out in various kinds of kinds of furnaces to try to simulate the formation of chondrules textures in a single heating/cooling cycle. These furnace experiments have been used to infer the cooling rates of chondrules during the temperature range at which olivine crystallized from the melt. Most of these inferred values are in the range 0.01-1 K per second. These low cooling rates are problematical because there is no long-term nebular environment that yields such values. In transparent regions chondrules would cool at rates orders of magnitude higher, whereas in an opaque nebular disk the cooling rates would be many orders of magnitude lower. And these latter conditions are not suitable locations for chondrule formation because such high temperatures would cause the complete evaporation of chondrules (which have melting temperatures about 600 K higher than their evaporation temperatures). During the last several years several kinds of petrographic evidence indicating rapid chondrule cooling have been recognized. These include thin overgrowths on relict grains, clusters of small crystals that grew following the most recent melting event, and fragments preserving shardlike shapes that would have been hidden had several tens of micrometers of new growth occurred following the last melting event. And there are other indications of rapid cooling including the preservation of volatiles such as FeS and one case where the modeling of O-isotopic and FeO/(FeO+MgO) gradients indicated a high cooling rate.