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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Stochastic self-tuning hybrid algorithm for reaction-diffusion systems

Many biochemical phenomena involve reactants with vastly different concentrations, some of which are amenable to continuum-level descriptions, while the others are not. We present a hybrid self-tuning algorithm to model such systems. The method combines microscopic (Brownian) dynamics for diffusion with mesoscopic (Gillespie-type) methods for reactions and remains efficient in a wide range of regimes and scenarios with large variations of concentrations. Furthermore, its accuracy, robustness, and versatility are balanced by redefining propensities and optimizing the mesh size and time step. We use a bimolecular reaction to demonstrate the potential of our method in a broad spectrum of scenarios: from almost completely reaction-dominated systems to cases where reactions rarely occur or take place very slowly. The simulation results show that the number of particles present in the system does not degrade the performance of our method. This makes it an accurate and computationally efficient tool to model complex multireaction systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Skeletal Kinetics Reduction for Astrophysical Reaction Networks

A novel methodology is developed to extract accurate skeletal reaction models for nuclear combustion. Local sensitivities of isotope mass fractions with respect to reaction rates are modeled based on the forced optimally time-dependent (f-OTD) scheme. These sensitivities are then analyzed temporally to generate skeletal models. The methodology is demonstrated by conducting skeletal reduction of constant density and temperature burning of carbon and oxygen relevant to Type Ia supernovae (SNe Ia). The 495-isotopes Torch model is chosen as the detailed reaction network. A map of maximum production of 56 Ni in SNe Ia is produced for different temperatures, densities, and proton-to-neutron ratios. The f-OTD simulations and the sensitivity analyses are then performed with initial conditions from this map. A series of skeletal models are derived and their performances are assessed by comparison against currently existing skeletal models. Previous models have been constructed intuitively by assuming the dominance of α-chain reactions. The comparison of the newly generated skeletal models against previous models is based on the predicted energy release and 44 Ti and 56 Ni abundances by each model. The consequences of ye ≠ 0.5 in the initial composition are also explored where ye is the electron fraction. The simulated results show that 56 Ni production decreases by decreasing ye as expected, and that the 43 Sc is a key isotope in proton and neutron channels toward 56 Ni production. It is shown that an f-OTD skeletal model with 150 isotopes can accurately predict the 56Ni abundance in SNe Ia for ye ≲ 0.5 initial conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

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

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

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

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↗
DOE OSTI2020

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

Ensemble cure kinetics network (ECK-Net): A method to derive cure kinetics of thermosetting resin

This paper introduces an Ensemble Cure Kinetics Network (ECK-Net), a neural network (NN)–based framework for modeling the cure kinetics of thermosetting resins within a phenomenological context. ECK-Net replaces traditional analytic models, which require extensive chemical insight and multiple isothermal/non-isothermal experiments, with a data-driven surrogate that maps nonlinear relationships between temperature, degree of cure, and reaction rate from differential scanning calorimetry data. The proposed approach predicts input-dependent kinetic coefficients of a generalized nth-order reaction equation rather than reaction rates directly, enabling a single unified model to represent various epoxy systems without relying on iso-conversional analysis or predefined functional forms. To ensure robustness, multiple independently trained networks under different random initializations are blended through an ensemble strategy, effectively mitigating the stochastic variability inherent to neural networks. The framework is validated using experimental datasets from multiple resin systems, including aerospace-grade materials (Toray 3900-2, Cycom 5320-1, and Hexcel 8552) and a windmill-grade resin (RIMR 035c). The model accurately reproduces the temporal evolution of the degree of cure under manufacturers’ recommended cure cycles across all tested resins systems, yielding Pearson’s correlation coefficients of 0.992, 0.994, 0.993, 0.997, respectively. To demonstrate process-level applicability, the trained network was implemented within the Abaqus environment to simulate out-of-autoclave (OOA) curing process of the CFRP panel composed of Toray T830H-6K/3900-2D prepreg. The simulation results showed excellent agreement with experimental temperature response (maximum peak temperature, simulation: 189.6 °C, experiment: 188.5 °C) and the final degree of cure (simulation: 0.948, experiment: 0.960 ± 0.013), confirming ECK-Net’s capability as a reliable alternative to conventional cure kinetics modeling methods.

Composite curing↗
NASA NTRS1982

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

PolyODENet: Deriving mass-action rate equations from incomplete transient kinetics data

Kinetics of a reaction network that follows mass-action rate laws can be described with a system of ordinary differential equations (ODEs) with polynomial right-hand side. However, it is challenging to derive such kinetic differential equations from transient kinetic data without knowing the reaction network, especially when the data are incomplete due to experimental limitations. We introduce a program, PolyODENet, toward this goal. Based on the machine-learning method Neural ODE, PolyODENet defines a generative model and predicts concentrations at arbitrary time. As such, it is possible to include unmeasurable intermediate species in the kinetic equations. Importantly, we have implemented various measures to apply physical constraints and chemical knowledge in the training to regularize the solution space. Here, using simple catalytic reaction models, we demonstrate that PolyODENet can predict reaction profiles of unknown species and doing so even reveal hidden parts of reaction mechanisms.

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

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

Suppression of Composition g -modes in Chemically Equilibrating Warm Neutron Stars

We investigate the impact of chemical equilibration and the resulting bulk viscosity on nonradial oscillation modes of warm neutron stars at temperatures up to T ≈ 5 MeV, relevant for protoneutron stars and neutron star postmerger remnants. In this regime, the relaxation rate of weak interactions becomes comparable to the characteristic frequencies of composition g-modes in the core, resulting in resonant damping. To capture this effect, we introduce the dynamical sound speed, a complex, frequency-dependent generalization of the adiabatic sound speed that encodes both the restoring force and the dissipative effects of bulk compression. Using realistic weak reaction rates and three representative equations of state, we compute the complex frequencies of composition g-modes with finite-temperature profiles. We find that bulk viscous damping becomes increasingly significant with temperature and can completely suppress composition g-modes. In contrast, the f-mode remains largely unaffected by bulk viscosity due to its nearly divergence-free character. Our results highlight the sensitivity of g-mode behavior to thermal structure, weak reaction rates, and the equation of state, and establish the dynamical sound speed as a valuable descriptor characterizing oscillation properties in dissipative neutron star matter.

Neutron stars↗
NASA NTRS2013

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↗
DOE OSTI2021

Quarkonium semiclassical transport in quark-gluon plasma: factorization and quantum correction

We study quarkonium transport in the quark-gluon plasma by using the potential nonrelativistic QCD (pNRQCD) effective field theory and the framework of open quantum systems. We argue that the coupling between quarkonium and the thermal bath is weak using separation of scales, so the initial density matrix of the total system factorizes and the time evolution of the subsystem is Markovian. We derive the semiclassical Boltzmann equation for quarkonium by applying a Wigner transform to the Lindblad equation and carrying out a semiclassical expansion. We resum relevant interactions to all orders in the coupling constant at leading power of the nonrelativistic and multipole expansions. The derivation is valid for both weakly coupled and strongly coupled quark-gluon plasmas. We find reaction rates in the transport equation factorize into a quarkonium dipole transition function and a chromoelectric gluon distribution function. For the differential reaction rate, the definition of the momentum dependent chromoelectric gluon distribution function involves staple-shaped Wilson lines. For the inclusive reaction rate, the Wilson lines collapse into a straight line along the real time axis and the distribution becomes momentum independent. The relation between the two Wilson lines is analogous to the relation between the Wilson lines appearing in the gluon parton distribution function (PDF) and the gluon transverse momentum dependent parton distribution function (TMDPDF). The centrality dependence of the quarkonium nuclear modification factor measured by experiments probes the momentum independent distribution while the transverse momentum dependence and measurements of the azimuthal angular anisotropy may be able to probe the momentum dependent one. We discuss one way to indirectly constrain the quarkonium in-medium real potential by using the factorization formula and lattice calculations. The leading quantum correction to the semiclassical transport equation of quarkonium is also worked out. The study can be easily generalized to quarkonium transport in cold nuclear matter, which is relevant for quarkonium production in eA collisions in the future Electron-Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The Reaction of o-Benzyne with Vinylacetylene: An Unexplored Way to Produce Naphthalene

Herein, the mechanisms and kinetics of the reaction of ortho-benzyne with vinylactylene have been studied by ab initio and density functional CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G(d,p) calculations of the pertinent potential energy surface combined with Rice-Ramsperger-Kassel-Marcus - Master Equation calculations of reaction rate constants at various temperatures and pressures. Under prevailing combustion conditions, the reaction has been shown to predominantly proceed by the biradical acetylenic mechanism initiated by the addition of C4H4 to one of the C atoms of the triple bond in ortho-benzyne by the acetylenic end, with a significant contribution of the concerted addition mechanism. Following the initial reaction steps, an extra six-membered ring is produced and the rearrangement of H atoms in this new ring leads to the formation of naphthalene, which can further dissociate to 1- or 2-naphthyl radicals. The o-C 6 H 4 + C 4 H 4 reaction is highly exothermic, by ~143 kcal/mol to form naphthalene and by 31-32 kcal/mol to produce naphthyl radicals + H, but features relatively high entrance barriers of 9-11 kcal/mol. Although the reaction is rather slow, much slower than the reaction of phenyl radical with vinylacetylene, it forms naphthalene and 1- and 2-naphthyl radicals directly, with their relative yields controlled by the temperature and pressure, and thus represents a viable source of the naphthalene core under conditions where ortho-benzyne and vinylacetylene are available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

On the mechanism of soot nucleation. II. E-bridge formation at the PAH bay

A recently proposed mechanism of soot nucleation based upon the formation of a rotationally-activated dimer in the collision of an aromatic molecule and a radical leading to a stable, doubly-bonded E-bridge between them, rooted in the existence of a five-membered ring on the molecule edge, has been further investigated with a focus on the 5-6 E-bridge forming in the reaction of a PAH cyclopenta group with a bay site of another PAH. As a prototype reaction of this kind, we examined the reaction between 4-phenanthrenyl and acenaphthylene and, to project these results to larger aromatic structures, we also explored the size effect of the E-bridge forming reactions by computing the 1-naphthyl + acenaphthylene system and comparing these results with the prior results for pyrenyl + acepyrene. The two systems have been studied through high-level G3(MP2,CC)//B3LYP/6-311G(d,p) ab initio calculations of their potential energy surfaces combined with the RRKM/Master Equation calculations of reaction rate constants. With PAH monomers of similar sizes involved, the formation of E-bridge structures at the bay radical sites appeared to be faster and lead to increased nucleation rates as compared to the zigzag-forming ones. Furthermore, a model combining both the bay and zigzag rotationally-induced formation of E-bridges successfully reaches the level of nucleation fluxes comparable to those of the irreversible pyrene dimerization, thus affirming the rotationally-activated dimerization as a feasible mechanism for soot particle nucleation.

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

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↗
NASA NTRS1993

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