Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface Modeling”

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

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

936 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Geospatial modeling of near subsurface temperatures of the contiguous United States for assessment of materials degradation

Abstract Understanding subsurface temperature variations is crucial for assessing material degradation in underground structures. This study maps subsurface temperatures across the contiguous United States for depths from 50 to 3500 m, comparing linear interpolation, gradient boosting (LightGBM), neural networks, and a novel hybrid approach combining linear interpolation with LightGBM. Results reveal heterogeneous temperature patterns both horizontally and vertically. The hybrid model performed best achieving a root mean square error of 2.61 °C at shallow depths (50–350 m). Model performance generally decreased with depth, highlighting challenges in deep temperature prediction. State-level analyses emphasized the importance of considering local geological factors. This study provides valuable insights for designing efficient underground facilities and infrastructure, underscoring the need for depth-specific and region-specific modeling approaches in subsurface temperature assessment.

Science & Technology - Other Topics↗
DOE OSTI2025

Carbon Utilization and Storage Partnership of the Western United States

This technical report documents research conducted under DOE Award No. DE-FE0031837 focused on evaluating the feasibility of carbon capture, utilization, and storage (CCUS) systems in the central and western United States. The project integrated geologic characterization, reservoir simulation, infrastructure modeling, and economic analysis to assess CO₂ storage potential near industrial sources and develop strategies for transport and sequestration. The work included subsurface modeling, risk assessment, monitoring and verification (MRV) planning, and evaluation of regulatory pathways such as EPA Underground Injection Control (UIC) Class VI permitting and IRS 45Q tax credit eligibility. Results demonstrate the viability of multiple storage approaches, including saline formations, enhanced coalbed methane recovery, and basalt mineralization, supported by data-driven workflows and regional analyses. The project also produced permitting templates, technology transfer activities, and stakeholder engagement efforts to support deployment readiness. These findings contribute to the development of scalable, economically viable CCUS systems and provide a repeatable framework for future carbon management projects.

20 FOSSIL-FUELED POWER PLANTS↗
NASA NTRS2026

Synchronized Eruptions on Io: Possible Evidence of Interconnected Subsurface Magma Reservoirs

On 27 December 2024, Juno's JIRAM (Jovian InfraRed Auroral Mapper) instrument observed an unprecedented volcanic event in Io's southern hemisphere, covering a vast region of ∼65,000 km2, near 73°S, 140°E. Within the imaged region, only one hot spot was previously known (Pfd454). This feature was earlier estimated to cover an area of 300 km2 with a total power output of 34 GW. JIRAM results show that the region produces a power output of 140–260 TW, over 1,000 times higher than earlier estimates and likely exceeding the brightest eruption ever recorded on Io, that of Surt in 2001 (∼80 TW). Three adjacent hot spots also exhibited dramatic power increases: P139, PV18, and an unnamed feature south of the main one that surged to ∼1 TW, placing all of them among the top 10 most powerful hot spots observed on Io. A temperature analysis of the features supports the simultaneous onset of these brightenings and suggests a single eruptive event propagating beneath the surface across hundreds of kilometers; this is the first time this has been observed on Io. This in turn would imply a connection among the hotspots' magma reservoirs, while other nearby hotspots that have been known to be active in the recent past, such as Kurdalagon Patera, appear unaffected. The simultaneity supports models of massive, interconnected magma reservoirs. The topology of these regional magma systems may resemble that of a large-scale sponge, in which the massive reservoirs are the pores, interconnected through a largely solid outer shell.

A Mura↗
DOE OSTI2026

Decadal-Scale Hydrothermal Alteration of Basalt: Implications for Long-Term Reactivity in Subsurface Technologies

Global demands to secure critical mineral supply chains and managing water resources require innovative resource management strategies. Subsurface basalt reservoirs offer a dual-purpose application, where CO2 enhanced critical metal recovery can be coupled with water remediation strategies for the extraction of economically important and environmentally relevant elements. This study utilizes static batch experiments conducted on Columbia River Basalt Group samples to evaluate long-term fluid–rock interactions relevant to coupled subsurface applications. Samples were reacted under elevated temperature and pressure conditions, for 296 and 4,061 days, providing a unique dataset to assess geochemical evolution, element partitioning, and aqueous chemistry over decadal timescales. The alteration was dominated by silicate, aluminosilicate, and oxide phases. Progressive hydrothermal alteration increased the cation exchange capacity through the formation of secondary minerals (e.g. clays and zeolites), influencing solute retention, ion exchange behavior, and aqueous chemistry. Grain-scale and bulk analyses revealed selective cation release and retention patterns, establishing a mechanistic baseline for predicting element mobility. These geochemical constraints inform the prospective use of CO2-enhanced metal recovery for critical resources (e.g. Co, Ni, Mn, Mg, Ti, Fe) and associated water-quality impacts, providing a foundation for the design and management of engineered fluid–rock systems in basalt reservoirs.

Gibson, Marie N.↗
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 NTRS2003

Global Characterization of Polygonally Fractured Terrain on Venus and Implications for a Climate Change Origin

On Earth, polygonal fractures with spacings of 1-100 cm form in lava flows that cool slowly under isotropic stress conditions. On Venus, polygonal fractures are observed with typical fracture spacings of 1-2 km, covering regions 10s to 100s of km across. These features have been proposed to form via cooling of lava flows or above subsurface intrusions, or due to surface cooling caused by climate change. A watershed algorithm was modified to search the entire Magellan image database autonomously for terrains with polygonal fractures. Here we characterize the 204 polygon terrains found using this method (including over 170 newly identified regions) with respect to size, morphology, stratigraphy, and geologic setting. These results are used to examine different models of origin.

S. E. Smrekar↗
DOE OSTI2026

Carbon Tetrachloride Degradation Results for 200-ZP-1 Operable Unit

Carbon tetrachloride (CT) contamination in the 200-ZP-1 Operable Unit (OU) at the Hanford Site originated from large-volume discharges to the subsurface during plutonium production operations between 1955 and 1973. Contamination migrated through more than 70 meters of unsaturated sediment to reach the underlying unconfined aquifer, where it persists as a large and complex groundwater plume. The 200-ZP-1 OU Record of Decision (ROD) requires that groundwater CT concentrations be reduced to 3.4 µg/L within 125 years. Current groundwater modeling projections estimate that the existing pump-and-treat, even when combined with monitored natural attenuation (specifically hydrolysis), will not achieve this target within the designated timeframe. A fundamental contributor to this shortfall is the extremely slow rate of CT hydrolysis under Hanford aquifer conditions, which has been estimated to have a half-life of 630 years. If faster-acting biotic and abiotic degradation processes are operating within the aquifer, their contribution to CT mass reduction could have a meaningful impact. However, site-specific measurements of these processes and their rates have not previously been performed. This report documents the results of a two-phase laboratory investigation designed to characterize and quantify the capacity of site-specific 200-ZP-1 OU sediments and groundwater to support natural attenuation of CT through biotic and abiotic pathways. In this context, degradation capacity is defined as the intrinsic potential of the subsurface matrix to transform CT under optimized, controlled conditions. System capacity is evaluated in two ways: (1) as rate-limited capacity, which establishes the maximum kinetic velocity of CT transformation and is measured using half-lives and first order rate constants; and (2) as mass limited capacity, which defines the total contaminant mass the batch experimental system can degrade before reactants are exhausted, representing the maximum amount of contaminant the microbial community and reactive mineral phases can transform under the experimental conditions.

abiotic degradation↗
DOE OSTI2024

Mechanochemically responsive polymer enables shockwave visualization

Abstract Understanding the physical and chemical response of materials to impulsive deformation is crucial for applications ranging from soft robotic locomotion to space exploration to seismology. However, investigating material properties at extreme strain rates remains challenging due to temporal and spatial resolution limitations. Combining high-strain-rate testing with mechanochemistry encodes the molecular-level deformation within the material itself, thus enabling the direct quantification of the material response. Here, we demonstrate a mechanophore-functionalized block copolymer that self-reports energy dissipation mechanisms, such as bond rupture and acoustic wave dissipation, in response to high-strain-rate impacts. A microprojectile accelerated towards the polymer permanently deforms the material at a shallow depth. At intersonic velocities, the polymer reports significant subsurface energy absorption due to shockwave attenuation, a mechanism traditionally considered negligible compared to plasticity and not well explored in polymers. The acoustic wave velocity of the material is directly recovered from the mechanochemically-activated subsurface volume recorded in the material, which is validated by simulations, theory, and acoustic measurements. This integration of mechanochemistry with microballistic testing enables characterization of high-strain-rate mechanical properties and elucidates important insights applicable to nanomaterials, particle-reinforced composites, and biocompatible polymers.

Science & Technology - Other Topics↗
NASA NTRS2026

The Role of Disruptive Impacts on Ocean Generation and Longevity in Icy Moons

Among icy moons of the outer Solar System, subsurface water oceans and large collisions both seem to be common, but the impact of the latter on the presence and persistence of the former is unclear and has rarely been investigated. Here we interface a smoothed-particle hydrodynamics model to simulate collisions with a thermal-structural evolution model to simulate the evolution of moons pre-collision, post-collision and without a collision. Overall, even such large-scale collisions affect only the ocean thickness or longevity, and the presence or absence of an ocean is affected for only a part of a moon’s history. In reaccreted moons, the ocean survives the impact and becomes much thicker inside larger moons with radius near 1,000 km, whereas an ocean that would otherwise arise inside moons with radius near 500 km is absent because the collision promotes ice–rock differentiation. Our simulations have not yielded an ocean developed post-impact—whether directly via collisional or reaccretional heating or indirectly through tidal heating due to collision-induced orbital changes—in a moon that would otherwise have remained frozen. The ocean-enhancing effect is pronounced only for late disruptive impacts onto large, 1,000-km-class targets, which are unlikely in recent solar system history.

Marc Neveu↗
DOE OSTI2026

Upper ocean mixing, surface heat fluxes, and heat content variability in the upper 150 m during Hurricane Laura (2020)

Observations (from thermistor sensors and a hurricane glider) and model-based temperature were used to answer questions including the following: 1) What physical and environmental conditions explain the rapid intensification of Hurricane Laura? 2) How did the pre-existing warm mixed layer at Stone mooring (StM) influence the degree of cooling in the mixed layer? 3) How did vertical mixing, surface heat fluxes, and advective processes collectively shape the mixed layer heat evolution at StM? Observations measured by thermistor sensors showed that Hurricane Laura induced a cooling of -1.2oC in the mixed layer when Hurricane Laura’s wind speed increased to 47 m/s on 26 August. This contrasts with the -1.04oC mixed layer temperature change estimated from model-based temperature. At StM, the presence of a 31oC warm mixed layer and elevated heat content (60-80kJ/cm2) effectively preconditioned the upper ocean ahead of Hurricane Laura’s passage. After Hurricane Laura passed, turbulent mixing associated with the storm transported the heat anomalies downward into subsurface layers. In addition, this study shows that a surface heat flux of 5.04 kJ/cm2 supported the intensity of Hurricane Laura as it traversed over the StM. Using a mixed layer heat balance model, this study shows that entrainment flux and surface flux mostly contributed to the observed change in the mixed layer temperature.

Tsai, Senam↗
NASA NTRS1999

DS-2 Mars Microprobe Battery

In January of 1999 the NM DS-2 Mars microprobe will be launched to impact on Mars in December. The technical objectives of the missions are to demonstrate: key technologies, a passive atmospheric entry, highly integrated microelectronics which can withstand both low temperatures and high decelerations, and the capability to conduct in-situ, surface and subsurface science data acquisition. The scientific objectives are to determine if ice is present below the Martian surface, measure the local atmospheric pressure, characterize the thermal properties of the martian subsurface soil, and to estimate the vertical temperature gradient of the Martian soil. The battery requirements are 2-4 cell batteries, with voltage of 6-14 volts, capacity of 550 mAh at 80C, and 2Ah at 25C, shelf life of 2.5 years, an operating temperature of 60C and below, and the ability to withstand shock impact of 80,000 g's. The technical challenges and the approach is reviewed. The Li-SOCL2 system is reviewed, and graphs showing the current and voltage is displayed, along with the voltage over discharge time. The problems encountered during the testing were: (1) impact sensitivity, (2) cracking of the seals, and (3) delay in voltage. A new design resulted in no problems in the impact testing phase. The corrective actions for the seal problems involved: (1) pre weld fill tube, (2) an improved heat sink during case to cover weld and (3) change the seal dimensions to reduce stress. To correct the voltage delay problem the solutions involved: (1) drying the electrodes to reduce contamination by water, (2) assemblage of the cells within a week of electrode manufacture, (3) ensure electrolyte purity, and (4) provide second depassivation pulse after landing. The conclusions on further testing were that the battery can: (1) withstand anticipated shock of up to 80,000 g, (2) meet the discharge profile post shock at Mars temperatures, (3) meet the required self discharge rate and (4) meet environmental requirements.

Frank, H.↗
DOE OSTI2026

Enhanced Recovery of Critical Minerals and Geological Hydrogen in The Mine of the Future

invited perspective article: The rapid demand for critical minerals (CMs) and hydrogen (H2) necessitates innovative solutions beyond conventional mining. This study explores enhanced mineral recovery (EMR) and geological hydrogen production (GeoH2) by leveraging the Earth’s subsurface as a reactive platform. Using ultramafic rocks and engineered fluids, the approach simultaneously mobilizes critical resources and produces H2 through natural processes like serpentinization. This paradigm offers a transformative pathway to secure essential materials and diversify energy resources, setting a foundation for the Terrestrial Mine of the Future

Yan, Keju↗
DOE OSTI2026

Value of Geothermal Energy Storage for Supply-Side and Demand-Side Applications

This report presents the results of a study examining the value potential for geothermal energy storage (GES), a long-duration energy storage resource that stores thermal and/or geomechanical energy in the subsurface. GES could benefit the overall U.S. power system by temporally shifting electricity generation (supply-side) or meeting building heating and cooling load (demand-side). This report analyzes supply-side and demand-side opportunities independently because of differences in applications and models. Currently there is significant uncertainty about the development costs for GES, with only a limited number of demonstration plants for electric energy storage and building heating and cooling storage developments. In this report, we estimate the value of supply-side and demand-side GES to the bulk power system in the contiguous United States. Because of the significant uncertainty about GES development costs, this analysis does not consider GES deployment costs but instead focuses on the value of GES to the U.S. electricity system. The estimated values of GES provide reference points for economically competitive commercial cost targets. Supply-side GES is modeled as part of an enhanced geothermal system (EGS) generation plant in NREL's Regional Energy Deployment System (ReEDS) capacity expansion model (Ho et al. 2021). In contrast to conventional geothermal plants, which generate constant power, EGS plants have unique features that may allow for in-reservoir energy storage for flexible generation. Demand-side GES for heating and cooling, including seasonal hot and cold storage and short-duration heat pump storage, is incorporated into a price-taker model using Cambium electricity marginal cost projections. To establish an upper bound for the value of GES, analysis focused on favorable scenarios for storage with high generation from zero marginal cost, variable renewable energy resources. High penetrations of variable renewable energy generation can increase hourly electricity price variability, which increases the value of temporal energy arbitrage for storage technologies like GES.

15 GEOTHERMAL ENERGY↗
NASA NTRSDate not supplied

MEAD In-situ Sampling and Testing at Haughton Crater

The Mars Exploration through Analog-site Drilling (MEAD) project is designed to demonstrate the feasibility and scientific value of robotic drilling missions to Mars in the search for evidence of past or extant life. MEAD conducts high-fidelity field simulations integrating drilling, autonomous operations, and in-situ life-detection and mineralogical instruments. The project tests the Signs of Life Detector (SOLID), proposed for the 2019 Discovery Mars Icebreaker mission, and ARIA (Astronaut Raman Instrument for ISRU and Astrobiology), which provides mineral identification and detects trace organics and volatiles, including indicators of organic thermal maturity. Because current and planned Mars drilling missions operate largely “blind,” MEAD also evaluates the use of drill-induced vibrations recorded by deployed geophones to enable local subsurface mapping for improved targeting of drilling and sampling. In August 2025, MEAD completed its first-year field deployment to Haughton Crater in the Canadian High Arctic, an established Mars impact-crater terrestrial analog site, bringing a Honeybee Robotics TRIDENT drill with applied automation, alongside ARIA, SOLID, and FleetSpace geophones.

Haughton Crater↗
DOE OSTI2026

Role of Fluid and Temperature in Fracture Mechanics and Coupled THMC Processes for Enhanced Geothermal Systems

The ability to sustain high energy extraction efficiency from Enhanced Geothermal Systems (EGS) is affected by the ability to measure, characterize, and predict the effect of different perturbations that arise from fluid injection rates, fluid temperature, and shut-in conditions on existing and stimulated fracture systems throughout a subsurface reservoir’s lifecycle. The difficulty arises from limited knowledge of the interaction of fluid and temperature driven fractures with frictional interfaces that govern local deformation and frictional behavior in the surrounding rock. Thermo-poro-mechanical coupled processes tend to dominate this interaction and strongly affect the permeability and local stress distributions that impact efficiency and prevent optimal stimulation conditions over multi-decadal time frames.

Pyrak-Nolte, Laura [Purdue Univ., West Lafayette, ↗
DOE OSTI2023

Radiolytically reworked Archean organic matter in a habitable deep ancient high-temperature brine

Abstract Investigations of abiotic and biotic contributions to dissolved organic carbon (DOC) are required to constrain microbial habitability in continental subsurface fluids. Here we investigate a large (101–283 mg C/L) DOC pool in an ancient (>1Ga), high temperature (45–55 °C), low biomass (10 2 −10 4 cells/mL), and deep (3.2 km) brine from an uranium-enriched South African gold mine. Excitation-emission matrices (EEMs), negative electrospray ionization (–ESI) 21 tesla Fourier-transform ion cyclotron resonance mass spectrometry (FT-ICR MS), and amino acid analyses suggest the brine DOC is primarily radiolytically oxidized kerogen-rich shales or reefs, methane and ethane, with trace amounts of C 3 –C 6 hydrocarbons and organic sulfides. δ 2 H and δ 13 C of C 1 –C 3 hydrocarbons are consistent with abiotic origins. These findings suggest water-rock processes control redox and C cycling, helping support a meagre, slow biosphere over geologic time. A radiolytic-driven, habitable brine may signal similar settings are good targets in the search for life beyond Earth.

Science & Technology - Other Topics↗