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 163 records · Page 9

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

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

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↗

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

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↗

Theory of Dust Voids in Plasmas

Dusty plasmas in a gas discharge often feature a stable void, i.e., a dust-free region inside the dust cloud. This occurs under conditions relevant to both plasma processing discharges and plasma crystal experiments. The void results from a balance of the electrostatic and ion drag forces on a dust particle. The ion drag force is driven by a flow of ions outward from an ionization source and toward the surrounding dust cloud, which has a negative space charge. In equilibrium the force balance for dust particles requires that the boundary with the dust cloud be sharp, provided that the particles are cold and monodispersive. Numerical solutions of the one-dimensional nonlinear fluid equations are carried out including dust charging and dust-neutral collisions, but not ion-neutral collisions. The regions of parameter space that allow stable void equilibria are identified. There is a minimum ionization rate that can sustain a void. Spatial profiles of plasma parameters in the void are reported. In the absence of ion-neutral collisions, the ion flow enters the dust cloud's edge at Mach number M = 1. Phase diagrams for expanding or contracting voids reveal a stationary point corresponding to a single stable equilibrium void size, provided the ionization rate is constant. Large voids contract and small voids expand until they attain this stationary void size. On the other hand, if the ionization rate is not constant, the void size can oscillate. Results are compared to recent laboratory and microgravity experiments.

Goree, J.↗

Iron Redox Systematics in Martian Mantle Melts and Identification of Primary Liquids

Many MgO-rich shergottites have olivine phenocrysts and a natural question to ask is if any of these represent primary melts from the martian mantle. To answer this requires knowledge of olivine liquid equilibrium in appropriate magma compositions. Because olivine-phyric shergottites represent either primary liquids or liquids that have been modified by fractionation or accumulation of olivine (or other phases), the equilibrium (and its Kd) 2MgO (liq) + Fe2SiO4 (olivine) = 2FeO (liq) + Mg2SiO4 (olivine), can be used to evaluate whether the olivine-bearing rocks represent liquids.

olivine↗

Electrochemical measurements and thermodynamic calculations of redox equilibria in pallasite meteorites - Implications for the eucrite parent body

The intrinsic oxygen fugacity (IOF) of olivine separates from the Salta, Springwater, and Eagle Station pallasites was measured between 850 and 1150 C using oxygen-specific solid zirconia electrolytes at 100,000 Pa. Thermodynamic calculations of redox equilibria involving equalibrium pallasite assemblages are in good agreement with the experimental results and provide a lower limit to pallasite redox stability; others involving disequilibrium assemblages, suggest that pallasites experienced localized, late-stage oxidation and reduction effects. Consideration of the redox buffer metal-olivine-orthopyroxene utilizing calculated Eucrite Parent Body (EPB) mantle phase compositions indicates that small redox gradients may have existed in the EPB. Such gradients may have produced strong compositional variation within the EPB. In addition, there is apparently significant redox heterogeneity in the source area of Eagle Station Trio pallasites and Bocaiuva iron meteorites.

Righter, Kevin↗

Magnetotail dynamics under isobaric constraints

Using linear theory and nonlinear MHD simulations, we investigate the resistive and ideal MHD stability of two-dimensional plasma configurations under the isobaric constraint dP/dt = 0, which in ideal MHD is equivalent to conserving the pressure function P = P(A), where A denotes the magnetic flux. This constraint is satisfied for incompressible modes, such as Alfven waves, and for systems undergoing energy losses. The linear stability analysis leads to a Schroedinger equation, which can be investigated by standard quantum mechanics procedures. We present an application to a typical stretched magnetotail configuration. For a one-dimensional sheet equilibrium characteristic properties of tearing instability are rediscovered. However, the maximum growth rate scales with the 1/7 power of the resistivity, which implies much faster growth than for the standard tearing mode (assuming that the resistivity is small). The same basic eigen-mode is found also for weakly two-dimensional equilibria, even in the ideal MHD limit. In this case the growth rate scales with the 1/4 power of the normal magnetic field. The results of the linear stability analysis are confirmed qualitatively by nonlinear dynamic MHD simulations. These results suggest the interesting possibility that substorm onset, or the thinning in the late growth phase, is caused by the release of a thermodynamic constraint without the (immediate) necessity of releasing the ideal MHD constraint. In the nonlinear regime the resistive and ideal developments differ in that the ideal mode does not lead to neutral line formation without the further release of the ideal MHD constraint; instead a thin current sheet forms. The isobaric constraint is critically discussed. Under perhaps more realistic adiabatic conditions the ideal mode appears to be stable but could be driven by external perturbations and thus generate the thin current sheet in the late growth phase, before a nonideal instability sets in.

Birn, Joachim↗

Intraseasonal variability in a barotropic model with seasonal forcing

It has recently been suggested that oscillatory topographic instability could contribute to low-frequency variability over the Northern Hemisphere midlatitudes. A barotropic potential vorticity model, with a hierarchy of forcing and topography configurations on the sphere, is used to investigate the nature of low-frequency oscillations induced by such instabilities. Steady-state solutions of the model include multiple unstable equilibria that sustain oscillatory instabilities with periods of 10 - 15 days, and 150 - 180 days, for a realistic forcing pattern, Time-dependent solutions exhibit chaotic behavior with episodic oscillations, featuring both the intraseasonal (35 - 50 day) and biweekly (10 - 15 day) modes. The former is dominated by standing spatial patterns, the latter by traveling wave patterns. The phases of the intraseasonal oscillation are robust for all cases, exhibiting a clear oscillatory exchange of atmospheric angular momentum with the solid earth via mountain troque. It is demonstrated, through linear stability analysis on the sphere, that the intraseasonal oscillations are induced by topographic instabilities. The role of the seasonal cycle is studied by prescribing an annual cycle in the forcing. In this case, the winter forcing is more favorable than the summer for the occurrence of episodic intraseasonal oscillations. Recent observations are consistent with this model result.

Strong, Christopher M.↗

High pressure cosmochemistry applied to major planetary interiors: Experimental studies

The measurement of equilibria in binary fluid-solid systems in diamond anvil cells, represents a major advance of the art of high-pressure experimentation. Vibrational spectroscopy, direct visual observations, and X-ray diffraction crystallography of materials confined in externally heated cells are the primary experimental probes being used. Adiabats in these systems are being measured in order to constrain models of heat flow in these bodies and to detect phase transitions by thermal anomalies. Other studies are directed toward interpreting high pressure reactions in these systems that are suggested by shockwave measurements, and developing methods for reaching high temperatures and high pressures of planetary interest in diamond cells. The overall objective of this project is to determine the properties of the H2-He-H2O-HN3-CH4 system and related small-molecule systems that are needed to constrain theoretical models of the interiors of the major planets.

Nicol, M. F.↗

Iron Redox Systematics of Martian Magmas

Martian magmas are known to be FeO-rich and the dominant FeO-bearing mineral at many sites visited by the Mars Exploration rovers (MER) is magnetite [1]. Morris et al. [1] propose that the magnetite appears to be igneous in origin, rather than of secondary origin. However, magnetite is not typically found in experimental studies of martian magmatic rocks [2,3]. Magnetite stability in terrestrial magmas is well understood, as are the stability of FeO and Fe2O3 in terrestrial magmas [4,5]. In order to better understand the variation of FeO and Fe2O3, and the stability of magnetite (and other FeO-bearing phases) in martian magmas we have undertaken an experimental study with two emphases. First we document the stability of magnetite with temperature and fO2 in a shergottite bulk composition. Second, we determine the FeO and Fe2O3 contents of the same shergottite bulk composition at 1 bar and variable fO2 at 1250 C, and at variable pressure. These two goals will help define not only magnetite stability, but pyroxene-melt equilibria that are also dependent upon fO2.

Righter, K.↗

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

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↗

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

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

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

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↗