Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsurface flow”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Pore-Scale Analysis of the Permeability Damage and Recovery during Cyclic Freshwater and Brine Injection in Porous Media Containing Non-Swelling Clays

Permeability damage in subsurface porous media caused by clay mobilization is encountered in many engineering applications, such as geothermal energy, water disposal, oil recovery, and underground CO2 storage. During the freshwater injection into rocks containing brine, the sudden decrease in salinity causes native clay fines to detach and clog pore throats, leading to a significant decline in permeability. The clay fines detach due to weakened net-attractive forces binding them to each other and the grain. Past experiments link this permeability damage on the immediate history of the salinity and the direction of flow. To better understand this phenomenon, we conducted pore-scale simulations of cyclic injection of freshwater and brine into sandstone containing Kaolinite clay. Our simulations establish a link between the clay-fine trajectory and the permeability trend observed by Khilar and Fogler (1983). For a uniform clay size of 3 microns, we observe a permeability decline by two orders of magnitude during freshwater injection with respect to brine injection. Increasing salinity and simultaneously reversing flow direction restores the permeability. The permeability restoration upon reversing the brine flow direction is attributed to the unblocking of pore throats in the reverse direction by the movement of the clay particles along the grain surfaces by the hydrodynamic force and the strong net-attractive force under high salinity.

Bhuvankar, Pramod (ORCID:000000027739973X)↗
NASA NTRS1991

Large area QNDE inspection for airframe integrity

Quantitative Nondestructive Evaluation (QNDE) technology is being developed to provide new options for cost effective inspection of airframes. An R&D effort based on five NDE technologies is addressing questions of structural bonding assessment, corrosion detection, multisite damage detection, and fatigue characterization. The research/applications are being conducted by prioritized focussing and staging of the following technologies: (1) thermal NDE; (2) ultrasonic NDE; (3) coherent optical NDE; (4) magnetic imaging NDE; and (5) radiographic NDE. The focus here is on the most recent applications of thermal NDE technology to large area inspection of lap-joint and stiffener bonds. The approach is based on pulsed radiant heating of the airframe and measurement of the surface temperature of the structure with an infrared imager. Several advantages of the technique are that it is noncontacting, inspects one square meter area in a period of less than 2 minutes and has no difficulty inspecting typical curvatures of the fuselage. Numerical models of heat flow in these geometries are used to determine appropriate techniques for reduction of the infrared images, thereby delineating regions of disbonds. These models are also used to determine the optimum heating and measurement times for maximizing the contrast between bonded and unbonded structures. Good agreement is found between these results and experimental measurements, and a comparison of the two are presented. Also presented are results of measurements on samples with fabricated defects which show the technique is able to clearly indicate regions of disbonds. Measurements on an airframe also clearly image subsurface structure.

Winfree, William P.↗
DOE OSTI2024

Lagrangian Decomposition of the Atlantic Ocean Heat Transport at 26.5°N

Abstract The Atlantic Meridional Overturning Circulation (AMOC) plays a critical role in the global climate system through the redistribution of heat, freshwater and carbon. At 26.5°N, the meridional heat transport has traditionally been partitioned geometrically into vertical and horizontal circulation cells; however, attributing these components to the AMOC and Subtropical Gyre (STG) flow structures remains widely debated. Using water parcel trajectories evaluated within an eddy‐rich ocean hindcast, we present the first Lagrangian decomposition of the meridional heat transport at 26.5°N. We find that water parcels recirculating within the STG account for 37% (0.36 PW) of the total heat transport across 26.5°N, more than twice that of the classical horizontal gyre component (15%). Our findings indicate that STG heat transport cannot be meaningfully distinguished from that of the basin‐scale overturning since water parcels cooled within the gyre subsequently feed the northward, subsurface limb of the AMOC.

Tooth, Oliver J.↗
DOE OSTI2024

Modeled estimates of permafrost hydrology and related fields for pan-Arctic region over the period 1980-2100

This dataset contains gridded estimates of total runoff (mm/day); subsurface runoff (mm/day); snow water equivalent (SWE, mm/day); vertical upward water flux of combined surface evaporation, evapotranspiration, and snow sublimation (mm/day); and soil temperature (C) across the pan-Arctic basin from simulations of the Permafrost Water Balance Model (PWBM v4). These data were generated to allow researchers to advance understanding of the changes unfolding in the Arctic's terrestrial water cycle under the influences of climate warming, hydrological cycle intensification, and permafrost thaw. The dataset will enable researchers to address important questions spanning multiple disciplines. For example, how will rising river flows exported to the coast impact critical zone ecosystems, near-shore ice dynamics, and ocean biogeochemical processes. Spatial resolution is 25x25 km with N=39926 grid cells spanning the domain. Three forcing datasets were used to drive simulations at an implicit daily time step over the period 1980-2100. Daily air temperature, precipitation, and wind speed data from the WATCH Forcing Data methodology applied to ERA5 reanalysis (W5E5) were used for meteorological forcings over the period 1980-2019. Simulations for the period 1980-2100 used meteorological data from two global climate models; IPSL-CM6A-LR and MPI-ESM1-2-HR. The PWBM hydrological outputs are organized by grid and month in gzipped ASCII text files within each archived tar file. Soil temperature is organized by grid cell and day of year from DOY 205 to DOY 260. No special software is required to work with these data. The package contains spatial domain and river basin information files. These data are most useful for analyses of the spatial and temporal variations in terrestrial hydrology and soil thermal dynamics. A manuscript describing the model simulations and associated analysis has been accepted for publication in The Cryosphere (Rawlins and Karmalkar, 2024).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Mechanisms for Microseismicity Occurrence Due to CO 2 Injection at Decatur, Illinois: A Coupled Multiphase Flow and Geomechanics Perspective

Here, we numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO 2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologically consistent model that honors key stratigraphic horizons and 3D fault surfaces interpreted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO 2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO 2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection interval to the basement is essential to explain the destabilization of the regions where microseismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO 2 sequestration operations.

58 GEOSCIENCES↗
DOE OSTI2024

Creating microfluidic channels functionalized with micro- and nano-scale features via femtosecond laser surface processing

We demonstrate the ability to create functionalized microfluidic channels using femtosecond laser surface processing (FLSP). FLSP is an emerging advanced manufacturing technology used to modify the surface properties of materials directly and permanently by producing self-organized quasi-periodic micro- and nano-scale surface features along with surface and subsurface chemical and grain structure changes. We demonstrate on Hastelloy X that by controlling the laser fluence and pulse count, the depth of the microchannels and height of the FLSP microstructures within the microchannels can be controlled independently.

Ultrashort pulse laser applications, micro- and na↗
DOE OSTI2018

Devising a framework to Couple Electromagnetic Geophysical Measurements with Reservoir Simulations to Monitor CO2 Plume Movement

Monitoring CO2 plume movement and extent is critical to any CO2 utilization and storage (CCUS) project but can present significant challenges. Downhole measurements (such as cross-well seismic) can give the resolution needed to track the CO2 plume. That resolution comes at a significant expense in equipment and lost production time, and the resolution is only in the vicinity of that well. Surface seismic is less expensive but has decreasing resolution with depth and may not be sensitive to changes in the CO2 plume. Conversely, electrical conductivity has been shown to be correlated with changes in phase saturation associated with CCUS operations. A novel approach for subsurface monitoring that can be applied with minimal expense and operation disruption is charged wellbore casing controlled-source electromagnetics (CWC-CSEM). This technique injects electrical current into two legacy wells. Surface measurements of the induced electric and magnetic fields are obtained with distributed data loggers. Time-lapse measurements are sensitive to changes in CO2 plume movement, but inversion requires the application of prior knowledge from prior characterization efforts. We propose a methodology where simulations are used to inform CWC-CSEM inversion while simultaneously being validated. As such, uncertainty is reduced, and more robust flow models can be devised. A three-year study is underway at the Bell Creek Field in southeast Montana (DOE-NETL Project DE-FE0028320) to evaluate the effectiveness of CWC-CSEM on an active CO2-EOR operation. Determining the relationship between electrical resistivity and phase saturations, and coupling that to a reservoir model is the thrust of this project. Initially, the relative permeability/capillary pressure relationship and the compositional oil model were identified as model parameters with high uncertainty. To elucidate the impact of these uncertain parameters, simulations were run for ten years with a single five-spot pattern. Final water saturation distributions were converted into resistivity using Archie’s Law and available fluid sampling data to provide an exploration of the time lapse electromagnetic anomaly which can be anticipated using CWC-CSEM methods. This technique has promise as an inexpensive and accurate way to monitor CO2 injection in a CCUS site.

Moodie, Nathan↗
DOE OSTI2018

Detailed analysis of geothermal potential in Crescent Valley, north-central Nevada

This study has synthesized geological, geochemical, and geophysical data sets to elucidate the structural architecture and history of the Crescent basin in order to assess its geothermal potential. The Crescent basin is a complex composite basin, consisting primarily of multiple east-tilted half grabens and horst blocks bounded by several major faults of variable strike. Fault kinematic data indicate a present-day west-northwesttrending least principal stress, which in turn suggests that north-northeast-striking normal faults have relatively high slip and dilation tendency. However, a protracted history of changing Cenozoic extension directions has produced a mosaic of fault patterns, which has generated both segmented fault geometries and multiple fault intersections. These structural complexities have ultimately primed the basin for geothermal activity. Detailed analyses and application of the play fairway methodology have respectively denoted and ranked 12 favorable structural settings capable of hosting geothermal fluids in the Crescent basin. Of the 12 favorable structural settings, five were identified as having a higher likelihood of hosting a relatively high-temperature (>130°C) geothermal system, including Hand-Me-Down Creek, Hot Springs Point, and three other unnamed areas. The Hand-Me-Down Creek hybrid structural setting along the Crescent Valley fault is the only documented relatively high temperature geothermal system identified in the basin. It contains an active 87°C hot spring actively depositing sinter and geothermometry greater than 170°C (possibly >180°C). Future exploration efforts at Hand-Me-Down Creek would benefit from focusing on targeting high permeability zones at depth utilizing densely spaced electromagnetic surveys, fluid flow modeling, and selective drilling that targets specific fault intersections within the greater step-over at depth. The Hot Springs Point hybrid structural setting along the Dry Hills fault hosts a geothermal system, but geothermometry, spring temperatures, and well temperatures indicate that fluids do not exceed ~100°C. Three other relatively promising areas include a fault intersection in the central part of the basin and two step-overs (relay ramps) along the Crescent Valley fault, but additional exploration is needed at these sites, such as 2-m temperature surveys in the dry season, measuring deeper subsurface temperatures through temperature-gradient drilling, and collecting water samples through geoprobe drilling to estimate geothermometry.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Corrosion of 316 L stainless steel under the natural circulation of molten NaCl-MgCl 2 salt

The corrosion behavior of 316 L stainless steel (SS) was studied via the natural circulation of molten eutectic NaCl-MgCl 2 salt through a microloop. The post-corrosion tested 316 L SS microloop sections were characterized with microscopy techniques to determine the microstructural and microchemical changes that occurred at the alloy/salt interface. It was found that 316 L SS showed heterogeneous dissolution at the hot-leg, whereas deposition of corrosion products occurred at the cold-leg. For the first time, experimentally obtained molten salt flow-induced corrosion of 316 L SS results were combined with computational thermodynamic-kinetic models to validate the dissolution and deposition in terms of elemental compositional changes at the alloy/salt interface. The thermodynamic-kinetic modeling predicted that the heterogeneous dissolution of Cr and Fe from the hot-leg section of 316 L SS persisted throughout the salt circulation. The model also estimated that, despite Cr deposition starting earlier than Fe, the total redeposition of Fe is expected to be significantly greater than that of Cr over the circulation of salt. Furthermore, the modeling accurately predicted the subsurface enrichment of Mo which is attributed to the reduced Cr activity and the relatively higher diffusion rate of Mo within the alloy matrix. Here, the agreement between modeling and experimental results confirms that Fe chlorides dissolve at the hot-leg and subsequently deposit at the cold-leg due to activity changes driven by the thermal gradient. By contrast, Cr was not detected in the cold-leg deposits, which is attributed to its weaker temperature dependence on activity, limiting its redeposition under these conditions.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University: Preprint

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ("CUBO") was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗
DOE OSTI2025

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ('CUBO') was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗
DOE OSTI2024

Selecting Appropriate Model Complexity: An Example of Tracer Inversion for Thermal Prediction in Enhanced Geothermal Systems

Abstract A major challenge in the inversion of subsurface parameters is the ill‐posedness issue caused by the inherent subsurface complexities and the generally spatially sparse data. Appropriate simplifications of inversion models are thus necessary to make the inversion process tractable and meanwhile preserve the predictive ability of the inversion results. In this study, we investigate the effect of model complexity on fracture aperture inversion and thermal performance prediction in a field‐scale EGS model. Principal component analysis was used to map the aperture field to a low‐dimensional latent space. The complexity of the inversion model was quantitatively represented by the percentage of total variance in the original aperture fields preserved by the latent space. Tracer, pressure and flow rate data were used to invert for fracture aperture through an ensemble‐based inversion method, and the inferred aperture field was used to predict thermal performance. With an over‐simplified aperture model, ensemble collapse occurred. The inverted aperture models failed to resolve necessary flow and transport features, leading to a biased thermal performance prediction. A complex aperture model involved excessive features and was prone to overinterpreting the inversion data. Both the tracer/pressure/flow rate data reproduction and thermal prediction showed significant uncertainties, making it difficult to properly estimate long‐term thermal performance. Fortunately, our results indicate that there exists an appropriate model complexity which can simultaneously match inversion data and predict thermal performance with an acceptable uncertainty. The quality of the fit of tracer data appears to be a useful indicator of such an appropriate model complexity.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

A Boundary Element Model for Assessing Large‐Scale Pressurization in Faulted Geological Storage Systems

Assessing large-scale pressurization at the regional scale—a possible outcome of large subsurface storage applications such as wastewater injection and geological carbon sequestration—presents significant computational challenges. These challenges are particularly pronounced when accounting for complex geologic structures with multiple reservoir and caprock layers, fault zones, and wells. This study introduces a computationally efficient model that integrates single-phase semi-analytical solutions with a boundary element (BE) approach. The model simulates pressure propagation in multilayered 3D systems, including vertical faults, caprock, basement, and confining units. We apply this new model to a representative scenario involving CO 2 injection near a partially sealing fault with verification against an independent two-phase flow model. Results demonstrate that our model accurately captures far-field pressure responses and that, outside the CO 2 plume zone, pressure predictions from single-phase and two-phase models are nearly identical. This supports the use of single-phase models like ours for efficient estimation of far-field pressure changes. Additionally, we demonstrate its effectiveness at a large scale, incorporating multiple wells and faults. With its ability to represent multiple wells, fault zones, and geological heterogeneity, our model is well suited for assessments of basin-scale pressurization. Its computational efficiency also makes it a promising tool for integration with optimization frameworks aimed at designing and managing injection strategies in faulted storage systems.

Cihan, A. [Lawrence Berkeley National Laboratory (↗
DOE OSTI2025

Hydrologic response of artificially drained agricultural watersheds: insights from high-resolution integrated surface/subsurface simulations

Artificial drainage systems comprising subsurface networks of perforated pipes (tile drains) and engineered surface ditches are widely used to remove excess water from poorly drained agricultural regions. Artificial drainage lowers the water table by design but also has inadvertent effects on the watershed-scale hydrologic response with important implications for flood risk and nutrient exports. We investigated the effects of tile drains on watershed-scale hydrologic response in the Portage River, OH, Watershed using a high-resolution physically based integrated surface/subsurface hydrology model with recently developed capabilities to represent artificial drainage. Tile drains were found to enhance streamflow during times of low flow, generally consistent with previous studies. Streamflow flashiness was found to have a non-monotonic dependence on tile spacing with a minimum at intermediate spacings (∼50 m). Flashiness and the event hydrographs for small tile spacing were similar to the situation with no tiles, but flow paths from farm to stream were very different for those two end member cases, emphasizing the limitations of the stream hydrograph in characterizing hydrologic response. For typical tile spacings, peak flow can either be enhanced or attenuated by the presence of tiles, depending on the size of the event and the antecedent meteorological conditions. Tiles enhance peak flow when the event is below a threshold of ∼25 mm or when events arrive in dry conditions. Peak flow is reduced by tiles when events are large and arrive in conditions that are not overly dry. The dependence on event size and antecedent conditions is explained by differences in available storage and flow paths to the streams. These results provide additional insights into how tile drainage modulates event-scale hydrologic response, an important control on flood generation mechanisms and nutrient exports.

agricultural watershed↗
DOE OSTI2018

Above zone pressure interpretation for leaky well characterization and its identification from leaky caprock/fault

Pressure transient analysis has long been used for reservoir characterization. Above-zone (AZ) pressure has been recently investigated for inferring leakage pathway characteristics in leakage events from subsurface injection operations. The recorded pressure in the AZ should be purely related to leakage and therefore it can be safely inverted to deduce leakage characteristics. It is crucial to evaluate fluid leakage through abandoned wells to plan for further measurements of leakage prevention. However, the recorded AZ leakage signal may not be related to leaky well(s). Therefore, identification and spatial investigation of well leakage is required for leakage evaluation. Here in this paper, we propose a pressure interpretation method for early detection of leaky pathways, applying two observation points in the AZ. We distinguish leaky well, fault and caprock based on their corresponding flow regime identification. We show that the pressure difference of the two observation wells can be applied as a proxy for unknown leakage rate, which is crucial for leakage identification as well as characterization. Results show that the estimated location of the leaky well, leakage coefficient, and the leakage rate are in good agreement with the actual values. The estimated leakage coefficient of the leaky well can be used to evaluate well leakage in multiphase systems such as carbon dioxide leakage in deep saline aquifers.

02 PETROLEUM↗
DOE OSTI2020

Reactive chemical transport simulations of geologic carbon sequestration: Methods and applications

Chemical reaction simulations are considerably used to quantitatively assess the long-term geologic carbon sequestration (GCS), such as CO 2 sequestration capacity estimations, leakage pathway analyses, enhanced oil recovery (EOR) efficiency studies, and risk assessments of sealing formations (caprocks), wellbores, and overlying underground water resources. All these require a deep understanding of the CO 2 - associated chemical reactions. To ensure long-term, safe CO 2 sequestration in the intended formations, modeling is the only way to plausibly assess the CO 2 flow, reaction, and transport over thousands of years. This review summarizes the multiple methodologies for describing homogeneous and heterogeneous chemical reaction patterns and multiscale application examples, the recent progress and current status of chemical reaction simulations for GCS, and the impact of such simulations on geological CO2 sequestration performance. Technical gaps and future challenges are also discussed for further study. The trends and challenges of such studies include: (1) the combination of coupled chemical, mechanical, and transport processes with calibrated experiments and associated uncertainty/risk assessments; (2) enhancement of the ability to simulate detailed geophysical and geochemical equations to mimic in situ conditions; and (3) characterization of multiscale subsurface systems with detailed conceptual models and assignment of suitable boundary conditions for field-scale sequestration fields. One major gap remaining is the current lack of accurate (and scale-justified) kinetic and equilibrium chemical reaction parameters under reservoir conditions. Advanced models that couple chemical, mechanical, and transport processes with scale-justified parameters, from lab to field-scale experiments, are required for quantitative assessments of sequestration capacity and the long-term safety of GCS projects.

58 GEOSCIENCES↗
NASA NTRS2015

Subscale Validation of the Subsurface Active Filtration of Exhaust (SAFE) Approach to NTP Ground Testing

Brief History of NTP: Project Rover Began in 1950s by Los Alamos Scientific Labs (now Los Alamos National Labs) and ran until 1970s Tested a series of nuclear reactor engines of varying size at Nevada Test Site (now Nevada National Security Site) Ranged in scale from 111 kN (25 klbf) to 1.1 MN (250 klbf) Included Nuclear Furnace-1 tests Demonstrated the viability and capability of a nuclear rocket engine test program One of Kennedys 4 goals during famous moon speech to Congress Nuclear Engines for Rocket Vehicle Applications (NERVA) Atomic Energy Commission and NASA joint venture started in 1964 Parallel effort to Project Rover was focused on technology demonstration Tested XE engine, a 245-kN (55-klbf) engine to demonstrate startup shutdown sequencing. Hot-hydrogen stream is passed directly through fuel elements potential for radioactive material to be eroded into gaseous fuel flow as identified in previous programs NERVA and Project Rover (1950s-70s) were able to test in open atmosphere similar to conventional rocket engine test stands today Nuclear Furance-1 tests employed a full scrubber system Increased government and environmental regulations prohibit the modern testing in open atmosphere. Since the 1960s, there has been an increasing cessation on open air testing of nuclear material Political and national security concerns further compound the regulatory environment

Test Facilities↗