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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Projecting the Performance of a Continuous Bosch Carbon Dioxide Reduction System

Kinetic studies using a subscale reactor were used to develop power-law rate equations for the two primary reactions contributing to the Bosch carbon dioxide reduction process: the reverse water-gas shift reaction and the carbon monoxide hydrogenation reaction. An alloy catalyst in bead form was used in these studies, which investigated both the dependence on gas composition and the dependence on time as the beads gradually disintegrate during carbon formation. Bead disintegration impacts both the rates and the relative activities of the two primary reactions. Although the kinetic testing was conducted in a single-pass mode, the target Bosch application requires operation in a recycle mode to achieve high conversion. Further, the proposed design concept for a full-scale carbon dioxide reduction system would allow carbon to be removed continuously in contrast to previous Bosch processes that required replacement of expendable catalyst cartridges. The developed power-law kinetic rate equations were used in conjunction with a reactor model to project the operating conditions, performance, and sizing of a continuous Bosch system.

carbon dioxide reduction↗
NASA NTRS1979

Simplified model of a volumetric direct nuclear pumped He-3-Ar laser

A physical model for the volumetric nuclear pumped He-3-Ar laser is described and the source terms arising from the He-3(n, p)H-3 reaction as well as the rate equations are given. Steady state solutions are derived for the densities of the species vs the fraction of argon. The dominant pumping mechanism is found to be collisional radiative recombination of the argon atomic ions and subsequent cascading into the upper laser level. The reactions following recombination are described and an expression is derived for the difference in densities of the argon 3d and 4p levels. Comparison with experiments shows good agreement.

Harries, W. L.↗
NASA NTRS1989

Boost-phase discrimination research activities

Theoretical research in two areas was performed. The aerothermodynamics research focused on the hard-body and rocket plume flows. Analytical real gas models to describe finite rate chemistry were developed and incorporated into the three-dimensional flow codes. New numerical algorithms capable of treating multi-species reacting gas equations and treating flows with large gradients were also developed. The computational chemistry research focused on the determination of spectral radiative intensity factors, transport properties and reaction rates. Ab initio solutions to the Schrodinger equation provided potential energy curves transition moments (radiative probabilities and strengths) and potential energy surfaces. These surfaces were then coupled with classical particle reactive trajectories to compute reaction cross-sections and rates.

Cooper, David M.↗
NASA NTRS2008

Simplified Modeling of Oxidation of Hydrocarbons

A method of simplified computational modeling of oxidation of hydrocarbons is undergoing development. This is one of several developments needed to enable accurate computational simulation of turbulent, chemically reacting flows. At present, accurate computational simulation of such flows is difficult or impossible in most cases because (1) the numbers of grid points needed for adequate spatial resolution of turbulent flows in realistically complex geometries are beyond the capabilities of typical supercomputers now in use and (2) the combustion of typical hydrocarbons proceeds through decomposition into hundreds of molecular species interacting through thousands of reactions. Hence, the combination of detailed reaction- rate models with the fundamental flow equations yields flow models that are computationally prohibitive. Hence, further, a reduction of at least an order of magnitude in the dimension of reaction kinetics is one of the prerequisites for feasibility of computational simulation of turbulent, chemically reacting flows. In the present method of simplified modeling, all molecular species involved in the oxidation of hydrocarbons are classified as either light or heavy; heavy molecules are those having 3 or more carbon atoms. The light molecules are not subject to meaningful decomposition, and the heavy molecules are considered to decompose into only 13 specified constituent radicals, a few of which are listed in the table. One constructs a reduced-order model, suitable for use in estimating the release of heat and the evolution of temperature in combustion, from a base comprising the 13 constituent radicals plus a total of 26 other species that include the light molecules and related light free radicals. Then rather than following all possible species through their reaction coordinates, one follows only the reduced set of reaction coordinates of the base. The behavior of the base was examined in test computational simulations of the combustion of heptane in a stirred reactor at various initial pressures ranging from 0.1 to 6 MPa. Most of the simulations were performed for stoichiometric mixtures; some were performed for fuel/oxygen mole ratios of 1/2 and 2.

Bellan, Josette↗
DOE OSTI2022

Approximation rates of DeepONets for learning operators arising from advection–diffusion equations

Here we present the analysis of approximation rates of operator learning in Chen and Chen (1995) and Lu et al. (2021), where continuous operators are approximated by a sum of products of branch and trunk networks. In this work, we consider the rates of learning solution operators from both linear and nonlinear advection–diffusion equations with or without reaction. We find that the approximation rates depend on the architecture of branch networks as well as the smoothness of inputs and outputs of solution operators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A theory of phonon-induced friction on molecular adsorbates

In this manuscript, we provide a general theory for how surface phonons couple to molecular adsorbates. Our theory maps the extended dynamics of a surface’s atomic vibrational motions to a generalized Langevin equation, and by doing so captures these dynamics in a single quantity: the non-Markovian friction. The different frequency components of this friction are the phonon modes of the surface slab weighted by their coupling to the adsorbate degrees of freedom. Using this formalism, we demonstrate that physisorbed species couple primarily to acoustic phonons while chemisorbed species couple to dispersionless local vibrations. We subsequently derive equations for phonon-adjusted reaction rates using transition state theory and demonstrate that these corrections improve agreement with experimental results for CO desorption rates from Pt(111).

Science & Technology - Other Topics↗
DOE OSTI2020

Neutrino signal from proto-neutron star evolution: Effects of opacities from charged-current–neutrino interactions and inverse neutron decay

In this work, we investigate the impact of charged-current–neutrino processes on the formation and evolution of neutrino spectra during the deleptonization of proto-neutron stars. To this end we develop the full kinematics of these reaction rates consistent with the nuclear equation of state, including weak magnetism contributions. This allows us to systematically study the impact of inelastic contributions and weak magnetism on the $ν_e$ and $\bar{ν}_e$ luminosities and average energies. Furthermore, we explore the role of the inverse neutron decay, also known as the direct Urca process, on the emitted spectra of $\bar{ν}_e$. This process is commonly considered in the cooling scenario of cold neutron stars but has so far been neglected in the evolution of hot proto-neutron stars. We find that the inverse neutron decay becomes the dominating opacity source for low-energy $\bar{ν}_e$. Accurate three-flavor Boltzmann neutrino transport enables us to relate the magnitude of neutrino fluxes and spectra to details of the treatment of weak processes. This allows us to quantify the corresponding impact on the conditions relevant for the nucleosynthesis in the neutrino-driven wind, which is ejected from the proto-neutron star surface during the deleptonization phase.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

An improved flux-split algorithm applied to hypersonic flows in chemical equilibrium

An explicit, finite-difference, shock-capturing numerical algorithm is presented and applied to hypersonic flows assumed to be in thermochemical equilibrium. Real-gas chemistry is either loosely coupled to the gasdynamics by way of a Gibbs free energy minimization package or fully coupled using species mass conservation equations with finite-rate chemical reactions. A scheme is developed that maintains stability in the explicit, finite-rate formulation while allowing relatively high time steps. The codes use flux vector splitting to difference the inviscid fluxes and employ real-gas corrections to viscosity and thermal conductivity. Numerical results are compared against existing ballistic range and flight data. Flows about complex geometries are also computed.

Palmer, Grant↗
DOE OSTI2020

An extended methodology for automated calculations of non-Boltzmann kinetic sequences: H + C 2 H 2 + X and combustion impact

It is generally assumed in phenomenological kinetic models that bimolecular reactions only occur between species whose rovibrational energy follows a Boltzmann (thermal) distribution. That is, any complexes initially formed in non-Boltzmann distributions are assumed to be thermalized by energy-transferring collisions prior to bimolecular reactions. Given the high mole fractions of reactive species, X, in combustion environments, reactive collisions of the complexes with X often occur on the same timescale as energy-transferring collisions – yielding sequences proceeding through non-Boltzmann intermediates across multiple potential energy surfaces. Recent studies have shown that such non-Boltzmann kinetic sequences can have substantial impact on the global reactivity in combustion systems. Simulations of these non-Boltzmann reaction sequences, which can be described in phenomenological kinetic models via chemically termolecular reactions, require that rovibrational excitation from one potential energy surface be carried over to the next. Here this paper presents an extended theoretical and computational methodology that couples multiple master equations and derives rate constants for phenomenological reactions describing the conversion of thermal reactants to thermal products for use in phenomenological kinetic schemes. The methodology is then implemented using in-house scripts for non-Boltzmann sequences involving C2H$^*_3$ + X (with X = O 2 , H, and OH) where C 2 H$^*_3$ is formed via H + C 2 H 2 association – which were identified as having strong potential for influencing combustion predictions in a recent study. The results reveal that non-Boltzmann reaction sequences for X = O 2 (the primary focus of this paper) significantly alters the total conversion rate from H + C 2 H 2 to products and product branching fractions from those of thermal sequential pathways. Furthermore, the present results demonstrate that non-Boltzmann reaction sequences have significant impact – as high as an order of magnitude – on predicted ignition delay times. Similarly, they yield significantly different dependence of ignition delay times with temperature and O 2 mole fraction – yielding signatures that are likely observable experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1976

Periodic solutions to systems of reaction-diffusion equations

Necessary and sufficient conditions are derived for the existence of temporally periodic 'dissipative structure' solutions in weak diffusion with the reaction rate terms dominant in a generic system of reaction-diffusion differential equations. The enumerator index i of the equations denotes the density or concentration of the ith participating molecular or biological species, and D sub i is the diffusivity constant for the ith species while Q sub i (c), an algebraic function of the n-tuple c, expresses the local rate of production of the ith species due to chemical reactions or biological interactions.

Rosen, G.↗
NASA NTRS1979

Analysis of turbulent free-jet hydrogen-air diffusion flames with finite chemical reaction rates

A numerical analysis is presented of the nonequilibrium flow field resulting from the turbulent mixing and combustion of an axisymmetric hydrogen jet in a supersonic parallel ambient air stream. The effective turbulent transport properties are determined by means of a two-equation model of turbulence. The finite-rate chemistry model considers eight elementary reactions among six chemical species: H, O, H2O, OH, O2 and H2. The governing set of nonlinear partial differential equations was solved by using an implicit finite-difference procedure. Radial distributions were obtained at two downstream locations for some important variables affecting the flow development, such as the turbulent kinetic energy and its dissipation rate. The results show that these variables attain their peak values on the axis of symmetry. The computed distribution of velocity, temperature, and mass fractions of the chemical species gives a complete description of the flow field. The numerical predictions were compared with two sets of experimental data. Good qualitative agreement was obtained.

Sislian, J. P.↗
DOE OSTI2020

Steady detonation propagation in thin channels with strong confinement

We examine asymptotically the dynamics of two-dimensional, steady detonation wave propagation and failure for a strongly confined high explosive (HE), in which the width of the explosive is small relative to the reaction zone length. An energy balance equation is derived, which shows how the longitudinal acceleration of subsonic flow behind the detonation shock is influenced both by chemical reaction and by the effects of HE boundary streamline deflection, specifically via the induced rate of change of mass flux through the detonation wave. The latter serves to either counteract or reinforce the acceleration of longitudinal flow, depending on the sign of the gradient of the boundary streamline deflection at the detonation shock. The analysis is valid for general equations of state and chemical reaction rates in the HE. The asymptotically derived form of the energy equation represents an eigenvalue problem for the determination of the steady detonation propagation speed, solved via a shooting method. We explore specific results for ideal and stiffened equations of state, along with a pressure-dependent reaction rate for which changes in the pressure exponent and reaction order are also studied. We consider the influences of both straight and curved HE boundary streamline shapes. The asymptotic analysis reveals significant physical insights into how detonation propagation and failure are affected by strong confinement.

42 ENGINEERING↗
DOE OSTI2024

Open Circuit Potential of a Au Catalyst during Selective Oxidation of Glycerol

It has been of great interest to understand relationships between the potential of a metal catalyst and its thermochemical catalytic activity, especially in an aqueous phase environment. In the literature, there are correlations of open circuit potential with reaction rates or surface concentrations of reaction intermediates. In this study, we measured the open circuit potential (OCP) of a Au gauze during selective oxidation of glycerol to glyceric acid and hydrogen peroxide in a basic solution, as a function of hydroxide and alcoholate concentrations. It is found that applying a potential to the Au catalyst has no influence on the reaction rate. Although a rough correlation appears to exist between OCP and reaction rate, the data are better fit to an equation which assumes that the potential of the metal (i.e. the OCP) is in equilibrium with the electrochemical potential of the solution, which is defined by the thermodynamic activities of the oxidizing and reducing species. The equation: OCP = constant + Σ ox [(RT/FZ ox )ln[Ox]] - Σ red [(RT/FZ red )ln[Red]] implicitly assumes that the Au metal functions as a probing electrode. It is further found that this equation also applies to the Au-catalyzed H 2 O 2 oxidation/decomposition reaction in a basic medium, and possibly to formic acid dehydrogenation. Here, we postulate that the apparent correlation between OCP and reaction rate is due to the fact that the reaction rate is proportional to the concentrations of the reaction products and/or reactants that define the electrochemical potential of the solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1973

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

Doppler Broadening and Other Temperature Effects

In this paper I attempt to document what I have learned and still remember about temperature effects, particularly Doppler broadening, on nuclear data, as it applies to ENDF formatted data. My focus is on the SIGMA1 method of Doppler broadening that I developed over 50 years ago. However, I do realize that there are currently many different methods used in computer codes, and I have tried to keep my discussion general as it applies to ALL of these methods. In particular, ALL of the MYTHS I describe below apply to all Doppler broadening methods: they all solve the parabolic diffusion equation in spherical geometry in (speed, reaction rate) versus temperature. These methods differ only in how the cross section is represented. So, at 85 years old here is what I still remember.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

CompLaB v1.0: a scalable pore-scale model for flow, biogeochemistry, microbial metabolism, and biofilm dynamics

Abstract. Microbial activity and chemical reactions in porous media depend on the local conditions at the pore scale and can involve complex feedback with fluid flow and mass transport. We present a modeling framework that quantitatively accounts for the interactions between the bio(geo)chemical and physical processes and that can integrate genome-scale microbial metabolic information into a dynamically changing, spatially explicit representation of environmental conditions. The model couples a lattice Boltzmann implementation of Navier–Stokes (flow) and advection–diffusion-reaction (mass conservation) equations. Reaction formulations can include both kinetic rate expressions and flux balance analysis, thereby integrating reactive transport modeling and systems biology. We also show that the use of surrogate models such as neural network representations of in silico cell models can speed up computations significantly, facilitating applications to complex environmental systems. Parallelization enables simulations that resolve heterogeneity at multiple scales, and a cellular automaton module provides additional capabilities to simulate biofilm dynamics. The code thus constitutes a platform suitable for a range of environmental, engineering and – potentially – medical applications, in particular ones that involve the simulation of microbial dynamics.

58 GEOSCIENCES↗
NASA NTRS1985

On convergence of computation of chemically reacting flows

The computational problems associated with high-temperature flows undergoing finite-rate ionization reactions is investigated. The conservation equations governing chemical species and vibrational and electron energies are solved simultaneously with those for overall mass, momentum, and energy for a one-dimensional subsonic flow, through a constant-area duct, originating behind a normal shock wave, using an implicit time-marching technique. Boundary conditions are imposed in the form of characteristic wave variables accounting for the effects of chemical reactions on the speed of sound. Converging solutions are obtained for cases in which chemical reactions are weak, but difficulty is encountered in other cases. The cause of the difficulty is investigated and shown to be the sharp pressure disturbances produced by such reactions.

Park, C.↗
NASA NTRS1978

Odd nitrogen production by meteoroids

The process by which odd nitrogen species (atomic nitrogen and nitric oxide) are formed during atmospheric entry of meteoroids is analyzed theoretically. An ablating meteoroid is assumed to be a point source of mass with a continuum regime evolving in its wake. The amounts of odd nitrogen species, produced by high-temperature reactions of air in the continuum wake, are calculated by numerical integration of chemical rate equations. Flow properties are assumed to be uniform across the wake, and 29 reactions involving five neutral species and five singly ionized species are considered, as well as vibrational and electron temperature nonequilibrium phenomena. The results, when they are summed over the observed mass, velocity, and entry-angle distribution of meteoroids, provide odd-nitrogen-species annual global production rates as functions of altitude. The peak production of nitric oxide is found to occur at an altitude of about 85 km; atomic nitrogen production peaks at about 95 km. The total annual rate for nitric oxide is 40 million kg; for atomic nitrogen it is 170 million kg.

Park, C.↗