SEARCH · Engineering Papers

Results for “Enthalpy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Thermocatalytic Heat Pipes for Geothermal Resource Recovery

Heat pipes are an important technology that allow orders of magnitude faster heat transfer than simple conduction. However, operating principles in heat pipes place fundamental bounds on their performance (critical heat flux and efficiency). Conventional heat pipe functionality is inherently tied to vaporization and condensation of the working fluid charged in the heat pipe. These fluids each have different operating temperature ranges based on the capillary, entrainment, sonic, and boiling limits of the heat pipe design. These limits, typically the capillary limit, dictate the maximum heat flux a heat pipe can carry, and most importantly for geothermal systems, the distance over which the pipes can operate (100 to 200 m maximum under optimum conditions). A thermocatalytic heat pipe breaks the inherent limitations of phase change thermo- and hydrodynamics and can transform heat pipe technology as a potentially more efficient means of extracting heat from a geothermal resource. The thermocatalytic heat pipe uses a working fluid to transport both sensible and chemical heat. An endothermic chemical reaction at depth removes heat from the reservoir and produces reactive intermediates, which are transported to the surface and used to run a reverse exothermic reaction that releases heat for use in power generation or other useful purposes. This technology offers two distinct advantages over conventional geothermal heat recovery technologies: (1) lower heat loss to the rock outside of the geothermal reservoir, and (2) higher heat transfer rates to the well field within the geothermal reservoir. Both advantages offer opportunity to reduce risks and lower costs of geothermal energy recovery. In this report, we discuss an initial effort to assess the efficacy and limitations of this technology for extracting heat from both porous/permeable and nominally impermeable geothermal reservoirs. Numerical simulation capabilities of the STOMP-GT code were enhanced to enable simulations of thermochemical heat pipes traversing geothermal reservoirs. An array of potential thermochemical reaction systems was evaluated and screened. Of these, an ethanol dehydration reaction was most promising in the vapor-liquid reaction set. A solid-phase dehydration reaction (CuSO4·5H2O) showed the highest reaction enthalpy per unit volume but would require development of a nonaqueous carrier fluid to implement it in a heat pipe. Subsurface reservoir simulations predicted long-term performance of the heat pipes for each geothermal reservoir type. The performance of U-shaped wells and coaxial wells was evaluated for a suite of reactions for both hydrothermal and hot dry rock reservoirs and was compared with a baseline case of simply pumping water through the wells. The heat pipe technology was additionally evaluated for an enhanced geothermal system (EGS) with an injection borehole, production borehole, and intervening hydraulically conductive fracture. All reservoir types showed significant improvement in heat recovered over a 20-year operating period ranging from a 1.8X increase for the hot dry rock case to more than 2.5X more energy recovered for the EGS case.

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Two-Phase Equation of State for Al 2 O 3

Two-phase equation of state (EOS) for AI 2 O 3 (alumina) is calibrated using the experimental/simulation calibration data uploaded to the caldata repository. The two phases are solid α-Al 2 O 3 (corundum, sapphire, ruby) and its melt. The enthalpy of formation is set to that available from the NIST-JANAF database at T = 0K, P =1 bar. The expected applicability range of the EOS is from the ambient conditions to temperatures of ~ 5000 K and pressures of ~ 50-100 GPa.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Constitutive Model Development for Aging Polymer Encapsulants (ASC P&EM FY2021 L2 Milestone 7836)

This SAND report fulfills the completion requirements for the ASC Physics and Engineering Modeling Level 2 Milestone 7836 during Fiscal Year 2021. The Sandia Simplified potential energy clock (SPEC) non-linear viscoelastic constitutive model was developed to predict a whole host of polymer glass physical behaviors in order to provide a tool to assess the effects of stress on these materials over their lifecycle. Polymer glasses are used extensively in applications such as electronics packaging, where encapsulants and adhesives can be critical to device performance. In this work, the focus is on assessing the performance of the model in predicting material evolution associated with long-term physical aging, an area that the model has not been fully vetted in. These predictions are key to utilizing models to help demonstrate electronics packaging component reliability over decades long service lives, a task that is very costly and time consuming to execute experimentally. The initiating hypothesis for the work was that a model calibration process can be defined that enables confidence in physical aging predictions under ND relevant environments and timescales without sacrificing other predictive capabilities. To test the hypothesis, an extensive suite of calibration and aging data was assembled from a combination of prior work and collaborating projects (Aging and Lifetimes as well as the DoD Joint Munitions Program) for two mission relevant epoxy encapsulants, 828DGEBA/DEA and 828DGEBA/T403. Multiple model calibration processes were developed and evaluated against the entire set of data for each material. A qualitative assessment of each calibration's ability to predict the wide range of aging responses was key to ranking the calibrations against each other. During this evaluation, predictions that were identified as non-physical, i.e., demonstrated something that was qualitatively different than known material behavior, were heavily weighted against the calibration performance. Thus, unphysical predictions for one aspect of aging response could generate a lower overall rating for a calibration process even if that process generated better quantitative predictions for another aspect of aging response. This insurance that all predictions are qualitatively correct is important to the overall aim of utilizing the model to predict residual stress evolution, which will depend on the interplay amongst the different material aging responses. The DSC-focused calibration procedure generated the best all-around aging predictions for both materials, demonstrating material models that can qualitatively predict the whole host of different physical aging responses that have been measured. This step forward in predictive capability comes from an unanticipated source, utilization of calorimetry measurements to specify model parameters. The DSC-focused calibration technique performed better than compression-focused techniques that more heavily weigh measurements more closely related to the structural responses to be predicted. Indeed, the DSC-focused calibration procedure was only possible due to recent incorporation of the enthalpy and heat capacity features into SPEC that was newly verified during this L2 milestone. Fundamentally similar aspects of the two material model calibrations as well as parametric studies to assess sensitives of the aging predictions are discussed within the report. A perspective on the next steps to the overall goal of residual stress evolution predictions under stockpile conditions closes the report.

36 MATERIALS SCIENCE↗

DOE OSTI2022

An Open-Source Framework for the Computational Analysis and Design of Autothermal Chemical Processes

The main project object was to develop software tools for simulating non-equilibrium autothermal processes, improving the prospects for identifying and designing such systems. The project demonstrates the use of these tools to simulate autothermal pyrolysis, a process recently developed at the pilot-plant scale at Iowa State University. In such process, instead of externally heating a reactor to pro-vide the enthalpy of pyrolysis, sufficient oxygen in the form of air is introduced into the reactor to support partial oxidation of reactants and products with the exothermic energy released supporting endothermic pyrolysis reactions. A fluidized bed is used to assure good mixing of biomass and oxidant and provide an isothermal reaction environment. The amount of oxygen required depends upon the kind of biomass being pyrolyzed and parasitic heat losses from the reactor. For example, for woody biomass pyrolyzed under conditions that simulate adiabatic operation, equivalence ratios can be as low as 0.06, compared to 0.20 or higher for autothermal gasifiers. By removing the heat transfer bottleneck of conventional pyrolysis, operation in autothermal mode allowed a significant increase in reactor throughput process, approaching five times the throughput of the conventionally operated pyrolyzer. Different simulation strategies were considered and developed: a zero-dimensional chemistry model was used to verify the applicability of kinetic schemes to predict biomass fast pyrolysis in autothermal conditions. Conventional chemical reactor models such as the plug flow reactor and the partially stirred reactor were used to investigate the role of mixing in the fluidized bed pyrolyzer and to establish the impact of mixing time on the gas-phase reactions. A comprehensive multiphase computational fluid dynamics (mCFD) framework, including polydisperse granular phase modeling and detailed chemical kinetics was formulated and used to model the experimental setup for autothermal biomass fast pyrolysis at ISU. Multiphase CFD was also used to investigate the role of biomass feed positioning on the mixing of biomass in the pyrolizer. Finally, a reduced order model (ROM), suitable to be implemented in process simulators was obtained. Both the mCFD and the ROM were validated against experiments.

09 BIOMASS FUELS↗

DOE OSTI2022

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Solar Receiver with Integrated Thermal Storage for a Supercritical Carbon Dioxide Power Cycle

This project endeavors to design and demonstrate an integrated CSP system that meets the Department of Energy’s 6¢/kW-hr LCOE goal [1]. This design is based on the 10 MWe Supercritical Transformational Electric Power (STEP) engine specifications, but adaptable to any sCO 2 power module up to 10 MW e . The enabling technology to be tested is a new class of metal hydride based thermal energy storage materials which can store thermal energy in the hydrogen bonds between the metal hydride and hydrogen. This thermochemical energy storage solution provides a much greater energy density than sensible and latent energy storage enabling an integrated system design that closely linked the metal hydride to the power block. This allows the entire system to be mounted up tower to provide cost reductions. A novel metal hydride, capable of operating above 760 °C, was developed during the course of this project and thoroughly tested at a lab scale. The material’s hydrogen capacity, reaction enthalpy, thermal conductivity, and cyclic stability were measured. A demonstration test article was built that can store 2kWh of thermal energy and was coupled to Brayton Energy’s CO 2 loop where the performance of the metal hydride under commercial conditions was measured. An integrated model was created that uses the inputs from all the different sub-models in order to return an overall system LCOE using SAM. This project greatly advanced metal hydride based thermochemical energy storage technology and demonstrated its use. Previously there were no metal hydrides which were capable of operating at the proposed temperature limits (760°C), there is now a high temperature metal hydride with documented performance/stability and a path for low-cost production. Additionally, the metal hydride was tested at a scale >650x larger than what had previously been done in the lab. While not at commercial scale, it was an excellent demonstration of the technology and provided knowledge that can be leveraged for other applications.

08 HYDROGEN↗

DOE OSTI2022

Liquid and Gamma-BCC Phase Equilibria in the Uranium-Niobium System: A Simple Thermodynamic Model for Phase Diagram and Other Predictions

A simple thermodynamic model for the uranium-niobium (U-Nb) system was created using a Calphad-type approach, circa 2002. Only the liquid and gamma-BCC phases were modeled, with their subregular solution mixing terms optimized by trial and error. The custom (user-defined) database file (.tdb) and experimental data file (.exp) are listed, both of which are compatible with Thermo-Calc software. This model is used to calculate the solidus, liquidus, miscibility gap on the U-Nb phase diagram and also the enthalpy of melting in the U-6Nb alloy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

M4SF-22LL010301052-SUPCRTNE Thermodynamic Database Integration

This progress report (Level 4 Milestone Number M4SF-22LL010301052) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Argillite Activity Number SF-22LL01030105. SUPCRTNE (starting with v. 1.0) is being developed as the primary engine for future thermodynamic database development. It replaces SUPCRT92 (Johnson et al., 1992), which was earlier used on the Yucca Mountain Project (Wolery and Jove-Colon, 2007). Various changes to SUPCRT92 were and are presently being made. A major change is replacing the original water properties module H2O92D. For this, we have developed a modularized form of the program H2OI95. H2OI95 (Wolery, 2018, 2020) evaluates a newer equation of state for water (IAPWS-95, (Wagner and Pruß, 2002)) that is the present international standard. Thermochemical properties including Gibbs energy, enthalpy, and entropy are now tied to the CODATA standard (Cox et al., 1989).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

In-situ thermodynamics measurements at metal oxides-solution interfaces using Flow Adsorption Microcalorimetry.

Mineral-fluid interfaces are the principal sites of geochemical processes near Earth’s surface, hosting chemical reactions that play a fundamental role in (bio)geochemical cycles and in the fate and transport of anthropogenic contaminants, and that effectively control the compositions of soil and water environments. These complex interfaces are critical for our energy and environmental future. The mineral-fluid interface has been studied in an unprecedented level of detail with both experimental and computational approaches, separately and in combination. However, conspicuously missing from studies of the mineral-fluid interface are direct measurements of the energies of ion sorption and exchange. The literature on energetics and enthalpies of exchange, adsorption, dissolution, precipitation, and surface protonation reactions, especially those directly supported by experimental data, remains scarce despite their fundamental nature. The overarching goal of this project is to complete a systematic study of the thermodynamics properties of interfacial reactions at four MO surfaces (Rutile (α-TiO2), Quartz (SiO2), boehmite (γ-AlOOH) and goethite (α-FeOOH)) through the application and construction of novel flow adsorption microcalorimetry techniques and instrumentations. These unique and specialized microcalorimeters will operate at various temperatures and solution chemical compositions allowing for in-situ measurements across metal oxides and ligands of various characteristics. The overall research goal will be accomplished by completing the following three specific objectives (O): O1) Determine the energetics of surface protonation and deprotonation, ion exchange and ligand sorption reactions; O2) Investigate the surface charge thermodynamic properties under a range of temperature and solution chemical compositions; and O3) Develop predictive trends of the interplay between MO structure, surface coverage and surface reactivity. In addition to key thermodynamics parameters, calorimetric measurements provide a wealth of mechanistic information about reactions energetics and kinetics, surface charge characteristics, and structure-reactivity or selectivity relationships, all obtained in-situ and in real-time. This report includes science highlights from various projects completed over the performance period of the project. Also listed are dissemination opportunities, people supported on the grant and the impact on available physical resources and the discipline as a whole.

58 GEOSCIENCES↗

DOE OSTI2023

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Overlapping Domain Coupling of Multidimensional and System Codes in NEAMS - Pronghorn and SAM

This report describes the implementation and testing of domain-overlapping coupling of the Pronghorn and SAM codes. Both Pronghorn and SAM are codes developed by the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program for analysis of advanced nuclear reactors. Pronghorn focuses on analyzing multi-dimensional core flow conditions, while SAM focuses on analysis of the entire system including, among other components, piping, pumps, and heat exchangers. It is desirable for many advanced reactor thermal-hydraulics simulations to couple Pronghorn and SAM to obtain self-consistent solutions in the core and the system. The domain-overlapping coupling approach provides a robust and numerically efficient candidate for this coupling. We describe the algorithm for obtaining self-consistent solutions between a SAM network of flow channels and an overlapped multi-dimensional Pronghorn domain. The algorithm is tested for a variety of simple canonical test problems and a model of the Molten Salt Reactor Experiment (MSRE). Results demonstrate that the domain overlapping coupling approach provides consistent pressure, mass flow rate, enthalpy, and passive scalar distributions between Pronghorn and SAM and converges reliably and efficiently (at most 15 iterations, usually much less than 10) for the considered test problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Bio-fuel effects and seasonal fuel property variations related to abnormal combustion occurrences observed in the field

Since 2015, the automotive industry has regularly observed sharp increases in customer complaints of poor driveability and increased warranty claims during the months of September and October. These persistent seasonal complaint increases have been tracked to occur from increased abnormal combustion events that coincide with the seasonal transition from summer to winter fuel volatility, beginning on September 15th annually. The Reid Vapor Pressure (RVP) of gasoline during this seasonal transition has been found to be the primary property attributed to this annual issue. It’s hypothesized that an increase in volatile butanes and pentanes from advanced petroleum recovery techniques and hydrofracking production has hastened the transition from summer to winter fuel volatility in the marketplace, rather than the gradual increase of the RVP over time as historically intended. At worst, mid-winter fuel, with an RVP ranging from 12 to 16 psi, depending on location in the U.S., enters the market while ambient temperatures can still be above average for the winter season. This combination of high volatility fuel with unseasonably warm weather will increase fuel system flash boiling, injector spray cavitation, and spray collapse inside the engine’s combustion chamber, ultimately linking fuel property variations to abnormal engine combustion processes. This CRADA program will investigate these known knowledge gaps by coupling with two national laboratories, Oak Ridge National Laboratory, and the National Renewable Energy Laboratory to employ both unique tools and expertise throughout the industry to provide the insight, knowledge, and data critical to understand the observed interactions. Overall findings of this work are that fuel wall retention increased SPI proportionally, and fuels with an increased distillation (i.e. energy fraction) above the wall temperature increased fuel retention. Likewise fuels with increased enthalpy of vaporization (HoV) increased fuel retention proportional to the temperature reduction from the HoV (~13K reduction in this work). Thus, fuel with increased HoV and increased fuel distillation above the linear temperature also exhibited increased SPI rates. Fuel with high RVP did exhibit increased fuel spray collapse form spray imaging measurements, but in SPI testing the high volatility of these fuels at the tested conditions did not necessarily increase SPI rates despite exhibiting increased spray collapse in injector imaging and increased spray-wall impingement measured in fired engine testing. It was also found that spray imaging showed collapse was possible at throttled conditions, but this was not relevant to high boost pressures relevant to SPI conditions, however, under transient load and or speed encountered in real-world operation fuel spray collapse could affect fuel retention and real-world SPI that were not probed in the automated high load test sequence of this work.

09 BIOMASS FUELS↗

DOE OSTI2024

Predicted thermophysical properties of UN, PuN and (U,Pu)N

Molecular dynamics and density functional theory simulations are used to predict the lattice and electronic contributions of thermophysical properties for UN, PuN, and mixed (U,Pu)N systems. The properties predicted include the lattice parameter, linear thermal expansion, enthalpy, and specific heat capacity, as a function of temperature. The simulation predictions for high temperature specific heat capacity are compared against experimental measurements to understand the behavior, and why differences in the experimental measurements are observed. The influence of adding U vacancies, N interstitials, and Pu to UN is also examined. For this, a new PuN potential parameter set is developed and used with the Kocevski UN potential, enabling the dynamics of mixed (U,Pu)N systems to be studied. How defects impact the thermophysical properties is important for understanding fuel behavior under different reactor conditions, and these mechanistic predictions can be used to support fuel performance codes where data is scarce.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Impact of intermolecular interactions on the spectroscopic signals, energetics, and redox behavior of high valent 237 Np (Final Report)

Our studies have led to definitive spectral assignments for neptunyl-neptunyl interactions in aqueous solutions that will be valuable for predicting Np(V) speciation in solution. We have definitively demonstrated how actinyl-hydrogen and actinyl-cation interactions can influence spectral features for uranyl and neptunyl system. In addition, we have developed DFT methodologies to predict enthalpies of formation for uranyl and neptunyl systems and linked stability to chemical descriptors (hydrogen interaction energies or the electrostatic attraction energies). The counterion influenced the rate of Np(V) oxidation to Np(VI) neptunates and [Np(VII)O 4 OH] 3- complexes. Oxidation to Np(VII) during ozonolysis takes place through the presence of hydroxyl radicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Density functional theory (DFT) study of UF 6 hydrolysis: reaction pathways, spectroscopy, and chemical kinetics

Depleted uranium hexafluoride (UF 6 ), a stockpiled byproduct of the nuclear fuel cycle, reacts readily with atmospheric humidity, but the gas-phase reaction mechanism and associated chemical kinetics are poorly understood. During the performance period we undertook development of a state-of-the-art ab initio gas-phase chemical kinetics simulation workflow to model the hydrolysis of uranium hexafluroride (UF 6 ). In doing so, we addressed several outstanding issues in the theoretical treatment of uranium-containing systems. At the outset it was unclear how to generate accurate estimates of kinetic and thermodynamic data for U-containing chemical reactions. Generation of such data has been made routine. Prior to our work, the literature associated with UF 6 hydrolysis were disparate and inaccurate. This body of work provides a modern and comprehensive theoretical assessment of the reaction mechanism, molecular clustering towards deposition, and chemical kinetics. New methodological implementations and software integrations resulting from this work are also highlighted. As much as possible, our predictions were validated against experimental data including particle morphologies, vibrational spectroscopy, atomization enthalpies, and kinetic rate constants. Nevertheless, we were unable to reconcile kinetic measurements with high-accuracy simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Metal Organic Frameworks for Noble Gas Isotope Harvesting at FRIB (Final Technical Report)

This project was a collaborative effort between Lawrence Livermore National Laboratory (LLNL) and Michigan State University (MSU) to investigate the use of promising metal organic frameworks (MOFs) for radioactive noble gas capture, with a focus on harvesting exotic radiokryptons from FRIB. After screening several candidate materials, two MOFs were selected for testing: SIFSIX-3Cu and SBMOF-1. Further evaluation showed that although SIFSIX-3Cu has a high selectivity for Kr, SBMOF-1 is less sensitive to the humidity that is present in the FRIB harvesting system and is more readily integrated into the harvesting infrastructure. SBMOF-1 was then evaluated for temperature-dependent Kr and Xe uptake in order to determine the sorption enthalpy. The SBMOF-1 data led to the design of a noble gas capture system that will be fabricated and put into service for isotope harvesting at FRIB as part of a separate project.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Noble Gas Transport in the MSRE

This study explores the relationship between the physicochemical parameter known as Henry’s gas constant and gas transport across a two-layer film interface. The investigation utilized the Gibbs free energy, incorporating surface and volume terms to elucidate trends in enthalpy and entropy. Notably, our findings align with experimental data and offer predictive insights into the Henry’s gas constants for helium and krypton, which hold significance for future experiments and theoretical developments. Furthermore, this study enhances the Gibbs free energy theory pertaining to the liquid–gas interface. It underscores the substantial contribution of noble gases in this region to volumetric energy as temperature increases. Additionally, we employed Monte Carlo simulations to analyze the effective thermal neutron multiplication factor, denoted as k eff . Our analysis reveals a linear correlation between graphite density and uniform density as a function of temperature. For the 1D Molten Salt Reactor Experiment (MSRE) system, we employed the Mole code to conduct heat and mass transfer calculations. These computations enable us to ascertain the distribution of fuel temperature based on coefficients and thermal properties. We also studied delayed neutron precursors during fuel cycling, taking into account the drift of cycling fuel through Mole–Griffin coupling. Our model represents k eff and β eff across various volume flow rates and salinity compositions. Finally, this study leveraged xenon-135 for continuous on-line monitoring of fuel salts to investigate the impact of steady-state xenon-135 on the MSRE and to better understand its distribution. These efforts build upon previous research related to removal processes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Measurement of Thermal Conductivity in a Supercooled Hydrogel-Salt Complex Near Its Phase Transition

Solid-liquid phase transitions, i.e. solidification processes, have applications in data storage, development of novel thermoelectric materials, cooling of microelectronic substrates and air conditioning condensers. Standard analyses of solidification (Stefan problem) assume constant thermal properties of the solid and liquid sides. It is not known how these properties change across the spatial transition interface, though most studies report a discontinuity in the solid and liquid properties through the transition temperature [1]. When the phase transition releases enthalpy, recent research has shown if the phonon or electron transport time is of the same order of magnitude as the time scale of the atomic transformation, this increases the heat capacity of the solid material at temperatures near the phase transition temperature [2]. A fundamental understanding of the phase transition may help shed light on the molecular origins of supercooling and spontaneous nucleation, which will help with applications involving them. We report studies of supercooling and nucleation of sodium sulfate decahydrate, a salt hydrate which is of recent interest in thermal storage, and of a hydrogel-sodium sulfate complex which shows limited supercooling. We will also report the measurement of their thermal properties during the phase transition process. This will be done with a hot-wire setup which produces small temperature changes of the order of ~1 C.

crystallization, phase transformation, thermal con↗