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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Antecedent Hydrologic Conditions Reflected in Stream Lithium Isotope Ratios During Storms

Antecedent hydrological conditions are recorded through the evolution of dissolved lithium isotope signatures (δ 7 Li) by juxtaposing two storm events in an upland watershed subject to a Mediterranean climate. Discharge and δ 7 Li are negatively correlated in both events,but mean δ 7 Li ratios and associated ranges of variation are distinct between them. We apply a previously developed reactive transport model (RTM) for the site to these event-scale flow perturbations, but observed shifts in stream δ 7 Li are not reproduced. To reconcile the stability of the subsurface solute weathering profile with our observations of dynamic stream δ 7 Li signatures, we couple the RTM to a distribution of fluid transit times that evolve based on storm hydrographs. The approach guides appropriate flux-weighting of fluid from the RTM over a range of flow path lengths, or equivalently fluid residence times. This flux-weighted RTM approach accurately reproduces dynamic storm δ 7 Li-discharge patterns distinguished by the antecedent conditions of the watershed.

58 GEOSCIENCES↗
DOE OSTI2025

Ubiquity and Causes of Soil Water Preferential Flow Across 17 Ecoregions

Abstract Preferential flow (PF) in soil causes the rapid transport of water, nutrients, and contaminants into the subsurface, influencing groundwater recharge and streamflow. Data scarcity has hindered the quantification of PF occurrence and the identification of its drivers across diverse ecoregions. We address this gap by analyzing high‐frequency, multi‐depth soil moisture data across 17 ecoregions in the USA, using ∼1,500 sensors at 40 sites. We discovered that PF is widespread, with sites experiencing PF in up to 60% of rainfall events ≥2 mm. Multiple approaches consistently show that PF is more likely to occur with increased peak rainfall intensity, finer textured material, low soil moisture variability, humid climate, and higher net primary productivity. This suggests that PF patterns could shift with projected climate changes, increasing uncertainty in predictions of groundwater recharge, water quality, and streamflow generation. Plain Language Summary Water can bypass part of the soil's matrix through a process called preferential flow (PF). This quick transport of water through the soil brings with it nutrients and contaminants and eventually makes it to groundwater and streams. To ensure ample amounts of good quality groundwater and surface water we need to understand when and where PF occurs. We inferred when PF occurred at 40 different sites across 17 ecoregions in the USA using soil moisture and rainfall data. We found that PF happened at all sites and in up to 60% of rainfall events ≥2 mm. Preferential flow was most likely at sites with high rainfall intensities, high clay content in soils, low variability in soil moisture, and high vegetation productivity. As rainfall intensities are predicted to increase due to climate change and vegetation becomes more productive, PF occurrence becomes more important for predicting groundwater recharge, water quality, and streamflow generation. Key Points Preferential flow (PF) is ubiquitous across the USA and occurs in up to 60% of all rainfall events ≥2 mm Rainfall intensity, soil texture, and antecedent soil moisture emerge as critical in generating PF across diverse ecoregions Two different PF detection approaches show similar relationships between key drivers and occurrence of PF

Li, Bonan↗
DOE OSTI2017

Analytical and variational numerical methods for unstable miscible displacement flows in porous media

The miscible displacement of one fluid by another in a porous medium has received considerable attention in subsurface, environmental and petroleum engineering applications. When a fluid of higher mobility displaces another of lower mobility, unstable patterns – referred to as viscous fingering – may arise. Their physical and mathematical study has been the object of numerous investigations over the past century. In writing this paper the objective is to present a review of these contributions with particular emphasis on variational methods. These algorithms are tailored to real field applications thanks to their advanced features: handling of general complex geometries, robustness in the presence of rough tensor coefficients, low sensitivity to mesh orientation in advection dominated scenarios, and provable convergence with fully unstructured grids.

42 ENGINEERING↗
DOE OSTI2026

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)↗
DOE OSTI2020

The Effects of Contaminant-Aging on Decontamination Efficacy for Rapid Remediation of Concrete Surfaces - 20412

Understanding the relationship between decontamination efficacy and contaminant aging is imperative for developing effective remediation strategies following a large-scale contamination event. Non-destructive decontamination of concrete following a contamination incident becomes more difficult with increased time between contamination and decontamination because contaminants may chemically bind to concrete and/or penetrate into the concrete subsurface. In this work, we evaluated the decontamination efficacy of two decontamination methods on concrete samples contaminated with soluble Cesium-137 (Cs-137) and silica particles with 0.5 μm and 2 μm diameters. Concrete samples were aged between 1 and 59 days, with half of the coupons receiving 1 mL of artificial rainwater about once every three days. After aging, coupons were either decontaminated or analyzed to determine contaminant penetration depths. Coupons were decontaminated with 0.1 M Potassium Chloride (KCI) solution applied either by flowing solution across the contaminated coupon face or non-destructive power washing. Contaminant depth profiles were created by removing the top surface of the coupon at least 15 times, measuring the activity in each removed layer, normalizing the total activity removed to the measured decontamination efficacy following the last layer removal, and determining the penetration thickness using the coupon dimensions, bulk density, and mass of material removed for each layer. Pressurized washing effectively removed particles, but not soluble cesium, despite the majority of Cs-137 being located within the first millimeter of the subsurface. Soluble cesium removal by flowing 0.1 M KCI solution across the contaminated surface dropped to less than 12% within the first ten days. This drop in Cs-137 removal efficacy occurred within the first five days for coupons experiencing artificial rainfall events. Flowing 0.1 M KCI solution across the contaminated face was ineffective at removing particles, likely because particles had localized in surface depressions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Aromatic amino acid metabolism and active transport regulation are implicated in microbial persistence in fractured shale reservoirs

Abstract Hydraulic fracturing has unlocked vast amounts of hydrocarbons trapped within unconventional shale formations. This large-scale engineering approach inadvertently introduces microorganisms into the hydrocarbon reservoir, allowing them to inhabit a new physical space and thrive in the unique biogeochemical resources present in the environment. Advancing our fundamental understanding of microbial growth and physiology in this extreme subsurface environment is critical to improving biofouling control efficacy and maximizing opportunities for beneficial natural resource exploitation. Here, we used metaproteomics and exometabolomics to investigate the biochemical mechanisms underpinning the adaptation of model bacterium Halanaerobium congolense WG10 and mixed microbial consortia enriched from shale-produced fluids to hypersalinity and very low reservoir flow rates (metabolic stress). We also queried the metabolic foundation for biofilm formation in this system, a major impediment to subsurface energy exploration. For the first time, we report that H. congolense WG10 accumulates tyrosine for osmoprotection, an indication of the flexible robustness of stress tolerance that enables its long-term persistence in fractured shale environments. We also identified aromatic amino acid synthesis and cell wall maintenance as critical to biofilm formation. Finally, regulation of transmembrane transport is key to metabolic stress adaptation in shale bacteria under very low well flow rates. These results provide unique insights that enable better management of hydraulically fractured shale systems, for more efficient and sustainable energy extraction.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2016

HORIZONTAL FLOWS IN ACTIVE REGIONS FROM RING-DIAGRAM AND LOCAL CORRELATION TRACKING METHODS

Continuous high-cadence and high spatial resolution Dopplergrams allow us to study subsurface dynamics that may be further extended to explore precursors of visible solar activity on the surface. Since the p-mode power is absorbed in the regions of high magnetic field, the inferences in these regions are often presumed to have large uncertainties. In this paper, using the Dopplergrams from space-borne Helioseismic Magnetic Imager, we compare horizontal flows in a shear layer below the surface and the photospheric layer in and around active regions. The photospheric flows are calculated using the local correlation tracking (LCT) method, while the ring-diagram technique of helioseismology is used to infer flows in the subphotospheric shear layer. We find a strong positive correlation between flows from both methods near the surface. This implies that despite the absorption of acoustic power in the regions of strong magnetic field, the flows inferred from the helioseismology are comparable to those from the surface measurements. However, the magnitudes are significantly different; the flows from the LCT method are smaller by a factor of 2 than the helioseismic measurements. Also, the median difference between the direction of corresponding vectors is 49°.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Effect of Shear Displacement and Stress Changes on Fracture Hydraulic Aperture and Flow Anisotropy

Fluid flow in fractures has the potential to drastically change the economic, environmental, and safety risks associated with a subsurface operation. This work focuses on the analysis of experimental shear fracturing and subsequent shear displacement of a shale specimen that contains multiple pre-existing natural fractures. A triaxial direct shear apparatus with concurrent µCT imaging allowed continuous monitoring of sample permeability and changes in fracture geometry during fracture creation, displacement, and changes in stress. Steady-state, low Reynolds number fluid flow simulations were used to determine the changes of each individual fracture’s hydraulic aperture due to displacement and effective stress. The development of anisotropy in the hydraulic aperture of a generated shear fracture was determined for shear displacements of 0.38, 0.74, and 2.04 mm. At high confining stress (30 MPa), shear displacement significantly reduced the hydraulic aperture of the shear fracture to below flow detectable limits and X-ray imaging resolution (23 × 23 × 23 µm 3 ). The shear generated fracture was found to be more sensitive to changes in effective stress compared to the pre-existing natural fractures. Existing models derived from fluid flow equations accounting for stochastic fracture roughness were able to partially fit the simulated flow properties of the pre-existing fractures, but no model was found to fit the flow properties of the newly formed shear fracture. These results are the first report of simultaneous changes in the geometry and flow properties of surrounding, parallel natural fractures as a shear fracture is generated with varying stress and shear displacement conditions.

58 GEOSCIENCES↗
DOE OSTI2021

Significance and complications of drilling fluid rheology in geothermal drilling: A review

The harsh downhole conditions of high pressure and high temperature (HPHT) encountered in geothermal wells make the drilling operation challenging. Drilling in such environments requires a special drilling mud formulation with high thermal stability and good rheological properties to fulfill the drilling fluid functions. Therefore, great efforts should be put into selecting the suitable drilling fluid, optimize and monitor the drilling fluid properties throughout drilling operations, and predicting its performance under downhole conditions. Rheological properties significantly impact many drilling parameters such as hole cleaning, fluid and wellbore stability, wellbore hydraulics, torque and drag, and other drilling issues. This paper discusses water-based drilling fluids' flow behavior under HPHT conditions and highlights the significance of fluid rheology in geothermal drilling. The common challenges and complications related to fluid rheology encountered in geothermal drilling are addressed in this paper, such as hole cleaning, wellbore hydraulics, and drilling fluid stability. This article also reviews the recent advances in drilling mud systems, rheology enhancement, and rheological properties measurements at surface and subsurface conditions. Furthermore, the rheology models of drilling fluid at elevated temperatures are reviewed to fully understand their flow behavior and establish a method for drilling engineers to optimize fluid formulations for geothermal drilling.

58 GEOSCIENCES↗
DOE OSTI2020

A review of spatial Markov models for predicting pre-asymptotic and anomalous transport in porous and fractured media

Heterogeneity across a broad range of scales in geologic porous media often manifests in observations of non-Fickian or anomalous transport. While traditional anomalous transport models can successfully make predictions in certain geological systems, increasing evidence suggests that assumptions relating to independent and identically distributed increments constrain where and when they can be reliably applied. A relatively novel model, the Spatial Markov model (SMM), relaxes the assumption of independence. The SMM belongs to the family of correlated continuous time random walks and has shown promise across a wide range of transport problems relevant to natural porous media. It has been successfully used to model conservative as well as more recently reactive transport in highly complex flows ranging from pore scales to much larger scales of interest in geology and subsurface hydrology. In this review paper we summarize its original development and provide a comprehensive review of its advances and applications as well as lay out a vision for its future development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Enhanced delivery of engineered Fe-Mn binary oxides in heterogeneous porous media for efficient arsenic stabilization

Heterogeneity in sediment and aquifer is universal, resulting in preferential flows of injected materials in the high permeability regions and forming flow by-passed zones in the low permeability regions during in-situ subsurface remediation. This adverse effect can considerably delay the completion of remedial operations and significantly increase the cost. Column experiments were designed and conducted to study the transport of starch- and starch-xanthan gum modified Fe-Mn binary oxide particles (SFM and SXFM) in saturated heterogeneous porous media and to reveal the particles’ arsenic (As) stabilization performance. Fine-in-Coarse (FIC) and Coarse-in-Fine (CIF) patterns of heterogeneous packings were set up in the columns. Testing results demonstrated that starch-xanthan gum dual treatment on Fe-Mn binary oxides successfully improved the particles’ migration capability in heterogeneous porous media and their distribution uniformity attributed to the profound shear thinning behavior of xanthan gum solution. The addition of xanthan gum to the system increased the viscosity and shear thinning property of SXFM suspension, making it a better candidate for delivery. Both SFM and SXFM stabilized As in heterogeneously packed sediment collected from a contaminated site, with SXFM showing better stabilization performance than SFM. The stabilization effects of SXFM were 90.7-97.0%, compared to 82.0-95.2% of SFM.

Yan, Xiulan↗
DOE OSTI2020

Morphological and nanomechanical changes in tungsten in high heat flux conditions

Abstract Morphological and nanomechanical alteration of tungsten in extreme environments, like those in edge localized modes in nuclear fusion environments, up to 46.3 GWm −2 heat fluxes were experimentally simulated using electrothermal plasma. Surface and subsurface damage to the tungsten is seen mainly in the form of pore formation, cracks, and resolidified melt instabilities. Mirco voids, rosette-type microfeatures, core-shell structure, particle enrichment, and submicron channels all manifest in the damaged subsurface. The formation of voids in the subsurface was determined to originate from the ductile fracture of hot tungsten by plastic flow but not developed to cracking. The voids were preferentially settled in grain boundaries, interfaces. The directionality of elongated voids and grains is biased to the heat flow vector or plasma pathway, which is the likely consequence of the thermally driven grain growth and sliding in the high-temperature conditions. The presence of a border between the transient layer and heat-affected zone is observed and attributed to plasma shock and thermal spallation of fractural tungsten at high temperature. Plasma peening-like hardening effects in tungsten were observed in the range of 22.7–46.3 GWm −2 but least in the case of the lowest heat flux, 12.5 GWm −2 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Increasing freshwater and dissolved organic carbon flows to Northwest Alaska’s Elson lagoon

Manifestations of climate change in the Arctic are numerous and include hydrological cycle intensification and permafrost thaw, both expected as a result of atmospheric and surface warming. Across the terrestrial Arctic dissolved organic carbon (DOC) entrained in arctic rivers may be providing a carbon subsidy to coastal food webs. Yet, data from field sampling is too often of limited duration to confidently ascertain impacts of climate change on freshwater and DOC flows to coastal waters. This study applies numerical modeling to investigate trends in freshwater and DOC exports from land to Elson Lagoon in Northwest Alaska over the period 1981-2020. While the modeling approach has limitations, the results point to significant increases in freshwater and DOC exports to the lagoon over the past four decades. The model simulation reveals significant increases in surface, subsurface (suprapermafrost), and total freshwater exports. Significant increases are also noted in surface and subsurface DOC production and export, influenced by warming soils and associated active-layer thickening. The largest changes in subsurface components are noted in September, which has experienced a ~50% increase in DOC export emanating from suprapermafrost flow. Direct coastal suprapermafrost freshwater and DOC exports in late summer more than doubled between the first and last five years of the simulation period, with a large anomaly in September 2019 representing a more than fourfold increase over September direct coastal export during the early 1980s. These trends highlight the need for dedicated measurement programs that will enable improved understanding of climate change impacts on coastal zone processes in this data sparse region of Northwest Alaska.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Transient flow modeling in fractured media using graphs

In this work, we describe a method to simulate transient fluid flows in fractured media using an approach based on graph theory. Our approach builds on past work where the graph-based approach was successfully used to simulate steady-state fluid flows in fractured media. We find a mean computational speedup of the order of 1400 from an ensemble of a 100 discrete fracture networks in contrast to the O(10 4 ) speedup that was obtained for steady-state flows earlier. However, the transient flows considered here involve an additional degree of complexity that was not present in the steady-state flows considered previously with a graph-based approach, that of time marching and solution of the flow equations within a time-stepping scheme. We verify our method with an analytical test case and demonstrate its use on a practical problem related to fluid flows in hydraulically fractured reservoirs. By enabling the study of transient flows, we create an opportunity for a wide set of possibilities where a steady-state approximation is not sufficient, such as the example motivated by hydraulic fracturing that we present here. This work validates the concept that graphs are able to reliably capture the topological properties of the fracture network and serve as effective surrogates in an uncertainty-quantification framework.

58 GEOSCIENCES↗
DOE OSTI2025

Experimental Investigation of Buoyant Flow in Realistic Bedforms With Heterogeneous Wettability

Submeter-scale geologic heterogeneity greatly affects CO 2 plume migration and retention. In this work, we present meter-scale laboratory experiments that can capture the impact of realistic submeter-scale geologic heterogeneity on multiphase flow and trapping. We produce realistic sedimentary formations consisting of ripple deposits with varying grain size contrast and wettability in a meter-scale slab chamber. Then, we conduct multiphase flow experiments with analog fluids through these structures and measure the saturation patterns, capillary heterogeneity trapping (CHT), and overall trapping performance. When we alter the ripple bedform architecture, variations in trapped saturation and CHT (10–20%) increment are exhibited. Similar growth in trapping performance is also observed when grain size contrast increases. Finally, wettability changes (water- to oil-wet) can increase nonwetting saturation and CHT up to 5% and 10–20%, respectively. These results emphasize the importance of correctly characterizing the impact of small-scale heterogeneities and wettability changes. We believe this is the first time that multiphase flow experiments were conducted in meter-scale domains with realistic ripple bedforms and heterogeneous wettability to investigate plume migration and trapping.

58 GEOSCIENCES↗
DOE OSTI2015

Southeast Regional Carbon Sequestration Partnership (SECARB): Completed History Match: Subsurface Plume Development, Early Test

In this paper, we compare 4D seismic interpretations of CO 2 plume evolution with fluid flow numerical simulation results for the SECARB Early test at Cranfield, Mississippi. Historic pressure trends, oil and gas production rates and current CO 2 -EOR production data from the field were history matched and tuned model used for predictive simulations. For CO 2 -EOR operations the tuned model results showing CO 2 plume distribution and CO 2 first arrival (break-through) times in production wells were compared to the available field data. Reasonable matches were observed with breakthrough times. However 4D seismic interpretations of saturated region distributions and simulation results mismatch in some significant ways, likely showing limitations in 4D seismic interpretations for Cranfield. The tuned model then used to predict reservoir response and storage capacity in different field development scenarios. We compared hypothetical scenarios where the operator decided to transition into CO 2 injection without oil production (CO 2 -EORT) when oil production is not economical versus a scenario of continuing with CO 2 -EOR. Our results show that CO 2 -EOR can store more CO 2 and operations will last longer whereas if switched to CO 2 -EORT field has to be abandoned earlier because of out-of-area movement of the CO 2 plume. However, the amount of CO2 stored per year during the operational period is larger for CO 2 -EORT.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2019

Mathematical Modeling on the Influence of Casting Parameters on Initial Solidification at the Meniscus of Slab Continuous Casting

In the current study, the two-dimensional mold model was applied to investigate the influence of different casting parameters on the initial solidification at the meniscus of the continuous casting mold. The profile of the meniscus and the slag rim, and the heat flux on the hot face of the copper plate, were compared with different locations of the mold and different casting parameters. The depth of the oscillation marks and solidified meniscus through simulation was measured with different casting parameters. With the mold oscillating downward from the peak to the valley, the maximum heat flux on the hot face of the copper plate near the meniscus increased from 4.3 to 5.5 MW/m{sup 2}, the thickness of the slag rim decreased from 4.7 to 4.2 mm, and the height of the curved meniscus also decreased. With the increase of casting speed from 1.2 to 1.45 m/min, the maximum speed of the upper backflow increased from 0.6 to 0.8 m/s, and the lowest location of the meniscus at approximately 310 mm from the copper plate decreased from 22 to 12 mm below the meniscus near the nozzle, where the slag entrapment easily occurred. With the increase of the casting speed, casting superheat, and oscillation frequency, and the decrease of the water flow rate in the mold and oscillation amplitude, the thickness of the slag rim decreased. With the increase of casting speed and superheat, and the decrease of the water flow rate in the mold, oscillation frequency, and oscillation amplitude, the depth of the solidified meniscus decreased. The regression equation was fitted to estimate the depth of hooks at the subsurface of slabs.

36 MATERIALS SCIENCE↗