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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS2013

The W-W02 Oxygen Fugacity Buffer at High Pressures and Temperatures: Implications for f02 Buffering and Metal-silicate Partitioning

Oxygen fugacity (fO2) controls multivalent phase equilibria and partitioning of redox-sensitive elements, and it is important to understand this thermodynamic parameter in experimental and natural systems. The coexistence of a metal and its oxide at equilibrium constitutes an oxygen buffer which can be used to control or calculate fO2 in high pressure experiments. Application of 1-bar buffers to high pressure conditions can lead to inaccuracies in fO2 calculations because of unconstrained pressure dependencies. Extending fO2 buffers to pressures and temperatures corresponding to the Earth's deep interior requires precise determinations of the difference in volume (Delta) V) between the buffer phases. Synchrotron x-ray diffraction data were obtained using diamond anvil cells (DAC) and a multi anvil press (MAP) to measure unit cell volumes of W and WO2 at pressures and temperatures up to 70 GPa and 2300 K. These data were fitted to Birch-Murnaghan 3rd-order thermal equations of state using a thermal pressure approach; parameters for W are KT = 306 GPa, KT' = 4.06, and αKT = 0.00417 GPa K-1. Two structural phase transitions were observed for WO2 at 4 and 32 GPa with structures in P21/c, Pnma and C2/c space groups. Equations of state were fitted for these phases over their respective pressure ranges yielding the parameters KT = 190, 213, 300 GPa, KT' = 4.24, 5.17, 4 (fixed), and αKT = 0.00506, 0.00419, 0.00467 GPa K-1 for the P21/c, Pnma and C2/c phases, respectively. The W-WO2 buffer (WWO) was extended to high pressure by inverting the W and WO2 equations of state to obtain phase volumes at discrete pressures (1-bar to 100 GPa, 1 GPa increments) along isotherms (300 to 3000K, 100 K increments). The slope of the absolute fO2 of the WWO buffer is positive with increasing temperature up to approximately 70 GPa and is negative above this pressure. The slope is positive along isotherms from 1000 to 3000K with increasing pressure up to at least 100 GPa. The WWO buffer is at a higher fO2 than the IW buffer at pressures lower than 40 GPa, and the magnitude of this difference decreases at higher pressures. This qualitatively indicates an increasingly lithophile character for W at higher pressures. The WWO buffer was quantitatively applied to W metal-silicate partitioning by using the WWO-IW buffer difference in combination with literature data on W metal-silicate partitioning to model the exchange coefficient (KD) for the Fe-W exchange reaction. This approach captures the pressure dependence of W metal-silicate partitioning using the WWO-IW buffer difference and models the activities of the components in the silicate and metallic phases using an expression of the Gibbs excess energy of mixing. Calculation of KD along a peridotite liquidus predicts a decrease in W siderophility at higher pressures that supports the qualitative behavior predicted by the WWO-IW buffer difference, and agrees with findings of others. Comparing the competing effects of temperature and pressure on W metal-silicate partitioning, our results indicate that pressure exerts a greater effect.

Shofner, G. A.↗

NASA NTRS1975

Origin of titaniferous lunar basalts

Delineation of low pressure phase equilibria in the composition space relevant to titaniferous lunar basalts demonstrates a significant degree of control by those equilibria on the compositions of the basalts. The existence of two distinct chemical groups of basalts (high and low K) which cannot be related one to the other by fractional crystallization at any pressure, suggests that melting is responsible for the two groups. Consideration of the pressure shift required to produce the differences between groups constrains magma segregation to have occurred in the outer 150 km of the moon. It is difficult to relate low-Ti and high-Ti basalts to the same source region. The preferred source region of high-Ti basalts, based on phase equilibrium considerations, is a late ilmenite-rich cumulate produced from the residual liquid of the primordial differentiation of the outer portions of the moon. This ilmenite-rich layer is sandwiched between the lunar feldspathic crust and a complementary mafic cumulate.

Walker, D.↗

DOE OSTI2023

Mixture Model for Refrigerant Pairs R-32/1234yf, R-32/1234ze(E), R-1234ze(E)/227ea, R-1234yf/152a, and R-125/1234yf

In this work, thermodynamic models based on the corresponding states framework with departure terms are developed for the refrigerant pairs R-32/1234yf, R-32/1234ze(E), R-1234ze(E)/227ea, R-1234yf/152a, and R-125/1234yf. These models are based on new measurements of density, speed of sound, and phase equilibria, combined with the data available in the literature. The model for R-32/1234yf is most comprehensive in its data coverage, with speed of sound deviations within 1%, density deviations within 0.1%, and bubble- and dew-point pressure deviations within 1%. In conclusion, the other mixtures have generally more limited data availability but a similar goodness of fit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Thermodynamic modeling of CsF with LiF-NaF-KF for molten fluoride-fueled reactors

Gibbs energy models were developed to describe the thermochemical behavior of CsF in molten FLiNaK (46.5LiF-11.5NaF-42KF mol%), a proposed molten salt reactor (MSR) fuel solvent and coolant, as cesium is of concern due to its high radiotoxicity and volatility. Initially, it was necessary to obtain a more accurate Gibbs energy function for CsF which required fitting parameters to reported vapor pressures over condensed phase CsF. The pseudo-binary systems CsF-LiF, CsF-NaF and CsF-KF were then evaluated utilizing phase equilibria and enthalpy of mixing (Δ mix H) values, together with original differential scanning calorimetry (DSC) measurements performed for the CsF-LiF and CsF-KF systems. The CsF-LiF-NaF, CsF-LiF-KF and CsF-NaF-KF pseudo-ternary system representations were obtained by interpolation of the constituent pseudo-binary systems, with DSC measurements performed for the CsF-LiF-NaF system to corroborate the calculated liquidus temperature. Ultimately, the pseudo-ternary systems were interpolated to obtain Gibbs energy models for the pseudo-quaternary CsF-LiF-NaF-KF system, supported by DSC measurements at low CsF compositions (1–10 mol%), yielding computed equilibria and cesium-containing vapor pressures that compare favorably with reported values. In conclusion, the Molten Salt Thermal Properties Database – Thermochemical (MSTDB-TC) was subsequently expanded to include these Gibbs energy models allowing description of the thermochemical behavior of the CsF-LiF-NaF-KF system.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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 OSTI2025

Strength, deformation, and the fcc–hcp phase transition in condensed Kr and Xe to the 100 GPa pressure range

The rare gas solids exhibit systematic differences in crystal structure, phase transition conditions, bond strength, and other physical properties. The physical properties of heavy rare gas solids krypton and xenon are modified by the martensitic phase transition from face-centered cubic to hexagonal close packed structure over a broad pressure range. Crystal structure, strength, and plastic deformation of krypton and xenon have been investigated at 300 K using compression in the diamond-anvil cell with synchrotron angle-dispersive x-ray diffraction and complementary ruby fluorescence spectroscopy for Xe. Stacking faults indicative of the fcc–hcp phase transition are observed at pressures at and above 1.23 ± 0.05 and 1.9 ± 0.6 GPa in Kr and Xe, respectively. The transition remains incomplete in both solids to pressures greater than 100 GPa. Strength determined from stress measurements in Pt and ruby standards at pressures up to 111 GPa and complemented by observations of strain and texture measurements obtained by x-ray diffraction in the radial geometry to 100 GPa indicates similar or higher strength than Ar at all conditions, with significant stiffening at 15–20 GPa. Radial diffraction data reveal the persistence of broad highly textured fcc diffraction lines to 101 GPa in Xe, suggesting that the axial measurements may underestimate the metastable persistence of the fcc phase due to biased sampling of hcp crystallites resulting from preferred crystallite orientation. Kr and Xe are compared with He, Ne, and Ar for a systematic understanding of physical properties and phase equilibria of rare gas solids.

Compressive stress↗

DOE OSTI2024

Computational thermodynamics-guided alloy design and phase stability in CoCrFeMnNi-based medium-and high-entropy alloys: An experimental-theoretical study

A computational thermodynamics approach has been employed to design CoCrFeMnNi-based medium- and high-entropy alloys (M/HEAs) with systematically varied compositions (Co(( 80-X )/ 2 )Cr(( 80-X )/ 2 )Fe X Mn 10 Ni 10 with x = 30, 40, and 50 at.%) and phase stability. Since the formation of sigma phase, usually brittle and undesirable, is a common concern, when this class of alloys is subjected to elevated temperatures (600–1000 °C), predicting its formation becomes essential. Thus, its formation and the phase equilibria were studied using the CALPHAD method, and two empirical methods, namely, valence electron concentration (VEC) and paired sigma-forming element (PSFE). Isothermal aging treatments at 900–1100 °C for 20 h were performed, since CALPHAD and VEC/PSFE predictions diverged. Both prediction methods were compared with experimental characterization by a combination of scanning electron microscopy and high-energy synchrotron X-ray diffraction. In conclusion, the predictions from the VEC/PSFE and CALPHAD calculations (depending on the database used) were shown to be quite accurate.

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

Phase Transitions and Melting in the Venusian Basaltic Crust: Implications for Crustal Recycling

Understanding phase transitions and melting in the Venusian crust and the associated changes in density are critical to constrain crustal thickness, recycling, and remelting processes. For example, average surface conditions on Venus of 92 bars and 460 °C correspond to low-grade metamorphic conditions on Earth equivalent to hornfels- and greenschist-facies. Further, the primary igneous mineralogy in the Venusian crust should be replaced by metamorphic minerals over time and potentially melt, as temperatures and pressures increase with depth along geothermal gradients. Here, we use Perple_X (Connolly, 2005), a Gibbs free energy minimization program to calculates table phase equilibria over a range of pressure and temperature conditions for whole-rock compositions of a dry basalt, alkali basalt, and peridotite. Further, we extract mineral abundancies along five geotherms of 5, 10, 15, 20 and 25 °C/km and calculate the rock density after the extraction of melt. On the coldest geotherm of 5 °C/km, melting would start at ~ 197 km and ~121 km in the basalt and alkali basalt, respectively, and therefore at depths that significantly exceed the estimated thickness of the Venusian crust (8-25 km). A density cross-over, where the crustal density would exceed that of the mantle could potentially induce delamination and is estimated to occur at ~ 40 km for both compositions. On the 10 °C/km geotherm, melting of the basalt starts at 65 km and melt extraction only causes a gradual densification due the formation of initially lower melt proportions. In contrast, the alkali basalt starts to melt at a depth of ~ 57 km and the subsequent densification of the residual composition could trigger delamination at ~ 60km depth. For hotter geotherms (i.e., 25 °C/km gradient), melting would occur at shallower depths of 21 – 23 km. However, the basalt will not get significantly denser than the mantle while the residuum of the alkali basalt would reach a density cross-over at 30 km. Hence, differences in Venusian crustal compositions can significantly influence the thickness of the crust on individual geotherms. Our results demonstrate that phase transitions and melting are strongly dependent on the assumed geotherm and could, combined with geodynamic models, further constrain crustal parameters.

Julia Semprich↗

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 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

Phase Stability Through Machine Learning

Understanding the phase stability of a chemical system constitutes the foundation of materials science. Knowledge of the equilibrium state of a system under arbitrary thermodynamic conditions provides valuable information about the types of phases that are likely to be synthesized and how to get there. Accessing the phase diagram in a materials system provides one with the information necessary to design materials and microstructures with optimal properties. While the materials science community has long been focused on exploiting this knowledge to navigate the materials space, recent advances in machine learning (ML) and artificial intelligence (AI) have provided the community with novel ways of interrogating the materials thermodynamics space. Furthermore, this work presents some of the most recent advances in ML/AI applied to phase stability and thermodynamics of materials. Prof. John Morral always had a passion for understanding and teaching the fundamental characteristics of phase diagrams. This review is written to honor his memory.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Effect of alloying with Sc, Nb and Zr on reduction-diffusion synthesis of magnetically hard Sm(Fe,Co,Ti) 12 -based monocrystalline powders

Powders of Sm(Fe,Co) 11.2 Ti 0.8 alloys modified with Sc, Nb and Zr, as well as with additional Ti were prepared by reducing mechanically activated raw oxides with Ca metal in the furnace preheated to 990–1250 °C. Expansion of the crystal lattice upon introduction of Nb or additional Ti implies that atoms of these elements replace the smaller Fe atoms in the tetragonal ThMn 12 -type structure. On the other hand, contraction of the lattice upon introduction of Sc or Zr was smaller than what was expected for replacement of the Sm atoms, which suggests that the Sc and Zr atoms replace both the Sm and Fe atoms. Washing away the reduction byproducts expands the crystal lattice of the 1:12 particles and increases their coercivity. The lattice expansion associated with the washing is believed to be caused by interstitial H atoms; more research, however, is needed to establish the mechanism(s) of the washing effect on the coercivity. The earlier reported development of a high coercivity in zirconium-modified monocrystalline particles achieved by increasing the reduction annealing temperature to ≈1200 °C was similarly characteristic of the particles modified with Sc (the coercivity reaches 11.5 kOe) and Nb (8.1 kOe), but not for the particles prepared with additional Ti where the maximum coercivity of 8.3 kOe develops for a lower annealing temperature. Furthermore, it is concluded that Sc, Nb and Zr modify the high-temperature phase equilibria of the Sm(Fe,Co) 11.2 Ti 0.8 alloys allowing for an effective high-temperature processing, whereas the alloy coercivity increases with the synthesis temperature through a different, still unknown mechanism which may involve suppression of the defects specific to the 1:12 crystals.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Thermodynamic assessment of lithium halide reciprocal salt systems for energy applications

Lithium halides are of interest because of their unique applications as fuel or coolants in molten salt reactors (MSRs) and electrolytes in thermally activated batteries. Because iodine is an important fission product for MSRs, the Molten Salt Thermal Properties Database-Thermochemical (MSTDB-TC) is being expanded to include relevant iodide reciprocal salt systems. In this work, we assess the thermochemical properties of lithium halide pseudobinary systems based on reported experimental phase diagrams, eutectic temperatures, and enthalpies of mixing. Within the framework of the modified quasi-chemical model in the quadruplet approximation (MQMQA), extrapolations to the related pseudoternary and pseudoquaternary systems have been performed. The pseudoquaternary representations of phase equilibria, enthalpy, entropy, and heat capacity are represented in the resultant models. Of interest for generating source terms for accident analysis is the ability to compute the vapor pressures of the dominant iodide species for the systems at any composition. Here, an example of such calculations is provided for the LiF-LiI system at x LiI = 0.01 and x LiI = 0.70, in this case indicating that the LiI vapor species are predominant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Modeling aqueous association constants and mineral solubilities at subcritical and supercritical temperatures

The need for sustainable power generation has increased interest in the use of hydrothermal fluids for industrial applications. New high-enthalpy geothermal systems and biowaste-to-fuel processes are two relevant examples that employ supercritical fluids which require an in-depth understanding of complex chemical reactions occurring near the supercritical temperature of water (374 °C). As these processes operate in thermodynamic regimes that are not currently covered by a standard molar Gibbs energy of formation model, only empirical fits for single reaction systems are available which limit the use of multi-component phase equilibria calculations that are standard practice for less extreme environments. Here, in this paper, we advance a standard molar Gibbs energy of formation model able to operate in these otherwise inaccessible thermodynamic states to include species needed for key mineral solubility systems and ion association reactions. This work extends a model based on molecular statistical thermodynamics (MST) into four new systems (Na 3 PO 4 -H 2 O, LiOH-H 2 O, KOH-H 2 O, and BaSO 4 -H 2 O) by extending the model to cover 10 new species. For each of these systems, model predictions were consistently within the experimental uncertainties for the new systems covered. A breakdown of MST contributions to the model revealed that electrostatic and hard sphere contributions were key to reproducing density dependencies of standard molar Gibbs energy of formation values around the critical point of water.

42 ENGINEERING↗

DOE OSTI2022

Leveraging computational thermodynamics to guide SiC-ZrC chemical vapor deposition process development

Here using the CALPHAD approach to understand zirconium carbide deposition, a series of phase equilibria were calculated from a custom thermodynamic database based on a literature source, and the equilibria were used to explore the potential chemical vapor deposition (CVD) processing space in the ZrCl 4 -CH 3 SiCl 3 -CH 4 -H 2 system as a function of pressure, temperature, and gas composition. Several gas ratios were considered. At a given ZrCl 4 :CH 3 SiCl 3 ratio within the range studied, the most important factor was found to be the ratios of CH 4 :ZrCl 4 , wherein the nature of the composition – carbide vs. silicide – could be controlled. A pure binary composition of ZrC and SiC is expected to form by increasing the initial amount of methane and decreasing the amount of hydrogen from values predicted purely based on thermodynamic equilibrium. Rietveld analysis of the x-ray diffractograms from corresponding experimental depositions confirmed that increasing the CH 4 :ZrCl 4 ratio increased the fraction of carbon-containing species (SiC, ZrC) and decreased the fraction of non-carbides (ZrSi, ZrSi 2 , etc.), as predicted from the CALPHAD results.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Electrically Powered High-Salinity Brine Separation Using Dimethyl Ether

Dewatering highly saline aqueous streams, from mining and geothermal leachates to industrial wastewater, is essential for effective resource recovery and safe disposal. Membraneless water extraction (MWE) uses a low-polarity solvent to separate water from concentrated aqueous solutions. In this study, we design a new MWE that uses dimethyl ether (DME) to selectively extract water from high-salinity brines, leveraging the volatility of DME to achieve rapid solvent recovery. By separating water and dissolved salts at a liquid–liquid interface, MWE minimizes the deleterious effects of scaling on vulnerable membrane and heat exchanger surfaces, reducing the need for extensive pretreatment and expensive materials. We begin by developing a computational framework for a multistage counterflow liquid–liquid contactor, which extracts water into DME, coupled with a multistage solvent regenerator that uses vapor compression to efficiently separate the desalinated water from the DME extractant. Excess Gibbs free energy and equation of state frameworks are used to model fluid phase equilibria in water–DME–sodium chloride (NaCl) mixtures, with interaction parameters estimated from experimental data. Incorporating equilibrium calculations into a system-scale computational model, we examine the performance of MWE using DME for the first time. Our analysis demonstrates that MWE can concentrate seawater desalination brine (>1.0 mol NaCl kg –1 ) to zero-liquid discharge salinities, with an energy consumption of under 50 kW h per m 3 of water extracted with a solvent recovery ratio greater than 99.9%. We highlight the importance of staging the vapor compression process to simultaneously minimize energy consumption while enabling brine concentration and product water solvent contamination. Finally, the thermodynamic framework developed here allows for the robust evaluation of new MWE solvents and systems for critical brine concentration and fractional precipitation applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Dislocation-induced stop-and-go kinetics of interfacial transformations

Most engineering materials are based on multiphase microstructures produced either through the control of phase equilibria or by the fabrication of different materials as in thin-film processing. In both processes, the microstructure relaxes towards equilibrium by mismatch dislocations (or geometric misfit dislocations) across the heterophase interfaces (1-5) . Despite their ubiquitous presence, directly probing the dynamic action of mismatch dislocations has been unachievable owing to their buried nature. In this work, using the interfacial transformation of copper oxide to copper as an example, we demonstrate the role of mismatch dislocations in modulating oxide-to-metal interfacial transformations in an intermittent manner, by which the lateral flow of interfacial ledges is pinned at the core of mismatch dislocations until the dislocation climbs to the new oxide/metal interface location. Together with atomistic calculations, we identify that the pinning effect is associated with the non-local transport of metal atoms to fill vacancies at the dislocation core. These results provide mechanistic insight into solid-solid interfacial transformations and have substantial implications for utilizing structural defects at buried interfaces to modulate mass transport and transformation kinetics.

25 ENERGY STORAGE↗

DOE OSTI2025

MCCCS-MN

The MCCCS‒MN (Monte Carlo for Complex Chemical Systems‒Minnesota) software is developed by the Siepmann research group at the University of Minnesota. MCCCS‒MN allows for the simulation of multi-component molecular systems in the canonical, isobaric-isothermal (including constant stress for solids), grand-canonical, semi-grand, and Gibbs (NVT, NPT, and more than two simulation boxes) ensembles. It uses the configurational-bias Monte Carlo method to efficiently sample phase space for linear, branched and cyclic chain molecules, the adiabatic nuclear and electronic sampling Monte Carlo method to treat many-body polarization effects, and the aggregation-volume-bias Monte Carlo algorithm to efficiently sample the spatial distribution of associating molecules. MCCCS-MN employs a molecular representation of the system where force fields contain bonded and non-bonded terms. Funding for the development of MCCCS-MN through grants from the National Science Foundation (simulation of fluid phase equilibria and chromatography) and the Department of Energy (simulation of adsorption equilibria) is gratefully acknowledged.

Siepmann, J.Ilja [University of Minnesota - Twin C↗