Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic phase”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Effects of surface chemical potentials on cation segregation

Surface cation segregation on perovskite-type electrodes is one of the major issues limiting the durability of high-temperature solid-oxide electrochemical cells, and this process is strongly dependent on temperature, the external gas-environment, and impurities in air, such as CO 2 . Here cation segregation on La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-δ (LSCF) is systematically evaluated under a set of chemical potentials to determine the origin of surface segregation and the dominant factors that govern the segregation process. Temperature is the main driving force for surface strontium segregation as the thermodynamic stability of each phase varies, and the SrO particles appear primarily in a specific operating window. In addition to the required thermal energy, the presence of gaseous oxygen-containing molecules as reactants help drive the precipitation of SrO. Oxygen partial pressure (pO 2 ) controls the defect chemistry of LSCF, leading to promotion or suppression of surface segregation. The presence of CO 2 promotes the nucleation process and suppresses the surface migration step, significantly altering the surface morphology. We also show that A-site deficiency can limit the SrO segregation in certain conditions but shows no effect in others. This study reveals the impact of gas–solid interactions on surface segregation and highlights the subtle relationship between multiple segregation driving forces.

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

Development of uncertainty-aware equation-of-state models: Application to copper

Sophisticated hydrodynamic codes are commonly used to understand and predict events relevant to natural and applied sciences. The degree to which these simulations reflect reality, however, is dependent on how well we understand the materials and underlying physics involved. These research communities need material models that communicate the uncertainty in the physical properties, which at their basest form comes from the uncertainty in the underlying experimental measurements. Additionally, we have constructed a new framework for using both experimental measurements and the associated experimental uncertainties to build equation-of-state models that reflect not only current best measurements but also the accuracy of those measurements. This method had been used to construct an ensemble of equation-of-state models for copper that communicates the experimental uncertainties in the data through the equation-of-state model, which is available for application in any simulation metric of interest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Effect of atom substitutions on the magnetic properties in Ce 2 Fe 17 : Toward permanent magnet applications

Due to the rapidly developing technologies and huge market demand, there has been increasing interest internationally in exploring permanent magnet formulations in addition to the well-known Nd 2 Fe 14 B and SmCo 5 /Sm 2 Co 17 . Given Fe's low materials cost and generally high magnetization, Fe-rich rare earth binaries such as Ce 2 Fe 17 comprise a rich “hunting ground” for such new materials. While this compound suffers from a low ordering point and is a heli magnet, these difficulties are easily remedied by the substitution of appropriate amounts of cobalt for Fe, with room-temperature saturation magnetization as high as 1.5 T. In this work, we try to switch the all-important magnetic anisotropy from planar to uniaxial behavior in Ce 2 Fe 17 via 18 h - and 6 c -type atom substitutions with Si, Ir, and numerous other atoms. The uniaxial magnetocrystalline anisotropy is successfully achieved in the 6 c -site-substituted Ce 2 Fe 15 Ir 2 systems, along with large magnetization. We find that iridium substitution, in particular, induces a substantial uniaxial anisotropy of 11.25 MJ/m 3 , which is comparable to most of the current rare earth permanent magnets. Although the iridium substitution is costly, the finding of Ir-triggered uniaxial magnetic anisotropy indicates the potential of Ce–Fe-based alloys for permanent magnets.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Borderline first-order phase transition and large cryogenic magnetocaloric effect in PrNdIn

We report a large cryogenic magnetocaloric effect stemming from an unconventional borderline first-order magnetic phase transition with negligibly small thermomagnetic hysteresis in a rare-earth-based intermetallic compound PrNdIn. The sample exhibits maximum magnetic field-induced entropy change as large as –10 J/Kg K (for 20 kOe magnetic field change) near the boiling point of oxygen. Magnetocaloric properties of PrNdIn are comparable to those of other known potential magnetocaloric materials with operating temperatures ranging between 50 and 125 K. Here, the magnetic properties of the present sample are qualitatively reminiscent of those of the binary Pr 2 In and Nd 2 In, including the emergence of a second low-temperature anomaly in the temperature dependence of magnetization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

An analytic and complete equation of state for condensed phase materials

Analytic equations of state (EOS) are intended to reproduce theoretical and experimental data in a single phase portion of the thermodynamic space. We devise a complete and thermodynamically consistent model with four distinct features: (1) a reference isotherm that remains thermodynamically stable, (2) a flexible specific heat model based on a fourth-order rational polynomial, (3) a Grüneisen parameter that depends on specific volume and temperature, and (4) pressure and internal energy functions that can be inverted analytically in temperature. The model aims to improve the accuracy of existing equations of state while remaining computationally efficient. To demonstrate its features, we include calibrations for single-crystal pentaerythritol tetranitrate (PETN), liquid nitromethane (NM), and hexagonal close-packed beryllium (Be) metal. The parameter optimization uses the specific heat capacity, Grüneisen parameter, and static compression curves obtained from density functional theory for the crystalline solids and molecular dynamics simulations for liquid NM. We also present a velocity autocorrelation function that yields accurate phonon densities of states for the EOS calibration from the molecular dynamics trajectories. Each of the three calibrations is constrained to enforce the ambient state from experimental measurements and validated against experimental Hugoniot data from multiple sources. We also include one-dimensional hydrodynamic simulations of the isentropic compression experiments for beryllium conducted at the Z facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Noncollinear polar magnet Fe 2 (SeO 3 ) 3 (H 2 O) 3 with inequivalent Fe 3+ sites

The emergence of novel magnetic states becomes more likely when the inversion symmetry of the crystal field, relative to the center between two spins, is broken. We propose that placing magnetic spins in inequivalent sites in a polar lattice can promote a realization of nontrivial magnetic states and associated magnetic properties. To test our hypothesis, we study Fe 2 (SeO 3 )(H 2 O) 3 as a model system that displays two distinct Fe(1) and Fe(2) magnetic sites in a polar structure (R3c space group). At low fields μ 0 H≤ 0.06 T, the material undergoes an antiferromagnetic ordering with T N1 = 77 K and a second transition at T N2 ≈ 4 K. At μ 0 H≥ 0.06 T and 74 K ≤T≤ 76 K, a positive entropy change of ∼0.12 mJ mol −1 K −1 can be associated with a metamagnetic transition to possibly nontrivial spin states. At zero field, Fe(1) is nearly fully ordered at T≈ 25 K, while Fe(2) features magnetic frustration down to T = 4 K. The magnetic ground state, a result corroborated by single-crystal neutron diffraction and 57 Fe Mössbauer spectroscopy, is a noncollinear antiparallel arrangement of ferrimagnetic Fe(1)–Fe(2) dimers along the c-axis. The results demonstrate that placing distinct magnetic sites in a polar crystal lattice can enable a new pathway to modifying spin, orbital, and lattice degrees of freedom for unconventional magnetism.

Crystal lattices↗
DOE OSTI2025

High temperature melting of dense molecular hydrogen from machine-learning interatomic potentials trained on quantum Monte Carlo

We present results and discuss methods for computing the melting temperature of dense molecular hydrogen using a machine learned model trained on quantum Monte Carlo data. In this newly trained model, we emphasize the importance of accurate total energies in the training. We integrate a two phase method for estimating the melting temperature with estimates from the Clausius–Clapeyron relation to provide a more accurate melting curve from the model. We make detailed predictions of the melting temperature, solid and liquid volumes, latent heat, and internal energy from 50 to 180 GPa for both classical hydrogen and quantum hydrogen. At pressures of roughly 173 GPa and 1635 K, we observe molecular dissociation in the liquid phase. Here, we compare with previous simulations and experimental measurements.

08 HYDROGEN↗
DOE OSTI2025

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗
DOE OSTI2020

Weak scaling of the contact distance between two fluctuating interfaces with system size

A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers-of-mass remain separated by a nonzero distance ℓ. An example of such a situation is the motion of interfaces between two phases at conditions of thermodynamic coexistence, and in particular the annihilation of domain wall pairs under periodic boundary conditions. Here we present a general approach to calculate the probability distribution of the contact distance ℓ and determine how its most likely value ℓ* depends on the surfaces' lateral size L. Using the Edward-Wilkinson equation as a model for interfaces, we demonstrate that ℓ* scales weakly with system size, i.e., the dependence of ℓ* on L for both (1+1)- and (2+1)-dimensional interfaces is such that lim L→∞ (ℓ*/L) = 0. In particular, for (2+1)-dimensional interfaces ℓ* is an algebraic function of logL, a result that is confirmed by computer simulations of slab-shaped domains formed under periodic boundary conditions. Overall, this weak scaling implies that such domains remain topologically intact until ℓ becomes very small compared to the lateral size of the interface, contradicting expectations from equilibrium thermodynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Ab initio free energies of liquid metal alloys: Application to the phase diagrams of Li-Na and K-Na

Comparison of free energies between different phases and different compositions underlies the prediction of alloy phase diagrams. To allow direct comparison, consistent reference points for the energies or enthalpies are required, and the entropy must be placed on an absolute scale, yielding absolute free energies. Here we derive absolute free energies of liquids from ab-initio molecular dynamics (AIMD) by combining the directly simulated enthalpies with an entropy derived from simulated densities and pair correlation functions. Additionally, as an example of the power of this method we calculate the phase diagrams of two binary alkali metal alloys, Li-Na and K-Na, revealing a critical point and liquid-liquid phase separation in the former case, and a deep eutectic in the latter. Good agreement with experimental data demonstrates the power of this simple method.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Chemomechanics in alloy phase stability

We describe a first-principles statistical mechanics method to calculate the free energies of crystalline alloys that depend on temperature, composition, and strain. The approach relies on an extension of the alloy cluster expansion to include an explicit dependence on homogeneous strain in addition to site occupation variables that track the degree of chemical ordering. The method is applied to the Si-Ge binary alloy and is used to calculate free energies that describe phase stability under arbitrary epitaxial constraints. We find that while the incoherent phase diagram (in which coexisting phases are not affected by coherency constraints) hosts a miscibility gap, coherent phase equilibrium predicts ordering and negative enthalpies of mixing. Instead of chemical instability, the chemomechanical free energy exhibits instabilities along directions that couple the composition of the alloy with a volumetric strain order parameter. Furthermore, this has fundamental implications for phase field models of spinodal decomposition as it indicates the importance of gradient energy coefficients that couple gradients in composition with gradients in strain.

Materials Science↗
DOE OSTI2025

Superstructural phase transitions in polymer-grafted nanooctahedra

Superlattices of polyhedral nanocrystals exhibit emergent properties defined by their structural arrangements, but native nanocrystal ligands often limit their programmability. Polymeric ligands address this limitation by enabling tunable nanocrystal softness through modifications of polymer molecular weight and grafting density. Here, we investigate phase transitions in polymer-grafted nanooctahedra by varying polymer length, nanocrystal size, truncation, and ligand density. In two-dimensional superlattices, longer polymers or smaller nanooctahedra induce a transition from orientationally ordered to hexagonal rotator lattices. In three-dimensional superlattices, increasing polymer length drives transitions from Minkowski to body-centered cubic and plastic hexagonal close-packed phases, while higher grafting densities further enable transitions to simple hexagonal phases. Polymer brush and thermodynamic perturbation theories, supported by Monte Carlo simulations, uncover the entropic and enthalpic forces that govern these transitions. This work highlights the versatility of polymer-grafted anisotropic nanocrystals as building blocks for designing hierarchical superstructures and metamaterials with customizable properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Reconsidering short-range order in complex concentrated alloys

Abstract The seemingly contradictory state of research on short-range order in many-component alloys is addressed through a critical review of the characterization of face-centered-cubic 3 d systems. Despite the paucity of direct observations, the ordering of many widely studied alloys is argued to be primarily interesting for its potential ubiquity. To clarify this situation, future research directions are proposed with reference to historical results, including a review of the fundamental principles of ordering and clustering. Graphical abstract

36 MATERIALS SCIENCE↗
DOE OSTI2021

Developing Practical Models of Complex Salts for Molten Salt Reactors

Molten salt reactors (MSRs) utilize salts as coolant or as the fuel and coolant together with fissile isotopes dissolved in the salt. It is necessary to therefore understand the behavior of the salts to effectively design, operate, and regulate such reactors, and thus there is a need for thermodynamic models for the salt systems. Molten salts, however, are difficult to represent as they exhibit short-range order that is dependent on both composition and temperature. A widely useful approach is the modified quasichemical model in the quadruplet approximation that provides for consideration of first- and second-nearest-neighbor coordination and interactions. Its use in the CALPHAD approach to system modeling requires fitting parameters using standard thermodynamic data such as phase equilibria, heat capacity, and others. A shortcoming of the model is its inability to directly vary coordination numbers with composition or temperature. Another issue is the difficulty in fitting model parameters using regression methods without already having very good initial values. The proposed paper will discuss these issues and note some practical methods for the effective generation of useful models.

Besmann, Theodore M. (ORCID:0000000155980550)↗
DOE OSTI2021

Visualizing temperature-dependent phase stability in high entropy alloys

High entropy alloys (HEAs) contain near equimolar amounts of five or more elements and are a compelling space for materials design. In the design of HEAs, great emphasis is placed on identifying thermodynamic conditions for single-phase and multi-phase stability regions, but this process is hindered by the difficulty of navigating stability relationships in high-component spaces. Traditional phase diagrams use barycentric coordinates to represent composition axes, which require ( N – 1) spatial dimensions to represent an N -component system, meaning that HEA systems with N > 4 components cannot be readily visualized. Here, we propose forgoing barycentric composition axes in favor of two energy axes: a formation-energy axis and a ‘reaction energy’ axis. These Inverse Hull Webs offer an information-dense 2D representation that successfully captures complex phase stability relationships in N ≥ 5 component systems. We use our proposed diagrams to visualize the transition of HEA solid-solutions from high-temperature stability to metastability upon quenching, and identify important thermodynamic features that are correlated with the persistence or decomposition of metastable HEAs.

Chemistry↗
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, θ'', θ' 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 θ', as those of the other phases were negligible. The simulation predictions of the phase boundaries between α-Al and either θ'', θ' 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 θ'' at low temperatures. They also demonstrate that θ' 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 θ'. 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↗