Search NASA⌕ Search

SEARCH · Search NASA

Results for “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.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Comparing π‐complexation capabilities of ionic liquids containing silver(I) and copper(I) ions by headspace single drop microextraction in combination with high‐performance liquid chromatography

Selective π‐complexation capabilities of silver(I) and copper(I) ions can be effectively facilitated in ionic liquids. To understand the effects of environmental factors that influence the π‐complexation of these metal ions with analytes, techniques that employ small volumes of ionic liquid that can be readily analyzed are desired. In this study, headspace single drop microextraction coupled with HPLC is used to investigate a diverse set of environmental factors on the metal ion‐mediated complexation with aromatic compounds in ionic liquid media. Silver(I) and copper(I) bis[(trifluoromethyl)sulfonyl]imide salts were both studied by dissolving them in the 1‐decyl‐3‐methylimidazolium bis[(trifluoromethyl)sulfonyl]imide ionic liquid and employing the mixture as extraction media for aromatic compounds. Water and acetonitrile within the sample solution were observed to interfere with the complexation of silver(I) ions and aromatic compounds, while ethylene glycol and triethylene glycol did not. The temperature and extraction times were optimized to fully facilitate the π‐complexation capabilities of metal ions in ionic liquid media. Partition coefficients between the sample headspace and metal ion were determined using a three‐phase equilibria model. Although no discernable difference in analyte partitioning between the headspace and ionic liquid solvent was observed, analyte partition coefficients to silver(I) ion tended to be greater compared to copper(I) ion.

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

Submodular optimization problems and greedy strategies: A survey

The greedy strategy is an approximation algorithm to solve optimization problems arising in decision making with multiple actions. How effective is the greedy strategy compared to the optimal solution? In this survey, we mainly consider two classes of optimization problems where the objective function is submodular. The first is set submodular optimization, which is to choose a set of actions to optimize a set submodular objective function, and the second is string submodular optimization, which is to choose an ordered set of actions to optimize a string submodular function. Our emphasis here is on performance bounds for the greedy strategy in submodular optimization problems. Specifically, we review performance bounds for the greedy strategy, more general and improved bounds in terms of curvature, performance bounds for the batched greedy strategy, and performance bounds for Nash equilibria.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A Regularized Variance-Reduced Modified Extragradient Method for Stochastic Hierarchical Games

We consider an N -player hierarchical game in which the i th player’s objective comprises of an expectation-valued term, parametrized by rival decisions, and a hierarchical term. Such a framework allows for capturing a broad range of stochastic hierarchical optimization problems, Stackelberg equilibrium problems, and leader-follower games. We develop an iteratively regularized and smoothed variance-reduced modified extragradient framework for iteratively approaching hierarchical equilibria in a stochastic setting. We equip our analysis with rate statements, complexity guarantees, and almost-sure convergence results. We then extend these statements to settings where the lower-level problem is solved inexactly and provide the corresponding rate and complexity statements. Our model framework encompasses many game theoretic equilibrium problems studied in the context of power markets. We present a realistic application to the study of virtual power plants, emphasizing the role of hierarchical decision making and regularization. Preliminary numerics suggest that empirical behavior compares well with theoretical guarantees.

Tikhonov regularization↗
DOE OSTI2020

Experimental and Modeling Review of the Plutonium-Zirconium (Pu-Zr) System: Lost in Translation and Over Time?

This article presents a review of the phase equilibria of the Pu-Zr system and discusses the contradictory reports on the experimental phase diagram and misinterpretations that have led to confusion over time. In addition, a review of the few Pu-Zr CALPHAD (CALculation of PHAse Diagrams) assessments is presented with emphasis on the heat of formation of the bcc phase, highlighting a disagreement between CALPHAD and ab initio calculations. Based on the information gathered in this study, a critical re-assessment of the Pu-Zr system involving the ζ (Pu 28 Zr), θ (Pu 4 Zr), kinetically hindered κ (PuZr 2 ), and once observed ι phases is called for.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermodynamic Modeling of the Al-Ce-Cu-Mg-Si System and Its Application to Aluminum-Cerium Alloy Design

Recently discovered AlCe alloys have shown promise in a number of applications, but the propensity of Ce to react with Al and other alloying elements can complicate the phase equilibria and design approach. To solve this, the CALPHAD method is used to explore an alloy within the quinary Al-Ce-Cu-Mg-Si system by developing a thermodynamic database with self-consistent parameters. The database includes a description of all 10 binary systems and 8 ternary systems consisting of: (i) 6 Al-containing ternaries (Al-Ce-Cu, Al-Ce-Mg, Al-Ce-Si, Al-Cu-Mg, Al-Cu-Si and Al-Mg-Si); and (ii) 2 additional ternaries that include Mg and Si (i.e., Ce-Mg-Si and Cu-Mg-Si). The thermodynamic description for the Al-Ce-Mg and Al-Mg-Si systems were reassessed to ensure consistency with the binary systems and the Ce-Mg-Si system is presented for the first time and compared to theoretical data from DFT (Density Functional Theory). In addition to the ternary interactions, the quaternary compound Al3Cu2Mg9Si7 and solid solution extending from the ternary Al2CuMg phase (Al,Si)2CuMg are incorporated. The CALPHAD method is employed and leveraged through the use of a Materials Design Simulator (MDS) to accelerate the design of novel aluminum-cerium-based alloys. The combination of a CALPHAD-based framework with experimental efforts and industrial insight permits the development of three new Al-Ce alloys: Al-3.5Ce-0.4Mg-7Si (Ce-modified A356), Al-5Ce-1Cu-0.5Mg-10Si and Al-19Ce-0.9Mg-1.1Si.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Thermodynamics and Kinetics of High-Entropy Alloys

Here, this editorial introduces a special issue of the Journal of Phase Equilibria and Diffusion that is comprised of 19 high-quality papers that report the state-of-the-art frontier research in various high entropy alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Interactive Exploration of High-Dimensional Phase Diagrams

High-dimensional thermodynamic phase stability databases are becoming increasingly common due to the convergence of three recent trends: (i) the widespread interest in so-called “high-entropy” alloys, (ii) the availability of high-throughput computational assessments of phase stability in broad composition spaces and (iii) the ongoing development of ever-increasingly broad, multicomponent, multiphase CALPHAD databases. Although automated computational tools can readily process such high-dimensional data, scientists are often unable to visualize the relevant phase relations, an ability that is crucial to gaining an intuitive understanding of the stability constraints governing materials design. The present work addresses this need by providing algorithms that enable the interactive exploration of phase equilibria in high-dimensional spaces. These algorithms concentrate the complex nonlinear nonsmooth optimization needed into a preprocessing step that generates a large number of high-dimensional yet elementary graphical primitives. Furthermore, these primitives can then be cross-sectioned to yield 3-dimensional views in a computationally efficient manner that enables an interactive exploration of high-dimensional spaces. All of these operations are highly parallelizable, thus facilitating scaling of this method to large data sets.

36 MATERIALS SCIENCE↗
DOE OSTI2024

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.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effect of magneto-mechanical synergism in the process-structure correlation in Fe–C alloys: A phase-field modeling approach

Applied magnetic fields can alter phase equilibria and kinetics in steels; however, quantitatively resolving how magnetic, chemical, and elastic driving forces jointly influence the microstructure remains challenging. We develop a quantitative magneto-mechanically coupled phase-field model for the Fe–C system that couples a CALPHAD-based chemical free energy with demagnetization-field magnetostatics and microelasticity. Here, the model reproduces single- and multi-particle evolution during the α → γ inverse transformation at 1023 K under external fields up to 20 T, including ellipsoidal morphologies observed experimentally at 8 T. Chemically driven growth is isotropic; a magnetic interaction introduces an anisotropic driving force that elongates γ precipitates along the field into ellipsoids, while elastic coherency promotes faceting, yielding elongated cuboidal or “brick-like” particles under combined magneto-elastic coupling. Growth kinetics increase with C content, and decrease with field strength and misfit strain. Multi-particle simulations reveal dipolar interaction-mediated coalescence for field-parallel neighbors and ripening for field-perpendicular neighbors. Incorporating field-dependent diffusivity from experiment slows kinetics as expected; a first-principles-motivated anisotropic diffusivity correction is estimated to be small (<2%). These results establish a process-structure link for magnetically assisted heat treatments of Fe–C alloys and provide guidance for microstructure control via chemo-magneto-mechanical synergism.

Magnetic field↗
DOE OSTI2023

Thermodynamic modeling of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems

Molten salt electrolysis can recover the actinides from spent nuclear fuels, and it involves a eutectic LiCl-KCl in molten form as an electrolyte. During reprocessing, the concentration of fission products such as La, Nd, Pr, etc., increases and affects the recovery efficiency of the electrolyte. In this work, thermodynamic modeling of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems was carried out using the CALPHAD (Calculation of Phase Diagrams) approach for the first time. The thermodynamic functions for the pure salts were taken from the SGTE (Scientific Group Thermodata Europe) Substances (SSUB) database. The experimental thermochemical and phase equilibria data available in the literature were used as input for the assessment of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems. The model parameters for the KCl-LiCl system were adjusted to include the new Gibbs energy descriptions for the pure salts. In addition, the sublattice model for the liquid phase in the LiCl-PrCl 3 system was modified and reassessed to ensure the model compatibility for higher-order extrapolation. There is a good agreement between the experimental and calculated thermochemical and phase diagram data for all the optimized constituent binaries and the ternary system. Furthermore, this work will be beneficial for determining the solubility limit of PrCl 3 in molten LiCl-KCl electrolytes and their thermodynamic properties for improving the efficiency of the pyrochemical process.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Thermodynamic modeling of the Nb-Ni system with uncertainty quantification using PyCalphad and ESPEI

Here, the Nb–Ni system is remodeled with uncertainty quantification (UQ) using software tools of PyCalphad and ESPEI (the Extensible, Self-optimizing Phase Equilibria Infrastructure) with the presently implemented capability of modeling site fraction based on Wyckoff positions. The five- and three-sublattice models are used to model the topologically close pack (TCP) μ-Nb 7 Ni 6 and δ-NbNi 3 phases according to their Wyckoff positions. The inputs for CALPHAD-based thermodynamic modeling include the thermochemical data as a function of temperature predicted by first-principles and phonon calculations based on density functional theory (DFT), ab initio molecular dynamics (AIMD) simulations, together with phase equilibrium and site fraction data in the literature. In addition to phase diagram and thermodynamic properties, the CALPHAD-based predictions of site fractions of Nb in μ-Nb 7 Ni 6 agree well with experimental data. Furthermore, the UQ estimation using the Markov Chain Monte Carlo (MCMC) method as implemented in ESPEI is applied to study the uncertainty of site fraction in μ-Nb 7 Ni 6 and enthalpy of mixing (ΔH mix ) in liquid.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Time-dependent water vapor desorption isotherm model of hardened cement paste

Highlights: • Short-term and long-term desorption isotherms of hardened cement paste are reviewed. • Water content at RH depends on both thermodynamic and microstructural equilibria. • An extended GAB equation is proposed to fit the cavitation drop during desorption. • A dynamic time-dependent sorption isotherm model is shown and used for simulations. • Numerical results capture anomalous drying in the same RH as seen experimentally. This study aims to enhance moisture transport modeling by elucidating the so-called anomalous water diffusion in cementitious materials. Water desorption isotherms are measured from samples at various drying stages using techniques with different durations to quantify the impact of the drying duration on the water sorption ability. A single water sorption isotherm does not solely give a relationship between the water content and relative humidity but also represents the state of microstructure. The continuous evolution of a desorption isotherm due to drying-induced microstructural rearrangement is demonstrated. For numerical modeling, the microstructural alteration can be explicitly considered through a dynamic desorption isotherm model, which governs the local thermodynamic equilibrium at the capillary meniscus. This approach is implemented into a multiphase transport code, whose ability to predict drying is validated using literature data. Finally, the effect of prolonged drying on the colloidal nature of the calcium-silicate-hydrate gel is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

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.

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

Major element composition controls rare earth element solubility during leaching of coal fly ash and coal by-products

Coal combustion ash and pre-combustion coal refuse are currently under consideration as potential sources of rare earth elements. One of the early steps in recovering REEs from coal by-products is often acid leaching, which can result in low pH leachates with complex aqueous chemistry. The aim of this work was to understand the connection between REE solubility, pH, and major elemental components of leachates for coal by-products. To accomplish this, we investigated the effects of solids concentration (i.e., pulp density) and pH adjustment on REE solubility in acid leachates of coal fly ashes from the Powder River Basin (PRB) and Appalachian Basin in the United States, and a coal processing refuse from the Southwestern U.S. For PRB ashes, the concentrations of soluble REEs generally increased with increasing pulp density; however, at pulp density values above 80–100 g/L, the soluble REE concentrations in the leachates were markedly lower. Similarly, the soluble concentrations of other major solutes (Fe, Al, Si) that leached from PRB fly ashes were also non-linear with pulp density. These major elements tended to reach maximum concentration values at 60–70 g/L pulp density. In contrast, for the Appalachian fly ashes and the coal by-product, soluble concentrations of REE and major elements in leachates increased linearly with pulp density. Chemical equilibria calculations of mineral saturation indices indicated that trends in soluble REE concentrations could be explained by saturation conditions for Fe and Al-(hydr)oxides and possibly sulfate minerals, but not lanthanide hydroxides. Furthermore, pH adjustment of the acid leachates showed that REEs and many major solutes were removed from solution at pH values above 4.5, also consistent with Fe- and Al-(hydr)oxide precipitation. These results highlight the importance of understanding the chemical composition of leachates when designing REE recovery processes for low-grade geologic feedstocks and that precipitation of hydr(oxide) or sulfate minerals of major elements rather than discreet formation of REE mineral phases could be used for process optimization.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Comparison of Na 2 SO 4 , K 2 SO 4 and Na 2 SO 4 -K 2 SO 4 deposit induced hot corrosion of a β-NiAl coating

Here, hot corrosion behaviour of a β-NiAl coating induced by Na 2 SO 4 , K 2 SO 4 , and Na 2 SO 4 -K 2 SO 4 deposits was investigated at 700 °C in air + 300 ppm SO 2 . To provide an insight into the corrosion mechanisms, phase equilibria of the Na 2 SO 4 -K 2 SO 4 -NiO-SO 3 system were calculated using an in-house developed thermodynamic database. The attack rates and morphologies could be correlated with the deposit compositions and their reactions with oxides. K 2 SO 4 and Na 2 SO 4 -K 2 SO 4 deposits are more corrosive than Na 2 SO 4 . This is attributed to a reduction of the melting temperature of the deposit and the formation of Ni-containing mixed sulphate phases in the presence of K 2 SO 4 .

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

Adaptive Hybridizable Discontinuous Galerkin discretization of the Grad–Shafranov equation by extension from polygonal subdomains

In this work, we propose a high-order adaptive numerical solver for the semilinear elliptic boundary value problem modeling magnetic plasma equilibrium in axisymmetric confinement devices. In the fixed boundary case, the equation is posed on curved domains with piecewise smooth curved boundaries that may present corners. The solution method we present is based on the hybridizable discontinuous Galerkin method and sidesteps the need for geometry-conforming triangulations thanks to a transfer technique that allows to approximate the solution using only a polygonal subset as computational domain. Moreover, the solver features automatic mesh refinement driven by a residual-based a posteriori error estimator. As the mesh is locally refined, the computational domain is automatically updated in order to always maintain the distance between the actual boundary and the computational boundary of the order of the local mesh diameter. Numerical evidence is presented of the suitability of the estimator as an approximate error measure for physically relevant equilibria with pressure pedestals, internal transport barriers, and current holes on realistic geometries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Solving the Grad–Shafranov equation using spectral elements for tokamak equilibrium with toroidal rotation

The Grad–Shafranov equation is solved using spectral elements for tokamak equilibrium with toroidal rotation. The Grad–Shafranov solver builds upon and extends the NIMEQ code (Howell and Sovinec, 2014) previously developed for static tokamak equilibria. Both geometric and algebraic convergence are achieved as the polynomial degree of the spectral-element basis increases. A new analytical solution to the Grad–Shafranov equation is obtained for Solov’ev equilibrium in presence of rigid toroidal rotation, in addition to a previously obtained analytical solution for a different set of equilibrium and rotation profiles. The numerical solutions from the extended NIMEQ are benchmarked with the analytical solutions, with good agreements. Besides, the extended NIMEQ code is benchmarked with the FLOW code (Guazzotto et al., 2004). The modification of pressure profile induced by toroidal flow is investigated. The relative change of pressure profile is found significant around the edge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Hermes-3: Multi-component plasma simulations with BOUT++

A new open source tool for fluid simulation of multi-component plasmas is presented, based on a flexible software design that is applicable to scientific simulations in a wide range of fields. Hermes-3 is built on plasma simulation framework BOUT++, consolidating earlier SD1D and Hermes models into a single code that can be configured at run-time to solve plasma models in 1D, 2D or 3D, either for transport (steady-state) or turbulent (time-evolving) problems, with an arbitrary number of ion and neutral species. Here, we describe the improved numerical algorithms and software design that have been implemented in Hermes-3. To demonstrate the capabilities of this tool, applications relevant to the boundary of tokamak plasmas are presented: 1D simulations of diveror plasmas evolving equations for all charge states of neon and deuterium; 2D transport simulations of tokamak equilibria in single-null X-point geometry with plasma ion and neutral atom species; and simulations of the time-dependent propagation of plasma filaments (blobs). Hermes-3 is publicly available on Github under the GPL-3 open source license. The repository includes documentation and a suite of unit, integrated and convergence tests.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗