Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic Modeling”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Coupled neutronics, thermochemistry, corrosion modeling and sensitivity analyses for isotopic evolution in molten salt reactors

This study presents a computational methodology for analyzing isotopic evolution and associated uncertainties in molten salt reactors (MSRs), focusing on both fluoride- and chloride-based fuel salts. The primary goal is to enhance the understanding of isotopic behavior in MSRs and provide data to support future experimental efforts. The methodology integrates transport-coupled depletion calculations using OpenMC, equilibrium thermodynamics modeling with Thermochimica, and a corrosion model. Sensitivity analyses are performed to evaluate the impact of power density, air ingress, and humidity content on isotopic evolution in MSR concepts. This study examines representative F- and Cl-based MSR designs, highlighting the dominant influence of power density on isotopic composition, which significantly affects isotope production and depletion rates, accounting for approximately 76% of the observed variance in element concentration. Air ingress and humidity content also affect the redox potential, solubility of heavier elements, and corrosion rates, thereby altering the expected isotopic evolution in the reactor. On average, air ingress accounts for around 17% of the variance in element concentrations, while humidity explains the remaining 7%. These variances differ significantly from element to element, depending on the element’s role in depletion, redox potential evolution, and galvanic corrosion. The findings indicate that power density, air ingress, and humidity content are all critical factors for optimizing reactor design and operational strategies. Furthermore, the study provides expected ranges for key impurities in the fuel salt, which are crucial for guiding future experimental studies and refining MSR designs. Finally, this study demonstrates the importance of modeling depletion coupled with the evolution of redox potential and chemical interactions in MSR fuel salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

An ab initio molecular dynamics exploration of associates in Ba-Bi liquid with strong ordering trends

Fictive associates are widely used to describe and model liquid phases with strong ordering trends. However, little evidence is known about the assumed associates in most cases. In the present work, an ab initio molecular dynamics (AIMD) study is employed to investigate the characters of the Ba-Bi liquid, in which associates have been assumed in existing thermodynamic modeling. It is found that in the Ba rich melt, the Bi atoms are almost completely surrounded by Ba atoms. The Bi-centered coordination polyhedrons are strongly associated to crystalline structures of Ba 5 Bi 3 and Ba 4 Bi 3 with a longer lifetime than other polyhedrons during the AIMD simulations. In addition, these Bi-centered polyhedrons in Ba rich melt connect with each other through vertex, edge, face, and/or bipyramid sharing to form medium range orders (MRO). In the Bi rich melt, the Ba-centered polyhedrons also form MROs, but they are both structurally and compositionally diverse with a shorter lifetime. These findings from AIMD study provide evidences that there exist a strongly ordering Ba 4 Bi 3 associate and a weakly ordering BaBi 3 associate in the Ba-Bi liquid. Here, the predicted enthalpy of mixing in the liquid agrees well with the results by the CALPHAD modeling in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Studies on the Encapsulation of Brine Generated from a Process Using Selective Electrodialysis Reversal

This topical report summarizes the laboratory testing, geochemical modeling, and mineralogical testing conducted to examine the efficacy of wastewater encapsulation for the brine generated from the Saltworks Flex electrodialysis reversal technology. Within the laboratory testing, the short- and long-term physical properties of seven separate mix designs were examined along with the leaching of each mix design. Additionally, geochemical thermodynamic modeling was used along with mineralogical analyses to help gain a fundamental understanding of the chemistry involved in the encapsulation process. Encapsulation is a process that results in the immobilization of constituents of interest by creating an engineered material that contains the waste within its matrix, both physically and chemically, in such a manner that the receiving environment is protected from adverse effects. Within the context of this study, encapsulation mix designs involved mixing the Saltworks brine with fly ash and various additives, principally lime (calcium oxide) and cement, to form a paste that eventually sets and hardens into a monolith upon placement in a landfill. Moisture conditioning fly ash with wastewater sources for on-site landfill disposal is a process that is in use at some sites and has been presented as a viable option to dispose of liquid wastes. However, because chemical reactions that can lead to sequestration of constituents of interest are not optimized, leaching can occur. Therefore, in this study, mixes were created to mimic moisture conditioning, as well as engineered pastes of varying mix designs. The study found that engineered encapsulation mixes provide significant resistance to leaching compared to moisture conditioning alone. The best performing encapsulation mix used lime as the additive, which proved superior to the use of cement. After 90 days of curing, this encapsulated mix was 5.5 times stronger than the moisture conditioned mix. Additionally, this mix exhibited an approximate 2.5 order of magnitude drop in hydraulic conductivity compared to the moisture conditioned mix. When comparing leaching, nearly all the chloride was found to evolve from the moisture conditioned mix, whereas only 6.7% of the chlorine contained in the best performing encapsulated mix leached out. Mineralogical analysis via X-ray diffraction (XRD) showed a significant concentration of the chloride sequestering mineral, Friedel’s salt, in the best performing encapsulation mix, whereas none was observed in the moisture conditioned mix, pointing to it being the primary chemical mechanism for chloride retention, along with the physical impermeability of the material. With respect to the modeling study, the model was successful in predicting the presence of most mineral phases observed in the XRD analysis.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS2020

Stability of Actinolite on Venus

Venus currently has a hostile surface environment with temperatures of ~460 ºC, pres-sures near 92 bars, and an atmosphere composed of super critical CO2 hosting a myriad of other potentially reactive gases (e.g., SO2, HCl, HF). However, it has been proposed that its surface may not have always been so harsh. Models suggest there may have been billions of years of clement conditions allowing an Earth-like environment with liquid water oceans. If such conditions existed, it is possible Venus formed a similar array of hydrous or aqueous minerals as seen on other planets with liquid surface water (e.g., Mars, Earth). Based on thermodynamic modeling, many of these phases would not be stable under the current atmospheric conditions on Venus, dehydrating due to the high temperatures and low concentration of H2O in the atmosphere. However, the rate of decomposition of these phases may allow them to remain present on the surface over geologic time. For example, experiments on the reaction rate of tremolite (Ca2Mg5Si8O22(OH)2) show a 50% decomposition time of 2.7 Gyr for micrometer sized grains in unreactive atmospheres (i.e., without SO2) at 740 K, and a 50% decomposition time of 70 Gyr for crystals several millimeters to centimeters in size. If hydrous minerals can remain on the surface of Venus over geologic time, it has implications for our detection of evidence of these past environments, and also for the overall water budget of the planet. If after surficial dehydration the planet was able to still store water in its crust, possible processes such as subduction or metamorphism could still have operated using stored water long after liquid surface water evaporated. Several previous studies have focused on experimental investigations of mineral stability on Venus. In particular, the works of studied the decomposition rate of tremolite under conditions relevant to Venus. As their focus was on decomposition of the mineral due to lack of water in the atmosphere, their experiments were undertaken using only CO2 or N2 gas at atmospheric pressure. Re-cent experiments have examined reactivity of other minerals with the Venusian atmosphere using more complex gas compositions at similar pressures to those seen on Venus. These studies show reaction of silicate minerals with atmospheric components on relatively short timescales (i.e., on the order of days). The reported reactions of silicate materials in both studies produced iron oxides, Ca sulfates, and Na sulfates. These ions are present in many amphiboles, and Ca was proposed by Johnson and Fegley to potentially have an important role in the decomposition mechanism for tremolite, with the Ca-O bond being the first to break during decomposition. The potential involvement of Ca in both processes raises the question of whether or not the reaction to form a secondary mineral phase will influence the rate of amphibole break-down (e.g., discussion in for tremolite). Additionally, reaction of Ca with atmospheric gases may result in a different secondary mineral assemblage than simple amphibole decomposition, which will need to be recognized when searching for evidence of past hydrated minerals on the Venusian surface. In order to understand the effect of this reaction on the overall preservation potential of amphibole on the surface of Venus, we are conducting experiments in both reactive and nonreactive atmospheres using the mineral actinolite (Ca2(Mg,Fe)5Si8O22(OH)2), an amphibole with similar crystal structure to tremolite that contains both Ca and Fe.

Santos, A. R.↗
DOE OSTI2021

Investigation of aluminum nitrate as a set retarder of magnesium potassium phosphate cement: Mechanisms involved in diluted suspension

In this study, the influence of aluminum nitrate on the hydration process of magnesium potassium phosphate (MKP) cement was investigated experimentally and through thermodynamic modelling using diluted systems with Mg/PO{sub 4} molar ratio and water-to-cement ratio of 1 and 100 respectively. Three processes contributed to retard cement hydration: the acidic nature of aluminum nitrate which reduced the initial pH of the solution by two units, the early precipitation of a taranakite-like amorphous potassium aluminophosphate, which reduced the aqueous phosphate and potassium concentrations, and a salt effect due to the presence of soluble nitrate. Precipitation of potassium-containing magnesium phosphate hydrates, like Mg{sub 2}KH(PO{sub 4}){sub 2}⋅15H{sub 2}O and K-struvite, was prohibited, whereas precipitation of newberyite and cattiite was delayed. Furthermore, the modelled phase development with ongoing cement hydration agreed rather well with the experimental findings.

36 MATERIALS SCIENCE↗
DOE OSTI2023

FWP FEAA149: “Next Generation Environmental Barrier Coatings”

Environmental barrier coatings (EBCs) are required coatings for utilization of SiC/SiC ceramic matrix composite (CMC) components in gas turbines, where the EBC represents the life-limiting factor for such components. EBC/CMC systems have shown success in aero-engine applications with increased turbine inlet temperatures and improved efficiencies, which are achieved through higher temperature stability, lower density, and decreased reliance on cooling air compared to traditional superalloys. While industrial gas turbines (IGTs) do not currently utilize SiC/SiC CMCs, the current shift towards low-carbon or carbon-free fuel sources for power generation could result in a need for EBC/CMC components with higher temperature capabilities. In this work, three tasks were outlined to improve understanding of EBC lifetimes to encourage use in IGTs with carbon-free fuel such as hydrogen: 1. Define the bond coating oxidation kinetics and EBC failure criteria, 2. Measure thermal expansion coefficients of each layered material, and 3. Perform advanced characterization and modeling to assess EBC lifetimes. Cyclic steam oxidation tests were conducted on various EBC/Si/SiC chemistries and EBC/SiC architectures to define substrate oxidation kinetics and EBC failure modes. An open-source code was developed to quantify the undulating thermally grown oxide thickness with thousands of measurements from specimen cross-section images. Bond coating oxidation kinetics were determined and used to develop a kinetic and thermodynamic model for predicting EBC lifetimes. High-temperature Raman spectroscopy was utilized for determining the SiO 2 thermally grown oxide phase transformation as the life-limiting feature for EBCs. Model efforts supported the claim that the SiO 2 phase transformation causes elevated stress during thermal cycling with associated cracking that decreases the adhesion strength of the EBC, eventually resulting in coating spallation. The finite element model subroutine will be made publicly available upon internal review. Further development of an EBC lifetime model for IGTs involves definition of a critical SiO 2 thickness for EBC spallation and must also consider both environmental (gas velocity, pressure, etc.) and specimen (EBC dopants, layer architectures, etc.) effects into predicted bond coating oxidation kinetics for long-lifetime components.

36 MATERIALS SCIENCE↗
NASA NTRS1983

Numerical simulation of the Weddell Sea pack ice

In order to determine the degree to which a coupled dynamic thermodynamic model can reproduce the seasonal cycle of the Antarctic Sea ice in the Weddell Sea, a series of two-year simulations of the Weddell Sea ice pack were carried out. The simulations employed Hibler's (1979) model and used a one-day time step on an 18 x 15 grid with a resolution of 222 km. Daily atmospheric data from 1979 were used to drive the simulations. The simulations yielded a seasonal cycle of ice with maximum and minimum extents close to that observed. Except for portions of the western Weddell, the advance of the ice is found to be primarily thermodynamic in nature, while the rapid decay depends critically on the presence of both leads and lateral ice advection. In early summer the ice motion causes a residual tongue of ice to extend eastward from the Antarctic Peninsula in agreement with observations. Mean ice thicknesses are consistent with observations, as are the mean drift rates of about 5 km/day.

Hibler, W. D., III↗
NASA NTRS2013

Development of a Phasor Diagram Creator to Visualize the Piston and Displacer Forces in an Advanced Stirling Convertor

The steady state, nearly sinusoidal behavior of the components in a Free Piston Stirling Engine allows for visualization of the forces in the system using phasor diagrams. Based on Newton's second law, F=ma, any phasor diagrams modeling a given component in a system should close if all of the acting forces have been considered. Since the Advanced Stirling Radioisotope Generator (ASRG), currently being developed for future NASA deep space missions, is made up of such nearly sinusoidally oscillating components, its phasor diagrams would also be expected to close. A graphical user interface (GUI) has been written in MATLAB by taking user input data, passing it to Sage, a 1-D thermodynamic modeling program used to model the Stirling convertor, running Sage and then automatically plotting the phasor diagrams. Using this software tool, the effect of varying different Sage inputs on the phasor diagrams was determined. The parameters varied were piston amplitude, hot end temperature, cold end temperature, operating frequency, and displacer spring constant. By using these phasor diagrams, better insight can be gained as to why the convertor operates the way that it does.

Saha, Dipanjan↗
NASA NTRS2013

Development of a Phasor Diagram Creator to Visualize the Piston and Displacer Forces in an Advanced Stirling Convertor

The steady-state, nearly sinusoidal behavior of the components in a free-piston Stirling engine allows for visualization of the forces in the system using phasor diagrams. Based on Newton's second law, F = ma, any phasor diagrams modeling a given component in a system should close if all of the acting forces have been considered. Since the Advanced Stirling Radioisotope Generator (ASRG), currently being developed for future NASA deep space missions, is made up of such nearly sinusoidally oscillating components, its phasor diagrams would also be expected to close. A graphical user interface (GUI) has been written in MATLAB (MathWorks), which takes user input data, passes it to Sage (Gedeon Associates), a one-dimensional thermodynamic modeling program used to model the Stirling convertor, runs Sage, and then automatically plots the phasor diagrams. Using this software tool, the effect of varying different Sage inputs on the phasor diagrams was determined. The parameters varied were piston amplitude, hot-end temperature, cold-end temperature, operating frequency, and displacer spring constant. These phasor diagrams offer useful insight into convertor operation and performance.

Saha, Dipanjan↗
DOE OSTI2020

MODELING OF ENHANCED AIR DEHUMIDIFICATION THROUGH ELECTRICALLY CHARGED VAPOR CAPTURING ELECTROSTATIC DROPLETS

This paper investigated the potential of a new electrostatic-based water vapor separation system for humidity control in buildings. The Electrospray Vortical Flow eXchanger (EVFX) consists of small electrically charged water droplets released inside air vortexes to attract water vapor molecules, effectively wiping out the humidity from the air. Larger droplets are then removed at the outer wall of the vortex. If a water molecule is placed in a gradient electric field, such a polar molecule experiences a dielectrophoresis force, which moves the molecules and produces a vapor concentration gradient. When the vapor density exceeds the saturation level, nucleation occurs. A thermodynamic model that captures this phenomenon was described in this paper. Electrosprayed droplets were initially charged close to their maximum electrical charge dictated by the Rayleigh limit; these droplets were used as nucleation centers. The high electric charge decreased the pressure of the vapor, which was in equilibrium with the droplets. This enhanced the charged droplets' nucleation and growth that depleted the vapor phase near the droplets, which was compensated for by the dielectrophoresis flow and diffusion. Significant water vapor condensation was predicted for small droplets at a high electrical charge. Beyond a specific droplet size, while the condensation of water vapor was still theoretically possible, it was somewhat limited because the further growth of the large charged droplets was inhibited due to their high inherent water vapor pressure near the surface. A case study illustrated the potentials and limitation of the EVFX for droplets charged to the Rayleigh limit, for which the effective range of diameters resulted between 10 nm up to few micron.

Morcelli, Stefano↗
DOE OSTI2026

A wide-range thermodynamically consistent constitutive model for dynamic loading of partially saturated porous media

SUMMARY This paper presents a novel thermodynamically consistent constitutive model for partially saturated porous rocks across a wide range of conditions. The material states generated behind the shock wave from an explosive source can vary significantly, ranging from crushed and melted rock near the source to a poroelastic medium in the far field. In the model, rock strength is determined by the effective pressure, which is calculated using two independent equations of state: one for the solid rock and another for the pore fluid. The model accounts for shock-induced liquefaction resulting from fluid pressure buildup in the pore spaces near the explosive source. Simultaneously, it describes the increase in wave propagation speed due to elastic pore contraction in both dry and partially saturated rocks. This model is applied to investigate how fluid saturation affects the amplitude and shape of the generated waves, as well as the residual stress surrounding the cavity formed by spherical explosions.

Vorobiev, Oleg [Lawrence Livermore National Labora↗
DOE OSTI2026

Optimal Operation of Residential High Performance Water Heater for Reduction of Electricity Cost and Peak Demand Through Field Validation

Water heating accounts for about 18% of a typical US home’s energy use. Modern water heaters have enabled control options through APIs, offering customers the opportunity to reduce their energy cost and peak demand by dynamically adjusting settings. A water heater’s capacity to store energy using its storage tank makes it an asset for peak demand reduction and energy cost savings. For this reason, a mixed-integer linear programming model is proposed to minimize the energy cost of a high-performance water heater while also reducing the peak demand of the residential household under a time-of-use utility rate by dynamically changing the water heater’s running mode. Specifically, a multi-objective optimization model is formulated to determine the mode settings of the water heater considering hot water use, time-of-use rate, and peak demand limit of the residential household. The mode settings are associated with different dead bands of water temperature for triggering on/off action of the heat pump and heating element. A 66-gal hybrid electric high performance water heater was used for numerical simulation and practical experiments. The simulation results were well aligned with measurements of practical experiments, validating the soundness of the thermodynamic model. In addition, reductions of energy cost, enabling affordability, and reducing peak demand are demonstrated. The research team also developed a software framework with dashboards to automatically and continuously monitor and manage devices.

Liu, Guodong [ORNL] (ORCID:0000000213498608)↗
NASA NTRS2014

Liquid Oxygen Thermodynamic Vent System Testing with Helium Pressurization

This report presents the results of several thermodynamic vent system (TVS) tests with liquid oxygen plus a test with liquid nitrogen. In all tests, the liquid was heated above its normal boiling point to 111 K for oxygen and 100 K for nitrogen. The elevated temperature was representative of tank conditions for a candidate lunar lander ascent stage. An initial test series was conducted with saturated oxygen liquid and vapor at 0.6 MPa. The initial series was followed by tests where the test tank was pressurized with gaseous helium to 1.4 to 1.6 MPa. For these tests, the helium mole fraction in the ullage was quite high, about 0.57 to 0.62. TVS behavior is different when helium is present than when helium is absent. The tank pressure becomes the sum of the vapor pressure and the partial pressure of helium. Therefore, tank pressure depends not only on temperature, as is the case for a pure liquid-vapor system, but also on helium density (i.e., the mass of helium divided by the ullage volume). Thus, properly controlling TVS operation is more challenging with helium pressurization than without helium pressurization. When helium was present, the liquid temperature would rise with each successive TVS cycle if tank pressure was kept within a constant control band. Alternatively, if the liquid temperature was maintained within a constant TVS control band, the tank pressure would drop with each TVS cycle. The final test series, which was conducted with liquid nitrogen pressurized with helium, demonstrated simultaneous pressure and temperature control during TVS operation. The simultaneous control was achieved by systematic injection of additional helium during each TVS cycle. Adding helium maintained the helium partial pressure as the liquid volume decreased because of TVS operation. The TVS demonstrations with liquid oxygen pressurized with helium were conducted with three different fluid-mixer configurations-a submerged axial jet mixer, a pair of spray hoops in the tank ullage, and combined use of the axial jet and spray hoops. A submerged liquid pump and compact heat exchanger located inside the test tank were used with all the mixer configurations. The initial series without helium and the final series with liquid nitrogen both used the axial jet mixer. The axial jet configuration successfully demonstrated the ability to control tank pressure; but in the normal-gravity environment, the temperature in the upper tank region (ullage and unwetted wall) was not controlled. The spray hoops and axial jet combination also successfully demonstrated pressure control as well as temperature control of the entire tank and contents. The spray-hoops-only configuration was not expected to be a reliable means of tank mixing because there was no direct means to produce liquid circulation. However, surprisingly good results also were obtained with the sprayhoops- only configuration (i.e., performance metrics such as cycle-averaged vent flowrate were similar to those obtained with the other configurations). A simple thermodynamic model was developed that correctly predicted the TVS behavior (temperature rise or pressure drop per TVS cycle) when helium was present in the ullage. The model predictions were correlated over a range of input parameters. The correlations show that temperature rise or pressure drop per cycle was proportional to both helium mole fraction and tank heat input. The response also depended on the tank fill fraction: the temperature rise or pressure drop (per TVS cycle) increased as the ullage volume decreased.

Thermodynamic vent system↗
NASA NTRS1993

An experimental study of trace element partitioning between perovskite, hibonite and melt: Equilibrium values

The presence of perovskite (CATiO3) and hibonite (Ca Al12O19) within different regions of Calcium-, Aluminum-rich Inclusions (CAI) and the trace element concentrations of these minerals in each circumstance, constrain models of precursor formation, nebular condensation, the thermal history of inclusions with relict perovskite and hibonite, and the formation of the Wark-Lovering rim. At present mineral/melt partition coefficient data for hibonite are limited to a few elements in simple experimental systems, or to those derived from hibonite-glass pairs in hibonite/glass microspherules. Similarly, there is only limited data on perovskite D that are applicable to meteorite compositions. Apart from the importance of partitioning studies to meteorite research, D values also are invaluable in the development of thermodynamic models, especially when data is available for a large number of elements that have different ionic charge and radii. In addition, study of the effect of rapid cooling on partitioning is crucial to our understanding of meteorite inclusions. To expand our knowledge of mineral/melt D for perovskite and hibonite, a study was instituted where D values are obtained in both equilibrium and dynamic cooling experiments. As an initial phase of this study mineral/melt D was measured for major elements (Ca, Mg, Al, Ti, and Si), 15 rare earth elements (La-Lu) and 8 other elements (Ba, Sr, U, Th, Nb, Zr, Hf, and Ge) in perovskite and hibonite grown under equilibrium conditions, in bulk compositions that are respectively similar to Compact Type A (CTA) CAI and to a hibonite/glass microspherule. Experimental mixes were doped with REE at 20-50x chondritic (ch) abundances, Ba at 50 ppm, Sr, Hf, Nb, and Zr at 100 ppm and, U and Th at 200 ppm. Trace element abundances were measured with the PANURGE ion microprobe. Major element compositions were obtained by electron microprobe analysis.

Kennedy, A. K.↗
NASA NTRS2016

Metabolic Potential and Activity in Fluids of the Coast Range Ophiolite Microbial Observatory, California, USA

Metabolic potential and activity associated with hydrogen and carbon monoxide were characterized in fluids sampled from the the Coast Range Ophiolite Microbial Observatory (CROMO). CROMO consists of two clusters of science-dedicated wells drilled to varying depths up to 35m in the actively serpentinizing, Jurassic-age Coast Range Ophiolite of Northern California, along with a suite of pre-existing monitoring wells at the same site. Consistent with the fluid chemistry observed in other serpentinizing systems, CROMO fluids are highly alkaline, with pH up to 12.5, high in methane, with concentrations up 1600 micromolar, and low in dissolved inorganic carbon (DIC), with concentrations of 10's to 100's of micromolar. CROMO is conspicuous for fluid H2 concentrations that are consistently sub-micromolar, orders of magnitude lower than is typical of other systems. However, higher H2 concentrations (10's -100's of micromolar) at an earlier stage of fluid chemical evolution are predicted by, or consistent with: thermodynamic models for fluid chemistry based on parent rock composition equivalent to local peridotite and with water:rock ratio constrained by observed pH; the presence of magnetite at several wt% in CROMO drill cores; and concentrations of formate and carbon monoxide that would require elevated H2 if formed in equilibrium with H2 and DIC. Calculated Gibbs energy changes for reaction of H2 and CO in each of several metabolisms, across the range of fluid composition encompassed by the CROMO wells, range from bioenergetically feasible (capable of driving ATP synthesis) to thermodynamically unfavorable. Active consumption relative to killed controls was observed for both CO and H2 during incubation of fluids from the pre-existing monitoring wells; in incubations of freshly cored solids, consumption was only observed in one sample set (corresponding to the lowest pH) out of three. The specific metabolisms by which H2 and CO are consumed remain to be determined.

Hydrogen↗
DOE OSTI2025

Characterization of kerogen nanopores using 2D NMR relaxation and MD simulations

The characterization of kerogen nanopores is crucial for predicting the geostorage capacity and recoverability of natural gas in unconventional gas shale reservoirs. Towards this end, a powerful technique is presented which integrates 2D NMR T 1 -T 2 relaxation measurements with molecular dynamics (MD) simulations of hydrocarbons confined in the nanopores of kerogen. The integrated NMR-MD technique is demonstrated using T 1 -T 2 measurements of kerogen isolates and organic-rich chalks saturated with heptane, together with MD simulations of heptane completely dissolved in a realistic kerogen model. The NMR-MD results are used to extract the swelling ratio and nanopore size distribution of kerogen as a function of depth in the reservoir. The effects of organic nanoconfinement on the T 1 relaxation dispersion and T 2 residual dipolar coupling of heptane are investigated, as well as the effect of downhole effective stress on the kerogen nanopore size as a function of depth and compaction. Potential applications in partially depleted gas shale reservoirs are discussed, including CO 2 utilization/geostorage, geostorage of green H 2 , and integration of the NMR-MD technique with thermodynamic models for predicting the competitive sorption of gas mixtures in kerogen.

Compaction↗
NASA NTRS1985

A self-consistent linear-mode model of stellar convection

A normal-mode expansion of the linearized fluid equations in terms of small subset of spherical harmonics can provide a foundation for a physically motivated, self-consistent description of a solar-type convection zone. In the absence of dissipation, a second-order differential equation governs the radial dependence of the modes, so that interpretation of the effects on convection quantities of the normal-form 'potential well' is straightforward. The philosophy is quite different from the more recent work of Narasimha and Antia (1982): all envelopes presented here differ substantially from MLT envelopes, and therefore, from theirs, which are constructed to be consistent with MLT. The amplitude of all modes is set by a Kelvin-Helmholtz-('shear'-) instability argument unrelated to solar observations, with the result that the convection description may be considered to arise from 'first-hueristic-principles'. The thermodynamics modelled vaguely resemble the sun's, and more vigorously convective envelopes show some phenomena qualitatively like solar observations (e.g., atmospheric velocity spectra).

Macauslan, J.↗