Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reaction rate constants”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electrolytes in Multiple-Phase Hydrogen Storage Reactions

Multiple-phase hydrogen storage materials such as metal alanates and borohydrides, and destabilized systems offer the possibility of high hydrogen storage capacity with favorable thermodynamics. However, the multiphase nature of these materials intrinsically limits the kinetics due to the required transport of species between phases, which are typically in dry powder form. To address this limitation, the influence of added electrolytes is explored. This approach is motivated by analogy with similar multiphase battery reactions that show reduced kinetic limitations while necessarily containing electrolytes. Previous experimental results showing improved kinetics for MgH2/Sn (using a LiBH4/KBH4 eutectic electrolyte) and NaAlH4 (using a diglyme electrolyte) are further analyzed in terms of this analogy. The results show that the analogy is useful and rate constants are increased. Importantly, the inclusion of an electrolyte also appears to alleviate the continuously decreasing rates with the extent of reaction, which is characteristic of many multiphase hydrides. Instead, reaction rates are approximately constant until near completion. Together, these effects can lead to >10× shorter overall reaction times. In addition, new results are presented for the hydrogenation of MgB2 using Li/K/CsI and Li/K/CsCl eutectic electrolytes, where >60% conversion to Mg(BH4)2 was demonstrated at 350 bar.

08 HYDROGEN↗
DOE OSTI2020

Water Catalysis of the Reaction of Methanol with OH Radical in the Atmosphere is Negligible

Abstract Faced with the contradictory results of two recent experimental studies [Jara‐Toro et al., Angew. Chem. Int. Ed . 2017 , 56 , 2166 and Chao et al., Angew. Chem. Int. Ed . 2019 , 58 , 5013] of the possible catalytic effect of water vapor on CH 3 OH + OH reaction, we report calculations that corroborate the conclusion made by Chao et al. and extend the rate constant evaluation down to 200 K. The rate constants of the CH 3 OH + OH reaction catalyzed by a water molecule are computed as functions of temperature and relative humidity using high‐level electronic structure and kinetics calculations. The Wuhan–Minnesota Scaling (WMS) method is used to provide accurate energetics to benchmark a density functional for direct dynamics. Both high‐frequency and low‐frequency anharmonicities are included. Variational and tunneling effects are treated by canonical variational transition state theory with multidimensional small‐curvature tunneling. And, most significantly, we include multistructural effects in the rate constant calculations. Our calculations show that the catalytic effect of water vapor is not observable at 200–400 K.

Wu, Junjun↗
DOE OSTI2021

Probing the antiknock effect of anisole through an ignition, speciation and modeling study of its blends with isooctane

In order to unravel the reaction pathways relevant to anisole co-oxidation within a fuel blend, a detailed study of isooctane/anisole blends was performed with the ULille RCM. Ignition delays as well as mole fraction profiles were measured during a two-stage ignition delay using sampling and GC techniques. These results are used to validate a kinetic model developed from ab initio calculations for the most relevant rate constants which included H-atom abstraction reactions from anisole, and reactions on the potential energy surfaces of methoxyphenyl + O 2 and anisyl + O 2 . Pressure dependent rate constants were computed for the methoxyphenyl + O 2 and anisyl + O 2 reactive systems using master equation code analysis. The new kinetic model shows good agreement with the experimental data. Dual brute-force sensitivity analysis was performed, on both first- and second-stages of ignition, allowing the identification of the most important reactions in the prediction of both ignition delays. It was observed that while pure anisole does not show NTC behavior, a 60/40 isooctane/anisole blend displays such behavior, as well as two-stage ignition. This suggests anisole addition may not be as beneficial to knock resistance as expected from its high octane number. The kinetic modeling results demonstrate the importance of H-abstraction reactions both from the methoxy group and from the aryl ring in ortho-position and the addition of the resultant radicals to O2, mostly leading to the formation of polar or non-aromatic products.

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

Development of an Enzyme-Inhibitor Reaction Using Cellular Retinoic Acid Binding Protein II for One-Pot Megamolecule Assembly

This paper presents an enzyme building block for the assembly of megamolecules. The system is based on the inhibition of the human-derived cellular retinoic acid binding protein II (CRABP2) domain. We synthesized a synthetic retinoid bearing an arylfluorosulfate group, which uses sulfur fluoride exchange click chemistry to covalently inhibit CRABP2. We conjugated both the inhibitor and a fluorescein tag to an oligo(ethylene glycol) backbone and measured a second-order rate constant for the protein inhibition reaction of approximately 3,600 M -1 s -1 . We used this new enzyme-inhibitor pair to assemble multi-protein structures in one-pot reactions using three orthogonal assembly chemistries to demonstrate exact control over the placement of protein domains within a single, homogeneous molecule. This work enables a new dimension of control over specificity, orientation, and stoichiometry of protein domains within atomically precise nanostructures.

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

Solvent effects on the heterogeneous growth of TiO 2 nanostructure arrays by solvothermal synthesis

One-dimensional titanium dioxide (TiO 2 ) nanostructure arrays (nanoarrays) are important metal oxide nanomaterials that can be synthesized via facile solvothermal methods. The properties of the organic solvents can impose significant effects on the microstructures and properties of the final products. However, the discussions were limited to the homogeneously nucleated TiO 2 nanomaterials in free-standing powder form, while the solvent effects are less understood during the heterogeneous growth of TiO 2 nanoarrays on a substrate surface. In this work, six organic compounds, namely 2-butanone, n-decane, n-hexane, toluene, ethylene glycol and ethanol, were selected as the solvents for the solvothermal synthesis of TiO 2 nanoarrays on the cordierite monolithic substrates. Special attentions are paid to the morphology, crystallinity, specific surface area, and porosity of the samples. The heterogeneous growth of TiO 2 nanoarrays on substrate surfaces is found to favor the solvents with moderate dielectric constants, which can be partially dissolved in the aqueous solution and modulate the reaction rate during the solvothermal synthesis. Organic solvents with low dielectric constants may result in a complete separation between the precursors and aqueous solution, and therefore slow down the overall reaction, causing the insufficient growth of the nanoarrays. The TiO 2 nanoarrays are obtained with optimum morphology from the combination of 2-butanone and titanium (IV) butoxide as solvent and precusor, respectively, with a high specific surface area up to 56 m 2 /g including cordierite substrate given a micron thickness. When loaded with Pt catalyst, the TiO 2 nanoarray-based monolithic catalysts show excellent low-temperature catalytic activity and hydrothermal stability for the CO and hydrocarbon oxidation under the simulated exhausted conditions. This article shall shed light on a better understanding of the growth mechanism and rational design of TiO 2 nanoarrays for high-performance catalytic converters.

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

Quantum and anharmonic effects in non-adiabatic transition state theory

Quantitative descriptions of non-adiabatic transition rates at intermediate temperatures are challenging due to the simultaneous importance of quantum and anharmonic effects. Here, in this paper, the interplay between quantum effects—for motion across or along the seam of crossing—and anharmonicity in the seam potential is considered within the weak coupling limit. The well-known expression for quantized 1-D motion across the seam (i.e., tunneling) in the linear terms approximation is derived in the thermal domain using the Lagrangian formalism, which is then applied to the case when tunneling is distributed along the seam of crossing (treating motion along the seam classically). For high-frequency quantum modes, a vibrationally adiabatic (VA) approach is developed that introduces to the non-adiabatic rate constant a factor associated with high-frequency wavefunction overlap; this approach treats the high-frequency motion along the seam quantum mechanically. To test these methodologies, the reaction N 2 O ↔ N 2 + O( 3 P) was chosen. CCSD(T)-F12b/cc-pVTZ-F12 explorations of the 3 A'- 1 A' seam of N 2 O revealed that seam anharmonicity has a strong effect on the rate constant (a factor of ~20 at 2000 K). Several quantum effects were found to be significant at intermediate/lower temperatures, including the quantum N–N vibration that was coupled with seam anharmonicity using the VA approach. Finally, a 1-D approximation to non-adiabatic instanton theory is presented to estimate the validity limit of the linear terms model at low temperatures (⁠~ 250 K for N 2 O). We recommend that the assumptions built into many statistical theories for non-adiabatic reactions—harmonic behavior, classical motion, linear terms, and weak coupling—should be verified on a case-by-case basis.

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

Computing Reaction Kinetics with MC-PDFT–OPESf: Combining Multireference Electronic Structure Theory and Enhanced Sampling

Accurate rate constants are crucial for understanding and optimizing catalytic reactions mediated by enzymes, metalloproteins, and heterogeneous catalysts. These systems frequently present a dual computational challenge. Multiconfigurational reaction sites require multireference techniques for the accurate treatment of the electronic structure, and high activation barriers prevent efficient sampling of unbiased reactive transitions. In this work, we combine multiconfiguration pair-density functional theory (MC-PDFT) as an accurate and efficient multireference electronic structure method with on-the-fly probability-enhanced sampling flooding (OPESf) as an enhanced sampling method capable of accelerating reactive transitions. We demonstrate the approach on the Diels–Alder [4+2] cycloaddition between cis-butadiene and ethene as a reaction characterized by a large activation barrier and multireference character. MC-PDFT–OPESf provides reaction rates in agreement with experiments at a fraction of the computational cost required by conventional unbiased ab initio calculations. Here, we propose MC-PDFT–OPESf as an efficient approach for computing kinetics in strongly correlated molecular systems.

Chemical calculations↗
DOE OSTI2025

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The role of energy transfer and competing bimolecular reactions in characterizing the unimolecular dissociations of allylic radicals

Resonantly stabilized radicals (RSR's) such as allylic radicals typically have bond energies exceeding 45 kcal/mol and are therefore reasonably long-lived up to quite high temperatures. Consequently, kinetics of the unimolecular dissociations and bimolecular reactions of these long-lived RSR's assume significance in high temperature chemically reacting systems such as in flames. Our recent analyses (J. Cho et al., Proc. Comb. Inst., 2022) of the uncertainty propagated by the energy transfer parameters for 1-methylallyl (1MA) dissociation to the laminar flame speed, indicates an intricate coupling between the kinetics of 1MA dissociation (chain propagation) and its reaction with H-atoms (chain-termination). Here, this work extends such analyses to other C 3 – C 5 allylic radicals: allyl, 2-methylallyl, 1,1-dimethylallyl, 1,2-dimethylallyl, and 1,3-dimethylallyl. Theoretical kinetics for the C 3 – C 4 allylic radicals were taken from prior literature studies (J.A. Miller et al., J. Phys. Chem., 2008; J. Cho et al., Proc. Comb. Inst., 2022; R.S. Tranter et al., Proc. Comb. Inst., 2017). Potential energy surfaces were newly computed for the dissociations of the larger C5 allylic radicals and the reactions of these radicals with H-atoms. These PES's were then used to calculate rate constants for the relevant unimolecular and bimolecular reactions. Energy transfer has a crucial role to play in these pressure-dependent reactions and therefore the present kinetics calculations relied on a-priori predictions of collisional energy transfer parameters to best characterize the branching between chain propagating unimolecular dissociations of these allylic radicals and chain terminating bimolecular reactions of these RSR's with H-atoms. Prompt dissociations are now established to be a universal feature of free radical dissociations. Here we have also characterized the propensity for prompt dissociations in these more-stable allylic RSR's and their role in combustion simulations.

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

An experimental and kinetic modeling study of the pyrolysis of isoprene, a significant biogenic hydrocarbon in naturally occurring vegetation fires

Isoprene dominates the carbon flux emitted by vegetation and constitutes 40% of non-methane biogenic emissions worldwide. Despite pyrolysis experiments at temperatures above 1000 K showing a link between isoprene combustion and aromatic species formation, comprehensive mechanistic research on isoprene is scarce in the literature. Here, in this work, we carry out an experimental and theoretical study to build, for the first time, a chemical kinetic model describing isoprene pyrolysis. The formation of polycyclic aromatic hydrocarbon (PAH) precursor species, often observed in vegetation fire plumes, is partially explained by isoprene pyrolysis experiments and theoretical modeling. Molecular dynamics (MD) simulations unveil reaction pathways from allylic isoprenyl radicals to allene and cyclopentadiene (CPD) intermediates, two relevant species detected in the experiments. Rate constants for these identified pathways are calculated using variational transition state theory to update the kinetic model, which is validated against single-pulse shock tube (SPST), and jet-stirred reactor (JSR) experimental data in the temperature range of 850–1690 K. The kinetic model presents satisfactory agreement with the SPST experimental data, and a reaction pathway analysis shows that association of propargyl radicals results in benzene formation. The JSR pathway analysis also identifies the prominent reactions for CPD, benzene, styrene, and toluene formation. Our model does not reproduce the CPD experimental profiles, indicating that additional studies are necessary. Overall, our findings advance the understanding of isoprene pyrolysis and its related atmospheric pollutants in naturally occurring vegetation fires where smoldering and oxygen-deficient combustion processes are present.

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

Bipolar Organic Cathodes for Stable and Sustainable Na‐Ion and Rechargeable Al Batteries

Bipolar organic materials have emerged as promising cathode materials for rechargeable batteries because of their high voltage and high capacity. However, they suffer from poor cyclic stability and slow reaction kinetics. In this work, we designed and synthesized two bipolar organic cathode materials, containing carbonyl (n‐type) and amine (p‐type) functional groups, as well as extended conjugation structures, for Na‐ion batteries (NIBs) and rechargeable aluminum batteries (RABs). As universal electrode materials, bipolar organic materials exhibited exceptional electrochemical performance in terms of high capacity, high voltage, long cycle life, and fast rate capability. Further, the extended conjugation structures in backbones of the bipolar organic materials facilitate the π–π stacking with graphene, playing a critical role in the high performance. Furthermore, the formation of a stable and robust NaF‐rich cathode electrolyte interphase was shown to stabilize the bipolar organic cathode in NIBs. Electrochemical kinetic measurements reveal that both functional groups undergo reversible redox reactions. Specifically, the electron transfer rate constant of the p‐type amine group is one order of magnitude higher than that of the n‐type carbonyl group. These results highlight the efficacy of developing bipolar organic materials for achieving high‐performance organic cathode in NIBs and RABs.

25 ENERGY STORAGE↗
DOE OSTI2022

Atomistic thermodynamics and kinetics of dicalcium silicate dissolution

Low-temperature cement manufacturing has garnered academic and industrial attention for its low environmental footprints. However, the sluggish hydration kinetics of the resultant cement affects their early-age strength development. This motivates fundamental studies to unravel the mechanistic picture of the dissolution process and discover science-informed pathways to accelerate hydration. Standard atomistic simulations seldomly exceed a microsecond making them impractical to study slow dissolution processes. We are using rare event sampling techniques, we provide the mechanistic picture of Ca 2+ ion dissolution from a kink site on the dicalcium silicate surface. The Ca 2+ ion dissolution is comprised of two sequential stages: breaking restraints from the kink sites to form a ledge adatom and detaching from the ledge/terrace adatom sites into the solution. The first and second stages feature free energy barriers of ~63 kJ/mol and ~ 29 kJ/mol respectively, making the first stage the rate-limiting step of the entire Ca 2+ dissolution kinetics. Using the reactive flux method, the rate and equilibrium constants for each reaction step are calculated, which yield the Ca 2+ ion activity of ~1.03 × 10 -5 . The diffusion calculations indicate that the surface effects lower the self-diffusion coefficient of Ca 2+ ions at the solid-water interface.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development of a kinetic-thermodynamic model for lime-stabilization of Na-bentonite

This study presents the first kinetic model to predict the solid and pore solution composition of Na-bentonite clay reacting with slaked lime over a period of 720 days. The model successfully accounts for most experimental data using a single kinetic rate constant. The following sequence of reactions was predicted by the model: initial rapid dissolution of portlandite within the first 7 days, leading to a decrease in pH and dissolved calcium, and concurrent formation of calcium silicate hydrates (C-S-H: jennite), calcium aluminate hydrate (C-A-H: C₄AH₁₃), calcium aluminosilicate hydrates (stratlingite) and hydrotalcite. After 7 days, jennite and stratlingite are predicted to transform into tobermorite-II, contributing to strength development up to 28 days. From 28 to 90 days, continued montmorillonite dissolution is predicted, along with minor formation of ettringite, partial tobermorite-II dissolution, and precipitation of secondary phases such as albite and talc. Experimentally, portlandite dissolution was confirmed by TGA and XRD and found to be complete within 7 days, in agreement with model predictions. However, other predicted solid-phase transformations (e.g., tobermorite-II formation and dissolution, ettringite, albite, and talc formation) could not be conclusively verified through experimental techniques. Aqueous phase measurements confirmed that the pH and Ca trends in solution, and that equilibrium was reached by 90 days.

Chemical kinetics↗
DOE OSTI2023

Proton transfer kinetics of transition metal hydride complexes and implications for fuel-forming reactions

Proton transfer reactions involving transition metal hydride complexes are prevalent in a number of catalytic fuel-forming reactions, where the proton transfer kinetics to or from the metal center can have significant impacts on the efficiency, selectivity, and stability associated with the catalytic cycle. Here, this review correlates the often slow proton transfer rate constants of transition metal hydride complexes to their electronic and structural descriptors and provides perspective on how to exploit these parameters to control proton transfer kinetics to and from the metal center. A toolbox of techniques for experimental determination of proton transfer rate constants is discussed, and case studies where proton transfer rate constant determination informs fuel-forming reactions are highlighted. Opportunities for extending proton transfer kinetic measurements to additional systems are presented, and the importance of synergizing the thermodynamics and kinetics of proton transfer involving transition metal hydride complexes is emphasized.

08 HYDROGEN↗
DOE OSTI2020

Theoretical study of the reaction mechanism and kinetics of the phenyl + propargyl association

Potential energy surface for the phenyl + propargyl radical recombination reaction has been studied at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G** level of theory for the closed-shell singlet species and at the triplet–singlet gap CASPT2/cc-pVTZ-CCSD(T)-F12/cc-pVTZ-f12//CASSCF/cc-pVTZ level of theory for the diradical species. High-pressure limit rate constants for the barrierless channels were evaluated with variable reaction coordinate transition state theory (VRC-TST). Rice–Ramsperger–Kassel–Marcus Master Equation (RRKM-ME) calculations have been performed to assess temperature- and pressure-dependent phenomenological rate constants and product branching ratios. Here, the entrance channels of the radical association reaction produce 3-phenyl-1-propyne and phenylallene which can further dissociate/isomerize into a variety of unimolecular and bimolecular products. Theoretical evidence is presented that, at combustion relevant conditions, the phenyl + propargyl recombination provides a feasible mechanism for the addition of a second five-member ring to the first six-member aromatic ring producing the prototype two-ring species indene and indenyl. Rate expressions for all important reaction channels in a broad range of temperatures and pressures have been generated for kinetic modeling.

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

Substitution Reactions in the Pyrolysis of Acetone Revealed through a Modeling, Experiment, Theory Paradigm

The development of high-fidelity mechanisms for chemically reactive systems is a challenging process that requires the compilation of rate descriptions for a large and somewhat ill-defined set of reactions. The present unified combination of modeling, experiment, and theory provides a paradigm for improving such mechanism development efforts. Here we combine broadband rotational spectroscopy with detailed chemical modeling based on rate constants obtained from automated ab initio transition state theory-based master equation calculations and high-level thermochemical parametrizations. Broadband rotational spectroscopy offers quantitative and isomer-specific detection by which branching ratios of polar reaction products may be obtained. Using this technique, we observe and characterize products arising from H atom substitution reactions in the flash pyrolysis of acetone (CH 3 C(O)CH 3 ) at a nominal temperature of 1800 K. The major product observed is ketene (CH 2 CO). Minor products identified include acetaldehyde (CH 3 CHO), propyne (CH 3 CCH), propene (CH 2 CHCH 3 ), and water (HDO). Literature mechanisms for the pyrolysis of acetone do not adequately describe the minor products. The inclusion of a variety of substitution reactions, with rate constants and thermochemistry obtained from automated ab initio kinetics predictions and Active Thermochemical Tables analyses, demonstrates an important role for such processes. Furthermore, the pathway to acetaldehyde is shown to be a direct result of substitution of acetone's methyl group by a free H atom, while propene formation arises from OH substitution in the enol form of acetone by a free H atom.

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

Electrochemical Corrosion of SIMFUEL: Effects of dissolved H2 and noble metal particles

This report describes the results from electrochemical corrosion tests that were conducted using simulated spent fuel (SIMFUEL) materials comprised of UO2 and surrogate fission products. Two SIMFUEL compositions were tested to quantify the effect of noble metal inclusions and dissolved H2 on the UO2 dissolution rate. One material consisted of UO2 with added lanthanide oxides (UO2 N) and the other consisted of UO2 with added lanthanide oxides and noble metals (UO2-H) at concentrations to simulate high burnup fuel. The electrochemical corrosion tests on the SIMFUEL materials were conducted in aqueous electrolyte solutions (pH 10) that were either saturated with air or purged with a H2/Ar gas mixture to maintain a constant dissolved H2 concentration. The results from electrochemical corrosion tests on SIMFUEL can be applied to qualitatively understand and also quantify the separate effects of water chemistry and fuel composition on the degradation behavior of the UO2 matrix. Open circuit potential (OCP) measurements on a SIMFUEL material of known composition (i.e., known fraction of NMPs at the fuel surface) immersed in a known solution chemistry (i.e., pH, Eh, [O2], [H2]) at known temperature provide qualitative insight into the degradation behavior of the UO2 matrix. If the OCP is above the threshold potential at which the oxidative dissolution of U(IV) to U(VI) occurs, the SNF is expected to degrade by oxidative processes under the experimental conditions. The net currents that are measured during potentiostatic tests, during which the surface potential of the SIMFUEL material is fixed by a potentiostat, can be used to quantify the UO2 degradation rate and optimize the rate constant values used in the Fuel Matrix Degradation Model (FMDM) for half reactions that occur on the SNF surface. Specifically, electrochemical measurements enable the estimation of the total anodic current at ECORR–which is the surface potential at which the total anodic and total cathodic currents are equal–so that the rate constant values for key reactions can be calculated. The open circuit potential measurements, potentiostatic tests, surface property measurements (scanning electron microscopy images and electrochemical impedence spectroscopy plots), and solution elemental composition analyses are being performed to update and optimize the FMDM. A case study is presented herein to show how electrochemical corrosion test results and accompanying characterization results for the UO2 N material in air-saturated solution can be used to validate the SNF surface reaction module of the FMDM. Future electrochemical tests will be conducted to provide quantitative information on the effects of NM content (burnup), H2 concentration, water chemistry, temperature, and galvanic couples with cladding and EBS alloys on UO2 degradation that can be used to improve the accuracy and functionality of the FMDM.

Thomas, Sara↗