Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “THERMODYNAMICS”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Thermodynamic Activity Measurements with Knudsen Cell Mass Spectrometry

Coupling the Knudsen effusion method with mass spectrometry has proven to be one of the most useful experimental techniques for studying the equilibrium between condensed phases and complex vapors. The Knudsen effusion method involves placing a condensed sample in a Knudsen cell, a small "enclosure", that is uniformly heated and held until equilibrium is attained between the condensed and vapor phases. The vapor is continuously sampled by effusion through a small orifice in the cell. A molecular beam is formed from the effusing vapor and directed into a mass spectrometer for identification and pressure measurement of the species in the vapor phase. Knudsen cell mass spectrometry (KCMS) has been used for nearly fifty years now and continues to be a leading technique for obtaining thermodynamic data. Indeed, much of the well-established vapor specie data in the JANAF tables has been obtained from this technique. This is due to the extreme versatility of the technique. All classes of materials can be studied and all constituents of the vapor phase can be measured over a wide range of pressures (approximately 10(exp -4) to 10(exp -11) bar) and temperatures (500-2800 K). The ability to selectively measure different vapor species makes KCMS a very powerful tool for the measurement of component activities in metallic and ceramic solutions. Today several groups are applying KCMS to measure thermodynamic functions in multicomponent metallic and ceramic systems. Thermodynamic functions, especially component activities, are extremely important in the development of CALPHAD (Calculation of Phase Diagrams) type thermodynamic descriptions. These descriptions, in turn, are useful for modeling materials processing and predicting reactions such as oxide formation and fiber/matrix interactions. The leading experimental methods for measuring activities are the Galvanic cell or electro-motive force (EMF) technique and the KCMS technique. Each has specific advantages, depending on material and conditions. The EMF technique is suitable for lower temperature measurements, provided a suitable cell can be constructed. KCMS is useful for higher temperature measurements in a system with volatile components. In this paper, we briefly review the KCMS technique and identify the major experimental issues that must be addressed for precise measurements. These issues include temperature measurements, cell material and cell design and absolute pressure calibration. The resolution of these issues are discussed together with some recent examples of measured thermodynamic data.

Copland, Evan H.↗
NASA NTRS2004

Thermodynamic Vent System Test in a Low Earth Orbit Simulation

A thermodynamic vent system for a cryogenic nitrogen tank was tested in a vacuum chamber simulating oxygen storage in low earth orbit. The nitrogen tank was surrounded by a cryo-shroud at -40 F. The tank was insulated with two layers of multi-layer insulation. Heat transfer into cryogenic tanks causes phase change and increases tank pressure which must be controlled. A thermodynamic vent system was used to control pressure as the location of vapor is unknown in low gravity and direct venting would be wasteful. The thermodynamic vent system consists of a Joule-Thomson valve and heat exchanger installed on the inlet side of the tank mixer-pump. The combination is used to extract thermal energy from the tank fluid, reducing temperature and ullage pressure. The system was sized so that the tank mixer-pump operated a small fraction of the time to limit motor heating. Initially the mixer used sub-cooled liquid to cool the liquid-vapor interface inducing condensation and pressure reduction. Later, the thermodynamic vent system was used. Pressure cycles were performed until steady-state operation was demonstrated. Three test runs were conducted at tank fills of 97, 80, and 63 percent. Each test was begun with a boil-off test to determine heat transfer into the tank. The lower tank fills had time averaged vent rates very close to steady-state boil-off rates showing the thermodynamic vent system was nearly as efficient as direct venting in normal gravity.

VanOverbeke, Thomas J.↗
NASA NTRS2019

The Role of Thermodynamic Phase Shifts in Cloud Optical Depth Variations with Temperature

We present a novel method that identifies the contributions of thermodynamic phase shifts and processes governing supercooled liquid and ice clouds to cloud optical depth variations with temperature using MODIS observations. Our findings suggest that thermodynamic phase shifts outweigh the net influence of processes governing supercooled liquid and ice clouds in causing increases in mid-latitudinal cold cloud optical depth with temperature. Cloud regime analysis suggests that dynamical conditions appear to have less influence on the contribution of thermodynamic phase shifts to cloud optical depth variations with temperature. Thermodynamic phase shifts also contribute more to increases in cloud optical depth during colder seasons due to the enhanced optical thickness contrast between liquid and ice clouds. The results of this study highlight the importance of thermodynamic phase shifts in explaining cold cloud optical depth increases with temperature in the current climate and may elucidate their role in the cloud optical depth feedback.

cloud optical depth↗
NASA NTRS2021

Estimation of Thermodynamic and Dynamic Contributions to Sea Ice Growth in the Central Arctic Using ICESat-2 and MOSAiC SIMBA Buoy Data

The fine spatial resolution of the ICESat-2 (IS2) satellite altimeter allows monitoring the evolution of sea ice thickness with detailed dynamic information (e.g. ridges and leads). In this study, we first assess the ability of IS2 to estimate thermodynamic ice growth and dynamic thickening during the ice-growing season in the central Arctic Ocean. As an indicator of the thermodynamic ice growth, we use 10 thermistor string-based sea ice mass balance array (SIMBA) buoys deployed at a scale of ~ 50 km from the Icebreaker Polarstern during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition. We collect IS2 data within 20 km buffer distance from the individual buoys, and calculate the mode, median, and mean of the IS2-derived ice thickness. The IS2 modal thickness shows the least bias (−0.169 m) with the buoy ice thickness, representing level ice thickness. In addition, the increasing rate of the IS2 modal thickness is close to the thermodynamic ice growth with a small bias of −0.054 cm/day. However, the increasing rates of the IS2 median and mean thickness are greater than the thermodynamic ice growth by about 0.114 cm/day and 0.198 cm/day, respectively, because they also include ice growth caused by thickness redistribution during dynamic deformation. The dynamic contributions may account for 26.1 ± 10.3% and 34.4 ± 10.1% of the total increase of the IS2 median and mean thickness, respectively. Within a ~ 50 km radius area from the MOSAiC Central Observatory, IS2 measurements exhibit that the ridge fraction increased from <2% in November to ~ 4% in March (~0.029%/day of average increasing rate) and ridge height increased about 0.047 cm/day during the same period. However, lead formation does not show significant contributions to the dynamic ice thickening because leads are temporary features lasting only 2–3 days. Although there are considerable uncertainties in IS2 ice thickness estimation and IS2-buoy thickness comparison, this study emphasizes the importance of combining measurements by IS2 and SIMBA buoys to explain the regional sea ice mass balance with separating the thermodynamic and dynamic contributions.

satellite altimeter↗
NASA NTRS2023

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗
DOE OSTI2025

Inverse Thermodynamics: Designing Interactions for Targeted Phase Behavior

The traditional goal of inverse self-assembly is to design interactions that drive particles toward a desired target structure. However, achieving successful self-assembly also requires tuning the thermodynamic conditions under which the structure is stable. In this work, we extend the inverse design paradigm to explicitly address this challenge by developing a framework for inverse thermodynamics, i.e., the design of interaction potentials that realize specific thermodynamic behavior. As a step in this direction, using patchy particle mixtures as a model system, we demonstrate how precise control over both bonding topology and bond energetics enables the programming of targeted phase behavior. In particular, we establish design principles for azeotropic demixing and show how to create mixtures that exhibit azeotropy at any prescribed composition. Our predictions are validated through Gibbs-ensemble simulations [Panagiotopoulos, Mol. Phys. 1987, 61, 813−826]. These results highlight the necessity of coupling structural design with thermodynamic engineering, and provide a blueprint for controlling complex phase behavior in multicomponent systems.

Azeotropes↗
DOE OSTI2024

On Gibbs Equilibrium and Hillert Nonequilibrium Thermodynamics

During his time at Royal Institute of Technology (Kungliga Tekniska högskolan) in Sweden, the present author learned nonequilibrium thermodynamics from Mats Hillert. The key concepts are the separation of internal and external variables of a system and the definitions of potentials and molar quantities. In equilibrium thermodynamics derived by Gibbs, the internal variables are not independent and can be fully evaluated from given external variables. While irreversible thermodynamics led by Onsager focuses on internal variables though often mixed with external variables. Hillert integrated them together by first emphasizing their differences and then examining their connections. His philosophy was reflected by the title of his book “Phase Equilibria, Phase Diagrams and Phase Transformations” that puts equilibrium, nonequilibrium, and internal processes on equal footing. Here, in the present paper honoring Hillert, the present author reflects his experiences with Hillert and his work in last 40 years and expresses his gratitude for all the wisdom and support from him in terms of “Hillert nonequilibrium thermodynamics” and discusses some recent topics that the present author has been working on.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A thermodynamically consistent discretization of 1D thermal-fluid models using their metriplectic 4-bracket structure

Thermodynamically consistent models in continuum physics, i.e. models which satisfy the first and second laws of thermodynamics, may be expressed using the metriplectic formalism. In this work, we leverage the structures underlying this modeling formalism to preserve thermodynamic consistency in discretizations of a fluid model. The procedure relies (1) on ensuring that the spatial semi-discretization retains certain symmetries and degeneracies of the Poisson and metriplectic 4-brackets, and (2) on the use of an appropriate energy conserving time-stepping method. Here, the minimally simple yet nontrivial example of a one-dimensional thermal-fluid model is treated. It is found that preservation of the requisite symmetries and degeneracies of the 4-bracket is relatively simple to ensure in Galerkin spatial discretizations, suggesting a path forward for thermodynamically consistent discretizations of more complex fluid models using more specialized Galerkin methods.

Hamiltonian structure↗
DOE OSTI2024

Thermodynamically consistent Cahn–Hilliard–Navier–Stokes equations using the metriplectic dynamics formalism

Cahn–Hilliard–Navier–Stokes (CHNS) systems describe flows with two-phases, e.g., a liquid with bubbles. Obtaining constitutive relations for general dissipative processes for such systems, which are thermodynamically consistent, can be a challenge. We show how the metriplectic 4-bracket formalism (Morrison and Updike, 2024) achieves this in a straightforward, in fact algorithmic, manner. First, from the noncanonical Hamiltonian formulation for the ideal part of a CHNS system we obtain an appropriate Casimir to serve as the entropy in the metriplectic formalism that describes the dissipation (e.g. viscosity, heat conductivity and diffusion effects). General thermodynamics with the concentration variable and its thermodynamics conjugate, the chemical potential, are included. Having expressions for the Hamiltonian (energy), entropy, and Poisson bracket, we describe a procedure for obtaining a metriplectic 4-bracket that describes thermodynamically consistent dissipative effects. The 4-bracket formalism leads naturally to a general CHNS system that allows for anisotropic surface energy effects. Furthermore, this general CHNS system reduces to cases in the literature, to which we can compare.

Cahn–Hilliard↗
DOE OSTI2025

Proton, Electron, and Hydrogen-Atom Transfer Thermodynamics of the Metal–Organic Framework, Ti-MIL-125, Are Intrinsically Correlated to the Structural Disorder

Interfacial charge transfer reactions involving protons and/or electrons are fundamental to heterogeneous catalysis and many other reactions relevant to energy, chemical, and biological sectors. Metal–organic frameworks (MOFs) with redox-active metal-oxo nodes have emerged as candidate materials to examine these reactions with near-atomic-level precision, given their crystalline nature. Here, we employed a colloidally stable, Ti-based MOF, Ti-MIL-125, with different crystal sizes to examine catalytically relevant charge transfer thermodynamics. The Ti 8 (μ 2 -O) 8 (μ 2 -OH) 4 nodes structurally mimic TiO 2 , which has shown some PCET reactivity toward reactions of H 2 , O 2 , and others. In this report, we have demonstrated that a change in crystal size induces different amounts of structural disorder to the Ti-oxo node, further changing the thermodynamics of proton/electron/hydrogen-atom transfer reactions. Using electrochemical open-circuit potential (E OCP ) measurements, we have determined that all crystallites undergo a 1H + /1e – redox reaction, which, given the stoichiometry, can be considered as a net H atom transfer (HAT) reaction. The thermodynamics of this HAT reaction, the Ti 3+ O–H bond dissociation free energy (BDFE), was dependent on the crystal size of the MOF, as the decrease in crystal size induced more structural disorder. Our computational calculations have indicated that this difference in BDFE is due to a local change in the geometry of Ti cations, rather than the commonly invoked defects, such as the “missing-linker” defect sites. Individual proton/electron transfer (PT/ET) thermodynamics were also highly dependent on the crystal sizes. These were probed using pK a or band gaps (E g ), respectively. These findings suggest that, particularly when MOFs are nanosized with a large amount of structural disorder, they should no longer be considered “true” single-site catalysts; this is an implicit, but widespread assumption within the MOF-based catalysis field. Implications of these findings will be contrasted with structurally similar metal oxides like TiO 2 and other redox-active MOFs.

Bond dissociation free energy↗
DOE OSTI2024

Pyrophosphate-free glycolysis in Clostridium thermocellum increases both thermodynamic driving force and ethanol titers

Background: Clostridium thermocellum is a promising candidate for production of cellulosic biofuels, however, its final product titer is too low for commercial application, and this may be due to thermodynamic limitations in glycolysis. Previous studies in this organism have revealed a metabolic bottleneck at the phosphofructokinase (PFK) reaction in glycolysis. In the wild-type organism, this reaction uses pyrophosphate (PP i ) as an energy cofactor, which is thermodynamically less favorable compared to reactions that use ATP as a cofactor. Previously we showed that replacing the PP i -linked PFK reaction with an ATP-linked reaction increased the thermodynamic driving force of glycolysis, but only had a local effect on intracellular metabolite concentrations, and did not affect final ethanol titer. Results: In this study, we substituted PP i -pfk with ATP-pfk, deleted the other PPi-requiring glycolytic gene pyruvate:phosphate dikinase (ppdk), and expressed a soluble pyrophosphatase (PPase) and pyruvate kinase (pyk) genes to engineer PP i -free glycolysis in C. thermocellum. We demonstrated a decrease in the reversibility of the PFK reaction, higher levels of lower glycolysis metabolites, and an increase in ethanol titer by an average of 38% (from 15.1 to 21.0 g/L) by using PP i -free glycolysis. Conclusions: By engineering PP i -free glycolysis in C. thermocellum, we achieved an increase in ethanol production. These results demonstrate that optimizing the thermodynamic landscape through metabolic engineering can enhance product titers. While further increases in ethanol titers are necessary for commercial application, this work represents a significant step toward engineering glycolysis in C. thermocellum to increase ethanol titers.

09 BIOMASS FUELS↗
DOE OSTI2025

Thermodynamic Cloud Phase Classifications Using Machine Learning at NSA and ANX

Vertically resolved thermodynamic cloud phase classifications are essential for studies of atmospheric cloud and precipitation processes. The Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) THERMOCLDPHASE Value-Added Product (VAP) uses a multi-sensor approach to classify thermodynamic cloud phase by combining lidar backscatter and depolarization, radar reflectivity, Doppler velocity, spectral width, microwave radiometer-derived liquid water path, and radiosonde temperature measurements. The measured voxels are classified as ice, snow, mixed-phase, liquid (cloud water), drizzle, rain, and liq_driz (liquid+drizzle). We use this product as the ground truth to train three machine learning (ML) models to predict the thermodynamic cloud phase from multi-sensor remote sensing measurements taken at the ARM North Slope of Alaska (NSA) observatory: a random forest (RF), a multilayer perceptron (MLP), and a convolutional neural network (CNN) with a U-Net architecture. Evaluations against the outputs of the THERMOCLDPHASE VAP with one year of data show that the CNN outperforms the other two models, achieving the highest test accuracy, F1-score, and mean Intersection over Union (IOU). Analysis of ML confidence scores shows ice, rain, and snow have higher confidence scores, followed by liquid, while mixed, drizzle, and liq_driz have lower scores. Feature importance analysis reveals that the mean Doppler velocity and vertically resolved temperature are the most influential datastreams for ML thermodynamic cloud phase predictions. The ML models’ generalization capacity is further evaluated by applying them at another Arctic ARM site in Norway using data taken during the ARM Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) field campaign. Finally, we evaluate the ML models’ response to simulated instrument outages and signal degradation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Insights into controlling bacterial cellulose nanofiber film properties through balancing thermodynamic interactions and colloidal dynamics

In recent years, nanocellulose has emerged as a sustainable and environmentally friendly alternative to traditional petroleum-derived structural polymers. Sourced either from plants, algae, or bacteria, nanocellulose can be processed into colloid, gel, film and fiber forms. However, the required fundamental understanding of process parameters that govern the morphology and structure–property relationships of nanocellulose systems, from colloidal suspensions to bulk materials, has not been developed and generalized for all forms of cellulose. This further hinders the more widespread adoption of this biopolymer in applications. Our study investigates the dispersion of cellulose nanofibers (CNFs) produced by a bacterial–yeast co-culture, in solvents, highlighting the role of thermodynamic interactions in influencing their colloidal behavior. By adjusting Hansen solubility parameters, we controlled the thermodynamic relationship between CNFs and solvents across various concentrations, studying the dilute to semi-dilute regimes. Rheological measurements revealed that the threshold at which a concentration-based regime transition occurs is distinctly solvent-dependent. Complementing rheological analysis with small angle X-ray scattering and zeta potential measurements, our findings reveal that enhancing CNF–solvent interactions increases excluded volume in the dilute regime, emphasizing the importance of the balance between fiber–fiber and fiber–solvent interactions. Moreover, we investigated the transition from colloidal to solid state by creating films from dispersions with varying interaction parameters in semi-dilute regimes. Through mechanical testing and scanning electron microscopy imaging of the fracture surfaces, we highlight the significance of electrokinetic effects in such transitions, as dispersions with higher electrokinetic stabilization gave rise to stronger and tougher films despite having less favorable thermodynamic interaction parameters. Finally, our work provides insights into the thermodynamic and electrokinetic interplay that governs bacterial CNF dispersion, offering a foundation for future application and a deeper understanding of nanocellulose's colloidal and structure-property relationships.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Computing chemical potentials with machine-learning-accelerated simulations to accurately predict thermodynamic properties of molten salts

The successful design and deployment of next-generation nuclear technologies heavily rely on thermodynamic data for relevant molten salt systems. However, the lack of accurate force fields and efficient methods has limited the quality of thermodynamic predictions from atomistic simulations. Here we propose an efficient free energy framework for computing chemical potentials, which is the central free energy quantity behind many thermodynamic properties. We accelerate our simulations without sacrificing accuracy by using machine learning interatomic potentials trained on density functional theory (DFT) data. Using lithium chloride as our model system, we compute chemical potentials with DFT-accuracy for solid and liquid phases by transmuting ions into noninteracting particles. Notably, in the liquid phase, we demonstrate consistency whether we transmute one ion pair or the entire system into ideal gas particles. By locating the temperature where the chemical potential of solid and liquid phases cross, we predict a melting point of 880 ± 18 K for lithium chloride, which is remarkably close to the experimental value of 883 K. With this successful demonstration, we lay the foundation for high-throughput thermodynamic predictions of many properties that can be derived from the chemical potentials of the minority and majority components in molten salts.

Gibson, Luke D. [Oak Ridge National Laboratory (OR↗
DOE OSTI2024

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗
DOE OSTI2025

Thermodynamics of calcium binding to heparin: Implications of solvation and water structuring for polysaccharide biofunctions

Heparan sulfates are found in all animal tissues and have essential roles in living systems. This family of biomacromolecules modulates binding to calcium ions (Ca 2+ ) in low free energy reactions that influence biochemical processes from cell signaling and anticoagulant efficacy to biomineralization. Despite their ubiquity, the thermodynamic basis for how heparans and similarly functionalized biomolecules regulate Ca 2+ interactions is not yet established. Using heparosan (Control) and heparins with different positions of sulfate groups, we quantify how SO 3 − and COO − content and SO 3 − position modulate Ca 2+ binding by isothermal titration calorimetry. The free energy of all heparin-Ca2+ interactions (ΔG rxn ) is dominated by entropic contributions due to favorable water release from polar, hydrophilic groups. Heparin with both sulfate esters (O-SO 3 − ) and sulfamides (N-SO 3 − ) has the strongest binding to Ca 2+ compared to heparosan and to heparin with only O-SO 3 − groups (~3X). By linking Ca 2+ binding thermodynamics to measurements of the interfacial energy for calcite (CaCO 3 ) crystallization onto polysaccharides, we show molecule-specific differences in nucleation rate can be explained by differences in water structuring during Ca 2+ interactions. A large entropic term (-TΔS rxn ) upon Ca 2+ –polysaccharide binding correlates with high interfacial energy to CaCO 3 nucleation. Combining our measurements with literature values indicates many Ca 2+ –polysaccharide interactions have a shared thermodynamic signature. The resulting enthalpy–entropy compensation relationship suggests these interactions are generally dominated by water restructuring involving few conformational changes, distinct from Ca 2+ –protein binding. Our findings quantify the thermodynamic origins of heparin-specific interactions with Ca 2+ and demonstrate the contributions of solvation and functional group position during biomacromolecule-mediated ion regulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

M3SF-24LL010302062-NEA-TDB Management and International Collaborations in Sorption and Thermodynamic Modeling

This progress report (Level 3 Milestone Number M3SF-24LL010302062) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Crystalline International Collaborations Work Package Number SF-24LL01030206. The activity is focused on our long-term commitment to engaging our partners in international nuclear waste repository research. This includes participation in the Nuclear Energy Agency Thermochemical Database (NEA-TDB) Project and development of methodologies for integrating US and international thermodynamic databases for use in SFWST Generic Disposal System Assessment (GDSA) efforts. A continuing focus for FY24 efforts is to support the US participation in the NEA-TDB effort. The focus of FY24 activities was the development of an agreement for a Phase 7 activity that will start in Q1 of 2025. Mavrik Zavarin is now the US representative on both the Management Board and the Executive group to the NEA-TDB. He is also the POC for the Cements State of the Art Report that is undergoing peer review in FY24. In FY24, we used our position on the NEA-TDB MB and EG to facilitate the integration of NEA-TDB thermochemical data with LLNL’s SUPCRTNE thermodynamic database that supports the SFWST GDSA activities. This effort is coordinated with the Argillite work package SUPCRTNE database development efforts (Wolery, 2024). The goal is to provide a downloadable database that will be hosted on LLNL’s thermodynamics website which incorporates NEA-TDB data into the LLNL database where appropriate. We also began engagement with the EURAD2 program that was initiated in FY24 by our European collaborators at the Karlsruhe Institute of Technology (KIT), Germany. The primary focus of the engagement is with WP20: DITUSC Thermodynamic database evaluation program. A kickoff meeting for this activity is planned for early FY25. Finally, we have been selected to co-host (with Clemson University) the International Conference on Chemistry and Migration Behaviour of Actinides and Fission Products in the Geosphere in 2025 (Migration2025). The meeting will be held September 21-26, 2025, in New Orleans, Louisiana, and will focus on international efforts to understand the risks of radionuclide releases into the environment. This central focus of this conference is on international efforts to develop safe disposal options for nuclear wastes. As such, we are developing a theme focused on US underground nuclear waste repository science.

58 GEOSCIENCES↗
NASA NTRS1972

Thermodynamic properties of oxygen and nitrogen III

The final equation for nitrogen was determined. In the work on the equation of state for nitrogen, coefficients were determined by constraining the critical point to selected critical point parameters. Comparisons of this equation with all the P-density-T data were made, as well as comparisons to all other thermodynamic data reported in the literature. The extrapolation of the equation of state was studied for vapor to higher temperatures and lower temperatures, and for the liquid surface to the saturated liquid and the fusion lines. A new vapor pressure equation was also determined which was constrained to the same critical temperature, pressure, and slope (dP/dT) as the equation of state. Work on the equation of state for oxygen included studies for improving the equation at the critical point. Comparisons of velocity of sound data for oxygen were also made between values calculated with a preliminary equation of state and experimental data. Functions for the calculation of the derived thermodynamic properties using the equation of state are given, together with the derivative and integral functions for the calculation of the thermodynamic properties using the equations of state. Summary tables of the thermodynamic properties of nitrogen and oxygen are also included to serve as a check for those preparing computer programs using the equations of state.

Stewart, R. B.↗