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 397 records · Page 22

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Imaging Complex Subsurface Structures for Geothermal Exploration at Pirouette Mountain and Eleven-Mile Canyon in Nevada

Accurate imaging of subsurface complex structures with faults is crucial for geothermal exploration because faults are generally the primary conduit of hydrothermal flow. It is very challenging to image geothermal exploration areas because of complex geologic structures with various faults and noisy surface seismic data with strong and coherent ground-roll noise. In addition, fracture zones and most geologic formations behave as anisotropic media for seismic-wave propagation. Properly suppressing ground-roll noise and accounting for subsurface anisotropic properties are essential for high-resolution imaging of subsurface structures and faults for geothermal exploration. We develop a novel wavenumber-adaptive bandpass filter to suppress the ground-roll noise without affecting useful seismic signals. This filter adaptively exploits both characteristics of the lower frequency and the smaller velocity of the ground-roll noise than those of the signals. Consequently, this filter can effectively differentiate the ground-roll noise from the signal. We use our novel filter to attenuate the ground-roll noise in seismic data along five survey lines acquired by the U.S. Navy Geothermal Program Office at Pirouette Mountain and Eleven-Mile Canyon in Nevada, United States. We then apply our novel anisotropic least-squares reverse-time migration algorithm to the resulting data for imaging subsurface structures at the Pirouette Mountain and Eleven-Mile Canyon geothermal exploration areas. The migration method employs an efficient implicit wavefield-separation scheme to reduce image artifacts and improve the image quality. Our results demonstrate that our wavenumber-adaptive bandpass filtering method successfully suppresses the strong and coherent ground-roll noise in the land seismic data, and our anisotropic least-squares reverse-time migration produces high-resolution subsurface images of Pirouette Mountain and Eleven-Mile Canyon, facilitating accurate fault interpretation for geothermal exploration.

15 GEOTHERMAL ENERGY↗

Organizational Principles of Hyporheic Exchange Flow and Biogeochemical Cycling in River Networks Across Scales

Hyporheic zones increase freshwater ecosystem resilience to hydrological extremes and global environmental change. However, current conceptualizations of hyporheic exchange, residence time distributions, and the associated biogeochemical cycling in streambed sediments do not always accurately explain the hydrological and biogeochemical complexity observed in streams and rivers. Specifically, existing conceptual models insufficiently represent the coupled transport and reactivity along groundwater and surface water flow paths, the role of autochthonous organic matter in streambed biogeochemical functioning, and the feedbacks between surface-subsurface ecological processes, both within and across spatial and temporal scales. While simplified approaches to these issues are justifiable and necessary for transferability, the exclusion of important hyporheic processes from our conceptualizations can lead to erroneous conclusions and inadequate understanding and management of interconnected surface water and groundwater environments. This is particularly true at the landscape scale, where the organizational principles of spatio-temporal dynamics of hyporheic exchange flow (HEF) and biogeochemical processes remain largely uncharacterized. This article seeks to identify the most important drivers and controls of HEF and biogeochemical cycling based on a comprehensive synthesis of findings from a wide range of river systems. We use these observations to test current paradigms and conceptual models, discussing the interactions of local-to-regional hydrological, geomorphological, and ecological controls of hyporheic zone functioning. This improved conceptualization of the landscape organizational principles of drivers of HEF and biogeochemical processes from reach to catchment scales will inform future river research directions and watershed management strategies.

54 ENVIRONMENTAL SCIENCES↗

Using passive, thermal remote sensing techniques for detecting subsurface gravel accumulations in vegetated, unconsolidated sedimentary terrains

Multiband radiometric data from an airborne imaging thermal scanner are being studied for use in finding buried gravel deposits. The techniques are based on measuring relative differences in the thermal properties between gravel-laden targets and the surrounding gravelless background. These properties are determined from modeling the spectral radiant emittance recorded over both types of surfaces in conjunction with ground measurements of the most significant heat flows above and below the surface. Thermodynamic properties of sampled materials from control sites are determined, and diurnal and annual subsurface heat waves are recorded. Thermal models that account for heat exchange at the surface, as well as varying levels of soil moisture, humidity, and vegetation, are needed for adaptation and modification to simulate the physical and radiative environments of this region.

Burns, Gregory S.↗

Development of a Compact, Deep-Penetrating Heat Flow Instrument for Lunar Landers: In-Situ Thermal Conductivity System

Geothermal heat flow is obtained as a product of the geothermal gradient and the thermal conductivity of the vertical soil/rock/regolith interval penetrated by the instrument. Heat flow measurements are a high priority for the geophysical network missions to the Moon recommended by the latest Decadal Survey and previously the International Lunar Network. One of the difficulties associated with lunar heat flow measurement on a robotic mission is that it requires excavation of a relatively deep (approx 3 m) hole in order to avoid the long-term temporal changes in lunar surface thermal environment affecting the subsurface temperature measurements. Such changes may be due to the 18.6-year-cylcle lunar precession, or may be initiated by presence of the lander itself. Therefore, a key science requirement for heat flow instruments for future lunar missions is to penetrate 3 m into the regolith and to measure both thermal gradient and thermal conductivity. Engineering requirements are that the instrument itself has minimal impact on the subsurface thermal regime and that it must be a low-mass and low-power system like any other science instrumentation on planetary landers. It would be very difficult to meet the engineering requirements, if the instrument utilizes a long (> 3 m) probe driven into the ground by a rotary or percussive drill. Here we report progress in our efforts to develop a new, compact lunar heat flow instrumentation that meets all of these science and engineering requirements.

Nagihara, S.↗

A detailed model for simulation of catchment scale subsurface hydrologic processes

A catchment scale numerical model is developed based on the three-dimensional transient Richards equation describing fluid flow in variably saturated porous media. The model is designed to take advantage of digital elevation data bases and of information extracted from these data bases by topographic analysis. The practical application of the model is demonstrated in simulations of a small subcatchment of the Konza Prairie reserve near Manhattan, Kansas. In a preliminary investigation of computational issues related to model resolution, we obtain satisfactory numerical results using large aspect ratios, suggesting that horizontal grid dimensions may not be unreasonably constrained by the typically much smaller vertical length scale of a catchment and by vertical discretization requirements. Additional tests are needed to examine the effects of numerical constraints and parameter heterogeneity in determining acceptable grid aspect ratios. In other simulations we attempt to match the observed streamflow response of the catchment, and we point out the small contribution of the streamflow component to the overall water balance of the catchment.

Paniconi, Claudio↗

Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications

This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.

42 ENGINEERING↗

Integration and Demonstration of Monitoring, Modeling, and Prediction of DV-1 Amendment Performance at the Bench Scale: DV-1 Amendment Demonstration

During fiscal years 2024 and 2025, the U.S. Department of Energy’s Hanford Field Office commissioned Pacific Northwest National Laboratory to conduct applied research aimed at reducing the cost, time, and uncertainty associated with in situ treatment of vadose zone contaminants at the Hanford Site. This report outlines the integration of three key research efforts into a meso-scale demonstration designed to advance field-scale solutions that aim to (1) optimize the delivery of chemical amendments to contaminated soils, (2) reduce uncertainty in amendment delivery performance assessment using advanced monitoring techniques, and (3) provide real-time insights into when and where amendment-induced precipitation reactions occur in the subsurface. To achieve these objectives, the tank-scale (~ 1 cubic meter) Geophysical Imaging of Flow and Transport (GIFT) system was developed. GIFT enables experimental testing of amendment delivery while incorporating automated multi-modal monitoring approaches, including pressure measurements, direct fluid sampling, and remote time-lapse geophysical imaging. The data generated from these monitoring techniques will serve as inputs for a generative artificial-intelligence-driven digital twin – a numerical simulation model designed to honor observed data while quantifying uncertainty in simulation accuracy. Using this simulator, researchers will refine an amendment injection strategy to maximize delivery efficiency within a low-permeability soil zone. Monitoring data will be interpreted through simulated outputs to enhance understanding of the injection process. The efficacy of this integrated approach will be evaluated through direct sampling at the conclusion of the experiment.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Versatile soil gas concentration and isotope monitoring: optimization and integration of novel soil gas probes with online trace gas detection

Abstract. Gas concentrations and isotopic signatures can unveil microbial metabolisms and their responses to environmental changes in soil. Currently, few methods measure in situ soil trace gases such as the products of nitrogen and carbon cycling or volatile organic compounds (VOCs) that constrain microbial biochemical processes like nitrification, methanogenesis, respiration, and microbial communication. Versatile trace gas sampling systems that integrate soil probes with sensitive trace gas analyzers could fill this gap with in situ soil gas measurements that resolve spatial (centimeters) and temporal (minutes) patterns. We developed a system that integrates new porous and hydrophobic sintered polytetrafluoroethylene (sPTFE) diffusive soil gas probes that non-disruptively collect soil gas samples with a transfer system to direct gas from multiple probes to one or more central gas analyzer(s) such as laser and mass spectrometers. Here, we demonstrate the feasibility and versatility of this automated multiprobe system for soil gas measurements of isotopic ratios of nitrous oxide (δ18O, δ15N, and the 15N site preference of N2O), methane, carbon dioxide (δ13C), and VOCs. First, we used an inert silica matrix to challenge probe measurements under controlled gas conditions. By changing and controlling system flow parameters, including the probe flow rate, we optimized recovery of representative soil gas samples while reducing sampling artifacts on subsurface concentrations. Second, we used this system to provide a real-time window into the impact of environmental manipulation of irrigation and soil redox conditions on in situ N2O and VOC concentrations. Moreover, to reveal the dynamics in the stable isotope ratios of N2O (i.e., 14N14N16O, 14N15N16O, 15N14N16O, and 14N14N18O), we developed a new high-precision laser spectrometer with a reduced sample volume demand. Our integrated system – a tunable infrared laser direct absorption spectrometry (TILDAS) in parallel with Vocus proton transfer reaction mass spectrometry (PTR-MS), in line with sPTFE soil gas probes – successfully quantified isotopic signatures for N2O, CO2, and VOCs in real time as responses to changes in the dry–wetting cycle and redox conditions. Broadening the collection of trace gases that can be monitored in the subsurface is critical for monitoring biogeochemical cycles, ecosystem health, and management practices at scales relevant to the soil system.

54 ENVIRONMENTAL SCIENCES↗

Assimilation of Altimeter Data into a Quasigeostrophic Model of the Gulf Stream System: Dynamical Considerations - Part 1

The dynamical consequences of constraining a numerical model with sea surface height data have been investigated. The model used for this study is a quasigeostrophic model of the Gulf Stream region. The data that have been assimilated are maps of sea surface height obtained as the superposition of sea surface height variability deduced from the Geosat altimeter measurements and a mean field constructed from historical hydrographic data. The method used for assimilating the data is the nudging technique. Nudging has been implemented in such a way as to achieve a high degree of convergence of the surface model fields toward the observations. The assimilation of the surface data is thus equivalent to the prescription of a surface pressure boundary condition. The authors analyzed the mechanisms of the model adjustment and the characteristics of the resultant equilibrium state when the surface data are assimilated. Since the surface data are the superposition of a mean component and an eddy component, in order to understand the relative role of these two components in determining the characteristics of the final equilibrium state, two different experiments have been considered: in the first experiment only the climatological mean field is assimilated, while in the second experiment the total surface streamfunction field (mean plus eddies) has been used. It is shown that the model behavior in the presence of the surface data constraint can be conveniently described in terms of baroclinic Fofonoff modes. The prescribed mean component of the surface data acts as a 'surface topography' in this problem. Its presence determines a distortion of the geostrophic contours in the subsurface layers, thus constraining the mean circulation in those layers. The intensity of the mean flow is determined by the inflow/outflow conditions at the open boundaries, as well as by eddy forcing and dissipation.

Capotondi, Antonietta↗

Electrochemical approaches for chemical and biological analysis on Mars

Obtaining in situ chemical data from planetary bodies such as Mars or Europa can present significant challenges. The one analytical technique that has many of the requisite characteristics to meet such a challenge is electroanalysis. Described here are three electroanalytical devices designed for in situ geochemical and biological analysis on Mars. The Mars Environmental Compatibility Assessment (MECA) was built and flight qualified for the now cancelled NASA Mars 2001 Lander. Part of MECA consisted of four "cells" containing arrays of electrochemical based sensors for measuring the ionic species in soil samples. A next-generation MECA, the Robotic Chemical Analysis Laboratory (RCAL), uses a carousel-type system to allow for greater customization of analytical procedures. A second instrument, proposed as part of the 2007 CryoScout mission, consists of a flow-through inorganic chemical analyzer (MICA). CryoScout is a torpedo-like device designed for subsurface investigation of the stratigraphic climate record embedded in Mars' north polar cap. As the CryoScout melts its way through the ice cap, MICA will collect and analyze the meltwater for a variety of inorganics and chemical parameters. By analyzing the chemistry locked in the layers of dust, salt, and ice, geologists will be able to determine the recent history of climate, water, and atmosphere on Mars and link it to the past. Finally, electroanalysis shows its abilities in the detection of possible microorganism on Mars or elsewhere in the solar system. To identify an unknown microorganism, one that may not even use Earth-type biochemistry, requires a detection scheme which makes minimal assumptions and looks for the most general features. Recent work has demonstrated that the use of an array of electrochemical sensors which monitors the changes in a solution via electrical conductivity, pH, and ion selective electrodes, can be used to detect minute chemical perturbations caused by the growth of bacteria and with the correct methodology provide unamibiguous detection of such life forms.

Non-NASA Center↗