Engineering Papers⌕ Search

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.

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

At least 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Nitrogen thermometry in an inductively coupled plasma torch using broadband nanosecond coherent anti-Stokes Raman scattering

The development of atmospheric hypersonic flight and re-entry capabilities requires the characterization of the thermo-chemical state of representative test environments. This study demonstrates the usage of multiplex nanosecond N 2 coherent anti-Stokes Raman scattering (CARS) to measure temperatures in an atmospheric, high-temperature (>6000K), air plasma plume, generated by an inductively coupled plasma torch. These are some of the highest temperatures ever accessed via gas-phase CARS, to our knowledge. Temperatures of N 2 in the equilibrium plasma plume are determined via theoretical fits to measured CARS spectra. We discuss the practical implementation of CARS at very high temperatures, including the scaling of the N 2 CARS signal strength from 300 to 6700 K, where the expected peak signal from the high-temperature plasma torch gases is two orders of magnitude less than commonly encountered in combustion environments. An intensified CCD camera enables single-laser-shot detection at temperatures as high as 6200 K, by increasing sensitivity and providing a time gate against intense background luminosity. We also discuss the impacts of unwanted two-beam CARS contributions from outside the nominal three-beam measurement volume. We present mean axial and radial temperature profiles, as well as time-series data derived from both single-laser-shot and accumulated CARS spectra. The single-laser-shot precision is 1.7%–2.6% at temperatures of 3500 to 6200 K. The presented results pave the way for the use of CARS at very high temperatures and the measurement of spatially resolved interface processes in high-enthalpy flows.

Fries, Dan (ORCID:0000000229513294)↗
DOE OSTI2025

Nitric oxide molecular tagging velocimetry of a free-flight model in a reflected shock tunnel

Nitric oxide molecular tagging velocimetry (NO MTV) is used to characterize the wake behind a free-flight spherical model in the test section of a free-piston reflected shock tunnel using a burst-mode laser operated at 100 kHz. A novel (to our knowledge) multi-delay timing scheme was implemented to measure velocity in a varying collisional environment in the wake of the free-flight model. Four simultaneous velocity profiles were measured in the wake of the model from -600 to 3600 m/s for flow enthalpies of 10.3 and 12.0 MJ/kg. Finally, the measured velocity distributions show good agreement when compared to computational fluid dynamics (CFD) modeling.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Wurtzite materials in alloys of rock salt compounds

Materials with crystal structures containing tetrahedral motifs are preferable for optoelectronic applications because they often have direct band gaps and low electron effective masses. However, crystal structures of manganese chalcogenides typically contain octahedral motifs, such as in rocksalt (RS) MnS and MnSe materials. Here we experimentally show that MnS 1-x Se x alloys with tetrahedrally-bonded wurtzite (WZ) structure can form between MnSe and MnS parent compounds with octahedral RS structures, at S-rich compositions (x<0.4) and low synthesis temperatures (~300°C). The calculated mixing enthalpies of MnS1-xSex alloys in RS and WZ structures cannot explain this experimental observation, so we hypothesize that WZ stabilization may be related to smaller structure density and lower surface energy compared to RS. The resulting WZ MnS 1-x Se x alloys have 3.0-3.2 eV optical absorption onset and lower electrical conductivity (<0.0001 S/cm) compared to the parent RS compounds. These experimental measurement result are consistent with computationally predicted band gaps and effective masses.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Theoretical Multinary Perovskite Oxides Dataset

Perovskite oxides (ternary chemical formula ABO3) are a diverse class of materials with applications including heterogeneous catalysis, solid-oxide fuel cells, thermochemical conversion, and oxygen transport membranes. However, their multicomponent chemical space is underexplored due to the immense number of possible compositions. To expand the number of computed compounds we report a dataset of 66,516 theoretical multinary oxides, 59,708 of which are perovskites. First, the AA'BB'O6 compositions were generated in the a-b+a- Glazer tilting mode using the computationally-inexpensive Structure Prediction and Diagnostic Software (SPuDS) program. Next, we optimized these structures with DFT using parameters compatible with the Materials Project (MP) database. Our dataset contains these optimized structures and their formation (ΔHf) and decomposition enthalpies (ΔHd) computed relative to MP tabulated elemental references and competing phases, respectively. This dataset can be mined, used to train machine learning models, and rapidly and systematically expanded by optimizing more SPuDS-generated AA'BB'O6 perovskite structures using MP-compatible DFT calculations

Bare, Zachary↗
DOE OSTI2022

Hydration-driven stabilization and volume collapse of grain boundaries in Mg2SiO4 forsterite predicted by first-principles simulations

Grain boundaries in mantle minerals are of critical importance to geophysical and geochemical processes of the Earth’s interior. One of the fundamental issues is to understand how the water (H 2 O) component influences the properties of grain boundaries in silicate materials. Here, we report the results of the structure and stability of several tilt grain boundaries in Mg 2 SiO 4 forsterite over the pressure range 0 to 15 GPa using density functional theory-based first-principles simulations. Here, the results suggest greater energetic stability and hydration-driven volume collapse (negative excess volume) at zero pressure for the majority of hydrous grain boundaries relative to the anhydrous (dry) ones. All the hydrous grain boundaries become increasingly favorable at elevated pressures as the calculated hydration enthalpy systematically decreases with increasing pressure. The hydrous components at the interfacial regions are predominantly in the hydroxyl form and, to a lesser extent, in the molecular H 2 O form. Their calculated ratio ranges from 1.6 to 8.7 among the different grain boundary configurations. Our structural analysis also reveals that the hydroxyls are bound to either both Mg and Si or to Mg only. In comparison, the molecular species are bound only to Mg sites. Besides direct oxygen-hydrogen bonding, intermolecular hydrogen bonding becomes important with compression. On the basis of our results, we suggest that local atomic rearrangements caused by dissociative adsorption of water facilitate efficient compaction of the boundary interfaces, which, in turn, results in greater relative stability of hydrous grain boundaries. This means that water prefers to be incorporated within the grain boundaries over the bulk of silicate materials.

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

Kinetics of TATB and LX-17 Decomposition by Thermal Analysis

The kinetics of TATB decomposition have been investigated by thermal analysis (thermogravimetric analysis, TGA, and differential scanning calorimetry, DSC) at both LLNL and the University of Alabama Birmingham. This report summarizes progress in understanding from those results within the context of identifying the global reaction mechanism, i.e., sequential, autocatalytic, competitive, or some mixture thereof. A particularly important issue is how pressure affects both the reaction rate and enthalpy. The effect of pressure helps identify dominant mechanisms and construct and constrain a global mechanism that will work over a wide range of conditions. This memorandum is simply a snapshot of current progress and not a definitive interpretation. I have written it to organize work of others and my interpretations in a systematic manner with the intention of stimulating further advances in understanding. SDT data comes from Jen Montgomery at LLNL, and high-pressure DSC comes from Victoria Stanford and Sergey Vyazovkin at University of Alabama at Birmingham. Major observations and conclusions, some new and some not, are (1) outgassing below 150 °C prior to decomposition releases adsorbed water vapor, (2) outgassing between 150 and 250 °C may come primarily from intercalated water (3) initial thermal decomposition above 250 °C or so involves dehydration of TATB to create furazans, (4) quickly heating to higher temperatures in open pans with flowing gas leads primarily to TATB evaporation between 270 and 390 °C, (5) using a pinhole crucible at atmospheric pressure largely eliminates TATB sublimation, (6) at typical thermal analysis heat rates, the main stage of decomposition is described by a single mass loss process, and decomposition of residual material is spread over a wide range of thermal severity, (7) at typical thermal analysis heat rates and constant temperatures above 300 °C, heat release occurs in two overlapping peaks, with the first one a little faster than mass loss and the second on a little slower than mass loss, (8) total heat release increases significantly with heating rate and substantially with increased pressure, and (9) isothermal experiments at 340 °C indicate that the first exotherm increases by more with LX-17 than TATB, which suggests a possible chemical role of Kel-F decomposition.

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

Fluorinated Covalent Organic Frameworks: A Novel Pathway to Enhance Hydrogen Sorption and Control Isosteric Heats of Adsorption; HyMARC Seed Project Final Report

Organic framework materials, as a unique sub-class of carbon-based sorbents, have gained increasing attention for promising attributes toward gas storage. Calculations of the hydrogen capacity of optimized frameworks show potential to achieve greater than 60 g/L storage of hydrogen, placing metal-organic frameworks (MOFs) and COFs near the top of the class of porous materials. However, at the ensemble level, poor stacking creates a quasi-amorphous material with low structural integrity and low effective surface area. The chemical versatility of COFs allows for additional methods for producing long-range order that specifically target the inter-layer interactions in COFs but that leave key pore-accessible linker sites open. These modifications (either during synthesis or post-synthetic) enable a host of strategies to both improve crystalline order for more stable and higher surface, as well as add metals with open coordination for enhanced H 2 binding enthalpy.

08 HYDROGEN↗
DOE OSTI2020

Ab Initio Phase Diagram of Tungsten

The phase diagram of tungsten (W) to a pressure (P) of 2500 GPa is investigated using a comprehensive ab initio approach that includes (i) the calculation of the zero temperature (T) free energies (enthalpies) of different solid structures, (ii) the quantum molecular dynamics simulation of the melting curves of different solid structures, (iii) the derivation of the analytic form for the solid-solid phase transition boundary, and (iv) the simulations of the solidification of liquid W into the final solid states on both sides of the solid-solid phase transition boundary, in order to confirm the corresponding analytic form. There are two solid structures confirmed to be present on the phase diagram of W, the ambient body-centered cubic (bcc) and the high-pressure double hexagonal close-packed (dhcp). At T = 0, the bcc-dhcp transition occurs at 1060 GPa, and the transition boundary has a positive slope dT/dP : the bcc-dhcp-liquid triple point is at (P, T) = (1675 GPa, 23680 K).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simulation of flame acceleration and deflagration-to-detonation transition with heat transfer in HE product and fractured HE

A simple heat transfer model has been developed for high explosives. For the gaseous high explosive (HE) product, the Smargorinsky eddy viscosity model has been adopted for estimating the turbulent thermal conductivity. Convective heat transfer in the fractured solid HE is modeled as the effective thermal conduction in the porous media, which is equivalent to the enthalpy flux carried by the velocity determined by Darcy’s law with Forchheimer term.Other features of the model include a multi-reaction chain, equation-of-state for each species involved, and a two-phase hotspot model. Test calculations are performed with various initial disturbances. Onset and progress of flame acceleration (FA) and chemical reactions are strongly affected by the heat transfer and initial disturbance. Staged development of FA also affirms the necessity of modeling based on a multi-reaction chain. When conditions warrant,e.g., stronger initial disturbance, FA evolves to deflagration-to-detonation transition (DDT),demonstrating automatic capture of the DDT.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
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↗