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 307 records · Page 17

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.↗

An Eulerian multimaterial framework for simulating high-explosive aquarium tests

Aquarium tests of cylindrical high-explosive charges provide optical data of the detonation front velocity and shape, propagation of the shock wave in the surrounding water, and expansion rates of the detonation products behind the front. Data from aquarium experiments is often used for calibration of reactive burn models based on phenomenological equations of state (EOS) and reaction rate laws. This paper presents a multimaterial numerical modeling framework to solve the 2D axisymmetric reactive Euler equations for high-explosive aquarium tests, in particular for ammonium nitrate - fuel oil (ANFO) explosives. An extension of the Ghost Fluid Method (GFM) is used to handle the dynamic material interfaces for the ANFO explosion products, the charge-confining material (polymethyl methacrylate PMMA), and the surrounding water. This study analyzes the sensitivity of calculations (both computational efficiency and numerical accuracy) to different algorithms for the material interface models including the original GFM versus Riemann solver-based strategies. A novel method for defining the left and right states in the interfacial Riemann problem eliminates the need for sorting or nodal interpolation during the projection along the material interface. Numerical tests indicate that populating the interface node values using the Riemann solution mitigate the overheating error observed in steady-state calculations. Solution convergence and computational efficiency are explored as a function of the spatial and temporal order of the schemes. Results from the computational model with analytical equations of state and fitted reaction rate parameters show very good quantitative agreement with experimentally observed detonation front velocity, reaction products expansion, and shock wave propagation in the surrounding water for a cylindrical ANFO charge encased in PMMA. Finally, the proposed modeling framework, in conjunction with experimental tests, provides a reliable tool to assess equations of state and reaction rate expressions for reactive burn models of confined high explosives.

42 ENGINEERING↗

Electrocatalytic Valorization into H2 and Hydrocarbons of an Aqueous Stream Derived from Hydrothermal Liquefaction

Electrocatalytic oxidation is an attractive process for valorizing the carbon and removal of nitrogen present in the aqueous waste streams at ambient conditions. Here, we evaluated the electrocatalytic oxidation as a function of applied potential over a commercial Pt electrode of aqueous waste stream generated via hydrothermal liquefaction. We validated the conversion of quantifiable carbon (e.g., carboxylic acids, alpha hydroxyacids, alcohols, ketones, amides) and unquantifiable carbon as well as nitrogen removal. The unquantifiable carbon was valorized to quantifiable short chain organic molecules (e.g., acetic). Model studies showed that carboxylic acids, alcohol and ketones converted following (non-) Kolbe chemistry. Amides and alpha hydroxyacids formed carboxylic acids that subsequently converted via (non-) Kolbe chemistry. Ammonia was oxidized into N2. The reaction rate varied for each molecule family with carboxylic acids and ammonia having the highest reaction rates. The main reaction products from the electrocatalytic oxidation at the anode were volatile, short chain hydrocarbons (i.e., olefins, paraffins) and CO generated via Kolbe chemistry, while H2 was generated at the cathode. The parallel denitrification, valorization of carbon, and H2 generation of the aqueous waste can replace and simplify the current unit operations used in a commercial hydrothermal liquefaction process. The high operation potential required to drive the electrocatalytic reaction is responsible for high operation costs; however, they can be either partially subsidized with sale of excess H2 at the current DOE goal of $2/kg H2. Thus, the present work shows how electrocatalytic oxidation can be used to valorize aqueous waste streams into volatile hydrocarbons and H2.

Electrocatalytic oxidation, Aqueous waste valoriza↗

Impact of Pycnonuclear Fusion Uncertainties on the Cooling of Accreting Neutron Star Crusts

Abstract The observation of X-rays during quiescence from transiently accreting neutron stars provides unique clues about the nature of dense matter. This, however, requires extensive modeling of the crusts and matching the results to observations. The pycnonuclear fusion reaction rates implemented in these models are theoretically calculated by extending phenomenological expressions and have large uncertainties spanning many orders of magnitude. We present the first sensitivity studies of these pycnonuclear fusion reactions in realistic network calculations. We also couple the reaction network with the thermal evolution code dStar to further study their impact on the neutron star cooling curves in quiescence. Varying the pycnonuclear fusion reaction rates alters the depth at which nuclear heat is deposited although the total heating remains constant. The enhancement of the pycnonuclear fusion reaction rates leads to an overall shallower deposition of nuclear heat. The impurity factors are also altered depending on the type of ashes deposited on the crust. These total changes correspond to a variation of up to 9 eV in the modeled cooling curves. While this is not sufficient to explain the shallow heat source, it is comparable to the observational uncertainties and can still be important for modeling the neutron star crust.

79 ASTRONOMY AND ASTROPHYSICS↗

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.↗

CO 2 hydrogenation over rhodium cluster catalyst nucleated within a manganese oxide framework

Rhodium-based manganese oxide frameworks were explored as a prototype for carbon dioxide reactive capture and conversion. Three-dimensional frameworks of MnOx were utilized as support structures to isolate Rh metal centers. V, Na, and Zn were introduced as counterions to stabilize the structure and for their beneficial effect as promoters. Here, with this multicomponent catalyst, Rh active centers with MnOxs and varied counterions, we were able to selectively tune the catalytic performance of the material via the choice of counterion and structure of the host material. With cryptomelane-type tunnel manganese oxides octahedral molecular sieve (OMS2), we found that Rh-V-OMS2 was highly stable even after 48 hours on stream with a reaction rate of around 1.5x10 -4 mol CO 2 /g Rh /s, surpassing the net reactivity of other initially more active combinations. Furthermore, during CO 2 hydrogenation, in situ XAFS showed that single Rh atoms nucleated into nanoparticles/ sub-nanometer clusters with a coordination number of 5.5 or less. Our finding of the correlation between the reaction rate and particle size offers the potential for enhanced control over the reaction rate by tuning particle size. Our activity study with control experiments demonstrates that the activities of the catalysts are proved due to the unique metal support interaction offered by the Rh-X-MnO.

10 SYNTHETIC FUELS↗

Oxidation and Associated Pore Structure Modification During Experimental Alteration of Granite

Weathering plays a crucial role in a number of environmental processes, and the microstructure and evolution of multi-scale pore space is a critically important factor in the weathering of geological formations. In igneous rocks the infiltration of meteoric water into initially relatively dry material can cause the rock to begin to disaggregate, increasing porosity and surface area, and allowing further disaggregation and weathering to occur. These processes, in turn, allow biota to colonize the rock, further enhancing the weathering rate. In some igneous rocks this may be driven by oxidation of primary minerals. One such mineral, biotite, has been repeatedly mentioned as a cause of cracking during oxidation. However, the scale-dependence of the processes by which this occurs are poorly understood. We cannot, therefore, accurately extrapolate laboratory reaction rates to the field in predictive numerical models. In order to better understand the effects of oxidation and test the hypothesis that fracture and disaggregation are initiated by swelling of oxidizing biotites, we reacted granite cores in a selenic acid-rich aqueous solution at 200°C for up to 438 days. Elevated temperatures and selenic acid were used to provide relatively fast reaction rates and highly oxidizing conditions in sealed reaction vessels. These experiments were analyzed using a combination of imaging, X-ray diffraction, Mössbauer spectroscopy, and small- and ultra-small angle neutron scattering to interrogate porosity and microfracture formation. The experimental results show little observable biotite swelling, but significant transport and growth of iron oxides and/or clays along grain boundaries throughout the sample. Significant increases in porosity were also observed at the sample rim, likely associated with feldspar alteration. Fractures and transport were observed throughout the core, suggesting that stresses due to crystallization pressures caused by the growing iron phases may be the initiating factors in granite weathering, possibly followed by biotite swelling after sufficient permeability is achieved.

Anovitz, Lawrence M.↗

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.↗

Accuracy of High Explosive Data used for Reactive Burn Model Calibration

Reactive burn models for a high explosive (HE) have three components: an equation of state (EOS) for the products, an EOS for the reactants and a reaction rate. For the products, we discuss the uncertainty in data for the planar detonation state and its release isentrope. For the reactants, the uncertainty in the principal shock locus data is discussed. For the reaction rate, we discuss uncertainties in shock initiation data (distance-of-run to detonation as a function of shock pressure), and curvature effect data (detonation speed as a function of front curvature). An important source of uncertainty for calibrating the reaction rate for a plastic-bonded explosive (PBX) results from variations in heterogeneities; that is the mesoscale structure of HE grains, binder and porosity. This shows up in the scatter of both shock initiation data and curvature effect data with small changes in initial density and initial temperature, and systematic changes with PBX lot. As illustrative examples, we use data for the conventional HMX based PBX 9501 and the insensitive TATB based PBX 9502 since these explosives have been extensively studied and high quality data are readily available in the literature.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

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↗

Estimation of the NiCu Cycle Strength and Its Impact on Type I X-Ray Bursts

Abstract Type I X-ray bursts (XRBs) are powered by thermonuclear burning on proton-rich unstable nuclides. The construction of burst models with accurate knowledge of nuclear physics is required to properly interpret burst observations. Numerous studies that have investigated the sensitivities of burst models to nuclear inputs have commonly extracted the strength of the NiCu cycle in the rp process, determined by the 59 Cu( p , α ) 56 Ni and 59 Cu( p , γ ) 60 Zn thermonuclear reaction rates, as critical in the determination of reaction flow in the burst. In this study, the strength of the cycle at the XRB temperature range was estimated based on published experimental data. The nuclear properties of the compound nucleus 60 Zn were evaluated for the 59 Cu( p , α ) 56 Ni and 59 Cu( p , γ ) 60 Zn reaction rate calculations. Monte Carlo rate calculations were conducted to include the large uncertainties of nuclear properties in the calculations. In the current work, a weak NiCu cycle is expected, whereas the rates adopted by the previous studies suggest a strong NiCu cycle. Model simulations were performed with the new rates to assess the impact on Type I XRBs. The results show that the estimated cycle strength does not strongly influence the model predictions of the burst light curve or synthesized abundances.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗