Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

The Evolution of Oblique Impact Flow Fields Using Maxwell's Z Model

Oblique impacts are the norm rather than the exception for impact craters on planetary surfaces. This work focuses on the excavation of experimental oblique impact craters using the NASA Ames Vertical Gun Range (AVGR). Three-dimensional particle image velocimetry (3D PIV) is used to obtain quantitative data on ejection positions, three dimensional velocities and angles. These data are then used to constrain Maxwell's Z Model and follow the subsurface evolution of the excavation-stage flow-field center during oblique impacts.

Anderson, J. L. B.↗
DOE OSTI2022

Voronoi Meshing to Accurately Capture Geological Structure in Subsurface Simulations

Mesh generation lies at the interface of geological modeling and reservoir simulation. Highly skewed or very small grid cells may be necessary to accurately capture the geometry of geological features, but the resulting poorly scaled or small grid cells can have a substantial negative impact on simulator accuracy and speed. One way to minimize numerical errors caused by gridding complex structures is to simulate on high-quality Voronoi meshes, which reduce grid orientation effects in fluid flow. This work presents a complete methodology to create Voronoi simulation grids, model fluid flow in complex geological systems, and visualize the results. A recently developed Voronoi meshing method that can automatically generate provably good unstructured meshes that conform to input surfaces creating closed volumes is used. Initially an analytical benchmark simulation is presented to validate the quality of the meshes and simulation results and demonstrate the superiority of simulation results using Voronoi meshes over flexed-hexahedral meshes on a domain with internal features. Next, meshes are created for test structures representing four of the most common geological features in the subsurface: layering, pinch-out, an interior lens that tapers to zero thickness on all sides and a fault with offset. Two benchmark flow simulations are run for each test structure. Finally, a realistic geological example for CO 2 injection into an anticline is simulated. Three realizations of the Voronoi mesh at the same resolution are generated for the simulations. Each mesh is highly refined near the injection wells and coarse in areas of less interest. These three meshes are used to model the CO 2 plume in the subsurface as it migrates to the top of the structure and then fills downward. Simulations on the meshes with randomly generated elements inside the input volumes each give slightly different fingering patterns for the viscous-unstable buoyant gas flow. The results presented in this work show a promising step towards utilizing fully automated Voronoi meshing for subsurface flow simulations in complex geology.

58 GEOSCIENCES↗
DOE OSTI2022

Confinement-Driven Heterogeneous Benzene Crystallization in Silica Nanopores

Nanoconfinement alters the thermodynamics, dynamics, and kinetics of fluids hosted in nanoscale solid nanopores to an extent that depends on the characteristics of the confining space and the chemistry of the confined fluids. Confinement-induced alterations in the phase behavior of confined energetic fluids under high pressure or low temperature are highly relevant to subsurface and subsea phenomena such as fluid flow in porous media, hydrate formation and dissociation, and gas storage capacity. Although extensive efforts have been directed toward understanding the phase behavior of confined fluids, the role of solid-liquid interfaces in the phase transitions of organic liquids has not been resolved yet. Here, we explore the onset and growth of benzene crystallization confined in 6 nm sized SBA-15 silica nanopores in the temperature range 300-200 K using in situ extended range small-angle and wide-angle X-ray scattering (SAXS/WAXS) measurements and atomistic classical molecular dynamics (MD) simulations. The crystallization onset of confined benzene depresses to 265 K compared to the freezing point of bulk benzene ( ~278 K), followed by the continuous growth of the emerged crystals in the pore space with complete crystallization at 200 K. The orientation of the emerged benzene crystals is dominated by parallel (π-π stacking) and perpendicular (T-shape stacking) orientations along the cylindrical pore radius and pore length, respectively. The onset of benzene crystals occurs heterogeneously on the pore surface and grows continuously toward the pore center. Further, confined benzene undergoes a dynamical crossover from fragile to strong dynamics behavior, inferred from the rotational and translational diffusion. The insights provided by this study have significant implications for the phase transitions of confined organic liquids that are relevant to a wide range of applications in the biological, geological, environmental, and chemical fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

NRAP-Open-IAM Multisegmented Wellbore Reduced-Order Model: Improvement and Quality Assurance

The multisegmented wellbore model (MSW) semi-analytically estimates the amount of CO 2 and brine leakage from a leaking legacy well by segmenting it into intervals to simulate site-specific stratigraphic and hydrogeologic properties. The model is a component of the National Risk Assessment Partnership Open-Source Integrated Assessment Model (NRAP-Open-IAM), which was developed to perform risk assessment for geologic CO 2 storage. The new wellbore leakage model, which uses deep learning networks for a caprock segment, was developed to enhance the analytical MSW. The model was trained and validated using a synthetic data set of Subsurface Transport Over Multiple Phases (STOMP) multiphase flow simulations from various geological, well attribute, and operational conditions to ensure its quality. The results demonstrate that the model is more accurate than the existing model in predicting the transport of two-phase fluids (brine and injected CO 2 ) through the well. This report provides a detailed explanation of the model development and quality assurance.

58 GEOSCIENCES↗
DOE OSTI2023

Knowledge-informed deep learning for hydrological model calibration: an application to Coal Creek Watershed in Colorado

Abstract. Deep learning (DL)-assisted inverse mapping has shown promise in hydrological model calibration by directly estimating parameters from observations. However, the increasing computational demand for running the state-of-the-art hydrological model limits sufficient ensemble runs for its calibration. In this work, we present a novel knowledge-informed deep learning method that can efficiently conduct the calibration using a few hundred realizations. The method involves two steps. First, we determine decisive model parameters from a complete parameter set based on the mutual information (MI) between model responses and each parameter computed by a limited number of realizations (∼50). Second, we perform more ensemble runs (e.g., several hundred) to generate the training sets for the inverse mapping, which selects informative model responses for estimating each parameter using MI-based parameter sensitivity. We applied this new DL-based method to calibrate a process-based integrated hydrological model, the Advanced Terrestrial Simulator (ATS), at Coal Creek Watershed, CO. The calibration is performed against observed stream discharge (Q) and remotely sensed evapotranspiration (ET) from the water year 2017 to 2019. Preliminary MI analysis on 50 realizations resulted in a down-selection of 7 out of 14 ATS model parameters. Then, we performed a complete MI analysis on 396 realizations and constructed the inverse mapping from informative responses to each of the selected parameters using a deep neural network. Compared with calibration using observations covering all time steps, the new inverse mapping improves parameter estimations, thus enhancing the performance of ATS forward model runs. The Nash–Sutcliffe efficiency (NSE) of streamflow predictions increases from 0.53 to 0.8 when calibrating against Q alone. Using ET observations, on the other hand, does not show much improvement on the performance of ATS modeling mainly due to both the uncertainty of the remotely sensed product and the insufficient coverage of the model ET ensemble in capturing the observation. By using observed Q only, we further performed a multiyear analysis and show that Q is best simulated (NSE > 0.8) by including in the calibration the dry-year flow dynamics that show more sensitivity to subsurface characteristics than the other wet years. Moreover, when continuing the forward runs till the end of 2021, the calibrated models show similar simulation performances during this evaluation period as the calibration period, demonstrating the ability of the estimated parameters in capturing climate sensitivity. Our success highlights the importance of leveraging data-driven knowledge in DL-assisted hydrological model calibration.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1971

Observations of actively forming lava tubes and associated structures, Hawaii, part 2

A ground examination is made of lave tubes and channels. The surface morphology and the changes noted through lava flow activity are cited, and compared to earlier aerial observations. The lava activity was believed to be caused by a small lava lake exposed by the collapse of a crust covering it. Drainage of the lake was caused by a fissure erruption. New tubes or extensions of existing ones were noted from the flow. Molten lava was not seen in any tubes examined on the ground, but some of the flows were not sufficiently cooled to allow subsurface examination and survey of the tubes.

Greeley, R.↗
NASA NTRS1982

Morphology and network patterns of Martian valleys

The erosional origin of Martian valleys is briefly discussed. The absence of streamlined obstacles and interior channels, the steep and rugged walls and flat floors, and the presence of steep-walled, cuspate amphitheater terminations and sparse, open networks, argue against erosion by rainfall and surface flow and in favor of erosion by subsurface water.

Pieri, D. C.↗
NASA NTRS2003

The Evolution of Oblique Impact Flow Fields Using Maxwell's Z Model

Oblique impacts are the norm rather than the exception for impact craters on planetary surfaces. This work focuses on the excavation of experimental oblique impact craters using the NASA Ames Vertical Gun Range (AVGR). Three-dimensional particle image velocimetry (3D PIV) is used to obtain quantitative data on ejection positions, three-dimensional velocities and angles. These data are then used to test the applicability and limitations of Maxwell's Z Model in representing the subsurface evolution of the excavation-stage flow-field center during vertical and oblique impacts.

Anderson, J. L. B.↗
NASA NTRS2004

Radio Sounding Techniques for the Galilean Icy Moons and their Jovian Magnetospheric Environment

Radio sounding of the Earth's topside ionosphere and magnetosphere is a proven technique from geospace missions such as the International Satellites for Ionospheric Studies (ISIS) and the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE). Application of this technique to Jupiter's icy moons and the surrounding Jovian magnetosphere will provide unique remote sensing observations of the plasma and magnetic field environments and the subsurface conductivities, of Europa, Ganymede, and Callisto. Spatial structures of ionospheric plasma above the surfaces of the moons vary in response to magnetic-field perturbations from (1) magnetospheric plasma flows, (2) ionospheric currents from ionization of sputtered surface material, and (3) induced electric currents in salty subsurface oceans and from the plasma flows and ionospheric currents themselves. Radio sounding from 3 kHz to 10 MHz can provide the global electron densities necessary for the extraction of the oceanic current signals and supplements in-situ plasma and magnetic field measurements. While radio sounding requires high transmitter power for subsurface sounding, little power is needed to probe the electron density and magnetic field intensity near the spacecraft. For subsurface sounding, reflections occur at changes in the dielectric index, e.g., at the interfaces between two different phases of water or between water and soil. Variations in sub-surface conductivity of the icy moons can be investigated by radio sounding in the frequency range from 10 MHz to 50 MHz, allowing the determination of the presence of density and solid-liquid phase boundaries associated with oceans and related structures in overlying ice crusts. The detection of subsurface oceans underneath the icy crusts of the Jovian moons is one of the primary objectives of the Jupiter Icy Moons Orbiter (JIMO) mission. Preliminary modeling results show that return signals are clearly distinguishable be&een an ice crust with a thickness of 7 km on 1) an ocean and 2) a layer of bedrock. Knowledge of the ionospheric contributions to the time delay of the low-frequency subsurface radar is shown to be important in obtaining accurate depth information.

Green, James L.↗
DOE OSTI2024

An integrated experimental–modeling approach to identify key processes for carbon mineralization in fractured mafic and ultramafic rocks

Abstract Controlling atmospheric warming requires immediate reduction of carbon dioxide (CO2) emissions, as well as the active removal and sequestration of CO2 from current point sources. One promising proposed strategy to reduce atmospheric CO2 levels is geologic carbon sequestration (GCS), where CO2 is injected into the subsurface and reacts with the formation to precipitate carbonate minerals. Rapid mineralization has recently been reported for field tests in mafic and ultramafic rocks. However, unlike saline aquifers and depleted oil and gas reservoirs historically considered for GCS, these formations can have extremely low porosities and permeabilities, limiting storage volumes and reactive mineral surfaces to the preexisting fracture network. As a result, coupling between geochemical interactions and the fracture network evolution is a critical component of long-term, sustainable carbon storage. In this paper, we summarize recent advances in integrating experimental and modeling approaches to determine the first-order processes for carbon mineralization in a fractured mafic/ultramafic rock system. We observe the critical role of fracture aperture, flow, and surface characteristics in controlling the quantity, identity, and morphology of secondary precipitates and present where the influence of these factors can be reflected in newly developed thermo-hydro-mechanical–chemical models. Our findings provide a roadmap for future work on carbon mineralization, as we present the most important system components and key challenges that we are overcoming to enable GCS in mafic and ultramafic rocks.

58 GEOSCIENCES↗
DOE OSTI2021

Understanding Discrete Fracture Networks Through Spectral Graph Theory

Discrete Fracture Network models (DFNs) are used to simulate fluid flow and particle transport through fracture networks in low permeability rock. Understanding these processes are essential in many subsurface applications, such as environmental restoration of contaminated fractured media, CO 2 sequestration, detection of low-level nuclear tests, and hydrocarbon extraction. Compared with other models, DFNs allow for incorporation of a wider range of network characteristics but have substantially greater computation cost. These networks can be represented with graphs, allowing the use of graph theory tools to study the networks. I used Python to simulate flow and transport on a range of DFNs and analyzed these networks using methods from network analysis and spectral graph theory. My purpose was to find ways to gain insight about flow and transport on DFNs using these graph representations, bypassing the computationally intensive meshing typically required. My work is still in progress, but I have discovered several interesting trends and patterns that I believe could be useful towards my goal. If I am able to bring these results to fruition, they will aid subsurface geologists in extracting flow and transport information about fracture networks more efficiently.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2024

Emplacement of Shergottites in the Martian Crust Inferred From 3D Petrofabric and Crystal Size Distribution Analyses

Shergottites are mafic to ultramafic igneous rocks that represent a majority of known martian meteorites. They are subdivided into gabbroic, poikilitic, basaltic, and olivine-phyric categories based on differences in mineralogy and textures. Their geologic contexts are unknown so analyses of crystal sizes and preferred orientations have commonly been used to infer where shergottites solidified. Such environments range from subsurface cumulates to shallow intrusives to extrusive lava flows, which all have contrasting implications for interactions with crustal material, cooling histories, and potential in situ exposure at the surface. In this study, we present a novel three-dimensional (3D) approach to better understand the solidification environments of these samples and improve our knowledge of shergottites’ geologic contexts. Shape preferred orientations of most phases and crystal size distributions of late-forming minerals were measured in 3D using X-ray computed tomography (CT) on eight shergottites representing the gabbroic, poikilitic, basaltic, and olivine-phyric categories. Our analyses show that highly anisotropic, rod-like pyroxene crystals are strongly foliated in the gabbroic samples but have a weaker foliation and a mild lineation in the basaltic sample, indicating a directional flow component in the latter. Star volume distribution analyses revealed that most phases (maskelynite, pyroxene, olivine, and oxides/sulfides) preserve a foliated texture with variable strengths, and that the phases within individual samples are strongly to moderately aligned with respect to one another. In combination with relative cooling rates during the final stages of crystallization determined from interstitial oxide/sulfide crystal size distribution analyses, these results indicate that the olivine-phyric samples were emplaced as shallow intrusives (e.g., dikes/sills) and that the gabbroic, poikilitic, and basaltic samples were emplaced in deeper subsurface environments.

Martian meteorites↗
DOE OSTI2021

Predicting Thermal Performance of an Enhanced Geothermal System From Tracer Tests in a Data Assimilation Framework

Abstract Predicting the thermal performance of an enhanced geothermal system (EGS) requires a comprehensive characterization of the underlying fracture flow patterns from practically available data such as tracer data. However, due to the inherent complexities of subsurface fractures and the generally insufficient geological/geophysical data, interpreting tracer data for fracture flow characterization and thermal prediction remains a challenging task. The present study aims to tackle the challenge by leveraging a data assimilation method to maximize the utilization of information inherently contained in tracer data, and meanwhile maintain the flexibility to handle various uncertainties. A tracer data interpretation framework was proposed with the following three components integrated: (a) We use principal component analysis (PCA) to reduce the dimensionality of model parameter space. (b) We use ES‐MDA (ensemble smoother with multiple data assimilation) to invert for fracture aperture/flow fields and obtain posterior model ensembles for uncertainty quantification. Various data types are assimilated jointly to improve the predictive ability of the posterior ensemble. (c) The inverted fracture aperture fields are then incorporated into reservoir models to predict thermal performance. We developed a field‐scale EGS model to verify the ability of the framework to characterize highly heterogeneous fracture aperture/flow fields and predicting thermal performance. We also applied the framework to a mesoscale field experiment to demonstrate its potential application in real‐world geothermal reservoirs. The results indicate that the proposed framework can effectively retrieve fracture flow information from tracer data for thermal prediction and uncertainty quantification, and thus provide informative guidance for EGS optimization and risk management.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Strong temporal variability in methane fluxes from natural gas well pad soils

We measured methane and carbon dioxide fluxes at natural gas well pad soils and undisturbed soils in the Rocky Mountain and Gulf Coast regions of the United States, including producing and gas storage wells. We collected both short-term (15 min) and multi-day (between 3 and 8), continuous measurements at 47 well pads and two undisturbed locations. Methane fluxes varied by more than an order of magnitude over periods as short as 30 min (e.g., 19–593 mg m -2 h -1 in one instance), and diurnal and seasonal variability was also significant (e.g., spring-to-fall change from 509 to 14174 mg m -2 h -1 ). We hypothesize that short-term flux variability was caused by pulsed flow of methane during its migration through the subsurface. Barometric pressure and well conditions likely impacted fluxes, but we found only weak evidence for this. Bacterial methanotrophy appeared to impact methane flux magnitude and variability. We injected methane into the subsurface at one well, and we found that, while fluxes of methane and carbon dioxide, and combustible soil gas concentrations, increased in response to the injection, the response was not uniform, and fluxes exhibited high hourly-scale variability, in spite of a constant injection rate. Methane fluxes tended to be higher at well pad soils compared to background soils (often much higher), and fluxes tended to be higher at well pad locations closer to the well head.

03 NATURAL GAS↗
DOE OSTI2020

Operational and geological controls of coupled poroelastic stressing and pore-pressure accumulation along faults: Induced earthquakes in Pohang, South Korea

Coupled poroelastic stressing and pore-pressure accumulation along pre-existing faults in deep basement contribute to recent occurrence of seismic events at subsurface energy exploration sites. Our coupled fluid-flow and geomechanical model describes the physical processes inducing seismicity corresponding to the sequential stimulation operations in Pohang, South Korea. Simulation results show that prolonged accumulation of poroelastic energy and pore pressure along a fault can nucleate seismic events larger than M w 3 even after terminating well operations. In particular the possibility of large seismic events can be increased by multiple-well operations with alternate injection and extraction that can enhance the degree of pore-pressure diffusion and subsequent stress transfer through a rigid and low-permeability rock to the fault. This study demonstrates that the proper mechanistic model and optimal well operations need to be accounted for to mitigate unexpected seismic hazards in the presence of the site-specific uncertainty such as hidden/undetected faults and stress regime.

58 GEOSCIENCES↗
DOE OSTI2024

Model data for Flood Frequency Analysis using Stochastic Storm Transposition and an Integrated Surface-Subsurface Hydrological Model

This archived provides scripts and input files used for the implementation of a novel approach to conduct process-based Flood Frequency Analysis using a Stochastic Storm Transposition (SST) and an Integrated Surface-Subsurface Hydrological Model (ISSHM). As a proof-of-concept, this study uses the ISSHM, Advanced Terrestrial Simulator (Amanzi-ATS) model, and the SST model, RainyDay, to conduct flood frequency analysis by simulating the flood response to 5,000 annual synthetic storm events in a ~2000 km2 Southeast Texas watershed.The Watershed Workflow package is implemented in Python3. The Jupyter notebooks can be executed through multiple open-source tools, for example, Anaconda Jupyter Lab, VS Studio Code, etc. Other data files include TXT, CSV, DAT, SBATCH, SHP, TIF, NetCDF, and HDF5 files, which can be read through Python scripts. The input files for the ATS model and RainyDay model have .XML and .SST extensions, respectively, and can be edited in any commonly used text editors.This archive contains:* Scripts and data files essential for generating the ATS model input. It uses the Watershed Workflow package to produce both mesh and ATS input files. * Jupyter notebooks designated for the ATS model evaluation, covering both long-term simulations and 40 rainfall-runoff events.* Input files required to simulate SST storm events using RainyDay.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1986

Experimental forecasts of El Nino

A deterministic numerical model of the coupled evolution of the tropical ocean and atmosphere was used to forecast all El Nino/Southern Oscillation (ENSO) events from 1970 to 1986. More particularly, the model, originally developed for studying large-scale ocean-atmosphere interactions in the tropics, successfully predicted the characteristics of the spatial and temporal structure of ENSO observed in the study interval. The model indicated that rainfall moving eastward over the Pacific slackens the surface winds that would otherwise cool the eastern Pacific by drawing up cooler subsurface waters. The oceanic thermocline increases, a poleward flow of westerly flowing warn waters deplets the equatorial warm water reservoir, and sea surface temperatures decline. These ENSO conditions are statistically tractable with the model several months in advance, provided upper ocean layer thermal data are available.

Cane, M. A.↗
NASA NTRS2010

Time-Distance Helioseismology

Time-distance helioseismology is a method of ambient noise imaging using the solar oscillations. The basic realization that led to time-distance helioseismology was that the temporal cross correlation of the signals at two 'surface' (or photospheric) locations should show a feature at the time lag corresponding to the subsurface travel time between the locations. The temporal cross correlation, as a function of the location separation, is the Fourier transform of the spatio-temporal power spectrum of the solar oscillations, a commonly used function in helioseismology. It is therefore likely the characteristic ridge structure of the correlation function had been seen before without appreciation of its significance. Travel times are measured from the cross correlations. The times are sensitive to a number of important subsurface solar phenomena. These include sound speed variations, flows, and magnetic fields. There has been much interesting progress in the 17 years since the first paper on this subject (Duvall et al., Nature, 1993, 362, 430-432). This progress will be reviewed in this paper.

Duvall, Thomas L., Jr.↗