SEARCH · Engineering Papers

Results for “subsurface diffusion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Discrete Fracture Network Modeling to Estimate Upscaled Parameters for the Topopah Spring, Lava Flow, and Tiva Canyon Aquifers at Pahute Mesa, Nevada National Security Site

This report describes the results of Discrete Fracture Network (DFN) simulations for the Topopah Spring Aquifer (TSA), Lava Flow Aquifer, and Tiva Canyon Aquifer (TCA), at Pahute Mesa on the Nevada National Security Site (NNSS), formerly the Nevada Test Site. The research focuses on calculating upscaled groundwater flow and contaminant transport parameters using DFNs generated according to fracture characteristics observed in the TSA, LFA and TCA at Pahute Mesa. The highly fractured and heterogeneous nature of these aquifers makes them candidates for stochastic DFN modeling of radionuclide transport on a small scale with subsequent upscaling. One hundred independent DFN realizations are generated for each aquifer, and the upscaled parameters for continuum simulations of subsurface flow and transport in fractured media at Pahute Mesa are calculated. Our goal is to implement a modeling approach that can translate parameters to larger-scale models that account for local-scale flow and transport processes, such as channelization of flow and transport along a few well connected, large fractures. Additionally, to simulate advective and advective-diffusive transport through the fracture networks, the Time Domain Random Walk (TDRW) approach is applied to account for matrix diffusion into a finite half-space. Moreover, a novel approach to calculate dynamic (active) fracture surface area to reflect flow channeling is implemented. This work will improve the representation of radionuclide transport processes in largescale, regulatory-focused models by providing estimates of hard-to-measure flow and contaminant transport parameters at large scales. In this report, we (1) show recent results of flow and transport simulations on multiple DFN realizations of the TSA, LFA, TCA; (2) discuss the resulting distributions of estimated upscaled parameters; (3) describe the estimation of upscaled parameters for an equivalent parallel-plate continuum model and (4) present a comparison between simulated transport from the equivalent continuum model and an actual DFN.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Stochastic Ensemble Generation for Improved Characterization of Representing Geologic Variability in a Reservoir: IBDP Case Study for SMART Initiative

This document is a poster covering the findings from activities on training data generation, specifically geologic ensemble generation. The generated geologic realizations captured the range of possible permeability distributions of the subsurface at the Illinois Basin - Decatur Project (IBDP) site, based on available well log variabilities. The percentages of reservoirs and baffles in the injection zone and a truncation of baffle permeability led to more variance in the simulations. This will be used to build forward modeling, history matching, and optimization workflows. The geologic realizations were also ranked according to dynamic measures of hydraulic diffusivity, and simulations confirm a greater contrast between the reservoir and the baffles during injection.

stochastic ensemble generation↗

DOE OSTI2020

Helium effects on the surface and subsurface evolutions in single-crystalline tungsten

Tungsten (W) has been perceived as one of the most promising plasma facing materials (PFMs) for future fusion reactors. In the past decade, its behavior under irradiation and helium (He) plasma interaction has been extensively studied. However, some key knowledge gaps still exist, such as the influence of crystallographic orientation on the surface and subsurface evolutions. In this work, we focus on the He ion-beam irradiation damage effects in mirror-polished single-crystalline W samples with three different surface planes of {100}, {110} and {111}. Irradiation was performed at room temperature using 40 keV He+ to a fluence of 1 × 10 16 /cm 2 , followed by thermal desorption spectroscopy (TDS) up to ~1920 K. The microstructures of He-irradiated W before and after TDS heat treatment were characterized by scanning and transmission electron microscopy. Subsurface He bubbles were imaged in all irradiated samples, but newly formed <111>-oriented surface grains and surface blisters were only observed in W {100} and {110} starting orientations. These results reveal that radiation damage, He thermal desorption, and surface/subsurface evolution are all strongly dependent on crystallographic orientation. Underlying physical mechanisms are discussed based on ion channeling effects, He-vacancy interactions, and surface diffusion. In conclusion, these findings provide new insights into He effects in W.

36 MATERIALS SCIENCE↗

NASA NTRS1985

The rise and fall of sunspot group 18962 - A case of magnetic submergence

Normal sunspot groups emerge from below the surface as a big arch; the footpoints form a dipole that moves apart and the fields eventually diffuse and disappear. BBSO 18962 was a small dipole with elements which, after spreading apart in the first day of existence, drew back together and disappeared below the surface. It is proposed that the flux loop was pulled back down by magnetic tension and submerged. The circumstances suggest that this behavior may not be uncommon. It implies that some dipoles do not decouple from subsurface fields.

Zirin, H.↗

NASA NTRS1973

Radar signal return from near-shore surface and shallow subsurface features, Darien Province, Panama

The AN/APQ-97 radar imagery over eastern Panama is analyzed. The imagery was directed toward extraction of geologic and engineering data and the establishment of operational parameters. Subsequent investigations emphasized landform identification and vegetation distribution. The parameters affecting the observed return signal strength from such features are considered. Near-shore ocean phenomena were analyzed. Tidal zone features such as mud flats and reefs were identified in the near range, but were not detectable in the far range. Surface roughness dictated the nature of reflected energy (specular or diffuse). In surf zones, changes in wave train orientation relative to look direction, the slope of the surface, and the physical character of the wave must be considered. It is concluded that the establishment of the areal extent of the tidal flats, distributary channels, and reefs is practical only in the near to intermediate range under minimal low tide conditions.

Hanson, B. C.↗

DOE OSTI2021

The subsurface structure of abraded Al–Zn–Mg–Cu alloy

Surface abrasion has been shown to introduce heavy shear stress and produce a few hundred nanometers thick altered surface layer (ASL) on Al alloys. Such ASLs exhibit microstructural changes at room temperature and different corrosion resistance than the underlying substrate. Here, we report how sample dimensions and temperature affect the deformed subsurface microstructure evolution on an abraded Al–Zn–Mg–Cu alloy using transmission electron microscopy (TEM) and atom probe tomography. The ASL microstructure evolves differently in the bulk sample than in the thin TEM specimen. In the freshly abraded ASL, the original η'/η precipitates and grain boundaries were deformed toward the abrasion direction. After 7-day natural aging of the freshly abraded bulk sample, well-defined subgrains formed locally in the ASL, and unusual Al 2 Cu (θ) and Mg x Zn y phases precipitated at the subgrain boundaries. However, during natural aging of this TEM specimen taken from the abraded and aged bulk samples, much faster changes occurred, including the continued formation of subgrains and θ particles in the entire ASL, and the diffusion of Zn atoms out of Mg x Zn y particles to form a Zn particle. On the other hand, no obvious microstructural changes occurred in the ASL during 3-month natural aging of the TEM specimen immediately taken from the freshly abraded sample. Finally, this work advances the understanding of how Al alloy microstructures respond during and following any shear deformation, and complex and unexpected effects of sample dimensions on the behavior.

36 MATERIALS SCIENCE↗

NASA NTRS2003

Polar and Non-Polar Layers on Mars: A Single Mechanism for Formation?

The recent discovery of vast quantities of near-subsurface ice in both polar regions of Mars by the Mars Odyssey Gamma Ray Spectrometer (GRS) has presented us with an interesting quandary. On one hand, these deposits, found poleward of 60 deg in both hemispheres, are consistent with thermal models suggesting ice will be best protected in these regions during periods of high obliquity. On the other hand, the current paradigm regarding the placement of these deposits, i.e., diffusive deposition of water vapor, appears to be inconsistent with the large volume mixing ratios (approx. 90%) inferred from the GRS data. This incongruity argues that diffusion alone cannot be the primary mechanism for the creation of these reservoirs, and that an alternate, large-scale process should be considered.

Mischna, M. A.↗

DOE OSTI2025

Behaviour of cadmium isotopes in sulfidic waters and sediments of the Black Sea: Implications for global cadmium cycling and the application of cadmium isotopes as a paleo-oceanographic proxy

The cadmium isotope system has found use as a tracer for biological productivity, redox and organic carbon burial in the oceans. There are, however, very few observational constraints on Cd isotopic behaviour in modern sulfidic marine conditions, limiting our understanding of the modern Cd cycle, and our ability to use Cd isotope measurements of ancient sedimentary deposits as a paleoceanographic proxy. Here we study the behaviour of dissolved Cd in the water column of the Black Sea and its incorporation into sediments on the basin floor. The isotopic composition of dissolved Cd in the upper ∼50 m of the Black Sea water column is controlled by a combination of biological uptake and regeneration along with mixing of river water and Mediterranean seawater. Cadmium declines to <2 % of its peak subsurface concentration within the zone of nitrate-reduction but above the sulfide chemocline. The isotopic composition of dissolved Cd below the sulfide chemocline evolves in a manner that is consistent with experimental Cd sulfide fractionation factors, providing field-based support for these earlier studies. In contrast, the stability of the dissolved Cd isotopic signature above the chemocline, despite a large reduction in dissolved Cd concentrations, is explained by diffusion towards the Cd-sulfide sink below the chemocline. Sediments accumulating in the deep Black Sea below the sulfide chemocline have isotopic compositions that are similar to dissolved Cd in its near-surface aqueous concentration maximum. Sediments accumulating at shallower depths have lighter isotopic compositions that are likely affected by non-sulfidic burial pathways mediated by macrofaunal mixing and diagenesis of shallow burial phases. Our study shows that the isotopic composition of bulk sediments accumulating under sulfidic conditions is similar to dissolved Cd in the upper water column, supporting its use as a paleo-chemical proxy for ancient seawater as long as basin-scale controls on open ocean seawater sources are considered. We show that Cd burial in the deep Black Sea, predominantly as CdS, accounts for ∼0.5–15 % of total annual Cd burial in the global ocean.

Biogeochemical cycles↗

DOE OSTI2020

Surface effects on deuterium permeation through vanadium membranes

Dense vanadium-based membranes offer high permeability and perfect selectivity to hydrogen isotopes, maintain favorable neutronic properties, and are compatible with liquid metals such as PbLi. These properties make vanadium membranes a promising fusion fuel cycle technology for processes such as tritium extraction from PbLi and exhaust processing. Surface contamination has a deleterious effect on the gas-phase hydrogen permeation through vanadium, and the reported permeabilities range from 10 -14 to 10 -7 mol m -1 s -1 Pa -0.5 . Thin dense films of palladium applied to clean vanadium surfaces enable a consistently high hydrogen permeability. In this study, uncoated vanadium resulted in deuterium permeabilities ranging from 2.8 × 10 -11 to 6.4 × 10 -9 mol m -1 s -1 Pa -0.5 at 300 °C–700 °C, respectively. Post-test analysis revealed a VO x surface layer and VC x subsurface layer formed on the feed side, while the as-received surface oxide dissolved leaving a submonolayer oxide on the permeate surface. Furthermore, the Pd-coated V resulted in a maximum deuterium permeability of 2.1 × 10 -7 m -1 s -1 Pa -0.5 at 375 °C upon activation of the Pd surface by oxidation and reduction. The deuterium permeation declined upon heating to 500 °C due to intermetallic diffusion between the Pd and V. The Mo 2 C-coated V resulted in deuterium permeabilities ranging from 2.7 × 10 -10 to 1.8 × 10-9 at 500 °C–700 °C, respectively, and a post-test analysis found the carbon in the Mo 2 C layer had dissolved into the V near the interface.

13 HYDRO ENERGY↗

DOE OSTI2020

Physics-Informed Deep Neural Networks for Learning Parameters and Constitutive Relationships in Subsurface Flow Problems

In this paper, we present a physics informed deep neural network (DNN) method for estimating parameters and unknown physics (constitutive relationships) in partial differential equation (PDE) models. We use PDEs in addition to measurements to train DNNs to approximate unknown parameters and constitutive relationships as well as states. The proposed approach increases the accuracy of DNN approximations of partially known functions when a limited number of measurements is available and allows for training DNNs when no direct measurements of the functions of interest are available. We employ physics informed DNNs to estimate the unknown space-dependent diffusion coefficient in a linear diffusion equation and an unknown constitutive relationship in a non-linear diffusion equation. For the parameter estimation problem, we assume that partial measurements of the coefficient and states are available and demonstrate that under these conditions, the proposed method is more accurate than state-of-the-art methods. For the non-linear diffusion PDE model with a fully unknown constitutive relationship (i.e., no measurements of constitutive relationship are available), the physics informed DNN method can accurately estimate the non-linear constitutive relationship based on state measurements only. Finally, we demonstrate that the proposed method remains accurate in the presence of measurement noise.

42 ENGINEERING↗

DOE OSTI2020

High-temperature and high-pressure NMR investigations of low viscous fluids confined in mesoporous systems

Abstract In this contribution, the relaxation and diffusional behaviors of low viscous fluids, water and methanol confined into mesoporous silica and controlled size pore glass were investigated. The engineered porous systems are relevant to geologically important subsurface energy materials. The engineered porous proxies were characterized by Brunauer–Emmett–Teller (BET) surface analyzer, nuclear magnetic resonance (NMR) spectroscopy, and electron microscopy (EM) to determine surface area, pore-wall protonation and morphology of these materials, respectively. The confined behavior of the low viscous fluids was studied by varying pore diameter, fluid-to-solid ratio, temperature, and pressure, and then compared to bulk liquid state. Both relaxation and diffusion behaviors for the confined fluids showed increasing deviation from pure bulk fluids as the fluid-to-solid ratio was decreased, and surface-to-volume ratio ( S / V ) was varied. Variable pressure deuteron NMR relaxation of confined D 2 O and confined methanol, deuterated at the hydroxyl or methyl positions, were performed to exploit the sensitivity of the deuteron quadrupole moment to molecular rotation. The methanol results demonstrated greater pressure dependence than those for water only in bulk. The deviations from bulk liquid behavior arise from different reasons such as confinement and the interactions between confined fluid and the nano-pore wall. The results of the present report give insight into the behavior of low viscosity fluid in nano-confined geometries under different state conditions.

Chemistry↗

DOE OSTI2021

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRSDate not supplied

Improving Ocean Reanalyses and ENSO Forecasts By Assimilation of Rain Corrected Satellite Sea Surface Salinity Using the GMAO S2S Forecast System

During the past years, we have seen that the La Nina to El Nino transition has had a significant global socio-economic impact and so has been the key focus for improving coupled ocean-atmosphere forecasts. Assimilation of satellite altimetry and subsurface temperature and salinity from (mostly) Argo help to improve the initialization of the thermocline, while satellite Sea Surface Temperature (SST) aids in constraining surface heat-fluxes, leading to improved subseasonal to seasonal forecasts of the coupled system. However, few studies have focused on improving the fresh-water flux and near-surface density and mixing through assimilation of satellite sea surface salinity (SSS). For expediency, the few ocean models that do assimilate SSS bias-correct the satellite SSS data to normalize towards the near-surface Argo data. However, in rainy regions, where buoyant water sits as a fresh lens at the surface, this assumption is likely inadequate. In previous work, we have shown that adjusting SSS data to bulk salinity (Sb) using the Rain Impact Model (RIM) of Santos-Garcia et al., 2014 has improved the near-surface density and mixed layer depth, leading to deeper thermocline and improved the NINO3.4 SST forecasts. Now we utilize the Soil Moisture/Ocean Salinity, Rain Corrected (SMOS_RC) SSS product provided by the Centre Aval de Traitement des données SMOS (CATDS CPDC) to represent the Sb more accurately at first model layer (in our case 5 m). Rather than using a diffusivity model as with RIM, SMOS_RC relies on an observed relationship between the spatial heterogeneity of SMOS SSS and instantaneous rain rate (RR) (Supply et al., 2020). In order to test the impact of SMOS_RC versus SMOS, we compare two reanalyses over the period 2014 to 2021. For both reanalysis experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF scheme (Penny et al., 2013). One reanalysis additionally assimilates SMOS SSS data as is (i.e., with the fresh bias), and a separate reanalysis is performed assimilating the SMOS_RC data. We assess the impact for near-surface and subsurface dynamics within ocean reanalyses by validating against observations and explore how SSS assimilation (SMOS versus SMOS_RC) impacts dynamical ENSO forecasts using the NASA GMAO Sub-seasonal to Seasonal coupled forecast system (GEOS S2S-3, Molod et al., 2020, Hackert et al., 2023). We will show that improved SSS estimates and near-surface density and mixing led to more accurate coupled air/sea interaction and better ENSO forecasts.

Eric Hackert↗

NASA NTRSDate not supplied

The Stop-and-Go Mechanism: Towards an Integrated Approach to Model Seismicity, Outgassing, Deformation, and Thermal Unrest at Active Volcanoes

Connecting the geophysical and geochemical signals recorded at and above the surface of volcanoes with source mechanisms is fundamental to understand transitions from quiescence to eruption, and to integrate the behavior of volcanoes with their regional seismotectonic and hydrological context. Some of these signals include shallow volcanic tremor, a long-lasting (minutes-to-years) ground vibration detected in volcanic areas during unrest; volcanic outgassing, which shows intriguing periodicities over multiple scales, as revealed by space-borne instruments and by ground-based, high-frequency (~1 Hz), sampling methods; inflation and deflation of volcanic edifices, commonly associated with subsurface volume changes; and large-scale (from a few to tens of km2), low-temperature, thermal anomalies, a newly-discovered signal that has been observed to emerge on volcanic flanks from months-to-years prior to gas- and magma-driven eruptions. Traditionally, these geophysical and geochemical signals have been modelled independently, which limits our interpretation of subsurface processes and thus our assessment of unrest. In this work, we show through lumped-parameter models, numerical simulations, and preliminary laboratory experiments that many of the signals recorded around volcanoes can be integrated into a common “stop-and-go” mechanism. In particular, we found that they can arise spontaneously when taking into account the feedbacks between the temporary accumulation of gas and/or magma in the crust (“stop”), and the diffusive transfer of gas and heat towards the surface (“go”). For example, we find that the “stop-and-go” mechanism can explain: (i) the emergence of monochromatic, broadband, and harmonic tremor; (ii) the periodic components commonly observed in volcanic outgassing time series; (iii) the link between large-scale thermal anomalies and shallow hydrothermal systems; (iv) the reported time lags between large-scale thermal anomalies and deformation; and (v) the breathing of potentially destructive volcanic calderas. The “stop-and-go” mechanism provides a consistent and realistic framework to link multiple geophysical and geochemical signals with the processes leading to volcanic unrest and eruption.

volcanoes↗

DOE OSTI2024

Mechanisms for Microseismicity Occurrence Due to CO 2 Injection at Decatur, Illinois: A Coupled Multiphase Flow and Geomechanics Perspective

Here, we numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO 2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologically consistent model that honors key stratigraphic horizons and 3D fault surfaces interpreted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO 2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO 2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection interval to the basement is essential to explain the destabilization of the regions where microseismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO 2 sequestration operations.

58 GEOSCIENCES↗

DOE OSTI2025

Numerical simulation of flow and mixing in fracture intersections

Fluid transport through fractured geological formations is strongly influenced by the redistribution of solutes at fracture intersections. In this study, we perform detailed numerical simulations of flow and scalar transport within the intersection of two smooth, planar fractures. The analysis focuses on the mixing ratio, the proportion of solute flux exiting along the outlet branch aligned with the primary inlet flow direction, relative to the total solute flux at the outlets. We systematically investigate how the mixing ratio varies with four key parameters: Peclet number, Reynolds number, flow rate ratio between outlet branches, and fracture intersection angle. Results show that the mixing ratio decreases with increasing Peclet number and outlet flow rate ratio, consistent with reduced diffusive spreading and enhanced streamline routing. While low Reynolds numbers have minimal impact, inertial effects at higher Reynolds numbers significantly increase the mixing ratio. Additionally, acute and obtuse intersection angles alter flow partitioning and modify the solute distribution at the outlets. These findings provide a quantitative basis for incorporating physically realistic mixing behavior—intermediate between complete mixing and streamline-following assumptions—into network-scale transport models. The results have direct relevance to subsurface energy systems, including geothermal energy production, carbon sequestration, and contaminant remediation.

58 GEOSCIENCES↗

DOE OSTI2023

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2012

Astrobiological Journeys to and from the South Polar Sea of Enceladus - Bidirectional Interactions with the Saturn Magnetosphere

The spectacularly cryovolcanic moon Enceladus is a major source of plasma for the Saturn magnetosphere via ionization of the ejected molecular species and ice grains. Field-aligned plasma flows from the Enceladus environment visibly impact the moon's magnetic footpoint in Saturn's polar auroral region, while water group and other ions from the moon emissions diffuse radially throughout the magnetosphere and may be the dominant source of oxygen for Titan's oxygen-poor upper atmosphere. But the moon-magnetosphere interaction is bidirectional in the sense that the moon surface is globally exposed to constant irradiation by the returning magnetospheric ions and by energetic electrons from the field-aligned and radially diffusing populations. The returning ion source operates both on global scales of the magnetosphere and locally for highly reactive species produced in the ejecta plume. All of these sources likely combine to produce a highly oxidized global surface layer. Since plasma electrons and ions are cooled by interaction with neutral gas and E-ring ice grains from Enceladus, the moon emissions have a governing effect on the seed populations of energetic particles that irradiate the surface. The proposed subsurface polar sea and transient crustal overturn in the south polar region could bring the polar surface oxidants into contact with hydrocarbons and ammonia to make oxidation product gases contributing to the cryovolcanic jets, a process first proposed by Cooper et al. (Plan. Sp. Sci., 2009). As has been previously suggested for Europa, the 'oxidants could contribute to enhanced astrobiological potential of Enceladus, perhaps even higher than for Europa where organic hydrocarbons have not yet been directly detected. Unlike Europa, Enceladus shows no sign of an oxygen-dominated exosphere that could otherwise be indicative of extreme surface and interior oxidation inhibiting the detectable survival and evolution of complex organics.

Cooper, John F.↗