Engineering Papers⌕ Search

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.

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

At least 307 records · Page 17

The evolving role of geothermal energy for decarbonizing the United States

Geothermal energy is often referred to as a niche technology that is too localized, too small or too expensive to make much of a difference in how renewable energy will be supplied in a fully decarbonized future. As a result, geothermal energy has been undervalued in terms of what it could provide to complement, rather than compete with, electricity generation from wind, solar photovoltaic, concentrating solar power and other renewables. Geothermal energy systems are fully dispatchable and can provide baseload or load-following electric power or heat suitable for a wide range of applications including supplying district heating for communities and cities, and heating and cooling of individual buildings. The focus of our study is on the potential of utilizing geothermal energy for providing heat for buildings and industry at lower temperatures, a substitute for the combustion of fossil fuels. Because heating represents about 20% (20 EJ per year) of the annual primary energy consumption in the U.S.—with most of it coming from burning natural gas, oil and/or propane in furnaces—deploying geothermal heating on a national scale could have a significant impact on lowering carbon emissions. In heating-dominated states in the U.S. Northern Tier, heating often is among the largest contributors to the state's carbon footprint. This review begins with a discussion of the motivation and rationale behind considering geothermal as a key low-carbon heating option for the U.S. The study summarizes the U.S. geothermal resource and describes the applications and main engineering components of using geothermal energy for heating and cooling, electric power generation, and co-generation using district heating, geothermal heat pumps, and power conversion with steam flashing and organic Rankine plants. Environmental benefits and impacts are described. An extensive discussion of geologic and thermal-hydraulic aspects of the subsurface is included in the review because of their critical role in determining reservoir designs at specific sites to ensure sufficient productivity that is both safe and economically viable. Models for estimating levelized costs of district heating are used to show how costs are affected by reservoir performance, infrastructure capital costs, and financial parameters. Here, the review concludes with an assessment of technical subsurface issues associated with reservoir performance and the economic requirements for providing geothermal heating in district heating systems at a sufficient scale to have an impact on decarbonizing the U.S.

15 GEOTHERMAL ENERGY↗

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗

A Review of Sandia Energy Storage Research Capabilities and Opportunities (2020 to 2030)

Large-scale integration of energy storage on the electric grid will be essential to enabling greater penetration of intermittent renewable energy sources, modernizing the grid for increased flexibility security, reliability, and resilience, and enabling cleaner forms of transportation. The purpose of this report is to summarize Sandia's research and capabilities in energy storage and to provide a preliminary roadmap for future efforts in this area that can address the ongoing program needs of DOE and the nation. Mission and vision statements are first presented followed by an overview of the organizational structure at Sandia that provides support and activities in energy storage. Then, a summary of Sandia's energy storage capabilities is presented by technology, including battery storage and materials, power conversion and electronics, subsurface-based energy storage, thermal/thermochemical energy storage, hydrogen storage, data analytics/systems optimization/controls, safety of energy storage systems, and testing/demonstrations/model validation. A summary of identified gaps and needs is also presented for each technology and capability.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Conceptual Model Testing Related to SDU 6 Drainwell Observations

From its inception in the early 1950s through the end of the Cold War in the early 1990s, the Savannah River Site (SRS) produced nuclear materials for national defense in five reactors. Additionally, irradiated reactor fuel and target tubes were dissolved in nitric acid to recover plutonium and uranium using the PUREX (Plutonium Uranium Reduction EXtraction) process. Liquid waste from these chemical separations processes was then stored onsite in 51 underground tanks. Eight waste storage tanks have been operationally closed (i.e. cleaned and grouted) and the remaining tanks hold a mixture of liquids, insoluble solids, and precipitated salts (SRMC-LWP-2022-00001), the latter generated by evaporating water from the liquid waste. Waste is currently being retrieved from tanks and separated into 1) high-radioactivity, low-volume, and 2) low-radioactivity, high-volume components, principally through the Salt Waste Processing Facility (SWPF) (SRMC-LWP-2023-00001). The former waste stream is vitrified in the Defense Waste Processing Facility (DWPF), stored onsite, and destined for offsite disposal in a deep geologic repository. The latter stream is mixed with dry cementitious materials in the Saltstone Production Facility (SPF) and the wet slurry placed in onsite Saltstone Disposal Units (SDUs) within the Saltstone Disposal Facility (SDF), where it hardens into a cement waste form termed saltstone. A low-infiltration surface cover system will be placed over the SDF at closure, where SDUs will then be in the subsurface post-closure (SRR-CWDA-2019-00001).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Fracture Intersections under Stress: Laboratory Data and Code [Data set]

The connectivity of natural and induced fractures governs the injection and withdrawal of fluids from subsurface reservoirs. Connectivity depends on intersections that control how fluids mix and move through the entire system. Here, we present data sets from 3D X-ray microscopy measurements of simple fracture networks under stress. 3D printing was used to create prismatic blocks that formed fracture networks composed of 2 orthogonal fractures. The network orientation was either "x" or "+" relative to an applied vertical stress. 3D data sets were collected for normal loads of 25, 100 and 200 Newtons for samples with fracture surfaces with either correlated or uncorrelated asperity distributions. The file contains data from the 12 samples analyzed along with an example code used to extract the intersection geometry. Additional experimental details can be found in the manuscript "Geologic Stress Modulates Fluid Mixing at Fracture Intersections" (10.1038/s43247-026-03525-9)and supplemental information to appear in Communications Earth & Environment in 2026.

02 PETROLEUM↗

Destructive Examination Protocol for 3013-Container-Package Storing Plutonium-Bearing Materials

The 3013-container-package consists of a convenience, inner and outer container and is used for long-term storage of plutonium-bearing materials. A destructive examination (DE) protocol has been developed to examine the container package visually and with microscopic aid to find any corrosion conditions that could result in the loss of the integrity of the container package over its lifetime. The DE protocol contains three main steps: initial container examination, helium (He) leak testing, and detailed imaging and analysis of the inner container closure weld region (ICCWR). The ICCWR has been determined to be bounding, defined as exhibiting worst case conditions for stress corrosion cracking (SCC) of the inner container. To assess SCC in the ICCWR, the inner container lid is cut into quarters and the weld and He-leak testing gasket are removed. Then a citric acid wash is performed to remove adherent chlorides from the ICCWR. The wash is then sent for analysis to determine the concentration of chlorides in the ICCWR. While the analysis for chloride concentration is being performed, the quarter sections are further sectioned into 1/8th subsections by cutting each section in half. These subsections are washed using nitric acid to remove corrosion products. Then each subsection is imaged using a Wide Angle 3-D Measurement System (WAMS). After analysis of microscope images for potential SCC, additional imaging can be performed, including subsurface imaging. After review, a determination is made of whether the container integrity may potentially affect the safe storage of the material.

container↗

Modeling Geologic Waste Repository Systems Below Residual Saturation

The heat generated by high-level radioactive waste can pose numerical and physical challenges to subsurface flow and transport simulators if the liquid water content in a region near the waste package approaches residual saturation due to evaporation. Here, residual saturation is the fraction of the pore space occupied by liquid water when the hydraulic connectivity through a porous medium is lost, preventing the flow of liquid water. While conventional capillary pressure models represent residual saturation using asymptotically large values of capillary pressure, here, residual saturation is effectively modeled as a tortuosity effect alone. Treating the residual fluid as primarily dead-end pores and adsorbed films, relative permeability is independent of capillary pressure below residual saturation. To test this approach, PFLOTRAN is then used to simulate thermal-hydrological conditions resulting from direct disposal of a dual-purpose canister in unsaturated alluvium using both conventional asymptotic and revised, smooth models. Importantly, while the two models have comparable results over 100 000 years, the number of flow steps required is reduced by approximately 94%.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Insights into the year-round vertical distribution of chlorophyll concentration in high-latitude Arctic Ocean: implications for primary production

Climate-induced rapid changes in the Arctic Ocean, such as decreasing sea ice extent and increasing water temperature, are altering nutrient and light availability, profoundly impacting primary producer growth. However, access to the high-latitude Arctic Ocean is limited, and satellite data are primarily available only during summer, making continuous in-situ data collection challenging. We collected year-round chlorophyll-a (Chl-a) concentration data in high-latitude regions using a mooring system and performed a comparative analysis with reanalysis data. Unlike previous satellite-based studies, which typically rely on surface measurements, we used the annual vertical distribution of Chl-a. These data were applied to the vertically generalized production model to accurately estimate annual primary production. The moored Chl-a concentration data showed that phytoplankton exhibited a typical subsurface chlorophyll maximum (SCM) layer as sea ice retreated in June. Contrary to the gradually deepening SCM distribution predicted by model-based reanalysis data, the SCM layer persisted for approximately 4 months. This indicates that light and nutrient conditions within the SCM layer remained stable, sustaining continuous phytoplankton growth. Annual primary production, reflecting this vertical distribution of Chl-a concentration, was 6.85 gC m −2 yr −1 . This exceeded satellite-based estimates by at least two-fold, highlighting the significant underestimation of primary production by satellite approaches. Estimating primary production while accounting for the vertical distribution of phytoplankton and light is essential for improving ecological models to better understand carbon cycle and food web changes in the Arctic Ocean, with important implications for climate change predictions.

Arctic Ocean↗

De-Risking Exploration for Geothermal Plays in Magmatic Environments Through Open-Source Tools: An Open-Source Python Framework for 2D and 3D Play Fairway Analysis

The De-Risking Exploration for Geothermal Plays in Magmatic Environments (DEEPEN) project seeks to accelerate superhot geothermal development by reducing exploration risk through advanced open-source modeling tools. This work presents a novel Python-based framework, geoPFA, for conducting 2D and 3D play fairway analysis (PFA) tailored to superhot geothermal systems. Building on previous methodologies, the framework integrates thermo-hydro-mechanical-chemical simulation outputs from TReactMech, resulting in improved representation of subsurface properties that are critical to superhot resource producibility. 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 and suggest possible upflow from the Hengill volcanic system. The geoPFA library is publicly available, offering a scalable and reproducible approach to geothermal exploration across varied geological contexts.

15 GEOTHERMAL ENERGY↗

Dynamic Surface Restructuring in Cu(Au) Alloys Driven by Oxygen-Mediated Au Mobility

Alloying plays a crucial role in tuning the surface properties of metals, but the atomic-level mechanisms by which alloying elements influence surface structure dynamics under reactive conditions remain elusive. Using Cu(Au) in oxidizing environments as a model system, we reveal a dynamic oxygen-induced transformation of the topmost atomic layer into a periodically hill-and-valley morphology with reversible switching between undulated and flattened surface states. These interconversions are driven by the retreat of surface Au to the subsurface during oxygen adsorption and its resegregation to the surface upon oxygen desorption. This cyclical mobility establishes a feedback loop, allowing the surface to dynamically reconfigure in response to changes in the oxygen pressure. The results offer a broadly applicable framework for understanding atomic-scale surface restructuring in alloy systems, where differences in the chemical reactivity of alloying elements drive dynamic redistribution between surface and subsurface regions. As a result, this dynamic coupling has practical implications for designing corrosion-resistant coatings and metastable nanostructures with tunable catalytic properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Metallene-related materials for electrocatalysis and energy conversion

As a member of graphene analogs, metallenes are a class of two-dimensional materials with atomic thickness and well-controlled surface atomic arrangement made of metals or alloys. When utilized as catalysts, metallenes exhibit distinctive physicochemical properties endowed from the under-coordinated metal atoms on the surface, making them highly competitive candidates for energy-related electrocatalysis and energy conversion systems. Significantly, their catalytic activity can be precisely tuned through the chemical modification of their surface and subsurface atoms for efficient catalyst engineering. Here, this minireview summarizes the recent progress in the synthesis and characterization of metallenes, together with their use as electrocatalysts toward reactions for energy conversion. In the Synthesis section, we pay particular attention to the strategies designed to tune their exposed facets, composition, and surface strain, as well as the porosity/cavity, defects, and crystallinity on the surface. We then discuss the electrocatalytic properties of metallenes in terms of oxygen reduction, hydrogen evolution, alcohol and acid oxidation, carbon dioxide reduction, and nitrogen reduction reaction, with a small extension regarding photocatalysis. At the end, we offer perspectives on the challenges and opportunities with respect to the synthesis, characterization, modeling, and application of metallenes.

36 MATERIALS SCIENCE↗

Using AI to build a hydrobiogeochemical soil model

Soil water content is a function of inputs from precipitation and outputs via evaporation, transpiration, lateral flow, and vertical percolation, and is sensitive to biogeochemical processes. As such, soils serve as an ideal integrator of atmospheric, hydrological, and biogeochemical processes affecting the water cycle. In addition, soil water retention capacity, infiltration rates, and hydraulic conductivity can buffer or exacerbate the effects of extreme precipitation events (e.g., flooding, runoff, subsurface transport, erosion, greenhouse gas emissions) and mitigate the impact of droughts and heat waves on land systems (e.g., fire, crop failure). However, integrating water cycle measurements spanning different land atmosphere compartments across scales is a fundamental barrier for numerical model predictability. A significant challenge is that each domain (soil, hydrology, biology, and atmosphere) typically collects different sets of data at different temporal and spatial frequencies/scales, and even different dimensionalities (2D vs 3D). To implement soil as an integrator of the water cycle in land models, we suggest that novel machine learning (ML) tools can be developed to effectively simulate complex landscapes across various domains and scales, extended to regions with sparse or no data. The ultimate goals are to improve predictive understanding of land-atmosphere interactions and to extend the predictability of current Earth System Models (ESMs) through better integration of hydrological and biogeochemical data. We envision a framework in which: (1) ML-aided data reconstructions enable the merger of data sources into a unified geospatial product; (2) automated detection techniques are used to improve the knowledge of complex soil processes and interactions; and (3) this knowledge is leveraged and incorporated into models through AI-based emulators to distinctly connect the land and atmospheric compartments of the water cycle in models.

54 ENVIRONMENTAL SCIENCES↗

3D Geologic Framework Modelling of the Los Alamos National Laboratory Site and Pajarito Plateau: Integrating a realistic 3D fault network and modelling subsurface relationships in a sparsely sampled and complex geologic region

The subsurface geology beneath the Pajarito Plateau is critical to understanding the seismic hazard of the Pajarito Fault System, yet our understanding of this geology is relatively poor. While previous 3D geologic framework models of the area have been created for the purposes of understanding hydrogeologic flow, they are inadequate for the purposes of understanding the Pajarito Fault System. The specific challenges of using oil and gas software for this purpose include: (1) the geologic complexities resulting from volcanism and tectonism; (2) a need for a high level of stratigraphic detail over a large area; (3) a near complete lack of seismic data; and (4) sparse wellbore data. Presented here is a workflow that handles these challenges of adapting commercially available software used by the oil and gas industries to this seismic hazard problem.

58 GEOSCIENCES↗

Calcined Polyethyleneimine Coated Optical Fibers for Distributed pH Monitoring at High Pressures and Temperatures

In the oil and gas, CO2 sequestration, H2 subsurface storage, and geothermal energy sectors, subsurface pH measurements are critical for monitoring the geochemical conditions and estimating potential corrosion rates of wellbore systems. Real-time pH measurements in these conditions are vital for detecting and predicting corrosion deterioration of wellbore components that may jeopardize the safety and continued operation of wellbore systems. Building off of earlier TiO2 based pH sensors and the known pH sensitivity of amine-based polymers, a coating based on the secondary amine polymer polyethyleneimine (PEI) was developed. The sensor coating was treated with a high temperature (500 °C) calcination procedure in air to convert it into a more stable oxidized coating capable of withstanding hot aqueous solutions without dissolving while retaining linear pH sensitivity from pH values between 2 and 11. The sensor performance was measured using optical transmission measurements in solutions of various pHs and using optical backscatter reflectometry for distributed pH sensing demonstration in wellbore-relevant pressures (up to 1000 psi) and temperatures (80 °C). A calibration curve with strong differentiation between acidic and alkaline pH was developed for both transmission-based and distributed pH measurements, using fixed wavelength transmission and integrated linear amplitude of backscattered light for distributed measurements.

Shumski, Alexander↗

Calcined Polyethyleneimine-Coated Optical Fibers for Distributed pH Monitoring at High Pressures and Temperatures

In the oil and gas, CO2 sequestration, H2 subsurface storage, and geothermal energy sectors, subsurface pH measurements are critical for monitoring the geochemical conditions and estimating potential corrosion rates of wellbore systems. Real-time pH measurements in these conditions are vital for detecting and predicting corrosion deterioration of wellbore components that may jeopardize the safety and continued operation of wellbore systems. Previous tests using metal oxide-based coatings (TiO2) provided strong responses at elevated temperatures and moderate pressure stability but provided poor differentiation between acidic and alkaline solutions. Building off the pH responsiveness of the TiO2 surface and known pH sensitivity of amine-based polymers, a coating based on the secondary amine polymer polyethyleneimine (PEI) was developed. As the polymer itself is highly water soluble and easily removed by aqueous solutions, the sensor coating was treated with a high temperature (500 °C) calcination procedure in air to convert it into a more stable oxidized coating capable of withstanding hot aqueous solutions without dissolving while retaining linear pH sensitivity from pH values between 2 and 11. The sensor performance was measured using optical transmission measurements in solutions of various pHs and using optical backscatter reflectometry for distributed pH sensing demonstration in wellbore-relevant pressures (up to 1,000 psi) and temperatures (80 °C).

Shumski, Alexander↗

Multiscale Dynamics of Reactive Fronts in the Subsurface

Understanding and predicting flow and reactive transport in rocks (i.e. geologic porous media) is critical to many technologies at the heart of the energy transition, including CO 2 sequestration and H2 storage. However, accurate modeling and prediction of these systems is very complex because physico-chemical processes that occur at very small spatial scales, i.e. in the pores of the rocks, can dramatically control the system performance at the field scale (km). For example, precipitation reactions at the pore-scale can lead to large permeability changes at the field scale and dramatically alter the migration of stored gases. Properly accounting for multi-scale coupling effects is critical to achieve predictivity and confidence in model outputs, which then can guide design and optimization at the system-scale. This can be achieved through the development of rigorous mathematical models that can appropriately account for fine-scale effects at the large scale. The final report of the Early Career award DE-SC0019075 “Multiscale dynamics of reactive fronts in the subsurface” summarizes the mathematical, numerical and experimental advancements to study and predict reactive transport in geologic porous media across scales.

58 GEOSCIENCES↗

Stochastic fracture generation and thermo-hydro-mechanical modeling in an equivalent continuum framework for enhanced geothermal systems

Enhanced geothermal systems (EGS) involve fracturing low permeability material to establish well connectivity and then injecting and circulating fluid into the fractured subsurface for geothermal power production. Changes in fracture aperture from contraction of the cooling matrix rock may alter network connectivity and risk thermal short-circuiting. Thermo-hydro-mechanical (THM) models are a useful tool to study these processes. However, as fracture networks are complex, and data may be limited, fracture networks in THM models are often stochastically generated. Given reliance on stochastic fracture networks and THM modeling to represent the subsurface and assess productivity of EGS, increased understanding of the influence of such statistically derived fracture networks on flow and heat transport in THM models is needed. Here, a new fracture process model is developed in the reactive transport code PFLOTRAN to stochastically generate fracture families and simulate changes in fracture aperture over time due to temperature changes of the rock matrix. Sixty-four different fracture networks ranging from well to poorly-connected, are modeled in PFLOTRAN with and without mechanical processes (THM vs TH). Results indicate that for well-connected fracture networks, thermal short-circuiting is less of a concern due to the abundance of available alternative flowpaths. For poorly-connected fracture networks, inclusion of mechanical processes showed steep thermal drawdown coincident with increase in fracture aperture along developing colder flowpaths, demonstrating the risk of thermal short-circuiting. Simulations with additional, larger fractures engineered to establish connectivity in a poorly-fractured subsurface, indicate that while stochastic variation of fracture orientation of the background network had limited influence, such variation in the engineered fractures significantly affected flow and heat transport.

Discrete fracture networks (DFN)↗

A dynamic Bayesian optimized active recommender system for curiosity-driven partially Human-in-the-loop automated experiments

Optimization of experimental materials synthesis and characterization through active learning methods has been growing over the last decade, with examples ranging from measurements of diffraction on combinatorial alloys at synchrotrons, to searches through chemical space with automated synthesis robots for perovskites. In virtually all cases, the target property of interest for optimization is defined a priori with the ability to shift the trajectory of the optimization based on human-identified findings during the experiment is lacking. Thus, to highlight the best of both human operators and AI-driven experiments, here we present the development of a human–AI collaborated experimental workflow, via a Bayesian optimized active recommender system (BOARS), to shape targets on the fly with human real-time feedback. Here, the human guidance overpowers AI at early iteration when prior knowledge (uncertainty) is minimal (higher), while the AI overpowers the human during later iterations to accelerate the process with the human-assessed goal. We showcase examples of this framework applied to pre-acquired piezoresponse force spectroscopy of a ferroelectric thin film, and in real-time on an atomic force microscope, with human assessment to find symmetric hysteresis loops. It is found that such features appear more affected by subsurface defects than the local domain structure. This work shows the utility of human–AI approaches for curiosity driven exploration of systems across experimental domains.

36 MATERIALS SCIENCE↗