Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactor Kinetics”

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

Thermodynamic assessment of nonstoichiometric oxides for solar thermochemical fuel production

Two-step solar thermochemical cycling (STC) based on nonstoichiometric oxides is an ideal means of solar fuel (e.g., H 2 , CO) production. Screening of nonstoichiometric oxides with excellent thermodynamic performance is key to achieving high solar-to-fuel efficiency. However, application-driven materials assessment intended for reactor-level solar fuel production performance requires mimicking realistic operating conditions of on-sun tests in a laboratory setting, which makes oxides assessment and screening an onerous task to accomplish experimentally. In this work, a rapid assessment and screening model of nonstoichiometric oxides for two-step solar thermochemical cycling assuming fixed-bed flow pattern and quasi-equilibration of the solid with the flowing gas phase is developed, with solar-to-fuel efficiency being the target function of optimization. The model accounts for the thermodynamic parameters of oxide materials and typical operating conditions of experimental thermochemical cycling. This study employed the model to explore and compare two groups of typical nonstoichiometric oxides (CeO 2 - and (LaSr)MnO 3 -based) for their maximum efficiency under their uniquely optimized cycling conditions. The results show that CeO 2 can reach a maximum efficiency of 12.9% at reduction temperature of 1500 °C, which is superior to other candidate materials, including 20 mol% Zr-doped CeO 2 (10.1%), La 0.6 Sr 0.4 MnO 3 (2.5%) and La 0.8 Sr 0.2 MnO 3 (3.3%). Even when the reduction temperature is lowered to 1350 °C, ceria yields the highest efficiency amongst the candidate STC materials considered. The optimal cycling strategy depends on the inherent thermodynamic properties of the oxides. This approach serves as a framework for assessing the maximum efficiency and optimal conditions of candidate thermochemical materials within a range of constraints rather than comparing materials under arbitrary cycling conditions, which may inherently favor one material over another. For oxidation temperatures below 800–1000 °C, the model could be further improved by considering reaction kinetics.

14 SOLAR ENERGY↗
DOE OSTI2023

Delayed Critical and Subcritical Experiments with Polyethylene Moderated Unreflected Thin 15 in. Diameter HEU Metal Plates

The thin ~15 in. diameter highly enriched uranium (HEU) metal plates were assembled to delayed criticality at the Oak Ridge Critical Experiments Facility (ORCEF) in 1969 with various thicknesses of polyethylene (varying from 1/16 to 2$\frac{3}{8}$ inches) between uranium metal plates. The average 235 U enrichment was 93.27 wt. %. These unreflected critical configurations contained 4$\frac{2}{3}$ to 20$\frac{5}{6}$ thin 15 in. diameter HEU metal plates (on loan from Los Alamos National Laboratory [LANL] and shipped to Oak Ridge National Laboratory [ORNL] on June 3, 1969). Depending on the thickness of polyethylene, the enriched uranium masses varying from 28,053 to 135,148 grams. Fractional plate sections consisted of the appropriate number of 60° pie sections. In addition to the measurement at delayed criticality, subcritical measurements were also performed by the inverse kinetic rod drop method. Prompt neutron decay constant measurements were also performed by the Rossi alpha and randomly pulsed neutron method using a time-tagged spontaneous fission californium neutron source; these are briefly reported here. At the time of these measurements in 1969, the thin HEU metal plates were in near-pristine condition with extremely little oxidation, allowing better descriptions of the uranium plates than the use of these plates in a heavily oxidized and deteriorated condition in recent reflected benchmark experiments at the LANL facility at the Nevada Test Site with these same thin highly enriched uranium metal plates. This report documents the experimental information for the measurements performed so that later researchers can perform the required uncertainty and calculational analyses and documentation to use these data for an International Nuclear Criticality Safety Benchmark Evaluation Program (ICSBEP) or a Nuclear Energy Agency (NEA) benchmark. Data from the experiments described should be acceptable for use as criticality safety benchmark experiments for the ICSBEP and the NEA nuclear criticality safety benchmark program once the uncertainty analysis on the measured neutron multiplication factors is completed. Additional data—such as the dimensional inspection reports, uranium isotopic information, and other relevant particulars—should be retrieved from the Y-12 Plant or LANL and incorporated in the final ICSBEP benchmark. Based on previous ICSBEP benchmarks with this enriched uranium metal at ORCEF, the uncertainties in $k_{eff}$ could be as low as ± 0.0002 for some configurations. Other experiments with smaller-diameter than 15 in. diameter HEU metal plates have been benchmarked in HEU-METFAST-001. The prompt neutron time decay measurements could be the basis for an International Reactor Physics Benchmark Program. Preparation of the present report is part of an effort at ORNL to document more than 15 undocumented critical and subcritical experiments enumerated in ORNL/TM-2019/18 and performed by ORNL at ORCEF and other US Department of Energy critical experiments facilities using more than 500 operational days of critical facility time. This work for this report publication was supported by the Nuclear Criticality, Radiation Transport and Safety NCSP Program at ORNL.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Modeling Radiolysis and Chemical Reactions during Dry Storage of Aluminum-clad Spent Nuclear Fuel

After aluminum-clad spent nuclear fuel (ASNF) is removed from the reactor, it is initially stored in spent fuel pools, which are specially designed water-filled basins that provide temporary cooling to reduce the temperature of the fuel assemblies and provide radiation shielding. ASNF continues to generate heat due to the radioactive decay of elements within the fuel, which persists for many years post-shutdown as the residual radioactive products decay into more stable elements. During the wet storage period, an oxyhydroxide layer composed of boehmite/bayerite forms on the surfaces of the aluminum cladding from exposure to water in the pools. Road-ready packaging for long-term disposition of the ASNF involves dry storage in helium backfilled DOE standard canisters (DSCs). When the ASNF is removed from water storage and dried, most of the water is removed, but some physisorbed and chemisorbed water remains in the oxyhydroxide layers. This residual water can produce hydrogen when exposed to radiation from the ASNF during dry storage. Predicting hydrogen accumulation over time in the DSCs is critical for long-term storage considerations. Previous modeling efforts have developed coupled computational fluid dynamics (CFD)-chemical models to simulate temperature, pressure, and gas phase concentrations within the DSCs. These models use the thermal field predicted by CFD as input to a radiolysis model for the gas phase and the surface oxyhydroxide layer chemistry. Given the long storage period of the DSCs and the impracticality of long-term experiments, a simulation-based approach is necessary to assess chemical evolution within the canisters. This study advances the development of a modeling framework designed to simulate the chemical evolution of spent fuel canisters. Both thermal and radiation-driven reactions are considered, with radiation kinetics quantified using G-values. Sensitivity analysis identifies key parameters influencing species composition. Reaction pathway diagrams offer insight into dominant species formation routes, enabling more effective comparisons between model predictions and experimental observations, particularly regarding the production of hydrogen. Results show that the model predicts significant hydrogen gas production with minimal oxygen generation, primarily due to hydrogen formation via boehmite pathways. These findings underscore the importance of accurately characterizing surface-bound species and radiolysis kinetics. A deeper understanding of these mechanisms is critical for evaluating the long-term safety of nuclear waste storage.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2021

Nitrogen-doped char as a catalyst for wet oxidation of phenol-contaminated water

Catalytic wet oxidation (CWO) of aqueous effluents rich in organic compounds is a very promising technology for the treatment of liquid wastes from biomass conversion processes. CWO reactions occur through the formation of free radical species, produced in the presence of an oxidant, which act on organic contaminates in the effluent. Although the reaction is well known, there exists a lack of affordable catalysts to conduct this process at the lower temperatures and pressures in novel bioenergy processes. Here, this study assessed the catalytic effect of nitrogen-doped chars as such an option. Phenol in aqueous solution was used as a model waste effluent. Treatment was conducted at moderate temperatures (190 to 260 °C), oxygen partial pressure of 1 MPa, and reaction times of 15, 30, and 45 min in stainless steel and glass-lined tube reactors. High pressure liquid chromatography (HPLC) analyses of the products quantified phenol and by-product concentrations used in the calculation of reaction activation energy. The char catalyst was studied by X-ray photoelectron spectroscopy (XPS), transmission electron microscopy (TEM), and scanning electron microscopy (SEM) in order to gain insight into its structure and surface composition. The results indicate that nitrogen-doped char catalysts accelerate the oxidation of phenol by decreasing its reaction activation energy from 82.2 kJ/mol (non-catalyzed) to 40.4 kJ/mol (catalyzed). An analysis from first principles using density functional theory (DFT) was conducted to ascertain which N functional group has the most significant impact on free radical formation in the presence of oxygen. Among all the N functional groups studied, the dipyridinic functional groups showed the most promising characteristics to facilitate the formation of hydroxyl free radicals.

Advanced Oxidation Processes for Organic Pollutant↗
DOE OSTI2020

Report on the FY20 Characterization of Age Hardened Alloy 709 Materials

This report provides the results on microstructural characterization of the first commercial heat of Alloy 709 plate material undergone additional age hardening treatment. In FY19, through thermodynamic-kinetic modeling and microstructural characterization, it is determined that the 775°C/10 h age-hardening treatment can promote desired precipitation on dislocations and has the potential to further improve the mechanical performance of the material. In FY20, the treatment protocol was applied to the ESR-1150°C condition, and the microstructure was then compared with the ESR-1100°C condition. It was found that the 1150°C solution annealing (SA) dramatically reduced the dislocation density compared to the 1100°C SA, resulting a smaller number density of Nb-rich fine precipitates which grow on dislocations during the age hardening treatment. Other features are very similar between the two. The effect of age hardening on creep performance is discussed. In addition, TEM characterization of a creep-ruptured sample (600°C/330 MPa) was performed. Specimens from the grip region and the gauge region was studied. The grip region specimen was effectively aged at 600°C for extended period, while the gauge region specimen was mechanically deformed at 600°C. The microstructures are very different, implying the importance of a dynamical effect on precipitation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗
DOE OSTI2024

Tracking the reaction networks of acetaldehyde oxide and glyoxal oxide Criegee intermediates in the ozone-assisted oxidation reaction of crotonaldehyde

The reaction of unsaturated compounds with ozone (O 3 ) is recognized to lead to the formation of Criegee intermediates (CIs), which play a key role in controlling the atmospheric budget of hydroxyl radicals and secondary organic aerosols. The reaction network of two CIs with different functionality, i.e. acetaldehyde oxide (CH 3 CHOO) and glyoxal oxide (CHOCHOO) formed in the ozone-assisted oxidation reaction of crotanaldehyde (CA), is investigated over a temperature range between 390 K and 840 K in an atmospheric pressure jet-stirred reactor (JSR) at a residence time of 1.3 s, stoichiometry of 0.5 with a mixture of 1% crotonaldehyde, 10% O 2 , at an fixed ozone concentration of 1000 ppm and 89% Ar dilution. Molecular-beam mass spectrometry in conjunction with single photon tunable synchrotron vacuum-ultraviolet (VUV) radiation is used to identify elusive intermediates by means of experimental photoionization energy scans and ab initio threshold energy calculations for isomer identification. Addition of ozone (1000 ppm) is observed to trigger the oxidation of CA already at 390 K, which is below the temperature where the oxidation reaction of CA was observed in the absence of ozone. The observed CA + O 3 product, C 4 H 6 O 4 , is found to be linked to a ketohydroperoxide (2-hydroperoxy-3-oxobutanal) resulting from the isomerization of the primary ozonide. Products corresponding to the CIs uni- and bi-molecular reactions were observed and identified. A network of CI reactions is identified in the temperature region below 600 K, characterized by CIs bimolecular reactions with species like aldehydes, i.e., formaldehyde, acetaldehyde, and crotonaldehyde and alkenes, i.e., ethene and propene. The region below 600 K is also characterized by the formation of important amounts of typical low-temperature oxidation products, such as hydrogen peroxide (H 2 O 2 ), methyl hydroperoxide (CH 3 OOH), and ethyl hydroperoxide (C 2 H 5 OOH). Detection of additional oxygenated species such as alcohols, ketene, and aldehydes are indicative of multiple active oxidation routes. This study provides important information about the initial step involved in the CIs assisted oligomerization reactions in complex reactive environments where CIs with different functionalities are reacting simultaneously. It provides new mechanistic insights into ozone-assisted oxidation reactions of unsaturated aldehydes, which is critical for the development of improved atmospheric and combustion kinetics models.

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

Mechanistic pathways and role of oxygen in oxidative coupling of methane derived from transient kinetic studies

Oxidative coupling of methane (OCM) is a promising industrial process to upgrade natural gas to high value chemicals. In this study, Temporal Analysis of Products (TAP) and steady-state experiments were conducted to distinguish how the composition of surface and gas phase oxygen influence mechanistic details of the selective conversion of CH 4 to C 2 H 4 over the Mn-Na 2 WO 4 /SiO 2 catalyst. The results from TAP studies indicate that methane activation on this catalyst proceeds predominantly via a short-lived, transient surface oxygen species and there is a competition for this species to form either CO or methyl radicals on the surface. This active species has a total lifetime of 3 s and is identified to have a dioxygen (e.g. O 2 2- or O 2 -) form. We show that the concentration of the transient surface oxygen species significantly impacts the OCM performance. Oxygen attributed to the catalyst lattice (in a singular form e.g., O - ), is found to activate methane to a lesser degree, but exclusively forms CO 2 . Evidence for these surface pathways for methyl radical, CO and CO 2 formation identified by TAP are also validated through steady-state experiments. Finally, by distinguishing different catalyst oxygen species and their role in selective/nonselective pathways, important screening criteria have been identified for the advancement of superior catalyst formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High temperature gamma radiation-induced chromium redox chemistry via in situ spectroscopic measurements

Chromium ions can make their way into the primary coolant of nuclear power reactors from the corrosion of stainless-steel reactor components, decreasing the material's corrosion resistance and resulting in increased transport of further corrosion products. Despite these potential effects, the radiation-induced redox speciation of chromium ions in aqueous solution is not well understood, especially at the elevated temperatures experienced by reactor coolants. Here, in the present work, we report new experimental results demonstrating that in aerated aqueous solution, the radiolytic oxidation of Cr(III) to Cr(VI) occurs at pH 4, while the reduction of Cr(VI) to Cr(III) occurs at pH 2. The oxidation of Cr(III) is primarily attributed to the reaction of the hydroxyl radical (˙OH) with the Cr(OH) 2+ species, while the reduction of Cr(VI) is attributed to reactions involving the hydrated electron (e aq - ) and hydrogen atom (H˙). Additionally, the steady-state equilibrium yield of Cr(VI) from the gamma irradiation of pH 4 Cr(III) solutions decreased with increasing temperature (over a range of 37–195 °C). This observation indicates that the activation energy of the Cr(VI) reduction reactions is higher than that for the Cr(III) oxidation reactions, such that it becomes relatively more favorable at higher temperatures. Overall, these data are important for the development of complementary multiscale models for the prediction of metal ion speciation in high temperature radiation environments.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Biogas sequestration to carbon nanofibers via tandem catalytic strategies

Upgrading decentralized biogas represents a sustainable route to produce valuable products while mitigating two potent greenhouse gases, namely, methane (CH 4 ) and carbon dioxide (CO 2 ). Conventional dry reforming of CH 4 with CO 2 yields syngas with low H 2 /CO ratios (≤1) and requires high temperatures (>800 °C) to overcome equilibrium constraints and abate coke deposition, which limits commercial implementation. Here we demonstrate the conversion of biogas into value-added carbon nanofibers via reaction integration in tandem reactors, while reducing the reaction temperature, shifting equilibrium limits and yielding H 2 -enriched syngas (H 2 /CO = 2–3) as a byproduct. Experimental and theoretical insights reveal that potassium (K) modification enhances carbon nanofiber formation due to synergistic effects via a balanced interplay between KO x -induced cobalt facets and cobalt carbide species. In conclusion, the energy cost and CO 2 footprint analyses highlight the potential advantages of tandem processes for the sustainable upgrading of biogas into valuable solid carbon products.

09 BIOMASS FUELS↗
DOE OSTI2024

Synthesis, Processing, and Use of Isotopically Enriched Epitaxial Oxide Thin Films

Isotopic engineering has emerged as a key approach to study the nucleation, diffusion, phase transitions, and reactions of materials at an atomic level. It aims to uncover mass transport pathways, kinetics, and operational and failure mechanisms of functional materials and devices. Understanding these phenomena leads to deeper insights into important physical processes, such as the transport of ions in energy conversion and storage devices and the role of active sites and supports during heterogeneous catalytic reactions. Likewise, isotopic engineering is being pursued as a means of modifying functionality to enable future technological applications. In this Account, we summarize our recent work employing isotope labeling (e.g., 18 O 2 and 57 Fe) during thin film synthesis and postgrowth processing to reveal growth mechanisms, defect chemistry, and elemental diffusion under working and extreme conditions. Isotope-resolved analysis techniques with nanometer-scale spatial resolution, such as time-of-flight secondary ion mass spectrometry and atom probe tomography, facilitate the accurate quantification of isotopic placement and concentration in our well-defined heterostructures with precisely positioned, isotope-enriched layers. By measuring the nanometer-scale redistribution between natural abundance and isotopically enriched oxygen layers during the deposition of Fe 2 O 3 and Cr 2 O 3 by molecular beam epitaxy, we identified intermixing processes driven by surface adatoms occurring both at the film growth surface and within the first few layers below the surface. Further insights into synthesis mechanisms were gained by studying the tungsten oxide thin films grown by evaporating WO 3 powder in the presence of background 18 O 2 , revealing minimal incorporation of background oxygen during the film formation process. Thermal and radiation-enhanced diffusion in epitaxial Fe and Cr oxides were precisely tracked using 18 O and 57 Fe tracer layers incorporated into model epitaxial oxide thin films. This approach has allowed us to access thermal diffusion behavior at lower temperatures than previously measured, revealing a potential changeover in diffusion mechanism. Understanding radiation-enhanced diffusion in model oxides that represent the surface layers on the structural components of nuclear reactors informs our understanding of their corrosion behavior under irradiation. Isotopic labeling can also provide unique insights into the surface exchange reactions and defect chemistry of electrocatalysts. For instance, tracking the change in 18 O concentration at the surface of an epitaxial LaNiO 3 thin film after the electrocatalytic oxygen evolution reaction revealed the participation of lattice oxygen, confirming a hypothesis that had been proposed previously. Lastly, we highlight a new direction wherein we perform in situ processing studies utilizing isotopic tracers in conjunction with model epitaxial thin films within the atom probe tomography instrument. Additionally, this Account illustrates the great potential of isotopic engineering to enable fundamental mechanistic insights into physical processes and engineer functional properties in epitaxial films, heterostructures, and superlattices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Initial Assessment of TeF 6 Adsorption

Research and development that supports the management of off-gases from nuclear fuel reprocessing has historically been focused on the off-gas streams that arise from aqueous reprocessing technology. However, as Gen-IV reactor development pathways move toward deployment, alternative spent nuclear fuel (SNF) processing and disposition pathways have been considered more actively. This work is focused upon aspects of fluoride volatility (FV) processing. The versatility of this method for deployment against multiple types of spent fuel encourages the continued advancement of both the primary separations processes and the secondary processes, including waste treatment, material control and accountability, and engineering designs. Recent work noted that TeF 6 , the most highly volatile fluorinated compound produced during FV processing, did not have a clear abatement technology recommended in the literature. Thus, this work performed scoping tests on activated alumina and copper shot to assess whether they could be used to remove TeF 6 from gas streams that bear F 2 . Previous work on this topic was not well described in the literature and was not performed with excess F 2 in the stream as would be typical of spent fuel processing via FV. To support an understanding of the concentration of TeF 6 contacting the adsorbent beds, a series of preliminary testing identified the TeF 6 production rate and the equilibrium concentration of TeF 6 in the gas stream contacting the adsorbent. Excess F 2 was determined to not affect the ability of activated alumina to remove TeF 6 from the gas stream quickly and completely. A determination of whether excess F 2 affected the distribution depth of Te in the sorbent bed is still pending analysis of the used sorbent. Literature suggests that when activated alumina is near saturation, TeF 6 could migrate from the sorbent bed. Future testing should investigate this possibility. Copper metal did not adsorb TeF 6 in the presence of F 2 across the sorbent temperature range of 50 to 335°C. F 2 was fully removed by the copper bed. Although preliminary thermodynamics would indicate adsorption to be energetically favorable, other factors that impact adsorption (e.g., slow kinetics, excess fluorine on the copper surface, an unfavorable transition state) are likely preventing the adsorption of TeF 6 at an easily measurable rate. The work also investigated TeF 6 production from multiple forms of Te in the temperature range of 100 to 250°C. No previous study had assessed the initial reaction rates for this process. A carefully designed study allowed determination of the activation energy for the production of TeF 6 from Te metal. The substantial amount of data collected during this study merits analysis beyond what is described here. A more in-depth kinetic analysis will be pursued. Additional analytical results will provide the ability to benchmark adsorption coefficients for TeF 6 , understand the distribution of TeF 6 within the alumina bed, and better understand the effect of F 2 partial pressure on Te fluorination. The data from this report, as supplemented by these additional analyses, will be submitted to a peer-reviewed journal.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Integration of Omics into a New Comprehensive Rate Law for Competitive Terminal Electron-Accepting Processes in Reactive Transport Models: Application to N, Fe, S, and Contaminant Transformations in Stream and Wetland Sediments

Surface waters represent important sources of alternate energy and drinking water in the United States, and characterizing the biogeochemical processes that affect surface water quality is relevant to the DOE-BER mission. Sediment biogeochemical processes regulate the release of carbon (C), nutrients, and contaminants to surface waters and thus influence water quality. Sediment biogeochemical processes are dynamic and affected by the deposition and remobilization of solid material and changes in environmental conditions driven by water discharge variations. Wetlands are important natural filters of surface waters which may either trap, metabolize, or mobilize nutrients and contaminants. Despite their importance, biogeochemical processes regulating nutrient and contaminant release and C transformation in stream and wetland sediments cannot be predicted accurately by current mathematical models. These reactive transport models largely rely on detectable changes in geochemical conditions to activate metabolic processes, do not accurately account for the competition between microbial processes, and poorly constrain effects of hydrological perturbations on biogeochemical processes. In this BER-SBR exploratory project, metagenomic and geochemical signatures were combined to identify microbially-mediated redox processes in anaerobic stream and wetland sediments from the Savannah River Site (SRS, ANL SFA) and East Fork Poplar Creek (EFPC, ORNL SFA) that play important roles in C, uranium (U), and mercury (Hg) transformations. In addition, sediment incubations were conducted to examine the competition between anaerobic respiration processes Finally, new rate laws were developed for reactive transport models that rely on complementary metagenomic and geochemical signatures to identify the underlying anaerobic microbial processes in stream and wetland sediments, describe the competition between the dominant metabolic processes involved in nutrient release and U and Hg mobilization, and more accurately quantify carbon transformation and the response of microbial processes to changes in redox conditions associated with hydrological forcing. These rate laws were optimized in batch reactors with SRS wetland sediments, where iron and sulfate reduction dominate. Anaerobic carbon remineralization processes followed the expected thermodynamic sequence of microbial respiration with depth in both sediments, except that geochemical signals indicated that sulfate reduction was inactive in EFPC sediments and moderate in SRS wetland sediments. Estimates indicated that microbial iron reduction contributed to at least half of the production of reduced iron in these sediments. Incubations demonstrated that nitrate reduction, denitrification, and dissimilatory nitrate reduction to ammonium were active in the natural EFPC sediment and activated upon nitrate amendment in these nitrate-rich sediments. In turn, these processes were outcompeted by the addition of either iron oxides or sulfate as alternative terminal electron acceptors. Although geochemical products of sulfate reduction were not detected in the incubations, the abundance of sulfate reduction genes increased with depth in the sediment and was equally more pronounced in treatments amended with either iron oxides or sulfate. Simultaneously, anaerobic sulfide oxidizing bacteria coupling sulfide oxidation to DNRA (and not conventional denitrification) were apparently enriched over time, regardless of the treatments. These findings indicate that sulfate reduction is important in freshwater stream sediments and probably catalyzed by a cryptic sulfur cycle involving nitrogen species, in which the sulfur products from sulfate reduction are immediately removed by side reactions and not detectable by geochemical measurements alone. Similar experiments in SRS sediments, however, demonstrated little interaction between nitrogen and sulfur cycling microorganisms. Sulfate reduction was impacted by the addition of more thermodynamically favorable electron acceptors, suggesting either that iron-reducing microorganisms outcompeted sulfate-reducing microorganisms for organic substrate, depleting the stock of electron donor available for sulfate reduction, or that the cryptic sulfur cycle was shunted by the precipitation of FeS generated as a result of the abiotic reduction of iron oxides by dissolved sulfide. A diagnostic modeling exercise was conducted to further investigate the competition between terminal electron accepting processes. As conventional kinetic models typically do not account for cryptic cycles and used inaccurate formulations to describe competition between microbial communities, new metabolic rate laws were developed that explicitly express the electron acceptor-specific enzyme of each energetically favorable metabolic process based on gene abundance detected in the incubations. The model was tested with the sediment slurry incubation data to determine whether substrate competition could explain the decrease in sulfate reduction observed in the presence of iron oxide competitor. The model was able to reproduce geochemical concentrations really well in each treatment once the model was calibrated with the unamended control, suggesting that microbial competition was indeed driven by thermodynamic considerations. Overall, carbon remineralization processes and rates will be reproduced much more realistically with the new metabolic rate laws.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Extended Modeling of DOE Sealed Canisters with Updated Chemistry Models

Road-ready and final disposition packaging configurations for the advanced test reactor (ATR) fuel currently specifies storage within helium backfilled DOE sealed standard canisters. The aluminum cladding of the ATR fuel contains an oxyhydroxide layer of boehmite/bayerite that generates hydrogen when subjected to irradiation. Understanding the effect of this hydrogen buildup over time to important for long term storage considerations. Previous modeling efforts have built a coupled CFD-chemical model to simulate the temperature gas phase concentrations within the DOE sealed standard canisters. These models have coupled the temperature conditions to both the gas phase radiolysis chemistry and the surface chemistry associated with the oxyhydroxide layer. The previous iteration of the model utilized constant G-values for the hydrogen generation a 50-year period. This new iteration of the model utilizes new experimental data to update the hydrogen generation rate, as well as increase the simulated time to a 200-year period. Given the half-life assumed for the primary Cs-137 isotope responsible for gamma radiation in the ATR spent fuel, the 200-year period decreases the decay heat and dose rate of the spent fuel by a factor of 100, which combined with the updated chemistry substantially decreases new hydrogen generation. Continued experimental work has identified trends for aluminum surrogate samples with oxyhydroxide layers for tests done at higher dose rates. At low initial doses a fast generation rate of hydrogen occurs which starts to roll over to a lower generation as the total dose applied increase. A small-scale chemical model was built to replicate as mini-canister surrogate system at SRNL as well as for the smaller capsule tests performed at INL. A variety of chemical models to capture this effect were tested, and a back reaction of H radical absorbing onto the surface, or inhibition of the reaction by significant H 2 cover gas were both able to fit both the mini-canister and the small capsule test data. Both kinetic fits were used to generate data from the new 200-year simulation. As additional experiments continue, the kinetic fits may be adjusted to adapt to new data. However, updated data shows that hydrogen atmosphere has little effect on actual generation data, so the model was reverted to use a star-stepped G-value for low-dose and high-dose regions. For the undried fuel case, the three models differ significantly with an end concentration of 1% for back reaction, 4.2% for inhibition, and 13.1% for constant G-value for the nominal scenario. For the dried fuel case, the nominal cases showed end concentrations of 0.23% for back reaction, 1.9% for inhibition, and 4.2% for const G-value for the nominal scenario. For the constant G-value case that is less conservative than the other two, the pressure for undried fuel increases to 1.94 atm over 200 years for the nominal case and 2.11 atm for the high decay heat case. The total hydrogen concentration after 200 years is 4.92% for the low decay heat case and 27.4% for the high decay heat case for undried fuel and is 1.4% and 9.6% for low and high decay heat fuel for the dried fuel case. If small amount of residual air is present, the potential for nitric acid formation of 595, 1568, and 2816 ppm for the lower, nominal, and upper fuel decay heat can occur. At long-term timeframes no significant shift in major species present occurs, so no appreciable amount of oxygen is present in the system. The continued rate of hydrogen generation in the canister occurs at an increase of 0.04% by mole over the final 10 years. Increasing the model range out to 1000 years continues to drop the hydrogen generation rate through decreasing dose rate.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

FY23 Progress Report on Viscosity and Thermal Conductivity Measurements of Molten Salts

As presented in this report, thermal conductivity and viscosity measurements were performed on key chloride pseudo-binary molten salt systems of relevance to molten salt reactor developers. Thermal conductivity measurements were conducted with a variable gap technique, in which a known heat flux is driven across a molten salt specimen and the temperature difference is measured, allowing calculation of the thermal conductivity. This is achieved by establishing a small gap between the bottom of a cylindrical inner containment, which houses electrical heating elements, and an outer containment, which houses cooling channels; the gap size can be varied by compression of a formed bellows. A new calibration scheme was developed herein, involving a correction factor to the heat flux based on He measurements at various temperatures. Furthermore, the data processing methodology was improved to minimize the impact of radiative heat transfer in the calculation of salt specimen thermal conductivity from the temperature difference measurements. Viscosity measurements were conducted with a rolling ball viscometer, in which a ball rolls some known distance in an angled tubular crucible, and the terminal velocity can be used to calculate the viscosity of the salt. The measurement can be performed in a quartz crucible, with which a standard camera can be used to track the ball, or in a metal crucible, with which x-ray radiography is required to track the ball. A new custom x-ray system was made and dedicated to the rolling ball viscometer to enable high throughput automated measurements with salts which require containment with metal. Both systems have been integrated with a ventilation stack which allows for off-gassing of radioactive material, enabling future measurements with U-bearing salts. The thermal conductivity measurements performed herein were with NaCl-KCl (44 mol% NaCl). This salt system was measured in the previous fiscal year, however the thermal conductivity values obtained were comparatively low, and so the system was remeasured with the aforementioned calibration scheme and improved post-processing techniques. The newly obtained thermal conductivity values for NaCl-KCl (44 mol% NaCl) indicate good agreement with kinetic theory and ab-initio models (within 5–10 %). The viscosity measurements performed herein were with three different compositions of NaCl-KCl: 75, 50, and 25 mol% NaCl; a new calibration scheme was applied to account for variable flow effects in the laminar regime. The results show reasonable agreement with literature (5–20 %, depending on the temperature, composition, and study); however, literature values are likely higher than true values based on pure end-member measurements performed in the comparative studies. The measurements conducted herein do show a trend such that viscosity increases with increasing NaCl concentration, which agrees well with one of two comparative studies. In general, the measurements conducted herein gives confidence in the capability to use these systems to accurately measure thermal conductivity and viscosity of actinide-bearing salts within the next fiscal year.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

The role of ion-scale micro-turbulence in pedestal width of the DIII-D wide-pedestal QH mode

The low-edge rotation, intrinsically ELM-free, and improved confinement wide-pedestal quiescent H-mode (QH-mode), discovered in DIII-D tokamak, has pedestal widths exceeding the EPED-kinetic-ballooning mode (KBM) model scaling typically by at least 25%. Ion-scale ($k_yp_s$ < 1) microturbulence and its role in setting the pedestal structure is investigated using the radially local δ$f$ gyrokinetic code CGYRO. The electromagnetic trapped electron mode (TEM) is unstable at the pedestal top, while plasma beta (β$_e$) is ∼60% below the KBM onset threshold and the electron temperature gradient mode is found to be unstable in the peak gradient region. Nonlinear simulation reveals that the ion-scale turbulence could produce electron energy flux consistent with the flux inferred from power balance at the pedestal top, with a reasonable variation of the local shearing rate; and the local neoclassical transport from NEO is dominant over the simulated turbulent transport in the ion energy flux channel. The simulated ion-scale turbulence produces much lower electron energy flux than inferred from experiment in the pedestal peak gradient region. A correction to the EPED-KBM pedestal width scaling is obtained based on the two-dimensional scan of pedestal top plasma beta (β$_e$) and normalized electron density and temperature scale lengths,${a}$/${L_n}_e$, ${a}$/${L_T}_n$using CGYRO linear simulations. Mode transitions among TEM, micro-tearing mode, ion-temperature gradient mode and KBM, are observed in the 2D scan at the pedestal top. A fixed normalized growth rate for these drift-type modes is taken to determine the pedestal width scaling, which shows good consistency with the QH experimental database on pedestal heights and widths. The onset of KBM instabilities and the local E × B shear suppression criterion set the lower and upper limit for the pedestal width of standard QH-mode, wide-pedestal QH-mode and type-I ELMy H mode. A potentially higher and wider pedestal is expected from the new scaling of pedestal width. This work presents an improved understanding of the ion-scale micro-turbulence of wide-pedestal QH-mode and sheds light on a promising scenario for future reactors, including ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Numerical Investigation of Fuel Property Effects on Mixed-Mode Combustion in a Spark-Ignition Engine

In this research, lean mixed-mode combustion is numerically investigated using computational fluid dynamics (CFD) in a spark-ignition engine. A new E30 fuel surrogate is developed using a neural network model with matched octane numbers. A skeletal mechanism is also developed by automated mechanism reduction and by incorporating a NO x submechanism. A hybrid approach that couples the G-equation model and the well-stirred reactor model is employed for turbulent combustion modeling. The developed CFD model is shown to well predict pressure and apparent heat release rate (AHRR) traces compared with experiment. Two types of combustion cycles (deflagration-only and mixed-mode cycles) are observed. The mixed-mode cycles feature early flame propagation and subsequent end-gas auto-ignition, leading to two distinctive AHRR peaks. The validated CFD model is then employed to investigate the effects of NO x chemistry. The NO x chemistry is found to promote auto-ignition through the residual gas, while the deflagration phase remains largely unaffected. Sensitivity analysis is finally performed to understand effects of fuel properties, including heat of vaporization (HoV) and laminar flame speed (S L ). An increased HoV tends to suppress auto-ignition through charge cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased S L is found to significantly promote both flame propagation and end-gas auto-ignition. The promoting effect of S L on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end-gas.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗