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.

337 records · Page 19

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↗

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

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

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↗

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↗

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↗

Two-fluid and kinetic transport physics of Kelvin–Helmholtz instabilities in nonuniform low-beta plasmas

Hall-magnetohydrodynamic (Hall-MHD) theory, two-fluid simulations, and kinetic simulations are used in this work to investigate the cross-field transport properties of Kelvin–Helmholtz instabilities in nonuniform low-beta collisionless plasmas. Hall-MHD analysis shows how the linear properties of the instability are modified by density gradients and magnetization. High-order accurate two-fluid and kinetic simulations, with complete dynamics of finite-mass electrons and ions, are applied to a suite of parameter cases to systematically assess the effects of diamagnetic drift, magnetization, charge separation, and finite Larmor motion. Initialization of exact two-species kinetic equilibria facilitates the study of isolated physical effects and enables detailed cross-comparisons between two-fluid and kinetic simulations, including for cases where ion gyroradii are comparable to gradient scale lengths. For nonuniform plasmas with significant space charge, the results of two-fluid and kinetic simulations are found to disagree with Hall-MHD predictions. Kelvin–Helmholtz instability growth rates, per unit shear, are shown to be smaller when ion diamagnetic drift and $E \times B$ drift are parallel and larger when the two drifts are antiparallel. The effect is attributed to polarization drift in the shear layer, which leads to redistribution of charge, alters the electric field that drives plasma advection, and consequently modifies growth rates. Instability-induced mass transport for different parameters is characterized in terms of the flux across the shear layer and a simplified diffusion model. Distribution functions from kinetic simulations are shown to deviate substantially from Maxwellian reconstructions, indicating the importance of kinetic physics during the nonlinear phase of the instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Phosphate adsorption kinetics and equilibria on natural iron and manganese oxide composites

We report that it is well known that phosphate retention in soils and sediments is strongly influenced by binding to secondary iron oxides, there have been relatively few studies examining its adsorption/desorption behavior on multicomponent particles of realistic natural complexity. In this study, natural Mn-rich limonite (LM), was used to prepare naturally complex Fe- and Mn-oxide composite materials to examine phosphate adsorption/desorption. To clarify the role of the Mn-oxides, results for the LM sample were compared to those for an acid treated version (LAT), in which the acid-extractable Mn-oxide fraction has been selectively eliminated while leaving the Fe-oxide fraction intact. The saturated adsorption capacity on LAT was almost double that on LM, suggesting that phosphate adsorption to the iron oxides is strongly occluded by the Mn-oxide fraction. This result is reinforced by the comparing the pH dependence and fits to adsorption isotherms, and by desorption experiments and STEM-EDS mapping showing that phosphate loading on Mn-oxides was limited. Hence, although the collective results confirm that phosphate uptake and strong binding is selectively controlled by the Fe-oxide fraction, our study reveals that the Mn-oxide fraction strongly interferes with this process. Therefore, phosphate uptake behavior on metal oxides cannot be predicted solely on the basis of the Fe-oxide fraction present, but instead must take into account the deleterious impacts of other intimately associated phases. For co-diagenetic Fe/Mn-oxide composites in particular, Mn-oxides appear to severely limit phosphate uptake on the Fe-oxide fraction, either by hindering access to binding sites on the Fe-oxide or by lowering their affinity for P.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Detailed Compositional Comparison of Hydrogenated Vegetable Oil with Several Diesel Fuels and Their Effects on Engine-Out Emissions

The Coordinating Research Council (CRC) is actively involved in developing and applying advanced analytical techniques to the chemical characterization of transportation fuels. Here, this article complements a 2017 CRC project to quantify and compare the effects of a commercially available renewable diesel fuel (hydrotreated vegetable oil [HVO]) and an ultralow sulfur diesel (ULSD) fuel on engine-out gaseous and particulate matter (PM) emissions from a light-duty vehicle. Results showed that the combustion of HVO fuel had an advantage over ULSD in terms of lowering engine-out emissions (THC, CO, NO x , etc.). Furthermore, this advantage is strongly related to the fuel composition. This article summarizes the results of advanced and comprehensive analytical tests on the same ULSD and HVO fuels and attempts to connect some of the engine-out emissions results to fuel composition and specific chemical structures. A variety of test methods, generally unavailable in combination, were employed, such as one-dimensional (1D) and two-dimensional (2D) gas chromatography (GC), nuclear magnetic resonance spectroscopy (NMR), and high-pressure solid-liquid phase transition experiments. In summary, the ULSD sample was found to have representation across the expected set of hydrocarbon classes typical for the sample type. Interestingly, a high content of cycloparaffins (>50 wt%) and a very low content of diaromatics (~2 wt%) were present. While not without precedent, these are higher and lower, respectively, than typically found for commercial ULSD compositions. In contrast, HVO was found to consist of only two hydrocarbon classes: n-paraffins (~10 wt%) and iso-paraffins (~90 wt%), both predominantly in a narrow carbon atom number range (i.e., C14–C18). HVO engine-out emissions results for the LA-92 and steady-state testing can be tracked to the narrow carbon atom number range of the n-paraffins and iso-paraffins, which result in a high cetane number fuel having a narrow distillation range. Previously, the low-temperature operability of HVO has been a concern, but that appears not be the case for this particular HVO. HVO and ULSD were evaluated at pressures up to ~275 MPa and found to have comparable solid-liquid equilibria despite significant compositional differences.

09 BIOMASS FUELS↗

An abiotic origin for hydrocarbons in the Allan Hills 84001 martian meteorite through cooling of magmatic and impact-generated gases

Thermodynamic calculations of metastable equilibria were used to evaluate the potential for abiotic synthesis of aliphatic and polycyclic aromatic hydrocarbons (PAHs) in the martian meteorite Allan Hills (ALH) 84001. The calculations show that PAHs and normal alkanes could form metastably from CO, CO2, and H2 below approximately 250-300 degrees C during rapid cooling of trapped magmatic or impact-generated gases. Depending on temperature, bulk composition, and oxidation-reduction conditions, PAHs and normal alkanes can form simultaneously or separately. Moreover, PAHs can form at lower H/C ratios, higher CO/CO2 ratios, and higher temperatures than normal alkanes. Dry conditions with H/C ratios less than approximately 0.01-0.001 together with high CO/CO2 ratios also favor the formation of unalkylated PAHs. The observed abundance of PAHs, their low alkylation, and a variable but high aromatic to aliphatic ratio in ALH 84001 all correspond to low H/C and high CO/CO2 ratios in magmatic and impact gases and can be used to deduce spatial variations of these ratios. Some hydrocarbons could have been formed from trapped magmatic gases, especially if the cooling was fast enough to prevent reequilibration. We propose that subsequent impact heating(s) in ALH 84001 could have led to dissociation of ferrous carbonates to yield fine-grain magnetite, formation of a CO-rich local gas phase, reduction of water vapor to H2, reequilibration of the trapped magmatic gases, aromatization of hydrocarbons formed previously, and overprinting of the synthesis from magmatic gases, if any. Rapid cooling and high-temperature quenching of CO-, H2-rich impact gases could have led to magnetite-catalyzed hydrocarbon synthesis.

Non-NASA Center↗

The life cycle of the Interstellar Medium in other galaxies

Gas in spiral galaxies cycles between the diffuse and dense phases as clouds collapse, form stars and are dispersed back into the ISM. Far infrared observations of continuum emission from interstellar dust and line emission from interstellar gas have revealed a wealth of information on the state of the ISM in galaxies of different morphological types. The analysis of these observations gives us information about the processes of star formation and about the evolution of the ISM. Star formation rates vary widely from galaxy to galaxy, with the rates in starburst galaxies being 10 - 100 times those in quiescent spiral galaxies. Far infrared spectroscopy of star-forming galaxies shows that the interstellar pressure increases with star formation rate. The structure of the interstellar medium in starburst galaxies is quite different from that of quiescent galaxies - much of the mass and volume are in HII regions and photodissociation regions. The size distribution of dust grains seems to depend on environment; small grains are abundant in the diffuse interstellar medium but not in dense molecular star forming regions. Quiescent spiral and elliptical galaxies contain a significant population of small grains, but starburst galaxies do not. Dwarf irregular galaxies also seem to contain few small grains; this may be the result of the higher UV flux in these galaxies. The star forming regions in dwarf irregulars also have a higher ratio of atomic to molecular gas than do those in the Galaxy. These results show that the ISM in galaxies of different morphological types reaches different equilibria, resulting in different modes of star formation and global galaxy evolution.

Knapp, G. R.↗

Analysis of Direct Samples of Early Solar System Aqueous Fluids

Over the past three decades we have become increasingly aware of the fundamental importance of water, and aqueous alteration, on primitive solar-system bodies. Some carbonaceous and ordinary chondrites have been altered by interactions with liquid water within the first 10 million years after formation of their parent asteroids. Millimeter to centimeter-sized aggregates of purple halite containing aqueous fluid inclusions were found in the matrix of two freshly-fallen brecciated H chondrite falls, Monahans (1998, hereafter simply "Monahans") (H5) and Zag (H3-6) (Zolensky et al., 1999; Whitby et al., 2000; Bogard et al., 2001) In order to understand origin and evolution of the aqueous fluids inside these inclusions we much measure the actual fluid composition, and also learn the O and H isotopic composition of the water. It has taken a decade for laboratory analytical techniques to catch up to these particular nanomole-sized aqueous samples. We have recently been successful in (1) measuring the isotopic composition of H and O in the water in a few fluid inclusions from the Zag and Monahans halite, (2) mineralogical characterization of the solid mineral phases associated with the aqueous fluids within the halite, and (3) the first minor element analyses of the fluid itself. A Cameca ims-1270 equipped with a cryo-sample-stage of Hokkaido University was specially prepared for the O and H isotopic measurements. The cryo-sample-stage (Techno. I. S. Corp.) was cooled down to c.a. -190 C using liquid nitrogen at which the aqueous fluid in inclusions was frozen. We excavated the salt crystal surfaces to expose the frozen fluids using a 15 keV Cs+ beam and measured negative secondary ions. The secondary ions from deep craters of approximately 10 m in depth emitted stably but the intensities changed gradually during measurement cycles because of shifting states of charge compensation, resulting in rather poor reproducibility of multiple measurements of standard fluid inclusions of +/- 90 0/00(2 sigma) for delta D, and +/- 29 0/00 (2 sigma) for delta O-18. On the other hand, the reproducibility of Delta O-17 is plus or minus 8 /00 (2 sigma ) because the observed variations of isotope ratios follow a mass dependent fractionation law. Variations of delta D of the aqueous fluids range over sog,a 330(90; 2 sigma ) to +1200(90) 0/00 for Monahans and delta 300(96) 0/00 to +90(98)0/00 for Zag. Delta O-17 of aqueous fluids range over delta 16(22) 0/00 to +18(10) 0/00 for Monahans and +3(10) 0/00 to +27(11) 0/00 for Zag. These variations are larger than the reproducibility of standard analyses and suggest that isotope equilibria were under way in the fluids before trapping into halite. The mean values of delta D and Delta O-17 are +290 0/00 and +9 0/00, respectively. The mean values and the variations of these fluids are different from the representative values of ordinary chondrites, verifying our working hypothesis that the fluid inclusion-bearing halites were not indigenous to the H chondrite parent-asteroid but rather represent exogenous material delivered onto the asteroid from a separate cryovolcanically-active body. This initial isotopic work has demonstrated the feasibility of the measurements, but also revealed sample processing and analytical shortcomings that are now being addressed. Examination of solid mineral inclusions within Monahans and Zag halite grains by confocal Raman spectroscopy at the Carnegie Geophysical Laboratory has revealed them to be metal, magnetite, forsteritic olivine (Fo.98), macromolecular carbon (MMC), pyroxenes, feldspar with Raman spectral affinity to anorthoclase and, probably, fine-grained lepidocrocite (FeO(OH)). In addition, one inclusion features aliphatic material with Raman spectral features consistent with a mixture of short-chain aliphatic compounds. We have initiated analyses of the bulk composition of the fluids within the inclusions in Zag and Monahans halites at Virginia Tech by LA ICPMS using angilent 7500ce quadrupole ICPMS and a Lambda Physik GeoLas 193 nm excimer laser ablation system. Preliminary results reveal that the inclusion aqueous fluids contain highly charged cations of Ca, Mg and Fe. The minerals and compounds discovered thus far within Monahans/Zag halites are indicative of an originating body at least partly composed of unequilibrated anhydrous materials (high Fo olivine, pyroxenes, feldspars, possibly the metal) which were subjected to aqueous alteration (the halite parent brine) and containing a light organic component (the short-chain aliphatic compounds). This material was ejected from the originating body with little or no disruption, as evidenced with the presence of fluid inclusions. An actively geysering body similar to modern Enceladus (Postberg et al., 2011) may be a reasonable analogue in this respect. Also, the originating body should have been within close proximity to the H chondrite parent in order to generate the number of halite grains seen in Monahans and Zag. Other candidates for Monahans/Zag halite parent bodie(s) may include a young Ceres with its possible liquid ocean, or Main Belt comets.

Zolensky, Michael E.↗

Corrosion of 316 SS in direct contact with Zircaloy-4 in repository relevant conditions

The effects of interfacial contact between 316 stainless steel (316 SS) waste canister and spent fuel cladding material (Zircaloy-4) in case of water intrusion into a deep geological repository are investigated. Hydrothermal corrosion experiments of 316 SS and Zircaloy-4 in direct contact were performed at elevated temperatures in MilliQ water and simulated Mont Terri groundwater. Hematite, magnetite, chromite, and baddeleyite form as secondary phases. pH of the Mont Terri solution increases with corrosion duration and magnetite forms preferentially instead of hematite. Thermodynamic modeling shows that iron oxidation in 316 SS controls the pH and redox potential of the aqueous environment and is in good agreement with experimental Mont Terri results, showing that experimental Mont Terri results trend to equilibrium. All oxides except hematite are predicted to exist at thermodynamic equilibrium based on calculations, implying that experimentally formed hematite will be reduced or transformed to magnetite at equilibrium. The thermodynamically predicted chromite was observed to form only at the interface where 316 SS was in contact with Zircaloy-4 at 250°C, showing that contact with Zircaloy-4 seems to accelerate the 316 SS corrosion to thermodynamic equilibrium. Finally, Mont Terri groundwater readily trends towards thermodynamic equilibrium, leading to a pH and redox buffered system that is dominated by iron in solution coming from the accelerated corrosion of 316 SS caused by Zircaloy-4 contact.

316 stainless steel↗