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.
The TChem open-source software is a toolkit for computing thermodynamic properties, source term, and source term’s Jacobian matrix for chemical kinetic models that involve gas and surface reactions.
To date, approximately 120 planet-sized objects have been discovered around other stars, mostly through the radial-velocity technique. This technique can provide information about a planet s minimum mass and its orbital period and distance; however, few other planetary data can be obtained at this point in time unless we are fortunate enough to find an extrasolar giant planet that transits its parent star (i.e., the orbit is edge-on as seen from Earth). In that situation, many physical properties of the planet and its parent star can be determined, including some compositional information. Our prospects of directly obtaining spectra from extrasolar planets may improve in the near future, through missions like NASA's Terrestrial Planet Finder. Most of the extrasolar giant planets (EGPs) discovered so far have masses equal to or greater than Jupiter's mass, and roughly 16% have orbital radii less than 0.1 AU - extremely close to the parent star by our own Solar-System standards (note that Mercury is located at a mean distance of 0.39 AU and Jupiter at 5.2 AU from the Sun). Although all EGPs are expected to have hydrogen-dominated atmospheres similar to Jupiter, the orbital distance can strongly affect the planet's temperature, physical, chemical, and spectral properties, and the abundance of minor, detectable atmospheric constituents. Thermochemical equilibrium models can provide good zero-order predictions for the atmospheric composition of EGPs. However, both the composition and spectral properties will depend in large part on disequilibrium processes like photochemistry, chemical kinetics, atmospheric transport, and haze formation. We have developed a photochemical kinetics, radiative transfer, and 1-D vertical transport model to study the atmospheric composition of EGPs. The chemical reaction list contains H-, C-, O-, and N-bearing species and is designed to be valid for atmospheric temperatures ranging from 100-3000 K and pressures up to 50 bar. Here we examine the effect of stellar distance (e.g., incident ultraviolet flux, atmospheric temperature) on the chemical properties of EGPs. The model is applied to two generic Class II and III intermediate temperature EGPs located at 3.3 and 0.27 AU from a solar-like parent star, and the results are compared with a model for Jupiter at 5.2 AU.
A CFD analysis of a low pressure nuclear rocket concept is presented with the use of an advanced chemical kinetics, Navier-Stokes code. The computations describe the flow field in detail,including gas dynamic, thermodynamic and chemical properties, as well as global performance quantities such as specific impulse. Computational studies of several rocket nozzle shapes are conducted in an attempt to maximize hydrogen recombination. These Navier-Stokes calculations, which include real gas and viscous effects, predict lower performance values than have been reported heretofore.
A CFD analysis of a low pressure nuclear rocket concept is presented with the use of an advanced chemical kinetics, Navier-Stokes code. The computations describe the flow field in detail, including gas dynamic, thermodynamic and chemical properties, as well as global performance quantities such as specific impulse. Computational studies of several rocket nozzle shapes are conducted in an attempt to maximize hydrogen recombination. These Navier-Stokes calculations, which include real gas and viscous effects, predict lower performance values than have been reported heretofore.
In chemical kinetics research, kinetic models containing hundreds of species and tens of thousands of elementary reactions are commonly used to understand and predict the behavior of reactive chemical systems. Reaction Mechanism Generator (RMG) is a software suite developed to automatically generate such models by incorporating and extrapolating from a database of known thermochemical and kinetic parameters. Here, we present the recent version 3 release of RMG and highlight improvements since the previously published description of RMG v1.0. Most notably, RMG can now generate heterogeneous catalysis models in addition to the previously available gas- and liquid-phase capabilities. For model analysis, new methods for local and global uncertainty analysis have been implemented to supplement first-order sensitivity analysis. The RMG database of thermochemical and kinetic parameters has been significantly expanded to cover more types of chemistry. The present release includes parallelization for faster model generation and a new molecule isomorphism approach to improve computational performance. RMG has also been updated to use Python 3, ensuring compatibility with the latest cheminformatics and machine learning packages. Overall, RMG v3.0 includes many changes which improve the accuracy of the generated chemical mechanisms and allow for exploration of a wider range of chemical systems.
A simplified chemical kinetic reaction mechanism for the combustion of a hydrocarbon fuel is presented and discussed. The observed kinetic behavior of propane combustion as determined in shock tube studies and the calculated kinetic behavior of propane/air mixtures as derived from a detailed propane combustion mechanism were used as a data base for constructing and refining the mechanism. The mechanism contains 13 chemical species and 26 chemical reactions. Numerical studies revealed that the simplified mechanism can reproduce the observed kinetic behavior of propane/air mixtures, including laminar flame speeds, over a wide range of temperatures and pressures. It is anticipated that the mechanism can be used in various combustion applications and provide a more realistic description of the combustion process than can be achieved using global reaction mechanisms.
The role of chemical kinetics in the formation of nitric oxide during the combustion of methane was examined analytically by means of a detailed chemical mechanism for the oxidation of methane, for the reaction between hydrocarbon fragments, and for the formation of nitric oxide. By comparing predicted nitric oxide levels with values reported in the literature from jet-stirred combuster experiments, it was determined that the nitric oxide levels observed in fuel-rich flames cannot be described by a mechanism in which the rate of nitric oxide formation is controlled solely by the kinetics of oxygen atom formation. A proposed mechanism for the formation of nitric oxide in methane-rich flames reproduces the observed levels. The oxidation of hydrogen cyanide appears to be an important factor in nitric oxide formation.
In a study of methane ignition sensitized by NO2, the induction period chemical kinetics were examined by monitoring reactants and intermediates (CH4, NO2, OH, NO, CO, CO2, and H2O) prior to ignition. Data were obtained at pressure of 2-4 atm, temperatures of 1300-1900 K, equivalence ratios of 0.5-1.0, and NO2 up to 3.4 mole %. A tentative NO2 sensitization mechanism was deduced from the analysis of photometric observations and the comparison of computer modeling with observations.
This computational singular perturbation library (CSPlib) software is useful for analysis of chemical kinetic systems. The code is used to analyze and describe the dynamical landscape underlying the time evolution of a chemical system. Given a database of state vectors, it allows the analysis of the source term and Jacobian matrix of the chemical system, the identification and decoupling of fast and slow processes, and the analysis of cause-and-effect relationships in the coupled system dynamics. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. SAND2021-2391 O
Approximate thermochemical tables are presented for some C-H and C-H-O species and for some ionized species, supplementing the JANAF Thermochemical Tables for application to finite-chemical-kinetics calculations. The approximate tables were prepared by interpolation and extrapolation of limited available data, especially by interpolations over chemical families of species. Original estimations have been smoothed by use of a modification for the CDC-6600 computer of the Lewis Research Center PACl Program which was originally prepared for the IBM-7094 computer Summary graphs for various families show reasonably consistent curvefit values, anchored by properties of existing species in the JANAF tables.
Work on understanding the mechanisms of mixing and reaction in high speed turbulent reacting flows was continued. Efforts, in particular, were concentrated on taking advantage of modern computational methods to simulate high speed turbulent flows. In doing so, two methodologies were used: large eddy simulations (LES) and direct numerical simulations (DNS). In the work related with LES the objective is to study the behavior of the probability density functions (pdfs) of scalar properties within the subgrid in reacting turbulent flows. The data base obtained by DNS for a detailed study of the pdf characteristics within the subgrid was used. Simulations are performed for flows under various initializations to include the effects of compressibility on mixing and chemical reactions. In the work related with DNS, a two-dimensional temporally developing high speed mixing layer under the influence of a second-order non-equilibrium chemical reaction of the type A + B yields products + heat was considered. Simulations were performed with different magnitudes of the convective Mach numbers and with different chemical kinetic parameters for the purpose of examining the isolated effects of the compressibility and the heat released by the chemical reactions on the structure of the layer. A full compressible code was developed and utilized, so that the coupling between mixing and chemical reactions is captured in a realistic manner.
Soot formation rates measured in fuel-rich premixed flat flames are frequently used to calibrate or validate chemical kinetic models of soot formation. Unfortunately, these flames feature an extended region of soot precursor particle inception and carbonization that complicates interpretation of soot measurements and leads to a fundamental inconsistency in the nature of the soot material that is modeled versus what is being measured when using non-intrusive, optical techniques. In the work presented here, previously reported data on two canonical sooting ethylene-air premixed flames at 1 atm pressure are interpreted via a new analysis approach that combines soot optical dispersion coefficient measurements with soot carbonization kinetics. This analytical approach explicitly accounts for the production of poorly ordered soot precursor particle mass and its carbonization over time in the flames, providing a clear distinction between the formation rate of precursor particles and their transformation into ordered, solid soot particulate mass. In particular, the results of the analysis show that the precursor particles form much earlier in the flame than the majority of the carbonized soot and their formation rate is two to three times faster than that of ordered soot. The results also show that particle agglomeration begins when the particles are at an intermediate state of carbonization. In conclusion, these results offer a valuable new interpretation of these important datasets and should lead to substantial improvements in the development and calibration of quantitative soot models.
Transient response of thermal-protection material to heat applied to surface calculated using CHAP III computer program. Used to analyze chemical kinetics of pyrolysis gas in detail and examine pyrolysis reactions-in-depth. Analysis includes deposition of solid products produced by chemical reactions in gas phase. Uses modeling technique that approximates wide range of ablation problems. CHAP III written in FORTRAN IV.
The growing interest in using ammonia as a carbon-neutral replacement fuel has prompted numerous recent research efforts towards understanding fundamental combustion characteristics and chemical kinetics of ammonia oxidation. Despite these recent efforts, there remains a large deficiency in extinction limit data for ammonia flames, particularly under elevated pressure and temperature conditions relevant to gas turbines. Even the literature reaction mechanisms developed specifically for ammonia have been untested under gas turbine relevant conditions, and existing kinetic models have shown large variations in laminar flame speed, extinction strain rate, and speciation predictions for NH3–H2 mixtures at ambient conditions. Recent work by Thomas and co-workers, for example, compared different ammonia kinetic models in the literature for predicting their extinction strain rate measurements for non-premixed NH3-H2 counterflow flames, and again have highlighted the discrepancies in model predictions. Here, our work has undertaken to further highlight important reaction pathways for ammonia oxidation that may be responsible for these discrepancies for extinction strain rate predictions, as well as elucidate sensitive reaction steps with large uncertainties that require further attention. This work points to a need for an improved understanding of ammonia combustion chemistry as well as new datasets at elevated pressure and temperature conditions to better constrain rate expressions for key reaction steps.
Computational fluid dynamics (CFD), chemical kinetics, and aerosol dynamics were combined to model the formation of uranyl fluoride particles from the hydrolysis of uranium hexafluoride gas. This chemical process has been studied for many decades, however, detailed measurements of aerosol formation have only become available in the past few years which provide a basis for model comparison. Here, CFD simulations predicted complicated flow patterns in the impinging jet gas reactor. Aerosol formation simulations also predicted higher mass concentrations than were observed experimentally. This suggested that experimental data provided a partial representation which was subsequently enhanced by modeling and simulation.
Some of the most interesting chemistry in the Solar System involves changes in the oxidation state of the simple carbon species. The chemical pathways for the conversion of CH4 to CO and CO2 are for the most part known. The reverse process, the reduction of CO to CH4, is, however, poorly understood. This is surprising in view of the importance of the reduction process in the chemistry of the Solar System. Recently we investigated the chemical kinetics of a hitherto unsuspected reaction. It is argued that the formation of the methoxy radical (CH3O) from H+H2CO may play an essential role in the reduction of CO to CH4. The rate coefficient for this reaction has been estimated using the approximate theory of J. Troe and transition state theory. We will discuss the implications of this reaction for the chemistry of CO on Jupiter, in the solar nebula, for interpreting the laboratory experiments of A. Bar-Nun and A. Shaviv and A. Bar-Nun and S. Chang, and for organic synthesis in the prebiotic terrestrial atmosphere. The possible relation of CO reduction in the solar nebula and polyoxymethylene observed in comet Halley will be discussed.
Four years ago this month, a hitherto unknown comet in loose orbit around Jupiter passed so near the giant planet that it was torn apart into 20 fragments by tides. One orbit later, two years ago this month, the fragments of doomed comet P/Shoemaker Levy (SL9) fell into Jupiter. The enormous energies of these impacts (the largest fragments were nearly 1 km across and, hitting at 60 km/s, released some 2-4 x 10(exp 27) ergs) produced enormous explosions. Several of the ejecta plumes were imaged towering 3000 km above Jupiter's limb. The heat released when the plumes fell was considerable and easily observed on Earth. The impacts produced strong shocks, both promptly at the impact site and again, later, and over thousands of kilometers, when the ejecta plume reentered the atmosphere. The focus of this talk will be to discuss what the SL9 impacts taught us about impact shock chemistry - the processes, the ingredients, the results - and what inferences we may draw for impacts on early Earth. Shock chemistry generates a suite of molecules not usually seen on Jupiter. The most surprising report was of a huge amount of diatomic sulfur S2 at the site of the G impact. Other unusual products include CS, CS2, OCS, H2S, SO2, HCN, CO, and H2O; although H2S and H2O are doubtless abundant below the visible clouds. Hot or enhanced CH4 and NH3 were also detected. A general rule of shock chemistry is that CO forms until either C or O is exhausted. If O greater than C, the other products are oxidized, and excess O goes to H2O. If C greater than O, the other products are reduced, and excess C goes to HCN, C2H2, and a wide variety of more complicated organics. Ultimately, given time, the carbon would react all the way to graphite, but in practice the reactions are incomplete. The dark ejecta debris were probably composed in part of carbonaceous particles generated by the shocks. In a sense, the SL9 impacts performed the famous Miller-Urey experiment on a grand scale, with one result being the production of a lot of complex brown organic solids (called "tholins"). We use, a straightforward chemical kinetics model for the H, N, C, O, S system to follow the nonequilibrium chemistry behind the shocks. The model traces the evolving chemical composition of a parcel of gas by directly integrating the web of chemical reactions. Pressure and temperature histories of the parcels are patterned after those calculated by numerical hydrodynamic simulations of the ejecta plume. A given plume parcel is generally shocked twice; t.e a parcel shocked near the impact site is ejected at high velocity and is shocked again when it reenters the atmosphere. The final state of the gas depends mostly on the second shock, provided that the latter is hot enough. The chemical evidence is ambiguous, but most indications are that C greater than O in the shocked, reacting gas. Telltale signatures of abundant oxygen - SO2, SO, CO2, O2 - were not seen, while signatures of abundant carbon - CS, CS2, and HCN - were. On the other hand, abundant H2O would appear to require O greater than C, and two other observed sulfur species, S2 and OCS, appear to form more easily in a somewhat oxidized gas, presumable vaporized from the comet itself. Since on general principles one expects the -comet to have had a more-or-less cosmic composition, i.e. O greater than C, the production of CS, CS2, and HCN probably requires C greater than O in the shocked jovian air. This in turn implies that even the largest fragments released the bulk of their energy above the jovian water table, in all likelihood above 5 bars . There is no evidence in favor of the proposition that a significant amount of wet jovian air was shocked strongly enough to coax water to react; i.e. wet jovian air saw only temperatures significantly below 2000 K.
The recent focus of this research project has been to model the combustion of isolated metal droplets and, in particular, to couple the existing theories and formulations of phenomena such as condensation, reaction kinetics, radiation, and surface reactions to formulate a more complete combustion model. A fully transient, one-dimensional (spherical symmetry) numerical model that uses detailed chemical kinetics, multi-component molecular transport mechanisms, condensation kinetics, and gas phase radiation heat transfer was developed. A coagulation model was used to simulate the particulate formation of MgO. The model was used to simulate the combustion of an Mg droplet in pure O2 and CO2. Methanol droplet combustion is considered as a test case for the solution method for both quasi-steady and fully transient simulations. Although some important processes unique to methanol combustion, such as water absorption at the surface, are not included in the model, the results are in sufficient agreement with the published data. Since the major part of the heat released in combustion of Mg, and in combustion of metals in general, is due to the condensation of the metal oxide, it is very important to capture the condensation processes correctly. Using the modified nucleation theory, an Arrhenius type rate expression is derived to calculate the condensation rate of MgO. This expression can be easily included in the CHEMKIN reaction mechanism format. Although very little property data is available for MgO, the condensation rate expression derived using the existing data is able to capture the condensation of MgO. An appropriate choice of the reference temperature to calculate the rate coefficients allows the model to correctly predict the subsequent heat release and hence the flame temperature.