SEARCH · Engineering Papers

Results for “subsurface system”

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 1,189 records · Page 66

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Numerical evaluation of photosensitive tracers as a strategy for separating surface and subsurface transient storage in streams

For this work, we numerically evaluated photosensitive tracers as a potential strategy for separating the effects of surface and hyporheic storage zones (SSZs and HSZs, respectively) on stream corridor transport. Correctly separating HSZ and SSZ effects is critical to estimating the hydro-biogeochemical function of a stream because HSZs and SSZs expose solutes to significantly different biogeochemical conditions, like sunlight exposure, microbial processes, and oxygen availability. Our numerical experiments used a multiscale river-corridor transport model implemented in the ATS code, which accommodates multiple storage zones with distinct travel time distributions and biogeochemical reactions. For parameter inferences, we used Bayesian inverse modeling. We found that breakthrough curves for photo-decaying tracers from day and night injection can delineate surface and hyporheic transient storage contributions, but only when interpreted jointly through a two-storage zone model. Numerical experiments that used only daytime injection or interpreted breakthrough curves with a single storage zone model yielded good fit to breakthrough curves, but parameter estimates were biased and controlling processes misattributed, examples of good model fits for the wrong reasons. Using those biased parameter estimates in reactive transport simulations resulted in significantly different projections of denitrification, which underscores the potential for stream function to be mischaracterized if tracer tests are interpreted through an inappropriately simplified model for transient storage. More generally, this study highlights the role of modeling in evaluating the experimental design and identifying the potential of system mischaracterization and its implications.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Estimating biofuel contaminant concentration from 4D ERT with mixing models

Here, we present the results of a lab-scaled feasibility study to assess the performance of electrical resistivity tomography for detection, characterization, and monitoring of fuel grade ethanol releases to the subsurface. Further, we attempt to determine the concentration distribution of the ethanol from the electrical resistivity tomography data using mixing-models. Ethanol is a renewable fuel source as well as an oxygenate fuel additive currently used to replace the known carcinogen methyl tert-butyl ether; however, ethanol is preferentially biodegraded and a cosolvent. When introduced to areas previously impacted by nonethanol-based fuels, it will facilitate the persistence of carcinogenic fuel compounds like benzene and ethylbenzene, as well as remobilize them to the ground water. These compounds would otherwise be retained in the soil column undergoing active or passive remediation processes such as soil vapor extraction or natural attenuation. Here, we introduce ethanol to a saturated Ottawa sand in a tank instrumented for four-dimensional geoelectrical measurements. Forward model results suggest pure phase ethanol released into a water saturated silica sand should present a detectable target for electrical resistivity tomography relative to a saturated silica sand only. We observe the introduction of ethanol to the closed hydraulic system and subsequent migration over the duration of the experiment. One-dimensional and three–dimensional temporal data are assessed for the detection, characterization, and monitoring of the ethanol release. Results suggest one-dimensional geoelectrical measurements may be useful for monitoring a release, while three-dimensional geoelectrical field imaging would be useful to characterize, monitor, and design effective remediation approaches for an ethanol release, assuming field conditions do not preclude the application of geoelectrical methods. We then attempt to use predictive mixing models to calculate the distribution of ethanol concentration within the measurement domain. For this study we examine four different models: a nested parallel mixing model, a nested cubic mixing model, the complex refractive index model (CRIM), and the Lichtenecker-Rother (L-R) model. The L-R model, modified to include an electrical formation factor geometry term, provided the best agreement with expected EtOH concentrations.

09 BIOMASS FUELS↗

DOE OSTI2025

Unraveling Hydrogen Induced Geochemical Reaction Mechanisms through Coupled Geochemical Modeling and Machine Learning

Underground hydrogen storage (UHS) provides a promising large-scale, long-term energy storage solution. A reasonable recovery of stored hydrogen is critical for a successful storage scheme. However, in subsurface reservoirs hydrogen is subject to active geochemical reactions that might result in hydrogen loss. In this study, we implemented a geochemical modeling approach coupled with an unsupervised machine learning technique called non-negative matrix factorization (NMF) to unravel the complex brine-rock-H 2 geochemical processes responsible for hydrogen losses, with particular focus on sulfate reduction reactions. NMF is applied to modeled mineral evolution and fluid component profiles to retrieve profiles that can be interpreted to more easily assess competing processes. NMF decouples simulated competing equilibrium reactions. This facilitates separation of overlapping reaction profiles from redox processes, dissolution fronts, and secondary precipitation while considering the effects of simulation parameters such as salinity, temperature, and total H 2 pressure. NMF successfully discriminates these competing effects in nonlinear ways, allowing robust interpretation. In addition, NMF reveals subtle coupled mineral associations and reaction fronts that are invisible to conventional model analysis. This integrated approach strengthens the conceptual understanding of complex nonlinear hydrogen-brine-rock interactions and advances geochemical research on UHS systems to resolve complexities in modeled geochemical systems without the need for direct experiments or prior knowledge. Furthermore, this study highlights the efficacy of combining geochemical modeling with machine learning techniques to enhance the interpretability of the intricate geochemical simulation output through deciphering the overlapping reaction path that cannot be achieved only using conventional analysis of geochemical models alone.

08 HYDROGEN↗

DOE OSTI2023

Quantifying model prediction sensitivity to model-form uncertainty

Computational and mathematical models are essential to understanding complex systems and phenomena. However, when developing such models, limited knowledge and/or resources necessitates the use of simplifying assumptions. It is therefore crucial to quantify the impact of such simplifying assumptions on the reliability and accuracy of resulting model predictions. This work develops a first-of-its-kind approach to quantify the impact of physics modeling assumptions on predictions. Here, we leverage the emerging field of model-form uncertainty (MFU) representations, which are parameterized modifications to modeling assumptions, in combination with grouped Sobol’ indices to quantitatively measure an assumption’s importance. Specifically, we compute the grouped Sobol’ index for the MFU representation’s parameters as a single importance measure of the assumption for which the MFU representation characterizes uncertainty. To ensure this approach is robust to the subjective choice of how to parameterize a MFU representation, we establish bounds for the difference between sensitivity results for two different MFU representations based on differences in model prediction statistics. The capabilities associated with this approach are demonstrated on three exemplar problems: an upscaled subsurface contaminant transport problem, ablation modeling for hypersonic flight, and nuclear waste repository modeling. We found that our grouped approach is able to assess the impact of modeling assumptions on predictions and offers computational advantages over classical Sobol’ index computation while providing more interpretable results.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Time-series dissolved oxygen, other bigeochemically-relevant analytes, and pressure gradients associated with the manuscript “Dissolved oxygen sensor in an automated hyporheic sampling system reveals biogeochemical dynamics”

This dataset contains time-series data from a vertical profile within the bed and banks of the Columbia river near Richland, WA. Water was sampled through 3 small tubes embedded in the sediment at 50,100, and 200 cm below the sediment-water interface. The goal of this study was to observe the correlations between hydraulic drivers and biogeochemical responses. The results of this study are published in the manuscript “Dissolved oxygen sensor in an automated hyporheic sampling system reveals biogeochemical dynamics”. The file types included in the data package are all time-series spreadsheet data, including hydraulic head gradients, physical parameters (temperature, pressure, SpC (specific conductivity)), and biogeochemical parameters (dissolved oxygen, pH, NO3 (nitrate) and ORP (oxidation-reduction potential)).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Optimizing Alabama’s CO 2 Storage in Shelby County (Project OASIS): Task 8 Deliverable 8.0 - Commercialization Plan

The “Optimizing Alabama’s CO 2 Storage in Shelby County: Project OASIS” CarbonSAFE Phase II Project seeks to build on regional data sets that demonstrate that the subsurface within Shelby County, Alabama has the potential to store commercial volumes of CO 2 safely, permanently, and economically. This work builds on the initiatives of the Southeast Regional Carbon Utilization and Storage Acceleration Partnership (SECARB-USA, DE-FE0031830) that identified nearly 500 million metric tonnes of CO 2 emitted on an annual basis from industrial facilities that are not collocated with prospective storage geology (i.e., the Coastal Plain of the Southeastern U.S. in this context). This observation suggests costly investments in connective infrastructure (e.g., pipelines) or exploratory well drilling campaigns are required to identify CO 2 storage opportunities in under-developed areas. While not traditionally thought of for saline storage, these studies suggest that storage prospects in the Valley and Ridge Province of the Appalachian Fold and Thrust Belt occur in relatively flat lying structural panels between thrust faults. For the Project OASIS region, available geologic studies related to hydrocarbon exploration suggest that Cambro-Ordovician carbonates and Cambrian clastic units offer multiple potential storage intervals, and that regional confining systems are present, such as the tectonically thickened Floyd-Parkwood Shale.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2022

The ParFlow Sandtank: An interactive educational tool making invisible groundwater visible

Physical aquifer models are a highly effective teaching tool for hydrology education, however they come with inherent limitations that include the high cost to purchase, the static configuration of the model materials, the time required to visualize hydrogeological phenomena, and the effort to reset and clean them over time. To address these and other limitations, we have developed an interactive computer simulation of a physical aquifer model called the ParFlow Sandtank. In this gamified interface, users run the simulation using a familiar web-app like interface with sliders and buttons while learning real hydrologic concepts. Our user interface allows participants to dive into the world of hydrology, understanding assumptions about model parameters such as hydraulic conductivity, making decisions about inputs to groundwater aquifer systems such as pumping rates, visualizing outputs such as stream flow, transport, and saturation, and exploring various factors that impact real environmental systems such as climate change. The ParFlow Sandtank has already been used in a variety of educational settings with more than 9,000 users per year, and we feel this emerging educational tool can be used broadly in educational environments and can be scaled-up to provide greater accessibility for students and educators. Here we present the capabilities and workflow of the ParFlow Sandtank, two use cases, and additional tools and custom templates that have been developed to support and enhance the reach of the ParFlow Sandtank.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

The stable isotopes of natural waters at the Marcell Experimental Forest

Here, we present a new data set from the Marcell Experimental Forest (MEF) that compiles water isotope measurements from multiple research catchments, some of which have been studied since the 1960s. The MEF is located in northern Minnesota, USA, and is home to heavily studied and monitored forests, streams, bogs, and fens. Peat-forming systems (bogs and fens) are an important component of the MEF landscape and have a profound impact on the water cycle in these catchments. Within the last decade, analysis of stable isotopes of water (expressed as δD and δ 18 O) has been implemented to characterize the different components of the water budget, and to allow researchers to look at catchment and peatland-specific hydrologic effects in the watershed. This δD and δ 18 O data set of natural waters from MEF catchments is primarily composed of measurements from three peatlands (S1, S2, S6) during an 11-year period. More recently collection and analysis were expanded to also include samples from the Spruce and Peatland Responses Under Changing Environments (SPRUCE) project in the S1 bog, peatlands S3, S4, S5, as well as nearby lakes. We establish a local meteoric water line by analyzing the isotopic composition of precipitation, which fills a void in regional meteoric water lines for Minnesota. Furthermore, we establish baseline isotopic composition for bog outlet streams, bog porewater, aquifer groundwater, overland flow, subsurface stormflow, and snowpack, as well as runoff from the SPRUCE experimental chambers. These data are publicly available and will be expanded upon in the future.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

In Situ Monitoring the Nucleation and Growth of Nanoscale CaCO 3 at the Oil–Water Interface

Interfaces can actively control the nucleation kinetics, orientations, and polymorphs of calcium carbonate (CaCO 3 ). Prior studies have revealed that CaCO 3 formation can be affected by the interplay between chemical functional moieties on solid–liquid or air–liquid interfaces as well as CaCO 3 ’s precursors and facets. Yet little is known about the roles of a liquid–liquid interface, specifically an oil–liquid interface, in directing CaCO 3 mineralization which are common in natural and engineered systems. Here, in this study, by using in situ X-ray scattering techniques to locate a meniscus formed between water and a representative oil, isooctane, we successfully monitored CaCO 3 formation at the pliable isooctane–water interface and systematically investigated the pivotal roles of the interface in the formation of CaCO 3 (i.e., particle size, its spatial distribution with respect to the interface, and its mineral phase). Different from bulk solution, ∼5 nm CaCO 3 nanoparticles form at the isooctane–water interface. They stably exist for a long time (36 h), which can result from interface-stabilized dehydrated prenucleation clusters of CaCO 3 . There is a clear tendency for enhanced amounts and faster crystallization of CaCO 3 at locations closer to isooctane, which is attributed to a higher pH and an easier dehydration environment created by the interface and oil. Our study provides insights into CaCO 3 nucleation at an oil–water interface, which can deepen our understanding of pliable interfaces interacting with CaCO 3 and benefit mineral scaling control during energy-related subsurface operation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

HFTS-1 Natural Joint Data and Engineering Summary

This report provides an analysis of engineering and geologic data collected from the Hydraulic Fracturing Test Site #1 (HTSF-1) project in the southern Midland Basin, Reagan County, Texas. The site is being studied as part of the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative at the National Energy Technology Laboratory (NETL). The data collected is intended to provide a basis of understanding of the site and to construct simulation models. This report provides a summary analysis of the engineering and geologic data collected from the project to provide a basis of understanding of the site and to construct simulation models. The report examines various possible correlations in engineering properties based on natural fracture data from the program, which consisted of 11 horizontal wells, one vertical well, and one slant well. The horizontal wells are in two horizons: 1) Upper and 2) Middle Wolfcamp formations. Program data include: fracture frequency and fracture orientation data from four core runs in a slant well; fracture spacing and orientation from a vertical pilot well; laboratory triaxial testing and mineralogical determinations; and porosity results from magnetic resonance analyses from various wells, together with observations based on the data collection. In addition, available references were reviewed on the site for additional insights. As the focus of the report is on the natural system, hydraulic fracture data from the site were not examined in detail in this report. Data variability is the chief observation in examination of the database. Fracture frequency in the slant well can range from sections with values as high as five fractures per ft to sections up to 100+ ft in length with no natural observed fractures. Fracture spacing across is typically less than 10 ft, but can range up to hundreds of feet. Apparent fracturing shows the trends in two predominate orientations, E-SW and WNW-ESE, but minor variations exist. The rock units vary across the site from siliceous mudstones to calcareous mudstones, showing a general layering with depth. The laboratory properties such as strength and modulus show no apparent trend with depth, but appear to correlate with rock mineralogy with high strength and modulus values where calcium content is high. In addition, an attempt to examine variability and mineralogy on a larger scale was made using a color-coded system based on gamma ray measurements. A staged colored approach was adopted, presuming that lower gamma ray values indicate higher value of calcium content (blue scale) and that higher gamma ray values indicate higher clay mineral content (orange scale). As provided in report appendices, the system correlated well with visual examination of the slant core and the petrofabric analyses of the vertical pilot well. The results showed large variability in mineral content along the horizontal plane across the site. The change in mineral content was also rapid, on a scale less than that of the average hydraulic fracture stage length of about 180 ft.

42 ENGINEERING↗

DOE OSTI2023

Application of Electromagnetic Methods for Reservoir Monitoring with Emphasis on Carbon Capture, Utilization, and Storage

The Controlled-Source ElectroMagnetic (CSEM) method provides crucial information about reservoir fluids and their spatial distribution. Carbon dioxide (CO2) storage, enhanced oil recovery (EOR), geothermal exploration, and lithium exploration are ideal applications for the CSEM method. The versatility of CSEM permits its customization to specific reservoir objectives by selecting the appropriate components of a multi-component system. To effectively tailor the CSEM approach, it is essential to determine whether the primary target reservoir is resistive or conductive. This task is relatively straightforward in CO2 monitoring, where the injected fluid is resistive. However, for scenarios involving brine-saturated (water-wet) or oil-wet (carbon capture, utilization, and storage—CCUS) reservoirs, consideration must also be given to conductive reservoir components. The optimization of data acquisition before the survey involves analyzing target parameters and the sensitivity of multi-component CSEM. This optimization process typically includes on-site noise measurements and 3D anisotropic modeling. Based on our experience, subsequent surveys tend to proceed smoothly, yielding robust measurements that align with scientific objectives. Other critical aspects to be considered are using magnetotelluric (MT) measurements to define the overall background resistivities and integrating real-time quality assurance during data acquisition with 3D modeling. This integration allows the fine tuning of acquisition parameters such as acquisition time and necessary repeats. As a result, data can be examined in real-time to assess subsurface information content while the acquisition is ongoing. Consequently, high-quality data sets are usually obtained for subsequent processing and initial interpretation with minimal user intervention. The implementation of sensitivity analysis during the inversion process plays a pivotal role in ensuring that the acquired data accurately respond to the target reservoirs’ expected depth range. To elucidate these concepts, we present an illustrative example from a CO2 storage site in North Dakota, USA, wherein the long-offset transient electromagnetic method (LOTEM), a variation of the CSEM method, and the MT method were utilized. This example showcases how surface measurements attain appropriately upscaled log-scale sensitivity. Furthermore, the sensitivity of the CSEM and MT methods was examined in other case histories, where the target reservoirs exhibited conductive properties, such as those encountered in enhanced oil recovery (EOR), geothermal, and lithium exploration applications. The same equipment specifications were utilized for CSEM and MT surveys across all case studies.

Barajas-Olalde, César↗

DOE OSTI2021

Surface Meteorology at Teller Mile 47 Watershed, Seward Peninsula, Alaska, Ongoing from 2018

Meteorological data are currently being collected at one location at the Teller Mile 47 (TL_MM47) Research Basin Site, Seward Peninsula (N64 58' 36.918", W166 12' 32.67", 67 meters above sea level). The site was installed and initial measurements started in September 2018 and it has operated continuously since then. The meteorological station is co-located with a continuous snow depth sensor and two soil pits for subsurface temperature and moisture measurements. These data are being collected to better understand the energy dynamics above the active layer and permafrost. They complement in-situ snow and soil measurements also at this location. The data could also be used as supporting measurements for other research and modeling activities. There are 35 comma separated value format (*.csv) files provided, where each file contains the full data for an individual parameter (e.g. air temperature at 1.5 meters above the ground surface (teller_m47_air_temperature_150cm_ags_Avg.csv) or soil temperature 20 centimeters below ground surface (teller_m47_dry_soil_pit_temperature_20cmbgs_Avg.csv)) plus the time in Universal Coordinated Time (UTC) and Alaska Standard Time (UTC time minus nine hours). The site was installed and initial measurements started in September 2018. It has been operated continuously since. Primary data gaps are due to battery failure or sensor failure. These data are being collected to better understand the surface energy dynamics above the active layer and permafrost. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 10-year research effort (2012-2022) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy?s Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy?s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Computed Tomography Scanning and Geophysical Measurements of UW Enterprises LP 1-250512-129 Well in Southwestern Indiana

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia, were used to characterize Illinois Basin core of the Upper Devonian-Early Mississippian New Albany Shale Formation from Posey County, Indiana. The primary impetus of this work is a collaboration between Indiana Geological and Water Survey at Indiana University Bloomington, NETL, and Woolsey Operating Company LLC to characterize and make publicly available core information from the New Albany Shale of the Illinois Basin. Core characterization of this unconventional oil/gas well will aid in understanding the lithology changes and the fracture complexity of the New Albany Shale. There is a potential for this formation to be developed in the future for critical mineral and rare earth element (CM/REE) extraction. The resultant datasets are presented in this report and can be accessed from NETL's Energy Data eXchange (EDX) online system using the following link: https://edx.netl.doe.gov/dataset/uw-enterprises. All equipment and techniques used were non-destructive, enabling future examinations and analyses to be performed on this core. None of the equipment used was suitable for direct visualization of the pore space in the fine-grained structures studied; however, fractures, discontinuities, and millimeter-scale features were readily detectable with the methods tested. Imaging with the NETL medical CT scanner was performed on the entire core. Targeted higher resolution CT scanning of select sections was performed with NETL’s industrial and micro-CT scanner. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), P-wave, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for more detailed analysis. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core; the resulting macro and micro descriptions are relevant to many subsurface energy related examinations traditionally performed at NETL.

58 GEOSCIENCES↗

DOE OSTI2025

Structural Evolution of the Hogback Monocline and Its Tectonic Significance in the San Juan Basin

The San Juan Basin is recognized as a Laramide foreland basin. It is located within the Colorado Plateau, a broad tectonic province characterized by a thick sedimentary sequence that was segmented into smaller sub basins during the Late Cretaceous to Paleogene Laramide orogeny. The Hogback Monocline lies along the northwestern margin of the San Juan Basin and is considered a Laramide-age structure formed in response to compressional stress. In this study, we interpret surface and subsurface datasets to construct a structural geological model and evaluate its tectonic significance. Through seismic data, we identify key fault and fold geometries at depth. The seismic dataset used in this study was reprocessed in depth and constrained with well log velocity data to enhance seismic imaging quality. Additionally, we performed well log correlations to identify formation tops and assess variations in basin infill and thickness geometry. A series of structural cross-sections, constructed using seismic data and a high density of boreholes, are presented to evaluate geometric variations along the structure and its evolution during basin development. Furthermore, kinematic restoration and forward modeling analyses were conducted to validate our structural interpretation. This work suggests that the Hogback Monocline formed through fault-propagation folding and flexural slip affecting the pre-Laramide sedimentary sequence under compressional stresses associated with the Laramide orogeny. This structure is interpreted as a high-angle reverse fault that influenced the geometry of the late basin infill. Additionally, monocline bending along the structure may have been controlled by fault relay systems and, in some cases, influenced by strike-slip faulting.

Reyes, Martin [New Mexico Bureau o fGeology and Mi↗

DOE OSTI2026

geoPFA: A Python-Based Open-Source Software for 3D Geothermal PFA

This work presents a novel Python-based framework, geoPFA, for conducting 3D play fairway analysis (PFA) tailored to superhot geothermal systems. The workflow has been applied to the Nesjavellir field in Iceland, a candidate site for the third Iceland Deep Drilling Project's superhot production scenarios. This application demonstrates the value of modular, transparent, and extensible workflows for integrating geological, geophysical, and simulation-derived datasets in high-enthalpy environments. Preliminary results indicate favorable zones consistent with known hydrothermal activity. The geoPFA library will soon be publicly available, offering a scalable and reproducible approach to geothermal exploration across varied geological contexts.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

Massively parallel modeling and inversion of electrical resistivity tomography data using PFLOTRAN

Abstract. Electrical resistivity tomography (ERT) is a broadly accepted geophysical method for subsurface investigations. Interpretation of field ERT data usually requires the application of computationally intensive forward modeling and inversion algorithms. For large-scale ERT data, the efficiency of these algorithms depends on the robustness, accuracy, and scalability on high-performance computing resources. In this regard, we present a robust and highly scalable implementation of forward modeling and inversion algorithms for ERT data. The implementation is publicly available and developed within the framework of PFLOTRAN, an open-source, state-of-the-art massively parallel subsurface flow and transport simulation code. The forward modeling is based on a finite-volume discretization of the governing differential equations, and the inversion uses a Gauss–Newton optimization scheme. To evaluate the accuracy of the forward modeling, two examples are first presented by considering layered (1D) and 3D earth conductivity models. The computed numerical results show good agreement with the analytical solutions for the layered earth model and results from a well-established code for the 3D model. Inversion of ERT data, simulated for a 3D model, is then performed to demonstrate the inversion capability by recovering the conductivity of the model. To demonstrate the parallel performance of PFLOTRAN's ERT process model and inversion capabilities, large-scale scalability tests are performed by using up to 131 072 processes on a leadership class supercomputer. These tests are performed for the two most computationally intensive steps of the ERT inversion: forward modeling and Jacobian computation. For the forward modeling, we consider models with up to 122 ×106 degrees of freedom (DOFs) in the resulting system of linear equations and demonstrate that the code exhibits almost linear scalability on up to 10 000 DOFs per process. On the other hand, the code shows superlinear scalability for the Jacobian computation, mainly because all computations are fairly evenly distributed over each process with no parallel communication.

58 GEOSCIENCES↗

NASA NTRS1996

Design and Laboratory Validation of a Capacitive Sensor for Measuring the Recession of Thin-Layered Ablator

Flight vehicles are typically instrumented with subsurface thermocouples to estimate heat transfer at the surface using inverse analysis procedures. If the vehicle has an ablating heat shield, however, temperature time histories from subsurface thermocouples no longer provide enough information to estimate heat flux at the surface. In this situation, the geometry changes and thermal energy leaves the surface in the form of ablation products. The ablation rate is required to estimate heat transfer to the surface. A new concept for a capacitive sensor has been developed to measure ablator depth using the ablator's dielectric effect on a capacitor's fringe region. Relying on the capacitor's fringe region enables the gage to be flush mounted in the vehicle's permanent structure and not intrude into the ablative heat shield applied over the gage. This sensor's design allows nonintrusive measurement of the thickness of dielectric materials, in particular, the recession rates of low-temperature ablators applied in thin (0.020 to 0.060 in. (0.05 to 0.15 mm)) layers. Twenty capacitive gages with 13 different sensing element geometries were designed, fabricated, and tested. A two-dimensional finite-element analysis was performed on several candidate geometries. Calibration procedures using ablator-simulating shims are described. A one-to-one correspondence between system output and dielectric material thickness was observed out to a thickness of 0.055 in. (1.4 mm) for a material with a permittivity about three times that of air or vacuum. A novel method of monitoring the change in sensor capacitance was developed. This technical memorandum suggests further improvements in gage design and fabrication techniques.

Gregory K Noffz↗

NASA NTRS2016

Sylph - A SmallSat Probe Concept Engineered to Answer Europa's Big Question

Europa, the Galilean satellite second from Jupiter, contains a vast, subsurface ocean of liquid water. Recent observations indicate possible plume activity. If such a plume expels ocean water into space as at Enceladus, a spacecraft could directly sample the ocean by analyzing the plume’s water vapor, ice, and grains. Due to Europa’s strong gravity, such sampling would have to be done within 5 km of the surface to sample ice grains larger than 5 µm, expected to be frozen ocean spray and thus to contain non-volatile species critical to a biosignature-detection mission. By contrast, the planned Europa Multiple Flyby Mission’s closest planned flyby altitude is 25 km. Sylph is a concept for a SmallSat free-flyer probe that, deployed from the planned Europa Mission, would directly sample the large grains by executing a single ~2-km altitude plume pass. The 40-kg probe would be deployed by the Europa mission just before it executes a plume fly-through. Within the probe’s 16-hour lifespan, it would autonomously navigate to perform a parallel, simultaneous pass at the lower altitude. The Sylph flight system design concept combines SmallSat technologies with robust traditional components and advanced manufacturing technologies. Its payload would comprise Mini-SUDA (SUrface Dust Mass Analyzer), a dual-channel, miniature impact ionization mass spectrometer. Sylph represents a novel type of SmallSat concept: purpose-built configuration, optimized for the harsh environment at Europa and for planetary-protection requirements, and hybridized from both mainstream and SmallSat components.

Imken, Travis↗