SEARCH · Search NASA

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 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS2017

The Substorm Cycle as Reproduced by Global MHD Models

Recently, Gordeev et al. (2015) suggested a method to test global MHD models against statistical empirical data. They showed that four community-available global MHD models supported by the Community Coordinated Modeling Center (CCMC) produce a reasonable agreement with reality for those key parameters (the magnetospheric size, magnetic field, and pressure) that are directly related to the large-scale equilibria in the outer magnetosphere. Based on the same set of simulation runs, here we investigate how the models reproduce the global loading-unloading cycle. We found that in terms of global magnetic flux transport, three examined CCMC models display systematically different response to idealized2 h north then 2 h south interplanetary magnetic field (IMF) Bz variation. The LFM model shows a depressed return convection and high loading rate during the growth phase as well as enhanced return convection and high unloading rate during the expansion phase, with the amount of loaded unloaded magnetotail flux and the growth phase duration being the closest to their observed empirical values during isolated substorms. Two other models exhibit drastically different behavior. In the BATS-R-US model the plasma sheet convection shows a smooth transition to the steady convection regime after the IMF southward turning. In the Open GGCM a weak plasma sheet convection has comparable intensities during both the growth phase and the following slow unloading phase. We also demonstrate potential technical problem in the publicly available simulations which is related to post processing interpolation and could affect the accuracy of magnetic field tracing and of other related procedures.

Gordeev, E.↗

NASA NTRS1987

Particle paths and phase plane for time-dependent similarity solutions of the one-dimensional Vlasov-Maxwell equations

The phase trajectories of particles in a plasma described by the one-dimensional Vlasov-Maxwell equations are determined qualitatively, analyzing exact general similarity solutions for the cases of temporally damped and growing (sinusoidal or localized) electric fields. The results of numerical integration in both untransformed and Lie-group point-transformed coordinates are presented in extensive graphs and characterized in detail. The implications of the present analysis for the stability of BGK equilibria are explored, and the existence of nonlinear solutions arbitrarily close to and significantly different from the BGK solutions is demonstrated.

Roberts, Dana Aaron↗

NASA NTRS1988

Computer modeling of the mineralogy of the Martian surface, as modified by aqueous alteration

Mineralogical constraints can be placed on the Martian surface by assuming chemical equilibria among the surface rocks, atmosphere and hypothesized percolating groundwater. A study was made of possible Martian surface mineralogy, as modified by the action of aqueous alteration, using the EQ3/6 computer codes. These codes calculate gas fugacities, aqueous speciation, ionic strength, pH, Eh and concentration and degree of mineral saturation for complex aqueous systems. Thus, these codes are also able to consider mineralogical solid solutions. These codes are able to predict the likely alteration phases which will occur as the result of weathering on the Martian surface. Knowledge of the stability conditions of these phases will then assist in the definition of the specifications for the sample canister of the proposed Martian sample return mission. The model and its results are discussed.

Zolensky, M. E.↗

NASA NTRS2004

The Relationship Between Clinopyroxene Fe3+ Content and Oxygen Fugacity

The intrinsic oxygen fugacity (fO2) imposed on a magma has the ability to influence the crystallization sequence of the melt, as well as the composition of the resulting minerals. fO2 is an easily controlled parameter in the lab, either through gas-mixing equilibria or with a solid-state buffer assemblage. In nature, the fO2 of a closed system is imposed on the system internally through multivalent equilibria involving the phenocryst-melt assemblage. This results in a characteristic oxidation state. The physical parameter used to quantify oxidation state is oxygen fugacity. Iron is the only major rock forming element in basaltic melts to exist in multiple valence states and, therefore, it is commonly used to assess fO2. Traditional methods to quantify fO2 utilize the ferric content of glasses or coexisting Fe-Ti oxides. However, many rocks, such as the Martian meteorites, do not contain the necessary phases or have oxides which have suffered reequilibration, thereby rendering them unmeasureable by current techniques. For these rocks, new methods, utilizing other phases are needed. Mafic minerals have Fe(3+)/SigmaFe ratios that are a function of two factors: 1) crystal chemistry and 2) their intrinsic fO2 during crystallization. Olivine and orthopyroxene, for example, have steric constraints on the extent to which Fe(3+) can be incorporated in their structures, and may not record changes in magmatic fO2 in a way that can easily be measured. The chemistry of clinopyroxene, however, allows for extensive incorporation of Fe(3+) in its crystal structure, making it a potentially useful oxybarometer. To date, there have been few, if any, systematic experimental studies of the variation of the Fe(3+)/SigmaFe ratio as a function of fO2 in clinopyroxene. This study seeks to address this lack of data.

McCanta, M. C.↗

NASA NTRS2004

Oxygen Fugacity of the Martian Mantle from Pigeonite/Melt Partitioning of Samarium, Europium and Gadolinium

This study is part of an ongoing effort to calibrate the pyroxene/melt Eu oxybarometer for conditions relevant to the martian meteorites. There is fairly good agreement between a determinations using equilibria between Fe-Ti oxides and the estimates from Eu anomalies in shergottite augites in tenns of which meteorites are more or less oxidized. The Eu calibration was for angrite composition pyroxenes which are rather extreme. However, application of a calibration for martian composition augites 113 does not significantly reduce the discrepancy between the two methods. One possible reason for this discrepancy is that augites are non-liquidus. The use of pigeonite rather than augite as the oxy-barometer phase is considered. We have conducted experiments on martian composition pigeonite/melt REE partitioning as a function of fO2.

Musselwhite, S.↗

DOE OSTI2025

The Temperature Dependence of the Langmuir Adsorption Model for a Single-Site Metal–Organic Framework

The single-site Langmuir adsorption model, also known as the Langmuir isotherm equation, is one of the simplest possible descriptions of adsorption phenomena, and yet finds widespread applicability across a range of disciplines. In its simplest form, it is deployed to treat adsorption equilibria at constant temperature (i.e., along isotherms); however, at the heart of its derivation is a more general class of models that each incorporate an explicit temperature dependence, subject to assumptions about the spatial/translational degrees of freedom of the adsorbed species. In this work, measurements of the temperature dependence of supercritical adsorption of H 2 on a single-site metal-organic framework (MOF) are presented, and fitted using a range of Langmuir models with distinct treatments of degrees of freedom in the adsorbed phase. Surprisingly, all of the models can be used to adequately represent the measured data, despite yielding significantly different values for binding energy and the temperature dependence of the isosteric enthalpy of adsorption (i.e., the isosteric heat, qst). Furthermore, a critical finding is that the mean-temperature isosteric enthalpy of adsorption remains consistent across all models within experimental error, highlighting its reliability for evaluating adsorption thermodynamics.

08 HYDROGEN↗

NASA NTRS1987

The evolution of a coronal streamer prior to mass ejection

A model describing the quasi-static evolution of a coronal helmet streamer as it is inflated with excess mass was developed. The model produces a sequence of magnetostatic equilibria for a magnetic field configuration that includes a current sheet in the equatorial plane as well as distributed volume currents. The equilibrium solutions are given as Legendre series, and the height at which the current sheet starts is adjusted iteratively to maintain constant flux in the closed-field region. The model was fit to SMM coronagraph/polarimeter observations showing the slow growth of a coronal streamer prior to its disruption in association with an eruptive prominence and mass ejection on August 18, 1980. The results suggest that the early phase of this event is consistent with the quasi-static evolution of the corona in response to the slow addition of mass to the closed-field region of the streamer.

Wolfson, Richard↗

NASA NTRS1975

The system anorthite-forsterite-fayalite-silica to 2 kbar with lunar petrologic applications

Near-liquidus phase relations in the magnesium-rich portion of the system anorthite-forsterite-fayalite-silica have been determined from synthesis experiments at 1 atm and 2 kbar at oxygen fugacities near the iron-wuestite buffer. The major-element chemistries of KREEP basalts are similar to peritectic liquids, L/ol,px,an/, produced experimentally. The reaction L/ol,px,an/ is well-located at low pressures. It intersects the join anorthite-enstatite-ferrosilite at Fe/(Fe+Mg) equals 0.45 (atomic ratio), and is peritectic only for Fe/(Fe+Mg) not greater than 0.45. Although positions of the equilibria, L/an,ol/ and L/ol,px/, differ significantly at 1 atm from their positions at 2 kbar, the position of the piercing point, L/ol,px,an/, on the anorthite-enstatite-ferrosilite join is not sensitive to pressure. The composition of liquids produced by partial melting of troctolites and pyroxenites in the lunar crustal environment will be affected by pressure, whereas liquids derived from gabbroic assemblages will not.

Merrill, R. B.↗

NASA NTRS1990

Liquid-Vapor Equilibrium of Multicomponent Cryogenic Systems

Liquid-vapor and solid-vapor equilibria at low to moderate pressures and low temperatures are important in many solar system environments, including the surface and clouds of Titan, the clouds of Uranus and Neptune, and the surfaces of Mars and Triton. The familiar cases of ideal behavior are limiting cases of a general thermodynamic representation for the vapor pressure of each component in a homogeneous multicomponent system. The fundamental connections of laboratory measurements to thermodynamic models are through the Gibbs-Duhem relation and the Gibbs-Helmholtz relation. Using laboratory measurements of the total pressure, temperature, and compositions of the liquid and vapor phases at equilibrium, the values of these parameters can be determined. The resulting model for vapor-liquid equilibrium can then conveniently and accurately be used to calculate pressures, compositions, condensation altitudes, and their dependencies on changing climatic conditions. A specific system being investigated is CH4-C2H6-N2, at conditions relevant to Titan's surface and atmosphere. Discussed are: the modeling of existing data on CH4-N2, with applications to the composition of Titan's condensate clouds; some new measurements on the CH4-C2H6 binary, using a high-precision static/volumetric system, and on the C2H6-N2 binary, using the volumetric system and a sensitive cryogenic flow calorimeter; and describe a new cryogenic phase-equilibrium vessel with which we are beginning a detailed, systematic study of the three constituent binaries and the ternary CH4-C2H6-N2 system at temperatures ranging from 80 to 105 K and pressures from 0.1 to 7 bar.

Thompson, W. Reid↗

NASA NTRS1998

Rare Potassium-Bearing Mica in Allan Hills 84001: Additional Constraints on Carbonate Formation

There have been presented several intriguing observations suggesting evidence of fossil life in martian orthopyroxenite ALH 84001. These exciting and controversial observations have stimulated extensive debate over the origin and history of ALH 84001, but many issues still remain unresolved. Among the most important is the question of the temperature at which the carbonates, which host the putative microfossils, formed. Oxygen- isotopic data, while showing that the carbonates are generally out of isotopic equilibria with the host rock, cannot constrain their temperature of formation. Both low- and high-temperature scenarios are plausible depending on whether carbonate growth occurred in an open or closed system. Petrographic arguments have generally been used to support a high-temperature origin but these appear to be suspect because they assume equilibrium between carbonate compositions that are not in contact. Some observations appear to be consistent with shock mobilization and growth from immiscible silicate-carbonate melts at high temperatures. Proponents of a low-temperature origin for the carbonates are hampered by the fact that there is currently no evidence of hydrous phases that would indicate low temperatures and the presence of a hydrous fluid during the formation of the carbonates. However, the absence of hydrous phases does not rule out carbonate formation at low temperatures, because the carbonate forming fluids may have been extremely CO2 rich, such that hydrous phases would not have been stabilized. In this study, I have carried out additional Transmission electron microscopy (TEM) studies of ALH-84001 and have found evidence of very rare phyllosilicates, which appear to be convincingly of pre-terrestrial origin. At present these observations are limited to one occurrence: further studies are in progress to determine if the phyllosilicates are more widespread.

Brearley, A. J.↗

DOE OSTI2025

Simulations of saturated MHD activity in the HBT-EP tokamak

The NIMROD code is used to perform simulations of MHD activity in the HBT-EP tokamak, including the effect of a resistive wall. Linear simulations are conducted to scan current and pressure profiles to establish self-consistent baseline equilibria that result in MHD instabilities with low error between experimental diagnostics and synthetic equivalents. A family of such equilibria is found with relatively peaked current and pressure profiles. Nonlinear simulations are performed using the optimal equilibrium, which is linearly unstable to a Resistive Wall Mode—stable in the ideal wall limit and with an Alfvénic growth rate in the no-wall limit. Using resistive wall boundary conditions, nonlinear Hall MHD simulations yield a rotating, saturated n=1 mode similar to the type observed experimentally in HBT-EP discharges. An n=1 mode around the q=2 surface mediates a localized axisymmetric perturbation that flattens the plasma current from a linearly unstable state to a 3D, rotating, stable state. During the saturated period, periodic changes in inter-diagnostic phase lag allow for qualitative estimations of effective, macroscopic transport coefficients used in MHD closure models and could be extended to model post-disruption MHD structure rotation after thermal quenches seen in HBT-EP.

Arnold, David A. (ORCID:0009000059535761)↗

NASA NTRS1988

Iron-base superalloys - A phase analysis of the multicomponent system (Fe-Mn-Cr-Mo-Nb-Al-Si-C)

In the course of studies on the iron-rich multicomponent system Fe-Mn-Cr-Mo-Nb-Al-Si-C, work was concentrated on pertinent quinary and six-component combinations namely Fe-Mn-Al-Si-C, Fe-Cr-Al-Si-C and Fe-Mn-Cr-Al-Si-C which had been elaborated at 65, 72, and 80 wt pct Fe. Manganese acts as a strong stabilizer for the cementite carbide. Chromium seems to stabilize the iron aluminide Fe2Al5 which forms in a considerable amount within an alloy of nominal composition Fe(65)Mn(15)Cr(12)Al(5)Si(2)C(1) (percent by weight). Although the Mn3AlC carbide is, like Fe3AlC, a perovskite carbide, manganese does not appear to favor the formation of the perovskite carbide. Because of the relatively low sintering temperature (700 C), for al large portion of the samples equilibria conditions are not always reached.

Gupta, H.↗

DOE OSTI2024

Non-Electricity Based Renewable Fuels: Theory and Computation for Solar Thermochemical Hydrogen

Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox active defects start to interact with each other, thereby impeding the formation of additional defects. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.

first-principles calculations↗

NASA NTRS1993

Molecular equilibria and condensation sequences in carbon rich gases

Chemical equilibria in stellar atmospheres have been investigated by many authors. Lattimer, Schramm, and Grossman presented calculations in both O rich and C rich environments and predicted possible presolar condensates. A recent paper by Cherchneff and Barker considered a C rich composition with PAH's included in the calculations. However, the condensation sequences of C bearing species have not been investigated in detail. In a carbon rich gas surrounding an AGB star, it is often assumed that graphite (or diamond) condenses out before TiC and SiC. However, Lattimer et al. found some conditions under which TiC condenses before graphite. We have performed molecular equilibrium calculations to establish the stability fields of C(s), TiC(s), and SiC(s) and other high temperature phases under conditions of different pressures and C/O. The preserved presolar interstellar dust grains so far discovered in meteorites are graphite, diamond, SiC, TiC, and possibly Al2O3.

Sharp, C. M.↗

NASA NTRS1988

The role of boron and fluids in high temperature, shallow level metamorphism of the Chugach Metamorphic Complex, Alaska

The possible role of boron (B) involvement in granite equilibria and generation of melts during crustal metamorphism has been a focus of speculation in recent literature. Most of the evidence for such involvement derives from experimental data which implies that the addition of B will lower the temperature of the granite solidus. Also the presence of tourmaline has a minor effect on the temperature of the solidus. Further indirect evidence that B may be involved in partial melting processes is the observation that granulites are commonly depleted in B, whereas the B content of low grade metapelites can be high (up to 2000 ppm). Researchers' measurements of the whole-rock B contents of granulites from the Madras region, India are low, ranging from 0.4 to 2.6 ppm. Ahmad and Wilson suggest that B was mobilized in the fluid phase during granulite facies metamorphism of the Broken Hill Complex, Australia. Thus, it appears that during the amphibolite to granulite transition, B is systematically lost from metasediments. The B that is released will probably partition into the vapor phase and/or melt phase. Preliminary measurements imply that the boron content of rocks in the Chugach Metamorphic Complex is not sufficient to influence the processes of partial melting at low pressures.

Sisson, V. B.↗

DOE OSTI2025

Integrated modeling of advanced tokamaks for tearing mode avoidance on DIII-D

Understanding the physics of the onset and evolution of tearing modes (TMs) in tokamak plasmas is crucial for high-performance steady-state operations. The onset condition of (m,n) = (2,1) tearing stability, $Δ′$ > $Δ_c^′$> 0, has been studied in DIII-D steady-state hybrid discharges with accurate equilibrium reconstruction and well-measured plasma profiles. Here, m is the poloidal mode number and n is the toroidal mode number. The onset of n = 1 TMs is carefully identified using the Mirnov array to determine in-phase signals for the onset time and using electron cyclotron emission measurements to trace back from the saturated island width in the experiments for both initial island width and its growth rate. The tearing stability index $Δ′$ at TM onset is calculated in multiple ways, using the modified Rutherford equation and using the PEST-III or resistive DCON magnetohydrodynamic codes, and the results are in reasonable agreement with each other. The calculated $Δ′$ is compared with an analytical formula of the tearing stability threshold $Δ_c^′$ to determine the condition for the mode onset, $Δ′$ > $Δ_c^′$. Finally, the local gradient effect of the plasma current and the pressure profiles on the tearing stability is investigated in terms of the difference between $Δ′$ and $Δ_c^′$ through a series of modeled equilibria with the parametric variations. A database of experimental mode instability onsets is compatible with the simulated regions of instability obtained from the parametric variations.

Linear stability analysis↗

DOE OSTI2024

Thermochemical Modeling of Radionuclide Vapor-Liquid Equilibria in Sodium Pools for SFR Mechanistic Source Term Analysis

Mechanistic source term (MST) analysis of sodium fast reactors (SFR) requires understanding of various radionuclide (RN) transport phenomena influencing potential releases from the fuel to the environment. One such phenomenon includes the retention or release of RNs from the sodium coolant pool, representing the step after possible fuel failures and influencing transport to the cover gas region. Thermodynamic vapor-liquid equilibria (VLE) calculations were performed on systems representing SFR sodium pools containing oxygen impurities and radionuclide (RN) inventories. First, an assessment and recreation of a previously developed thermodynamic database was completed, including updates to thermodynamic parameters. The RN inventories used in VLE calculations represented hypothetical source terms that might be released to the pool during previously analyzed fuel failure scenarios. The calculations were performed for all possible combinations of sodium pool size (i.e., total oxygen) and number of failed fuel pins (i.e., total RNs). In this way, multiple ratios of the RN relative to the oxygen impurity (RN:O) were compared for their impact on RN volatility, which is discussed in terms of the vapor fraction (VF), defined as the fraction of the RN that is calculated to exist in the vapor phase at equilibrium above condensed phases of that element. Similar trends in VLE behavior are seen in the equilibrium calculation results for elements of similar chemistry, and for some element types, it was found that the RN:O ratio can be important due to oxide formation, which typically exist in the condensed phase.

Shahbazi, Shayan↗

DOE OSTI2024

Variable-spectrum mode control of high poloidal beta discharges

DIII-D experiments demonstrate that high pressure, broad current profile equilibria can be accessed in the high poloidal beta regime by optimizing the MHD mode control poloidal spectrum. A novel, variable spectrum (VS) magnetic feedback scheme implemented using the DIII-D internal non-axisymmetric coils (I-coils) facilitated access to reduced internal inductance $l$ i operation above the no-wall beta limit compared with both no feedback and fixed spectrum feedback. In addition, the VS feedback helped avoid beta collapses caused by marginally unstable resistive wall mode activity. The lower and upper I-coil rows were configured in two independent feedback loops, allowing the feedback field's poloidal spectrum to vary and track changes in the plasma mode structure as the edge safety factor q 95 varied from 11 to 6 during the discharges. The q 95 dependence of the measured phase difference between the lower and upper I-coil rows during VS feedback is qualitatively compatible with ideal MHD simulations of the least-stable plasma kink mode and with plasma response simulations that included kinetic modifications to ideal MHD. The VS feedback approach is a straightforward way to improve resilience to variations in mode structure that occur as plasma parameters change. The demonstrated expansion of the operating space to lower $l$ i is expected to improve the coupling of the plasma kink mode to external fields and beneficial wall eddy currents, and is compatible with high bootstrap fraction operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗