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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

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.↗
DOE OSTI2021

Passive deconfinement of runaway electrons using an in-vessel helical coil

A helical coil designed to passively generate non-axisymmetric fields during a plasma disruption is shown (via electromagnetic analysis, linear MHD modeling, and relativistic drift orbit tracing) to be effective at deconfining runaway electrons (REs) on a time scale significantly faster than the plasma current quench. Magnetic equilibria from DIII-D RE-producing scenarios are used to calculate the toroidal electric field generated during the current quench phase of a disruption, which in turn drives current in the proposed n = 1 in-vessel helical coil, without the need for any external power supplies or disruption detection or prediction techniques. Simulations of the plasma evolution using the TokSys GS Evolve code predict the inductive coupling of coil currents up to 12% of the pre-disruption plasma current into the helical coil. The coil geometry is parametrically varied to maximize both the non-resonant and resonant components of the 3D magnetic perturbation, resulting in δB/B ≈ 10 –2 and a vacuumisland overlapwidth of up to 0.7ψ N . The REORBIT module of the MARS-F code is used to model the full non-axisymmetric magnetic field and trace RE drift orbits to determine the effect on RE deconfinement, with up to 70% of the RE orbits lost after 0.2 ms. A two-stage evolution of the RE orbit loss fraction is observed to be caused by resonant trapping between multiple magnetic island chains. Finally, electromagnetic and thermal stresses on the coil are calculated to be within operational limits for installation in DIII-D, and scale favorably to a reactor-size device. Furthermore, these findings motivate future experimental study of the helical coil concept in DIII-D or other tokamaks.

3D coil↗
NASA NTRSDate not supplied

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↗
NASA NTRS1988

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.↗
NASA NTRS2023

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↗
DOE OSTI2021

Neoclassical Tearing Mode Seeding by Nonlinear Three-Wave Interactions in Tokamaks

We report the experimental observation of seed magnetic island formation by nonlinear three-wave coupling of magnetic island triplets. In this experiment, disruptive 2,1 islands are seeded by the coupling of 4,3 and 3,2 tearing modes to a central 1,1 sawtooth precursor. Three-wave interactions between these modes are conclusively identified by bispectral analysis, indicating fixed phase relationships in agreement with theory. Furthermore, this new observation of this seeding mechanism has important implications for future reactors that must operate in stable plasma equilibria, free of disruptive 2,1 islands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

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↗
NASA NTRS2004

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.↗
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

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↗
NASA NTRS2004

The Crystallization of Canavalin as a Function of pH and NaCl Concentration

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 crystalline protein can be readily dissolved by weakly basic solution, which has been proposed to result 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.4 to 9.6 range. Fluorescence anisotropy, light scattering, and gel filtration experiments show that the solutions are primarily trimers, with association to form larger species occurring as a function of protein concentration.

Forsythe, Elizabeth L.↗
NASA NTRS1977

Solubility of K in Fe-S liquid, silicate-K/Fe-S/liq equilibria, and their planetary implications

Potassium has been found to have extremely limited absolute solubility in Fe-S liquid in the pressure-temperature range of 18 to 40 kbars, 1050 to 1150 C, and fO2 within the field of metallic iron. It also partitioned into a certain silicate phase highly in preference to Fe-S liquid at 30 kbar and 1100 C. The dependence of the partitioning of K between solid silicate and Fe-S liquid on fO2 and compositions of mineral solid solutions have been analyzed. These experimental data, along with those of others, limit the amount of K that could fractionate in Fe-S liquid layers or a core in the early history of the moon and, thus, act as localized heat sources in its thermal history models; the data also seem to argue against a chondritic abundance of potassium for earth. The question of fractionation of enough K-40 in an Fe-S liquid outer core of earth to provide the necesary thermal energy for the geomagnetic dynamo remains unresolved.

Gangully, J.↗
DOE OSTI2021

High-fidelity kinetic modeling of instabilities and gyromotion physics in nonuniform low-beta plasmas

A fourth-order accurate continuum kinetic Vlasov solver and a systematic method for constructing customizable kinetic equilibria are demonstrated to be powerful tools for the study of nonuniform collisionless low-beta plasmas. The noise-free methodology is applied to investigate two gradient-driven instabilities in 4D (x,y,vx,vy) phase space: the Kelvin–Helmholtz instability and the lower hybrid drift instability. Nonuniform two-species configurations where ion gyroradii are comparable to gradient scale lengths are explored. The approach sheds light on the evolution of the pressure tensor in Kelvin–Helmholtz instabilities and demonstrates that the associated stress tensor deviates significantly from the gyroviscous stress tensor. Even at high magnetization, first-order approximations to finite-gyromotion physics are shown to be inadequate for the Kelvin–Helmholtz instability, as shear scales evolve to become on par with gyromotion scales. The methodology facilitates exploring transport and energy partitioning properties associated with lower hybrid drift instabilities in low-beta plasma configurations. Distribution function features are captured in detail, including the formation of local extrema in the vicinity of particle-wave resonances. Finally, the approach enables detailed targeted investigations and advances kinetic simulation capability for plasmas in which gyromotion plays an important role.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRS2005

Optimal Sampling of a Reaction Coordinate in Molecular Dynamics

Estimating how free energy changes with the state of a system is a central goal in applications of statistical mechanics to problems of chemical or biological interest. From these free energy changes it is possible, for example, to establish which states of the system are stable, what are their probabilities and how the equilibria between these states are influenced by external conditions. Free energies are also of great utility in determining kinetics of transitions between different states. A variety of methods have been developed to compute free energies of condensed phase systems. Here, I will focus on one class of methods - those that allow for calculating free energy changes along one or several generalized coordinates in the system, often called reaction coordinates or order parameters . Considering that in almost all cases of practical interest a significant computational effort is required to determine free energy changes along such coordinates it is hardly surprising that efficiencies of different methods are of great concern. In most cases, the main difficulty is associated with its shape along the reaction coordinate. If the free energy changes markedly along this coordinate Boltzmann sampling of its different values becomes highly non-uniform. This, in turn, may have considerable, detrimental effect on the performance of many methods for calculating free energies.

Pohorille, Andrew↗
DOE OSTI2023

Chaotic and integrable magnetic fields in one-dimensional hybrid Vlasov–Maxwell equilibria

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗
NASA NTRS1990

Chemical heterogeneity and imperfect mixing in the solar nebula

The 'standard' accretion-disk model of the solar nebula is examined together with its basic equations, which are solved for the steady-state case. It is shown that, for this model of the solar nebula, the efficiency of mixing the products of thermochemical processing from small to large disk radii depends not on the magnitude of the eddy diffusivity but only on the ratio k between the eddy diffusivity and the eddy viscosity. In most cases, this ratio k is probably of order of unity or less, implying that most of the solar nebula was not contaminated by the consequences of thermochemical equilibria that were established at 'small' disk radii. It is noted that, while this condition is almost certainly satisfied if the physical process responsible for viscosity is thermal convection or waves, it has uncertain validity during the early phases of disk evolution where accretion-induced shear instabilities may dominate.

Stevenson, David J.↗
DOE OSTI2023

Geochemical Modeling Using OLI Systems Mixed Solvent Electrolyte Model in Support of the Waste Isolation Pilot Plant

Aqueous electrolyte thermodynamic models were developed for synthetic Salado (GWB) and Castile (ERDA-6) brines to inform the experimental design and to facilitate laboratory data interpretation in geochemical studies supporting disposal in the Waste Isolation Pilot Plant (WIPP). The thermodynamic models, which include the relevant vapor-liquid-solid equilibria for the oxidation-reduction conditions of interest, enable computation of system parameters such as activity-based pH, E h [vs. standard hydrogen electrode (SHE)], ionic strength, total dissolved solids, and solution density; the chemical speciation of all phases; saturation tendencies for thermodynamically feasible solids; and changes in brine chemistry because of the addition of H 2 (g), argon, MgO, Fe(0), and other waste form components. The impacts of oxidation-reduction potential (ORP) and pH on aqueous speciation and dominant solid phases were also assessed with the aid of real-solution Pourbaix diagrams for key elements of interest. In addition to the OLI modeling results, limited measurement data for pH and ORP standards are provided, analyzed, and discussed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1991

The high-pressure phase diagram of Fe(0.94)O - A possible constituent of the earth's core

Electrical resistivity measurements to pressures of 83 GPa and temperatures ranging from 300 K to 4300 K confirm the presence of both crystalline and liquid metallic phases of FeO at pressures above 60-70 GPa and temperatures above 1000 K. By experimentally determinig the melting temperature of FeO to 100 GPa and of a model-core composition at 83 GPa, it is found that the solid-melt equilibria can be described by complete solid solution across the Fe-FeO system at pressures above 70 GPa. The results indicate that oxygen is a viable and likely candidate for the major light alloying element of the earth's liquid outer core. The data suggest that the temperature at the core-mantle boundary is close to 4800 K and that heat lost out of the core accounts for more than 20 percent of the heat flux observed at the surface.

Knittle, Elise↗