Search NASA⌕ Search

SEARCH · Search NASA

Results for “predictive 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Three-Dimensional, Thermodynamically and Variationally Consistent, Fully Coupled, Electro-Chemo-Thermo-Mechanical Model of Solid-State Batteries

Here we present a theoretical framework for high fidelity modeling of three-dimensional electro-chemo-thermo-mechanical coupled systems such as solid-state batteries. The framework provides a useful means of utilizing molecular level information to predict the overall performance of the system at the continuum level. The central statement of the theoretical framework is a thermodynamically consistent inf-sup problem. The total potential of the coupled system consists of the Helmholtz free energy, the electrical kinetic potential, the chemical potential, the Fourier potential, the chemical reaction potential, and the external power expenditure. All governing equations derive from variational principles. We also show that the resulting initial-boundary-value problem (IBVP) is consistent with existing, widely adopted equations for modeling batteries. Moreover, the IBVP can be naturally discretized using the (variational) finite element method to simulate practical solid-state batteries with arbitrary geometries. Several examples illustrate multiphysics coupling between pairs of the physical field equations, culminating in a comprehensive three-dimensional model of charge-discharge in a solid state battery.

25 ENERGY STORAGE↗
DOE OSTI2021

Evolution of the Thermodynamic Properties of Clusters of Galaxies out to Redshift of 1.8

The thermodynamic properties of the hot plasma in galaxy clusters retains information on the processes leading to the formation and evolution of the gas in their deep, dark matter potential wells. These processes are dictated not only by gravity but also by gas physics, e.g. AGN feedback and turbulence. Here, we study the thermodynamic properties, e.g. density, temperature, pressure, and entropy, of the most massive and the most distant (z > 1.2) SPT-selected clusters, and compare them with those of the nearby clusters (z < 0.1) to constrain their evolution as a function of time and radius. We find that thermodynamic properties in the outskirts of high redshift clusters are remarkably similar to the low redshift clusters, and their evolution follows the prediction of the self-similar model. Their intrinsic scatter is larger, indicating that the physical properties that lead to the formation and virialization of cluster outskirts show evolving variance. On the other hand, thermodynamic properties in the cluster cores deviates significantly from self-similarity indicating that the processes that regulate the core are already in place in these very high redshift clusters. This result is supported by the unevolving physical scatter of all thermodynamic quantities in cluster cores.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Evolution of the Thermodynamic Properties of Clusters of Galaxies out to Redshift of 1.8

The thermodynamic properties of the hot plasma in galaxy clusters retains information on the processes leading to the formation and evolution of the gas in their deep, dark matter potential wells. These processes are dictated not only by gravity but also by gas physics, e.g. AGN feedback and turbulence. In this work, we study the thermodynamic properties, e.g. density, temperature, pressure, and entropy, of the most massive and the most distant ($z > 1.2$) SPT-selected clusters, and compare them with those of the nearby clusters ($z<0.1$) to constrain their evolution as a function of time and radius. We find that thermodynamic properties in the outskirts of high redshift clusters are remarkably similar to the low redshift clusters, and their evolution follows the prediction of the self-similar model. Their intrinsic scatter is larger, indicating that the physical properties that lead to the formation and virialization of cluster outskirts show evolving variance. On the other hand, thermodynamic properties in the cluster cores deviates significantly from self-similarity indicating that the processes that regulate the core are already in place in these very high redshift clusters. This result is supported by the unevolving physical scatter of all thermodynamic quantities in cluster cores.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Data‐Driven Insights into Rare Earth Mineralization: Machine Learning Applications Using Functional Material Synthesis Data

Understanding rare‐earth element (REE) mineralization mechanisms is essential for developing efficient separation strategies. Although the geochemical pathways that generate REE deposits are qualitatively known, quantitative links between specific conditions and mineralization outcomes remain limited. Herein, the repurpose laboratory REE hydrothermal synthesis data—originally collected for functional‐materials fabrication—as a surrogate for studying mineralization with data‐driven methods. The compiled 1,200+ hydrothermal reaction records and trained three machine‐learning models—K‐nearest neighbors (KNN), random forest (RF), and extreme gradient boosting (XGB)—to predict product elements and phases from precursors, additives, reaction conditions, and engineered features. Validation shows XGB achieves the highest accuracy. Feature importance indicates thermodynamic properties of cations and anions dominate model decisions. Correlations reveal positive relationships among precursor concentration, reaction time, pH, and temperature, consistent with classical crystallization behavior. XGB‐based regressors are built to predict crystallization temperature and pH from precursor/product attributes. Performance is strongest when similar training examples exist, while accuracy declines for underrepresented reactions, notably REE carbonates and heavy‐REE systems. Overall, the study shows that functional‐materials datasets can illuminate REE mineralization and provide priors for exploration and processing. Expanding datasets with less‐studied chemistries and conditions will improve generality and support deposit discovery and more efficient REE recovery.

feature importance analysis↗
DOE OSTI2025

Charting the chemical space of Zintl phases with graph neural networks and bonding insights

A large number of Zintl phases have been discovered by solid-state chemists driven by empirical knowledge, chemical intuition and in some cases, through serendipitous accidents. These discoveries have only scratched the surface, given the vast compositional and structural diversity that Zintl phases can accommodate. The large chemical space of Zintl phases, as well as intermetallic compounds in general, remain under-explored. Here, we use graph neural networks and the upper bound energy minimization approach to efficiently scan a large chemical space of >90 000 hypothetical Zintl phases and accurately discover 1810 new thermodynamically stable phases with 90% precision, as validated with first-principles calculations. We show that our approach is more than 2× more accurate in predicting DFT stability than M3GNet (40% precision) on the same dataset. Using a random forest model and SHAP analysis, we demonstrate the critical role of ionic bonding in the thermodynamic stability of Zintl phases. Our results not only expand the known chemical landscape of Zintl phases but also highlight the efficacy of machine learning frameworks combined with domain knowledge in uncovering chemically meaningful insights across complex intermetallics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A thermodynamically consistent finite strain phase field approach to ductile fracture considering multi-axial stress states

Phase field models for ductile fracture have gained significant attention in the last two decades due to their ability in implicitly tracking the nucleation and propagation of cracks. However, most crack phase field formulations for elastoplastic solids focus only on the effects of plastic deformation, and do not consider the different multi-axial stress states that may arise in practical designs. In this work, a thermodynamically consistent phase field approach coupled with finite strain plasticity, considering multi-axial stress states is presented. In order to account for the coupling between plasticity and stress states, the Stress-Weighted Ductile Fracture Model (SWDFM) is utilized. The SWDFM represents a criterion for predicting ductile crack initiation under both monotonic and cyclic loadings based on histories of an internal plastic variable, stress triaxiality, and the Lode angle parameter. The excellent performance of the SWDFM for predicting ductile crack initiation motivates for its incorporation into a phase field approach for predicting both crack initiation and propagation through degradation of the fracture toughness. Moreover, based on the second law of thermodynamics, exact requirements are imposed on the rate at which the fracture toughness can evolve. Here, a novel function for degrading the plastic yield surface during the evolution of damage is introduced. This function, in line with experimental observations, leads to an accumulation of plastic deformation in damaged regions of a solid, and avoids numerical instabilities arising from concentrations of large plastic deformations in severely damaged regions. For validating the proposed model, results of computational simulations are compared to data from selected tests considering different multi-axial stress states. Comparisons of the numerical results with data from laboratory experiments demonstrate the capabilities of the proposed framework.

42 ENGINEERING↗
DOE OSTI2024

Efficient data-driven models for prediction and optimization of geothermal power plant operations

Increasing the capacity of geothermal energy as a renewable resource calls for development and deployment of efficient control and optimization technologies for geothermal power plants. A data-driven prediction and optimization model is presented as a cost-effective and efficient alternative to physics-based approach. The model predicts power output and operational cost by propagating the influence of control and disturbance variables within an artificial neural network (ANN). Numerical experiments with simulated and field data from a real geothermal power plant are first used to demonstrate the prediction performance of the ANN model. The model is then adopted to maximize the net predicted power production by automatically adjusting the working fluid circulation rate. The optimization performance of the model in evaluated using a thermodynamic flowsheet simulation model. The workflow is applied to model and control the effect of ambient temperature on an air-cooled binary cycle power plant, which is complex and costly to perform using a physics-based predictive model. As a result, the performance of the method is demonstrated by applying it to both simulated and field datasets from a binary cycle geothermal power plant.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Strain phase equilibria and phase‐field method of ferroelectric polydomain: A case study of monoclinic K x Na 1 − x NbO 3 thin films

Abstract Knowledge of the thermodynamic equilibria and domain structures of ferroelectrics is critical to establishing their structure–property relationships that underpin their applications from piezoelectric devices to nonlinear optics. Here, we establish the strain condition for strain phase separation and polydomain formation and analytically predict the corresponding domain volume fractions and wall orientations of, relatively low symmetry and theoretically more challenging, monoclinic ferroelectric thin films by integrating thermodynamics of ferroelectrics, strain phase equilibria theory, microelasticity, and phase‐field method. Using monoclinic K x Na 1 − x NbO 3 (0.5 < x < 1.0) thin films as a model system, we establish the polydomain strain–strain phase diagrams, from which we identify two types of monoclinic polydomain structures. The analytically predicted strain conditions of formation, domain volume fractions, and domain wall orientations for the two polydomain structures are consistent with phase‐field simulations and in good agreement with experimental results in the literature. The present study demonstrates a general, powerful analytical theoretical framework to predict the strain phase equilibria and domain wall orientations of polydomain structures applicable to both high‐ and low‐symmetry ferroelectrics and provide fundamental insights into the equilibrium domain structures of ferroelectric K x Na 1 − x NbO 3 thin films that are of technology relevance for lead‐free dielectric and piezoelectric applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Co-Based superalloy morphology evolution: A phase field study based on experimental thermodynamic and kinetic data

Cobalt-based superalloys with two phase $\gamma$/$\gamma '$ microstructures offer great promise as candidates for next-generation high-temperature alloys for applications, such as turbine blades. It is essential to understand the thermodynamic and kinetic factors that influence the microstructural evolution of these alloys in order to optimize the alloy compositions and processing steps with a goal to improve their coarsening, creep and rafting behavior. We are using a continuum phase field approach to study the diffusion process and to predict the equilibrium shapes of Co-Al-W $\gamma '$ precipitates. In order to obtain quantitatively predictive capabilities, we extract chemical free energies for the $\gamma$/$\gamma '$ phases based on CALculation of PHAse Diagrams (CALPHAD) thermodynamic data and diffusion mobilities for Co alloys based on CALPHAD kinetic data. We also use experimental or first-principles data for other quantities, such as misfit strain and interface information, for the parameterization of our model. A particular focus of our study is to understand how different energy balances, misfit strain and kinetics affect the coarsening and rafting behavior of $\gamma '$ precipitates, and the sensitivity of the final precipitate shape to materials parameters. Here, we find that the equilibrium shape of the precipitate results from a delicate competition between chemical, interfacial, and elastic energies, and it is very sensitive to changes in model parameters. Here, we examine how modeling input parameters affect the equilibrium shape of precipitates and relate these parameters to experimentally available values.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Resolving local ordering and structure in Mn x Ge 1- x Te alloys through thermodynamic ensembles of pair distribution functions

Characterizing local bonding environments in complex materials is essential for understanding and optimizing their properties. Equally as important is the ability to predict local motifs as a function of synthesis conditions, enhancing chemists’ ability to design properties into materials. In this study, we present an approach to leverage statistical mechanics to generate temperature- and energy-informed ensemble averaged pair distribution functions (PDFs). This method, which we have named Thermodynamic Ensemble Averages of PDFs for Ordering and Transformations (TEAPOT), utilizes density functional theory (DFT) to relax supercells while incorporating energetic penalties for local order, enabling accurate and computationally efficient analysis of local structure. We apply this method to the neutron PDF measurements of the pseudobinary MnTe–GeTe (MGT) alloy, demonstrating its capability to resolve complex local distortions and chemical ordering. Our results reveal detailed insights into phase transformations and local distortions driven by Mn substitution. For compositions that globally present as rock salt, our analysis reveals that Ge coordination geometry is heavily impacted by synthesis temperature. We propose that high temperature synthesis conditions promote a lowered Ge polyhedra distortion, promoting high charge carrier mobility due to the alignment of local and global structure. Incorporating statistical mechanics and computation into experimental analysis thus guides synthesis of tailored local structure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Prediction of the Al-rich part of the Al-Cu phase diagram using cluster expansion and statistical mechanics

In this work, the thermodynamic properties of α-Al and other phases (GP zones, θ&#x27;&#x27;, θ&#x27; and θ) in the Al-rich part of the Al-Cu system have been obtained by means of the cluster expansion formalism in combination with statistical mechanics. This information was used to build the Al-rich part of the Al-Cu phase-diagram taking into account vibrational entropic contributions for θ&#x27;, as those of the other phases were negligible. The simulation predictions of the phase boundaries between α-Al and either θ&#x27;&#x27;, θ&#x27; or θ phases as a function of temperature are in good agreement with experimental data and extend the phase boundaries to a wider temperature range. The DFT calculations reveal the presence of a number of metastable Guinier-Preston-zone type configurations that may coexist with α-Al and θ&#x27;&#x27; at low temperatures. They also demonstrate that θ&#x27; is the stable phase below 550K but it is replaced by θ above this temperature due to the vibrational entropic contribution to the Gibbs energy of θ&#x27;. This work shows how the combination of cluster expansion and statistical mechanics can be used to expand our knowledge of the phase diagram of metallic alloys and to provide Gibbs free energies of different phases that can be used as input in mesoscale simulations of precipitation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Capturing Surface Coverage Effects in Heterogeneous Catalysis

Adsorbate–adsorbate lateral interactions at relevant surface coverages have a significant effect on chemical kinetics, thereby influencing the activity of a heterogeneous catalyst. Coverage-dependent kinetic and thermodynamic parameters therefore must be included in studies of such complex systems to properly predict the turnover frequencies and kinetic trends. Thus, it becomes extremely important to accurately capture the strength of lateral interactions between neighboring species under realistic reaction conditions. In this Perspective, we discuss the various existing computational and experimental methods for determining adspecies coverage and configurational effects. The choice of the tools and methods employed in such studies depends on factors such as time, length scales, computational cost, the presence of solvents, and reaction conditions. The applications of each method and the respective challenges are also discussed here. As a result, we discuss the recent developments and future of the state-of-the-art for inclusion of surface coverage and configuration into a holistic picture for accurate predictions of catalytic behavior.

09 BIOMASS FUELS↗
DOE OSTI2021

Thermodynamic stability versus kinetic accessibility: Pareto fronts for programmable self-assembly

A challenge in designing self-assembling building blocks is to ensure the target state is both thermodynamically stable and kinetically accessible. These two objectives are known to be typically in competition, but it is not known how to simultaneously optimize them. We consider this problem through the lens of multi-objective optimization theory: we develop a genetic algorithm to compute the Pareto fronts characterizing the tradeoff between equilibrium probability and folding rate, for a model system of small polymers of colloids with tunable short-ranged interaction energies. We use a coarse-grained model for the particles' dynamics that allows us to efficiently search over parameters, for systems small enough to be enumerated. For most target states there is a tradeoff when the number of types of particles is small, with medium-weak bonds favouring fast folding, and strong bonds favouring high equilibrium probability. Additionally, the tradeoff disappears when the number of particle types reaches a value m *, that is usually much less than the total number of particles. This general approach of computing Pareto fronts allows one to identify the minimum number of design parameters to avoid a thermodynamic–kinetic tradeoff. However, we argue, by contrasting our coarse-grained model's predictions with those of Brownian dynamics simulations, that particles with short-ranged isotropic interactions should generically have a tradeoff, and avoiding it in larger systems will require orientation-dependent interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Improved reference system for the corrected rigid spheres equation of state model

The Corrected Rigid Spheres (CRIS) EOS model, developed from fluid perturbation theory using a hard sphere reference system, has been successfully used to calculate the EOS of many materials, including gases and metals. The radial distribution function (RDF) plays a pivotal role in choosing the sphere diameter, through a variational principle, as well as the thermodynamic response. Despite its success, the CRIS model has some shortcomings in that it predicts too large a temperature for liquid-vapor critical points, can break down at large compression, and is computationally expensive. In this paper, we first demonstrate that an improved analytic representation of the hard sphere RDF does not alleviate these issues. Relaxing the strict adherence of the RDF to hard spheres allows an accurate fit to the isotherms and vapor dome of the Lennard-Jones fluid using an arbitrary reference system. Second order correction are eliminated, limiting the breakdown at large compression and significantly reducing the computation cost. The transferrability of the new model to real systems is demonstrated on argon, with an improved vapor dome compared to the original CRIS model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Computational and Experimental Investigation of Thermal-Mechanical- Chemical Mechanisms of High-burnup Spent Nuclear Fuel (SNF) Processes at Elevated Temperatures and Degradation Behavior in Geologic Repositories

The overarching goal of the combined computational and experimental R&D activities proposed in this project is to enhance understanding of the mechanisms and thermal-mechanical-chemical (TMC) parameters controlling the instant release fraction (IRF) and matrix dissolution of high-burnup (HB; burnup) spent nuclear fuels (SNFs) and the subsequent formation, stability, and phase transformations of SNF alteration products under long-term storage and geological disposal conditions. Uranium dioxide may undergo oxidative corrosion/alteration, and the IRF may be increased for HB SNF, both of which may affect environmental systems associated with SNF long-term storage and disposal. The oxidative matrix dissolution may form various complex uranyl-based phases, including a rich variety of oxides, silicates, carbonates and other secondary minerals in varied geological environments (e.g., studtite, metastudtite, amorphous uranyl peroxide, uranium trioxide, triuranium octoxide, schoepite, dehydrated schoepite, metaschoepite, becquerelite, soddyite, rutherfordine,...). These uranyl phases generally have higher mobility UO 2 +2 species than less soluble U 4+ phases. However, limited information on the thermodynamic properties and formation kinetics of these uranyl-bearing phases is available to predict explicitly paragenesis under the conditions relevant to long-term storage or disposal. The proposed project draws on complementary expertise and research backgrounds from the team members: (i) to apply a combined ab initio modeling (UNLV/UTEP and SNL) and experimental (UNLV) strategy investigating the high-temperature TMC mechanisms of alteration of SNF under α-radiolysis conditions; (ii) to investigate the mechanistic of phase transformations in UNF degradation products under various conditions expected in long-term storage systems (e.g. (UO 2 )O 2 (H 2 O) 4 → (UO 2 )O 2 (H 2 O) 2 → U 2 O 7 → UO 3 → U 3 O 8 ); (iii) to determine high-accuracy TMC parameters for complex uranyl-based phases formed in storage or geological disposal environments (e.g. UO 3 (H 2 O) 2 , Ca[(UO 2 ) 6 O 4 (OH) 8 ] 8 H 2 O, (UO 2 ) 2 (SiO 4 ) 3 2H 2 O,…). The unforeseen COVID-19 pandemic led to the laboratory/campus closure since March 2020, that resulted in a significant delay in reaching milestones in a satisfactory manner, due to (i) the statewide recommendation from stop-working to later limited work in the lab and work-from-home (WFH), (ii) no in-person interactions, and (iii) a hiring freeze at UNLV. Therefore, a no cost extension (10/01/2021- 9/30/2022) was requested to help make up the time we lost during the global pandemic in 2020-2021, leading to paradigm shifts in the focus of the project in the following three main tasks: Task 1 (Computational), Task 2 (Experimental), and Task 3 (Final report, due on 12/29/2022).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Time-dependent water vapor desorption isotherm model of hardened cement paste

Highlights: • Short-term and long-term desorption isotherms of hardened cement paste are reviewed. • Water content at RH depends on both thermodynamic and microstructural equilibria. • An extended GAB equation is proposed to fit the cavitation drop during desorption. • A dynamic time-dependent sorption isotherm model is shown and used for simulations. • Numerical results capture anomalous drying in the same RH as seen experimentally. This study aims to enhance moisture transport modeling by elucidating the so-called anomalous water diffusion in cementitious materials. Water desorption isotherms are measured from samples at various drying stages using techniques with different durations to quantify the impact of the drying duration on the water sorption ability. A single water sorption isotherm does not solely give a relationship between the water content and relative humidity but also represents the state of microstructure. The continuous evolution of a desorption isotherm due to drying-induced microstructural rearrangement is demonstrated. For numerical modeling, the microstructural alteration can be explicitly considered through a dynamic desorption isotherm model, which governs the local thermodynamic equilibrium at the capillary meniscus. This approach is implemented into a multiphase transport code, whose ability to predict drying is validated using literature data. Finally, the effect of prolonged drying on the colloidal nature of the calcium-silicate-hydrate gel is discussed.

36 MATERIALS SCIENCE↗