Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reaction Rates”

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 163 records · Page 9

Solutions of the chemical kinetic equations for initially inhomogeneous mixtures.

Following the recent discussions by O'Brien (1971) and Donaldson and Hilst (1972) of the effects of inhomogeneous mixing and turbulent diffusion on simple chemical reaction rates, the present report provides a more extensive analysis of when inhomogeneous mixing has a significant effect on chemical reaction rates. The analysis is then extended to the development of an approximate chemical sub-model which provides much improved predictions of chemical reaction rates over a wide range of inhomogeneities and pathological distributions of the concentrations of the reacting chemical species. In particular, the development of an approximate representation of the third-order correlations of the joint concentration fluctuations permits closure of the chemical sub-model at the level of the second-order moments of these fluctuations and the mean concentrations.

Hilst, G. R.↗

Measurements of HO2 chemical kinetics with a new detection method

In this research program, HO2 was detected by the OH(A-X) photofragment from dissociative excitation of HO2 at 147 nm. This detection method was applied to measure the reaction rate constant of HO2 + O3. This reaction rate constant is needed for the understanding of stratospheric chemistry. Since C12 was used in the flow system, photoexcitation of C12 may produce fluorescence to interfere with the measurements. Thus, the photoexcitation process of C12 in the vacuum ultraviolet region was also examined in this research period using synchrotron radiation as a light source. The research results are summarized.

Lee, L. C.↗

Effects of the low Earth orbital environment on spacecraft materials

It is evident from space flights during the last three years that the low Earth orbital (LEO) environment interacts with spacecraft surfaces in significant ways. One manifestation of these interactions is recession of, in particular, organic-polymer-based surfaces presumably due to oxidation by atomic oxygen, the major component of the LEO environment. Three experiments have been conducted on Space Shuttle flights 5, 8 and 41-G to measure reaction rates and the effects of various parameters on reaction rates. Surface recession on these flights indicates reaction efficiencies approximately 3 x 10(-24) cu cm/atoms for unfilled organic polymers. Of the metals, silver and osmium are very reactive. Effects on spacecraft or experiment surfaces can be evaluated using the derived reaction efficiencies and a definition of the total exposure to atomic oxygen. This exposure is obtained using an ambient density model, solar activity data and spacecraft parameters of altitude, attitude and operational date. Oxygen flux on a given surface is obtained from the ambient density and spacecraft velocity and can then be integrated to provide the total exposure or fluence. Such information can be generated using simple computational programs and can be converted to various formats. Overall, the extent of damage is strongly dependent on the type of surface and total exposure time.

Leger, L. J.↗

Experiments concerning the laser-enhanced reaction between vibrationally excited O3 and NO

The enhancement in reaction rate between O3 and NO is studied for the case of O3 vibrationally excited by a CO2 laser. Chemiluminescence observations of a vibrationally excited and an electronically excited nitrogen dioxide reaction product provide information on the separate contributions to the overall reaction rate of these two reactive channels. The contribution of the stretching and bending modes of O3 to the reaction rate enhancement is also discussed. In addition, consideration is given to the nonreactive vibrational deactivation of vibrationally excited O3.

Hui, K.-K.↗

Progress in the development of PDF turbulence models for combustion

A combined Monte Carlo-computational fluid dynamic (CFD) algorithm was developed recently at Lewis Research Center (LeRC) for turbulent reacting flows. In this algorithm, conventional CFD schemes are employed to obtain the velocity field and other velocity related turbulent quantities, and a Monte Carlo scheme is used to solve the evolution equation for the probability density function (pdf) of species mass fraction and temperature. In combustion computations, the predictions of chemical reaction rates (the source terms in the species conservation equation) are poor if conventional turbulence modles are used. The main difficulty lies in the fact that the reaction rate is highly nonlinear, and the use of averaged temperature produces excessively large errors. Moment closure models for the source terms have attained only limited success. The probability density function (pdf) method seems to be the only alternative at the present time that uses local instantaneous values of the temperature, density, etc., in predicting chemical reaction rates, and thus may be the only viable approach for more accurate turbulent combustion calculations. Assumed pdf's are useful in simple problems; however, for more general combustion problems, the solution of an evolution equation for the pdf is necessary.

Hsu, Andrew T.↗

Effect of photochemical models on calculated equilibria and cooling rates in the stratosphere

The determination of the relaxation time of a temperature perturbation in the stratosphere must take into account the effects of the absorption of solar energy by ozone, while the ozone density itself is dependent on temperature. A photochemical model, consisting of continuity equations for each of 29 constituents with reaction rates and adjustment times, is used to obtain the vertical distribution of ozone. The altitude range from 35 to 60 km in non-Arctic regions is shown to be in approximate joint radiative-photochemical equilibrium. The temperature and ozone distributions are thus well buffered. Modest changes in ozone and temperature were calculated on the basis of this model from large changes in cooling rates and reaction rates, and the results are shown to be more in line with actual observations. The importance of estimating the mixing ratios for NO and H2O is emphasized.

Blake, D.↗

Flame propagation under partially-premixed conditions

This study concentrates on developing a better understanding of triple flames. We relax the assumption of zero heat release, address the issue of stabilization, and, in order to investigate the role that heat release plays in flame propagation in partially premixed combustion, we return to a simple flow field and investigate the behavior of flames in a laminar environment. We solve the compressible Navier-Stokes equations in a two-dimensional domain. At the boundaries, we use an inflow boundary condition on the left and nearly-perfect reflective boundary conditions, required to avoid pressure drift, at the outflow and sides. After the flow and flame are initialized, the mixture fraction is varied at the inlet from its uniform stoichiometric value to a tanh profile varying from zero to one. As the mixture fraction gradient reaches the flame surface only the centerline is exposed to the stoichionetric mixture fraction and locally maintains the planar flame speed and reaction rate. Above this point the mixture is fuel rich, and below fuel lean. As a result, these regions of non-unity equivalence ratio burn less, the reaction rate drops, and the local flame speed is reduced. The excess fuel and oxidizer then combine behind the premixed flame along the stoichiometric surface and burn in a trailing diffusion flame. Thus the 'triple' flame refers to the fuel-rich premixed flame, the fuel-lean premixed flame, and the trailing diffusion flame. Due to heat release, the normal velocity across the flame is increased, whereas the tangential component remains unchanged. Far-field flame speed, local flame speed, and their differences are shown as a function of the local mixing thickness. It was also determined that the lateral position of the flame affects stabilization, and the distribution of the reaction rate along the premixed wings of triple flames affects flame propagation.

Ruetsch, Gregory R.↗

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗

EUPDF: Eulerian Monte Carlo Probability Density Function Solver for Applications With Parallel Computing, Unstructured Grids, and Sprays

The success of any solution methodology used in the study of gas-turbine combustor flows depends a great deal on how well it can model the various complex and rate controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as convective and radiative heat transfer and other phenomena. The phenomena to be modeled, which are controlled by these processes, often strongly interact with each other at different times and locations. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. The influence of turbulence in a diffusion flame manifests itself in several forms, ranging from the so-called wrinkled, or stretched, flamelets regime to the distributed combustion regime, depending upon how turbulence interacts with various flame scales. Conventional turbulence models have difficulty treating highly nonlinear reaction rates. A solution procedure based on the composition joint probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices (such as extinction, blowoff limits, and emissions predictions) because it can account for nonlinear chemical reaction rates without making approximations. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on the PDF method to unstructured grids, parallel computing, and sprays. EUPDF, which was developed by M.S. Raju of Nyma, Inc., was designed to be massively parallel and could easily be coupled with any existing gas-phase and/or spray solvers. EUPDF can use an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements. The application of the PDF method showed favorable results when applied to several supersonic-diffusion flames and spray flames. The EUPDF source code will be available with the National Combustion Code (NCC) as a complete package.

Raju, M. S.↗

Summary of Simplified Two Time Step Method for Calculating Combustion Rates and Nitrogen Oxide Emissions for Hydrogen/Air and Hydrogen/Oxygen

A simplified single rate expression for hydrogen combustion and nitrogen oxide production was developed. Detailed kinetics are predicted for the chemical kinetic times using the complete chemical mechanism over the entire operating space. These times are then correlated to the reactor conditions using an exponential fit. Simple first order reaction expressions are then used to find the conversion in the reactor. The method uses a two time step kinetic scheme. The first time averaged step is used at the initial times with smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, temperature, and pressure. The second instantaneous step is used at higher water concentrations (greater than l x 10(exp -20)) moles per cc) in the mixture which gives the chemical kinetic time as a function of the instantaneous fuel and water mole concentrations, pressure and temperature (T(sub 4)). The simple correlations are then compared to the turbulent mixing times to determine the limiting properties of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates are used to calculate the necessary chemical kinetic times. This time is regressed over the complete initial conditions using the Excel regression routine. Chemical kinetic time equations for H2 and NOx are obtained for H2/Air fuel and for H2/O2. A similar correlation is also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium temperature (T(sub 4)) as a function of overall fuel/air ratio, pressure and initial temperature (T(sub 3)). High values of the regression coefficient R squared are obtained.

Marek, C. John↗

Simplified Two-Time Step Method for Calculating Combustion Rates and Nitrogen Oxide Emissions for Hydrogen/Air and Hydorgen/Oxygen

A simplified single rate expression for hydrogen combustion and nitrogen oxide production was developed. Detailed kinetics are predicted for the chemical kinetic times using the complete chemical mechanism over the entire operating space. These times are then correlated to the reactor conditions using an exponential fit. Simple first order reaction expressions are then used to find the conversion in the reactor. The method uses a two-time step kinetic scheme. The first time averaged step is used at the initial times with smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, temperature, and pressure. The second instantaneous step is used at higher water concentrations (> 1 x 10(exp -20) moles/cc) in the mixture which gives the chemical kinetic time as a function of the instantaneous fuel and water mole concentrations, pressure and temperature (T4). The simple correlations are then compared to the turbulent mixing times to determine the limiting properties of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates are used to calculate the necessary chemical kinetic times. This time is regressed over the complete initial conditions using the Excel regression routine. Chemical kinetic time equations for H2 and NOx are obtained for H2/air fuel and for the H2/O2. A similar correlation is also developed using data from NASA s Chemical Equilibrium Applications (CEA) code to determine the equilibrium temperature (T4) as a function of overall fuel/air ratio, pressure and initial temperature (T3). High values of the regression coefficient R2 are obtained.

Molnar, Melissa↗

Application of mixing-controlled combustion models to gas turbine combustors

Gas emissions were studied from a staged Rich Burn/Quick-Quench Mix/Lean Burn combustor were studied under test conditions encountered in High Speed Research engines. The combustor was modeled at conditions corresponding to different engine power settings, and the effect of primary dilution airflow split on emissions, flow field, flame size and shape, and combustion intensity, as well as mixing, was investigated. A mathematical model was developed from a two-equation model of turbulence, a quasi-global kinetics mechanism for the oxidation of propane, and the Zeldovich mechanism for nitric oxide formation. A mixing-controlled combustion model was used to account for turbulent mixing effects on the chemical reaction rate. This model assumes that the chemical reaction rate is much faster than the turbulent mixing rate.

Nguyen, Hung Lee↗

Reduced Equations for Calculating the Combustion Rates of Jet-A and Methane Fuel

Simplified kinetic schemes for Jet-A and methane fuels were developed to be used in numerical combustion codes, such as the National Combustor Code (NCC) that is being developed at Glenn. These kinetic schemes presented here result in a correlation that gives the chemical kinetic time as a function of initial overall cell fuel/air ratio, pressure, and temperature. The correlations would then be used with the turbulent mixing times to determine the limiting properties and progress of the reaction. A similar correlation was also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium concentration of carbon monoxide as a function of fuel air ratio, pressure, and temperature. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates and the values obtained from the equilibrium correlations were then used to calculate the necessary chemical kinetic times. Chemical kinetic time equations for fuel, carbon monoxide, and NOx were obtained for both Jet-A fuel and methane.

Molnar, Melissa↗

Direct numerical simulation of turbulent H2-O2 combustion using reduced chemistry

Results of direct numerical simulations of hydrogen-oxygen combustion using a partial-equilibrium chemistry scheme in constant density, decaying, isotropic turbulence are reported. The simulations qualitatively reproduce many features of experimental results, such as superequilibrium radical species mole fractions, with temperature and major species mole fractions closer to chemical equilibrium. It was also observed that the peak reaction rates occur in narrow zones where the stoichiometric surface intersects regions of high scalar dissipation, as might be expected for combustion conditions close to chemical equilibrium. Another finding was that high OH mole fraction correspond more closely to the stoichiometric surface than to areas of high reaction rate for conditions of the simulations. Simulation results were compared to predictions of the Conditional Moment Closure model. This model was found to give good results for all quantities of interest when the conditionally averaged scalar dissipation was used in the prediction. When the nonconditioned average dissipation was used, the predictions compared well to the simulations for most of the species and temperature, but not for the reaction rate. The comparison would be expected to improve for higher Reynolds number flows, however.

Montgomery, Christopher J.↗

The rate of pyrite decomposition on the surface of Venus

We report the results of a detailed experiment study of the kinetics and mechanism of pyrite (FeS2) chemical weathering under Venus surface conditions. Pyrite is thermodynamically unstable on the surface of Venus and will spontaneously decompose to pyrrhotite (Fe7S8) because the observed S2 partial pressure in the lower atmosphere of Venus is lower than the S2 vapor pressure over coexisting pyrite and pyrrhotite. Pyrite decomposition kinetics were studied in pure CO2 and CO2 gas mixtures along five isotherms in the temperature range 390-531 C. In all gas mixtures studied, pyrite thermally decomposes to pyrrhotite (Fe7S8), which on continued heating loses sulfur to form more Fe-rich pyrrhotites. During this process the pyrrhotites are also being oxidized to form magnetite (Fe3O4), which converts to maghemite (gamma-Fe2O3), and then to hematite (alpha-Fe2O3). The reaction rates for pyrite thermal decomposition to pyrrhotite were determined by measuring the weight loss. The thickness of the unreacted pyrite in the samples provided a second independent reaction rate measurement. Finally, Mossbauer spectra done on 42 of the 115 experimental samples provided a third set of independent reaction rate data. Pyrite decomposition follows zero-order kinetics and is independent of the amount of pyrite present. The rate of pyrite decomposition is apparently independent of the gas compositions used and of the CO2 number density over a range of a factor of 40. The derived activation energy of approximately 150 kJ/mole is the same in pure CO2, two different CO-CO2 mixtures, and a ternary CO-SO2-CO2 mixture. Based on data for a CO-CO2-SO2 gas mixture with a CO number density approximately 10 times higher than at the surface of Venus and a SO2 number density approximately equal to that at the surface of Venus, the rate of pyrite destruction on the surface of Venus varies from about 1225 +/- 238 days/cm at the top of Maxwell Montes (approximately 660 K) to about 233 +/- 133 days/cm in the plains of Venus (approximately 740 K). These lifetimes are very short on a geological time scale and show that pyrite cannot exist on the surface of Venus for any appreciable length of time.

Fegley, B., Jr.↗

Determination of Chemical Kinetic Rate Constants of a Model for Carbothermal Processing of Lunar Regolith Simulant Using Methane

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.

Balasubramaniam, R↗

Rate of reaction between molecular hydrogen and molecular oxygen

The shock tube data of Jachimowski and Houghton were rigorously analyzed to obtain rate constants for the candidate initiation reactions H2 + O2 yields H + HO2, H2 + O2 yields H2O + O, and H2 + O2 yields OH + OH. Reaction (01) is probably not the initiation process because the activation energy obtained is less than the endothermicity and because the derived rates greatly exceed values inferred in the literature from the reverse of reaction (01). Reactions (02) and (03) remain as possibilities, with reaction (02) slightly favored on the basis of steric and statistical considerations. The solution of the differential equations is presented in detail to show how the kinetics of other ignition systems may be solved.

Brokaw, R. S.↗

New Reduced Two-Time Step Method for Calculating Combustion and Emission Rates of Jet-A and Methane Fuel With and Without Water Injection

A simplified kinetic scheme for Jet-A, and methane fuels with water injection was developed to be used in numerical combustion codes, such as the National Combustor Code (NCC) or even simple FORTRAN codes that are being developed at Glenn. The two time step method is either an initial time averaged value (step one) or an instantaneous value (step two). The switch is based on the water concentration in moles/cc of 1x10(exp -20). The results presented here results in a correlation that gives the chemical kinetic time as two separate functions. This two step method is used as opposed to a one step time averaged method previously developed to determine the chemical kinetic time with increased accuracy. The first time averaged step is used at the initial times for smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, initial water to fuel mass ratio, temperature, and pressure. The second instantaneous step, to be used with higher water concentrations, gives the chemical kinetic time as a function of instantaneous fuel and water mole concentration, pressure and temperature (T4). The simple correlations would then be compared to the turbulent mixing times to determine the limiting properties of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates were then used to calculate the necessary chemical kinetic times. Chemical kinetic time equations for fuel, carbon monoxide and NOx were obtained for Jet-A fuel and methane with and without water injection to water mass loadings of 2/1 water to fuel. A similar correlation was also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium concentrations of carbon monoxide and nitrogen oxide as functions of overall equivalence ratio, water to fuel mass ratio, pressure and temperature (T3). The temperature of the gas entering the turbine (T4) was also correlated as a function of the initial combustor temperature (T3), equivalence ratio, water to fuel mass ratio, and pressure.

Molnar, Melissa↗