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

Fractional Modeling in Action: A Survey of Nonlocal Models for Subsurface Transport, Turbulent Flows, and Anomalous Materials

Modeling of phenomena such as anomalous transport via fractional-order differential equations has been established as an effective alternative to partial differential equations, due to the inherent ability to describe large-scale behavior with greater efficiency than fully-resolved classical models. In this review article, we first provide a broad overview of fractional-order derivatives with a clear emphasis on the stochastic processes that underlie their use. We then survey three exemplary application areas – subsurface transport, turbulence, and anomalous materials – in which fractional-order differential equations provide accurate and predictive models. For each area, we report on the evidence of anomalous behavior that justifies the use of fractional-order models, and survey both foundational models as well as more expressive state-of-the-art models. We also propose avenues for future research, including more advanced and physically sound models, as well as tools for calibration and discovery of fractional-order models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

On the emergent scale of bedrock groundwater contribution to headwater mountain streams

We investigated the contribution of bedrock groundwater to streamflow as a function of catchment scale in a headwater stream. Synoptic surveys were conducted during hydrologically important periods of the year using multiple environmental tracers in stream water, soil water, and bedrock groundwater, along a first-order montane stream, in west-central Montana. Sampled analytes included 222 Rn, used to constrain total subsurface flux, and major and minor elements, used in end-member mixing analysis (EMMA) to identify the contributions of soil and bedrock groundwater to the stream. Partitioning between soil-derived and bedrock-derived groundwater was then analyzed as a function of the incremental and accumulated sub-catchment sizes. Radon results indicated that subsurface water contributions accounted for the majority of streamflow at all surveyed times. EMMA results revealed that the bedrock groundwater contribution to streamflow varied between 26% during peak snowmelt and 44% during late summer. Streamflow generation was dominated by soil groundwater contribution along the entire reach, but the bedrock groundwater contribution increased consistently with accumulated sub-catchment size. However, groundwater contributions were not well-correlated with incremental sub-catchment size. The scale at which increased bedrock groundwater discharge can be correlated with sub-catchment size appears to be >1 km 2 for our study. Our results are consistent with a conceptual model where streamflow is predominantly generated by a 3D subsurface nested flow system. Local subsurface heterogeneities control the stream source at local scales but begin to average out at scales >2 km 2 . Our study indicates that, while soil groundwater is the dominant source, bedrock groundwater remains an important and predictable contributor to streamflow throughout the year, even in a snow-dominated, mountainous headwater catchment.

environmental tracers↗
DOE OSTI2023

TOUGH3/EOS7MG (EOS7MG) v1.1

EOS7MG is a new TOUGH3 module for modeling water, brine, and multiple gas mixtures in the subsurface multiphase flow system. The module allows mixture of as more as three gas components selecting from a list of predefined gases (CH4, C2H6, C3H8, H2S, CO2, N2, NO2, H2, and AIR). EOS7MG was developed based on EOS7C and inherited the general modeling capabilities of the TOUGH3 platform. It uses cubic equation of states for gas density, viscosity, and enthalpy calculation. The module provides the selection options for Peng-Robinson, Redlich-Kwong, or Soave-Redlich-Kwong equations of state. The subsurface flow and transport of gas and aqueous phase mixtures over a wide range of pressures and temperatures is modeled with a multiphase Darcy's Law and an accurate chemical equilibrium approach as originally implemented in EOS7C is used for calculation of the solubilities of different mass components. EOS7MG was designed to allow flexible number of components based on the number of gases included in the model. The potential application of this module may include subsurface gas storage, energy storage, CCS, enhanced gas recovery, and more.

Zhang, Keni↗
DOE OSTI2025

The impact of capillary heterogeneity on CO 2 flow and trapping across scales

Capillary heterogeneity has been identified over the last decade as a key control on subsurface CO 2 flow behavior during geological CO 2 sequestration. These heterogeneities can be formed in all sedimentary rocks, ranging from slight variations in the sand grain sizes to extensive sequences of interbedded sands, shales, and limestones. Capillary heterogeneity has been largely, although not entirely, overlooked in subsurface flow modeling because it is assumed to only directly influence fluid redistribution over scales of centimeters to meters. However, even small-scale fluid movements can result in dramatic impacts on the mobility and trapping of the CO 2 over kilometers. Therefore, neglecting capillary heterogeneity at multiple scales could potentially lead to errors in modeling and predicting field-scale plume migration. In this review paper, we aim to provide a consistent overview to (1) establish that capillary heterogeneity can have a major impact on CO 2 plume migration, (2) establish the respective length scales at which capillary heterogeneity matters, and (3) provide guidance for numerical modeling. This review covers pertinent literature and extracts key observations from the core to the field scales. Experimental studies have shown that millimeter-decimeter scale capillary heterogeneity can cause the so-called capillary heterogeneity trapping in addition to pore-scale residual trapping. Even at such a small scale, capillary heterogeneity can already lead to complex upscaled constitutive relationships, such as flow-rate dependent and anisotropic relative permeability, which affects field-scale CO 2 migration even when field-scale heterogeneities are present. Under gravity-dominated flow regimes, centimeter-meter scale capillary heterogeneity can entrap a significant amount of CO 2 at field scale, not just after imbibition but also during drainage. In certain cases, the presence of capillary heterogeneity can even completely stop the vertical movement of the CO 2 plume, hence greatly reducing leakage risks. At meter-kilometer scale, the influence of capillary heterogeneity is more pronounced and can hinder or redirect CO 2 migration in both lateral and vertical directions. The impact of capillary heterogeneity across multiple spatial scales poses a great challenge in modeling CO 2 migration at field scale, because it is practically impossible to build a field-scale earth model with grid blocks at millimeter scale. We recommend a hierarchical modeling approach to address this challenge. At field scale, earth models are built to capture geological features and heterogeneities in high but still practical grid resolutions. For each facies or rock type of the field-scale model, high- resolution meter-scale “conceptual” models are built with millimeter-scale grid blocks to capture representative fine-scale bedding geometries and heterogeneities in various environments of deposition, bridging the gap from subcore scale to the size of a field-scale simulation grid block. Upscaling is then used to preserve the smaller-scale flow dynamics of various rock types in field-scale simulations. Here, future work is needed to (1) refine, improve, and validate the hierarchical modeling approach; (2) build libraries of fine-scale bedding models for facies in various environments of deposition; (3) quantify multiscale capillary heterogeneity effects under subsurface uncertainties; (4) gain learning from different storage formations; and (5) establish best practices that balance accuracy and computational speed.

Capillary heterogeneity↗
DOE OSTI2020

Relating the depth of the water table to the depth of weathering

Weathering of bedrock creates and occludes permeability, affecting subsurface water flow. Often, weathering intensifies above the water table. On the other hand, weathering can also commence below the water table. To explore relationships between weathering and the water table, a simplified weathering model for an eroding hillslope was formulated that takes into account both saturated and unsaturated subsurface water flow (but does not fully account for changes in dissolved gas chemistry). The phreatic line was calculated using solutions to mathematical treatments for both zones. In the model, the infiltration rate at the hill surface sets both the original and the eventual steady-state position of the water table with respect to the weathering reaction front. Depending upon parameters, the weathering front can locate either above or below the water table at steady state. Erosion also affects the water table position by changing porosity and permeability even when other hydrological conditions (e.g. hydraulic conductivity of parent material, infiltration rate at the surface) do not change. The total porosity in a hill (water storage capacity) was found to increase with infiltration rate (all else held constant). Furthermore, this effect was diminished by increasing the erosion rate. We also show examples of how the infiltration rate affects the position of the water table and how infiltration rate affects weathering advance.

58 GEOSCIENCES↗
DOE OSTI2022

Physically rigorous reduced-order flow models of fractured subsurface environments without explosive computational cost

Fractured media models comprise discontinuities of multiple lengths (e.g. fracture lengths and apertures, wellbore area) that fall into the relatively insignificant length scales spanning millimeter-scale fractures to centimeter-scale wellbores in comparison to the extensions of the field of interest, and challenge the conventional discretization methods imposing highly-fine meshing and formidably large numerical cost. By utilizing the recent developments in the finite element analysis of electromagnetics that allow to represent material properties on a hierarchical geometry, this project develops computational capabilities to model fluid flow, heat conduction, transport and induced polarization in large-scale geologic environments that possess geometrically-complex fractures and man-made infrastructures without explosive computational cost. The computational efficiency and robustness of this multi-physics modeling tool are demonstrated by considering various highly-realistic complex geologic environments that are common in many energy and national security related engineering problems.

42 ENGINEERING↗
DOE OSTI2025

Changes in soil water content and lateral flow exert large effects on soil thermal dynamics across Alaskan landscapes

Both lateral surface and subsurface water flow affect soil moisture dynamics, yet most land surface models only solve subsurface water movement vertically. Here, we use a 3D ecosystem model that considers both land surface and subsurface hydrologic processes to simulate soil moisture, which is then used to drive a 1-D vertical soil thermal model to simulate the soil moisture effects on soil thermal dynamics in central Alaska. Our coupled model improves soil temperature (ST) estimates by 43.5% in comparison with observational data. Soil moisture has little effect on ST during the wet season (-1.5%) and a substantial influence during the dry season (60%). Spatially, water lateral flow has significant impacts on both soil moisture and ST, causing model estimates for thawed areas in the transition season to increase by ~10% in the study area. Our results highlight the importance of considering dynamical soil moisture, as well as lateral flow effects, on soil thermal dynamics in permafrost regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The Complexity of Nonlinear Flow and non-Fickian Transport in Fractures Driven by Three-Dimensional Recirculation Zones

Subsurface fluid flow and solute transport are typically described by first-order or linear rate laws. However, deviations from the first-order rate laws (or anomalous behavior) are typical and lead to nonlinear flow and non-Fickian transport phenomena. The shared underlying mechanisms for these “anomalous” flow and transport phenomena, which are actually typical of fractures, are seldom analyzed. We investigated the co-occurrence of nonlinear flow and non-Fickian transport through numerical simulation experiments of pore-scale flow and solute transport processes through three-dimensional natural fractures. In these fractures, 3-D recirculation zones (RZs) developed and grew with increasing flow rates. Flow nonlinearity resulted from the RZ development and growth, which shrank the main flow channel. Moreover, the same RZs also captured solutes from the main flow channel and released them back later. This retardation resulted in non-Fickian transport. Thus, nonlinear flow and non-Fickian transport are intertwined via their shared dependence on RZs.

58 GEOSCIENCES↗
DOE OSTI2024

Development of inter-grid-cell lateral unsaturated and saturated flow model in the E3SM Land Model (v2.0)

Abstract. The lateral transport of water in the subsurface is important in modulating terrestrial water energy distribution. Although a few land surface models have recently included lateral saturated flow within and across grid cells, it is not a default configuration in the Climate Model Intercomparison Project version 6 experiments. In this work, we developed the lateral subsurface flow model within both unsaturated and saturated zones in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2.0). The new model, called ELMlat, was benchmarked against PFLOTRAN, a 3D subsurface flow and transport model, for three idealized hillslopes that included a convergent hillslope, divergent hillslope, and tilted V-shaped hillslope with variably saturated initial conditions. ELMlat showed comparable performance against PFLOTRAN in terms of capturing the dynamics of soil moisture and groundwater table for the three benchmark hillslope problems. Specifically, the mean absolute errors (MAEs) of the soil moisture in the top 10 layers between ELMlat and PFLOTRAN were within 1 %±3 %, and the MAEs of water table depth were within ±0.2 m. Next, ELMlat was applied to the Little Washita experimental watershed to assess its prediction of groundwater table, soil moisture, and soil temperature. The spatial pattern of simulated groundwater table depth agreed well with the global groundwater table benchmark dataset generated from a global model calibrated with long-term observations. The effects of lateral groundwater flow on the energy flux partitioning were more prominent in lowland areas with shallower groundwater tables, where the difference in simulated annual surface soil temperature could reach 0.3–0.4 ∘C between ELMv2.0 and ELMlat. Incorporating lateral subsurface flow in ELM improves the representation of the subsurface hydrology, which will provide a good basis for future large-scale applications.

58 GEOSCIENCES↗
DOE OSTI2022

Capillary Pressure to Relative Permeability: Task 7.3.4 Relative Permeability Analysis Final Report

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability relationship can be determined. Relative permeability is a numerical description of the interplay between two or more fluids and the porous media they flow through. It is a critical parameter for various tools used to characterized subsurface multiphase flow systems, such as numerical simulation for oil and gas development, carbon sequestration, and groundwater contamination remediation. Therefore, it is critical to get a good statistic distribution of relative permeability in the porous media under study. Empirical relationships for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility required to match laboratory derive relative permeability relationships. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability relationships. Existing laboratory-measured relative permeability data along with mercury intrusion capillary (MICP) data coupled with interfacial tension and contact angle measurements were used to determine the efficacy of this approach to relative permeability curve creation. The relative permeability relationships determined with this method were fit to existing laboratory data to elucidate common fitting parameters that were then used to create relative permeability relationships from MICP data that does not have an associated laboratory-measured relative permeability relationship.

58 GEOSCIENCES↗
DOE OSTI2020

Preliminary Study of Potential Utilization of Triassic Rift Basins for Long-Term Carbon Sequestration

Triassic rift basins of the eastern United States present a potential option for long–term carbon dioxide (CO 2 ) sequestration. Exposed and buried basins are located near large point CO 2 sources. Because of their similar origins, comparable fill successions indicate a level of reproducibility that can be conveyed between basins. The thick rock successions of the exposed Culpeper and Gettysburg basins served as proxies for understanding of the basin sequences. Their study allowed recognition of five mappable assemblages of rock types, herein termed lithofacies associations. These lithofacies associations were formed by alluvial fan, braided and meandering streams and marginal and distal lake depositional processes.Stratigraphic architecture of the lithofacies associations within exposed basins suggests a vertical succession that consists of an initial coarse-grained fluvial succession that is progressively replaced upward and basinward by finer grained lacustrine deposits. Along the faulted margins of the basins alluvial fans transitioning to fluvial delta deposits reflect lateral progradation and contemporaneous lateral components of the basin filling.To test the model devised in exposed basins, facies associations were developed for the buried Taylorsville Basin. This basin preserves over 8,000 feet of Triassic rocks that are concealed beneath more than 2,000 feet of Cretaceous and Tertiary Coastal Plain sediments. Composition and stratigraphic architecture of these rocks are similar to those observed in exposed basins. Infilling of the basin was the result of vertical aggradation and lateral progradation of coarse-grained to fine-grained facies. Based upon study of outcrops in the Culpeper and Gettysburg basins, groups of recurring lithologic facies were identified as the fundamental constructs of the basin’s sediment infilling. These groups of facies, termed lithofacies associations, represent an amalgamation of lithologic components from broadly similar depositional systems. These lithologic associations were then extrapolated into the exposed and buried portions of the Taylorsville basins. This effort determined that Triassic rift basins provide several avenues for potential study of geologic sequestration of carbon. The characteristic rift basin succession in the basins presents possible reservoir targets within the marginal fluvial deposits. Furthermore, intrusive and extrusive mafic bodies present a potential source of fracture porosity that could serve as CO 2 reservoirs. Porosity values from thin section analysis of exposed basin samples are higher in alluvial fan and braided fluvial lithofacies. Porosity is highest in samples that show less compaction, due to a lack of ductile lithic fragments and/or higher stratigraphic position. Log data from the Taylorsville Basin indicates that porosity and permeability values in basin marginal fluvial strata are elevated. Near the center of the basin, porosity and permeability values are greatly reduced, primarily owing to the fine-grained character of the lacustrine deposits. Some Triassic basins also contain thick intervals of concordant extrusive and intrusive mafic igneous rocks. These igneous bodies are potential CO 2 reservoirs for several reasons. Firstly, extrusive lava flows provide potential storage in the layers of primary porosity that occur in the vesicules formed at the top of lava flows. Secondly, both lava flows and subsurface igneous sills exhibit extensive fracture porosity produced by the rapid cooling of the flows and intrusions. Thirdly, these igneous rocks are mafic in composition, and studies have shown that iron- and magnesium-rich mafic rocks provide sequestration opportunities through carbonate remineralization. Lastly, extrusive and intrusive igneous rocks are invariably preserved within fine-grained lake deposits. These lacustrine sediments can serve as a fine-grained confining layer that encases the igneous rocks both above and below.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Flow Directions and Ages of Subsurface Water in a Salt Marsh System Constrained by Isotope Tracing

Salt marshes are dynamic hydrologic systems where terrestrial groundwater, terrestrial surface water, and seawater mix due to bi-directional flows and pressure gradients. Due to the counteracting terrestrial and marine forcings that control these environments, we do not comprehensively understand water fluxes in these complex coastal systems. To understand the water sources, flow directions, and velocities in salt marsh porewater, we employed a combination of geochemical tracers and analytical models across a hillslope-to-salt marsh continuum in a salt marsh experiencing daily inundation of estuarine surface water (SW) from tides and mixing of fresh seasonal groundwater. We used tritium ( 3 H) as a hydrologic tracer to assess porewater ages and stable water isotope ($\delta$ 2 H and $\delta$ 18 O) analyses to separate isotopically distinct estuarine and terrestrial groundwater across different depths and landscape positions in the study transect. We employed electrical conductivity to constrain the role of source mixing and evapotranspiration in salt marsh hydrology. Salinity and stable isotopes revealed that transpiration, rather than evaporation, increased subsurface water salinity to concentrations above estuarine SW during summer. Elevated salinity at depth indicated that salt marsh subsurface water is recharged during the dry growing season. Seasonal recharge patterns drive long-term deep subsurface water dynamics across the salt marsh, with 3 H ages of 3–7 years, and daily tidal cycles drive short-term shallow porewater dynamics with 3 H ages of 0 ± 3.6 years. Furthermore, our conceptual understanding of the spatiotemporal changes in SW-subsurface water interactions at the terrestrial-marine interface quantifies the hydrological constraints we are missing to improve our understanding of biogeochemical cycles within the salt marsh.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

4D Electrical Resistivity Imaging of Stress Perturbations Induced During High–Pressure Shear Stimulation Tests

Fluid flow through fractured media is typically governed by the distribution of fracture apertures, which are in turn governed by stress. Consequently, understanding subsurface stress is critical for understanding and predicting subsurface fluid flow. Although laboratory–scale studies have established a sensitive relationship between effective stress and bulk electrical conductivity in crystalline rock, that relationship has not been extensively leveraged to monitor stress evolution at the field scale using electrical or electromagnetic geophysical monitoring approaches. In this paper we demonstrate the use time–lapse 3–dimensional (4D) electrical resistivity tomography to image perturbations in the stress field generated by pressurized borehole packers deployed during shear–stimulation attempts in a 1.25 km deep metamorphic crystalline rock formation.

58 GEOSCIENCES↗
DOE OSTI2024

Fusing 4D coded thermoacoustic, electromagnetic, and acoustic/seismic wavefields for subsurface characterization and imaging of fluid flow in porous media

This research program aims to address critical challenges in subsurface exploration and monitoring of anthropogenic CO 2 storage by advancing quantitative dynamic sensing and imaging technologies. Two primary applications—petrophysical assessment of hydrocarbon reservoirs (C1) and CO 2 storage monitoring (C2)—require significant innovation to overcome three barriers (B1-B3). These barriers include the need for enhanced understanding of interactions between physical wave fields (thermoacoustic, electromagnetic, acoustic/seismic, and X-ray) with fluid-filled porous media, the development of multi-sensor data fusion for real-time imaging, and upscaling microscopic quantum effects for macroscopic observations. The project’s objective was to establish a unified 4D coded sensing and imaging approach, integrating EM, AC/S, and TA fields for multi-scale and multi-physics material characterization and subsurface imaging.

58 GEOSCIENCES↗