Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracturing”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Fe Oxidation and Species Distribution at the Rock–Fluid Interface of Marcellus Shale Reacted with Hydraulic Fracturing Fluid

Hydraulic fracturing of shale reservoirs resulted in significant opportunity for increased oil and gas production in the United States. Rock-fluid interactions can cause mineral dissolution and precipitation reactions that lead to permeability changes in the shale matrix, which ultimately may affect transport pathways and hydrocarbon production. Understanding the distribution of secondary precipitates, such as barite and Fe(III) (hydro)oxides, and cation leaching at the rock-fluid interface is an important step to further investigate how these geochemical processes can change permeability and transport pathways. In this study, thin sections of the fracture-matrix interface were made from reacted Marcellus shale cores. The thin sections were characterized using synchrotron X-ray fluorescence imaging and synchrotron X-ray absorption spectroscopy. Fe species with different oxidation states were identified in the maps, together with barite and Ca distribution. The results show that ferrihydrite, as newly formed Fe(III)-bearing precipitates, aligned well with the border of the Ca (e.g., calcite) leaching region in the reaction front. Some Fe-containing clay also dissolved, but the dissolution region for the clay was not as deep as the calcite. Further, the reaction front is about three times deeper in the direction parallel to the shale bedding than that perpendicular to the bedding. The Ca leaching region can be an index for reaction front detection for Marcellus shale. Reactive transport modeling was conducted and the predicted Ca leaching boarder align well with ferrihydrite precipitation, consistent with the experimental observation. The carbonate mineral dissolution can be crucial to promote fluid access into the shale matrix. Together with our previous study on the shale reactive surface, this follow-up study showed similar Ca leaching region and Fe(III) precipitates distribution in the reaction front regardless of barite precipitation on the surface, indicating that the barite coatings on the surface may not pose a significant impact on reactive transport at the shale-fluid interface.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2020

Chemical Speciation and Stability of Uranium in Unconventional Shales: Impact of Hydraulic Fracture Fluid

Uranium and other radionuclides are prominent in many unconventional oil/gas shales and is a potential contaminant in flowback/produced waters due to the large volumes/types of chemicals injected into the subsurface during stimulation. In this work, to understand the stability of U before and after stimulation, a geochemical study of U speciation was carried out on three shales (Marcellus, Green River, and Barnett). To do so, two types of samples for each shale were subjected to sequential chemical extractions: unreacted and shale-reacted with a synthetic hydraulic fracture fluid. Furthermore, a significant proportion of the total U (20–57%) was released from these three shales after reaction with fracture fluid, indicating that U is readily leachable. The total U released exceeds labile water-soluble and exchangeable fractions in unreacted samples, indicating that fluids leach more recalcitrant phases in the shale. Radiographic analysis of unreacted Marcellus shale thin sections shows U associated with detrital quartz and the clay matrix in the shale. Detrital zircon and TiO2 identified by an electron microprobe could account for the hot spots. This study shows that significant proportions of U in three shales are mobile upon stimulation. In addition, the extent of mobilization of U depends on the U species in these rocks.

58 GEOSCIENCES↗
DOE OSTI2020

Anomalous Transport in Three-Dimensional Discrete Fracture Networks: Interplay Between Aperture Heterogeneity and Injection Modes

In this work, we study how the interplay between fracture aperture heterogeneity and tracer injection mode controls fluid flow and tracer transport in three-dimensional (3D) discrete fracture networks (DFNs). The direct 3-D DFN simulations show that tracer injection mode has substantial effects on tracer spreading across all levels of aperture heterogeneity. The key controlling factor for effective transport is the initial Lagrangian velocity distribution, which is determined by the interplay between injection mode and aperture heterogeneity. The fundamental difference between initial Lagrangian velocity distribution and domain-scale Eulerian velocity distribution plays a vital role in determining anomalous transport. We effectively capture the observed anomalous transport using an upscaled transport model that incorporates initial velocity distribution, stationary velocity distribution, velocity correlation length, and average advective tortuosity. With the upscaled transport model, we accurately capture the evolution of Lagrangian velocity distribution and predict longitudinal spreading in 3-D DFN.

58 GEOSCIENCES↗
DOE OSTI2022

Phase-field framework with constraints and its applications to ductile fracture in polycrystals and fatigue

Abstract Modeling of ductile fracture in polycrystalline structures is challenging, since it requires integrated modeling of cracks, crystal plasticity, and grains. Here we extend the typical phase-field framework to the situations with constraints on the order parameters, and formulate two types of phase-field models on ductile fracture. The Type-I model incorporates three sets of order parameters, which describe the distributions of cracks, plastic strain, and grains, respectively. Crystal plasticity is employed within grain interiors accommodated by J 2 plasticity at grain boundaries. The applications of the Type-I model to single crystals and bicrystals demonstrate the influences of grain orientations and grain boundaries on crack growth. In the Type-II model, J 2 plasticity is assumed for the whole system and grain structures are neglected. Taking advantage of the efficiency of the fast Fourier transform, our Type-II model is employed to study low cycle fatigue. Crack closure and striation-like patterning of plastic strain are observed in the simulations. Crack growth rate is analyzed as a function of the J -integral, and the simulated fatigue life as a function of plastic strain agrees with the Coffin–Manson relation without a priori assumption.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
DOE OSTI2022

Atomistic measurement and modeling of intrinsic fracture toughness of two-dimensional materials

Quantifying the intrinsic mechanical properties of two-dimensional (2D) materials is essential to predict the long-term reliability of materials and systems in emerging applications ranging from energy to health to next-generation sensors and electronics. Currently, measurements of fracture toughness and identification of associated atomistic mechanisms remain challenging. Herein, we report an integrated experimental–computational framework in which in-situ high-resolution transmission electron microscopy (HRTEM) measurements of the intrinsic fracture energy of monolayer MoS 2 and MoSe 2 are in good agreement with atomistic model predictions based on an accurately parameterized interatomic potential. Changes in crystalline structures at the crack tip and crack edges, as observed in in-situ HRTEM crack extension tests, are properly predicted. Such a good agreement is the result of including large deformation pathways and phase transitions in the parameterization of the inter-atomic potential. The established framework emerges as a robust approach to determine the predictive capabilities of molecular dynamics models employed in the screening of 2D materials, in the spirit of the materials genome initiative. Moreover, it enables device-level predictions with superior accuracy (e.g., fatigue lifetime predictions of electro- and opto-electronic nanodevices).

Zhang, Xu↗
DOE OSTI2023

A probabilistic-phase field model for the fracture of brittle materials

Abstract We develop a computational method to determine the failure probability of brittle materials under general mechanical loading conditions. The method is a combination of two parts: (1) numerical simulations of materials with multiple cracks using phase field theory, where the complete fracture process is viewed as ‘damage percolation’ along critical paths or clusters of cracks, rather than the traditional weak-link failure mechanism of Weibull, and (2) an extension of the Batdorf statistical theory of fracture to finite domains, where it is implemented within the finite element framework. The results of phase-field simulations at the ‘percolation threshold’ are used as failure data in the Batdorf theory to determine the overall probability of failure. The input to this approach is the size distribution of cracks in a pristine material. An example is shown, where alumina samples that were previously tested by Abe and coworkers (Abeet al2003J. Am. Ceram. Soc.861019–21) in four-point loading are compared to the results of our numerical simulations. The approach developed here has the advantage of being extendable to more complex thermomechanical loading.

Materials Science↗
DOE OSTI2020

Quantum-well states in fractured crystals of the heavy-fermion material CeCoIn5

Quantum-well states appear in metallic thin films due to the confinement of the wave function by the film interfaces. Using angle-resolved photoemission spectroscopy, we unexpectedly observe quantum-well states in fractured single crystals of CeCoIn5. We confirm that confinement occurs by showing that these states’ binding energies are photon-energy independent and are well described with a phase accumulation model, commonly applied to quantum-well states in thin films. This indicates that atomically flat thin films can be formed by fracturing hard single crystals. For the two samples studied, our observations are explained by free-standing flakes with thicknesses of 206 and 101 Å. We extend our analysis to extract bulk properties of CeCoIn5. Specifically, we obtain the dispersion of a three-dimensional band near the zone center along in-plane and out-of-plane momenta. We establish part of its Fermi surface, which corresponds to a hole pocket centered at G. We also reveal a change of its dispersion with temperature, a signature that may be caused by the Kondo hybridization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Universality class for loopless invasion percolation models and a percolation avalanche burst model for hydraulic fracturing

Invasion percolation is a model that was originally proposed to describe growing networks of fractures. In this work, we describe a loopless algorithm on random lattices, coupled with an avalanche-based model for bursts. The model reproduces the characteristic b-value seismicity and spatial distribution of bursts consistent with earthquakes resulting from hydraulic fracturing (“fracking”). We test models for both site invasion percolation and bond invasion percolation. These have differences on the scale of site and bond lengths l. But since the networks are characterized by their large-scale behavior, l << L, we find small differences between scaling exponents. Though data may not differentiate between models, our results suggest that both models belong to different universality classes.

58 GEOSCIENCES↗
DOE OSTI2024

3D modeling of deep borehole electromagnetic measurements with energized casing source for fracture mapping at the Utah Frontier Observatory for Research in Geothermal Energy

Here, we present a 3D numerical modelling analysis evaluating the deployment of a borehole electromagnetic measurement tool to detect and image a stimulated zone at the Utah Frontier Observatory for Research in Geothermal Energy geothermal site. As the depth to the geothermal reservoir is several kilometres and the size of the stimulated zone is limited to several 100 m, surface-based controlled-source electromagnetic measurements lack the sensitivity for detecting changes in electrical resistivity caused by the stimulation. To overcome the limitation, the study evaluates the feasibility of using a three-component borehole magnetic receiver system at the Frontier Observatory for Research in Geothermal Energy site. To provide sufficient currents inside and around the enhanced geothermal reservoir, we use an injection well as an energized casing source. To efficiently simulate energizing the injection well in a realistic 3D resistivity model, we introduce a novel modelling workflow that leverages the strengths of both 3D cylindrical-mesh-based electromagnetic modelling code and 3D tetrahedral-mesh-based electromagnetic modelling code. The former is particularly well-suited for modelling hollow cylindrical objects like casings, whereas the latter excels at representing more complex 3D geological structures. In this workflow, our initial step involves computing current densities along a vertical steel-cased well using a 3D cylindrical electromagnetic modelling code. Subsequently, we distribute a series of equivalent current sources along the well's trajectory within a complex 3D resistivity model. We then discretize this model using a tetrahedral mesh and simulate the borehole electromagnetic responses excited by the casing source using a 3D finite-element electromagnetic code. This multi-step approach enables us to simulate 3D casing source electromagnetic responses within a complex 3D resistivity model, without the need for explicit discretization of the well using an excessive number of fine cells. We discuss the applicability and limitations of this proposed workflow within an electromagnetic modelling scenario where an energized well is deviated, such as at the Frontier Observatory for Research in Geothermal Energy site. Using the workflow, we demonstrate that the combined use of the energized casing source and the borehole electromagnetic receiver system offer measurable magnetic field amplitudes and sensitivity to the deep localized stimulated zone. The measurements can also distinguish between parallel-fracture anisotropic reservoirs and isotropic cases, providing valuable insights into the fracture system of the stimulated zone. Besides the magnetic field measurements, vertical electric field measurements in the open well sections are also highly sensitive to the stimulated zone and can be used as additional data for detecting and imaging the target. We can also acquire additional multiple-source data by grounding the surface electrode at various locations and repeating borehole electromagnetic measurements. This approach can increase the number of monitoring data by several factors, providing a more comprehensive dataset for analysing the deep-localized stimulated zone. The numerical analysis indicates that it is feasible to use the combination of the energized casing and downhole electromagnetic measurements in monitoring localized stimulated zone at large depths.

58 GEOSCIENCES↗
DOE OSTI2022

Specimen Size and Geometry Effects on the Master Curve Fracture Toughness Measurements of EUROFER97 and F82H Steels

EUROFER97 and F82H are two leading reduced-activation ferritic-martensitic (RAFM) steels for fusion blanket applications. Exposure to the harsh environment of fusion reactors can result in severe degradation of materials fracture toughness (FT). Thus, the post-irradiation evaluation of FT is critical to understanding the material behavior. Due to the space constraint of irradiation facilities, the development of small specimen test techniques (SSTT) is necessary to evaluate the performance of irradiated materials. In this study, we evaluated the specimen size and geometry effects on the ductile-to-brittle transition FT of EUROFER97 batch-3 and F82H-BA12 steels. The specimen thicknesses ranged from 1.65 to 12.7 mm and the geometries included 1.65 mm bend bar, 4 mm mini-compact tension (miniCT), and 0.5T compact tension (CT) specimens. Fracture toughness testing and evaluations were performed using the Master Curve method in the ASTM E1921-19 standard. After size correction to 1T size using the Master Curve method, no specimen size effect was observed between the 4 mm miniCT and 0.5T CT specimens for the Master Curve reference temperature T0Q, while the bend bars yielded a higher T0Q. A strong effect of fatigue precrack front straightness on T0Q for 0.5T CT specimens was observed. The minimum number of specimens needed for each specimen geometry has been determined.

Chen, Xiang↗
DOE OSTI2024

Code for Experiment in Publication “Determining the dominant factors for carbon-1 mineralization in three-dimensional fracture networks”

We plan to release the code used to perform the experiment described in our recent publication, entitled “Determining the dominant factors for carbon-1 mineralization in three-dimensional fracture networks.” This code fits a joint emulator to data from several Discrete Fracture Network (DFN) simulations, performed using the open-source software DFNworks (https://dfnworks.lanl.gov/). All code to be released implements existing methods; there are no novel algorithms nor any major innovations to existing software.

Murph, Alexander↗
DOE OSTI2024

Utah FORGE 5-2615: Determination and Analysis of Thermo-poromechanical Response of Fractured Rock - 2024 Annual Workshop Presentation

This is a presentation on the Determination and Modeling-Informed Analysis of Thermo-poromechanical Response of Fractured Rock for Application to FORGE by the University of Oklahoma, presented by Ahmad Ghassemi. This video presentation discusses how to improve understanding and control of coupled thermoporomechanical (or thermo-hydro-mechanical-THM) processes in reservoir development, and to demonstrate its role in interpretations of the fracture closure pressure. This presentation was featured in the Utah FORGE R&D Annual Workshop on August 13-15, 2024.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Fracture properties of high-entropy alloys

Abstract Since the concept of high-entropy alloys (HEAs) as materials with at least four or five principal elements in (near)-equiatomic composition was introduced in 2004, this new class of materials has penetrated essentially all materials science-related fields. The main reason for this is that some face-centered-cubic alloy compositions have been shown to exhibit truly outstanding mechanical properties with extraordinary combinations of strength, ductility, and fracture toughness, particularly at cryogenic temperatures, whereas certain body-centered-cubic refractory compositions display remarkable high-temperature strength. While significant efforts have been put into rapid screening and narrowing the compositional space of HEAs to a manageable scope, there are still only a few metallic alloys that push the limits of mechanical performance. Here, we review work on some of the most damage-tolerant HEAs discovered to date and discuss the fundamental reasons why their resistance to fracture and subsequent stable crack growth is so exceptional. Graphical abstract

Gludovatz, Bernd (ORCID:0000000224203879)↗
DOE OSTI2021

Deformation and fracture characteristics of zirconium plate produced via ultrasonic additive manufacturing

Abstract The microstructural evolution, deformation modes, and fracture mechanisms of zirconium plate produced using ultrasonic additive manufacturing (UAM) are presented. In addition to conventional tensile testing techniques, digital image correlation captured highly variable strain accumulation in specimens loaded perpendicular or parallel to the build height (Z). When tested in parallel to Z, delamination at prior foil/foil interfaces creates strain localization noticeable in strain rate maps, whereas specimens loaded perpendicular to Z illustrate conventional strain hardening until necking accelerates delamination. Although bond strengths are statistically and spatially variable, in situ electron backscattering diffraction tests illustrate the ability for grains near interfaces to accommodate strain with twinning and slip modes consistent with conventionally produced zirconium alloys. Finally, mixtures of ductile and delamination-induced fracture highlight the interface-driven failure modes of UAM zirconium plate in the as-built condition. Graphic abstract

36 MATERIALS SCIENCE↗
DOE OSTI2020

Potential Formation Damage: An Integrated Reservoir Characterization Study of the Naturally Fractured Carbonate Middle Duperow Formation at the Kevin Dome, Montana

In this study, we integrate geologic and engineering data of a naturally fractured carbonate reservoir at the Kevin Dome, Montana. Well test data are correlated with core description, geochemical and lithology study to determine the flow behavior and communication within the injection test interval and to the surrounding area. Based on a dual-continuum geologic model, numerical brine injection simulations are carried out to validate the interpretation results from our well test analytical models and forecast the probability of CO 2 injection success using current reservoir properties. As a result, our well test analytical models as well as lithology/core description suggest that fluid flow may be mainly restricted to the injection interval and the assumption of radial (horizontal) flow may be appropriate. The well test models also indicate that there is potentially formation damage with a positive skin factor although prior to brine injection well tests, well stimulation through acid treatment was performed. Our numerical simulation results appear to confirm this formation damage by showing additional pressure buildup in the injection data during later test periods. To explain this, acid may have dissolved dolomite then dolomite or calcite may have been formed again further into the matrix/fracture system. Another possible explanation is mechanical clogging of the fractures due to acid dissolving dolomite and dislodging fine grains. Our work also predicts that if no additional well stimulation is performed, the project will have a lower probability of successfully injecting 1 million tons of CO 2 into the Middle Duperow formation over 4 years.

42 ENGINEERING↗
DOE OSTI2021

Rock Fracturing Using High-Pressure Ethylene/Nitrous Oxide Detonations

The present work investigates high initial pressure detonations of a stoichiometric mixture of ethylene and nitrous oxide (C2H4 + 6N2O) as a method of fracturing rock beneath the ground surface. These tests were conducted at a test site operated by the Energetic Materials Research and Testing Center (EMRTC), Socorro, New Mexico. The volume under the surface used for testing (called the Down Hole Assembly) consists of a 0.438 in. ID x 50 ft. long stainless-steel tube running down from the test site to a well bore which is 3 in. ID x 10 ft. long and the rock in the well bore is exposed to the propagating combustion wave. The testing carried out at Zucrow Laboratories in the smaller, alloy steel combustion vessel provided a scaling of pressures expected in the well bore. The combustion is initiated by energizing an EBW (Exploding Bridge Wire) above the ground surface. The experimental setup accommodates one high pressure (100,000 psia) transducer to measure the pressure peak and is placed approximately 5 ft. above the ground surface and 5 ft. downstream of the EBW. The focus of this series of experiments is to investigate the dependence of fracture to the rock beneath the surface on initial pressures of the mixture of ethylene and nitrous oxide. Experiments were carried out at initial pressures varying between 125 psia and 300 psi. The transducer recorded elevated pressures, which were 2.3 to 2.6 times in excess of the CJ values. The experimental results are discussed and explained in this report.

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

Fracture Strength Determination Methods for Ceramic Materials Applied to Uranium Dioxide

The Advanced Fuels Campaign is currently focusing on development of accident-tolerant fuels that possess a range of property modifications intended to improve fuel performance during accident and transient conditions as well as extend license limits to burnups beyond 62 MWd/kgU. Both of these drivers have identified understanding and mitigating fuel cracking as a key performance criterion. Small-scale cantilever beam testing has been identified as a plausible method by which to collect fracture data for UO 2 as a function of chemical and structural evolutions introduced either during fabrication or irradiation. While this method has been found to be capable of providing data that are in reasonable agreement with literature for unirradiated UO 2 , the inherently small sample volumes that can be sampled limit its ability to capture the statistical nature of mechanisms that govern the fracture of brittle ceramics. A biaxial flexure strength test was developed to be used for unirradiated UO 2 . This method is standard in the community, but no systems presently in use at national laboratories, universities, or private companies are available to be used for nuclear fuel materials. This report describes the operation and benchmarking of this system, which has been validated for a number of common oxides. Future work will extend this system to characterization of both doped UO 2 and other relevant microstructural modifications that can also be measured using small-scale cantilever beam testing. Comparison of datasets collected using both methods will allow the overall applicability of small-scale techniques to be assessed and provide confidence or bounds on their use for irradiated fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗