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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2020

NRC Multiphysics Analysis Capability Deployment FY2020 - Part 4

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Sup- port for Advanced Non-Light Water Reactors.” The tasks completed for this report are: • Task 3a: This task demonstrates the effectiveness of Griffin’s methods for interpolating cross sections and SPH correction factors as a function of control drum rotation. The Monte Carlo code Serpent is used to generate reference results. The accuracy of Griffin-calculated eigenvalues and reaction rates are studied as a function of control drum rotation for cross section and SPH factor libraries of varying fidelity. Higher-fidelity libraries require less interpolation and are therefore inherently more accurate; however, they require more computational resources to create. Overall, Griffin is able to accurately model eigenvalue and reaction rates over the full range of drum rotation using a library generated from three to five discrete drum rotations for a particular microreactor design. • Task 3c: This preliminary task successfully demonstrates the performance of a coupled neutron-photon transport calculation with Griffin using the spherical harmonics approximation. The Griffin solutions are compared to reference solutions from the continuous energy Monte Carlo codes MCNP and Serpent. Furthermore, this work also demonstrates how to prepare the neutron and photon libraries using NJOY, MCNP, and Serpent with the latest ENDF/B-VIII data for 235 U and graphite. The Griffin total energy deposition in the active core region is calculated within 0.3% of the reference, but it degrades in the reflector regions with a maximum difference of -1.5%. Overall, the results look promising and future analysis should use an improved data stream and focus on determining the potential cancellation of error within each energy group.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Direct numerical simulations of three-component Rayleigh–Taylor mixing and an improved model for multicomponent reacting mixtures

We present direct numerical simulations of a three-layer Rayleigh–Taylor instability (RTI) problem with a configuration based on the experiments of Suchandra & Ranjan ( J. Fluid Mech. , vol. 974, 2023, A35) and Jacobs & Dalziel ( J. Fluid Mech. , vol. 542, 2005, pp. 251–279). The problem consists of a layer of light fluid between two layers of heavy fluid with an Atwood number of 0.3. These simulations are first validated through comparison with available experimental data. The validated simulations are then utilized to analyse statistics in this three-component flow. First, length scales are examined utilizing spectra and two-point spatial correlations of velocity and species concentration fluctuations. Next, joint probability density functions (p.d.f.s) of species concentration are compared against several model p.d.f.s representing generalizations of the bivariate beta distribution. Notably, the joint p.d.f.s do not appear to be accurately described by a Dirichlet distribution, indicating the marginal distributions do not conform to a beta distribution. Finally, similarity of the present configuration to three-component mixing found in inertial confinement fusion (ICF) applications is exploited to develop and validate an improved model for the impact of multicomponent mixing on thermonuclear (TN) reaction rates. A single time instant from the present simulations is chosen for a TN burn calculation under the hypothetical assumption of ICF materials and temperatures. Total TN output from this second calculation is then compared against the prediction of the improved model. The new model is found to accurately predict TN reaction rates in both premixed and non-premixed configurations.

42 ENGINEERING↗

NASA NTRS1983

Rate processes in gas phase

Reaction-rate theory and experiment are given a critical review from the engineers' point of view. Rates of heavy-particle, collision-induced reaction in gas phase are formulated in terms of the cross sections and activation energies for reaction. The effect of cross section function shape and of excited state contributions to reaction both cause the slope of Arrhenius plots to differ from the true activation energy, except at low temperature. The master equations for chemically reacting gases are introduced, and dissociation and ionization reactions are shown to proceed primarily from excited states about kT from the dissociation or ionization limit. Collision-induced vibration, vibration-rotation, and pure rotation transitions are treated, including three-dimensional effects and conservation of energy, which have usually been ignored. The quantum theory of transitions at potential surface crossing is derived, and results are found to be in fair agreement with experiment in spite of some questionable approximations involved.

Hansen, C. F.↗

DOE OSTI2020

Urca nuclide production in Type-I X-ray bursts and implications for nuclear physics studies

ABSTRACT The thermal structure of accreting neutron stars is affected by the presence of urca nuclei in the neutron star crust. Nuclear isobars harbouring urca nuclides can be produced in the ashes of Type I X-ray bursts, but the details of their production have not yet been explored. Using the code MESA, we investigate urca nuclide production in a one-dimensional model of Type I X-ray bursts using astrophysical conditions thought to resemble the source GS 1826-24. We find that high-mass (A ≥ 55) urca nuclei are primarily produced late in the X-ray burst, during hydrogen-burning freeze-out that corresponds to the tail of the burst light curve. The ∼0.4–0.6 GK temperature relevant for the nucleosynthesis of these urca nuclides is much lower than the ∼1 GK temperature most relevant for X-ray burst light curve impacts by nuclear reaction rates involving high-mass nuclides. The latter temperature is often assumed for nuclear physics studies. Therefore, our findings alter the excitation energy range of interest in compound nuclei for nuclear physics studies of urca nuclide production. We demonstrate that for some cases this will need to be considered in planning for nuclear physics experiments. Additionally, we show that the lower temperature range for urca nuclide production explains why variations of some nuclear reaction rates in model calculations impacts the burst light curve but not local features of the burst ashes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Learning heterogeneous reaction kinetics from X-ray videos pixel by pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries and electrocatalysts. Experimental characterizations of such materials by operando microscopy produce rich image datasets, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Verification of EvaluateFLux Utility Program

The EvaluateFlux program is a post-processing utility program for DIF3D, specifically DIF3D-VARIANT which handles Cartesian and hexagonal geometries. The EvaluateFlux program was developed to allow users to obtain flux and power traverses through the geometry domain, and its initial purpose was to facilitate foil analysis by evaluating the flux solution from DIF3D-VARIANT and combining it with foil cross section data. The EvaluateFlux program can calculate the neutron flux, as well as the reaction rates, at any user provided evaluation point. It does this by identifying the spatial mesh associated with the evaluation point and then evaluates the polynomial based neutron flux moments stored in the NHFLUX file at that point. The output of EvaluateFlux varies depending on the input setup. The maximum output includes the neutron flux and microscopic and macroscopic reaction rates at each evaluation point. The purpose of this work is to verify the outputs of EvaluateFlux. Simple models that have hand calculatable results are first defined and used to verify the EvaluateFlux outputs. More complex cases are then added where a duplicate program of EvaluateFlux that uses PrintTables outputs of the binary files is used to verify the EvaluateFlux outputs. In those complex cases, hand calculations of selected evaluation points were also displayed to confirm the software verification. For all the tests done, the hand calculations agreed well with those calculated by EvaluateFlux. For the larger complex problems, the duplicate program that can process hundreds of evaluation points was able to identify that zero points within some meshes have large errors. This aspect was attributed to the truncation error on the input provided to the duplicate program and is not a concern for the accuracy of the EvaluateFlux software.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Verification of the EvaluateFlux Utility Program

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1976

Ionospheric composition in SAR-arcs

Theoretical ion and electron density profiles in the SAR-arc region are calculated using a model of the ionosphere based on the coupled continuity, momentum, and energy equations for O(+), NO(+), and O2(+). It is found that an increase in the reaction O(+) + N2 yields NO(+) + N, which results from enhanced N2 vibrational excitation due to the high electron temperatures found in SAR arcs, can cause a reduction in F-region electron densities by up to a factor of two. The increase in the O(+) + N2 reaction rate is shown to result in a marked change in the ion composition in SAR arcs, with NO(+) being an important ion up to altitudes of about 350 km at night. Since observed electron-density depressions in SAR arcs generally vary between factors of two and seven, it is concluded that the increase in the O(+) + N2 reaction rate cannot account for these depressions by itself.

Raitt, W. J.↗

DOE OSTI2023

Learning heterogeneous reaction kinetics from X-ray movies pixel-by-pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries 1 and electrocatalysts 2. Experimental characterizations of such materials by operando microscopy produce rich image datasets 3–6, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation 7. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

Chueh, William↗

NASA NTRS2012

Performance Evaluation of Staged Bosch Process for CO2 Reduction to Produce Life Support Consumables

Utilizing carbon dioxide to produce water and hence oxygen is critical for sustained manned missions in space, and to support both NASA's cabin Atmosphere Revitalization System (ARS) and In-Situ Resource Utilization (ISRU) concepts. For long term missions beyond low Earth orbit, where resupply is significantly more difficult and costly, open loop ARS, like Sabatier, consume inputs such as hydrogen. The Bosch process, on the other hand, has the potential to achieve complete loop closure and is hence a preferred choice. However, current single stage Bosch reactor designs suffer from a large recycle penalty due to slow reaction rates and the inherent limitation in approaching thermodynamic equilibrium. Developmental efforts are seeking to improve upon the efficiency (hence reducing the recycle penalty) of current single stage Bosch reactors which employ traditional steel wool catalysts. Precision Combustion, Inc. (PCI), with support from NASA, has investigated the potential for utilizing catalysts supported over short-contact time Microlith substrates for the Bosch reaction to achieve faster reaction rates, higher conversions, and a reduced recycle flows. Proof-of-concept testing was accomplished for a staged Bosch process by splitting the chemistry in two separate reactors, first being the reverse water-gas-shift (RWGS) and the second being the carbon formation reactor (CFR) via hydrogenation and/or Boudouard. This paper presents the results from this feasibility study at various operating conditions. Additionally, results from two 70 hour durability tests for the RWGS reactor are discussed.

Vilekar, Saurabh A.↗

NASA NTRS1996

Hydrocarbon photochemistry in the upper atmosphere of Jupiter

The hydrocarbon photochemistry in the upper atmosphere of Jupiter is investigated using a one-dimensional, photochemical-diffusive, and diurnally averaged model. The important chemical cycles and pathways among the major species are outlined and a standard model for the North Equatorial Belt region is examined in detail. It is found that several traditionally dominant chemical pathways among the C and C2 species are replaced in importance by cycles involving C-C4 species. The pressure and altitude profiles of mixing ratios for several observable hydrocarbon species are compared with available ultraviolet- and infrared-derived abundances. The results of sensitivity studies on the standard model with respect to variations in eddy diffusion profile, solar flux, atomic hydrogen influx, latitude, temperature, and important chemical reaction rates are presented. Measured and calculated airglow emissions of He at 584 angstroms and H at 1216 angstroms are also used to provide some constraints on the range of model parameters. The relevance of the model results to the upcoming Galileo mission is briefly discussed. The model is subject to considerable improvement; there is a great need for laboratory measurements of basic reaction rates and photodissociation quantum yields, even for such simple species as methylacetylene and allene. Until such laboratory measurements exist there will be considerable uncertainty in the understanding of the C3 and higher hydrocarbons in the atmospheres of the jovian planets.

NASA Discipline Number 52-20↗

NASA NTRS1994

In Situ Measurements of the NO2/NO Ratio for Testing Atmospheric Photochemical Models

Simultaneous in situ measurements of NO2, NO, O3, ClO, pressure and temperature have been made for the first time, presenting a unique opportunity to test our current understanding of the photochemistry of the lower stratosphere. Data were collected from several flights of the ER-2 aircraft at mid-latitudes in May 1993 during NASA's Stratospheric Photochemistry, Aerosols and Dynamics Expedition (SPADE). The daytime ratio of NO2/NO remains fairly constant at 19 km with a typical value of 0.68 and standard deviation of +/- 0.17. The ratio observations are compared with simple steady-state calculations based on laboratory-measured reaction rates and modeled NO2 photolysis rates. At each measurement point the daytime NO2/NO with its measurement uncertainty overlap the results of steady-state calculations and associated uncertainty. However, over all the ER-2 flights examined, the model systematically overestimates the ratio by 40% on average. Possible sources of error are examined in both model and measurements. It is shown that more accurate laboratory determinations of the NO + O3 reaction rate and of the NO2 cross-sections in the 200-220 K temperature range characteristic of the lower stratosphere would allow for a more robust test of our knowledge of NO(x) photochemistry by reducing significant sources of uncertainties in the interpretation of stratospheric measurements. The present measurements are compared with earlier observations of the ratio at higher altitudes.

Jaegle, L.↗

DOE OSTI2020

Grain detachment and transport clogging during mineral dissolution in carbonate rocks with permeable grain boundaries

In this paper, we explore the effect of permeable grain boundaries on carbonate rock dissolution in order to improve our understanding of grain detachment and migration during reactive flow. To do so we investigated the effects of grain size and enhanced reactivity along grain boundaries on global dissolution kinetics. Variations of permeability and porosity were calculated, and exponential relationships were observed in fractured rocks. Our model employed a reactive transport framework based on the Darcy-Brinkman equation to simulate calcite dissolution in carbonate rocks composed of microporous grains and, included fluid transport along grain boundaries. The model includes fluid flow, solute transport by advection–diffusion, heterogeneous reaction between fluid and minerals and grain detachment with subsequent grain transport in macropores. The migration of solid particles due to dissolution was based on cluster analysis and local movement, and the results show that grain detachment can lead to significant decreases in permeability due to the clogging of transport pathways. Microporous media with smaller grain sizes (fine grains in this study) showed a higher average reaction rate than those with coarser grains. In addition, as the overall rates of geochemical processes are commonly affected by the presence of texture heterogeneities such as fractures, a single fracture (macropore) introduced into a microporous matrix composed of permeable grains and grain boundaries was modelled as a large channel connecting the inlet and outlet of the simulation domain. It was found that macropore clogging by grain detachment temporary decreased permeability, lowering the long-term global reaction rates. Lastly, we observed that increases in flow rate can reduce detachment and fracture clogging by reducing local dissolution along grain boundaries.

58 GEOSCIENCES↗

DOE OSTI2021

Impact of Flux Wire Selection on Neutron Spectrum Adjustment

In reactor dosimetry, spectrum unfolding methods are used to assess the energy-dependent neutron flux following the irradiation of flux wires or foils. For some spectrum unfolding methods, an a priori, or initial guess, spectrum is necessary. However, if the a priori spectrum is not representative of the irradiation spectrum experienced by the flux wires, the spectrum unfolding method may fail or create a nonphysical adjusted spectrum. The unfolded spectrum is dependent on flux wire reaction rates, and therefore cross section. As a result, each selection or combination of flux wires will have a different sensitivity to the group-wise fluence difference between the a priori spectrum and the irradiation spectrum. Therefore, a selection of flux wires that are insensitive to errors in the a priori spectrum would be considered robust choices for reactor dosimetry. In this work, flux wire reaction rates are calculated using a simulated neutron spectrum from the High Flux Isotope Reactor (HFIR) and spectrum unfolding methods are used in each flux wire combinations. The flux wire combinations are compared against a traditional flux wire set of Ti, Fe, and Co, with a summary table of the deemed “good” wire combinations that are insensitive to number of a priori neutron energy groups, group-wise difference between the a priori and irradiation, and magnitude of difference.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2006

Development of Comprehensive Reduced Kinetic Models for Supersonic Reacting Shear Layer Simulations

Large-scale simulations of multi-dimensional unsteady turbulent reacting flows with detailed chemistry and transport can be computationally extremely intensive even on distributed computing architectures. With the development of suitable reduced chemical kinetic models, the number of scalar variables to be integrated can be decreased, leading to a significant reduction in the computational time required for the simulation with limited loss of accuracy in the results. A general MATLAB-based automated mechanism reduction procedure is presented to reduce any complex starting mechanism (detailed or skeletal) with minimal human intervention. Based on the application of the quasi steady-state (QSS) approximation for certain chemical species and on the elimination of the fast reaction rates in the mechanism, several comprehensive reduced models, capable of handling different fuels such as C2H4, CH4 and H2, have been developed and thoroughly tested for several combustion problems (ignition, propagation and extinction) and physical conditions (reactant compositions, temperatures, and pressures). A key feature of the present reduction procedure is the explicit solution of the concentrations of the QSS species, needed for the evaluation of the elementary reaction rates. In contrast, previous approaches relied on an implicit solution due to the strong coupling between QSS species, requiring computationally expensive inner iterations. A novel algorithm, based on the definition of a QSS species coupling matrix, is presented to (i) introduce appropriate truncations to the QSS algebraic relations and (ii) identify the optimal sequence for the explicit solution of the concentration of the QSS species. With the automatic generation of the relevant source code, the resulting reduced models can be readily implemented into numerical codes.

Zambon, A. C.↗

DOE OSTI2023

Towards the first plasma-electron screening experiment

The enhancement of fusion reaction rates in a thermonuclear plasma by electron screening of the Coulomb barrier is an important plasma-nuclear effect that is present in stellar models but has not been experimentally observed. Experiments using inertial confinement fusion (ICF) implosions may provide a unique opportunity to observe this important plasma-nuclear effect. Herein, we show that experiments at the National Ignition Facility (NIF) have reached the relevant physical regime, with respect to the density and temperature conditions, but the estimated impacts of plasma screening on nuclear reaction rates are currently too small and need to be increased to lower the expected measurement uncertainty. Detailed radiation hydrodynamics simulations show that practical target changes, like adding readily available high-Z gases, and significantly slowing the inflight implosion velocity, while maintaining inflight kinetic energy, might be able to push these conditions to those where plasma screening effects may be measurable. We also perform synthetic data exercises to help understand where the anticipated experimental uncertainties will become important. But challenges remain, such as the detectability of the reaction products, non-thermal plasma effects, species separation, and impacts of spatial and temporal gradients. This work lays the foundation for future efforts to develop an important platform capable of the first plasma electron screening observation.

inertial confinement fusion↗

DOE OSTI2022

A reaction mechanism for plasma electrolysis of AgNO 3 forming silver nanoclusters and nanoparticles

In plasma-driven solution electrolysis (PDSE), gas-phase plasma-produced species interact with an electrolytic solution to produce, for example, nanoparticles. An atmospheric pressure plasma jet (APPJ) directed onto a liquid solution containing a metallic salt will promote reduction of metallic ions in solution, generating metallic clusters that nucleate to form nanoparticles. In this article, results from a computational investigation are discussed of a PDSE process in which a radio-frequency APPJ sustained in helium impinges on a silver nitrate solution, resulting in growth of silver nanoparticles. A reaction mechanism was developed and implemented in a global plasma chemistry model to predict nanoparticle growth. To develop the reaction mechanism, density functional theory was used to generate probable silver growth pathways up to Ag 9 . Neutral clusters larger than Ag 9 were classified as nanoparticles. Kinetic reaction rate coefficients for thermodynamically favorable growth pathways were estimated based on an existing, empirically determined base reaction mechanism for smaller Ag particle interactions. These rates were used in conjunction with diffusion-controlled reaction rate coefficients that were calculated for other Ag species. The role of anions in reduction of Agn ions in forming nanoparticles is also discussed. Oxygen containing impurities or admixtures to the helium, air entrainment into the APPJ, and dissociation of saturated water vapor above the solution can produce additional reactive oxygen species in solution, resulting in the production of anions and O 2 – in particular. For a given molarity, delivering a sufficient fluence of reducing species will produce similar nanoparticle densities and sizes for all applied power levels. In conclusion, comparisons are made to alternate models for nanoparticle formation, including charged nanoparticles and use of direct current plasmas.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Summary of Structural Alloy Compatibility in Supercritical CO2 at 450°-800°C

Supercritical CO2 (sCO2) is of interest as a working fluid for several concepts including the direct- fired Allam cycle as a low-emission fossil energy power cycle. Over the past 10 years, laboratory exposures at 300 bar sCO2 have found reasonably good compatibility for Ni-based alloys at <800°C, including an assessment of the sCO2 impact on room temperature mechanical properties after 750°C exposures. However, initial screening tests at 1 and 20 bar CO2 at 900°-1100°C showed poor compatibility for Ni-based alloys. In an open cycle, the introduction of 1%O2 and 0.1- 0.25%H2O impurities at 300 bar increased the reaction rates ≥2X at 750°C. At lower temperatures, steels are susceptible to C ingress and embrittlement. Creep-strength enhanced ferritic steels may be limited to <550°C and conventional stainless steels to <600°C. Two strategies to increase those temperatures are higher Ni and Cr alloying additions and Al- or Cr-rich coatings. Alloy 709 (Fe- 20Cr-25Ni) shows some promising results at 650°C in sCO2 but reaction rates were accelerated with the addition of O2 and H2O impurities. Pack aluminized and chromized Gr.91 (Fe-9Cr-1Mo) and type 316H stainless steel show some promise at 600°-650°C but further coating optimization is needed.

Pint, Bruce↗