Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gas fugacity”

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.

The Reservoir Temperature Estimator (RTEst): A multicomponent geothermometry tool

The Reservoir Temperature Estimator (RTEst) is a multicomponent geothermometry tool for estimating reservoir geochemical parameters including reservoir temperature, CO 2 fugacity, mass of water lost or gained, and a reaction factor. It estimates these parameters and their associated uncertainties by minimizing an objective function that is the weighted sum of squares of the saturation indices of a user-selected set of minerals believed to be equilibrated with the reservoir fluid. RTEst accomplishes these estimates by combining the geochemical modeling capabilities of The Geochemist’s Workbench® with the optimization/parameter estimation resources of PEST®. An included interface aids the user in selecting plausible mineral phases to comprise the objective function and calculates their weighting factors. The working principles of RTEst are described and its efficacy is illustrated by presenting results of its application to various geothermal fields with known conditions. These examples show RTEst can account for the alteration of ascending reservoir fluid by mineral (calcite) re-equilibration with changes in temperature, reconstruct waters with CO 2 loss, correct for the deficit of water and other volatiles (CO 2,gas , H 2 S gas ) from boiling, and determine the amount of mixing of thermal and non-thermal waters. RTEst can use data with basis species below detection limit, missing, or unreliable either by assuming equilibrium with a controlling mineral (fixed-analyte method) or by treating the analyte concentration as an optimization parameter. The inverse of variance weighting method included in RTEst provides more representative results than either the normalization or unit weighting methods. Finally, the ability of RTEst to calculate reservoir temperatures, gas fugacity, and mixing fractions demonstrates its usefulness as a tool for evaluating geothermal systems.

15 GEOTHERMAL ENERGY↗

Community Based Data of Uranium Adsorption onto Quartz

This data upload includes a compilation of experiments for quantifying uranium adsorption onto quartz. The provided .csv file has been compiled from the literature in a findable, accessible, interoperable, reusable (FAIR) data format. This was accomplished using the Lawrence Livermore National Laboratory Surface Complexation Database Converter (SCDC) code written in the R programming language (free licensing available at https://ipo.llnl.gov/technologies/software/llnl-surface-complexation-database-converter-scdc). This constitutes all current data compiled on uranium-quartz interactions in the L-SCIE (LLNL Surface Complexation/Ion Exchange) database (as of 02/28/2022). This data was used to develop surface complexation models that fit the global community dataset (https://doi.org/10.1021/acs.est.1c07109). The FAIR-formatted dataset also enables the implementation of alternative machine-learning approaches that can be explored in the future.

54 ENVIRONMENTAL SCIENCES↗

Structure factors of the unitary gas under supernova conditions

Here, we compute with lattice field theory the vector and axial static structure factors of the unitary gas for arbitrary temperature above the superfluid transition and for fugacities 0.1 < $\mathcal{z}$ < 1.0 . Using the lattice formulation, we calculate beyond the validity of the virial expansion, a commonly used technique in many-body physics. We find qualitative differences in the behavior of the structure factors at high fugacity compared to the predictions of the virial expansion. Due to the large scattering length of neutrons, we expect the unitary gas structure factors to approximate the structure factors of hot neutron gases, and we therefore expect our calculations to be useful in supernova simulations, where neutron gas structure factors are needed to compute in-medium neutrino-neutron scattering rates.

79 ASTRONOMY AND ASTROPHYSICS↗

Pure- and mixed-gas transport properties of a microporous Tröger's Base polymer (PIM-EA-TB)

Polymers of intrinsic microporosity (PIMs) offer tantalizing combinations of high selectivity and permeability in initial gas permeation measurements. Here, we report characterization of pure- and mixed-gas permeation properties of a thick (~80 μm) PIM consisting of Tröger's Base (TB) and Ethanoanthracene (EA) films. The effects of feed pressure and temperature on pure-gas permeabilities of CH 4 , N 2 , O 2 , H 2 , and CO 2 were investigated. The physical aging behavior of the thick film was tracked via pure-gas O 2 , N 2 , and CH 4 permeability at 35 °C. Gas permeability decreased noticeably and selectivity increased as aging time increased. Particular attention was given to mixed-gas measurements of CO 2 and CH 4 (50/50) permeabilities at 35 °C and fugacities ranging from 2 to 18 atm to explore whether the rigid, bridged, bicyclic TB and EA units could resist CO 2 -induced plasticization. These results are presented along with pure-gas CO 2 and CH 4 results for membrane samples aged at different times. PIM-EA-TB aged for ~24 h did not show signs of plasticization across the fugacity range considered. Furthermore, dual-mode competitive sorption presumably caused the CO 2 /CH 4 mixed-gas selectivity to be slightly higher than its corresponding pure-gas selectivity. However, as aging time increased, aged films underwent progressively more rapid and extensive CO 2 -induced plasticization with increasing fugacity, suggesting a systematic relationship between physical aging and plasticization in PIMs. Consequently, physical aging caused less improvement in mixed-gas CO 2 /CH 4 selectivity than it did in pure-gas selectivity, due mainly to plasticization effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Asymptotic Freedom at the Berezinskii-Kosterlitz-Thouless Transition without Fine-Tuning Using a Qubit Regularization

We propose a two-dimensional hard-core loop-gas model as a way to regularize the asymptotically free massive continuum quantum field theory that emerges at the Berezinskii-Kosterlitz-Thouless transition. Without fine-tuning, our model can reproduce the universal step-scaling function of the classical lattice X Y model in the massive phase as we approach the phase transition. This is achieved by lowering the fugacity of Fock-vacuum sites in the loop-gas configuration space to zero in the thermodynamic limit. Some of the universal quantities at the Berezinskii-Kosterlitz-Thouless transition show smaller finite size effects in our model as compared to the traditional X Y model. Our model is a prime example of qubit regularization of an asymptotically free massive quantum field theory in Euclidean space-time and helps understand how asymptotic freedom can arise as a relevant perturbation at a decoupled fixed point without fine-tuning. Published by the American Physical Society 2024

Physics↗

Pressure-stable supported ionic liquid membranes using isoporous supports for evaluating pure- and mixed-gas light paraffin fractionation

Advances in horizontal drilling and hydraulic fracturing have spurred the growth of domestic U.S. energy production. Membranes, typically silicone rubbers, have found utility in shale gas treatment as fuel gas conditioning units to selectively remove C 2 + hydrocarbons at pressures up to 30 bar, producing clean CH 4 for gas engines. However, more selective materials could be beneficial for broader shale gas treatment applications, such as dew point control units. Supported ionic liquid membranes (SILMs) offer a potential opportunity for improving C 3 H 8 /CH 4 selectivity, but they lack pressure-stability. Here, we report C 3 H 8 /CH 4 selective and pressure-stable SILMs using isoporous supports. SILMs with supports that had minimal defects remained stable up to 15 bar of transmembrane pressure. This stability allowed for pure- and mixed-gas testing of the resulting membranes at elevated pressures. Furthermore, these tests revealed that low viscosity ILs (<100 cp) may display mixed-gas permeances and permselectivity nearly identical to pure-gas results. On the other hand, higher viscosity ILs may display increasing permeances and permselectivity with increasing C 3 H 8 fugacity, similar to rubbery polymers. Ultimately, the SILMs demonstrated relatively high pressure-stability due to the isoporous supports and competitive mixed-gas C 3 H 8 /CH 4 permselectivity compared to silicone rubber.

Rosenthal, Justin J. [University of Texas at Austi↗

Synthesis and Properties of the Helium Clathrate and Defect Perovskite [He 2– x $^\square$ x ][CaNb]F 6

The defect double perovskite [He 2– x $^\square$ x ][CaNb]F 6 , with helium on its A-site, can be prepared by the insertion of helium into ReO 3 -type CaNbF 6 at high pressure. Upon cooling from 300 to 100 K under 0.4 GPa helium, ~60% of the A-sites become occupied. Helium uptake was quantified by both neutron powder diffraction and gas insertion and release measurements. After the conversion of gauge pressure to fugacity, the uptake of helium by CaNbF 6 can be described by a Langmuir isotherm. The enthalpy of absorption for helium in [He 2– x $^\square$ x ][CaNb]F 6 is estimated to be ~+3(1) kJ mol –1 , implying that its formation is entropically favored. Helium is able to diffuse through the material on a time scale of minutes at temperatures down to ~150 K but is trapped at 100 K and below. The insertion of helium into CaNbF 6 reduces the magnitude of its negative thermal expansion, increases the bulk modulus, and modifies its phase behavior. On compressing pristine CaNbF 6 , at 50 and 100 K, a cubic (Fm$\overline{3}$m) to rhombohedral (R$\overline{3}$) phase transition was observed at <0.20 GPa. However, a helium-containing sample remained cubic at 0.4 GPa and 50 K. CaNbF 6 , compressed in helium at room temperature, remained cubic to >3.7 GPa, the limit of our X-ray diffraction measurements, in contrast to prior reports that upon compression in a nonpenetrating medium, a phase transition is detected at ~0.4 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vapor-Phase Aggregation of Cerium Oxide Nanoparticles in a Rapidly Cooling Plasma

Local conditions, such as temperature and oxygen availability, have a pronounced effect on the formation and evolution of fallout following a nuclear explosion. While the behavior of nuclear-relevant materials such as uranium has begun to be explored under a wider range of environments, little is known about the behavior of plutonium. Here, using cerium as a surrogate, we track the vapor-phase aggregation of cerium oxide nanoparticles created in a plasma flow reactor under conditions of controlled temperature at two different oxygen fugacities. In situ optical emission spectroscopy is used to measure the variations in the spectral intensity of atomic and molecular species with temperature and oxygen content. We find that the relative rate of gas-phase oxidation of cerium is highly dependent on both temperature and local redox conditions within the flow reactor, to the extent that doubling the oxygen availability effectively doubles the amount of vapor-phase cerium monoxide at high temperatures (>2000 K). Condensed cerium oxide nanoparticles are also collected and analyzed ex situ via transmission electron microscopy and grazing-incidence small-angle X-ray scattering to determine their elemental composition, crystal structure, and size distribution. The size and morphology of the condensed nanoparticles are independent of local redox conditions, forming the same crystal type with the same size distribution regardless of oxygen availability. Postcondensation particle evolution, however, is found to be predominantly driven by temperature, with the average particle size increasing as particles cool and subsequently aggregate. These results expand our understanding of the chemical and physical behavior of refractory oxides that form during the early stages of fallout formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental investigation of elemental and isotopic evaporation processes by laser heating in an aerodynamic levitation furnace

In this report we carried out evaporation experiments on a B-type calcium–aluminium-rich inclusion (CAI) melt in a gas-mixing aerodynamic levitation laser furnace, at 1873 K and an oxygen partial pressure of 10 -9.1 atm, for durations ranging from 60 to 600 s. Evaporation of SiO 2 and MgO follow the same trend as those observed in vacuum furnace experiments at the same temperature and starting composition, showing that their evaporation relative to one another from the melt is independent of pressure, oxygen fugacity, and hydrodynamical regime specific to the furnace. Isotopic ratios of Mg and Si in evaporation residues are used to derive fractionation factors of α 26/24 Mg vap-liq = 0.9906± 0.0004 and α 30/28 Si vap-lip = 0.9943±0.0003, which are both significantly closer to unity than those found for evaporation in a vacuum, which translates to less isotope fractionation. The residues are also less isotopically fractionated than expected for cases in which transport of the gas species away from the melt is diffusion-controlled at 1-atm. By analysing the flow regimes in our furnace, we find that advection by the levitating gas is the primary mode of mass transport away from the melt surface, as opposed to diffusion-limited transport in a vacuum or 1-atm tube furnace. A modified Hertz–Knudsen–Langmuir formulation accounts for this process, and shows that isotopic fractionation of both Si and Mg reflect a saturation factor (ratio of the pressure of the evaporating species to vapour saturation pressure) equal to 0.75. This is in perfect accord with recent measurements of Cu isotopic fractionation using a similar furnace. The fact that three elements (Mg, Si, Cu) with varying equilibrium vapour pressures, activity coefficients in the liquid, and diffusion coefficients in the gas have the same scaling behaviour to saturation pressure is a strong indication that the mechanism controlling evaporation is driven by the hydrodynamical regime imposed in the furnace. Therefore, this class of experiments can be used to constrain processes in which advection dominates over diffusion, such as (but not limited to) planetary ejecta, tektites, giant impacts, nebular condensation in a turbulent flow, or nuclear fallout material. Finally, the possibility to reach high temperatures (in excess of 3500 K) in this furnace allows it to be used to evaluate the activity coefficients of melt components in extreme conditions relevant to molten planetary interiors (i.e., magma oceans), with a specific focus on refractory elements.

58 GEOSCIENCES↗

Scale translation yields insights into gas adsorption under nanoconfinement

This work describes a scale-translating simulation framework to investigate gas adsorption behavior in nanoconfined pores. The framework combines molecular simulations (MSs), equation of state (EoS), and lattice Boltzmann (LB) simulations. MSs reveal the physics of methane adsorption in nano-sized pores, where input values of fugacity coefficients are optimized based on EoS predictions. Then, an LB free-energy model, which incorporates a viral EoS, upscales intermolecular forces and estimates adsorption behavior via a proposed fluid–wall interaction model. Armed with the values of the LB interaction parameter as a function of pressure, the LB model is used to predict fluid behavior in irregular nanopores, and the results are validated against reference MS data. The LB model is then used to study adsorption behavior at a continuum scale in representative organic shale nanopores based on finely characterized Vaca Muerta shale samples. Furthermore, the results show that methane adsorption could significantly increase contained fluids by 10%–25% in pores smaller than 20 nm. However, in larger pores (40 nm to 90 nm), adsorption's impact diminishes to 2%–3%, suggesting sorption's negligible role beyond a 40 nm pore size.

74 ATOMIC AND MOLECULAR PHYSICS↗

Thermal relaxation and the complete set of second-order transport coefficients for the unitary Fermi gas from kinetic theory

We compute the complete set of second-order transport coefficients of the unitary Fermi gas, a dilute gas of spin-1/2 particles interacting via an 𝑠 -wave interaction tuned to infinite scattering length. The calculation is based on kinetic theory and the Chapman-Enskog method at second order in the Knudsen expansion. We take into account the exact two-body collision integral. We extend previous results on second-order coefficients related to shear stress by including terms related to heat flow and gradients of the fugacity. We confirm that the thermal relaxation time is given by the simple estimate 𝜏 𝜅 = 𝜅⁢𝑚/(𝑐 𝑃 ⁢𝑇) even if the full collision kernel is taken into account. Furthermore, 𝜅 is the thermal conductivity, 𝑚 is the mass of the particles, 𝑐𝑃 is the specific heat at constant pressure, and 𝑇 is the temperature.

Kinetic theory↗

Bridging adsorption behavior of confined CH 4 -CO 2 binary mixtures across scales

An accurate understanding of the competitive adsorption of CH 4 -CO 2 binary mixtures in nano-confined systems is critical for engineering CO 2 storage in shale gas reservoirs. Due to difficulties in making reliable experimental observations in nano-scale, atomistic simulations (ASs), such as the Grand Canonical Monte Carlo (GCMC) method, provide a viable approach to studying the adsorption behavior of confined fluids. ASs are, however, limited in the size of the compositional domain due to the high computational cost. This work proposes a framework that combines AS and the lattice Boltzmann (LB) method to bridge the physics of confined fluids across scales. The Peng–Robinson equation of state (PR-EoS) produces fugacity coefficients, which serve as input for conducting multi-component GCMC simulations. These GCMC simulations explore the competitive adsorption behavior of CH 4 -CO 2 in nano-slits at various composition, pressure, and channel-width conditions. Both components generate adsorption layers with high densities near the walls with CO 2 preferentially adsorbing compared to CH 4 on the organic walls of carbon sheets. At the mesoscale, a pseudopotential model represents the intermolecular forces in multi-component, multiple-relaxation-time LB simulations. The LB simulations are in good agreement with the GCMC results, allowing us to obtain values for tunable LB parameters. We then extend the use of LB to simulate adsorption behavior in complex networks with nano-sized channels. The phase behavior and fluid properties in the complex geometries of nano-channels differ from nano-slits and bulk systems. Furthermore, the bridging of physics from GCMC (microscale) to LB (mesoscale) via the macroscale PR-EoS connects the adsorption behavior of binary systems across scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The ever-evolving active site: transformation of single atoms to extended structures during the Rh-catalyzed reverse water–gas shift reaction

At low temperatures (< 400°C), single atoms of Rh supported on rutile TiO 2 (rTiO 2 ) are responsible for the formation of CO during the reverse water gas shift (RWGS), while methane production is associated with the Rh-TiO 2 interface due to the correlation between methane formation rates and the volume-averaged Rh nanoparticle diameter. As the temperature is increased to >540°C, there is a notable increase in CO selectivity and methane production rates tend towards zero. The time to reach zero depends on the temperature but independent of the initial Rh structure (single atoms and/or nanoparticles), which is controlled by the catalyst preparation method (wetness impregnation versus colloidal nanoparticles). At 600°C and > 4 h time on stream, the catalytic behaviour becomes completely agnostic to the initial Rh structure as well as weight loading, and the catalysts are highly selective for the reverse water gas-shift reaction. Post-reaction HR-TEM image analysis confirms Rh nanoparticles crystallize/order during the reaction; at 400°C, most of the Rh particles are disordered, while at 600°C, they are more ordered (i.e., development of defined faceting). Infrared spectroscopy of CO adsorption on Rh nanoparticles confirms the appearance of defined facets after annealing in nitrogen at high temperatures. Annealing the Rh/rTiO 2 catalysts prior to the RWGS reaction demonstrates the structural transformation of Rh depends only on time and temperature and not on reactant or product fugacity. Sites responsible for stabilizing Rh single atoms are no longer competent at higher temperatures, enabling their integration into existent nanoparticles. As the reaction temperature is increased to temperatures >540°C, the dominant active site for CO production evolves from single atoms to extended Rh structures.

36 MATERIALS SCIENCE↗

An Experimental Study of the Solubility of Rare Earth Chloride Salts (La, Nd, Er) in HCl Bearing Water Vapor from 350 – 425 °C

In this work, the solubilities of the rare earth chlorides REECl 3 , where REE = (La, Nd, Er), were measured in HCl bearing water vapor from 350 – 425°C with water partial pressures ranging from 8 – 170 bar. Solubility data were fit to the Pitzer-Pabalan quasi-chemical model in order to extract thermodynamic parameters for the formation of the gaseous REE-chloride-water clusters REECl 3 (H 2 O) n . The data show that the solubility of the REE chlorides are orders of magnitude higher than salts such as NaCl or CuCl at low water fugacities, despite their sublimation energies being substantially higher. This enhanced solubility is likely due to the high enthalpy associated with binding a single water molecule to form the species REECl 3 (H 2 O), with derived enthalpies ranging from –378 to –465 kJ/mol. Addition of further water molecules to form higher order clusters (n > 1) involves enthalpy changes of ~ -20 kJ/mol, and are in effect thermodynamically suppressed over the temperature range 350 – 425°C. Despite the enhanced solubility of small REECl 3 water clusters, simulations of boiling processes demonstrate that the REE show highly conservative behavior, partitioning strongly into the dense aqueous phase. Not surprisingly, the presence of phosphates in the system makes this effect even more pronounced, completely immobilizing the REE. This would reduce transport in both the vapor and aqueous phase to negligible levels. However, we suggest that vapor phase transport of the REE may play a significant role in systems having a relatively low partial pressure of water (below the saturation point), where the relatively high stability of the first hydrated REE chloride clusters (REECl 3 (H 2 O) and REECl 3 (H 2 O) 2 ) will give a preference for gas transport of the REE relative to other elements. This can likely happen in systems involving a gas/melt exchange in fumarolic exhalations, where water vapor discharges at relatively low (close to atmospheric) pressures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding Pore Filling Processes and Adsorption/Desorption Hysteresis in Nanoporous Metal–Organic Frameworks: Insights from Grand Canonical Monte Carlo Simulations and Free Energy Calculations

Grand canonical Monte Carlo (GCMC) simulations were used to investigate pore filling and hysteresis in nanoporous metal-organic frameworks (MOFs). Adsorption and desorption isotherms were calculated for argon at 87 K in 1866 MOFs from the CoRE MOF database and for short n-alkanes in selected MOFs, keeping the adsorbent structure rigid. Analysis of the molecular configurations showed two different mechanisms and origins of hysteresis: one involving a transition of the adsorbate arrangement in the pores similar to a gas-to-liquid transition associated with a large change in the loading and one more similar to a liquid-to-solid transition associated with a relatively small change in the loading. Our GCMC simulations in MOFs with diverse pore topologies indicate exceptions to an empirical relationship for the minimum diameter of a cylindical pore required for hysteresis as a function of the adsorbate diameter and reduced temperature. The simulations reveal some structures where isotherms exhibit two steps in the adsorption branch and only one step in the desorption branch. Hysteresis loops with a different number of adsorption and desorption steps are not common. Here, to better understand why hysteresis is observed in the GCMC simulations, the concept of the transition probability for observing a step in the adsorption isotherm at a given pressure in a GCMC simulation is introduced. We used two different methods to calculate the transition probabilities and find that these yield comparable results. Furthermore, the transition probability provides a measure for the length of GCMC simulations to yield reliable results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗