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.
Estimate of chemical kinetics behind normal shock waves in carbon dioxide and nitrogen mixtures for Martian atmospheric entry
Combustion instability in hybrid rocket motors, noting chemical kinetics role and steady state regression rate
Martian entry condition analyzed as nonequilibrium flow field phenomenon, noting chemical kinetics and radiation behind shock wave in gas mixtures
Laminar boundary layers effect on chemical-kinetic measurements in shock tube, applying correction by computer program for subsonic diffuser
A general chemical kinetics code is described for complex, homogeneous ideal gas reactions in any chemical system. The main features of the GCKP84 code are flexibility, convenience, and speed of computation for many different reaction conditions. The code, which replaces the GCKP code published previously, solves numerically the differential equations for complex reaction in a batch system or one dimensional inviscid flow. It also solves numerically the nonlinear algebraic equations describing the well stirred reactor. A new state of the art numerical integration method is used for greatly increased speed in handling systems of stiff differential equations. The theory and the computer program, including details of input preparation and a guide to using the code are given.
A number of chemical-kinetic problems related to phenomena occurring behind a shock wave surrounding an object flying in the earth atmosphere are discussed, including the nonequilibrium thermochemical relaxation phenomena occurring behind a shock wave surrounding the flying object, problems related to aerobraking maneuver, the radiation phenomena for shock velocities of up to 12 km/sec, and the determination of rate coefficients for ionization reactions and associated electron-impact ionization reactions. Results of experiments are presented in form of graphs and tables, giving data on the reaction rate coefficients for air, the ionization distances, thermodynamic properties behind a shock wave, radiative heat flux calculations, Damkoehler numbers for the ablation-product layer, together with conclusions.
Thermoconductivity in a chemically reacting gas mixture with one finite reaction rate
The purpose of this study is to provide a capability for determining the chemical kinetic behavior of one family of pollutants, the nitrogen oxides, in the flow of a coal-fired MHD generator facility. The method used in the study allows the prediction of the nonequilibrium concentration of the minor NO(x) species in a flow otherwise assumed in equilibrium. Consideration is given to the effect of preheat, stoichiometry, and oxygen enrichment of the NO(x) concentration. The effect of preheat and stoichiometry is found to have a significant influence on the NO(x) concentration at the exit of the radiant boiler.
The development of chemical kinetics for the simulation of propane-air combustion with the use of computer code KIVA-II since 1989 is summarized here. In order to let readers understand the general feature well, a brief description of the KIVA-II code, specially related with the chemical reactions is also given. Then the results of recent work with 20 reaction mechanism is presented. It is also compared with the 5 reaction mechanism. It may be expected that the numerical stability of the 20 reaction mechanism is better as compared to that of 5 reaction mechanism, but the CPU time of the CRAY computer is much longer. Details are presented in the paper.
Turbulent boundary layers produced by nonuniform shock tube flow, considering effects on chemical kinetic measurements
We have previously developed a chemical conversion model of the carbothermal processing of lunar regolith using methane to predict the rate of production of carbon monoxide. In this carbothermal process, gaseous methane is pyrolyzed as it flows over the hot surface of a molten zone of lunar regolith and is converted to carbon and hydrogen. Hydrogen is carried away by the exiting stream of gases and carbon is deposited on the melt surface. The deposited carbon mixes with the melt and reacts with the metal oxides in it to produce carbon monoxide that bubbles out of the melt. In our model, we assume that the flux of carbon deposited is equal to the product of the surface reaction rate constant gamma and the concentration of methane adjacent to the melt surface. Similarly, the rate of consumption of carbon per unit volume in the melt is equal to the product of the melt reaction rate constant k and the concentrations of carbon and metal oxide in the melt. In this paper, we describe our effort to determine gamma and k by comparison of the predictions from our model with test data obtained by ORBITEC (Orbital Technologies Corporation). The concentration of methane adjacent to the melt surface is a necessary input to the model. It is inferred from the test data by a mass balance of methane, adopting the usual assumptions of the continuously-stirred-tank-reactor model, whereby the average concentration of a given gaseous species equals its exit concentration. The reaction rates gamma and k have been determined by a non-linear least-squares fit to the test data for the production of carbon monoxide and the fraction of the incoming methane that is converted. The comparison of test data with our model predictions using the determined chemical kinetic rate constants provides a consistent interpretation of the process over the full range of temperatures, pressures, and methane flow rates used in the tests, thereby increasing our confidence to use the model for scale-up purposes.
Chemical kinetic model for hydrocarbon fuels combustion and application to multidimensional finite difference mixing analyses in hypersonic engines and nozzles
Numerical analysis of chemical kinetics of hydrogen-air combustion in presence of carbon dioxide and water vapor
A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas phase chemical reactions at constant pressure, and include all three combustion regimes: induction, heat release, and equilibration. The comparisons show that LSODE is the most efficient code - in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.
A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas phase chemical reactions at constant pressure, and include all three combustion regimes: induction heat release, and equilibration. The comparisons show that LSODE is the most efficient code - in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.
The effect of switching between nonstiff and stiff methods on the efficiency of algorithms for integrating chemical kinetic rate equations was examined. Different integration methods were tested by application of the packaged code LSODE to four practical combustion kinetics problems. The problems describe adiabatic, and homogeneous gas phase combustion reactions. It is shown that selective use of nonstiff and stiff methods in different regimes of a typical batch combustion problem is faster than the use of either method for the entire problem. The implications which result in the development of fast integration techniques for combustion kinetic rate equations are discussed.
The effect of switching between nonstiff and stiff methods on the efficiency of algorithms for integrating chemical kinetic rate equations is presented. Different integration methods are tested by application of the packaged code LSODE to four practical combustion kinetics problems. The problems describe adiabatic, homogeneous gas-phase combustion reactions. It is shown that selective use of nonstiff and stiff methods in different regimes of a typical batch combustion problem is faster than the use of either method for the entire problem. The implications of this result to the development of fast integration techniques for combustion kinetic rate equations are discussed.
A CO-CH4-air reaction system is used to demonstrate the computational singular perturbation (CSP) method for deriving time-resolved simplified chemical kinetics models. CSP provides a programmable algorithm to group the given collection of elementary reactions into reaction groups which are ordered according to their speed. The concept of Importance Index k(m)exp s is introduced: k(m)exp s is defined to be a number between 0 and 1 which measures the importance of the m-th reaction group to the s-th reactant and can readily be computed from data generated by CSP. It is suggested that the robustness of the solutions of the reaction system can be qualitatively assessed by inspecting the Importance Index data.