Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid characterization”

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
DOE OSTI2025

Using principal component analysis to distinguish different dynamic phases in superconducting vortex matter

Vortices in type-II superconductors driven over random disorder are known to exhibit a remarkable variety of distinct nonequilibrium dynamical phases that arise owing to the competition between vortex-vortex interactions, the quenched disorder, and the drive. These include pinned states, elastic flows, plastic or disordered flows, and dynamically reordered moving crystal or moving smectic states. The plastic flow phases can be particularly difficult to characterize since the flows are strongly disordered. Here, we perform principal component analysis (PCA) on the positions and velocities of vortex matter moving over random disorder for different disorder strengths and drives. We find that PCA can distinguish the known dynamic phases as well as or better than previous measures based on transport signatures or topological defect densities. In addition, PCA recognizes distinct plastic flow regimes, a slowly changing channel flow and a moving amorphous fluid flow, that do not produce distinct signatures in the standard measurements. In conclusion, our results suggest that this position and velocity-based PCA approach could be used to characterize dynamic phases in a broader class of systems that exhibit depinning and nonequilibrium phase transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Final Report DE-FE0031785 Mohsen Ahmadian, Ph.D. Demonstration of Proof of Concept of a Multiphysics Approach for Real-Time Remote Monitoring of Dynamic Changes in Pressure and Salinity in Hydraulically Fractured Networks

Hydraulic fracturing has evolved into a multistep process with varying flow rates, carrier fluids (e.g., gel or slickwater), proppant loadings, and proppant grain sizes. As a result, primary recovery from a hydraulically fractured tight-oil reservoir is often a tiny fraction of the original oil in place, ranging between 5 and 10%. As stated in the FOA1990, “part of this problem is due to the inability of current well completion processes to effectively stimulate the entire reservoir volume in contact with the wellbore. Innovative technologies are needed that can help improve the effectiveness of reservoir completion methods, maximize stimulated reservoir volumes, and optimize recovery over the entire producing life span of a well”. We first need to enhance the current fracture diagnostic techniques to improve a well-completion design. However, detecting and delineating a subsurface hydraulic fracture is extremely difficult because the induced fracture network is only fractionally propped, and these propped fractures are generally very thin. Microseismic and tiltmeter monitoring techniques can provide information on the fracture extent but provide little or no information on the movement and final distribution of proppant or production fluids. On the other hand, electromagnetic (EM) imaging has shown the capability to monitor proppant distribution throughout the fracture area, especially in the presence of Electrically Active Proppants (EAPs). A previous EM survey of hydraulic fracturing at the Devine Fracture Pilot Site (DFPS) and subsequent EM code developments demonstrated this survey as a robust technique to remotely interrogate the extent of the EAP-filled hydraulic fracture during its propagation. The objectives of the project were threefold: (1) to capitalize on the material properties of an EAP to demonstrate remote monitoring of relative changes in pressure, pressure, and flow that are commonly encountered during production from a hydraulically fractured reservoir; (2) to evaluate EM imaging tools, to achieve Objective 1 in near real-time; and (3) to develop a multi-physics joint inversion approach to precisely predict flow patterns and physiochemical changes within an EAP-filled fracture network. This research project was built upon our previous work at the Devine Test Site managed by the Bureau of Economic Geology (BEG) at The University of Texas at Austin (UT-Austin). It also leveraged a significant investment from the Advanced Energy Consortium (AEC) to address the DOE's interest in subsurface flow, containment, and characterization by multiple signals. This three-year and three-month project succeeded in demonstrating the feasibility of a real-time dynamic fluid flow mapping technique at Technology Readiness Level 5 (TRL5) by utilizing a commercially available surface-based Controlled-Source Electromagnetic (CSEM) method (Objectives 1, 2). We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, the observed electric field in our study is influenced by bottomhole pressure, flow rate, and salinity, which is demonstrated by comparing these parameters with the electrical field potentials. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce the whole measured electric field magnitudes. Preliminary estimates showed that including Streaming Potential (SP) in our geophysical model is likely needed to reduce the simulation misfit.

02 PETROLEUM↗
DOE OSTI2020

Tuning the quantumness of simple Bose systems: A universal phase diagram

We present a comprehensive theoretical study of the phase diagram of a system of many Bose particles interacting with a two-body central potential of the so-called Lennard-Jones form. First-principles path-integral computations are carried out, providing essentially exact numerical results on the thermodynamic properties. The theoretical model used here provides a realistic and remarkably general framework for describing simple Bose systems ranging from crystals to normal fluids to superfluids and gases. The interplay between particle interactions on the one hand and quantum indistinguishability and delocalization on the other hand is characterized by a single quantumness parameter, which can be tuned to engineer and explore different regimes. Taking advantage of the rare combination of the versatility of the many-body Hamiltonian and the possibility for exact computations, we systematically investigate the phases of the systems as a function of pressure (P) and temperature (T), as well as the quantumness parameter. Here, we show how the topology of the phase diagram evolves from the known case of 4 He, as the system is made more (and less) quantum, and compare our predictions with available results from mean-field theory. Possible realization and observation of the phases and physical regimes predicted here are discussed in various experimental systems, including hypothetical muonic matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

One-dimensional, transient modeling of a fixed-bed regenerator as a replacement for recuperators in supercritical CO 2 power cycles

Fixed-bed regenerators may be a low cost alternative to recuperators for supercritical CO 2 power cycles. However, the transient nature of regenerators operating in sCO 2 power cycles has not been previously characterized. This study details a one-dimensional, transient model of a fixed-bed regenerator that captures phenomena such as local fluid properties and heat transfer coefficients, entrained fluid heat capacity, and the effects of conservation of mass during process switching. The model is validated with experimental results of a 10 kWth regenerator. This comparison shows that the model predicts temperature as a function of both time and space exceptionally well. The model predicts effectiveness, heat transfer, and discharge temperatures within 5%, and predicts pressure drop and fluid carryover within 10%. This model could be used in future studies to simulate regenerators within a supercritical CO 2 power cycle.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Simulation of the supercritical CO 2 recompression Brayton power cycle with a high-temperature regenerator

The supercritical carbon dioxide (sCO 2 ) recompression Brayton cycle promises higher efficiency and lower capital cost than traditional steam Rankine power cycles. However, achieving high efficiency requires large, highly effective recuperators. Regenerators may be a low-cost alternative to printed circuit and micro-tube heat exchangers for recuperation in sCO 2 power cycles. Regenerators are a periodic heat exchanger in which thermal energy is extracted from the hot stream, stored in solid media, and then released to the cold stream at a later time. Fixed bed regenerators with valves to direct fluid are the preferred method for implementing regenerators in power cycles, but the inherently transient nature of these systems has not been characterized for this application. This study presents the simulation of a high-temperature regenerator within a 10 MWe sCO 2 recompression Brayton cycle. A transient, one-dimensional regenerator model presented in a previous study is used to simulate the regenerator. Dynamic heat exchanger models are also developed for the precooler, low-temperature recuperator, and primary heat exchanger, and the compressors and the turbine are modeled with off-design performance maps. We assess two regenerator-valve design options; one for fast switching, and one for reduced flow rate fluctuations. System simulation finds that both designs see significant fluctuations in turbomachinery flow rate, turbomachinery and system power, and regenerator discharge process outlet temperature. While designing the regenerator-valve subsystem for lower fluctuations is possible, the regenerator cold discharge temperature and net power still fluctuate by ±77.6°C and 6%, respectively. Increasing buffer volume is not effective at sufficiently reducing these fluctuations, but adding a packed bed in between the regenerator and the primary heat exchanger can reduce regenerator discharge process outlet temperature fluctuations to 6.4°C. Further reductions could be possible by increasing the size of this packed bed.

42 ENGINEERING↗
DOE OSTI2021

Modeling Nanosecond-Pulsed Spark Discharge and Flame Kernel Evolution

Dilute combustion, either using exhaust gas recirculation or with excess air, is considered a promising strategy to improve the thermal efficiency of internal combustion engines. However, the dilute air-fuel mixture, especially under intensified turbulence and high-pressure conditions, poses significant challenges for ignitability and combustion stability, which may limit the attainable efficiency benefits. In-depth knowledge of the flame kernel evolution to stabilize ignition and combustion in a challenging environment is crucial for effective engine development and optimization. To date, a comprehensive understanding of ignition processes that result in the development of fully predictive ignition models usable by the automotive industry does not yet exist. Spark-ignition consists of a wide range of physics that includes electrical discharge, plasma evolution, joule-heating of gas, and flame kernel initiation and growth into a self-sustainable flame. In this study, an advanced approach is proposed to model spark-ignition energy deposition and flame kernel growth. To decouple the flame kernel growth from the electrical discharge, a nanosecond-pulsed high-voltage discharge is used to trigger spark-ignition in an optically accessible small ignition test vessel with a quiescent mixture of air and methane. Initial conditions for the flame kernel, including its thermodynamic state and species composition, are derived from a plasma-chemical equilibrium calculation. The geometric shape and dimension of the kernel are characterized using a multi-dimensional thermal plasma solver. Here, the proposed modeling approach is evaluated using a high-fidelity computational fluid dynamics procedure to compare the simulated flame kernel evolution against flame boundaries from companion Schlieren images.

42 ENGINEERING↗
DOE OSTI2020

Application of an Embedded Fracture and Borehole Modeling Approach to the Understanding of EGS Collab Experiment 1

EGS Collab is a series of meso-scale experiments and associated numerical simulation activities being funded by the United States Department of Energy, Geothermal Technologies Office (GTO) to investigate enhanced geothermal system processes under in-situ stress and slightly elevated temperature conditions. This project is designed to provide scientists and engineers with immediate access to impermeable rock at scales larger than possible in the laboratory, but generally smaller than those for commercial production. Immediate access to rock is provided via the existing drifts of the former Homestake Gold Mine, now operated as the Sanford Underground Research Facility in Lead, South Dakota. The objectives of the EGS Collab project are to develop well controlled fracture networks between injector and producer boreholes using normal and shear stimulation for permeability enhancement. The first experimental site was located off the West Access drift on the 4850 Level (4850 feet below ground surface) in phyllite of the Precambrian Poorman formation, and involved the creation of a fracture network comprising a combination of hydraulic and natural fractures. A second experimental site is now being considered near the battery alcove on the 4100 Level in amphibolite and rhyolite of the Yates Unit. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO’s flagship EGS research effort. Experiment 1 has comprised a series of successful tests, including a long-term chilled-water circulation test, but the testbed has two atypical EGS elements. Active ventilation in the adjacent drift over a 50-year period cooled the testbed rock mass and hydraulic fracturing intersected monitoring boreholes making them conduits for fluid flow. Numerical simulations executed in support of the design of the EGS Collab Experiment 1, computed a radial temperature distribution and stress gradient orthogonal to the drift axis, resulting in the forecast of an oblong hydraulic fracture geometry, extended in the direction of the drift from the stimulation borehole. This paper describes the numerical simulation of long-term chilled-water circulation test with an embedded fracture and borehole modeling approach. The principal objective of the simulation work is show agreement between experimental observations in terms of production fluid temperatures, monitoring borehole temperatures, injection and production pressures, and tracer recoveries using a single conceptual model for the fracture network, using characterization and monitoring data generally available for EGS. A secondary objective is to improve the numerical simulation result comparisons with additional data available from the broader monitoring equipment within the testbed versus deeper, hotter, and more remote EGS. One complicating factor for comparisons between the numerical simulation results and experimental observations, is the Joule-Thomson heating associated with large pressure drops across the fracture network.

numerical simulation, embedded fractures and boreh↗
DOE OSTI2024

NORTH DAKOTA CARBONSAFE PHASE III: SITE CHARACTERIZATION AND PERMITTING OF GEOLOGIC STORAGE OF CARBON DIOXIDE

The Energy & Environmental Research Center (EERC), in partnership with Minnkota Power Cooperative Inc. (Minnkota), SLB, and Computer Modelling Group Ltd. (CMG), supported wide-scale deployment of carbon capture and storage (CCS) as part of the U.S. Department of Energy (DOE) National Energy Technology Laboratory Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative Phase III. This phase included the acquisition, analysis, and development of information to fully characterize two storage complexes to demonstrate viable storage resources for commercial volumes of CO2 (defined by DOE as a minimum of 50 million tonnes [MMt] of CO2 within a 30-year period) (National Energy Technology Laboratory, 2024). Phase III also involved the preparation, submission, and approval of North Dakota underground injection control (UIC) Class VI storage facility permits (SFPs)—required precursors to applications for Class VI injection well permits. The presumed viability of commercial-scale CCS, situated adjacent to Minnkota’s Milton R. Young Station (MRYS), is validated by Minnkota’s continued pursuit of Project Tundra—an initiative to build the world’s largest lignite-based CCS project in central North Dakota (www.projecttundrand.com). Project Tundra comprises two scopes of work, Tundra Capture (installation of postcombustion CO2 capture at MRYS) and Tundra SGS (secure geologic storage). The efforts of North Dakota CarbonSAFE Phase III, Site Characterization and Permitting, supported Tundra SGS. Extensive site-specific characterization activities included a successful multimeasurement geophysical approach and drilling a stratigraphic test well (J-ROC 1, subsequently renamed Liberty-1) adjacent to MRYS. Core collection and analyses, downhole testing and fluid sampling, and geophysical logging were performed on J-ROC 1 and on a nearby stratigraphic test well (known as J-LOC 1), which was drilled, cored, and tested under a complementary project funded by the North Dakota Lignite Research Program. The injection tests performed on J-LOC 1 positively impacted the CarbonSAFE project, resulting in fewer proposed injection wells and significant construction, operations, and monitoring cost savings. The characterization data collected and analyses performed were integrated into geologic models, and successive numerical simulations were run to determine CO2 plume extent and subsurface pressure buildup associated with the planned CO2 injection rate of nearly 4 MMt per year. The latter doubles the CarbonSAFE Initiative goal with an estimated 100 MMt of CO2 stored in 20 years. Application of the U.S. Environmental Protection Agency’s (EPA’s) method for estimating the Class VI Rule area of review (AOR) to the overpressurized Broom Creek Formation inspired an alternative method of calculation, called risk-based AOR delineation. This peer-reviewed method was applied for the first time during the storage facility-permitting process. The two SFP applications submitted in 2021 successfully resulted in North Dakota Industrial Commission (NDIC) orders in 2022 authorizing the creation of the storage facility areas and amalgamation of pore space as well as establishing financial responsibility requirements. After approval of the SFPs, Minnkota filed in 2022 applications for permits to reenter the J-ROC 1 well and to drill two new wells—all with the intended purpose to become Class VI injection wells. To establish eligibility under the Internal Revenue Code for Section 45Q tax incentives, a monitoring, reporting, and verification (MRV) plan was prepared and submitted by Minnkota to EPA in 2022, resulting in the first such plan approved in North Dakota. Also in 2022, under the National Environmental Policy Act (NEPA), Minnkota prepared and submitted an environmental information volume (EIV) describing the proposed CCS project and associated potential environmental impacts. Based on the EIV, DOE determined that the proposed construction project required an environmental assessment, and Minnkota submitted the first draft in 2023 and a revised draft in 2024. Both submissions were followed by a public comment period. Subsequently, DOE issued a finding of no significant impact (FONSI) on September 13, 2024. A successful outreach program, strongly based in the production, presentation, and dissemination of informational material, fostered an environment to aid stakeholders in making informed decisions regarding the planned project. Opportunities for public input were provided at various steps along the way, including at county planning and zoning meetings, before and during the SFP administrative hearing, and during environmental assessment public comment periods. In addition, land/pore space owners and mineral owners had various points of contact, including granting access rights, securing pore space leasing, and mineral owner notifications. Based upon the successful storage facility permitting issued by NDIC, approval of the MRV plan by EPA, and receipt of a FONSI under the NEPA, Minnkota is continuing its pursuit of Project Tundra. In December 2023, the Office of Clean Energy Demonstrations under its Carbon Capture Demonstrations Projects Program announced funding for the capture system (Office of Clean Energy Demonstrations, 2023) and a proposal for CarbonSAFE Phase IV: Construction funding was submitted in March 2024 for the storage project. A go/no-go decision to proceed with construction and operations in the Broom Creek Formation is anticipated in 2024. References National Energy Technology Laboratory, CarbonSafe Initiative, https://netl.doe.gov/carbon-management/carbon-storage/carbonsafe (accessed August 2024). Office of Clean Energy Demonstrations, 2023, OCED selects three projects in CA, ND, and TX to reduce harmful carbon pollution, create new economic opportunities, and advance carbon reducing technologies, December, www.energy.gov/oced/articles/oced-selects-three-projects-ca-nd-and-tx-reduce-harmful-carbon-pollution-create-new (accessed August 2024).

Peck, Wesley↗
DOE OSTI2025

SiC Receiver/Reactor by Additive Manufacturing for Concentrated Solar Thermocatalysis with Thermal Energy Storage (Final Technical Report - Public)

The direct use of solar thermal energy provides opportunities for low-cost heating sources for a variety of applications. Ultra-high temperatures around 1000°C are high value and highly useful for energy-demanding industries. Many materials cannot withstand these conditions. In the area of Sustainable Chemicals, further limitations on material stability exist. Combining state-of-the-art materials with new designs provides a promising pathway for harvesting solar thermal energy and performing high temperature chemical processes. However, conventional manufacturing limits the potential for design flexibility. In this project, Additive Manufacturing was combined with advanced materials and new chemical reactor designs. In addition, 24/7 energy is necessary for chemical processing, and designs for ultra-high temperature thermal storage were devised. Specifically, preliminary design of a novel solar thermal receiver was developed in this project and designed to work with thermocatalytic reactors for producing sustainable chemicals and fuels. An ultra-high temperature particle storage system and heat exchangers were proposed to transport ultra-hot air as thermal fluid for the system. On a broader scale, this system could be used to tap solar thermal energy for a centralized facility with capability of transferring that heat to various segments at a full range of temperatures to 1000°C. The project pushed the temperature boundaries past those in current use, and Additive Manufacturing was envisaged for fabricating the receiver to meet requirements of extreme environments. An extensive analysis of silicon carbide additive manufacturing was performed to compare the thermal and mechanical properties of complex geometries compared to conventional material and those manufactured via other methods. The Additive Manufacturing via Binder-Jet printing was optimized and characterized to provide high quality and reproducible components capable of withstanding the proposed extreme environments. The designs for the concentrating solar thermal cavity with ultra-hot air thermal fluid showed high performance in simulations, attributable to the complex optimized geometries of the 3D printed systems. The bright future of Additive Manufacturing with advanced materials developments should provide more options and even higher quality as the technology further develops. Current costs for Additive Manufacturing of advanced ceramics is relatively low, however post-processing of the materials for extreme environments is currently high. There is little industrial-scale infrastructure for these, but it is growing as niche applications become more mainstream. The results of the project can be translated into similar extreme environments for concentrating solar thermal energy as well as its integration with ultra-hot air thermal fluids. A number of industries that require ultra-high temperatures need to electrify or otherwise decarbonize for climate goals, and this project showed that theoretically there is a pathway to do so with direct concentrated solar thermal power.

10 SYNTHETIC FUELS↗
DOE OSTI2023

Characterizing different motility-induced regimes in active matter with machine learning and noise

Here we examine motility-induced phase separation (MIPS) in two-dimensional run-and-tumble disk systems using both machine learning and noise fluctuation analysis. Our measures suggest that within the MIPS state there are several distinct regimes as a function of density and run time, so that systems with MIPS transitions exhibit an active fluid, an active crystal, and a critical regime. The different regimes can be detected by combining an order parameter extracted from principal component analysis with a cluster stability measurement. The principal component-derived order parameter is maximized in the critical regime, remains low in the active fluid, and has an intermediate value in the active crystal regime. We demonstrate that machine learning can better capture dynamical properties of the MIPS regimes compared to more standard structural measures such as the maximum cluster size. The different regimes can also be characterized via changes in the noise power of the fluctuations in the average speed. In the critical regime, the noise power passes through a maximum and has a broad spectrum with a 1/f 1.6 signature, similar to the noise observed near depinning transitions or for solids undergoing plastic deformation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Actinide in Air (Rn-Progeny Rejected) Alpha Spectroscopy With Tensioned Metastable Fluid Detectors

Here, this article discusses outcome of research for deriving a methodology and apparatus for ascertaining for the presence of ultratrace-level actinides in air from their alpha emission signatures, while remaining blind to the relatively large (1000× higher activity) alpha emissions from Rn-progeny. Apparatus and techniques were developed to collect and characterize alpha-emitting nuclides of Rn-progeny and actinides in air on a polycarbonate 3 μm pore size continuous air monitor (CAM) filter. A wet-chemistry approach was developed and validated for successfully separating the Rn-progeny alpha emitting isotopes of Po-214 and Po-218, while extracting the actinides (U, Pu, Am) in a fluid mixture that is suitable for conduct of alpha spectroscopy with a centrifugally tensioned metastable fluid detector (CTMFD). The resulting α-TMFD technology was compared against the state-of-art “Alpha-Sentry™” Continuous Air Monitor (CAM) system commonly utilized world-wide. Results indicate that the α-TMFD technology can potentially offer complementary and superior performance in multiple performance categories, and ∼18× improvement in the time to detect (e.g., at 0.02 derived air concentration (DAC) within ∼3 h, versus ∼70 h for Alpha-Sentry) for actinides of interest while also remaining ∼100% blind to ∼10 3 × higher Rn-progeny background—with the added potential for offering few keV scale energy resolution without resorting to peak shape fitting, versus ∼300–400 keV for existing CAM systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Boron isotopes in blue diamond record seawater-derived fluids in the lower mantle

Boron is a quintessential crustal element but its conspicuous presence in diamond – a mantle mineral – raises questions about potential subduction pathways for boron and other volatiles. It has been a long-standing goal to characterize the isotopic composition of boron in blue, boron-bearing (Type IIb) diamonds to reveal its origin. Mineral inclusions indicate that Type IIb diamonds crystallize at transition zone to lower mantle depths, meaning that if the boron is subducted it would trace a pathway of volatile elements into the deep mantle. Here, using off-line laser ablation sampling, we present the first boron isotope compositions, along with trace element contents and carbon isotope compositions, of a suite of blue diamonds mainly sourced from the Cullinan diamond mine in South Africa. The ten analyzed blue diamonds have a wide range in δ 11 B, between -9.2 ± 2‰ and -0.5 ± 2‰, compared to the more restricted range for mid-ocean ridge basalts (-7.1 ± 0.9‰). Carbon isotope values for the blue diamonds range between -20.6 and -1.8‰, with a mode at -17‰, significantly more negative than the main mantle mode at -5‰. Combined, the boron and carbon isotope compositions require fluid input from subducted oceanic lithosphere down into deep mantle source regions of blue diamonds. This finding highlights a deep subduction pathway for volatiles to stoke the source regions of deeply-derived magmas.

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

A data-driven framework for permeability prediction of natural porous rocks via microstructural characterization and pore-scale simulation

Understanding the microstructure–property relationships of porous media is of great practical significance, based on which macroscopic physical properties can be directly derived from measurable microstructural informatics. However, establishing reliable microstructure–property mappings in an explicit manner is difficult, due to the intricacy, stochasticity, and heterogeneity of porous microstructures. In this paper, a data-driven computational framework is presented to investigate the inherent microstructure–permeability linkage for natural porous rocks, where multiple techniques are integrated together, including microscopy imaging, stochastic reconstruction, microstructural characterization, pore-scale simulation, feature selection, and data-driven modeling. A large number of 3D digital rocks with a wide porosity range are acquired from microscopy imaging and stochastic reconstruction techniques. A broad variety of morphological descriptors are used to quantitatively characterize pore microstructures from different perspectives, and they compose the raw feature pool for feature selection. Here high-fidelity lattice Boltzmann simulations are conducted to resolve fluid flow passing through porous media, from which reliable permeability references are obtained. The optimal feature set that best represents permeability is identified through a performance-oriented feature selection process, upon which a cost-effective surrogate model is rapidly fitted to approximate the microstructure-permeability mapping via data-driven modeling. This surrogate model exhibits great advantages over empirical/analytical formulas in terms of prediction accuracy and generalization capacity, which can predict reliable permeability values spanning four orders of magnitude. Besides, feature selection also greatly enhances the interpretability of the data-driven prediction model, from which new insights into the mechanism of how microstructural characteristics determine intrinsic permeability are obtained.

58 GEOSCIENCES↗
DOE OSTI2020

Addressing wellbore integrity and thief zone permeability using microbially-induced calcium carbonate precipitation (MICP): A field demonstration

Microbially-induced calcium carbonate precipitation (MICP) is an emerging biotechnology for wellbore integrity applications including sealing defects in wellbore cement and modifying the permeability of rock formations. The goal of this field demonstration was to characterize a failed waterflood injection well and provide proof of principle that MICP can reduce permeability in the presence of oil using conventional oilfield fluid delivery methods. Here we compared well logs performed at the time the well was drilled with ultrasonic logs, sonic cement evaluation, and temperature logs conducted after the well failed. Analysis of these logs suggested that, rather than entering the target waterflood formation, injectate was traveling through defects in the well cement to a higher permeability sandstone layer above the target formation. Sporosarcina pasteurii cultures and urea-calcium media were delivered 2290 ft (698 m) below ground surface using a 3.75 gal (14.2 L) slickline dump bailer to promote mineralization in the undesired flow paths. By Day 6 and after 25 inoculum and 49 calcium media injections, the injectivity [gpm/psi] had decreased by approximately 70%. This demonstration shows that 1) common well logs can be used to identify scenarios where MICP can be employed to reduce system permeability, remediate leakage pathways, and improve waterflood efficiency, and 2) MICP can occur in the presence of hydrocarbons.

42 ENGINEERING↗
DOE OSTI2022

Machine learning-assisted upscaling analysis of reservoir rock core properties based on micro-computed tomography imagery

Optimum solutions for geologic modeling and reservoir simulation in industries such as oil and gas recovery and carbon capture and storage require accurate characterization of reservoir properties, which are often heterogeneous. In this study, high-quality micro-computed tomography (CT) images (1.475-μm/pixel resolution) of a sandstone core acquired from the Bell Creek oil field, USA, were used to provide nondestructive analysis of pore- and core-scale heterogeneity across measurement scales of 94–566 μm. In addition to characterizing the as-received sample, the core sample was flooded with brine to evaluate the capacity of the core sample to receive injected fluids. The micro-CT images were systematically segmented into pore spaces and grains via machine learning (ML) steps including image preprocessing, label creation using a traditional ML method based on limited manual image annotation, and finally U-Net segmentation. The segmented image stacks were reconstructed into digital cubes of various scales of voxel lengths. The 3D porosity values were calculated for all the digital cubes, and the fractal dimensions of the cubes were estimated using a box-counting method. The results showed that smaller cubes had greater heterogeneity and that the porosity values could be accurately estimated by fractal dimension and voxel lengths using ML models. For the core sample with brine flooding, the ratio of pores filled by brine to the total pore space was related to the porosity and could also be accurately estimated by porosity, fractal dimension, and voxel lengths using ML models. In conclusion, the results of this study demonstrate that the concept of fractal dimension can be a useful vector to perform upscaling analysis of sandstone rock heterogeneity from the pore to core scale and that fractal dimensions can be used to estimate porosity values and pore space-filling capacity across those scales.

58 GEOSCIENCES↗
DOE OSTI2023

Motion of a viscous slug on heterogeneous surfaces: crossover from stick–slip to steady sliding

We present a theoretical study of viscous slug motion inside a microscopically rough capillary tube, where pronounced stick–slip motion can emerge at slow displacement rates. The mathematical description of this intermittent motion can be reduced to a system of ordinary differential equations, which also describe the motion of a pendulum inside a fluid-filled rotating drum. We use this analogy to show that the stick–slip motion transitions to steady sliding at high displacement rates. We characterize this crossover with a simple scaling relation and show that the crossover is accompanied by a shift in the dominant energy dissipation mechanisms within the system.

Mechanics↗
DOE OSTI2020

A novel, magnetically driven convergent Richtmyer–Meshkov platform

Here, we introduce a novel experimental platform for the study of the Richtmyer–Meshkov instability in a cylindrically converging geometry using a magnetically driven cylindrical piston. Magnetically driven solid liner implosions are used to launch a shock into a liquid deuterium working fluid and, ultimately, into an on-axis rod with a pre-imposed perturbation. The shock front trajectory is tracked through the working fluid and up to the point of impacting the rod through the use of on axis photonic Doppler velocimetry. This configuration allows for precise characterization of the shock state as it impacts the perturbed rod interface. Monochromatic x-ray radiography is used to measure the post-shock interface evolution and rod density profile. The ALEGRA MHD model is used to simulate the dynamics of the experiment in one dimension. We show that late in time the perturbation growth becomes non-linear as evidenced by the observation of high-order harmonics, up to n = 5. Two dimensional simulations performed using a combination of the GORGON MHD code and the xRAGE radiation hydrodynamics code suggest that the late time non-linear growth is modified by convergence effects as the bubbles and spikes experience differences in the pressure of the background flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗