Search NASA⌕ Search

SEARCH · Search NASA

Results for “EQUILIBRIUM DIAGRAM”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Equation of state and phase diagram of dense hydrogen

The equation of state of hydrogen was calculated for specific volumes ranging from 0.01 to 0.0001 cm3/mole and for temperatures ranging from 200 to 1 million K. Three phases are considered: the molecular solid, the metallic solid and the fluid. Chemical equilibrium between molecules, atoms, ions and electrons is considered in calculating the properties of the fluid phase. Transitions between the three phases will be discussed. The triple point, where the three phases coexist, is calculated to occur at 2.3 Mbar and 1679 K. At higher temperatures and pressures, the molecular solid is unstable.

Kerley, G. I.↗
DOE OSTI2024

Applications of Thermochemical Modeling in Molten Salt Reactors

The extensively evaluated and consistent thermodynamic database, the Molten Salt Thermal Properties Database—Thermochemical (MSTDB-TC), was used along with additional thermodynamic values from other sources as examples of ways to examine molten salt reactor (MSR) fuel behavior. Relative stability with respect to halide potential and temperature for likely fuel and fission product components were mapped in Ellingham diagrams for the chloride and fluoride systems. The Ellingham diagrams provide a rich, visual means for identifying halide-forming components in proposed fuel/solvent salt systems. Thermochemical models and values from MSTDB-TC and ancillary sources were used in global equilibrium calculations to provide compositions for a close analysis of the behavior of a possible Molten Chloride Salt Fast Reactor and a Molten Salt Reactor Experiment-type system at high burnup (100 GWd/t). The results illustrated the oxidative nature of burnup in MSRs and provided information about redox behavior and possible control.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Computing Free Energies with Fluctuation Relations on Quantum Computers

One of the most promising applications for quantum computers is the dynamic simulation of quantum materials. Current hardware, however, sets stringent limitations on how long such simulations can run before decoherence begins to corrupt results. The Jarzynski equality, a fluctuation theorem that allows for the computation of equilibrium free energy differences from an ensemble of short, non-equilibrium dynamics simulations, can make use of such short-time simulations on quantum computers. Here, we present a quantum algorithm based on the Jarzynski equality for computing free energies of quantum materials. We demonstrate our algorithm using the transverse field Ising model on both a quantum simulator and real quantum hardware. As the free energy is a central thermodynamic property that allows one to compute virtually any equilibrium property of a physical system, the ability to perform this algorithm for larger quantum systems in the future has implications for a wide range of applications including the construction of phase diagrams, prediction of transport properties and reaction constants, and computer-aided drug design.

Bassman, Lindsay↗
DOE OSTI2023

Effect of Mn on eutectic phase equilibria in Al-rich Al-Ce-Ni alloys

Microstructural analysis of additively manufactured (AM) Al-Ce-Ni-Mn alloys has identified phases not predicted from existing ternary liquidus projections in the Al-Ce-Ni system. Because the rapid cooling rate of AM is orders of magnitude above that of traditional casting, it is unclear if these additional phases arose from the non-equilibrium processing conditions of AM, a drastic shift in phase stability in the system due to the addition of 1 wt% Mn, or some combination of these two influences. The phases and microstructure of cast samples of Al-Ce-Ni and Al-Ce-Ni-Mn alloys were characterized for several annealing conditions which revealed the equilibrium phases at different temperatures. Phase analysis confirmed that minute levels of Mn substituted for Ni in the system drastically shifts the liquidus projection in the Al-rich corner of the ternary phase diagram such that the eutectic Al 3 Ni phase is suppressed in favor of the Al 23 Ni 6 (Ce,Mn) 4 phase. Further addition of Mn promotes the formation of Al 20 Mn 2 Ce and Al 10 Mn 2 Ce phases. The phase analysis data was then used to improve the CALPHAD modeling of the liquidus projection and isothermal sections for the Al-rich Al-Ce-Ni-Mn quaternary system. Thermodynamic modeling and experimental analysis on phases in the AM sample of Al-Ce-Ni with Mn confirmed that the phases present are consistent with Mn-containing Al-Ce-Ni cast samples. Here, this investigation demonstrates the potential for using secondary alloying elements to drastically alter phase stability and microstructure in alloy systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Chiral effective field theory calculation of neutrino reactions in warm neutron-rich matter

Neutrino scattering and absorption rates of relevance to supernovae and neutron star mergers are obtained from nuclear matter dynamical structure functions that encode many-body effects from nuclear mean fields and correlations. In this study, we employ nuclear interactions from chiral effective field theory to calculate the density, spin, isospin, and spin-isospin response functions of warm β-equilibrium nuclear matter. We include corrections to the single-particle energies in the mean field approximation as well as vertex corrections resummed in the random phase approximation (RPA), including both direct and exchange diagrams. We find that correlations included through the RPA redistribute the strength of the response to higher energy for neutrino absorption and lower energy for antineutrino absorption. This tends to suppress the absorption rate of electron neutrinos across all relevant energy scales. In contrast, the inclusion of RPA correlations enhances the electron antineutrino absorption rate at low energy and suppresses the rate at high energy. These effects are especially important at high-density and in the vicinity of the neutrino decoupling region. Implications for heavy element nucleosynthesis, electromagnetic signatures of compact object mergers, supernova dynamics, and neutrino detection from galactic supernovae are discussed briefly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Thermodynamic Modeling of Point Defects in Triple Conducting Perovskite Ba 0.95 La 0.05 FeO 3- $_δ$ with Incorporation of the Hydride Defect Formation Reaction for Solid Oxide Cells

Distinct from the proton defect, the hydride defect species may be present in certain perovskite materials in reducing environments such as in fuel electrodes of solid oxide cells (SOCs) or on the reducing side of ceramic membranes. A generalized defect thermodynamic model was developed for the triple-conducting perovskites (La,Ba)Fe 1-x M x O 3-δ (M = Y and Zr) to allow inclusion of the hydride defect formation reaction in addition to the other three main defect reactions, namely, the oxygen vacancy formation, hydration, and charge disproportionation reactions. This comprehensive defect model also allows the incorporation of polynomial functional forms of oxygen nonstoichiometry δ to describe the defect reaction energies and entropies and to enable refinements of the defect reaction equilibrium constants in the defect thermodynamic analysis. As a first step, the developed model is applied to the Ba 0.95 La 0.05 FeO 3-δ material as an illustrative system to obtain its Brouwer diagrams with both the proton and hydride defects in relevant SOC conditions, particularly for more reducing environments. In conclusion, the results provide direct guidance on the influence of electronic and ionic defect concentrations upon thermodynamic properties and ultimately on the performance of Ba 0.95 La 0.05 FeO 3-δ and potentially other (La,Ba)Fe 1-x M x O 3-δ perovskite materials involved in SOC applications.

25 ENERGY STORAGE↗
DOE OSTI2024

An ion specific continuum model on ionic surfactant's binary phase diagram in aqueous solution

Here, a continuum aggregation model is proposed to account for specific ion effects, enabling accurate phase diagrams calculations for the micellar, cylindrical and lamellar aggregates of sodium/potassium carboxylate surfactants in aqueous solution across a range of temperatures. Three groups of concentrations at distinctive temperatures are fitted to build empirical temperature dependence given the limited availability of relevant experimental measurements. The specific ion effects are manifested in the aggregates' surface tension as well as in the distributions of counter-ions' concentrations in the vicinity of the aggregates. The aggregates' geometric sizes are well-reproduced. The differential evolution algorithm is applied to address boundary conditions of the electrostatic potential, aggregate size optimization as well as the equilibrium of monomers transferring between the aggregate and aqueous region.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Geochemical Insights Into Volcanic and Lithospheric Evolution of Mars

Introduction: The Martian lithosphere plays a pivotal role in volcanic processes, altering surface features, climate shifts, and the planet's early habitability. Crucially, the rigid lithosphere significantly influences Mars' long-term thermal conditions. Studying various geological eons’ volcanic compositions unveils how the lithosphere has evolved, even though understanding the early Noachian volcanic chemistry remains challenging due to weathering and resurfacing. In this study, newly discovered Noachian volcanic terranes[1-2] along with Hesperian and Amazonian volcanic terranes[3] are considered(Figure 1), to infer the evolution of the Martian lithosphere. Figure 1:Mars topographic map derived from MOLA-HRSC data [4] highlights the study regions including the Noachian volcanic terranes, Arabia Terra (AT), Thaumasia Minor (TM), and the Southwest Margin of Icaria Planum (SWIP)[5].Additionally, the Hesperian volcanic terranes such as Syrtis Major (SM), Hesperia Planum (HP), and, the Amazonian volcanic terranes encompass Elysium Mons (EM), Alba Patera (AP), and Olympus Mons (OM) are marked. Methods and Datasets: We use the most recent elemental mass fraction maps derived from gamma spectroscopy data from the Mars Odyssey 2001mission [6]to ascertain the bulk composition of the Martian landscape. Utilizing the updated geochemical provinces dataset of Mars[7], we consider bulk chemistry variations within our study regions. The maps, at a resolution of 50×50, encompass elements like Al, Ca, Fe, Si, K, Th, H2O, Cl, and S, detected through Gamma-Ray Spectroscopy (GRS) with decimeter scale depth sensitivity[6].Additional major and minor elements like Mg, Na, Ti, P, and Mn, using mass balance methods [8]serve as inputs for melting and crystallization simulations. We employ the pMELTS model with FMQ-3to FMQ oxygen fugacity to estimate the Pressure(P)and melting degree(F)in various eons of volcanic terranes. We have also used conventional thermobarometric calculations to determine their formation P-T conditions, using silica activity and olivine-melt Mg-exchange thermometry[9-10]. Results: Chemical Composition of Studied Volcanic Terranes: To assess formation conditions, we have investigated if GRS-measured compositions signify primary igneous processes by evaluating weathering through the Chemical Index of Alteration (CIA), K/Th ratios, and ternary plots for aqueous alteration. Consistent K/Th ratios and low CIA (<50) across volcanic regions suggest minimal large-scale weathering, affirming well-preserved igneous material extending at decimeter depths. In addition, the calculated bulk compositions align with Martian mantle equilibrium, confirming studied volcanic terranes-Arabia Terra, Thaumasia Minor, and SWIP -as primary or less altered compositions, reinforcing their representativeness at a regional GRS scale on Mars. We have also estimated the bulk compositions of Hesperian and Amazonian volcanic terranes shown in Figure 1. All volcanic terranes represent the primary or minimally altered composition. Temporal Evolution of Martian Volcanic Terranes: Our results show P-T conditions ranging from 1.3-1.6 GPa and 1350-1390°C for Noachian terranes, while Hesperian volcanic terranes vary from 1.6-1.7 GPa and 1370-1395°C, and Amazonian volcanic terranes range between 1.9-2.8 GPa and 1380-1415°C. Our analyses correspond to lithospheric thicknesses/ depth of melting for Noachian which ranges from110-135 km, 120-235 km for Hesperian terranes, and 160-235 km for Amazonian terranes. Partial melting percentages (F) range from 8-12% for Noachian, 10-11% for Hesperian, and 10-12% for Amazonian. We have also calculated mantle potential temperatures (Tp) to show variations across epochs(Figure 2). Heat flux estimation shows variations from 51-65 mW/m2 for Noachian, 38-45 mW/m2for Hesperian, and 27-39 mW/m2for Amazonian terranes, implying temporal changes in lithospheric thickness and heat flow. Comparison of Noachian volcanic terranes with Hesperian and Amazonian reveals variations in lithospheric thickness and heat flux, indicating potential implications for Martian surface conditions, including volcanic activity, climate evolution, and early habitability. Figure 2: P-T diagram displays formation conditions of studied Martian volcanic regions using GRS data, alongside Noachian-age surface basalt samples from Gusev, Gale, and Jezero. Dashed lines represent Martian mantle potential temperatures, calculated by adding the latent heat of fusion to the equilibrium temperature at which melt and solid mantle coexist and subtracting the Martian adiabatic gradient corresponding to the depth of melt. The solid line shows the Martian mantle solidus [11].Discussion& Implication: The study explores Martian volcanic activity across eons. Uniform lithospheric thickness from the Noachian to Hesperian eons suggests persistent weak plumes until Hesperian, influencing heat flow variations. Noachian's high heat flux, likely due to concentrated radioactive elements, shows spatial heterogeneity, indicating early Mars' differentiated mantle. Volcanic outgassing plays avital role in the formation of the Martian climate, likely sustained water on the surface, influencing warm-wet or cold-dry conditions. Hesperian's favorable lithospheric heat flow (45 mW/m2) supports extended water existence, resembling Earth's conditions. Whereas, Amazonian's diminished heat flux (25 mW/m2) corresponds to water loss and arid climate. Hesperian's similar lithospheric thickness to mid-Noachian implies prolonged water retention, influencing climate and potential for astrobiological research, expanding perspectives for future Martian missions. Future work: As volcanism was active throughout its history. The emergence of Mars' extensive volcanic regions, known as large igneous provinces, which are linked to intermittent volcanic activity driven by mantle plumes, poses challenges to understanding the scale and persistence of these phenomena in Martian convection models. Therefore, understanding the regional-scale changes in eruptive processes within Martian volcanic provinces over time remains unclear and crucial for deciphering the evolution of the Martian interior without Earth-like plate tectonics. We will continue this work with petrological and thermoelastic analyses to study the compositional and spatiotemporal evolution of the Elysium Volcanic Province.

A Rani↗
DOE OSTI2026

Realization of two-dimensional discrete time crystals with anisotropic Heisenberg coupling

A discrete time crystal (DTC) is an out-of-equilibrium phase of matter that spontaneously breaks discrete time-translation symmetry. Previous studies have been limited to a set of models with Ising-like couplings - and mostly only in one dimension - thus precluding our understanding of the existence (or not) of DTCs in models with more realistic interactions. In this work, by combining the latest generation of IBM quantum processors with state-of-the-art tensor network methods, we demonstrate the existence of a DTC in a two-dimensional system governed by anisotropic Heisenberg interactions. We uncover a rich phase diagram encompassing spin-glass, ergodic, and time-crystalline phases, and identify the interplay of initialization, interaction anisotropy, and driving protocols in stabilizing the DTC phase. By extending the study of Floquet matter beyond simplified models, we lay the groundwork for exploring how driven systems bridge the gap between quantum coherence and emergent non-equilibrium thermodynamics.

Phase transitions and critical phenomena↗
DOE OSTI2022

Transferability of data-driven, many-body models for CO 2 simulations in the vapor and liquid phases

Here, extending on the previous work by Riera et al. [J. Chem. Theory Comput. 16, 2246–2257 (2020)], we introduce a second generation family of data-driven many-body MB-nrg models for CO 2 and systematically assess how the strength and anisotropy of the CO 2 –CO 2 interactions affect the models’ ability to predict vapor, liquid, and vapor–liquid equilibrium properties. Building upon the many-body expansion formalism, we construct a series of MB-nrg models by fitting one-body and two-body reference energies calculated at the coupled cluster level of theory for large monomer and dimer training sets. Advancing from the first generation models, we employ the charge model 5 scheme to determine the atomic charges and systematically scale the two-body energies to obtain more accurate descriptions of vapor, liquid, and vapor–liquid equilibrium properties. Challenges in model construction arise due to the anisotropic nature and small magnitude of the interaction energies in CO 2 , calling for the necessity of highly accurate descriptions of the multidimensional energy landscape of liquid CO 2 . These findings emphasize the key role played by the training set quality in the development of transferable, data-driven models, which, accurately representing high-dimensional many-body effects, can enable predictive computer simulations of molecular fluids across the entire phase diagram.

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

Activity coefficient acquisition with thermodynamics–informed active learning for phase diagram construction

This work explores the use of thermodynamics-informed Gaussian processes (GPs) and active learning (AL) to model activity coefficients and construct phase diagrams. Relying on synthetic data generated from an excess Gibbs energy model, GPs were found to accurately describe the activity coefficients of several binary mixtures across large composition and temperature ranges. Moreover, GPs could estimate their own uncertainty and identify composition/temperature regions where activity coefficient data provide the most information to the models. This was leveraged to build AL algorithms targeted at modeling phase equilibria. In many cases, a single active-learning-acquired data point was sufficient to describe the phase diagrams studied. Lastly, the ability of AL to greatly reduce the amount of data needed to obtain accurate models was further verified on experimental case studies, namely individual ion activity coefficients, the solid–liquid and vapor–liquid equilibrium of deep eutectic solvents, and phase equilibria in ternary mixtures.

25 ENERGY STORAGE↗
DOE OSTI2021

Constraining the nonanalytic terms in the isospin-asymmetry expansion of the nuclear equation of state

In this report we examine the properties of the isospin-asymmetry expansion of the nuclear equation of state from chiral two- and three-body forces. We focus on extracting the high-order symmetry energy coefficients that consist of both normal terms (occurring with even powers of the isospin asymmetry) as well as terms involving the logarithm of the isospin asymmetry that are formally nonanalytic around the expansion point of isospin-symmetric nuclear matter. These coefficients are extracted from numerically precise perturbation theory calculations of the equation of state coupled with a new set of finite difference formulas that achieve stability by explicitly removing the effects of higher-order terms in the expansion. We consider contributions to the symmetry energy coefficients from both two- and three-body interactions. It is found that the coefficients of the logarithmic terms are generically larger in magnitude than those of the normal terms from second-order perturbation theory diagrams, but overall the normal terms give larger contributions to the ground state energy. The high-order isospin-asymmetry terms are especially relevant at large densities where they affect the proton fraction in β-equilibrium matter, and in particular we find that at twice saturation density they can reduce the proton fraction by up to 0.02.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1988

The solubility of hen egg-white lysozyme

The equilibrium solubility of chicken egg-white lysozyme in the presence of crystalline solid state was determined as a function of NaCl concentration, pH, and temperature. The solubility curves obtained represent a region of the lysozyme phase diagram. This diagram makes it possible to determine the supersaturation of a given set of conditions or to achieve identical supersaturations by different combinations of parameters. The temperature dependence of the solubility permits the evaluation of Delta-H of crystallization. The data indicate a negative heat of crystallization for the tetragonal crystal form but a positive heat of crystallization for the high-temperature orthorhombic form.

Howard, Sandra B.↗
NASA NTRS2009

Long Term Measurement of the Vapor Pressure of Gold in the Au-C System

Incorporating the {Au(s,l) + graphite} reference in component activity measurements made with the multiple effusion-cell vapor source mass spectrometry (multicell KEMS) technique provides a fixed temperature defining ITS-90 (T(sub mp)(Au) = 1337.33K) and a systematic method to check accuracy. Over a 2 year period delta H sub(298)Au was determined by the 2nd and 3rd law methods in 25 separate experiments and were in the ranges 362.2 plus or minus 3.3 kJmol(sup -1) and 367.8 plus or minus 1.1 kJmol(sup -1), respectively. This 5 kJmol-1 discrepancy is transferred directly to the measured activities. This is unacceptable and the source of this discrepancy needs to be understood and corrected. Accepting the 2nd law value increases p(Au) by about 50 percent, brings the 2nd and 3rd law values into agreement and removes the T dependence in the 3rd law values. While compelling, there is no way to independently determine instrument sensitivities, S(sub Au), with T in a single experiment with KEMS. This lack of capability is stopping a deeper understanding of this problem. In addition, the Au-C phase diagram suggests a eutectic invariant reaction: L-Au(4.7at%C) = FCC-Au(0.08at%C) + C(graphite) at T(sub e) approximately 1323K. This high C concentration in Au(l) must reduce p(Au) in equilibrium with {Au(s,l) + graphite} and raises some critical questions about the Gibbs free energy functions of Au(s,l) and the Au fixed point (T(sub mp)(Au) = 1337.33K) which is always measured in graphite.

Copland, Evan H.↗
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↗
NASA NTRS1981

Momentum and energy balance in late-type stellar winds

Observations at ultraviolet and X-ray wavelengths indicate that the classical picture of a static stellar atmosphere containing a radiative equilibrium temperature distribution is inapplicable to the majority of late type stars. Mass loss and the presence of atmospheric regions characterized by gas temperatures in excess of the stellar effective temperature appear to be almost ubiquitous throughout the HR diagram. Evidence pertaining to the thermal and dynamical structure of the outer envelopes of cool stars is summarized. These results are compared with the predictions of several theoretical models which were proposed to account for mass loss from latetype stars. Models in which the outflow is thermally radiatively, or wave driven are considered for identification of the physical processes responsible for the observed wind properties. The observed variation of both the wind, thermal and dynamical structure as one proceeds from the supergiant branch toward the main sequence in the cool portion of the HR diagram give consideration to potential mechanisms for heating and cooling the flow from low gravity stars.

Macgregor, K. B.↗
DOE OSTI2026

Computational Design of Interlayers for Thermally Stable Compositionally Graded Coatings on Nickel Alloys

To extend the service life of Ni-based superalloys, refractory metal coatings are often used. However, direct bonding between metals with dissimilar crystal structure promotes brittle intermetallic phase formation. This work presents a computational thermodynamic framework for high throughput design of functionally graded interlayers to suppress deleterious phases that may form at the interlayer. The Thermo-Calc software package was used to screen candidate metallic interlayer elements based on stability of solid-solution phases. Vanadium was identified as a promising interlayer due to its consistent suppression of intermetallic phases. Temperature-dependent phase diagram mapping between 600 and 1000 °C guided selection of a compositional pathway that significantly reduced intermetallic formation compared to directly joining the Ni-based and Nb refractory alloys. Time–temperature–transformation analysis was performed to assess whether equilibrium-predicted phases are kinetically accessible along regions of the graded path where non-solid-solution phases are not fully suppressed. The methodology was further applied to additional Ni-based alloy and coating systems, illustrating its transferability as an approach for rapid computational design of graded interlayers in dissimilar high-temperature materials.

36 MATERIALS SCIENCE↗