Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Porosity gradient”

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

Segmentation and Classification of Fission as Pores in Reactor Irradiated Annular U–10Zr Metallic Fuel Using Machine Learning Models

Metallic fuels, particularly U—10Zr, are promising candidates for next-generation sodium-cooled fast reactors. Irradiation of nuclear fuels in reactors can lead to the formation of solid and gas fission product which subsequently forms microstructural pores, deteriorating fuel performance. Due to the massive amount of pores and complex phases formed, a quantitative description of fission gas pores is not yet available, preventing the development of microstructure-informed fuel performance modeling for fuel qualification. This paper applied a pre-trained deep learning model to ~10,260 high magnification scanning electron microscopy images. This method increased the accuracy of fission gas pore segmentation and allows statistical features to be extracted which cannot be achieved manually. A pre-trained decision tree model worked on the segemenation results and further classified the pores into different categories to produce a correlation between the pores, movement of lanthanides, and temperature gradient during irradiation. Finally, this paper emphasizes the potentials of machine learning models to accelerate fuel research, development, and qualification for advanced reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thermally Insulating Transparent Barrier (THINNER) coatings on single pane windows

Conventional silica aerogel monoliths can provide remarkable thermal insulation but the presence of large pores (> 30 nm) tends to scatter visible light and render the material opaque or translucent instead of transparent. In addition, they are prone to cracking during synthesis and handling which makes them difficult to integrate in products and in particular in window solutions. This project developed two new solgel synthesis methods using silica precursors or preformed silica nanoparticles and ambient drying to produce mesoporous organo-silica monoliths. The monoliths were (i) thermally insulating, (ii) optically transparent, (iii) flexible, and (iv) hydrophobic. They feature porosity ranging from 50% to 90% with narrow pore size distribution with pore less than 20 nm resulting in excellent optical clarity (haze < 2%) and very low thermal conductivity (< 30 mW/mK). Interestingly, not only porosity but also pore size and mass fractal dimension were found to affect the thermal conductivity of the mesoporous silica. The superior transparency of the monoliths was shown to be attributed to dependent scattering among silica nanoparticles. Flexibility was achieved through trimethylchlorosilane surface modification. Furthermore, process scale-up and integration of the ambigel monoliths into window solutions using optically clear adhesives were also demonstrated in 6”x6” double-pane windows. The aerogel and the window solution were shown to be durable to accelerated aging under UV, moisture, and/or temperature gradient and thermal cycling. Overall, the window solution achieved the technical performance and the cost target of $10/sqft set for the SHIELD program. However, scaling to industry relevant scale (> 10’x10’) remains a challenge due to requirements on the drying process and fume hood size and to the propensity of large aerogel slab to crack during drying.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1998

Rapid subsidence over oil fields measured by SAR

The Lost Hills and Belridge oil felds are in the San Joaquin Valley, California. The major oil reservoir is high porosity and low permeability diatomite. Extraction of large volumes from shallow depths causes reduction in pore pressure and subsequent compaction, forming a surface subsidence bowl. We measure this subsidence from space using interferometric analysis of SAR (Synthetic Aperture Radar) data collected by the European Space Agency Remote Sensing Satellites (ERS-1 and ERS-2). Maximum subsidence rates are as high as 40 mm in 35 days or > 400 mm/yr, measured from interferograms with time separations ranging from one day to 26 months. The 8- and 26-month interferograms contain areas where the subsidence gradient exceeds the measurement possible with ERS SAR, but shows increased detail in areas of less rapid subsidence. Synoptic mapping of subsidence distribution from satellite data powerfully complements ground-based techniques, permits measurements where access is difficult, and aids identification of underlying causes.

synthetic aperture radar remote sensing oil fields↗
NASA NTRS1985

Amorphous and porous ices in cometary nuclei

The heat flux into the cometary nucleus before and after perihelion and its influence on the nucleus are discussed. The surface temperature and the associated rate of evaporation are lower than usually calculated, and the effect is strongly dependent on the porosity of the nucleus. The surface temperature reaches a maximum after perihelion, in agreement with statistical studies of comae and tails. Delayed heating of the CO2-rich inclusions deep in the nucleus leads to outbursts and splittings as far as 9 AU after perihelion, also in agreement with observation. At sufficiently high temperatures, heat transport in porous cometary nuclei is controlled by vapor diffusion. The presence of CO2 vapor in the pores significantly increases the heat flow in mixed nuclei at temperatures above 140-150 C. Closed pores move up the radial thermal gradient so that the outside layer of a cometary nucleus should become denser.

Smoluchowski, R.↗
DOE OSTI2026

Tailoring Thermal and Mechanical Performance Through Multimaterial Laser Powder Directed Energy Deposition of Copper and 17-4PH Stainless Steel

This study investigates the additive manufacturing (AM) processing, microstructural evolution, and resulting mechanical and thermal properties of multimaterial components combining 17-4PH stainless steel and pure copper (Cu) fabricated via laser powder directed energy deposition (LP-DED). Conventional tooling steels exhibit limited thermal conductivity, significantly constraining production throughput in high-volume processes. Incorporating Cu, with its superior thermal conductivity, could significantly enhance tool performance, though Cu and steel present metallurgical incompatibilities when processed via AM. A systematic investigation was conducted across compositions ranging from 0 to 100 wt% Cu, revealing critical thresholds influencing solidification behavior, defect formation, microstructure, hardness, and thermal transport. Optical microscopy, electron backscatter diffraction (EBSD), hardness testing, and thermal conductivity measurements provided comprehensive process–structure–property correlations. Severe hot cracking occurred at low-Cu contents (6–25 wt%), aligning generally well with crack susceptibility modeling, with an unexpected discrepancy at 25 wt%. Porosity remained low (≥99% dense) throughout the compositional spectrum. EBSD analysis revealed a transformation from columnar martensitic structures at low-Cu contents to equiaxed FCC Cu-dominated structures at higher Cu concentrations, highlighting the complex microstructural transitions driven by Cu-induced changes in solidification and phase stability. Hardness decreased from 330 HV (pure 17-4PH) to 62 HV (pure Cu), consistent with microstructural changes. Concurrently, thermal conductivity improved substantially from 13.5 W/m K to 367.9 W/m K, emphasizing Cu’s dominant role in thermal transport. The findings highlight the feasibility of leveraging compositional gradients between 17-4PH and Cu to achieve tailored tooling with optimized thermal and mechanical performance.

17-4PH↗
DOE OSTI2022

Mechanical behavior of bimetallic stainless steel and gray cast iron repairs via directed energy deposition additive manufacturing

The utility of gray cast iron in engine components remains tied to the mechanical performance and cost. Repair and remanufacturing of castings offer economical and sustainable benefits; however, high thermal input from traditional fusion-based welding is unable to restore the original mechanical quality owing to brittle microstructures and porosity formed in situ. Directed energy deposition (DED) is an additive manufacturing method that has received considerable interest for repairs owing to the highly controllable nature of the process. Despite this, few works have connected the effect of DED parameters on actual interfacial strength. Consequently, distinct DED parameter combinations were identified to maximize the strength and fatigue life of the repaired cast iron. Further, high speed melt pool imaging and residual stress measurements are provided to aid in the understanding of the metallurgical quality and strength seen in these structures. In general, higher scanning speeds and lower thermal gradients promoted comparable tensile strength to that of the original gray cast iron. The results presented here provide a foundation to tune in the DED process to generate the required mechanical quality as a starting point for future process advancements.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multiscale mechanical design of the lightweight, stiff, and damage-tolerant cuttlebone: A computational study

Cuttlebone, the endoskeleton of cuttlefish, offers an intriguing biological structural model for designing low-density cellular ceramics with high stiffness and damage tolerance. Cuttlebone is highly porous (porosity ~93%) and lightweight (density less than 20% of seawater), constructed mainly by brittle aragonite (95 wt%), but capable of sustaining hydrostatic water pressures over 20 atmospheres and exhibits energy absorption capability under compression comparable to many metallic foams (~4.4 kJ/kg). Here, in this work, we computationally investigate how such remarkable mechanical efficiency is enabled by the multiscale structure of cuttlebone. Using the common cuttlefish, Sepia Officinalis, as a model system, we first conducted high-resolution synchrotron micro-computed tomography (µ-CT) and quantified the cuttlebone's multiscale geometry, including the 3D asymmetric shape of individual walls, the wall assembly patterns, and the long-range structural gradient of walls across the entire cuttlebone (ca. 38 chambers). The acquired 3D structural information enables systematic finite-element simulations, which further reveal the multiscale mechanical design of cuttlebone: at the wall level, wall asymmetry provides optimized energy absorption while maintaining high structural stiffness; at the chamber level, variation of walls (number, pattern, and waviness amplitude) contributes to progressive damage; at the entire skeletal level, the gradient of chamber heights tailors the local mechanical anisotropy of the cuttlebone for reduced stress concentration. Our results provide integrated insights into understanding the cuttlebone's multiscale mechanical design and provide useful knowledge for the designs of lightweight cellular ceramics.

36 MATERIALS SCIENCE↗
NASA NTRS1996

Improved Models for the Porous Surface with Passive Control

A computational investigation of afterbody flow using a passive control method is conducted. The passive control method consists of a porous surface placed over a plenum. The purpose of the passive control method is to exploit the adverse pressure gradient present in afterbody flow in an attempt to reduce boundary layer separation and afterbody drag. Four different porous wall models are used to model the transpiration velocity in the region of passive control. A three-dimensional, time-dependent, Reynolds-averaged, simplified Navier-Stokes solver, PAB3D, is used to simulate afterbody flow with and without passive control. Three afterbody configurations with boat-tail angles of 10, 20, and 30 deg. are used to obtain two-dimensional solutions with a freestream Mach number of 0.6 and nozzle pressure ratio of 6. The region of passive control was initially placed from 20-60% of the nozzle length. The effect of the porous placement and porous extent is also studied. Baseline (no porosity) two-dimensional solutions are qualitatively similar to experimental data but under-predict the magnitude of the pressure recovery. Results for the subsonic solutions show losses in the pressure recovery for some cases with passive control. Three-dimensional effects are also investigated and seen to be very significant. Three-dimensional baseline solutions, for both sub- and super-critical freestream Mach numbers, compare very favorably with the experimental data in comparison to the two-dimensional solution. Future work is required to examine three-dimensional afterbody flows with passive porosity.

Kraushaar, Sandra L.↗
NASA NTRS2001

Tortuosity measurement and the effects of finite pulse widths on xenon gas diffusion NMR studies of porous media

We have extended the utility of NMR as a technique to probe porous media structure over length scales of approximately 100-2000 microm by using the spin 1/2 noble gas 129Xe imbibed into the system's pore space. Such length scales are much greater than can be probed with NMR diffusion studies of water-saturated porous media. We utilized Pulsed Gradient Spin Echo NMR measurements of the time-dependent diffusion coefficient, D(t), of the xenon gas filling the pore space to study further the measurements of both the pore surface-area-to-volume ratio, S/V(p), and the tortuosity (pore connectivity) of the medium. In uniform-size glass bead packs, we observed D(t) decreasing with increasing t, reaching an observed asymptote of approximately 0.62-0.65D(0), that could be measured over diffusion distances extending over multiple bead diameters. Measurements of D(t)/D(0) at differing gas pressures showed this tortuosity limit was not affected by changing the characteristic diffusion length of the spins during the diffusion encoding gradient pulse. This was not the case at the short time limit, where D(t)/D(0) was noticeably affected by the gas pressure in the sample. Increasing the gas pressure, and hence reducing D(0) and the diffusion during the gradient pulse served to reduce the previously observed deviation of D(t)/D(0) from the S/V(p) relation. The Pade approximation is used to interpolate between the long and short time limits in D(t). While the short time D(t) points lay above the interpolation line in the case of small beads, due to diffusion during the gradient pulse on the order of the pore size, it was also noted that the experimental D(t) data fell below the Pade line in the case of large beads, most likely due to finite size effects.

Non-NASA Center↗
DOE OSTI2025

Optimization of direct air capture processes using reactive transport models of adsorption-desorption cycles

In this study, we develop and implement a reactive transport model in COMSOL Multiphysics® to address the challenges of direct air carbon capture. The model is validated against experimental data and used to simulate the cyclic steady state of the adsorption-desorption process. The optimization of this model is achieved through advanced trust-region methods integrated with Gaussian Processes. Key decision variables, including adsorption and desorption times, desorption temperature and pressure, input velocity, bed porosity, column length, and radius were optimized to minimize the capture cost. After optimization, a sensitivity analysis revealed the complex interplay between the decision variables and their effect on the specific energy and cost of removing the CO 2 . We optimized the capture cost while taking into account the trade-off between energy consumption and productivity. The resulting minimum capture cost was determined to be 265.2 $/t-CO 2 , which aligns with expected values reported in the literature. Numerical results suggest the effectiveness of the optimization strategies applied, and underscore the importance of simultaneous decision variable selection in improving the performance in direct air capture processes. We also extend the modeling approach to a 2D axisymmetric model to better visualize CO₂ uptake and temperature profiles, revealing significant radial gradients during the regeneration step. As a main drawback, this enhanced model comes with a computational cost approximately 40 times higher than that of the 1D model.

Adsorption-desorption process↗
NASA NTRS2013

Investigation of VEGGIE Root Mat

VEGGIE is a plant growth facility that utilizes the phenomenon of capillary action as its primary watering system. A cloth made of Meta Aramid fiber, known as Nomex is used to wick water up from a reservoir to the bottom of the plants roots. This root mat system is intended to be low maintenance with no moving parts and requires minimal crew interface time. Unfortunately, the water wicking rates are inconsistent throughout the plant life cycle, thus causing plants to die. Over-wicking of water occurs toward the beginning of the cycle, while under-wicking occurs toward the middle. This inconsistency of wicking has become a major issue, drastically inhibiting plant growth. The primary objective is to determine the root cause of the inconsistent wicking through experimental testing. Suspect causes for the capillary water column to break include: a vacuum effect due to a negative pressure gradient in the water reservoir, contamination of material due to minerals in water and back wash from plant fertilizer, induced air bubbles while using syringe refill method, and material limitations of Nomex's ability to absorb and retain water. Experimental testing will be conducted to systematically determine the cause of under and over-wicking. Pressure gages will be used to determine pressure drop during the course of the plant life cycle and during the water refill process. A debubbler device will be connected to a root mat in order to equalize pressure inside the reservoir. Moisture and evaporation tests will simultaneously be implemented to observe moisture content and wicking rates over the course of a plant cycle. Water retention tests will be performed using strips of Nomex to determine materials wicking rates, porosity, and absorptivity. Through these experimental tests, we will have a better understanding of material properties of Nomex, as well as determine the root cause of water column breakage. With consistent test results, a forward plan can be achieved to resolve the issue and give valuable insight for the next generation of VEGGIE.

Subbiah, Arun M.↗
DOE OSTI2022

In situ melt pool measurements for laser powder bed fusion using multi sensing and correlation analysis

Laser powder bed fusion is a promising technology for local deposition and microstructure control, but it suffers from defects such as delamination and porosity due to the lack of understanding of melt pool dynamics. To study the fundamental behavior of the melt pool, both geometric and thermal sensing with high spatial and temporal resolutions are necessary. This work applies and integrates three advanced sensing technologies: synchrotron X-ray imaging, high-speed IR camera, and high-spatial-resolution IR camera to characterize the evolution of the melt pool shape, keyhole, vapor plume, and thermal evolution in Ti–6Al–4V and 410 stainless steel spot melt cases. Aside from presenting the sensing capability, this paper develops an effective algorithm for high-speed X-ray imaging data to identify melt pool geometries accurately. Preprocessing methods are also implemented for the IR data to estimate the emissivity value and extrapolate the saturated pixels. Quantifications on boundary velocities, melt pool dimensions, thermal gradients, and cooling rates are performed, enabling future comprehensive melt pool dynamics and microstructure analysis. The study discovers a strong correlation between the thermal and X-ray data, demonstrating the feasibility of using relatively cheap IR cameras to predict features that currently can only be captured using costly synchrotron X-ray imaging. Such correlation can be used for future thermal-based melt pool control and model validation.

47 OTHER INSTRUMENTATION↗
NASA NTRSDate not supplied

Selection, Production, and Properties of Regolith Polymer Composites for Lunar Construction

NASA’s Artemis Program seeks to establish a long-term presence on the moon to enable scientific discoveries and utilization of lunar resources through public-private partnerships. Over the next decades, a lunar spaceport enabling regular transportation from and to Earth will need to be established to provide the services and facilities that are necessary to achieve this goal. Robotic construction technologies using in-situ materials must be developed to build up enabling infrastructure such as launch/landing pads, blast protection, power/communications towers, improved trafficability pathways and radiation protection shelters. Kennedy Space Center’s Granular Mechanics and Regolith Operations laboratory at Swamp Works has partnered with AI SpaceFactory to develop the architectural and structural design of an unpressurized shelter. The shelter, called Lina, is designed to protect astronauts and surface assets from radiation, meteoroid impact, moon quakes and thermal gradients. Lina’s structural design criteria and the resulting structure design are detailed in separate papers. A Fused Filament Fabrication (FFF) construction process using regolith polymer composites has been developed. This paper presents the material formulations and selection rationale for each of the composite components. Formulations include 70:30, 80:20, and 85:15 wt.% mixture ratios of lunar mare simulant Black Point-1 (BP-1): Polylactic Acid (PLA), 80:20 Lunar Highlands Simulant (LHS-1): PLA and an 80:20 BP-1:PLA formulation with a processing aid/compatibilizer additive. Test samples were printed in simulated lunar dirty thermal vacuum conditions (~-180 °C, ~10 -3 torr). The test environment evolved throughout the duration of the print process. A series of tests were performed to characterize the mixture ratios achieved for each formulation, the off-gassing products during vacuum printing, the strength properties, and porosity of printed products. The LHS-1: PLA formulation yielded an advantageous combination of properties and was used in a final test that additively constructed a sub-scaled Lina on regolith simulant in dirty vacuum conditions.

lunar infrastructure↗
NASA NTRSDate not supplied

Size Effects on Laser Powder Bed Fusion and Laser Powder Directed Energy Deposition GRCop-42 Alloy

Additive manufacturing (AM) plays a significant role in the fabrication of regeneratively cooled rocket engine combustion chambers, where internal channels are used to transport high pressure propellants to keep the engine from overheating. Cooling channels are designed with thin walls between the cooling and combustion sides, experiencing high thermal and mechanical stresses from large temperature gradients and internal pressures. Hence, selecting materials suited for this environment and maintaining a high-quality processing for this critical application is of utmost importance. GRCop-alloys (Cu-Cr-Nb) were developed as an ideal material for combustion chambers due to its high conductivity and high strength at temperature, being employed as a structural member for the hotwall with excellent heat transfer capabilities. GRCop-42 (Cu, 4 at% Cr, 2 at% Nb) has significantly matured using the Laser Powder Bed Fusion (L-PBF) process through microstructure, mechanical and thermophysical properties and evaluations using hot-fire testing. Additionally, Laser Powder Directed Energy Deposition (LP-DED) has demonstrated successful deposition and geometry samples using GRCop-42. While the maturation in properties has focused on bulk material properties from standard test specimens, these cannot be extrapolated to be used in thin walls, as different features (i.e., microstructure, porosity, and surface texture) have a greater effect on thin walls, especially on AM structures built using L-PBF and LP-DED. This presentation provides an overview of the effects introduced by reducing the thickness on the fabrication, microstructure, surface metrology and subsequent mechanical properties of GRCop-42 using L-PBF and LP-DED as well as future work planned for this research.

Gabriel Demeneghi↗
NASA NTRS2010

Radiation Environment and Surface Radiolytic Interactions at Mimas

Saturn's innermost principal moon Mimas shares the distinction with Europa at Jupiter of being the most irradiated icy moon in its respective planetary system, although the energetic electron energy flux at Mimas is forty times smaller than at Europa. High energy (> 10 MeV) proton fluxes are low in this moon's orbital corridor, likely since slowly diffusing protons from the weak but steady source of cosmic ray albedo neutron decay (CRAND) cannot accumulate without impacting the moon surface. Lower energy proton fluxes are also evidently suppressed in this orbital region. Plasma ion and electron fluxes are also low apparently due to cooling by interaction with E-ring dust and neutral gas from Enceladus. Due to energy-dependent effects of longitudinal gradient-curvature drift for the electrons, the trailing hemisphere is mainly irradiated by electrons at energies below 1 MeV that drift relative to Mimas in the prograde direction of orbital motion around Saturn, while higher energy electrons primarily impact the leading hemisphere. Plasma ions in the inner magnetosphere of Saturn are mainly pickup ions forming from the dissociation products of Enceladus plume water molecules, additionally including some contribution from photosputtering of the main rings, and do not introduce new elemental materials at Mimas via surface implantation from the corotating plasma. Thus the primary interaction at the surface is radiolytic chemistry induced in pure water ice by relatively deep penetration of the energetic electrons to millimeter and greater depths, as compared to the micron depths impacted by the corotating plasma ions. If surface erosion by sputtering from relatively low fluxes of the plasma and more energetic ions is indeed ineffective, then molecular products (OH, H2O2, 02, 03) of the radiolytic interactions may accumulate in the meters-deep impact regolith of the surface ices. An effect of regolith trapped gas accumulation could be to increase porosity and reduce thermal conductivity of the ice, potentially contributing to reported thermal anomalies from Cassini infrared map observations. Low amplitude of trailing-leading asymmetry in optical albedo and color maps at Mimas is suggestive of relative weakness of asymmetrical effects from low-energy ions. Greater induced asymmetries are expected and observed for the moons beyond Enceladus in the middle magnetospheric region of hot plasma ions at much greater fluxes than at Mimas. Low density (1.15 g/cc) of this moon indicates paucity of mineral salts and radiogenic heating to maintain subsurface liquids, so Enceladus-like cryovolcanism as a resurfacing process is unlikely despite closer proximity to Saturn, greater tidal forcing, and more intense surface irradiation than for Enceladus.

Cooper, J. F.↗
NASA NTRS2002

Diffusion NMR methods applied to xenon gas for materials study

We report initial NMR studies of (i) xenon gas diffusion in model heterogeneous porous media and (ii) continuous flow laser-polarized xenon gas. Both areas utilize the pulsed gradient spin-echo (PGSE) techniques in the gas phase, with the aim of obtaining more sophisticated information than just translational self-diffusion coefficients--a brief overview of this area is provided in the Introduction. The heterogeneous or multiple-length scale model porous media consisted of random packs of mixed glass beads of two different sizes. We focus on observing the approach of the time-dependent gas diffusion coefficient, D(t) (an indicator of mean squared displacement), to the long-time asymptote, with the aim of understanding the long-length scale structural information that may be derived from a heterogeneous porous system. We find that D(t) of imbibed xenon gas at short diffusion times is similar for the mixed bead pack and a pack of the smaller sized beads alone, hence reflecting the pore surface area to volume ratio of the smaller bead sample. The approach of D(t) to the long-time limit follows that of a pack of the larger sized beads alone, although the limiting D(t) for the mixed bead pack is lower, reflecting the lower porosity of the sample compared to that of a pack of mono-sized glass beads. The Pade approximation is used to interpolate D(t) data between the short- and long-time limits. Initial studies of continuous flow laser-polarized xenon gas demonstrate velocity-sensitive imaging of much higher flows than can generally be obtained with liquids (20-200 mm s-1). Gas velocity imaging is, however, found to be limited to a resolution of about 1 mm s-1 owing to the high diffusivity of gases compared with liquids. We also present the first gas-phase NMR scattering, or diffusive-diffraction, data, namely flow-enhanced structural features in the echo attenuation data from laser-polarized xenon flowing through a 2 mm glass bead pack. c2002 John Wiley & Sons, Ltd.

NASA Discipline Life Sciences Technologies↗
DOE OSTI2023

Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas

Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.

58 GEOSCIENCES↗
NASA NTRS2002

The narrow pulse approximation and long length scale determination in xenon gas diffusion NMR studies of model porous media

We report a systematic study of xenon gas diffusion NMR in simple model porous media, random packs of mono-sized glass beads, and focus on three specific areas peculiar to gas-phase diffusion. These topics are: (i) diffusion of spins on the order of the pore dimensions during the application of the diffusion encoding gradient pulses in a PGSE experiment (breakdown of the narrow pulse approximation and imperfect background gradient cancellation), (ii) the ability to derive long length scale structural information, and (iii) effects of finite sample size. We find that the time-dependent diffusion coefficient, D(t), of the imbibed xenon gas at short diffusion times in small beads is significantly affected by the gas pressure. In particular, as expected, we find smaller deviations between measured D(t) and theoretical predictions as the gas pressure is increased, resulting from reduced diffusion during the application of the gradient pulse. The deviations are then completely removed when water D(t) is observed in the same samples. The use of gas also allows us to probe D(t) over a wide range of length scales and observe the long time asymptotic limit which is proportional to the inverse tortuosity of the sample, as well as the diffusion distance where this limit takes effect (approximately 1-1.5 bead diameters). The Pade approximation can be used as a reference for expected xenon D(t) data between the short and the long time limits, allowing us to explore deviations from the expected behavior at intermediate times as a result of finite sample size effects. Finally, the application of the Pade interpolation between the long and the short time asymptotic limits yields a fitted length scale (the Pade length), which is found to be approximately 0.13b for all bead packs, where b is the bead diameter. c. 2002 Elsevier Sciences (USA).

NASA Discipline Life Sciences Technologies↗