Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase stability modeling”

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
DOE OSTI2022

Experimental and theoretical examination of shock-compressed copper through the fcc to bcc to melt phase transitions

Recent studies show a face-centered cubic (fcc) to body-centered cubic (bcc) transformation along the shock Hugoniot for several metals (i.e., Cu, Au, and Ag). Here, we combine laser-shock compression of Cu foils on nanosecond timescales with in situ x-ray diffraction (XRD) to examine the microstructural changes with stress. We study the fcc phase and the phase transition from fcc to bcc (pressures greater than 180 GPa). Textural analysis of the azimuthal intensities from the XRD images is consistent with transformation into the bcc phase through the Pitsch-distortion mechanism. We use embedded atom model molecular dynamics simulations to determine the stability of the bcc phase in pressure–temperature space. Finally, our results indicate that the bcc phase is stabilized only at high temperatures and remains stable at pressures greater than 500 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2024

FY-24 Progress on Computational Modeling of the Water Based NSTF

This report summarizes the computational modeling of the Natural Convection Shutdown Heat Removal Test Facility (NSTF) completed in FY24. This year’s modeling campaign focuses on the continual testing of the RELAP5 model against experimental data. Several fault cases were simulated with the RELAP5 model in FY23 to study the accuracy of the model under complex flow conditions. Similarly, in FY24, a case is simulated where the chimney outlet is throttled with a valve while the flow is undergoing two-phase oscillations. The purpose of this case is to study the effects of increasing pressure drop in the two-phase region on the system behavior. Even though the model is able to capture the experimental data qualitatively, it overpredicts the pressure drops experienced by the flow in the two-phase regime. Nevertheless, the model and the experiment show that oscillations are stabilized with increased pressure drop in the two-phase region.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Thermodynamics of monoclinic and tetragonal hafnium dioxide (HfO 2 ) at ambient pressure

Accurate and precise thermodynamic models of oxide compounds and their phases are important for calculating the phase stability of oxide materials. We develop and use a coupled quantum mechanical and molecular dynamics approach to create thermodynamic models of hafnia (HfO2) polymorphs from 0 K to 3000 K at ambient pressure. The approach is based on the quasi-harmonic approximation below the Debye temperature and on ab-initio molecular dynamics calculations above the Debye temperature to predict constant pressure heat capacities (C-p). A Bayesian model provides interpolated values between these regimes. As a case study, we develop thermodynamic models of monoclinic and tetragonal HfO2 polymorphs. The predicted heat capacities are in excellent agreement with experiment, and the predicted temperature of the monoclinic to tetragonal phase transition (2173 K) is in good agreement with the experimental value (2078K). These results provide a comprehensive and accurate thermodynamic model of the monoclinic and tetragonal phases of hafnia on a broad range of temperatures and can serve as input to CALPHAD assessment of multi-component hafnia-based phase diagrams.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Phase-Field Model for In-Space Manufacturing of Binary Alloys

The integrity of the final printed components is mostly dictated by the adhesion between the particles and phases that form upon solidification, which is a major problem in printing metallic parts using available In-Space Manufacturing (ISM) technologies based on the Fused Deposition Modeling (FDM) methodology. Understanding the melting/solidification process helps increase particle adherence and allows to produce components with greater mechanical integrity. We developed a phase-field model of solidification for binary alloys. The phase-field approach is unique in capturing the microstructure with computationally tractable costs. The developed phase-field model of solidification of binary alloys satisfies the stability conditions at all temperatures. The suggested model is tuned for Ni-Cu alloy feedstocks. We derived the Ginzburg-Landau equations governing the phase transformation kinetics and solved them analytically for the dilute solution. We calculated the concentration profile as a function of interface velocity for a one-dimensional steady-state diffuse interface neglecting elasticity and obtained the partition coefficient, k, as a function of interface velocity. Numerical simulations for the diluted solution are used to study the interface velocity as a function of undercooling for the classic sharp interface model, partitionless solidification, and thin interface.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Machine learning assisted derivation of minimal low-energy models for metallic magnets

Abstract We consider the problem of extracting a low-energy spin Hamiltonian from a triangular Kondo Lattice Model (KLM). The non-analytic dependence of the effective spin-spin interactions on the Kondo exchange excludes the use of perturbation theory beyond the second order. We then introduce a Machine Learning (ML) assisted protocol to extract effective two- and four-spin interactions. The resulting spin model reproduces the phase diagram of the original KLM as a function of magnetic field and single-ion anisotropy and reveals the effective four-spin interactions that stabilize the field-induced skyrmion crystal phase. Moreover, this model enables the computation of static and dynamical properties with a much lower numerical cost relative to the original KLM. A comparison of the dynamical spin structure factor in the fully polarized phase computed with both models reveals a good agreement for the magnon dispersion even though this information was not included in the training data set.

Sharma, Vikram (ORCID:0000000189105519)↗
DOE OSTI2022

Entanglement phase transitions in random stabilizer tensor networks

Here we explore a class of random tensor network models with “stabilizer” local tensors which we name random stabilizer tensor networks (RSTNs). For RSTNs defined on a two-dimensional square lattice, we perform extensive numerical studies of entanglement phase transitions between volume-law and area-law entangled phases of the one-dimensional boundary states. These transitions occur when either (a) the bond dimension D of the constituent tensors is varied or (b) the tensor network is subject to random breaking of bulk bonds, implemented by forced measurements. In the absence of broken bonds, we find that the RSTN supports a volume-law entangled boundary state with bond dimension D ≥ 3 where D is a prime number, and an area-law entangled boundary state for D = 2. Upon breaking bonds at random in the bulk with probability p, there exists a critical measurement rate p c for each D ≥ 3 above which the boundary state becomes area-law entangled. To explore the conformal invariance at these entanglement transitions for different prime D, we consider tensor networks on a finite rectangular geometry with a variety of boundary conditions, and extract universal operator scaling dimensions via extensive numerical calculations of the entanglement entropy, mutual information, and mutual negativity at their respective critical points. Our results at large D approach known universal data of percolation conformal field theory, while showing clear discrepancies at smaller D, suggesting a distinct entanglement transition universality class for each prime D. We further study universal entanglement properties in the volume-law phase and demonstrate quantitative agreement with the recently proposed description in terms of a directed polymer in a random environment.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Design of lightweight BCC multi-principal element alloys with enhanced hydrogen storage using a machine learning-driven genetic algorithm

Body-centered cubic (BCC) based multi-principal element alloy (MPEA) hydrides have demonstrated significant potential for compact and efficient hydrogen storage. In this work, we first leverage machine learning (ML) models to predict the hydrogen affinity, storage capacity and phase stability of BCC MPEAs, creating a unique hydrogen-to-metal (H/M) predictor for materials with unprecedented performance. We developed a metaheuristic optimizer high-throughput framework by interfacing ML models with a genetic algorithm for the accelerated search of {Mg, Al, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Nb, Mo} based lightweight BCC MPEAs with improved hydrogen storage characteristics. We report five new MPEAs with a predicted gravimetric hydrogen storage capacity of around 3.5 wt% or more, including Cr 0.09 Mg 0.73 Ti 0.18 (4.25 wt% H) and Cr 0.21 Nb 0.11 Ti 0.35 V 0.33 (3.5 wt% H). The electronic structure of the top-performing composition, Cr 0.09 Mg 0.73 Ti 0.18 , was analyzed using density functional theory (DFT) to understand the reasons for its improved hydrogen storage properties compared to TiFe (1.90 wt% H), LaNi 5 (1.37 wt% H) or BCC MPEAs like TiVNbCr (3.70 wt% H). Temperature-dependent molecular dynamics (MD) studies were further performed on optimized BCC MPEAs to qualitatively study hydrogen mobility and analyze the effect of different elemental composition on bulk hydrogen diffusion. Our findings demonstrate how a ML assisted genetic algorithm framework can be used for efficient search of stable, lightweight and cost-effective MPEAs while minimizing the need for expensive ab initio calculations.

DFT↗
DOE OSTI2023

Thermodynamic re-modeling of the Yb-Sb system aided by first-principles calculations

Here, the thermodynamic description of the Yb-Sb binary system is developed by means of the CALculations of PHAse Diagrams (CALPHAD) method by combining experimental data in the literature and predictions from first-principles calculations based on density functional theory (DFT) in the literature and the present work. Two pseudopotentials of Yb are compared in the present DFT-based calculations with 14 and 13 f-electrons frozen in the core, i.e., 5p 6 6s 2 and 5p 6 6s 2 5 d 1 electrons as valence electrons, termed Yb_2 and Yb_3, respectively. It is shown that the phonon spectrum of the YbSb phase calculated using the Yb_3 pseudopotential does not have imaginary phonon modes and is subsequently used to predict its temperature dependent thermodynamic properties by the DFT-based quasiharmonic phonon calculations. The present thermodynamic database includes the Yb 16 Sb 11 phase in addition to five intermetallic phases that were considered in previous modeling studies, i.e., YbSb 2 , YbSb, Yb 11 Sb 10 , Yb 4 Sb 3 , and Yb 5 Sb 3 . The high temperature orthorhombic structure of the Yb 5 Sb 3 phase is not considered in the present work as it was stabilized by hydrogen. The associate solution model is used to describe the short-range ordering behavior in the liquid phase. The calculations from the present thermodynamic model show good agreement with thermochemical and phase equilibrium data from both the present work and the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The u 3 si 2 -H system

U 3 Si 2 is of interest to the nuclear industry as a candidate fuel material due to its high uranium density and high thermal conductivity. However, it has been observed to react with hydrogen, resulting in material decrepitation. As a result, it is important to understand the thermodynamics of the U 3 Si 2 -H system. In this study, the thermodynamics of the hydrogen absorption reaction of USi were determined experimentally using Sievert’s gas absorption and related to crystallographic evolution with hydrogen content using X-ray diffraction. Experimentally-determined thermodynamic parameters were compared with results from density functional theory modeling. Results from this study were also compared with those determined in previous work. Sievert’s gas absorption results were used to develop the pressure-composition-temperature (PCT) curves of the U 3 Si 2 -H system. It was found that the hydride phase exhibited a maximum stoichiometry between U 3 Si 2 H 1.8 and U 3 Si 2 H 2 . The two-phase region for hydride formation from U 3 Si 2 exhibited a miscibility gap with a critical temperature between 623 and 673 K, as calculated from the PCT curves. Analysis of the PCT curves also showed that both the enthalpy and entropy of the hydrogen absorption reaction increased with hydrogen content but were lower than the values for uranium trihydride formation from uranium metal. The enthalpy of reaction for hydrogen absorption was calculated to range between -86.9 and -94.8 kJ mol -1 , while the entropy of reaction was calculated to range between 101.9 and 138.8 J mol -1 K -1 . Furthermore, DFT modeling of the thermoydnamic stability of the USi hydride phases yielded a decomposition temperature of U 3 Si 2 H 2 of approximately 600 K, which was consistent with the experimental results. Similarly, the DFT-calculated enthalpy and entropy of reaction to form USiH were determined to be -106.5kJ mol and 121.8J mol -1 K -1 , respectively, which were both in close agreement with the experimentally-determined values.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Computational design of high entropy alloy coating for hydrogen turbine applications

This project aims to develop novel high entropy alloy (HEA)-based coatings to protect critical components in hydrogen-fueled turbine power system. The HEA-coatings will demonstrate superior performance in hydrogen combustion environment to commercial NiCoCrAlY coating in current natural gas turbine system. The HEA coating facilitates the formation of a protective scale of alpha-alumina to slow down the inward diffusion of oxidizing species and the outward diffusion of metal elements, and possesses ultrahigh corrosion and spallation resistance to prolong the service lifetime of critical components in hydrogen turbine power system. Aimed to accelerate the discovery of novel HEA coating compositions, high throughput computational modeling including CALPHAD and density functional theory and machine learning are performed to predict phase stability, oxygen permeability, oxidation rate constant, coefficient of thermal expansion, and mechanical properties. Based on the modeling and machine learning prediction, experimental validation is performed. Preliminary results will be presented and approaches to minimize oxidation will be discussed.

alloy design↗
DOE OSTI2023

Stability of the interorbital-hopping mechanism for ferromagnetism in multi-orbital Hubbard models

The emergence of insulating ferromagnetic phase in iron oxychalcogenide chain system has been recently argued to be originated by interorbital hopping mechanism. However, the practical conditions for the stability of such mechanism still prevents the observation of ferromagnetic in many materials. Here, we study the stability range of such ferromagnetic phase under modifications in the crystal fields and electronic correlation strength, constructing a theoretical phase diagram. We find a rich emergence of phases, including a ferromagnetic Mott insulator, a ferromagnetic orbital-selective Mott phase, together with antiferromagnetic and ferromagnetic metallic states. We characterize the stability of the ferromagnetic regime in both the Mott insulator and the ferromagnetic orbital-selective Mott phase forms. We identify a large stability range in the phase diagram at both intermediate and strong electronic correlations, demonstrating the capability of the interorbital hopping mechanism in stabilizing ferromagnetic insulators. Our results may enable additional design strategies to expand the relatively small family of known ferromagnetic insulators.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Remote Influence of Andean Convection on Amazonian Rainfall and Its Mechanisms

Models from Coupled Model Intercomparison Project Phase 6 produce too much precipitation over the Andes but too little over the Amazon or the Wet Andes-Dry Amazon (WADA) bias pattern. Unlike the conventional view that convection parameterization and land model deficiencies can contribute to Amazonian rainfall biases, we approach this long-standing biased model behavior through the lens of Andean convection. Using Community Earth System Model v1.1 and focusing on the wet season, our mechanism-denial experiments demonstrate that Andean convection notably reduces precipitation over the Amazon during austral summer. The Andean forced Amazonian response operates on weather timescale. Furthermore, the reduction of Amazonian rainfall is detectable within a few hours after initial Andean forcing. The precipitation response is primarily driven by variations in the moisture budget and is moderated by changes in convective available potential energy over the Amazon. Changes in the total advection of moisture over the Amazon are dominated by the vertical advection term and can be attributed to discrepancies in the dynamic omega field. In the experiments, the Andean east flank region is scrutinized where the vertical velocity and moisture fields play an intermediary role for the Andean driven WADA connection. The Andean forcing induces descending anomalies on the Andean east flank. The disturbances of wind and geopotential fields over the Andean east flank propagate eastward via Kelvin waves. Over the Amazon, descending anomalies and advective drying lead to reduction of mid-to-high level cloud, increase of shortwave cloud forcing and surface net radiation, and enhancement of themodynamic stability and rainfall reduction.

58 GEOSCIENCES↗
DOE OSTI2026

BMINN: Learning chemical potentials and parameters from voltage data for multi-phase battery modeling

Free-energy landscapes and chemical potentials govern the dynamics of phase transitions, transport, and stability in functional materials, yet they remain experimentally inaccessible under realistic operating conditions. Here we introduce a Bayesian model-integrated neural network (BMINN) that embeds physics-based formulations of non-autonomous partial differential-algebraic equations into probabilistic learning. This approach reconstructs hidden thermodynamics directly from macroscopic current-voltage data, providing quantitative access to metastable states, staging transitions, and energy barriers without synchrotron probes. Demonstrated on lithium-graphite electrodes, BMINN recovers full Gibbs free-energy landscapes with fidelity validated against operando X-ray diffraction. The framework generalizes across dynamical regimes, enabling accurate voltage prediction, internal state estimation, and inference of governing parameters. Beyond batteries, BMINN exemplifies a broadly applicable strategy for learning missing physics in multiphase, non-equilibrium systems, offering a new pathway to uncover hidden thermodynamic functions across condensed matter and materials physics.

25 ENERGY STORAGE↗
DOE OSTI2020

Tuning the critical Li intercalation concentrations for MoX 2 bilayer phase transitions using classical and machine learning approaches

Transition metal dichalcogenides (TMDs) such as MoX 2 are known to undergo a structural phase transformation as well as a change in the electronic conductivity upon Li intercalation. These properties make them candidates for charge tunable ion-insertion materials that could be used in electro-chemical devices for neuromorphic computing applications. In this work we study the phase stability and electronic structure of Li-intercalated bilayer MoX 2 with X=S, Se or Te. Using first-principles calculations in combination with classical and machine learning modeling approaches we find that the energy needed to stabilize the conductive phase decreases with increasing atomic mass of the chalcogen atom X. A similar decreasing trend is found in the threshold Li concentration where the structural phase transition takes place. While the electronic conductivity increases with increasing ion concentration at low concentrations, we do not observe a conductivity jump at the phase transition point.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Competing magnetic states in transition metal dichalcogenide moiré materials

Small-twist-angle transition metal dichalcogenide (TMD) heterobilayers develop isolated flat moiré bands that are approximately described by triangular lattice generalized Hubbard models. In this paper we explore the metallic and insulating states that appear under different control conditions at a density of one electron per moiré period and the transitions between them. By combining fully self-consistent Hartree-Fock theory calculations with strong-coupling expansions around the atomic limit, we identify four different magnetic states and one nonmagnetic state near the model phase diagram’s metal-insulator phase transition line. Here, ferromagnetic insulating states, stabilized by nonlocal direct exchange interactions, are surprisingly prominent.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Predicting the recovery and grade of a rougher flotation circuit from liberation data

A flotation model has been validated against the survey data obtained from a copper rougher flotation circuit with a circulating load. Here, the model was developed using the hydrophobic force as a kinetic parameter, which made it possible to use the contact angles of the free and composite particles as input parameters so that the model can predict both recoveries and grades. The model also recognizes the role of bubble coarsening in determining the froth phase recoveries and grades. The model predictions are in good agreement with the survey data. The size-by-size recovery curve as predicted from the model shows an optimum particle size range that is consistent with most of the operating data reported in the literature. The model suggests that the upper particle size limit may be determined by the instability of the froth phase rather than the turbulence in the pulp phase of a flotation cell.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Polysulfide reduction and Li 2 S phase formation in the presence of lithium metal and solid electrolyte interphase layer

Lithium sulfur battery is an attractive next generation technology that can meet many demands of modern society. Unfortunately, the lithium sulfur battery faces unique issues related to the polysulfide shuttle effect, that is due to reduction products dissolving in the electrolyte and their subsequent reduction on the lithium metal electrode. This adds further problems to the already challenging needs of understanding and engineering a solid electrolyte interphase (SEI) layer with desired properties. One of the most important SEI properties is its passivation of lithium metal which is critically important to the overall battery performance. Passivation is difficult to measure experimentally without the influence of many factors. This study reports an investigation of the reduction of the intermediate Li 2 S 8 over lithium already passivated with Li 2 O, Li 2 CO 3 , LiOH, LiF and Li 2 S along with exploration of Li 2 S 8 reduction over pristine lithium nanoclusters using first principles computational models. Significant formation of Li 2 S phase nucleation is found to stabilize the reduction products of the Li 2 S 8 . Here, the formation of Li 2 S is explored in-depth with lithium nanocluster-based models determining a 2 V potential increase for the reduction of polysulfides due to the formation of Li 2 S. This investigation demonstrates passivation effects of important SEI components including Li 2 S.

25 ENERGY STORAGE↗