Engineering Papers⌕ Search

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.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Vertical transport of water in the Martian boundary layer

We are continuing our examination of the transport of H2O through the martian boundary layer, and we have written a one-dimensional numerical model of the exchange of H2O between the atmosphere and subsurface of Mars through the planetary boundary layer (PBL). Our goal is to explore the mechanisms of H2O exchange, and to elucidate the role played by the regolith in the local H2O budget. The atmospheric model includes effects of Coriolis, pressure gradient, and frictional forces for momentum, as well as radiation, sensible heat flux, and advection for heat. The model differs from Flasar and Goody by use of appropriate Viking-based physical constants and inclusion of the radiative effects of atmospheric dust. We specify the pressure gradient force or compute it from a simple slope model. The subsurface model accounts for conduction of heat and diffusion of H2O through a porous adsorbing medium in response to diurnal forcing. The model is initialized with depth-independent H2O concentrations (2 kg M(exp -3)) in the regolith, and a dry atmosphere. The model terminates when the atmospheric H2O column abundance stabilizes at 0.1 percent per sol.

Zent, Aaron P.↗
NASA NTRS2019

Subsurface Cl-bearing salts as potential contributors to recurring slope lineae (RSL) on Mars

We report laboratory experimental results that support a brine-related hypothesis for the recurring slope lineae (RSL) on Mars in which the subsurface Cl-salts, i.e., hydrous chlorides and oxychlorine salts (HyCOS) are the potential source materials. Our experiments revealed that within the observed RSL temperature window TRSL (250–300 K), the deliquescence of HyCOS could occur in relative humidity ranges (RH≥22%–46%) much lower than those for hydrous (Mg, Fe)-sulfates (RH≥75%–96%). In addition, we demonstrated that the RH values kept by common HyCOS and hydrous sulfates in enclosures have a general trend as RH(sulfates) > RH(perchlorates) > RH(chlorides) (with same type of cation) in wide T range. It means that the required RH range for a Cl-bearing salt to deliquescence can be satisfied by a co-existing salt of different type, e.g., in the subsurface layers of mixed salts on Mars. Furthermore, we found a strong temperature dependence of the deliquescence rates for all tested HyCOS, e.g., a duration of 1–5 sols for all HyCOS at the high end (300 K) of TRSL, and of 20–70 sols for all tested HyCOS (except NaClO4·H2O) at the low end (250 K) of TRSL, which is consistent with the observed seasonal behavior of RSL on Mars. From a mass-balance point of view, the currently observed evidences on Mars do not support a fully-brine-wetted track model, thus we suggest a brine-triggered granular-flow model for the most RSL. Considering the recurrence of RSL in consecutive martian years, our experimental results support the rehydration of remnant HyCOS layers during a martian cold season through H2O vapor-to-salt direct interaction. We found that the evidences of HyCOS rehydration under Mars relevant P-T-RH conditions are detectable in a few minutes by in situ Raman spectroscopy. This rehydration would facilitate the recharge of H2O back into subsurface HyCOS, which could serve as the source material to trigger RSL in a subsequent warm season. The major limiting factor for this rehydration is the H2O supply, i.e., the H2O vapor density carried by current Mars atmospheric circulation and the diffusion rate of H2O vapor into the salt-rich subsurface in a cold season. In a worst-case scenario, these H2O supplies can support a maximum increase of hydration degrees of two for totally dehydrated HyCOS, whereas the full rehydration of subsurface HyCOS layers can be easily reached during a>30° obliquity period that has H2O vapor density 10× to 20× times the value of current obliquity. Overall, our results imply the existence of a large amount of Cl-bearing salts in the subsurface at RSL sites.

Alian Wang↗
NASA NTRS1992

Regional variations in the stability and diffusion of water-ice in the Martian regolith

Geologic evidence suggests subsurface water-ice has played an important role in the formation of Martian landforms. Forms of mass-wasting such as debris aprons and flow patterns on valley floors suggest creep deformation of ice-laden soil, while thermokarst and chaotic terrain suggest once extensive deposits of ground ice that were later removed. The global distribution of ice-related morphology was mapped. The mapping showed regional variation, in both latitude and longitude, in the distribution of debris aprons, concentric fill craters, and 'softened' crater profiles.

Mellon, Michael T.↗
DOE OSTI2024

Utah FORGE: Geochemical Data for Cold Groundwaters and Produced Geothermal Fluids

Geochemical data for cold groundwaters and produced geothermal fluids around the Utah FORGE site. The data is compiled into four tables in the attached Excel File. Table 1 is a compilation of compositions (anions, cations, weak acids, oxygen, hydrogen, and carbon isotopes) for cold groundwaters and produced geothermal waters in the Milford valley, Utah. Table 2 is a compilation of noble gas (He, Ne, Ar) and He and Ne isotopic compositions for cold groundwaters and produced geothermal waters in the Milford valley, Utah. Table 3 provides values for calculated advective and diffusive fluxes of helium. Table 4 provides values of calculated subsurface stored heat between the Opal Mound fault and the Utah FORGE site, which are related to volumes of recently solidified magmatic heat sources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Time-Resolved Neutron Imaging for Hydrogen Uptake in Subsurface Lithologies

Geologic hydrogen production and underground storage are increasingly important for meeting rising energy demands while providing clean-combustion advantages. However, hydrogen’s high diffusivity and propensity for leakage through porous media necessitate direct evaluation of its transport behavior in subsurface materials. Whereas X-ray microcomputed tomography (μCT) studies often employ contrast agents or surrogate gases, this study leverages neutron transmission radiography/CT to observe hydrogen migration in situ. This work represents the first demonstration of real-time neutron radiography of hydrogen migration in reservoir and caprock lithologies. Cylindrical cores of Indiana limestone, Amherst Gray sandstone, and Tumey shale were subjected to constant-pressure hydrogen charging and scanned in real time using high-resolution neutron radiography. Results indicate immediate hydrogen infiltration in sandstone and limestone, with homogeneous distribution detected throughout their pore structure. In contrast, hydrogen remained largely absent from fine-grained shale under the same pressure, except in an apparently localized fracture zone, where neutron signatures confirmed the presence of hydrogen. Subsequent neutron CT of the sandstone sample, using image subtraction against an uncharged reference, corroborated hydrogen distribution patterns. Even under lowpressure, single-phase conditions, distinct neutron imaging signatures of hydrogen were achieved. These preliminary findings underscore the potential of neutron imaging for advancing subsurface hydrogen migration research.

08 HYDROGEN↗
NASA NTRS1993

Equatorial ground ice on Mars: Steady-state stability

Current Martian equatorial surface temperatures are too warm for water ice to exist at the surface for any appreciable length of time before subliming into the atmosphere. Subsurface temperatures are generally warmer still and, despite the presence of a diffusive barrier of porous regolith material, it has been shown by Smoluchowski, Clifford and Hillel, and Fanale et al. that buried ground ice will also sublime and be lost to the atmosphere in a relatively short time. We investigate the behavior of this subliming subsurface ice and show that it is possible for ice to maintain at a steady-state depth, where sublimation and diffusive loss to the atmosphere is balanced by resupply from beneath by diffusion and recondensation of either a deeper buried ice deposits or ground water. We examine the behavior of equatorial ground ice with a numercial time-marching molecular diffusion model. In our model we allow for diffusion of water vapor through a porous regolith, variations in diffusivity and porosity with ice content, and recondensation of sublimed water vapor. A regolith containing considerable amounts of ice can still be very porous, allowing water vapor to diffuse up from deeper within the ice layer where temperatures are warmer due to the geothermal gradient. This vapor can then recondense nearer to the surface where ice had previously sublimed and been lost to the atmosphere. As a result we find that ice deposits migrate to find a steady-state depth, which represents a balance between diffusive loss to the atmosphere through the overlying porous regolith and diffusive resupply through a porous icy regolith below. This depth depends primarily on the long-term mean surface temperature and the nature of the geothermal gradient, and is independent of the ice-free porosity and the regolith diffusivity. Only the rate of loss of ground ice depends on diffusive properties.

Mellon, Michael T.↗
DOE OSTI2026

Adsorption, charge transfer and a coverage-driven transition of alkali metals on rutile TiO 2 (110)

The interaction of alkali metals with metal oxide surfaces is central to tuning surface reactivity in heterogeneous catalysis and photocatalysis. Here we present a comprehensive DFT+U study of the adsorption of alkali metals (Li, Na, K, Rb, Cs) on the (110) surface of rutile TiO 2 . At low coverage (θ = 1/8), all alkali metals bind preferentially to bridging oxygen sites with adsorption energies in the range −4.06 to −3.33 eV, transferring nearly one full electron (0.92–0.99 |e|) to the substrate and inducing Ti 4+ → Ti 3+ reduction. The excess charge localizes preferentially at subsurface Ti sites in the form of small polarons. Diffusion barriers indicate facile motion along bridging-oxygen rows, whereas inter-row hopping is strongly hindered. Coverage effects were examined systematically for potassium: adsorption energy and charge transfer per K atom decrease monotonically with increasing θ. Strikingly, a sharp energy discontinuity occurs between θ = 4/8 and θ = 5/8 (ΔE ≈ 1 eV per atom), which we identify as a coverage-driven structural transition arising from steric packing constraints and enhanced K–K electrostatic repulsion once every (1×1) surface cell is occupied. This structural transition perfectly correlates with a dramatic drop in the work function down to an ultra-low minimum of 0.84 eV at θ=5/8, followed by a metallization- driven recovery at higher coverages. Ab initio molecular dynamics simulations confirm zigzag K arrangements at moderate coverage (θ = 1/3), while at high coverage (θ = 2/3) short-range K–K correlations emerge without long-range order. These results provide atomistic insight into the structure–activity relationships underlying alkali promotion effects on oxide-supported catalysts.

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

Electron Trapping and Ion Leaching at the Li-Modified Quartz–Water Interface

Electrons and ions at metal oxide-water interfaces have a critical role in many phenomena and applications, thus making their properties of considerable interest. We here investigate the behavior of a lithium impurity at the water interface with Li-doped SiO 2 using hybrid density functional based ab initio molecular dynamics simulations. We find that the excess electron donated by the lithium dopant localizes on a surface silicon atom and is partially solvated at the aqueous interface. While the position of the excess electron does not change during our ~60 ps simulation, the lithium cation diffuses from its initial position in the SiO 2 subsurface toward the interface and eventually leaches out of the oxide surface and becomes solvated by four water molecules, forming an aqueous Li+- electron complex at the interface. The degree of interaction between the localized electron and the lithium ion controls not only the energy level of the excess electron but also the extent of structural distortion and hydrogen-bonding at the interface. The results of our study can be relevant for understanding the initial stages of glass corrosion, where alkali ions are leached from the surface region by interactions with water.

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

Phenethylammonium Functionalization Enhances Near-Surface Carrier Diffusion in Hybrid Perovskites

Understanding semiconductor surface properties and manipulating them chemically are critical for improving their performance in optoelectronic devices. Hybrid halide perovskites have emerged as an exciting class of highly efficient solar materials; however, their device performance could be limited by undesirable surface properties that impede carrier transport and induce recombination. Here we show that surface functionalization of methylammonium lead iodide (MAPbI 3 ) perovskite with phenethylammonium iodide (PEAI), a commonly employed spacer cation in two-dimensional perovskites, can enhance carrier diffusion in the near-surface regions and reduce defect density by more than one order of magnitude. Using transient transmission and reflection microscopy, we selectively imaged the transport of the carriers near the (001) surface and the bulk carriers in single-crystal MAPbI 3 microplates. The surface functionalization increases the diffusion coefficient of the carriers in a 40-nm subsurface region from ~ 0.6 cm 2 s -1 to ~ 1 cm 2 s -1 , similar to the value for bulk carriers. These results suggest the PEA ligands are effective in suppressing surface defect and phonon scattering and shed light on the mechanisms for enhancing photophysical properties and improving solar cell efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

Dark Halo Craters and the Thickness of Grooved Terrain on Ganymede

Dark halo craters on grooved terrain on Ganymede represent potential probes of the subsurface geology. These craters are surrounded by a broad, diffuse low albedo annulus, or halo and have apparently lost the bright rim and ray deposits associated with very young craters. The recognition of dark halo craters, almost all greater than 12 km in diameter, indicates that material darker and redder than grooved terrain material forms a stratigraphy horizon at about 1 km depth the present surface of grooved terrain in Uruk Sulcus. This material is most likely downdropped cratered terrain material, which only larger craters have excavated into. This is most consistent with extensional tectonic models for grooved terrain formation whereby the lithosphere is stretched and blocks of ancient cratered terrain are downdropped along bounding faults, and subsequently resurfaced by a shallow layer of relatively clean icy material. Evidence for the preservation of cratered terrain material at shallow depth beneath grooved terrain material poses great difficulties for alternative grooved and smoothed terrain formation mechanisms.

Schenk, P. M.↗
DOE OSTI2020

Iodine effective diffusion coefficients through volcanic rock: Influence of iodine speciation and rock geochemistry

Accurate prediction of the subsurface transport of iodine species is important for the assessment of long-term nuclear waste repository performance, as well as monitoring compliance with the Comprehensive Nuclear-Test-Ban Treaty, given that radioiodine decays into radioxenon. However, the transport of iodine through intact geologic media is not well understood, compromising our ability to assess risk associated with radioiodine migration. The current study's goal is to quantify the matrix diffusion of iodine species through saturated volcanic rock, with particular attention paid to the redox environment and potential speciation changes. Here, diffusion experiments were run for iodide through lithophysae-rich lava, lithophysae-poor lava, and welded tuff, whereas iodate diffusion was studied through welded tuff. Iodine transport was compared with a conservative tracer, HDO, and effective diffusion coefficients were calculated. Likely due to a combination of size and anion exclusion effects, iodine species diffused more slowly than the conservative tracer through all rock types tested. Furthermore, oxidation of iodide to iodate was observed in the lithophysae-poor lava, affecting transport. Results provide much needed data for subsurface transport models that predict radioiodine migration from underground sources, and indicate the pressing need for geochemical and redox interactions to be incorporated into these models.

58 GEOSCIENCES↗
DOE OSTI2023

Machine Learning Predictions of Simulated Self-Diffusion Coefficients for Bulk and Confined Pure Liquids

Diffusion properties of bulk fluids have been predicted using empirical expressions and machine learning (ML) models, suggesting that predictions of diffusion also should be possible for fluids in confined environments. The ability to quickly and accurately predict diffusion in porous materials would enable new discoveries and spur development in relevant technologies such as separations, catalysis, batteries, and subsurface applications. Here in this work, we apply artificial neural network (ANN) models to predict the simulated self-diffusion coefficients of real liquids in both bulk and pore environments. The training data sets were generated from molecular dynamics (MD) simulations of Lennard-Jones particles representing a diverse set of 14 molecules ranging from ammonia to dodecane over a range of liquid pressures and temperatures. Planar, cylindrical, and hexagonal pore models consisted of walls composed of carbon atoms. Our simple model for these liquids was primarily used to generate ANN training data, but the simulated self-diffusion coefficients of bulk liquids show excellent agreement with experimental diffusion coefficients. ANN models based on simple descriptors accurately reproduced the MD diffusion data for both bulk and confined liquids, including the trend of increased mobility in large pores relative to the corresponding bulk liquid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

The distribution of ground ice on Mars

A wealth of geologic evidence indicates that subsurface water ice has played an important role in the evolution of Martian landforms. Theoretical models of the stability of ground ice show that in the near-surface regolith ice is currently stable at latitudes poleward of about +/- 40 deg and below a depth of a few centimeters to a meter, with some variations with longitude. If ice is not previously present at a particular location where it is stable, atmospheric water will diffuse into the regolith and condense as ice, driven by the annual subsurface thermal oscillations. The lower boundary of this ice deposit is found to occur at a depth (typically a few meters) where the annual thermal oscillations give way to the geothermal gradient. In the equatorial regions near-surface ice is currently not stable, resulting in the sublimation of any existing ice and subsequent loss to the atmosphere. However, subliming ice might be maintained at a steady-state depth, where diffusion and loss to the atmosphere are balanced by resupply from a possible deeper source of water (either deeper, not yet depleted, ice deposits or ground water). This depth is typically a few tens to hundreds of meters and depends primarily on the surface temperature and the nature of the geothermal gradient, being deeper for a higher surface temperature and a lower geothermal gradient. Such an equatorial deposit is characterized by the regolith ice content being low nearer the surface and increasing with depth in the deposit. Oscillations in the orbit will affect this picture of ground ice in two ways: by causing periodic changes in the pattern of near-surface stability and by producing subsurface thermal waves that may be capable of driving water ice deeper into the regolith.

Mellon, M. T.↗
NASA NTRS1991

Endotaxial growth of CoSi2 within (111) oriented Si in a molecular beam epitaxy system

A new mode of growth is reported in which buried metallic layers can be fabricated within a single-crystal semiconductor through preferential subsurface growth on previously-grown 'seed' regions. The deposition of Co at 800 C at a rate of 0.01 nm/s on (111) Si substrates containing buried CoSi2 columns 40-100 nm below the Si surface results in the growth and coalescence of these subsurface columns. The formation of a CoSi2 layer on the Si surface is suppressed by this growth mode. It is proposed that the high diffusion rate of Co at 800 C, coupled with the high growth rate of CoSi2 at the subsurface columns, is responsible for this preferred 'endotaxial' growth mode. This growth technique was used to produce a continuous buried single-crystal layer of CoSi2 under a single-crystal Si capping layer.

George, T.↗
DOE OSTI2025

Probing the Surface Chemistry of Lithium Nitridation

Chemical synthesis of Li 3 N through lithium nitridation has potential to advance rechargeable battery and nitrogen fixation technology. However, studies of the conditions for forming Li 3 N on the lithium surface via nitrogen gas exposure report contradictory findings, such as the spontaneous reaction of Li with pure N 2 , the impossibility of forming Li 3 N through pure Li and N 2 interaction, the requirement of trace H 2 O to catalyze the reaction, and evidence to the contrary. In this study, ambient pressure X-ray photoelectron spectroscopy (APXPS) was applied to evaluate the in situ chemical evolution of the lithium metal surface under nitrogen gas up to 800 mTorr. At pressures ≤10 mTorr, no Li 3 N was detected. At higher pressures, surface Li 3 N rapidly reacts with trace CO 2 . Additionally, because metallic lithium is readily oxidized by trace gases, the atomic nitrogen concentration of the lithium surface remains below 2%. When nitridation follows oxidation by O 2 gas, CO 2 gas, or H 2 O vapor, surface Li 3 N formation is inhibited. These results suggest that nitrogen gas can diffuse through the oxidized lithium metal surface to react with subsurface metallic lithium.

Binding energy↗
DOE OSTI2024

Subsurface Redox Interactions Regulate Ebullitive Methane Flux in Heterogeneous Mississippi River Deltaic Wetland

Abstract As interfaces connecting terrestrial and ocean ecosystems, coastal wetlands develop temporally and spatially complex redox conditions, which drive uncertainties in greenhouse gas emission as well as the total carbon budget of the coastal ecosystem. To evaluate the role of complex redox reactions in methane emission from coastal wetlands, a coupled reactive‐transport model was configured to represent subsurface biogeochemical cycles of carbon, nitrogen, and sulfur, along with production and transport of multiple gas species through diffusion and ebullition. This model study was conducted at multiple sites along a salinity gradient in the Barataria Basin at the Mississippi River Deltaic Plain. Over a freshwater to saline gradient, simulated total flux of methane was primarily controlled by its subsurface production and consumption, which were determined by redox reactions directly (e.g., methanogenesis, methanotrophy) and indirectly (e.g., competition with sulfate reduction) under aerobic and/or anaerobic conditions. At fine spatiotemporal scales, surface methane fluxes were also strongly dependent on transport processes, with episodic ebullitive fluxes leading to higher spatial and temporal variability compared to the gradient‐driven diffusion flux. Ebullitive methane fluxes were determined by methane fraction in total ebullitive gas and the frequency of ebullitive events, both of which varied with subsurface methane concentrations and other gas species. Although ebullition thresholds are constrained by local physical factors, this study indicates that redox interactions not only determine gas composition in ebullitive fluxes but can also regulate ebullition frequency through gas production.

03 NATURAL GAS↗
NASA NTRS2012

Modeling Radar Scatter from Icy and Young Rough Lunar Craters

For lunar orbital synthetic aperture radars, such as the Chandrayaan Mini-RF operating at S- band (13-cm) wavelength and the Lunar Reconnaissance Orbiter Mini-RF operating at S- band and X-band (3-cm) wavelengths, it is important to understand the radar backscattering characteristics of the icy and young, rough craters. Assuming a mixing model consisting of diffuse and quasi-specular scattering components, we have modeled the opposite-sense circular (OC) and same-sense circular (SC) backscattering characteristics. The specular component, consisting of only OC echoes, represents the echoes from the surface and subsurface layers that are oriented perpendicular to the radar's line-of-sight. The diffuse component, consisting of both SC and OC echoes, represents the echoes associated with either rocks or ice. Also, diffuse echoes have backscatter that is proportional to the cosine of the incidence angle. We modeled how these two (specular and diffuse) radar scattering components could be modulated by factors such as surface roughness associated with young craters. We also modeled how ice radar scattering components could be modulated by a thin regolith covering, and/or by the situation where ice occupies small patches within a larger radar pixel. We tested this modeling by examining 4 nonpolar craters and 12 polar craters using LRO Mini-RF data. Results indicate that icy and young rough craters can be distinguished based upon their SC enhancements (Alpha) and OC enhancements (Gamma). In addition, we also examined the craters that have unusual circular polarization ratios (CPRs) that likely result from a double bounce mode of scattering. Blocky fresh craters, icy craters, and craters exhibiting double bounce scattering can be separated based on the values of Alpha, Gamma, the ratio of Alpha/Gamma and the weighted sum of Alpha and Gamma.

lunar craters↗