Search NASA⌕ Search

SEARCH · Search NASA

Results for “reaction rate equations”

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.

At least 19 records

Effect of Wind Velocity on Flame Spread in Microgravity

A three-dimensional, time-dependent model is developed describing ignition and subsequent transition to flame spread over a thermally thin cellulosic sheet heated by external radiation in a microgravity environment. A low Mach number approximation to the Navier Stokes equations with global reaction rate equations describing combustion in the gas phase and the condensed phase is numerically solved. The effects of a slow external wind (1-20 cm/s) on flame transition are studied in an atmosphere of 35% oxygen concentration. The ignition is initiated at the center part of the sample by generating a line-shape flame along the width of the sample. The calculated results are compared with data obtained in the 10s drop tower. Numerical results exhibit flame quenching at a wind speed of 1.0 cm/s, two localized flames propagating upstream along the sample edges at 1.5 cm/s, a single line-shape flame front at 5.0 cm/s, three flames structure observed at 10.0 cm/s (consisting of a single line-shape flame propagating upstream and two localized flames propagating downstream along sample edges) and followed by two line-shape flames (one propagating upstream and another propagating downstream) at 20.0 cm/s. These observations qualitatively compare with experimental data. Three-dimensional visualization of the observed flame complex, fuel concentration contours, oxygen and reaction rate isosurfaces, convective and diffusive mass flux are used to obtain a detailed understanding of the controlling mechanism, Physical arguments based on lateral diffusive flux of oxygen, fuel depletion, oxygen shadow of the flame and heat release rate are constructed to explain the various observed flame shapes.

Prasad, Kuldeep↗

Development of a second order closure model for computation of turbulent diffusion flames

A typical eddy box model for the second-order closure of turbulent, multispecies, reacting flows developed. The model structure was quite general and was valid for an arbitrary number of species. For the case of a reaction involving three species, the nine model parameters were determined from equations for nine independent first- and second-order correlations. The model enabled calculation of any higher-order correlation involving mass fractions, temperatures, and reaction rates in terms of first- and second-order correlations. Model predictions for the reaction rate were in very good agreement with exact solutions of the reaction rate equations for a number of assumed flow distributions.

Varma, A. K.↗

Neutron radiation characteristics of plutonium dioxide fuel

The major sources of neutrons from plutonium dioxide nuclear fuel are considered in detail. These sources include spontaneous fission of several of the Pu isotopes, (alpha, n) reactions with low Z impurities in the fuel, and (alpha, n) reactions with O-18. For spontaneous fission neutrons a value of (1.95 + or - 0.07) X 1,000 n/s/g PuO2 is obtained. The neutron yield from (alpha, n) reactions with oxygen is calculated by integrating the reaction rate equation over all alpha-particle energies and all center-of-mass angles. The results indicate a neutron emission rate of (1.14 + or - 0.26) X 10,000 n/s/g PuO2. The neutron yield from (alpha, n) reactions with low Z impurities in the fuel is presented in tabular form for one part part per million of each impurity. The total neutron yield due to the combined effects of all the impurities depends upon the fractional weight concentration of each impurity. The total neutron flux emitted from a particular fuel geometry is estimated by adding the neutron yield due to the induced fission to the other neutron sources.

Taherzadeh, M.↗

Neutron Radiation Characteristics of Plutonium Dioxide Fuel

The major sources of neutrons from plutonium dioxide nuclear fuel are considered in detail. These sources include spontaneous fission of several of the Pu isotopes, reactions with low Z impurities in the fuel, and reactions with O-18. For spontaneous fission neutrons a value of (1.95 plus or minus 0.07) X 1,000 n/s/q PuO2 is obtained. The neutron yield from (alpha, neutron) reactions with oxygen is calculated by integrating the reaction rate equation over all alpha particle energies and all center-of-mass angles. The results indicate a neutron emission rate of (1.42 plus or minus 0.32) X 10,000 n/s/q PuO2. The neutron yield from (alpha, neutron) reactions with low Z impurities in the fuel is presented in tabular form for one part per million of each impurity. The total neutron flux emitted from a particular fuel geometry is estimated by adding the neutron yield due to the induced fission to the other neutron sources.

Taherzadeh, M.↗

A model for reaction rates in turbulent reacting flows

To account for the turbulent temperature and species-concentration fluctuations, a model is presented on the effects of chemical reaction rates in computer analyses of turbulent reacting flows. The model results in two parameters which multiply the terms in the reaction-rate equations. For these two parameters, graphs are presented as functions of the mean values and intensity of the turbulent fluctuations of the temperature and species concentrations. These graphs will facilitate incorporation of the model into existing computer programs which describe turbulent reacting flows. When the model was used in a two-dimensional parabolic-flow computer code to predict the behavior of an experimental, supersonic hydrogen jet burning in air, some improvement in agreement with the experimental data was obtained in the far field in the region near the jet centerline. Recommendations are included for further improvement of the model and for additional comparisons with experimental data.

Chinitz, W.↗

Dynamics of exploding magma chambers: Implications for K-T volcanism and mass extinctions

Although it is well known that unconfined chemical explosives may yield pressures to several megabars on detonation in air, the explosive literature has yet to be accessed by some contributors to the volcanological literature who've indicated that pressures in excess of the overburden and/or tensile cannot be obtained. Idealized ballistic assessments of pressures internal to volcanoes yield pressures in the hundreds of kilobar range upon correction by addition of friction, etc. Previous assessments of exploding magma chamber pressure have been made from the characteristics of the Mt. St. Helens explosion. A variety of methods yield pressures of similar value: at least hundreds of kilobars. Such results are consistent with free energy requirements for quench supersaturation explosion, a process occurring in solidifying industrial melts. Several reviews of geochemical literature emphasize the carbon event at the Cretaceous-Tertiary (K-T) boundary as being an indicator of a massive dump of CO2 derived from the mantle and entering the atmosphere by extensive global volcanism. Oxygen isotope data indicates extreme warming at the end of the Cretaceous which is consistent with a greenhouse effect attending the CO2 event. Reaction rate equations for the quench supersaturation explosion mechanism indicated, are consistent with the rise in pressure to 30 kbar on solidification of magmatic melts, these pressures limited by the strength of the experimental apparatus.

Rice, A. R.↗

Chemical Reaction and Flow Modeling in Fullerene and Nanotube Production

The development of processes to produce fullerenes and carbon nanotubes has largely been empirical. Fullerenes were first discovered in the soot produced by laser ablation of graphite [1]and then in the soot of electric arc evaporated carbon. Techniques and conditions for producing larger and larger quantities of fullerenes depended mainly on trial and error empirical variations of these processes, with attempts to scale them up by using larger electrodes and targets and higher power. Various concepts of how fullerenes and carbon nanotubes were formed were put forth, but very little was done based on chemical kinetics of the reactions. This was mainly due to the complex mixture of species and complex nature of conditions in the reactors. Temperatures in the reactors varied from several thousand degrees Kelvin down to near room temperature. There are hundreds of species possible, ranging from atomic carbon to large clusters of carbonaceous soot, and metallic catalyst atoms to metal clusters, to complexes of metals and carbon. Most of the chemical kinetics of the reactions and the thermodynamic properties of clusters and complexes have only been approximated. In addition, flow conditions in the reactors are transient or unsteady, and three dimensional, with steep spatial gradients of temperature and species concentrations. All these factors make computational simulations of reactors very complex and challenging. This article addresses the development of the chemical reaction involved in fullerene production and extends this to production of carbon nanotubes by the laser ablation/oven process and by the electric arc evaporation process. In addition, the high-pressure carbon monoxide (HiPco) process is discussed. The article is in several parts. The first one addresses the thermochemical aspects of modeling; and considers the development of chemical rate equations, estimates of reaction rates, and thermodynamic properties where they are available. The second part addresses modeling of the arc process for fullerene and carbon nanotube production using O-D, 1-D and 2-D fluid flow models. The third part addresses simulations of the pulsed laser ablation process using time-dependent techniques in 2-D, and a steady state 2-D simulation of a continuous laser ablation process. The fourth part addresses steady state modeling in O-D and 2-D of the HiPco process. In each of the simulations, there is a variety of simplifications that are made that enable one to concentrate on one aspect or another of the process. There are simplifications that can be made to the chemical reaction models , e.g. reduction in number of species by lumping some of them together in a representative species. Other simulations are carried out by eliminating the chemistry altogether in order to concentrate on the fluid dynamics. When solving problems with a large number of species in more than one spatial dimension, it is almost imperative that the problem be decoupled by solving for the fluid dynamics to find the fluid motion and temperature history of "particles" of fluid moving through a reactor. Then one can solve the chemical rate equations with complex chemistry following the temperature and pressure history. One difficulty is that often mixing with an ambient gas is involved. Therefore, one needs to take dilution and mixing into account. This changes the ratio of carbon species to background gas. Commercially available codes may have no provision for including dilution as part of the input. One must the write special solvers for including dilution in decoupled problems. The article addresses both ful1erene production and single-walled carbon nanotube (SWNT) production. There are at least two schemes or concepts of SWNT growth. This article will only address growth in the gas phase by carbon and catalyst cluster growth and SW T formation by the addition of carbon. There are other models that conceive of SWNT growth as a phase separation process from clusters me up carbon and metal catalyst, with the carbon precipitating from the cluster as it cools. We will not deal with that concept in this article. Further research is needed to determine the rates at which these composite clusters form, evaporate, and segregate.

Scott, Carl D.↗

A study of the Sabatier-methanation reaction kinetics

The kinetics of the Sabatier methanation reaction, the reduction of carbon dioxide with hydrogen to methane and water, was investigated for 58 percent nickel on kieselguhr catalyst and 20 percent ruthenium on alumina catalyst. Differential rate data from an experimental program were correlated with a power function rate equation both for forward and reverse reactions. The kinetic parameters of activation energy, frequency rate constant and reaction order were determined for the rate equation. The values of these parameters were obtained from an Arrhenius plot of the experimental differential rate data. Also the carbon monoxide side reaction effect was measured and included in the correlation of parameters. The reaction was found to fit the rate equation experimentally within the temperature range 421 K, where the reaction effectively begins, the 800 K where the reaction rate drops and departs from the rate equation form.

Verostko, C. E.↗

Algorithm Computes Hypersonic Flow Of Air

Algorithm simulates three-dimensional hypersonic flow of air or another mixture of gases around blunt body. Accounts for nonequilibrium thermochemical effects via two-temperature mathematical model of molecular excitations and equations for finite-rate chemical reactions. These equations fully coupled to fluid dynamical equations.

Palmer, Grant↗

Standard solar model

A set of solar models have been constructed, each based on a single modification to the physics of a reference solar model. In addition, a model combining several of the improvements has been calculated to provide a best solar model. Improvements were made to the nuclear reaction rates, the equation of state, the opacities, and the treatment of the atmosphere. The impact on both the structure and the frequencies of the low-l p-modes of the model to these improvements are discussed. It is found that the combined solar model, which is based on the best physics available (and does not contain any ad hoc assumptions), reproduces the observed oscillation spectrum (for low-l) within the errors associated with the uncertainties in the model physics (primarily opacities).

Guenther, D. B.↗

Standard solar model. II - g-modes

The paper presents the g-mode oscillation for a set of modern solar models. Each solar model is based on a single modification or improvement to the physics of a reference solar model. Improvements were made to the nuclear reaction rates, the equation of state, the opacities, and the treatment of the atmosphere. The error in the predicted g-mode periods associated with the uncertainties in the model physics is predicted and the specific sensitivities of the g-mode periods and their period spacings to the different model structures are described. In addition, these models are compared to a sample of published observations. A remarkably good agreement is found between the 'best' solar model and the observations of Hill and Gu (1990).

Guenther, D. B.↗

Test of a new theory for stellar convection using helioseismology

Two evolutionary models of the sun have been tested using helioseismological data. The two models use the same input microphysics (nuclear reaction rates, opacity, equation of state) and the same numerical evolutionary code, but differ in the treatment of turbulent convection. The first model employs the standard mixing - length theory of convection, while the second one employs a new turbulent convection model which overcomes some basic inconsistencies of the standard theory of convection. The test rests on the calculation of p-mode eigenfrequencies and on the comparison with the helioseismological data. The comparison shows an overall improvement of the eigenfrequencies calculated with the new model with respect to those calculated with the standard model, although it appears that both models still suffer from inaccuracies especially in the treatment of the surface layers.

Paterno, L.↗

Deuterium Exchange in the Systems of H2O+/H2O and H3O+/H2O

Using tandem mass spectrometry various water ion interactions were observed. These reactions consisted of a series of charge transfer, proton transfer, and isotopic exchange steps. The experimental data sets consist of variations of ion abundances over a neutral pressure range. An expected sequence of isotopic exchange reactions is given along with differential equation solutions & reaction rate data.

Spectroscopy Deuterium Proton Ion Atmospheric Scie↗

Investigation of the hydrochlorination of SiCl4

The hyrochlorination of SiC14 and m.g. silicon metal to produce SiHC13, was investigated. Reaction kinetic measurements were carried out to collect additional rate data at 525 C and 550 C. A theoretical study was carried out to provide a kinetic model and a rate equation for the hydrochlorination reaction. Results of this preliminary study show that the rate of formation of SiHC13 follows a pseudo first order kinetics. The rate constants were measured at three temperatures, 550 C, 500 C and 450 C, respectively. The activation energy was determined from the Arrhenius plot to give a value of 13.2 Kcal/mole. The design of a quartz reactor to measure reaction rates and equilibrium conversion of SiHC13 at reaction temperature up to 650 C was completed.

Mui, J. Y. P.↗

Incorporation of a Chemical Equilibrium Equation of State into LOCI-Chem

Renewed interest in development of advanced high-speed transport, reentry vehicles and propulsion systems has led to a resurgence of research into high speed aerodynamics. As this flow regime is typically dominated by hot reacting gaseous flow, efficient models for the characteristic chemical activity are necessary for accurate and cost effective analysis and design of aerodynamic vehicles that transit this regime. The LOCI-Chem code recently developed by Ed Luke at Mississippi State University for NASA/MSFC and used by NASA/MSFC and SSC represents an important step in providing an accurate, efficient computational tool for the simulation of reacting flows through the use of finite-rate kinetics [3]. Finite rate chemistry however, requires the solution of an additional N-1 species mass conservation equations with source terms involving reaction kinetics that are not fully understood. In the equilibrium limit, where the reaction rates approach infinity, these equations become very stiff. Through the use of the assumption of local chemical equilibrium the set of governing equations is reduced back to the usual gas dynamic equations, and thus requires less computation, while still allowing for the inclusion of reacting flow phenomenology. The incorporation of a chemical equilibrium equation of state module into the LOCI-Chem code was the primary objective of the current research. The major goals of the project were: (1) the development of a chemical equilibrium composition solver, and (2) the incorporation of chemical equilibrium solver into LOCI-Chem. Due to time and resource constraints, code optimization was not considered unless it was important to the proper functioning of the code.

Cox, Carey F.↗

Effects of Nonequilibrium Chemistry and Darcy-Forchheimer Pyrolysis Flow for Charring Ablator

The fully implicit ablation and thermal response code simulates pyrolysis and ablation of thermal protection materials and systems. The governing equations, which include energy conservation, a three-component decomposition model, and a surface energy balance, are solved with a moving grid.This work describes new modeling capabilities that are added to a special version of code. These capabilities include a time-dependent pyrolysis gas flow momentum equation with Darcy-Forchheimer terms and pyrolysis gas species conservation equations with finite rate homogeneous chemical reactions. The total energy conservation equation is also enhanced for consistency with these new additions. Two groups of parametric studies of the phenolic impregnated carbon ablator are performed. In the first group, an Orion flight environment for a proposed lunar-return trajectory is considered. In the second group, various test conditions for arcjet models are examined. The central focus of these parametric studies is to understand the effect of pyrolysis gas momentum transfer on material in-depth thermal responses with finite-rate, equilibrium, or frozen homogeneous gas chemistry. Results indicate that the presence of chemical nonequilibrium pyrolysis gas flow does not significantly alter the in-depth thermal response performance predicted using the chemical equilibrium gas model.

Chen, Yih-Kanq↗

Modeling of hypersonic reacting flows

The importance of thermochemical nonequilibrium in hypersonic flow is studied as well as the breakdown of a one-temperature kinetic model. Consideration is given to vibration relaxation, reaction rates, and conservation equations for chemical variables. Experimental data on the pitching moment characteristics of the raked-off blunted elliptic cone geometry proposed for use with an aeroassisted orbit transfer vehicle are presented.

Park, Chul↗

Helioseismological test of a new model for stellar convection

We compare the predictions of two solar models with the observed p-mode eigenfrequencies. The two models use the same input microphysics (nuclear reaction rates, opacity, and equation of a state) and the same numerical evolutionary code, but differ in the treatment of turbulent convection. The first model employs the standard mixing-length theory of convection (MLT), while the second model employs a new model of turbulent convection (CM) whose primary goal was that of accounting for the whole spectrum of turbulent eddies so as to avoid the MLT approximation that such a wide spectrum be represented by a single, large eddy. The main result of this paper is that the p-mode eigenfrequencies calculated with the CM model show an overall improvement with respect to those calculated with the standard MLT model.

Paterno, L.↗