Search NASA⌕ Search

SEARCH · Search NASA

Results for “Arrhenius equation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Unpacking model inadequacy: The quantification of silver release from TRISO fuel by considering empirical and mechanistic approaches

Increasing adoption of the proposed tristructural isotropic (TRISO) particle fuel for both advanced and existing reactors makes it critical to assess and address any uncertainties and inadequacies of TRISO fission product release models. Model inadequacy stems from simplifications made to the computational model when compared to the experiments. The modeling and simulation efforts conducted using the BISON fuel performance code, along with the experimental campaigns carried out under the Advanced Gas Reactor Fuel Development and Qualification Program, afford a unique opportunity to conduct a rigorous modeling inadequacy assessment within the Bayesian uncertainty quantification (UQ) framework. Here, this study compares the standard Bayesian framework against the Kennedy-O'Hagan (KOH) framework, which explicitly represents modeling inadequacy, in regard to UQ for TRISO silver release models. For this purpose, both the traditional Arrhenius equation fitted to experimental data and the more advanced lower-length-scale (LLS)-informed model, which considers microstructure information, are independently considered. Applying the inverse UQ process on the AGR-2 and -3/4 datasets revealed modeling inadequacy to be the most dominant source of uncertainty. Experimental noise uncertainty is also significant; however, model parameter uncertainty can be considered negligible. Interestingly, both the Arrhenius equation and the LLS-informed model demonstrated similar levels of modeling inadequacy. For the forward predictive UQ, the KOH framework improved both the accuracy and quality of quantified uncertainties in comparison to the standard Bayesian framework. This is true for both the Arrhenius equation and the LLS-informed model. In comparing these modeling approaches, both demonstrated similar performance at the engineering scale, while the LLS-informed model expectedly outperformed the Arrhenius equation at the mesoscale. These conclusions highlight the importance of explicitly accounting for modeling inadequacy in the UQ process, and reinforce the need for continuous refinement of physics-based models in order to address the modeling inadequacy.

Advanced reactors↗
DOE OSTI2025

Understanding Model Inadequacy in TRISO Nuclear Fuel Fission Products Release Models: Empirical and Mechanistic Approaches

The increasing use of tristructural isotropic (TRISO) particle fuel in both advanced and existing reactors necessitates a thorough evaluation of uncertainties and shortcomings in TRISO fission product release models. These inadequacies arise from the simplifications made in computational models compared to experimental data. Utilizing the BISON fuel performance code and experimental data from the Advanced Gas Reactor (AGR) program provides a unique chance to rigorously assess these inadequacies within a Bayesian uncertainty quantification (UQ) framework. This study contrasts the standard Bayesian framework with the Kennedy-O'Hagan (KOH) framework, which explicitly accounts for modeling inadequacies, in the context of UQ for TRISO silver release models. It examines both the traditional Arrhenius equation and a more advanced lower-length-scale (LLS)-informed model that incorporates microstructure information. The inverse UQ process applied to AGR-2 and AGR-3/4 datasets identified modeling inadequacy as the primary source of uncertainty, with experimental noise also being significant, while model parameter uncertainty was minimal. Both the Arrhenius and LLS-informed models showed similar levels of modeling inadequacy. For forward predictive UQ using the AGR-1 dataset, the KOH framework enhanced the accuracy and quality of quantified uncertainties by approximately 30% and 40%, respectively, compared to the standard Bayesian framework. This improvement was observed for both the Arrhenius and LLS-informed models. At the engineering scale, both models performed similarly, but the LLS-informed model outperformed the Arrhenius equation at the mesoscale. These findings underscore the importance of explicitly considering modeling inadequacy in the UQ process and highlight the need for ongoing refinement of physics-based models to address these shortcomings.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Water and Solute Activities Regulate CO 2 Reduction in Gas-Diffusion Electrodes

Electrolysis of CO 2 at gas-diffusion electrodes (GDEs) has typically been limited by the supply of gas to the electrocatalyst, overshadowing the importance of the supply of water. However, at high current densities that approach 1 A cm –2 , where the electrolyte becomes highly concentrated in the catalyst layer of a GDE, the activity of water and solutes deviate from their bulk dilute solution values, potentially slowing reaction rates and changing reaction equilibrium potentials. In addition, as flow plates for the gas stream are introduced to enable larger electrodes and high single pass conversion of CO 2 to product, variations in the gas composition will become important. By drawing upon literature for the oxygen reduction reaction (ORR), here we explain how to account for these effects in future modeling and experimental work, with particular attention to accurate use of the Nernst equation for electrode potentials and the Arrhenius equation for reaction rates. Specifically, using measurements of KOH solvent and solute activity reported in literature, and assuming the second protonation of CO 2 by water as the rate-determining step, we show the Nernst equation dilute-solution approximation of the CO 2 to CO equilibrium potential to be accurate below 5 M KOH, but it has a 74 mV error when increasing the concentration up to 10 M KOH. Finally, a simple one-dimensional model of a serpentine flow-field on a GDE demonstrated that a reactor with constant pressure of 1 bar and 1 A cm –2 at the inlet had only ~0.3 A cm –2 at the outlet for a conversion in CO 2 partial pressure from 0.90 to 0.48 bar, showing the significant practical implications of this work.

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

Performance Evaluation of AG-1 FC HEPA Filters and Medium in Nuclear Complex Facilities - 20434

High Efficiency Particulate Air (HEPA) filters are credited as the final barrier against the release of radioactive aerosol contamination in nearly every operating U. S. Department of Energy (DOE) and National Nuclear Security Agency (NNSA) nuclear facility. The Institute for Clean Energy Technology (ICET) at Mississippi State University maintains a research program that studies various aspects of these components of containment systems and seeks to answer key questions from across the industry. This programmatic overview will include results from recent and ongoing studies, including a study on filter design and performance envelope, degradation due to aging or fatigue, fire event impact on HEPA filters, and accelerated aging of medium and component parts. ICET has completed a study helping to define the threshold for combined elevated temperature and relative humidity resistance of separator style and separator-less style HEPA filters that are used in DOE complexes. The study provides experimental data that considers elevated temperature, elevated relative humidity (RH), and target differential pressure (dP) in an attempt to gain a better understanding or more comprehensive insight into the operating envelope of different configurations of ASME (American Society for Mechanical Engineers) AG-1 Section FC filters. Test variables including air temperature ranges of 48.88 deg. C, 54.44 deg. C, or 60 deg. C (120, 130, or 140 deg. F); air RH ranges including 60-70%, 80%, or 90+% with initial filter dP of either 497.68 Pa, 746.52 Pa, or 995.36 Pa (2 inches water column (in. w. c.), 3 in. w. c., or 4 in. w. c.) included for full bracketing of each set of conditions. A study examining the effects of aging and fatigue on nuclear grade HEPA filters and medium elucidates the physical properties of media along with testing of new and aged ASME AG-1 HEPA filters. In addition to performing autopsies of the tested, aged filters to help better understand the service life of their performance, newly designed accelerated aging chambers allow new medium to undergo accelerated aging and exposure treatments. Evaluation of degradation in physical properties and functionality requires a high population of aged filters of different ages, designs, manufacturers, and operational histories. The limited availability of this population and the impropriety of looking solely at formulations and components from past years requires experimental design that will provide foresight into future filter performance. This undertaking therefore involves not only properties analysis of aged medium from different operational histories, ages, designs, and manufacturers, but also prescience in analytical goals via accelerated aging studies of newly manufactured media. A more complete understanding of the mechanisms of degradation from past formulations as well as current and future formulations is possible. Accelerated Aging of newly manufactured media utilizes exposure treatments based on the Arrhenius equation to artificially age medium and therefore glimpse into the future. Evaluation of current media properties This testing allows for the comparison of performance and durability of new filters under upset or design basis conditions with aged filters that were in service under ambient conditions and other aged filters retained in storage. Susceptibility of HEPA filters to the effects of combustion byproducts, heat, and water is also being studied in an attempt to gain understanding on the prevention of filter failures during fire events in nuclear facilities. ICET seeks to determine the effect of filter performance due to smoke loading using characterized smoke from variable fuel compositions, geometries, and loads; the smoke capacity of various filters with different burning conditions, heat release rates and transport, as well as mass transport from the air stream will be studied in low flow containment systems. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Modeled sensitivity of multi-MA accelerator performance to electrode contaminant inventory

Significant particle-in-cell code development has enabled simulations of power flow in multi-MA accelerators to include the desorption of surface contaminants, their ionization into surface plasmas, and the impact of these plasmas on efficiency. The simulations base desorption on an Arrhenius equation, whose most significant unknown is the surface contaminant inventory. The sensitivity of power-flow simulations to this inventory is studied here using Sandia National Laboratories' Z accelerator with a 7-nH MagLIF load [Phys. Plasmas 17, 056303 (2010)]. Simulations are conducted in 3D cylindrical coordinates for the current-adder, or “convolute,” region of Z and in 2D for the final feed only. Simulated contaminant inventories are varied from 1 to 32 monolayers (MLs) in 2D, and 2 to 4 ML in 3D. The results reveal sensitivities to the local ratio of E/B⁠. The high B-field, low E-field region near the short-circuit load is insensitive to the contaminant inventory, where assumed values of 4–32 ML change the load current by ≤ 2%, and agree with experiment to within 2% at peak current. A 1-ML value is the outlier, increasing the load current by 5%, but still within measurement uncertainty. In contrast, the relatively higher E-field, lower B-field convolute region has slower contaminant desorption and higher-magnitude E-field penetration of the surface plasmas. The current loss in the convolute region does increase with contaminant inventory. The loss assuming 4 ML is 12% larger than for 2 ML, with 4 ML being the better match to experiment.

Arrhenius equation↗
DOE OSTI2021

Transformation of an Embedded Five-Membered Ring in Polycyclic Aromatic Hydrocarbons via the Hydrogen-Abstraction–Acetylene-Addition Mechanism: A Theoretical Study

Five-membered rings are constituents of many polycyclic aromatic hydrocarbons (PAHs) and their presence on the edges of large PAHs has been repeatedly observed experimentally. However, modern kinetic combustion models often do not consider the growth of PAHs through the transformation of the five-membered rings. Here, in connection with the above, we carried out a theoretical study of the mechanism of Hydrogen-Abstraction-aCetylene-Addition (HACA) transformation of an embedded five-membered ring on the armchair PAH edge to a six-membered ring, considering cyclopenta[d,e,f]phenanthrene (4,5-methylenephenanthrene) as a prototype system for this process. The potential energy surface for the reactions of cyclopenta[d,e,f]phenanthrenyl radicals produced by direct H abstractions from cyclopenta[d,e,f]phenanthrene with acetylene has been compiled at the G3(MP2,CC)//B3LYP/6-311G(d,p) level of theory including zero-point vibrational energy corrections. The computed energies and molecular parameters were then used to solve the Rice-Ramsperger-Kassel-Marcus Master Equation (RRKM-ME) in order to calculate reaction rate at various pressures and temperatures, which were fitted to modified Arrhenius equation for further kinetic modeling. The results show that the HACA transformation of the embedded five-membered ring to a six-membered ring is possible, albeit slow. The most viable reaction mechanism involves the R2 + C 2 H 2 reaction, where the acetylene molecules adds to a σ-radical in the six-membered ring adjacent to the five-membered ring via a low entrance barrier. The predominant product of R2 + C 2 H 2 is predicted to be 3 ethynyl 4H cyclopenta[def]phenanthrene Pr5 via immediate H elimination from the initial addition complex. Next, Pr5 undergoes H-assisted isomerization to 4aH-pentaleno[4,3,2,1-cdef]phenanthrene Pr4, and the latter adds an H atom eventually forming the 1-pyrenylmethyl radical Pr3: R2 + C 2 H 2 ⇌ 3 ethynyl 4H cyclopenta[def]phenanthrene (Pr5) + H ⇌ 4aH-pentaleno[4,3,2,1-cdef]phenanthrene (Pr4) + H Pr5 + H ⇌ Pr4 + H Pr4 + H → 1-pyrenylmethyl (Pr3) This HACA sequence may be competitive with the methyl radical addition to the R1 radical formed by H abstraction from the CH 2 group in the five-membered ring of cyclopenta[d,e,f]phenanthrene, which provides a pathway to pyrene following two H atom losses. Relative contributions of the two mechanisms of the five- to six-membered ring transformation would strongly depend on the branching ratios of the R1 and R2 radicals produced by the H abstractions and the available concentration of C 2 H 2 vs. CH 3 and hence differ in different flames.

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

Formation of phosphorus monoxide through the $$\mathbf {P}(^{4}S)+\mathbf {O}_{\mathbf {2}}(^{3}\Sigma ^{-})\rightarrow \mathbf {O}(^{3}P)+\mathbf {PO}(^{2}\Pi )$$ reaction

Phosphorus is a key and vital element for a diverse set of important biological molecules, being indispensable for life as we know. A deeper comprehension of its role in astrochemistry and atmospheric chemistry may aid in finding answers to how this element became available on Earth. The PO molecule is one of the main reservoirs of phosphorus in the interstellar medium (ISM), and a better understanding of the mechanisms and rate coefficients for its formation in the ISM is important for modelling its abundances. Here in this work, we perform multireference configuration interaction calculations on the formation of PO via the P( 4 S) + O 2 (3 Σ - ) reaction, analyzing its potential energy surface and rate coefficients for the global reaction on both doublet and quartet states. We also perform DFT (M06-2X) and CCSD(T) calculations, in order to compare the results. We found that the OPO system possesses a high multiconfigurational character, making DFT and CCSD methodologies not suitable for its potential energy landscape calculation. The rate coefficients have been calculated using the master equation system solver (MESS) package, and the results compared to recent experimental data. It is shown that the quartet state contributes for temperatures higher than 700K. The computed rate coefficient can be described by a modified Arrhenius equation [α(T/300) β exp (-γ/T)] with α = 1.44 x 10 -12 cm 3 s -1 , β= -1.66 and γ= 704 K.

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

Thermochemical Ablation Analysis of the Orion Heatshield

The Orion Multi-Purpose Crew Vehicle will one day carry astronauts to the Moon and beyond, and Orion's heatshield is a critical component in ensuring their safe return to Earth. The Orion heatshield is the structural component responsible for absorbing the intense heating environment caused by re-entry to Earth's atmosphere. The heatshield is primarily composed of Avcoat, an ablative material that is consumed during the re-entry process. Ablation is primarily characterized by two processes: pyrolysis and recession. The decomposition of in-depth virgin material is known as pyrolysis. Recession occurs when the exposed surface of the heatshield reacts with the surrounding flow. The Orion heatshield design was changed from an individually filled Avcoat honeycomb to a molded block Avcoat design. The molded block Avcoat heatshield relies on an adhesive bond to keep it attached to the capsule. In some locations on the heatshield, the integrity of the adhesive bond cannot be verified. For these locations, a mechanical retention device was proposed. Avcoat ablation was modelled in CHAR and the in-depth virgin material temperatures were used in a Thermal Desktop model of the mechanical retention device. The retention device was analyzed and shown to cause a large increase in the maximum bondline temperature. In order to study the impact of individual ablation modelling parameters on the heatshield sizing process, a Monte Carlo simulation of the sizing process was proposed. The simulation will give the sensitivity of the ablation model to each of its input parameters. As part of the Monte Carlo simulation, statistical uncertainties on material properties were required for Avcoat. Several properties were difficult to acquire uncertainties for: the pyrolysis gas enthalpy, non-dimensional mass loss rate (B´c), and Arrhenius equation parameters. Variability in the elemental composition of Avcoat was used as the basis for determining the statistical uncertainty in pyrolysis gas enthalpy and B´c. A MATLAB program was developed to allow for faster, more accurate and automated computation of Arrhenius reaction parameters. These parameters are required for a material model to be used in the CHAR ablation analysis program. This MATLAB program, along with thermogravimetric analysis (TGA) data, was used to generate uncertainties on the Arrhenius parameters for Avcoat. In addition, the TGA fitting program was developed to provide Arrhenius parameters for the ablation model of the gap filler material, RTV silicone.

Sixel, William↗
NASA NTRS2000

Accelerated Aging of Lead-Free Propellant

Following higher than expected 2-NDPA depletion rates in a lead-free doublebase formulation (RPD-422), an accelerated aging study was conducted to verify the depletion rates. A test plan was prepared to compare the aging characteristics of lead-free propellant and NOSIH-AA2. The study was also designed to determine which lead-free ballistic modifiers accelerated 2-NDPA depletion. The increased depletion rate occurred in propellants containing monobasic copper salicylate. Four lead-free propellants were then formulated to improved aging characteristics over previous lead-free propellant formulations. The new formulations reduced or replaced the monobasic copper salicylate. The new formulations had improved aging characteristics. Their burn rates, however, were unacceptable for use in a 2.75 inch rocket. To compare aging characteristics, stabilizer depletion rates of RPD-422, AA2, M28, and RLC 470/6A were measured or taken from the literature. The data were fit to a kinetic model. The model contained first and zero order terms which allowed the stabilizer concentration to go to zero. In the model, only the concentration of the primary stabilizer was considered. Derivatives beyond the first nitrated or nitroso derivative of 2-NPDA were not considered. The rate constants were fit to the Arrhenius equation and extrapolated to lower temperatures. The time to complete stabilizer depletion was estimated using the kinetic model. The four propellants were compared and the RPD-422 depleted faster at 45 C than both A22 and M28. These types of predictions depend on the validity of the model and on confidence in the Arrhenius relationship holding at lower temperatures. At 45 C, the zero order portion of the model dominates the depletion rate.

Furrow, Keith W.↗
DOE OSTI2024

Evidence of gas phase nucleation of nanodiamond in microwave plasma assisted chemical vapor deposition

The mechanism of ballas-like nanodiamond formation still remains elusive, and this work attempts to analyze its formation in the framework of activation energy (Ea) of nanodiamond films grown from a H2/CH4 plasma in a 2.45 GHz chemical vapor deposition system. The Ea was calculated from the Arrhenius equation corresponding to the thickness growth rate using substrate temperature (∼1000−1300 K) in all the calculations. While the calculated values matched with the Ea for nanodiamond formation throughout the literature, these values of ∼10 kcal/mol were lower compared to ∼15–25 kcal/mol for standard single crystal diamond (SCD) formation, concluding thus far that the energetics and processes involved were different. Further, the substrate preparation and sample collection method were modified while keeping the growth parameters constant. Unseeded Si substrate was physically separated from the plasma discharge by a molybdenum disk with a pinhole drilled in it. Small quantity of a sample substance was collected on the substrate. The sample was characterized by electron microscopy and Raman spectroscopy, confirming it to be nanodiamond, thus suggesting that nanodiamond self-nucleated in the plasma and flowed to the substrate that acted as a mere collection plate. It is hypothesized then, if nanodiamond nucleates in gas phase, gas temperature has to be used in the Arrhenius analysis. The Ea values for all the nanodiamond films were re-calculated using the simulated gas temperature (∼1500−2000 K) obtained from a simple H2/CH4 plasma model, giving new values within the range characteristic to SCD formation. Based on these findings, a unified growth mechanism for nanodiamond and SCD is proposed, concluding that the rate-limiting reactions for nanodiamond and SCD formation are the same.

Materials Science↗
DOE OSTI2020

Prediction of Climate-Specific Degradation Rate for Photovoltaic Encapsulant Discoloration

Encapsulant discoloration is a well-known field degradation mode of crystalline-silicon photovoltaic modules, particularly in the hot climate zones. The discoloration rate is influenced by several weathering factors, such as UV light, module temperature, and humidity, as well as the permeability of oxygen into the module. Here, a rate dependence model employing the modified Arrhenius equations to predict the degradation rate for encapsulant discoloration in different climates is presented. Two modeling approaches are introduced, which utilize the field and accelerated UV testing degradation data in conjunction with the field meteorological data to determine the acceleration factor for encapsulant browning. A novel method of accelerated UV stress testing at three simultaneous module temperatures in a single environmental chamber test run is implemented to estimate the activation energy for browning. The test was performed on three field-retrieved modules to capture the wear-out failure mechanism. The degradation in short-circuit current Isc rather than maximum power is used as a decisive parameter for the discoloration analysis. Furthermore, the developed model has been used to predict the Isc degradation rate for the Arizona field characterized by a hot and dry climate and is validated against the field-measured value. It has also been applied to other climate types, e.g., the cold and dry climate of New York.

14 SOLAR ENERGY↗
DOE OSTI2023

DISCOVR strain pipeline screening – Part I: Maximum specific growth rate as a function of temperature and salinity for 38 candidate microalgae for biofuels production

Here, to identify high productivity strains for microalgal biofuels generation, the maximum specific growth rate of 38 strains was measured as a function of salinity (i.e., 5, 15, and 35 PSU) and temperature (i.e., at 8 temperatures along a linear gradient from ca. 5 to 45°C) to determine the most suitable growth medium salinity and best growing season, respectively, for outdoor raceway pond cultivation. The following strains were evaluated: Agmenellum quadruplicatum UTEX 2268, Anabaena sp. ATCC 33081, Arthrospira fusiformis UTEX 2721, Arthrospira platensis UTEX 3086, Chlorella vulgaris NREL 4-C12, Chlorella autotrophica CCMP 243, Chlorella sorokiniana DOE1044, Chlorella sorokiniana DOE 1116, Chlorella sorokiniana DOE 1412 (UTEXB3016), Chlorella vulgaris LRB AZ-1201, Chlorococcum littorale UTEX 117, Chlorococcum sp. UTEX-B P7, Chloromonas reticulata CCALA 870, Coelastrella sp. DOE 0202, Cyanobacterium sp. AB1, Micractinium reisseri NREL 14-F2, Microchloropsis gaditana CCMP1894, Microchloropsis salina CCMP 1776, Monoraphidium sp. MONOR1, Monoraphidium minutum 26B-AM, Nannochloropsis oceanica CCAP 849/10, Oscillatoria cf. priestleyi CCMEE 5020.1-1, Picochlorum celeri TG2-WT-CSM/EMRE, Picochlorum oklahomensis CCMP 2329, Picochlorum renovo NREL 39-A8, Picochlorum soloecismus DOE 101, Porphyridium cruentum CCMP 675, Scenedesmus acutus LRB-AP-0401, Scenedesmus obliquus DOE 0152.z, Scenedesmus obliquus UTEX393,Scenedesmus rubescens NREL 46B-D3, Scenedesmus sp. IITRIND2, Stichococcus minor CCMP 819, Stichococcus minutus CCALA 727, Synechococcus elongatus UTEX2973.1, Tetraselmis striata LANL 1001, Tisochrysis lutea CCMP 1324, and Tribonema minus UTEXB3156. For each strain, the identity and the presence of bacterial cohorts was determined using 18S and 16S rDNA sequencing, respectively. The maximum specific growth rate versus temperature data were also used to determine the activation energies (Arrhenius equation) for most strains. For all strains, the measured salinity and temperature tolerance data were compared to those reported in the literature. The fastest growing strains were down-selected for subsequent biomass productivity measurements in climate-simulation photobioreactors, as reported in the next paper in the issue.

18S and 16S rDNA sequencing↗
DOE OSTI2023

Determining the oxidation behavior of matrix graphite

This work presents the oxidation behavior of matrix graphite in air. Matrix graphite, graphite powder/flakes bonded by a small amount of non-graphitic carbon, surrounds coated fuel particles in order to form cylindrical fuel compacts (in prismatic core designs) or spheres (in pebble-bed reactor designs). This work focuses on oxidation tests conducted on two matrix graphite materials, one provided by Kairos Power and the other A3 matrix graphite. Some of the tests followed American Society for Testing and Materials (ASTM) oxidation testing standards using a vertical furnace system and others were performed in a thermogravimetric analyzer (TGA). It was determined that, at temperatures of 450 °C–700 °C, the oxidation rate of the Kairos matrix graphite follows the Arrhenius equation. In comparison with A3 matrix graphite, the Kairos matrix graphite shows better oxidation resistance at high temperatures (≥550 °C), but also a higher oxidation rate at low temperatures. Both the A3 matrix graphite and the Kairos matrix graphite materials may experience preferential oxidation of the partially graphitized binder. An oxygen penetration gradient was also observed when using the three characterization methods (i.e., optical microscope, x-ray tomography [XCT], and density profile by the lathe) enlisted in this research. In conclusion, the oxygen penetration depth increases with decreasing isothermal oxidation temperature, while the center of the oxidized samples (10% weight loss) remains almost untouched even at 500 °C.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Understanding thermal and mechanical effects on lithium plating in lithium-ion batteries

Lithium plating is one of the main concerns for cell durability, as cyclable lithium and cell capacity are decreased during the process. Thermal and mechanical loads can influence cell behavior, including lithium plating driving forces, as suggested by the Arrhenius equation and the Bruggeman relation, respectively. Hence, a 1D and a 3D model are employed to study the thermal and mechanical effects on lithium plating in a lithium-ion battery with a high areal capacity (4 mAh/cm 2 ). The 1D model assumes uniform mechanical and temperature profile, while the 3D model is capable of capturing the nonuniformity of stress and temperature in the cell. Further, the effects of mechanical deformation and different temperatures on the process of lithium plating are investigated in this study. With the help of the 1D charging simulations, it is observed that lower temperatures increase the risk of lithium plating and the presence of mechanical deformations can accelerate its initiation. Furthermore, 3D charging simulations indicated that folds and boundaries of jellyrolls are more susceptible to lithium plating due to the higher stress concentration and lower temperature, respectively.

25 ENERGY STORAGE↗
DOE OSTI2021

Molecular-level insights into structure and dynamics in ionic liquids and polymer gel electrolytes

We report designing new electrolytes requires a better understanding of the correlation between their transport properties and their molecular structure. In this work, we present a detailed study of ionic liquids and polymer gel electrolytes probing their structure and dynamics by nuclear magnetic resonance (NMR) spectroscopy. In particular, ammonium- and phosphonium-based ionic liquids (ILs) are combined with different LiTFSI concentrations, and then with different ratios of poly(methylmethacrylate) polymer. The temperature dependence of self-diffusion coefficients of mobile species, D Li+ , D TFSI- and D P4441+ or D N4441+ , measured by pulsed field gradient (PFG) NMR spectroscopy obeys the Arrhenius equation. Diffusivity of [P4441][TFSI] ILs is found to be greater than those of the [N4441][TFSI] ILs in both liquid and gel electrolytes. Solid state NMR experiments including 13 C and 19 F MAS, 13 C{ 19 F}, 13 C{ 1 H} CPMAS probe the local structure and molecular-level interactions between ions and polymer in gel electrolytes, particularly for samples with higher PMMA content (≥25 wt%). Finally, the fast field cycling relaxometry has been used to unveil the rotational and translational dynamics of P4441 + or N4441 + and TFSI - by measuring 1 H and 19 F R 1 relaxation rate profiles at different temperatures. A comprehensive NMR analysis including relaxation studies at low magnetic field provides decisive new insights regarding the formation of ionic clusters and the interaction of ions with the polymer chain in the case of the gel electrolytes.

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

The P( 4 S) + NH( 3 Σ – ) and N( 4 S) + PH( 3 Σ – )reactions as sources of interstellar phosphorus nitride

Phosphorus nitride (PN) is believed to be one of the major reservoirs of phosphorus in the interstellar medium (ISM). For this reason, understanding which reactions produce PN in space and predicting their rate coefficients is important for modelling the relative abundances of P-bearing species and clarifying the role of phosphorus in astrochemistry. Here, in this work, we explore the potential energy surfaces of the $\textrm{P}(^4\textrm{S}) + \textrm{NH}(^3\Sigma^-)$ and $\textrm{N}(^4\textrm{S}) + \textrm{PH}(^3\Sigma^-)$ reactions and the formation of $\textrm{H}(^2\textrm{S}) + \textrm{PN}(^1\Sigma^+)$ through high accuracy ab initio calculations and the variable reaction coordinate transition state theory (VRC-TST). We found that both reactions proceed without an activation barrier and with similar rate coefficients that can be described by a modified Arrhenius equation ( $k(T)=\alpha\!\left( T/300 \right)^{\beta} \exp\!{(\!-\!\gamma/T)})$ with $\alpha=0.93\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=0.24\, \rm K$ for the $\textrm{P} + \textrm{NH} \longrightarrow \textrm{H} + \textrm{PN}$ reaction and $\alpha=0.88\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=1.01\, \rm K$ for the $\textrm{N} + \textrm{PH} \longrightarrow \textrm{H} + \textrm{PN}$ one. Both reactions are expected to be relevant for modelling PN abundances even in the cold environments of the ISM. Given the abundance of hydrogen in space, we have also predicted rate coefficients for the destruction of PN via H + PN collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Determining the Reaction Kinetics and Thermodynamics of a Diels–Alder Network Using Dynamic Gel Criteria

We undertook a detailed rheological investigation to evaluate the kinetic parameters of the forward and reverse Diels–Alder (DA) reactions of a model network cross-linked using a furan prepolymer and a common aromatic bismaleimide. At high temperature where the Winter–Chambon’s criterion of frequency-independence was more applicable, a multiwave technique permitted van’t Hoff analysis and calculation of the reaction thermodynamic parameters, specifically the enthalpy and entropy of the reaction: ΔH° = –38.3 ± 5.2 kJ mol –1 and ΔS° = –94.3 ± 13.4 J mol –1 . At mild temperatures where the G'–G" crossover point is experimentally convenient to measure gelation, isothermal tests were used to obtain reasonable fDA kinetic parameters from Eyring analysis such as the apparent activation enthalpy and entropy of ΔH$^{‡}_{fDA}$ = 76.8 ± 6.9 kJ mol –1 and ΔS$^{‡}_{fDA}$= –82.8 ± 22.2 J mol –1 K –1 . Comparable rheokinetic methods include cross-linking density measurements and stress relaxation tests to calculate effective kinetics, whereas the critical gel conversion was consistently applied here. As a result, rate data are fitted with the Arrhenius equation for comparison purposes and the Eyring equation to demonstrate its broader utility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗