Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase equilibria”

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.

48 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

First-principles thermodynamic assessment of Sr-containing secondary phase formation in strontium-substituted lanthanum manganites for solid oxide cell applications

Sr-secondary phase formation is a potentially significant degradation mode with direct impact upon solid-oxide cell (SOC) commercial viability. A first-principles based thermodynamic study was performed for La 1−x Sr x MnO 3±δ (LSM) perovskites to assess their stability against formation of different Sr-secondary phases, including SrO, SrCrO 4 , SrSO 4 , SrCO 3 , and Sr(OH) 2 , for SOC applications. The Sr-secondary phase formation reaction free energies were determined via a thermodynamic model by combining ab initio lattice dynamics calculations for the solid phases and ab initio thermodynamic data for the gas phases. The current approach expands the previously reported thermodynamic modeling studies by integrating first-principles based point defect equilibria into the thermodynamic analysis. The modeling results obtained using this new approach indicate an increased tendency to form the SrO oxide upon decreasing the oxygen partial pressure. Additionally, the enhancing factors to form the Sr-related secondary phase from the associated SrO activity in LSM are further quantified by considering the equilibrium of SrO reacting with the contaminant gas species as a function of temperature and pressure.

defect thermodynamics modeling↗
DOE OSTI2024

Generalized Brunauer–Emmett–Teller Isotherm for Mixed-Gas Multilayer Adsorption Equilibria

Here, this work presents a rigorous thermodynamic framework to calculate mixed-gas multilayer adsorption equilibria from single-component isotherms and vapor–liquid equilibria. Named the generalized Brunauer–Emmett–Teller (gBET) isotherm, the newly formulated isotherm considers adsorbent surface heterogeneity, competitive adsorption on the monolayer, condensation–evaporation on the subsequent layers, and adsorbed phase nonideality for the monolayer and the subsequent layers. The monolayer adsorbed phase nonideality is tracked by using an area-based adsorption nonrandom two-liquid activity coefficient model. The adsorbed phase composition and corresponding nonideality in the subsequent layers are calculated at the dew point condition of the mixed gas with either an equation of state or an activity coefficient model for the vapor–liquid equilibria. The proposed model is validated with six single and three binary multilayer adsorption equilibrium systems, and the model results are compared against those from the classical BET isotherm for single-component adsorption and ideal adsorbed solution theory for mixed-gas adsorption equilibria.

09 BIOMASS FUELS↗
DOE OSTI2026

Toroidal phase of post-disruption runaway electron loss to wall in presence of applied 3-D fields in DIII-D

The toroidal distribution of runaway electrons (REs) striking the centre post in DIII-D post-disruption RE final loss events is simulated utilizing the MARS-F code combined with the REORBIT module, by tracing the guiding-centre drift orbits of test REs in the presence of external (applied) plus internal (intrinsic instability) 3-D fields. To better recover the experimental results, three different equilibria with different safety factor profiles are adopted. Dominant resistive magnetohydrodynamic instabilities are found to be m/n=2/1, 3/2 and 1/1 modes (where m and n are the poloidal and toroidal mode numbers). Externally applied n=1 resonant magnetic perturbations (RMPs), with both even and odd parity configuration, are considered in each equilibrium. Locking of MHD instability to the applied RMPs is simulated by minimizing the total perturbed magnetic energy in the plasma. Here, the simulated toroidal locations, where most REs strike in the presence of an edge m/n= 2/1 instability and the n=1 RMPs in even parity, are found to agree well with experiments. Modeling also captures the measured n = 2 RE loss pattern on the wall for experiments where an edge m/n=3/2 instability is dominant instead of the m/n=2/1 mode. The key overall findings include (a) the measured global RE toroidal impact pattern in experiments can be explained by a locked resistive kink instability, (b) the toroidal phase of the locked mode is well predicted by minimizing the total perturbed magnetic energy, and (c) this edge locked mode determines the toroidal phase of the peak RE wall impact.

Bai, Xue [University of California San Diego, La J↗
DOE OSTI2024

Salt Hydrate Eutectic Thermal Energy Storage for Building Thermal Regulation (Final Technical Report)

Thermal energy storage is anticipated to play an important role in developing the power grid of the future - a power grid that meets increasing demands of users, is resistant to disruptions, but also allows for greater penetration of renewable resources. Specifically, thermal energy storage materials can be integrated into HVAC systems and building envelopes, where they can be used to shift power demands for building climate control from periods of peak demand to periods of low demand. Phase change materials (PCMs) are compelling as low-cost, high energy density thermal energy storage materials for building thermal management. However, there is a lack of high performance low-cost PCMs within the specific temperature ranges which would most effectively allow for power load shifting. Inorganic salt hydrates represent a promising class of PCMs, but their inherent limitations cause them to be currently unavailable for reliable building applications. The overarching goals of this research effort are to: 1) Discover low-cost, high volumetric density salt hydrate eutectic PCMs to store low-quality heat (10 to 40 °C); 2) Introduce a high thermal conductivity matrix to reduce the time constant for energy storage to ~0.1 to 1 hr, incorporate nucleation catalysts to decrease undercooling, and utilize microencapsulation and shape stabilization approaches, minimizing moisture loss/gain, mitigating phase separation, and maintaining stable melting behavior over the lifetime of the compounds; 3) Evaluate the impact these systems have on peak load shifting, and the potential for overall energy savings under different climatic scenarios and building configurations. These goals will be achieved by an integrated research program consisting of six cohesive research subtasks: 1) Materials discovery of eutectic salt hydrate PCMs by using computationally predicted thermodynamic equilibria, coupled with high-throughput experimental validation, 2) Rapid experimental screening of nucleation catalysts identified through robust computational databases, 3) Embedding salt hydrate PCM into a low cost and scalable high conductivity matrix, 4) Microencapsulation of salt hydrate microspheres using hybrid inorganic-polymer microencapsulation approach, 5) Shape stabilization by thermoreversible salt hydrate salogels, and 6) Analysis of end-use using thermal simulations, and characterization of mock-up energy storage finished components.

25 ENERGY STORAGE↗
DOE OSTI2026

Hydrogen defects in LaBi 2 O 4 X (X = Cl, Br, and I) Sillén oxyhalide phases and their impacts on ionic transport

Sillén oxyhalides have recently emerged as promising materials for both photocatalytic and ionic transport applications, yet the role of likely-ubiquitous hydrogen-related defects in these layered compounds remains largely unexplored. Here, we employ first-principles defect calculations to investigate incorporation energetics for hydrogen- and oxygen-related defects, as well as their migration barriers in LaBi 2 O 4 X (X = Cl, Br, I) phases. We find that hydrogen interstitials, particularly protonic species (H i + ), are readily accommodated within the open Bi–O layers. Protons compete with oxygen vacancy donors (V O 2+ ) and charge-compensate with oxygen interstitial acceptors (O i 2− ). By linking hydrogen defect formation to water- and oxygen-related redox equilibria, we reveal that V O 2+ facilitates H i + incorporation, while O i 2− promotes interstitial hydroxide formation, establishing a direct connection between proton and oxide-ion transport. Calculated migration barriers indicate that ionic diffusion is confined to Bi–O layers with low barriers of 0.20–0.25 eV for H i + and 0.14–0.25 eV for V O 2+ , suggesting that the materials contain intrinsic pathways for mixed ionic conduction. These results provide a microscopic picture of hydrogen behavior in Sillén oxyhalides and point to design strategies for integrating protonic and oxide-ion transport in layered oxyhalide electrolytes. Band-edge alignment analysis shows that LaBi 2 O 4 I provides the optimal combination of hydrogen solubility, oxygen defect stability, and mixed ionic conductivity, highlighting its potential for low-temperature electrochemical and energy-conversion applications. Overall, this work establishes the defect-driven origin of hydrogen transport in Sillén oxyhalides and expands their applicability beyond photocatalysis to mixed ionic conduction and hydrogen electrochemistry.

Energy - Conversion↗
DOE OSTI2025

The Temperature Dependence of the Langmuir Adsorption Model for a Single-Site Metal–Organic Framework

The single-site Langmuir adsorption model, also known as the Langmuir isotherm equation, is one of the simplest possible descriptions of adsorption phenomena, and yet finds widespread applicability across a range of disciplines. In its simplest form, it is deployed to treat adsorption equilibria at constant temperature (i.e., along isotherms); however, at the heart of its derivation is a more general class of models that each incorporate an explicit temperature dependence, subject to assumptions about the spatial/translational degrees of freedom of the adsorbed species. In this work, measurements of the temperature dependence of supercritical adsorption of H 2 on a single-site metal-organic framework (MOF) are presented, and fitted using a range of Langmuir models with distinct treatments of degrees of freedom in the adsorbed phase. Surprisingly, all of the models can be used to adequately represent the measured data, despite yielding significantly different values for binding energy and the temperature dependence of the isosteric enthalpy of adsorption (i.e., the isosteric heat, qst). Furthermore, a critical finding is that the mean-temperature isosteric enthalpy of adsorption remains consistent across all models within experimental error, highlighting its reliability for evaluating adsorption thermodynamics.

08 HYDROGEN↗
DOE OSTI2025

Simulations of saturated MHD activity in the HBT-EP tokamak

The NIMROD code is used to perform simulations of MHD activity in the HBT-EP tokamak, including the effect of a resistive wall. Linear simulations are conducted to scan current and pressure profiles to establish self-consistent baseline equilibria that result in MHD instabilities with low error between experimental diagnostics and synthetic equivalents. A family of such equilibria is found with relatively peaked current and pressure profiles. Nonlinear simulations are performed using the optimal equilibrium, which is linearly unstable to a Resistive Wall Mode—stable in the ideal wall limit and with an Alfvénic growth rate in the no-wall limit. Using resistive wall boundary conditions, nonlinear Hall MHD simulations yield a rotating, saturated n=1 mode similar to the type observed experimentally in HBT-EP discharges. An n=1 mode around the q=2 surface mediates a localized axisymmetric perturbation that flattens the plasma current from a linearly unstable state to a 3D, rotating, stable state. During the saturated period, periodic changes in inter-diagnostic phase lag allow for qualitative estimations of effective, macroscopic transport coefficients used in MHD closure models and could be extended to model post-disruption MHD structure rotation after thermal quenches seen in HBT-EP.

Arnold, David A. (ORCID:0009000059535761)↗
DOE OSTI2024

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

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

first-principles calculations↗
DOE OSTI2025

Integrated modeling of advanced tokamaks for tearing mode avoidance on DIII-D

Understanding the physics of the onset and evolution of tearing modes (TMs) in tokamak plasmas is crucial for high-performance steady-state operations. The onset condition of (m,n) = (2,1) tearing stability, $Δ′$ > $Δ_c^′$> 0, has been studied in DIII-D steady-state hybrid discharges with accurate equilibrium reconstruction and well-measured plasma profiles. Here, m is the poloidal mode number and n is the toroidal mode number. The onset of n = 1 TMs is carefully identified using the Mirnov array to determine in-phase signals for the onset time and using electron cyclotron emission measurements to trace back from the saturated island width in the experiments for both initial island width and its growth rate. The tearing stability index $Δ′$ at TM onset is calculated in multiple ways, using the modified Rutherford equation and using the PEST-III or resistive DCON magnetohydrodynamic codes, and the results are in reasonable agreement with each other. The calculated $Δ′$ is compared with an analytical formula of the tearing stability threshold $Δ_c^′$ to determine the condition for the mode onset, $Δ′$ > $Δ_c^′$. Finally, the local gradient effect of the plasma current and the pressure profiles on the tearing stability is investigated in terms of the difference between $Δ′$ and $Δ_c^′$ through a series of modeled equilibria with the parametric variations. A database of experimental mode instability onsets is compatible with the simulated regions of instability obtained from the parametric variations.

Linear stability analysis↗
DOE OSTI2024

Thermochemical Modeling of Radionuclide Vapor-Liquid Equilibria in Sodium Pools for SFR Mechanistic Source Term Analysis

Mechanistic source term (MST) analysis of sodium fast reactors (SFR) requires understanding of various radionuclide (RN) transport phenomena influencing potential releases from the fuel to the environment. One such phenomenon includes the retention or release of RNs from the sodium coolant pool, representing the step after possible fuel failures and influencing transport to the cover gas region. Thermodynamic vapor-liquid equilibria (VLE) calculations were performed on systems representing SFR sodium pools containing oxygen impurities and radionuclide (RN) inventories. First, an assessment and recreation of a previously developed thermodynamic database was completed, including updates to thermodynamic parameters. The RN inventories used in VLE calculations represented hypothetical source terms that might be released to the pool during previously analyzed fuel failure scenarios. The calculations were performed for all possible combinations of sodium pool size (i.e., total oxygen) and number of failed fuel pins (i.e., total RNs). In this way, multiple ratios of the RN relative to the oxygen impurity (RN:O) were compared for their impact on RN volatility, which is discussed in terms of the vapor fraction (VF), defined as the fraction of the RN that is calculated to exist in the vapor phase at equilibrium above condensed phases of that element. Similar trends in VLE behavior are seen in the equilibrium calculation results for elements of similar chemistry, and for some element types, it was found that the RN:O ratio can be important due to oxide formation, which typically exist in the condensed phase.

Shahbazi, Shayan↗
DOE OSTI2024

Variable-spectrum mode control of high poloidal beta discharges

DIII-D experiments demonstrate that high pressure, broad current profile equilibria can be accessed in the high poloidal beta regime by optimizing the MHD mode control poloidal spectrum. A novel, variable spectrum (VS) magnetic feedback scheme implemented using the DIII-D internal non-axisymmetric coils (I-coils) facilitated access to reduced internal inductance $l$ i operation above the no-wall beta limit compared with both no feedback and fixed spectrum feedback. In addition, the VS feedback helped avoid beta collapses caused by marginally unstable resistive wall mode activity. The lower and upper I-coil rows were configured in two independent feedback loops, allowing the feedback field's poloidal spectrum to vary and track changes in the plasma mode structure as the edge safety factor q 95 varied from 11 to 6 during the discharges. The q 95 dependence of the measured phase difference between the lower and upper I-coil rows during VS feedback is qualitatively compatible with ideal MHD simulations of the least-stable plasma kink mode and with plasma response simulations that included kinetic modifications to ideal MHD. The VS feedback approach is a straightforward way to improve resilience to variations in mode structure that occur as plasma parameters change. The demonstrated expansion of the operating space to lower $l$ i is expected to improve the coupling of the plasma kink mode to external fields and beneficial wall eddy currents, and is compatible with high bootstrap fraction operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Corrosion of 316 SS in direct contact with Zircaloy-4 in repository relevant conditions

The effects of interfacial contact between 316 stainless steel (316 SS) waste canister and spent fuel cladding material (Zircaloy-4) in case of water intrusion into a deep geological repository are investigated. Hydrothermal corrosion experiments of 316 SS and Zircaloy-4 in direct contact were performed at elevated temperatures in MilliQ water and simulated Mont Terri groundwater. Hematite, magnetite, chromite, and baddeleyite form as secondary phases. pH of the Mont Terri solution increases with corrosion duration and magnetite forms preferentially instead of hematite. Thermodynamic modeling shows that iron oxidation in 316 SS controls the pH and redox potential of the aqueous environment and is in good agreement with experimental Mont Terri results, showing that experimental Mont Terri results trend to equilibrium. All oxides except hematite are predicted to exist at thermodynamic equilibrium based on calculations, implying that experimentally formed hematite will be reduced or transformed to magnetite at equilibrium. The thermodynamically predicted chromite was observed to form only at the interface where 316 SS was in contact with Zircaloy-4 at 250°C, showing that contact with Zircaloy-4 seems to accelerate the 316 SS corrosion to thermodynamic equilibrium. Finally, Mont Terri groundwater readily trends towards thermodynamic equilibrium, leading to a pH and redox buffered system that is dominated by iron in solution coming from the accelerated corrosion of 316 SS caused by Zircaloy-4 contact.

316 stainless steel↗