Linear & Nonlinear Solvers for Simulating High-Temperature Multiphase Flow Within Large-scale Engineered Subsurface Systems.
Abstract not provided.
Park, Heeho↗
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.
Integrating the Interconnections Between Groundwater and Land Surface Processes Through the Coupled NASA Land Information System and ParFlow Environment
Understanding the interactions between the atmosphere, the land, and the subsurface is fundamental to hydrology and is critical for a better assessment of the impacts of climate change and human management on hydrological systems. However, many land surface models simplify the subsurface hydrology and thereby these interactions. In this study, we couple the land surface model Noah-MP included in the NASA Land Information System (LIS) with the integrated hydrologic model ParFlow (ParFlow-LIS) using the Earth System Modeling Framework (ESMF) and the National United Operational Prediction Capability (NUOPC). This coupling improves the simulation of water and energy cycle processes by adding the three-dimensional variably saturated and heterogeneous flow in the subsurface using sophisticated and nonlinear physics-based equations as well as the advances in satellite remote sensing-based data assimilation of the land surface, thereby benefiting the integrated hydrologic modeling and data assimilation community. We use the High Plains aquifer, located in the central United States, as a testbed to evaluate the coupled ParFlow-LIS system. The new ParFlow-LIS system accounts for the effects of topographically driven flows on the land surface, producing more fine-scale patterns of land surface states and fluxes than standalone LIS. In addition, ParFlow-LIS enables the consideration of the effect of subsurface water storage on evapotranspiration. This is particularly important in areas and times with dry soils, such as during drought conditions or in the presence of a cone of depression due to pumping.
land surface modeling↗
Understanding Subsurface Fracture Evolution Dynamics Using Time‐Lapse Full Waveform Inversion of Continuous Active‐Source Seismic Monitoring Data
Abstract Predicting the behavior, geometry, and flow properties of subsurface fractures remains a challenging problem. Seismic models that can characterize fractures usually suffer from low spatiotemporal resolution. Here, we develop a correlative double‐difference time‐lapse full waveform inversion of continuous active source seismic monitoring data for determining high‐spatiotemporal‐resolution time‐lapse V p models of in‐situ fracture evolution at a shallow contamination site in Wyoming, USA. Assisted by rock physics modeling, we find that (a) rapidly increasing pore pressure initializes and grows the fracture, increasing the porosity slightly (from ∼13.7% to ∼14.6%) in the tight clay formation, thus decreasing V p (∼50 m/s); (b) the fluid injection continues decreasing V p , likely through the introduction of gas bubbles in the injectate; and (c) final V p reductions reach over ∼150 m/s due to a posited ∼4.5% gas saturation. Our results demonstrate that high‐resolution V p changes are indicative of mechanical and fluid changes within the fracture zone during hydrofracturing.
58 GEOSCIENCES↗
Time-Distance Helioseismology Data-Analysis Pipeline for Helioseismic and Magnetic Imager Onboard Solar Dynamics Observatory (SDO-HMI) and Its Initial Results
The Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) provides continuous full-disk observations of solar oscillations. We develop a data-analysis pipeline based on the time-distance helioseismology method to measure acoustic travel times using HMI Doppler-shift observations, and infer solar interior properties by inverting these measurements. The pipeline is used for routine production of near-real-time full-disk maps of subsurface wave-speed perturbations and horizontal flow velocities for depths ranging from 0 to 20 Mm, every eight hours. In addition, Carrington synoptic maps for the subsurface properties are made from these full-disk maps. The pipeline can also be used for selected target areas and time periods. We explain details of the pipeline organization and procedures, including processing of the HMI Doppler observations, measurements of the travel times, inversions, and constructions of the full-disk and synoptic maps. Some initial results from the pipeline, including full-disk flow maps, sunspot subsurface flow fields, and the interior rotation and meridional flow speeds, are presented.
Sun: helioseismology↗
Relating land surface, water table, and weathering fronts with a conceptual valve model for headwater catchments
Knowing little about how porosity and permeability are distributed at depth, we commonly develop models of groundwater by treating the subsurface as a homogeneous black box even though porosity and permeability vary with depth. One reason for this depth variation is that infiltrating meteoric water reacts with minerals to affect porosity in localized zones called reaction fronts. We are beginning to learn to map and model these fronts beneath headwater catchments and show how they are distributed. The subsurface landscapes defined by these fronts lie subparallel to the earth-air interface but with lower relief. They can be situated above, below, or at the water table. These subsurface landscapes of reaction are important because porosity developed from weathering can control subsurface water storage. In addition, porosity often changes at the weathering fronts, and when this affects permeability significantly, the front can act like a valve that re-orients water flowing through the subsurface. We explore controls on the positions of reaction fronts under headwater landscapes by accounting for the timescales of erosion, chemical equilibration, and solute transport. One strong control on the landscape of subsurface reaction is the land surface geometry, which is in turn a function of the erosion rate. In addition, the reaction fronts, like the water table, are strongly affected by the lithology and water infiltration rate. We hypothesize that relationships among the land surface, reaction fronts, and the water table are controlled by feedbacks that can push landscapes toward an “ideal hill”. In this steady state, reaction-front valves partition water volumes into shallow and deep flowpaths. Furthermore, these flows dissolve low- and high-solubility minerals, respectively, allowing their reaction fronts to advance at the erosion rate. This conceptualization could inform better models of subsurface porosity and permeability, replacing the black box.
54 ENVIRONMENTAL SCIENCES↗
Containment strategy for subsurface H 2 storage based on time-dependent soft solids
Here, we propose an innovative containment strategy based on time-dependent soft solids, namely 2 wt% Laponite suspensions, to reinforce natural subsurface seals and engineer flow barriers, with an eye toward making H 2 subsurface storage scalable and geographically agnostic. This suspension can be injected at its initial low viscosity and elasticity into a porous medium, allowing for easy pumping and targeted delivery. Once inside the target zone, it matures into a soft solid with much higher viscosity and elasticity, acting as a potential flow barrier. We discuss the rheological properties of the suspensions and demonstrate that hydrogen does not adversely affect their microstructure, but rather increases the suspensions’ viscosity and elasticity. Moreover, the suspensions enhance the rock samples’ compressive strength, while hydrogen exposure increases their stiffness and ductility. The ability of rock samples saturated with the suspensions to contain higher injected gas pressures is enhanced by aging at higher temperatures.
08 HYDROGEN↗
Tomographic Imaging of the Suns Interior
A new method is presented of determining the three-dimensional sound-speed structure and flow velocities in the solar convection zone by inversion of the acoustic travel-time data recently obtained by Duvall and coworkers. The initial inversion results reveal large-scale subsurface structures and flows related to the active regions, and are important for understanding the physics of solar activity and large-scale convection. The results provide evidence of a zonal structure below the surface in the low-latitude area of the magnetic activity. Strong converging downflows, up to 1.2 km/s, and a substantial excess of the sound speed are found beneath growing active regions. In a decaying active region, there is evidence for the lower than average sound speed and for upwelling of plasma.
Kosovichev, A. G.↗
Influence of Clay Wettability Alteration on Relative Permeability
Abstract Understanding the wettability of porous materials is important to model fluid flow in the subsurface. One of the critical factors that influences wetting in real reservoirs is the composition of geologic materials. The wetting properties for clay minerals can have a particularly strong impact on flow and transport. In this work, we analyze the chemical composition of a Mt. Simon sandstone core to resolve the microscopic structure of clay regions and assess how alterations to the local wetting properties influence multiphase transport based on core flooding experiments and relative permeability simulations. We show that whichever fluid has greater affinity toward clay minerals will tend to accumulate within these high surface area regions, leading to dramatic shifts in the relative permeability. This work establishes the essential importance of the mineral composition and associated wetting properties in the modeling of flow and transport in reservoir‐scale systems.
Fan, Ming↗
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↗
Flow Channeling in Fracture Networks: Characterizing the Effect of Density on Preferential Flow Path Formation
Flow channelization is a commonly observed phenomenon in fractured subsurface media where the flow of fluids is restricted primarily to highly transmissive fracture networks surrounded by a low-permeability rock matrix. The multiscale structural heterogeneity of these networks results in multiscale flow channelization where preferential flow paths form at length scales ranging from the entire system down to the subfracture size. We present an analysis of how one of the largest scales in fractured media, the network density, influences the degree of flow channeling that occurs using an ensemble of semigeneric three-dimensional discrete fracture network (DFN) simulations. We construct 10 DFNs, whose fracture lengths follow a power law distribution, at four densities for a total of 40 networks. We characterize their structure in terms of the network topology and geometry. Eulerian and Lagrangian observations of the steady-state flow fields obtained within the networks are used to quantify the degree of flow channelization at the network scale. We introduce a measure for the importance-ranking/hierarchy of different flow paths in the network using graph-based analysis of Lagrangian transport by which the degree of flow channeling between networks is compared. These flow observations are then linked to the structural properties of the networks. In general, network-scale flow channeling decreases as the network density increases. However, at low densities, there is more uniform flow within the entire connected network than in high-density networks. We also demonstrate how standard transport observables can be used to infer the degree of flow channelization occurring within a fracture network.
58 GEOSCIENCES↗
Inferring Gas Flow and Fracture Network Damage After a Subsurface Detonation Using Explosive and Geogenic Noble Gases.
Abstract not provided.
Bauer, Stephen↗
Emergence of Unstable Focused Flow Induced by Variable‐Density Flows in Vertical Fractures
Abstract Fluids with different densities often coexist in subsurface fractures and lead to variable‐density flows that control subsurface processes such as seawater intrusion, contaminant transport, and geologic carbon sequestration. In nature, fractures have dip angles relative to gravity, and density effects are maximized in vertical fractures. However, most studies on flow and transport through fractures are often limited to horizontal fractures. Here, we study the mixing and transport of variable‐density fluids in vertical fractures by combining three‐dimensional (3D) pore‐scale numerical simulations and visual laboratory experiments. Two miscible fluids with different densities are injected through two inlets at the bottom of a fracture and exit from an outlet at the top of the fracture. Laboratory experiments show the emergence of an unstable focused flow path, which we term a “runlet.” We successfully reproduce the unstable runlet using 3D numerical simulations and elucidate the underlying mechanisms triggering the runlet. Dimensionless number analysis shows that the runlet instability arises due to the Rayleigh‐Taylor instability (RTI), and flow topology analysis is applied to identify 3D vortices that are caused by the RTI. Even under laminar flow regimes, fluid inertia is shown to control the runlet instability by affecting the size and movement of vortices. Finally, we confirm the emergence of a runlet in rough‐walled fractures. Since a runlet dramatically affects fluid distribution, residence time, and mixing, the findings in this study have direct implications for the management of groundwater resources and subsurface applications.
Environmental Sciences & Ecology↗
Helioseismic Observations of the Structure and Dynamics of a Rotating Sunspot Beneath the Solar Surface
Time-distance helioseismology is applied to study the subphotospheric structures and dynamics of an unusually fast-rotating sunspot observed by the Michelson Doppler Imager on bead SOH0 in 2000 August. The subsurface sound speed structures and velocity fields are obtained for the sunspot region at different depths from 0 to 12 Mm. By comparing the subsurface sound speed variations with the surface magnetic field, we find evidence for structural twists beneath the visible surface of this active region, which may indicate that magnetic twists often seen at the photosphere also exist beneath the photosphere. We also report on the observation of subsurface horizontal vortical flows that extend to a depth of 5 Mm around this rotating sunspot and present evidence that opposite vortical flows may exist below 9 Mm. It is suggested that the vortical flows around this active region may build up a significant amount of magnetic helicity and energy to power solar eruptions. Monte Carlo simulation has been performed to estimate the error propagation, and in addition the sunspot umbra is masked to test the reliability of our inversion results. On the basis of the three-dimensional velocity fields obtained from the time-distance helioseismology inversions, we estimate the subsurface kinetic helicity at different depths for the first time and conclude that it is comparable to the current helicity estimated from vector magnetograms.
Zhao, Junwei↗
Impacts of Permeability Heterogeneity and Background Flow on Supercritical CO 2 Dissolution in the Deep Subsurface
Motivated by CO 2 capture and sequestration (CCS) design considerations, we consider the coupled effects of permeability heterogeneity and background flow on the dissolution of a supercritical CO 2 lens into an underlying deep, confined aquifer. We present the results of a large-scale Monte Carlo simulation study examining the interaction of background flow rate and three parameters describing multi-Gaussian log-permeability fields: mean, variance, and correlation length. Hundreds of high-resolution simulations were performed using the PFLOTRAN finite volume software to model CO 2 dissolution in a kilometer-scale aquifer over 1,000 years. Predictive dimensionless scaling relationships relating CO 2 dissolution rate to heterogeneity statistics, Rayleigh (Ra) and Péclet (Pe) numbers were developed for both gravitationally dominated free convection to background flow-dominated forced convection regimes. An empirical criterion, Pe = Ra 3/4 , was discovered for regime transition. All simulations converged quickly to a quasi-steady, approximately linear dissolution rate. However, this rate displayed profound variability between permeability field realizations sharing the same heterogeneity statistics, even under mild permeability heterogeneity. In general, increased heterogeneity was associated with a lower mean and higher variance of dissolution rate, undesirable from a CCS design perspective. The relationship between dissolution rate and background flow was found to be complex and nonlinear. Dimensionless scaling relationships were uncovered for a number of special cases. Results call into question the validity of the Boussinesq approximation in the context of modest-to-high background flow rates and the general applicability of numerical simulations without background flow.
42 ENGINEERING↗
Multiscale Feature-feature Interactions Control Patterns of Hyporheic Exchange in a Simulated Headwater Mountain Stream
Recent predictions of hyporheic exchange at the basin-scale assume individual features control exchange independently of each other, which has been demonstrated in relatively uniform, low-gradient rivers. However, this assumption may not hold in steep catchments where both the type and size of individual features may vary over short distances, leading to irregular patterns of feature dominance on hyporheic exchange flows. Also, steep longitudinal gradients support substantial downvalley flows in the subsurface, which may create feedbacks between adjacent features. In this study, we test the extent to which features interact with one another and whether they can be aggregated to make reach-scale predictions in a headwater mountain stream. Using systematic manipulations of a 2-D stream centerline model and spectral analyses, we test for the presence of both feature-feature and multiscale interactions. Our results show that changing the height of individual step-pool features can alter hyporheic flow fields in neighboring, and sometimes distant, features. Spectral analyses revealed two scales of streambed topography—a local scale of single features and an intermediate scale that spanned multiple local-scale features. All features produced hyporheic exchange, but turnover of deeper hyporheic water only occurred at a few key locations where local- and intermediate-scale features amplified each other. Further, shallow bedrock increases the ratio of local- and intermediate-scale flowpaths to regional-scale flowpaths. Conceptual models portraying hyporheic exchange as a series of nested flowpaths should include the interactions among streambed topographic features in mountain streams. These results have implications for upscaling, field experiments, and stream restoration in steep catchments.
54 ENVIRONMENTAL SCIENCES↗
Role of Intersections in Fracture Connectivity
Networks of intersecting fractures often provide the flow paths through subsurface reservoirs. Assessing network connectivity is challenging because fracture intersections compose a vanishingly small fraction of the network void volume. In this paper, motivated by 3D X-ray imaging of the simplest element of fracture network, that is, two orthogonal fractures, we perform a percolation and finite-size scaling analysis to study the connectivity provided by fracture intersections. The conditions when an intersection enhances connectivity across a sample depend on spatial correlations in the fracture aperture distributions, on the stress state, and on the direction of flow. Here we consider three flow directions: (a) across intersections, (b) parallel to intersections and (c) around corners. For (a), intersections provide minimal enhancement of connectivity because they contribute little additional void area. For (b), intersections increase the probability of a connected path near threshold by enabling 3D connected pathways that are not possible in parallel fractures. Flow around corners, (c), is fundamentally the result of the intersection connecting two fractures in series and spatial correlations are broken around corners, suppressing the connectivity relative to (a). When the connected fractures are stressed equally, a joint percolation threshold emerges that continues to have scale invariance. However, when the fractures are stressed unequally, the system has mixed percolation without clearly defined percolation thresholds. In all cases, percolation probabilities are found to be scale dependent which has important consequences for the connectivity of larger fracture networks composed of the fundamental element studied here.
02 PETROLEUM↗
Groundwater Inflows to the Columbia River Along the Hanford Reach and Associated Nitrate Concentrations
Healthy river ecosystems require the interaction of many physical and biological processes to maintain their status. One physical process supporting biogeochemical cycling is hydrologic exchange (i.e., hydrologic exchange flows, HEFs) between relatively fast-flowing channel waters and slower-flowing surface and subsurface waters (lateral and vertical). Land uses adjacent to rivers have the potential to alter the water quality of off-channel surface and subsurface waters, and HEFs therefore have the potential to deliver solutes associated with river-adjacent land uses to rivers. HEFs can be nonpoint, diffuse sources of pollution, making the ultimate pollution source difficult to identify, especially in large rivers. Here, we seek to identify HEFs in the Columbia River near Richland, WA by looking for anomalies in temperature and electrical conductivity (EC) along the bed of the river in February, June, July, August, and November 2018. These are ultimately the “ends” of HEFs as they are locations of subsurface inflow to the river. We found these anomalies to be a combination of warmer or colder and higher (but not lower) EC than river water. We identified a majority of warm anomalies in February and July 2018, and majority cold anomalies in June, August, and November 2018. High-EC anomalies were found mostly in February, August, and November. Combined, we observe a shift from warm, high EC anomalies dominating in February to equivalent EC, warm anomalies in June, to equivalent EC, cool anomalies dominating July. In August, we also measured dissolved nitrate ($NO^-_3$) in-situ to determine whether anomalies were associated with increased $NO^-_3$ loading to the river, especially along the eastern shoreline, which is dominated by agricultural land use. Inflows along the eastern shoreline have greater concentrations of nitrate than river water (up to 10 mg $N–NO^-_3$/L). This research demonstrates that HEFs are temporally and spatially dynamic transferring heat and solutes to rivers.
54 ENVIRONMENTAL SCIENCES↗
Simplified Method for the In Situ Collection and Laboratory Analysis of Cosmogenic Tracers (Sulfur-35 and Sodium-22) to Determine Residence Time Distributions and Water Ages
The use of cosmogenically produced sulfur-35 (T 1/2 = 87 days) and sodium-22 (T 1/2 = 2.6 years) as intrinsic tracers can provide valuable information on catchment hydrology, flow paths, and subsurface storage. A new and straightforward method was created to determine the activities of both 35 S and 22 Na in various water sources by pumping large volumes (up to 1000 L) of water through cation- and anion-exchange resin columns in the field to collect sodium and sulfate ions and simple chemistry in the lab. Samples are counted for 35 S using liquid scintillation counting (LSC) and for 22 Na via γ spectroscopy. Our novel in situ method provides faster sample throughput as well as better counting statistics and lower detection limits. Both methods were successfully applied at the Southern Sierra Critical Zone Observatory.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗