Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamics of computing”

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 181 records · Page 10

Defect Diffusion Graph Neural Networks for Materials Discovery in High-Temperature Energy Applications

Here, the migration of crystallographic defects dictates material properties and performance for a plethora of technological applications. Density functional theory (DFT)-based nudged elastic band (NEB) calculations are a powerful computational technique for predicting defect migration activation energy barriers, yet they become prohibitively expensive for high-throughput screening of defect diffusivities. Without introducing hand-crafted (i.e., chemistry- or structure-specific) descriptors, we propose a generalized deep learning approach to train surrogate models for NEB energies of vacancy migration by hybridizing graph neural networks with transformer encoders and simply using pristine host structures as input. With sufficient training data, computationally efficient and simultaneous inference of vacancy defect thermodynamics and migration activation energies can be obtained to compute temperature-dependent vacancy diffusivities and to down-select candidates for more thorough DFT analysis or experiments. Thus, as we specifically demonstrate for potential water-splitting materials, candidates with desired defect thermodynamics, kinetics, and host stability properties can be more rapidly targeted from open-source databases of experimentally validated or hypothetical materials.

14 SOLAR ENERGY↗

Adsorbate Partition Functions via Phase Space Integration: Quantifying the Effect of Translational Anharmonicity on Thermodynamic Properties

A new method for computing anharmonic thermophysical properties for adsorbates on metal surfaces is presented. Classical Monte Carlo phase space integration is performed to calculate the partition function for the motion of a hydrogen atom on Cu(111). Here, a minima-preserving neural network potential energy surface is used within the integration routine. Two different sampling schema for generating the training data are presented, and two different density functionals are used. The results are benchmarked against direct state counting results by using discrete variable representation. The phase space integration results are in excellent quantitative agreement with the benchmark results. Additionally, both the discrete variable representation and the phase space integration results confirm that the motion of H on Cu(111) is highly anharmonic. The results were applied to calculate the free energy of dissociative adsorption of H 2 and the resulting Langmuir isotherms at 400, 800, and 1200 K in a partial pressure range of 0–1 bar. It shows that the anharmonic effects lead to significantly higher predicted surface site fractions of hydrogen.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ORNL/Thermo4PFM

Library to compute equilibrium compositions and equal chemical potential compositions in binary and ternary alloys. This library supports thermodynamic data formulated according to CALPHAD (Computer Coupling of Phase Diagrams and Thermochemistry) models, as well as dilute binary alloys. One major aim is to solve the Kim-Kim-Suzuki (KKS) equations to obtain local compositions given a local mixture of two phases. These quantities are then used to evaluate the phase-field model driving force (right handside of time-evolution equation), as well as the energy density

Fattebert, Jean-Luc↗

Efficiently predicting pressure-composition-temperature diagrams to discover low-stability metal hydrides

Quantitatively accurate computational predictions of metal hydride thermodynamics are challenging but critical for alloy performance optimization across a multitude of technological domains, including hydrogen storage, compression, purification, and getters. Recent machine learning approaches have demonstrated great success in this area, but can potentially suffer from several shortcomings since they rely on imbalanced experimental training data and can have poor out-of-distribution (ood) test performance. Here, in this study, we circumvent such pitfalls by developing a computationally efficient, first principles-based workflow for direct prediction of metal hydride phase equilibrium, i.e., the pressure-composition-temperature (PCT) diagram. We then demonstrate its utility on predicting low stability hydrides derived from compositionally complex C14 Laves phase AB2 alloys. Specifically, we computationally predict and then experimentally validate an AB 2 alloy series (z < 0.6 for Ti 2−z Zr z CrMnFeNi) with ideal hydriding thermodynamics for a two-stage metal hydride-based compressor for pressurizing boil off from liquefied hydrogen. Importantly, this study lays the groundwork for accurate and efficient discovery/optimization of ood, low-stability hydrides for which purely data-driven approaches lack sufficient accuracy.

08 HYDROGEN↗

Non-Electricity Based Renewable Fuels: Theory and Computation for Solar Thermochemical Hydrogen

Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox active defects start to interact with each other, thereby impeding the formation of additional defects. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.

first-principles calculations↗

Two-loop HTL-resummed thermodynamics for $\mathcal{N}$ = 4 supersymmetric Yang-Mills theory

We compute the two-loop hard-thermal-loop (HTL) resummed thermodynamic potential for \( \mathcal{N} \) = 4 supersymmetric Yang-Mills (SYM). Our final result is manifestly gauge-invariant and was renormalized using only simple vacuum energy, gluon mass, scalar mass, and quark mass counter terms. The HTL mass parameters m D , M D , and m q are then determined self-consistently using a variational prescription which results in a set of coupled gap equations. Based on this, we obtain the two-loop HTL-resummed thermodynamic functions of \( \mathcal{N} \) = 4 SYM. We compare our final result with known results obtained in the weak- and strong-coupling limits. We also compare to previously obtained approximately self-consistent HTL resummations and Padé approximants. We find that the two-loop HTL resummed results for the scaled entropy density is a quantitatively reliable approximation to the scaled entropy density for 0 ≤ λ ≲ 2 and is in agreement with previous approximately self-consistent HTL resummation results for λ ≲ 6.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Favorable Redox Thermodynamics of SrTi 0.5 Mn 0.5 O 3–δ in Solar Thermochemical Water Splitting

Two-step, solar thermochemical splitting of water using nonstoichiometric redox-active metal oxides has emerged as an intriguing approach for large-scale hydrogen production. Perovskites have been proposed as alternatives to state-of-the-art fluorite CeO 2–δ because of their potential for lowering reduction temperature while maintaining high fuel productivity. Guided by computational insights, we explore the thermodynamic properties and water splitting efficacy of the cubic perovskite SrTi 0.5 Mn 0.5 O 3–δ (STM55). Thermogravimetric analysis is performed under controlled oxygen partial pressures (pO 2 ) and temperatures up to 1500 °C, from which both the enthalpy and entropy of reduction as a function of oxygen nonstoichiometry are determined. STM55 provides an attractive combination of moderate enthalpy, 200–250 kJ (mol-O) –1 , and high entropy, with unusual δ dependence. Using a water splitting cycle in which the material is thermally reduced at 1350 °C (pO 2 , ~10 –5 atm) and subsequently exposed to steam at 1100 °C (steam partial pressure of pH 2 O = 0.4 atm), we demonstrate a hydrogen yield of 7.4 mL g –1 . Through both half-cycles, the material remains largely in quasi-equilibrium with the gas phase, as reflected in the agreement of the measured data with predicted profiles based on the thermodynamic data. This behavior indicates rapid surface and bulk diffusion kinetics. Cyclic operation showed the material to be free of degradation and always resulted in a 2:1 yield of H 2 /O 2 . Altogether, STM55 provides outstanding performance characteristics for thermochemical hydrogen production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Uncertainty-Driven Rapid Thermodynamic Assessment of Nb-Ta-Zr System and Effects of C impurities (L25GF9298S): Annual Progress Report

An integrated computational materials engineering (ICME) method is in development for rapid thermodynamic experimental investigation and high-fidelity computational modeling of refractory multi-principal element alloys (RMPEAs). These ultra-high-temperature (UHT) alloys are of interest for structural applications in extreme environments, but deficiency of reliable data, especially melting temperatures, impedes the prediction of alloys with favorable properties. The method leverages UHT capabilities and computational expertise of LLNL’s Materials Science Division and the McCormack Lab’s UHT conical nozzle levitation (CNL) system to iteratively map the Nb-Ta-Zr phase space, with focus on the liquidus surface, through targeted experiments selected by quantifying uncertainty in the thermodynamic model fitting parameters. This method will reduce the time to map uncharted RMPEA phase space and thereby accelerate discovery and development of advanced materials for applications in extreme environments.

36 MATERIALS SCIENCE↗

Effect of iodides on thermal behavior and phase partitioning in LiCl-KCl

Liquid-fueled molten salt reactors (MSRs) are designed to operate with fissile materials and, ultimately, fission products dissolved in the primary molten salt coolant. Understanding the speciation and transport of iodine—a high-yield fission product—is essential because this element’s accidental release poses significant environmental concerns due to its capacity to be readily absorbed by the human thyroid gland. Here, we report the impact of iodide species (LiI and KI) on phase transitions, phase distribution, and phase stability in LiCl–KCl-eutectic salt mixtures. The study employed a combination of computational and experimental techniques, including thermodynamic FactSage calculations, differential scanning calorimetry, and high-temperature X-ray diffraction. The results indicate that the presence of iodide (10–25 wt%) significantly alters the melting behavior of the LiCl–KCleutectic system. Adding 10 wt% LiI has a more-pronounced effect than 10 wt% KI, as LiI converts to KI, leading to formation of LiCl, thereby, altering the LiCl-KCl ratio which significantly affects the melting temperature of the mixture. Furthermore, the evolution of crystalline structure, solid-fraction composition, and the dynamics of mixed-halide solid–liquid partitioning as a function of temperature indicate the potential for selective iodide separation from chloride-salt mixtures via solid–liquid separation techniques. Overall, the presented findings provide valuable insights that are beneficial for the design and operation of MSRs, as well as for the safe handling and effective processing of used nuclear fuel using advanced pyrochemical techniques.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Integrated modeling methodology for ash agglomeration in poly-disperse fluidized beds using particle population framework

This article discusses a unique agglomeration modeling methodology developed based on binary collisions to combine the effects of heterogeneity in ash chemistry and granular physics. A simple population balance is defined to find changes in the particle size distribution (PSD) of a fluidized bed. Thermodynamic equilibrium calculations and a computational fluid dynamics (CFD) code, are used to obtain hydrodynamic parameters. A method to calculate and use a distribution of collision frequencies in fluidized beds, corresponding to a poly-disperse particle size distribution, was developed in order to incorporate the particle level heterogeneities. A distribution of collision frequencies obtained for poly-dispersed particles showed a three orders of magnitude higher collision frequency amongst the smaller particles than the coarser ones, at the initiation of agglomeration. Ash agglomeration occurred when the slag amount (binder) was less than 10 wt% at temperatures <850 °C for Pittsburgh No. 8 coal.

42 ENGINEERING↗

Solvent Isotherms and Structural Transitions in Nanoparticle Superlattice Assembly

Here, we introduce a Molecular Theory for Compressible Fluids (MOLT-CF) that enables us to compute free energies and other thermodynamic functions for nanoparticle superlattices with any solvent content, including the dry limit. Quantitative agreement is observed between MOLT-CF and united-atom molecular dynamics simulations performed to assess the reliability and precision of the theory. Among other predictions, MOLT-CF shows that the amount of solvent within the superlattice decreases approximately linearly with its vapor pressure and that in the late stages of drying, solvent-filled voids form at lattice interstitials. Applied to single-component superlattices, MOLT-CF predicts fcc-to-bcc Bain transitions for decreasing vapor pressure and for increasing ligand length, both in agreement with experimental results. We explore the stability of other single-component phases and show that the C14 Frank-Kasper phase, which has been reported in experiments, is not a global free-energy minimum. Implications for precise assembly and prediction of multicomponent nanoparticle systems are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Alkane Dehydrogenation and H/D Exchange by a Cationic Pincer-Ir(III) Hydride: Cooperative C–H Addition and β -H Elimination Modes Induce Anomalous Selectivity

We report that the cationic iridium complex ( iPr PCP)IrH + catalyzes the transfer-dehydrogenation of alkanes to give alkenes and hydrogen isotope exchange (HIE) of alkanes and arenes. Contrary to established selectivity trends found for C–H activation by transition metal complexes, strained cycloalkanes, including cyclopentane, cycloheptane, and cyclooctane, undergo C–H addition much more readily than n-alkanes, which in turn are much more reactive than cyclohexane. Aromatic C–H bonds also undergo H/D exchange much less rapidly than those of the strained cycloalkanes, but much more favorably than cyclohexane. The order of reactivity toward dehydrogenation correlates qualitatively with the reaction thermodynamics, but the magnitude is much greater than can be explained by thermodynamics. Accordingly, the cycloalkenes corresponding to the strained cycloalkanes undergo hydrogenation much more readily than cyclohexene, despite the less favorable thermodynamics of such hydrogenations. Here, computational (DFT) studies allow rationalization of the origin of reactivity and the unusual selectivity. Specifically, the initial C–H addition is strongly assisted by β-agostic interactions, which are particularly favorable for the strained cycloalkanes. Subsequent to α-C–H addition, the H atom of the β-agostic C–H bond is transferred directly to the hydride ligand of ( iPr PCP)IrH + to give a dihydrogen ligand.

02 PETROLEUM↗

Observation of laser ablation of silicon as a function of pulse length at constant fluence via time-resolved x-ray spectroscopy

We investigate the ablation of silicon as a function of laser pulse length at a constant fluence using time-resolved x-ray spectroscopy data obtained from OMEGA EP experiments at the University of Rochester's Laboratory for Laser Energetics. Our targets consisted of three-layer planar structures composed of Si (50 μm), Cu (25 μm), and SiO 2 (500 μm) layers. The Si layer was irradiated by a 351-nm laser with varying pulse widths of 250 ps, 500 ps, 1 ns, and 10 ns while maintaining a constant fluence of ~27.9 kJ/cm 2 . Electron temperatures and densities of the ablated plasma were determined by analyzing the time-resolved x-ray spectroscopy data through a comparison of experimental measurements with synthetic results obtained from Si atomic calculations in a steady state and non-local thermodynamic equilibrium. These calculations were computed using PrismSPECT. Additionally, radiation-hydrodynamics simulations with FLASH are used to generate simulated plasma-density and plasma-temperature profiles, which are then compared with the experimental measurements. Our analyses reveal that increasing the laser pulse length at a constant fluence results in a decrease in electron temperatures and densities. Furthermore, the longer pulses with lower intensities lead to deeper ablation regions before reaching the peak ablation but lower ionization balances in the silicon layer. Furthermore, these findings emphasize the critical role of laser pulse length in plasma ablation and shock generation for laser-impulse studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computational investigation of water glasses using machine-learning potentials

The molecular origins of water’s anomalous properties have long been a subject of scientific inquiry. The liquid–liquid phase transition hypothesis, which posits the existence of distinct low-density and high-density liquid states separated by a first-order phase transition terminating at a critical point, has gained increasing experimental and computational support and offers a thermodynamically consistent framework for many of water’s anomalies. However, experimental challenges in avoiding crystallization near the postulated liquid–liquid critical point have focused attention to water’s canonical glassy states: low-density and high-density amorphous ice. Here, we use two Deep Potential machine-learning models, trained on the Strongly Constrained and Appropriately Normed density functional and the highly accurate Many-Body Polarizable potential, to conduct an investigation of water’s glassy phenomenology based on quantum mechanical calculations. Despite not being explicitly trained on amorphous ices, both models accurately capture the structure and transformation of the water glasses, including their interconversion along different thermodynamic paths. Isobaric quenching of liquid water at various pressures generates a continuum of intermediate amorphous ices and density fluctuations increase near the liquid–liquid critical pressure. The glass transition temperatures of the amorphous ices produced at different pressures exhibit two distinct branches, corresponding to low-density and high-density amorphous ice behaviors, consistent with experiment and the liquid–liquid transition hypothesis. Extrapolating transformation pressures from isothermal compressions to experimental compression rates brings our simulations into excellent agreement with data. Our findings demonstrate that machine-learning potentials trained on equilibrium phases can effectively model nonequilibrium glassy behavior and pave the way for studying long-timescale, out-of-equilibrium processes with quantum mechanical accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Uranium Bed Design Parameters for Tritium Plants Supporting Fusion Reactors

Uranium hydride is commonly used to store hydrogen or its isotopes in a solid state. The Self-Assaying Tritium Accountancy and Containment Unit for ITER (STACI) is a 5.2-kg bed of depleted uranium (dU) capable of holding up to 33 mol of hydrogen or its isotopes. Here we summarize the data analysis of past experimental campaigns with STACI, with the aim of describing the kinetics and thermodynamics of the hydriding process. Computed tomography imaging was performed on STACI both before and after its experimental campaign, and a high degree of swelling was observed in the dU. Literature on studies in regard to the swelling of large (multikilogram) quantities of uranium hydride for storage applications was not identified during this study. Data from the experimental campaign, as well as data on the formation reaction, are presented. The authors hope to create an analytical model of STACI based on these data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Origin of mean-field behavior in an elastic Ising model

Simple elastic models of spin-crossover compounds are known empirically to exhibit classical critical behavior. Here, we demonstrate how the long-range interactions responsible for this behavior arise naturally upon integrating out mechanical fluctuations of such a model. A mean-field theory applied to the resulting effective Hamiltonian quantitatively accounts for both thermodynamics and kinetics observed in computer simulations, including a barrier to magnetization reversal that grows extensively with system size. For nanocrystals, which break translational symmetry, a straightforward extension of mean-field theory yields similarly accurate results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multimodal Ligand Non-Innocence for Reduction Catalysis (Final Technical Report)

The goal of this project was to develop multimodal ligand non-innocence in late transition metal complexes of aminophenolate ligands and to apply this ligand-centered reactivity to small-molecule reduction reactions relevant to energy-conversion chemistry. Specific aims were, 1) Prepare new coordination complexes of aminophenolate ligands and characterize them using a battery of structural, spectroscopic, and computational methods. 2) Benchmark the thermodynamics and kinetics of e - , H + , H • , and H - transfer reactivity to establish the baseline reactivity of new metal complexes containing aminophenolate ligands. 3) Apply new aminophenolate complexes to catalytic reduction reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗