Engineering Papers⌕ Search

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.

At least 145 records · Page 8

On the topological stability of magnetostatic equilibria

The topological stability of MHD equilibria is investigated by exploring the formal analogy, in the ideal MHD limit, between the topology of magnetic lines of force in coordinate space and the topology of integral surfaces of one- and two-dimensional Hamiltonian systems in phase space. It is demonstrated that in an astrophysical setting, symmetric magnetostatic equilibria satisfying the ideal MHD equations are exceptional. The principal result of the study is that previous infinitesimal perturbation theory calculations can be generalized to include finite-amplitude and symmetry-breaking effects. The effect of the ergodicity of perturbed symmetric equilibria on heat dispersal in magnetically dominated plasmas is discussed.

Tsinganos, K. C.↗

Determination of Systems Suitable for Study as Monotectic Binary Metallic Alloy Solidification Models

Succinonitrile-water and diethylene glycol-ethyl salicylate are two transparent systems which have been studied as monotectic binary metallic alloy solidification models. Being transparent, these systems allow for the direct observations of phase transformations and solidification reactions. The objective was to develop a screening technique to find systems of interest and then experimentally measure those systems. The succinonitrile-water system was used to check the procedures. To simulate the phase diagram of the system, two computer programs which determine solid-liquid and liquid-liquid equilibria were obtained. These programs use the UNIFAC method to determine activity coefficients and together with several other programs were used to predict the phase diagram. An experimental apparatus was developed and the succinonitrile-water phase diagram measured. The diagram was compared to both the simulation and literature data. Substantial differences were found in the comparisons which serve to demonstrate the need for this procedure.

Smith, J. E., Jr.↗

Laboratory studies of the nitric acid trihydrate - Implications for the south polar stratosphere

Vapor pressures of HNO3 and H2O have been measured over the trihydrate crystal, formed by vapor deposit on a glass surface. In the temperature range 190 to 205 K the two phase-equilibrium trihydrate/vapor was studied by adding and removing H2O. Coexistence equilibria vapor pressures of trihydrate/solid solutions of HNO3 in ice and of mono-/trihydrate were also measured. Results show that for typical mixing ratios of H2O and HNO3 found in the lower stratosphere (3 ppm H2O, 5 ppb HNO3) the trihydrate would start to form at temperatures about 7 K higher than the ice point. The pressure of atmospheric HNO3 would rapidly decrease as the atmosphere cools without large changes in partial pressures of H2O. These laboratory results provide information on the formation of polar stratospheric clouds containing H2O and HNO3.

Hanson, David↗

Machine Learning Applications to Metal-Silicate Equilibria and their Insights into Core Formation

An extensive number of studies have experimentally investigated how elements distribute between metal and silicate phases, to better constrain core-mantle chemical equilibrium. Here, we present a new database compiling all (to our knowledge) experimental data on liquid metal-silicate partitioning from 118 peer-reviewed publications. We applied various machine learning techniques to gain further insights into these partitioning equilibria and their dependencies. We performed a network analysis to investigate the relationship between experiments and partition coefficients, which enables visualizing gaps in the experimental dataset and biases related to varying experimental conditions and analytical setup. In addition, semi-empirical thermodynamic models are commonly used to extrapolate these chemical reactions to the wide range of pressure, temperature and compositional conditions of planetary differentiation. These models are based on linear regressions that assume continuous relationship between partition coefficients and experimental variables. Here, we considered random forest regressions, which are algorithms based on ensembles of decision trees and does not consider continuous effects of each variable. The application of this regression significantly improves the prediction of metal-silicate partitioning for several elements including Ni, Si and Cr. We will show how this new approach improves our understanding of elemental exchange between metal and silicate and their implications for the Earth’s core formation.

siderophile element↗

Liquid-solid equilibria involving spinel, ilmenite, and ferropseudobrookite in the system 'FeO'-Al2O3-TiO2 in contact with metallic iron

Phase relations in the liquidus temperature region of the system 'FeO'-Al2O3-TiO2 in contact with metallic iron, at a total pressure below 1 atm, have been determined by the quenching technique. Four invariant points have been located, with phase assemblages and temperatures as follows; wuestite, ulvoespinel, nercynite and liquid, 1306 C; ulvoespinel, ilmenite, ferropseudobrookite and liquid, 1340 C; ulvoespinel, hercynite, ferropseudobrookite and liquid, 1367 C; hercynite, ferropseudobrookite, corundum and liquid, 1465 C. The data obtained confirm the presence of a miscibility gap between titanate and aluminate spinels, and provide quantitative data for the effect of Al2O3 on mutual stability relations among spinel, ilmenite, and ferropseudobrookite phases in the presence of liquid at high temperatures and strongly reducing conditions. It is shown that Al2O3 has a strong stabilizing effect on the phase assemblage ferropseudobrookite and spinel relative to ilmenite.

Schreifels, W. A.↗

Thermodynamic models in cosmochemical systems.

Generalized computer methods are developed for inferring details of the formation of cosmochemical systems. Compositions of ideal gas mixtures existing in equilibrium with multicomponent solid and liquid phases are calculated. A comparison of computed results with experimental data is made for the ternary system MgO-FeO-SiO2. While the ideal-solution approximation is shown to be inaccurate in dealing with the silicate melts, the stable phases and compositions can be accurately calculated in a system where there are only solids and gas. A model system containing the elements H, O, Si, Mg, S, C, Cl, and F is investigated over a range of compositions involving the gas and ten solid phases, to show the power of the technique in dealing with complex gas-solid equilibria. Systems close to cosmic composition are next considered, both with and without iron.

Griffiths, P. R.↗

A thermodynamic model for subsolidus equilibria in the system CaO-MgO-Al2O3-SiO2

It is shown that the high temperature enthalpy of solution data for pure phases and solid solutions in the CMAS system are, for 'gabbroic' and 'peridotitic' compositions, consistent with available phase equilibrium data for the MAS, CAS and CMAS systems. A refined set of values of thermodynamic properties for these phases and solid solutions is tabulated. The small differences between the new data set and that of Helgeson et al. (1978) arise from new data on heat capacity and enthalpy being incorporated. The important constraints applied and the major difficulties in fitting which arose are summarized.

Wood, B. J.↗

Chassigny petrogenesis - Melt compositions, intensive parameters, and water contents of Martian (?) magmas

This paper considers petrogenesis of SNC meteorites by examining partially crystallized melt inclusions in cumulus olivine grains in Chassigny meteorite. First, all phases in the Chassigny inclusions were analyzed by electron microprobe, followed by an experimental investigation of the kaersutite/melt equilibria. Results yield the intensive parameters which stabilize igneous kaersutite and the melt compositions that coexist with kaersutite. The results are used to model Chassigny petrogenesis. It is argued that the existence in the Chassigny of hydrous melts and the conditions appropriate for amphibole crystallization suggest the existence on the Chassigny parent body (Mars) of SiO2-rich lavas.

Johnson, Marie C.↗

Predicting major element mineral/melt equilibria - A statistical approach

Empirical equations have been developed for calculating the mole fractions of NaO0.5, MgO, AlO1.5, SiO2, KO0.5, CaO, TiO2, and FeO in a solid phase of initially unknown identity given only the composition of the coexisting silicate melt. The approach involves a linear multivariate regression analysis in which solid composition is expressed as a Taylor series expansion of the liquid compositions. An internally consistent precision of approximately 0.94 is obtained, that is, the nature of the liquidus phase in the input data set can be correctly predicted for approximately 94% of the entries. The composition of the liquidus phase may be calculated to better than 5 mol % absolute. An important feature of this 'generalized solid' model is its reversibility; that is, the dependent and independent variables in the linear multivariate regression may be inverted to permit prediction of the composition of a silicate liquid produced by equilibrium partial melting of a polymineralic source assemblage.

Hostetler, C. J.↗

Spinel cataclasites in 15445 and 72435 - Petrology and criteria for equilibrium

The problem of establishing the existence of equilibrium among the coexisting phases in the rock is addressed by presenting petrographic and mineral chemistry data on a new spinel cataclasite from 15445 (clast H) and data more extensive than those previously available on two clasts in 72435. Criteria useful in reconstructing the original petrology of these and other spinel cataclasites are analyzed by considering equilibrium among the different phases, that is, the mono- or polymict nature of these cataclasized samples. Finally, the role of impact processes in disturbing the equilibria is discussed.

Baker, M. B.↗

Thermodynamic and Diffusion Model Estimates on Metamorphic Temperatures and Timescales for Basaltic Eucrite GRA 98098

Introduction: HED meteorites are thought to rep-resent igneous rocks from Vesta’s basaltic crust and preserve evidence of early crustal metamorphism. Determining the temperatures and timescales of thermal metamorphism is important for reconstructing crustal evolution in the early solar system. Here, we study basaltic eucrite Graves Nunataks (GRA) 98098, which has been identified as a highly metamorphosed eucrite [1]. We present new estimates on metamorphic temperatures determined via thermodynamic modeling as well as the initial results from diffusion models constraining timescales of thermal metamorphism. Sample Description: GRA 98098 is an unbrecciated eucrite with a granoblastic plagioclase and pyroxene mineralogy. Millimeter to cm-long lathes of tridymite cross-cut and poikilitically enclose plagioclase and pyroxene [1,this work]. Pyroxene grains have exsolved into Ca-rich (~Wo38En29Fs33) and Ca-poor (~Wo4.5En36Fs59.5) lamellae. Both unzoned and zoned plagioclase grains are observed. Unzoned plagioclase grains are found solely with tridymite laths. These grains have ~An92 compositions. The cores of the zoned plagioclase grains have the same composition and thin, relatively sodic rims (~An67), (Fig. 1). The bulk sample is unusually enriched in highly in-compatible elements and has one of the most fractionated REE patterns reported [1]. Maximum metamorphic temperatures of 985±78°C have been estimated using two-pyroxene thermometry [2]. Methods: Thermodynamic modeling. Thermodynamic models were constructed using the software Perple_X, which employs a Gibbs free energy minimization in order to determine the most stable phase assemblage for a given bulk rock composition [3]. Bulk composition was calculated using mineral com-positions acquired via EMPA (this study) and the observed abundancies present in the thin section. Two bulk compositions were estimated; 1) includes all phases present in the thin section, (assumes that all phases are present during metamorphism), 2) excludes tridymite from the bulk calculation (assumes that tridymite was not present during metamorphism). In order to determine whether metamorphic equilibria was achieved and estimate temperatures of metamorphism, we compared measured pyroxene compositions with thermodynamically predicted compositions [4, Fig. 2]. Diffusion Modeling. Several time-temperature de-pendent diffusion profiles were calculated in order to determine the best match for XAn chemical profiles observed at the edges of the zoned plagioclase (Fig. 3). We assumed that the start condition was a stepwise gradient at the plagioclase/pyroxene interface. We also assumed an average diffusion coefficient (D) and a constant temperature using the equation in [5]. D was determined for two temperatures (T = 1060ºC; near eucrite solidus [6] and T = 985ºC; metamorphism reported in [4]) and then XAn was calculated as a function of distance from plagioclase core to rim using an error function solution to Fick’s second law. Results: Thermodynamic model results are summarized in Fig. 2. For a bulk composition that includes all phases in the thin section, pyroxene endmember compositions plot in the following temperature ranges: Fs ~660-860ºC, En~1000ºC & 1150ºC, and Wo~760-900ºC (Fig. 2a). For a bulk composition that excludes tridymite from the peak metamorphic assemblage (i.e., the bulk composition minus the contribution from tridymite), a temperature range could not be determined for the Fs component of pyroxene. For Wo, T~760-900ºC and En, T~1000ºC & 1150ºC (Fig. 2b). Fig. 3 summarizes the diffusion model results. For T = 1060°C & 985°C, the most appropriate time interval was estimated based on which diffusion curve most matched (solid lines, Fig. 3) the EMPA data. For T = 1060°C, the best looking match was t = 500 ka. For T = 985°C, the best match was t = 7 Ma. Discussion and future work: Temperature estimates from thermodynamic models are not conclusive because the temperature ranges determined for pyroxene endmember stability do not overlap (colored fields in Fig. 2), thus implying that there is disequilibrium between pyroxene crystals and the bulk composition considered [4]. Thus, additional exploration is needed to define a metamorphically equilibrated do-main that accurately records peak temperature. The utility of defining metamorphically equilibrated do-mains to improve the accuracy and level of detail elucidated regarding the petrogenetic history of metamorphose samples has been demonstrated previously [4,7]. We suggest that in the case of Fig. 2a, the thin section composition is not representative of the length scales over which metamorphic equilibrium was achieve and in the case of Fig. 2b, the assumption that tridymite was not present during metamorphism was incorrect. However, results from thermodynamic models can provide insight into the relative timing of mineral and compositional textures. For example from texture alone, it is unclear whether tridymite was igneous in origin and represents the last bits of melt in a crystallizing magma chamber, or if it formed during (and possibly initiated) open system thermal metamorphism. The latter could be consistent with a partial melt hypothesis [8,9] while the former implies that simple fractional crystallization can yield the textures present in GRA 98098. The lack of coincidence be-tween pyroxene endmember compositions in Fig. 2b suggest that the bulk composition minus tridymite was not the assemblage in equilibrium with the pyroxene, suggesting that tridymite was present during metamorphism and formed during igneous crystallization. We conclude that the development of the Na-rich plagioclase rims likely occurred during or immediately after peak thermal metamorphism, because eucrites of similar metamorphic grade and texture have unzoned plagioclase (~An92) [2,4,8], and Na zoning is only observed in the plagioclase not included in the tridymite. This suggests that the zoning formed after tridymite formation, and therefore after igneous crystallization. Thus, the timescales calculated via diffusion modeling possibly represent the time interval over which thermal metamorphism occurred. Cooling rates approximated for the Vestan crust predict that the crust cooled below 300°C around 35-40 Ma after formation[10]. This is consistent with our modeling results that predict formation of the Na rich plagioclase rims occurring at higher temperatures over a period of 0.5 to 7 Ma years. Future work. Additional thermodynamic modeling work will focus on selecting an equilibrated bulk rock domain in which to elucidate metamorphic conditions. Diffusion models currently provide a minimum time-scale, since diffusion slows down as the system cools. Future work with will focus on integrating cooling into the diffusion models and constraining the depth at which thermal metamorphism occurs because it could be used to determine whether the range of time-scales calculated for thermal metamorphism are consistent with the geologic environment.

J S Gorce↗

Interfacial tension and vapor-liquid equilibria in the critical region of mixtures

In the critical region, the concept of two-scale-factor universality can be used to accurately predict the surface tension between near-critical vapor and liquid phases from the singularity in the thermodynamic properties of the bulk fluid. In the present work, this idea is generalized to binary mixtures and is illustrated using the data of Hsu et al. (1985) for CO2 + n-butane. The pressure-temperature-composition-density data for coexisting, near-critical phases of the mixtures are fitted with a thermodynamic potential comprised of a sum of a singular term and nonsingular terms. The nonuniversal amplitudes characterizing the singular term for the mixtures are obtained from the amplitudes for the pure components by interpolation in a space of thermodynamic 'field' variables. The interfacial tensions predicted for the mixtures from the singular term are within 10 percent of the data on three isotherms in the pressure range (Pc - P)/Pc of less than 0.5. This difference is comparable to the combined experimental and model errors.

Moldover, Michael R.↗

Equilibrium Selection via Current Sheet Relaxation and Guide Field Amplification

Although there is a continuous spectrum of current sheet equilibria, how a particular equilibrium is selected by a given system remains a mystery. Yet, only a limited number of equilibrium solutions are used for analyses of magnetized plasma phenomena. Here we present the exact process of equilibrium selection, by analyzing the relaxation process of a disequilibrated current sheet under a finite guide field. It is shown via phase-space analyses and particle-in-cell simulations that the current sheet relaxes in such a way that the guide field is locally amplified, yielding a mixed equilibrium from the spectrum. Comparisons to spacecraft observations and solar wind current sheet statistics demonstrate that such mixed equilibria are ubiquitous and exist as underlying local structures in various physical environments.

Young-Dae Yoon↗

Stability of Pseudobrookite-Type Titanium Oxides

Orthorhombic, (Bbmm), (Al, Fe, Cr, Ti)(sub 2) TiO5-(Mg, Fe)Ti2O5 solid solutions (pseudobrookites, s.l.) are found either as an oxidation product of ilmenite and/or spinel or a primary crystallizing phase in igneous and metamorphic rocks on Earth (e.g., basalt flows, crustal and mantle xenoliths, hornfels), and basaltic rocks on the Moon. Moreover, orthorhombic oxides are often part of the crystalline matrix in glass/ceramics with useful applications, and play a major role in the industrial production of TiO2. To fully exploit the potential of these compounds as petrogenetic indicators and/or useful materials we need to quantitatively understand the factors controlling their properties and stability, and thus, to extrapolate beyond the calibrating experiments. For that purpose, we need to combine thermochemistry, phase equilibrium, and in situ P-V-T-cation disorder experimental data that presently either are incomplete or lacking. Perhaps, the most complete data set is that for MgTi2O5 (karrooite) which allows the calibration of models for the Gibbs free energy of the MgTi2O5 as a function of pressure, temperature, and the Mg2+-Ti4+ distribution between the two nonequivalent octahedral sites. Consequently, the effect of cation disorder on MgTi2O5 stability, and the phase relations among MgTi2O5, other titanium oxides, and silicate minerals can be examined. Calculated phase relations in the Mg-Ti-Si-O system and phase equilibrium experiments in Fe-bearing compositions suggest that pseudobrookite-type oxides may be a more common in rocks than previously realized. However, homogeneous and heterogeneous equilibria, and crystallization paths likely affect their stability. For example, isobaric increases in temperature favor disordering and thus entropy-stabilization, in contrast, isothermal increases in pressure have the opposite effect. Although, currently, the potential effect of composition to cation disorder cannot be fully explored, it appears that enrichment in trivalent cations probably enhances entropy-stabilization and thus may increase the stability of (Al, Fe, Cr, Ti)-rich pseudobrookites relative to that of (Mg, Fe)-rich ones. In addition, high-temperature, nearly isothermal, decompression paths of olivine+orthopyroxene+oxide assemblages may favor pseudobrookites (s.l.) over rutile and/or ilmenite, in contrast, cooling at low pressures seems to favor ilmenite and/or rutile. In the case of crustal and mantle xenoliths, the presence or absence of orthorhombic oxides is probably controlled by reactions with olivine, orthopyroxene, ilmenite, and rutile. In oceanic mantle xenoliths such reactions may also involve a TiO2-enriched but not SiO2-enriched melt/fluid, because pseudobrookites (s.l.) would react with the SiO2-enriched melt/fluid to form orthopyroxene and rutile. Parenthetically, experiments and model calculations in the Mg-Ti-Si-O system suggest that low degree partial melting of low-TiO2 bulk compositions may produce Ti-enriched liquids in equilibrium with olivine, orthopyroxen ad=nd MgTi2O5, rutile or ilmenite.

Xirouchakis, Dimistrios↗

Solution-Phase Processes of Macromolecular Crystallization

We have proposed, for the tetragonal form of chicken egg lysozyme, that solution phase assembly processes are needed to form the growth units for crystal nucleation and growth. The starting point for the self-association process is the monomeric protein, and the final crystallographic symmetry is defined by the initial dimerization interactions of the monomers and subsequent n-mers formed, which in turn are a function of the crystallization conditions. It has been suggested that multimeric proteins generally incorporate the underlying multimers symmetry into the final crystallographic symmetry. We posed the question of what happens to a protein that is known to grow as an n-mer when it is placed in solution conditions where it is monomeric. The trypsin-treated, or cut, form of the protein canavalin (CCAN) has been shown to nucleate and grow crystals as a trimer from neutral to slightly acidic solutions. Under these conditions the solution is composed almost wholly of trimers. The insoluble protein can be readily dissolved by weakly basic solution, which results in a solution that is monomeric. There are three possible outcomes to an attempt at crystallization of the protein under monomeric (high pH) conditions: 1) we will obtain the same crystals as under trimer conditions, but at different protein concentrations governed by the self association equilibria; 2) we will obtain crystals having a different symmetry, based upon a monomeric growth unit; 3) we will not obtain crystals. Obtaining the first result would be indicative that the solution-phase self-association process is critical to the crystal nucleation and growth process. The second result would be less clear, as it may also reflect a pH-dependent shift in the trimer-trimer molecular interactions. The third result, particularly for experiments in the transition pH's between trimeric and monomeric CCAN, would indicate that the monomer does not crystallize, and that solution phase self association is not part of the crystal nucleation and growth path. Results are presented for crystallization experiments of CCAN over the pH 6.8 to 9.6 range.

Pusey, Marc L.↗

Calcium and Magnesium Distribution Between Sulfide and Silicate Melts: Thermodynamic Modeling and Insights into Mercury's Mantle-Crust Differentiation

Mercury stands out as the terrestrial planet with the lowest degree of oxidation among all planets of our solar system. It is characterized by a high abundance of sulfur, which is expected to replace oxygen in silicate melts. This sulfur enrichment could lead to sulfide saturation in its mantle, potentially resulting in the formation of sulfides rich in calcium and magnesium. Utilizing experimental data from the literature, we constructed a thermodynamic model to explore the partitioning of calcium and magnesium between sulfide and silicate melts. This model aims to improve our comprehension of the chemical equilibria driving the process of mantle-crust differentiation on Mercury. Our models indicate that increased temperature, decreased pressure, and low concentrations of iron, oxygen, and carbon within the sulfide phase, along with low sulfur concentrations in the silicate phase, favor calcium’s partitioning into the sulfide. Additionally, magnesium partitioning into the sulfide melt is favored with decreased pressure, decreased temperature, low concentrations of iron and calcium within the sulfide, and low concentrations of sulfur within the silicate. Our results suggest that up 9.4% and 5.6% Ca and Mg, respectively, can be present in sulfides equilibrated with sulfide-saturated silicate melts during mantle-crust differentiation. Further investigation will focus on the distribution of sulfide phases within the mantle and crust to better constrain Mercury’s internal structure.

Taylor McCombs↗

Aqueous Alteration on the Ryugu Parent Body Constrained by Chemical Equilibrium Models

The composition of solid samples returned from the asteroid 162173 Ryugu byHayabusa2 indicates aqueous alteration of rocks that formed a parent body of Ryugu [1, 2]. As on parent bodies of carbonaceous chondrites [3-5], Ryugu materials could have went through (1) accretion of anhydrous and reduced rocky grains together with water-rich ices, (2) melting of ices in the body’s interior through release of radiogenic heat [6], and (3)aqueous alteration of rocky and organic materials that produced secondary minerals, salt- and organic-bearing water solutions, gases, and altered organic matter. Although the bulk chemical composition and major mineralogy (phyllosilicates, magnetite, pyrrhotite, carbonates, phosphates, etc.) of Ryugu [1,2]are similar to those of CI carbonaceous chondrites [3-5, 7], the lack of sulfates and a deficiency of oxygen in the returned samples suggest a distinct alteration pathway. Chondritic materials [3-5] do not indicate complete chemical equilibration though low-temperature parent body processes. However, calculations of chemical equilibria in modeling rock-water-gas type systems[e.g., 8, 9], could assess deviations of chondritic materials from equilibrium conditions and constrain partitioning of chemical elements between solid, aqueous, and gas phase for different stages of alteration. Here we constrained speciation of solid, aqueous, and gaseous phases during aqueous alteration on the Ryugu parent body through chemical equilibrium calculations.

M Yu Zolotov↗