SEARCH · Engineering Papers
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.
Melting temperature, critical nucleus size, and interfacial free energy in single FCC metals — A Molecular Dynamics study of liquid–solid phase equilibria
Not Available
Deep Planet Chemistry from Atomistic Simulation: Phase Equilibria Without Mapping Free Energy
Explore the source record for details and available documents.
On Gibbs Equilibrium and Hillert Nonequilibrium Thermodynamics
During his time at Royal Institute of Technology (Kungliga Tekniska högskolan) in Sweden, the present author learned nonequilibrium thermodynamics from Mats Hillert. The key concepts are the separation of internal and external variables of a system and the definitions of potentials and molar quantities. In equilibrium thermodynamics derived by Gibbs, the internal variables are not independent and can be fully evaluated from given external variables. While irreversible thermodynamics led by Onsager focuses on internal variables though often mixed with external variables. Hillert integrated them together by first emphasizing their differences and then examining their connections. His philosophy was reflected by the title of his book “Phase Equilibria, Phase Diagrams and Phase Transformations” that puts equilibrium, nonequilibrium, and internal processes on equal footing. Here, in the present paper honoring Hillert, the present author reflects his experiences with Hillert and his work in last 40 years and expresses his gratitude for all the wisdom and support from him in terms of “Hillert nonequilibrium thermodynamics” and discusses some recent topics that the present author has been working on.
Large-Signal Stability of Phase-Balanced Equilibria in Single-Phase Grid-Forming Inverter Systems
This article explores the setup where large numbers of single-phase grid-forming inverters with droop control across distribution networks self-organize into a stable and balanced system with 120° phase offsets across aggregates in the absence of balanced three-phase generating resources or external communication. A suite of circuit- and system-theoretic notions are leveraged to derive a dynamical model for phase-angle differences across aggregates of inverters connected in the three phases. Focusing on this model, large-signal stability is established and the region of attraction of the phase-balanced equilibria is determined with the aid of a Lyapunov function. Experimental validation for a bench-top prototype network is included to support the analytical developments. Altogether, the effort supports the vision of facilitating balanced operation of distribution networks with grid-forming inverters during service disruptions at the bulk transmission network.
Equilipy
There is high demand for high-throughput calculations of phase equilibria based on the CALPHAD method (CALculation of PHAse Diagram). High-throughput calculations are possible through HPC, currently by a commercial program called Thermo-Calc. However, the number of nodes/processors are limited to the number of purchased license (16 processors per license). This program provides a toolkit for high-throughput calculations of phase equilibria by the CALPHAD method. The program is prepared for Python environment, so that it can be easily installed and used together with other open-source programs. The program can be run in supercomputer.
Development of a Thermodynamic Database for Corrosion in Chloride MSRs
Control of the salt redox condition is the salient method for corrosion inhibition in molten salts. For chloride salt-fueled molten salt reactors of primary concern is the U(IV)/U(III) ratio, influenced directly by the Cl2 activity in the salt. Although thermodynamic assessments of MSR-relevant salts abound, none have included consideration of UCl4 which is essential to correctly describe the salt redox condition. This work reports the assessment of the NaCl-KCl-MgCl2 base salt with UCl3-UCl4 and corrosion chlorides involving Cr, Fe, and Ni as minor components and the resultant thermodynamic models. The effort utilized available thermodynamic measurements and phase equilibria complemented by differential scanning calorimetry determinations of phase equilibria. Example calculations explore the temperature and composition space to understand its influence on corrosion. The foregoing models are incorporated into an update of the Molten Salt Thermal Properties Database – Thermochemical which is publicly available.
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.
Predicting Multicomponent Aqueous Phase Adsorption Equilibria of Organic Acids Using the Generalized Brunauer–Emmett–Teller Isotherm Model
Here, to support process development of adsorptive separation of organic acids, this work presents a thermodynamic modeling methodology to predict multicomponent aqueous phase organic acid adsorption equilibria from single-component adsorption isotherms using the generalized Brunauer–Emmett–Teller isotherm model (gBET). With the organic acid fugacities rigorously accounting for the aqueous phase organic acid solution chemistry and solution nonideality, gBET precisely represents the single-component adsorption isotherms and accurately predicts the binary and ternary mixed-acid adsorption equilibria for the ranges of initial pH (∼3–7), acid concentration (100–400 mmol/L), and temperature (298.15–328.15 K) with less than 10% average absolute relative deviation. In addition, gBET with pH-independent parameters provides insights into the underlying adsorption phenomena, including the adsorbate loadings and compositions in the monolayer and subsequent layers under varying initial pH, temperature, concentration, and composition. The gBET model predictions outperform the predictions from the overloading model and the Ideal Dilute Solution Theory.
The influence of cooling rates on strain phase diagrams and domain structures of ferroelectric thin films: A case study of PbTiO 3
Strain engineering has been established as an effective approach to control phase equilibria, domain configurations, and functional properties of ferroelectric thin films. Temperature-strain phase diagrams have been used as powerful tools for providing insights into strain engineering. However, almost all existing phase diagrams established using the phase-field approach assume quenching conditions without considering actual cooling rates during the post-deposition annealing process of ferroelectric thin films. Within this work, we systematically investigate the influence of cooling rates on domain structures and the strain-phase diagram of ferroelectric thin films using phase-field simulations, taking PbTiO 3 thin films as a model system. We found that both the position of phase boundaries in the strain phase diagrams and the domain morphology are significantly influenced by the cooling rates. It is revealed that while the paraelectric-ferroelectric phase boundary remains invariant, the phase boundaries between single-phase and multi-phase regions tend to shift toward the corresponding multi-phase region as the cool rate reduces. Slow cooling generally leads to more ordered domain structures with increased domain size. Using the obtained equilibrium domain structures, we calculated effective thermal conductivities and found significant variations that can be tuned by the cooling rates. In conclusion, this work reveals an underexplored yet critical impact of cooling rates on phase equilibria and domain structures in ferroelectric thin films, which may inspire further fine-tuning of domains and domain walls in low-dimensional ferroelectrics for multifunctional applications.
IDAES-PSE 2.5.0 Release
The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost, most environmentally sustainable solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.5.0 Release Highlights Upcoming Changes IDAES will be switching to the new Pyomo solver interface in the next release. Whilst this will hopefully be a smooth transition for most users, there are a few important changes to be aware of. The new solver interface uses a different version of the IPOPT writer (“ipopt_v2”) and thus any custom configuration options you might have set for IPOPT will not carry over and will need to be reset. By default, the new Pyomo linear presolver will be activated with ipopt_v2. Whilst are working to identify any bugs in the presolver, it is possible that some edge cases will remain. IDAES will begin deploying a new set of scaling tools and APIs over the next few releases that make use of the new solver writers. The old scaling tools and APIs will remain for backward compatibility but will begin to be deprecated. New Models, Tools and Features New diagnostics check for near-parallel variables and constraints. New diagnostics tools for identifying causes of infeasibility in models. New example for creating a custom model of a liquid-liquid extractor unit operation. Bug Fixes Fixed bug in Gibbs reactor that caused it to appear to have additional spurious degrees of freedom. Fixed bug in the Modular Property Framework that would cause errors when trying to use phase-based material balances with phase equilibria. Fixed bug in Modular Properties Framework that caused errors when initializing models with non-vapor-liquid phase equilibria. Testing and Robustness Deployed the IDAES Diagnostics Toolbox to confirm that there are no structural or numerical issues in the core model libraries. Additional robustness tests for core model, and some associated improvements in the converge tester class. Fixed a number of issues that were causing unexpected warnings to be emitted during testing. Deprecations and Removals Removed examples for RIPE tool which has not been supported for a number of releases.
IDAES-PSE 2.6.0 Release
The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.6.0 Release Highlights Upcoming Changes IDAES will be switching to the new Pyomo solver interface in the next release. Whilst this will hopefully be a smooth transition for most users, there are a few important changes to be aware of. The new solver interface uses a different version of the IPOPT writer (“ipopt_v2”) and thus any custom configuration options you might have set for IPOPT will not carry over and will need to be reset. By default, the new Pyomo linear presolver will be activated with ipopt_v2. Whilst are working to identify any bugs in the presolver, it is possible that some edge cases will remain. IDAES will begin deploying a new set of scaling tools and APIs over the next few releases that make use of the new solver writers. The old scaling tools and APIs will remain for backward compatibility but will begin to be deprecated. New Models, Tools and Features New Intersphinx extension automatically linking Jupyter notebook examples to project documentation New end-to-end diagnostics example demonstrated on a real problem New complementarity formulation for VLE with cubic equations of state, backward compatibility for old formulation New solver interface with presolve (ipopt_v2) in support of upcoming changes to the initialization and APIs methods, with default set to ipopt to maintain backwards compatibility; this will deprecate once all examples have been updated New forecaster and parameterized bidder methods within grid integration library Updated surrogates API and examples to support Keras 3, with backwards compatibility for older formats such as TensorFlow SavedModel (TFSM) Updated costing base dictionary to include the 2023 cost year index value Updated ProcessBlock to include information on the constructing block class Updated Flowsheet Visualizer to allow visualize() method to return value and functions Bug Fixes Fixed bug in the Modular Property Framework that would cause errors when trying to use phase-based material balances with phase equilibria. Fixed bug in Modular Properties Framework that caused errors when initializing models with non-vapor-liquid phase equilibria. Fixed typos flagged by June update to crate-ci/typos and removed DMF-related exceptions Minor corrections of units of measurement handling in power plant waste/transport costing expressions, control volume material holdup expressions, and BTX property package parameters Fixed throwing >7500 numpy deprecation warnings by replacing scalar value assignment with element extraction and item iteration calls Testing and Robustness Migrated slow tests (>10s) to integration, impacting test coverage but also yielding a nearly 30% decrease in local test runtime Pinned pint to avoid issues with older supported Python versions Pinned codecov versions to avoid tokenless upload behavior with latest version Bumped extensions to version 3.4.2 to allow pointing to non-standard install location Deprecations and Removals Python 3.8 is no longer supported. The supported Python versions are 3.9 through 3.12 The Data Management Framework (DMF) is no longer supported. Importing idaes.core.dmf will cause a deprecation warning to be displayed until the next release The SOFC Keras surrogates have been removed. The current version of the SOFC surrogate model in the examples repository is a PySMO Kriging model.
Thermodynamic evaluation and modeling of the CeF 3 -LaF 3 -LiF-NaF system
The thermodynamic properties of the LiF-NaF-CeF 3 -LaF 3 pseudo-quaternary system have been evaluated, and models for the relevant phases have been revised or newly developed. The solid solubility of NaF in LaF 3 and CeF 3 was described in the model based on reported phase-equilibria observations and results from density functional theory calculations. The Gibbs energy functions for the NaCeF 4 and NaLaF 4 phases have been updated with recently reported thermochemical measurements. The Gibbs energy functions and model relations for the LiF-NaF-LaF 3 and LiF-NaF-CeF 3 pseudo-ternary systems were successfully optimized based on thermochemical properties and phase equilibria observations.
Thermodynamic Modeling of Complex Solid Solutions in the Lu-H-N System via Graph Neural Network Accelerated Monte Carlo Simulations
Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., 𝑂(10 2 ) GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the 𝑂(10 0 ) GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va) 3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.
Thermodynamic properties of ZrCl 4 with LiCl, NaCl, KCl, CsCl, MgCl 2 , and UCl 3 for molten salt reactor applications
Here, a set of self-consistent Gibbs energy functions has been developed to describe the thermochemical behavior of the major reactive fission product ZrCl 4 with key components in chloride molten salt reactors (MSRs): LiCl–ZrCl 4 , NaCl–ZrCl 4 , KCl–ZrCl 4 , CsCl–ZrCl 4 , MgCl 2 –ZrCl 4 and UCl 3 –ZrCl 4 . Low ZrCl 4 concentration phase equilibria most relevant to MSR applications have been confirmed via differential scanning calorimetry for NaCl–ZrCl 4 , and X-ray diffraction analysis performed on equilibrated samples of NaCl–ZrCl 4 , KCl–ZrCl 4 , MgCl 2 –ZrCl 4 , and UCl 3 –ZrCl 4 . Within the framework of the modified quasi-chemical model in the quadruplet approximation, extrapolations of pseudo-binary models were generated to represent KCl–MgCl 2 –ZrCl 4 , KCl–NaCl–ZrCl 4 , and MgCl 2 –NaCl–ZrCl 4 , which show agreement with phase equilibria data where available without the use of ternary interaction parameters. The optimized thermodynamic descriptions for these systems and others are available in the open-source compendium Molten Salt Thermal Properties Database – Thermochemical (MSTDB–TC).
Probing ethane phase changes in bead pack via high-field NMR spectroscopy
Phase equilibria analysis in two-phase systems requires discernment of the individual phases as well as the interfacial region. The phase behavior of fluids in constraint geometries have been studied through several techniques. High-Field Nuclear Magnetic Resonance (HF-NMR) can be used to detect fluid phases, i.e. gas, liquid, or supercritical fluid, in the bulk or under confinement. In solution-state HF-NMR, each phase exhibits a distinct chemical shift that correlates to the phase density. In turn, the emergence of multimodal in NMR spectra reveals the coexistence of multiple phases in fluid systems. Our results for ethane show that HF-NMR can enable tracking phase changes in the bulk, but also in macro-scale (10–100's μm in pore size) porous systems at various pressure and temperature conditions, even when the system us out of equilibrium. Finally, these results creates potential opportunities to investigate phase changes influenced by confinement in porous systems, though additional steps are necessary.
Phase diagram to demarcate supercritical, transcritical, and continuous phase regimes for binary fluid equilibrium mixing relevant to combustion applications
Here, a robust methodology to develop phase diagrams of binary fluid mixtures at fixed thermo dynamic conditions (pressure, temperature, and mole fraction) as well as of two initially separated fluids undergoing mixing near critical conditions are presented for fluids and con ditions relevant to rockets, gas turbines, and diesel engine applications. Phase equilibria of mixtures is first examined to provide insight into the continuous-phase mixing behavior (including but not limited to supercritical behavior), and to develop a broadly applicable phase-diagram for binary fluid mixtures at fixed conditions. Next, adiabatic mixing theory and reduced Helmholtz equations of state are used to predict the thermodynamic conditions required to attain continuous-phase binary fluid mixing near critical conditions. Then, a 3D surface diagram (P,T fuel ,T amb ) separating single and two-phase regions is constructed by varying the ambient pressure and the initial temperatures of the two fluids. Polynomial fits of the 3D surfaces for 10 different binary mixtures are tabulated for nitrogen-alkane and methane-oxygen blends relevant to air-breathing and propellant based engines, respectively.
Theory-guided design of duplex-phase multi-principal-element alloys
Density-functional theory (DFT) is used to identify phase-equilibria in multi-principal-element and high-entropy alloys (MPEAs/HEAs), including duplex-phase and eutectic microstructures. Here, a combination of composition-dependent formation energy and electronic-structure-based ordering parameters were used to identify a transition from FCC to BCC favoring mixtures, and these predictions experimentally validated in the Al-Co-Cr-Cu-Fe-Ni system. A sharp crossover in lattice structure and dual-phase stability as a function of composition were predicted via DFT and validated experimentally. The impact of solidification kinetics and thermodynamic stability was explored experimentally using a range of techniques, from slow (castings) to rapid (laser remelting), which showed a decoupling of phase fraction from thermal history, i.e., phase fraction was found to be solidification rate-independent, enabling tuning of a multi-modal cell and grain size ranging from nanoscale through macroscale. Strength and ductility tradeoffs for select processing parameters were investigated via uniaxial tension and small-punch testing on specimens manufactured via powder-based additive manufacturing (directed-energy deposition). This work establishes a pathway for design and optimization of next-generation multiphase superalloys via tailoring of structural and chemical ordering in concentrated solid solutions.