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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Use of Neutron Flux Calculated by Shift in a Grizzly Reactor Pressure Vessel Fracture Simulation

The reactor pressure vessel (RPV) plays the critical role of containing the reactor and coolant in a light water reactor nuclear power plant, and must be able to safely perform this function under a variety of normal and off-normal transient loading conditions. The RPV is subjected to harsh environmental conditions (elevated temperature and high neutron flux) whenever the plant is operating, which leads to embrittlement of the steel over time. RPVs contain populations of flaws introduced during the manufacturing process, and the primary safety concern is that during a transient event, a fracture could initiate at the site of one of these flaws and lead to rupture of the RPV. RPVs are very rugged structures, and the probability of this occurring is very low, but as the RPV's steel becomes embrittled over time, its resistance to fracture decreases and this becomes more likely. Assessing the likelihood of RPV fracture during a transient event involves performing a probabilistic fracture mechanics analysis, in which random sampling procedures are used to assess the probability of fracture of a population of flaws. Grizzly offers unique capabilities for modeling the RPV using 1D, 2D, or 3D representations, taking advantage of parallel computing for this type of analysis. One of the important inputs to Grizzly is the spatial distribution of the neutron fluence, which has a large impact on the material embrittlement, and should thus be characterized as accurately as possible. The VERA core simulator developed by the CASL program takes advantage of high performance computing resources to perform high-fidelity multiphysics simulations of reactor physics, including the effects of neutron transport, fuel performance, and thermal hydraulics. Recently a capability was developed to also evaluate neutron fluxes away from the reactor core through the Shift code. This allows for a detailed 3D map of the spatial variation of the fluence in components such as the RPV to be computed, which is ideal for use with Grizzly. This report documents a procedure developed to use the fluences computed in Shift in a Grizzly simulation, which allows for improved realisim in the prediction of the effects of environmental exposure on propensity for reactor pressure vessel fracture. This capability demonstrated on RPVs with multiple types of configurations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Fracture Network Quantification during CO2 Injection

This is the presentation prepared for the ARMA 2025 (59th US Rock Mechanics/Geomechanics Symposium) Conference held in Santa Fe, New Mexico, June 8-11, 2025. Accurate mapping and quantification of these networks are essential to ensure the integrity of CO2 storage reservoirs, understand and reduce potential leakage, and maintain long-term environmental safety. This study presents a novel machine learning-driven approach, integrated with geomechanical analysis, to quantify fracture networks and assess their spatial distribution during CO2 injection. This paper combines microseismic monitoring data with principles of hydraulic diffusivity and geomechanical analysis to characterize reservoir scale fracture network. The novelty of our approach lies in its capacity to assimilate time-dependent pressure data and microseismicity into a cohesive framework, which not only identifies microseismic triggering fronts but also tracks fracture distribution during active injection. Besides, leveraging image log data and analysis our approach also provides another angle of the insights to solidate the fracture networks understanding and geomechanical impacts. Key results from our study include the detection of over 100 distinct fracture clusters across the injection site, with fracture orientations strongly correlated with the prevailing in-situ stress field.

CO2 storage and sequestration↗
DOE OSTI2024

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Simulating Fracture Sealing by Granular LCM Particles in Geothermal Drilling

Lost circulation occurs when the returned fluid is less than what is pumped into the well due to loss of fluid to pores or fractures. A lost-circulation event is a common occurrence in a geothermal well. Typical geothermal reservoirs are often under-pressured and have larger fracture apertures. A severe lost-circulation event is costly and may lead to stuck pipe, well instability, and well abandonment. One typical treatment is adding lost-circulation materials (LCMs) to seal fractures. Conventional LCMs fail to properly seal fractures because their mechanical limit is exceeded at elevated temperatures. In this paper, parametric studies in numerical simulations are conducted to better understand different thermal effects on the sealing mechanisms of LCMs. The computational fluid dynamics (CFDs) and the discrete element method (DEM) are coupled to accurately capture the true physics of sealing by granular materials. Due to computational limits, the traditional Eulerian–Eulerian approach treats solid particles as a group of continuum matter. With the advance of modern computational power, particle bridging is achievable with DEM to track individual particles by modeling their interactive forces between each other. Particle–fluid interactions can be modeled by coupling CFD algorithms. Fracture sealing capability is investigated by studying the effect of four individual properties including fluid viscosity, particle size, friction coefficient, and Young’s modulus. It is found that thermally degraded properties lead to inefficient fracture sealing.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

A New Modeling Framework for Multi-Scale Simulation of Hydraulic Fracturing and Production from Unconventional Reservoirs

This paper describes a new modeling framework for microscopic to reservoir-scale simulations of hydraulic fracturing and production. The approach builds upon a fusion of two existing high-performance simulators for reservoir-scale behavior: the GEOS code for hydromechanical evolution during stimulation and the TOUGH+ code for multi-phase flow during production. The reservoir-scale simulations are informed by experimental and modeling studies at the laboratory scale to incorporate important micro-scale mechanical processes and chemical reactions occurring within the fractures, the shale matrix, and at the fracture-fluid interfaces. These processes include, among others, changes in stimulated fracture permeability as a result of proppant behavior rearrangement or embedment, or mineral scale precipitation within pores and microfractures, at µm to cm scales. In our new modeling framework, such micro-scale testing and modeling provides upscaled hydromechanical parameters for the reservoir scale models. We are currently testing the new modeling framework using field data and core samples from the Hydraulic Fracturing Field Test (HFTS), a recent field-based joint research experiment with intense monitoring of hydraulic fracturing and shale production in the Wolfcamp Formation in the Permian Basin (USA). Below, we present our approach coupling the reservoir simulators GEOS and TOUGH+ informed by upscaled parameters from micro-scale experiments and modeling. We provide a brief overview of the HFTS and the available field data, and then discuss the ongoing application of our new workflow to the HFTS data set.

58 GEOSCIENCES↗
DOE OSTI2018

The distribution, orientation, and characteristics of natural fractures for Experiment 1 of the EGS Collab Project, Sanford Underground Research Facility

The EGS Collab project is focused on understanding and predicting permeability enhancement and evolution in crystalline rocks. To accomplish this, the project is creating a suite of intermediate-scale (~10-20 m) field test beds coupled with stimulation and interwell flow tests that will provide a basis to better understand the fracture geometries and processes that control heat transfer between rock and stimulated fractures. As part of the site characterization effort for the first experimental test bed, our team has worked on mapping the distribution, orientation, and nature of open and healed fractures exposed along the drift wall and within the eight bore holes drilled for this test bed. The fractures have been characterized through detailed description of continuous cores obtained from these boreholes, evaluation of televiewer logs, and mapping of fractures and seeps exposed along the drift wall. The fracture data are being compiled and interpreted for slip and dilation tendencies, and will be incorporated into coupled-process geomechanical flow and transport models to better constrain the planned flow and tracer tests.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Fracture Network Quantification during CO2 Injection

This is the conference paper accompanying an oral presentation at the ARMA 2025 (59th US Rock Mechanics/Geomechanics Symposium) Conference held in Santa Fe, New Mexico, June 8-11, 2025. Accurate mapping and quantification of these networks are essential to ensure the integrity of CO2 storage reservoirs, understand and reduce potential leakage, and maintain long-term environmental safety. This study presents a novel machine learning-driven approach, integrated with geomechanical analysis, to quantify fracture networks and assess their spatial distribution during CO2 injection. This paper combines microseismic monitoring data with principles of hydraulic diffusivity and geomechanical analysis to characterize reservoir scale fracture network. The novelty of our approach lies in its capacity to assimilate time-dependent pressure data and microseismicity into a cohesive framework, which not only identifies microseismic triggering fronts but also tracks fracture distribution during active injection. Besides, leveraging image log data and analysis our approach also provides another angle of the insights to solidate the fracture networks understanding and geomechanical impacts. Key results from our study include the detection of over 100 distinct fracture clusters across the injection site, with fracture orientations strongly correlated with the prevailing in-situ stress field.

CO2 storage and sequestration↗
DOE OSTI2019

EGS Collab Experiment 1: Common Discrete Fracture Network

This package includes data and models that support hydraulic fracture stimulation and fluid circulation experiments in the Sanford Underground Research Facility (SURF). A paper by Schwering et al. (2020) describes the deterministic basis for developing a "common" discrete fracture network (CDFN) model of significant natural fractures in EGS Collab Testbed 1 on the 4850-Level of SURF. The ReadMe for this model shows drift, wells, scanlines, fracture data, interpreted fractures, and geophysical visualizations. There is also a summary of the data that was used in this experiment and includes results from reviewing core, televiewer (TV) logs, core-TV depth/feature registration, and from mapping weeps in the 4850-Level drift. The CDFN is intended to be a baseline model of the pre-stimulated testbed (though some observations from stimulation helped inform the model).

15 GEOTHERMAL ENERGY↗
DOE OSTI2017

Frictional stability-permeability relationships for fractures in shales

There is wide concern that fluid injection in the subsurface, such as for the stimulation of shale reservoirs or for geological CO 2 sequestration (GCS), has the potential to induce seismicity that may change reservoir permeability due to fault slip. However, the impact of induced seismicity on fracture permeability evolution remains unclear due to the spectrum of modes of fault reactivation (e.g., stable versus unstable). As seismicity is controlled by the frictional response of fractures, we explore friction-stability-permeability relationships through the concurrent measurement of frictional and hydraulic properties of artificial fractures in Green River shale (GRS) and Opalinus shale (OPS). We observe that carbonate-rich GRS shows higher frictional strength but weak neutral frictional stability. The GRS fracture permeability declines during shearing while an increased sliding velocity reduces the rate of permeability decline. By comparison, the phyllosilicate-rich OPS has lower friction and strong stability while the fracture permeability is reduced due to the swelling behavior that dominates over the shearing induced permeability reduction. Furthermore, we conclude that the friction-stability-permeability relationship of a fracture is largely controlled by mineral composition and that shale mineral compositions with strong frictional stability may be particularly subject to permanent permeability reduction during fluid infiltration.

58 GEOSCIENCES↗
DOE OSTI2017

Three‐Phase CO 2 Flow in a Basalt Fracture Network

Abstract Geologic CO 2 sequestration in basalt reservoirs is predicated on permanent CO 2 isolation via rapid mineralization reactions. This process is supported by a substantial body of evidence, including laboratory experiments documenting rapid mineralization rates, regional storage estimates indicating large, accessible storage reservoirs, and two successful pilot‐scale studies. Nevertheless, there remains significant uncertainty in the behavior of CO 2 flow within basalt fracture networks, particularly in the context estimating physical trapping potential in early time and as CO 2 undergoes phase change. In this study, a Monte Carlo numerical model is designed to simulate a supercritical CO 2 plume infiltrating a low‐permeability flood basalt entablature. The fracture network model is based on outcrop‐scale LiDAR mapping of Columbia River Basalt, and CO 2 flow is simulated within fifty equally probable realizations of the fracture network. The spatial distribution of fracture permeability for each realization is randomly drawn from a basalt aperture distribution, and ensemble results are analyzed with e‐type estimates to compute mean and standard deviation of fluid pressure and CO 2 saturation. Results of this model after 10 years of simulation suggest that (1) CO 2 flow converges on a single dominant flow path, (2) CO 2 accumulates at fracture intersections, and (3) variability in permeability can account for a 1.6 m depth interval within which free CO 2 may change phase from supercritical fluid to subcritical liquid or gas. In the context of CO 2 sequestration in basalt, these results suggest that physical CO 2 trapping may be substantially enhanced as carbonate minerals precipitate within the basalt fracture network.

Gierzynski, Alec O.↗
DOE OSTI2016

Diffusion-Weighted MRI Assessment of Adjacent Disc Degeneration After Thoracolumbar Vertebral Fractures

ObjectiveThe purpose of this study was to assess, by the mean apparent diffusion coefficient (ADC), if a relationship exists between disc ADC and MR findings of adjacent disc degeneration after thoracolumbar fractures treated by anatomic reduction using vertebral augmentation (VAP).Materials and MethodsTwenty non-consecutive patients (mean age 50.7 years; range 45–56) treated because of vertebral fractures, were included in this study. There were 10 A3.1 and 10 A1.2 fractures (AO classification). Surgical treatment using VAP was applied in 14 cases, and conservative in 6 patients. MRI T2-weighted images and mapping of apparent diffusion coefficient (ADC) of the intervertebral disc adjacent to the fractured segment were performed after a mean follow-up of 32 months. A total of 60 discs, 3 per patient, were analysed: infra-adjacent, supra-adjacent and a control disc one level above the supra-adjacent.ResultsNo differences between patients surgically treated and those following a conservative protocol regarding the average ADC values obtained in the 20 control discs analysed were found. Considering all discs, average ADC in the supra-adjacent level was lower than in the infra-adjacent (1.35 ± 0.12 vs. 1.53 ± 0.06; p < 0.001). Average ADC values of the discs used as a control were similar to those of the infra-adjacent level (1.54 ± 0.06). Compared to surgically treated patients, discs at the supra-adjacent fracture level showed statistically significant lower values in cases treated conservatively (p < 0.001). The variation in the delay of surgery had no influence on the average values of ADC at any of the measured levels.ConclusionsADC measurements of the supra-adjacent discs after a mean follow-up of 32 months following thoracolumbar fractures, showed that restoration of the vertebral collapse by minimally invasive VAP prevents posttraumatic disc degeneration.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

A multi-resolution approach to hydraulic fracture simulation

Abstract We present a multi-resolution approach for constructing model-based simulations of hydraulic fracturing, wherein flow through porous media is coupled with fluid-driven fracture. The approach consists of a hybrid scheme that couples a discrete crack representation in a global domain to a phase-field representation in a local subdomain near the crack tip. The multi-resolution approach addresses issues such as the computational expense of accurate hydraulic fracture simulations and the difficulties associated with reconstructing crack apertures from diffuse fracture representations. In the global domain, a coupled system of equations for displacements and pressures is considered. The crack geometry is assumed to be fixed and the displacement field is enriched with discontinuous functions. Around the crack tips in the local subdomains, phase-field sub-problems are instantiated on the fly to propagate fractures in arbitrary, mesh independent directions. The governing equations and fields in the global and local domains are approximated using a combination of finite-volume and finite element discretizations. The efficacy of the method is illustrated through various benchmark problems in hydraulic fracturing, as well as a new study of fluid-driven crack growth around a stiff inclusion.

58 GEOSCIENCES↗
DOE OSTI2021

Characterizing Reactive Transport Behavior in a Three-Dimensional Discrete Fracture Network

While several studies have linked network and in-fracture scale properties to conservative transport behavior in subsurface fractured media, studies on reactive transport cases remain relatively underdeveloped. In this study, we explore the behavior of an irreversible kinetic reaction during the interaction of two solute plumes, one consisting of species A and the other species B. When the plumes converge, these species react kinetically to form a new species C via A+B→kC. This reactive system is studied using a three-dimensional discrete fracture network (DFN) model coupled with reactive Lagrangian particle tracking. We find that the interplay of network topology and chemical properties of the reactive solutes controls reactive transport processes. The network topology drives species A and B together, and the chemical properties dictate whether and how quickly a reaction occurs. Results demonstrate that reactions are most likely to occur in high-velocity fractures that make up the network backbone. The interplay between species’ chemical properties and transport is characterized by a non-dimensional Damköhler (Da) number. We show that the spatial distribution of reactions is sensitive to Da, which subsequently influences late-time tailing behavior in outlet breakthrough time distributions. The results of this study provide initial insights into how an irreversible reaction occurs during transport in a fracture network, using a methodology that can be applied to study reactive transport in a wide range of fractured media environments and contexts.

58 GEOSCIENCES↗
DOE OSTI2017

Effect of Normalizing Temperature on Fracture Characteristic of Tensile and Impact Tested Creep Strength-Enhanced Ferritic P92 Steel

The high-temperature Cr-Mo creep strength-enhanced ferritic (CSEF) steels are mainly used in nuclear and thermal power plants. In the present investigation, a systematic study on fracture surface morphologies of tensile and impact tested specimens and mechanical properties of cast and forged (C&F) P92 steel was performed for various heat treatment conditions. The heat treatment was carried out in normalizing temperature range of 950-1150 °C and then tempered to a fixed tempering temperature of 760 °C. The effect of varying normalizing temperatures before and after tempering on microstructure evolution, tensile properties, Vicker’s hardness and Charpy toughness was studied. The normalizing temperature before and after tempering was having a noticeable effect on mechanical properties of as-received P92 steel. The fracture surface of impact and tensile tested samples was also studied for various normalizing temperatures with or without tempering. Fracture surface morphology was affected by the presence of secondary phase carbide particles. The fraction area of cleavage facets on the tensile fracture surface was found to be increased with an increase in the normalizing temperature. The fractured tensile specimens were characterized by transgranular ductile dimples, tear ridges and transgranular cleavage facets for various heat treatments. The fracture mode of impact tested samples was more complex. It showed both quasi-cleavage facets and ductile dimple tearing for various normalizing temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Effect of the As-Forged and Heat-Treated Microstructure on the Room Temperature Anisotropic Ductile Fracture of Inconel 718

The purpose of the present work was to investigate the effect of primary carbides and the δ-phase on the anisotropic ductile fracture of Inconel 718 in the forging process. Inconel 718 alloys were prepared by VIM + VAR processes with various carbon contents (0.009 and 0.027 wt.%). Then, the alloys were forged and annealed at temperatures of 980 and 1030 °C. The room temperature mechanical anisotropy of the alloys was evaluated at the longitudinal direction (LD) and transverse direction (TD). Tensile and impact tests were used to characterize the mechanical properties of the specimens. The microstructural characterization and the fractography of the alloys were carried out by FE-SEM. The obtained results showed that the fracture strain and the impact energy in the TD were 30-50% lower than the LD. The fracture was accelerated by the δ-phase, leading to the reduction of impact energy in the longitudinal and the lateral directions up to 50%. The low-carbon alloy indicated similar characteristics in both the LD and the TD. Aligned carbides changed the fracture path from a zigzag path in the LD to a fibrous path in the TD, while the δ-phase created a flat fracture path. The shear lip area ratio in the tensile fracture cross section was decreased by reducing ductility.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Grain Boundary and Lattice Fracture Toughness of UO 2 Measured Using Small-Scale Mechanics

Small-scale mechanical testing has been used to isolate the properties of individual grain boundaries and regions of the lattice in undoped UO 2 polycrystals. Notched microcantilever specimens were prepared by focused ion beam milling and tested in bending. Fracture toughness values were calculated from finite element analysis of the sample geometry at fracture load using the J-integral method. Lattice and grain boundary fracture toughness values of 1.17 ± 0.16 MPa m 1/2 and 1.08 ± 0.08 MPa m 1/2 , respectively, were obtained. The similarity between the lattice and grain boundary fracture toughness values is consistent with the preference for transgranular room-temperature fracture in UO 2 . The grain boundary and lattice fracture toughness values also agree well with the toughness of dense polycrystals, which suggests that microstructural-scale toughening mechanisms are limited in such samples.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Vertebral Compression Fractures in Elderly: How to Recognize and Report

Purpose of Review: In elderly patients, when there is no obvious evidence of malignancy on plain radiography, it can be difficult to evaluate benign and malignant vertebral fractures, because these patients are predisposed to senile osteoporosis and malignant disease. The only complaint could be back pain, as both types of fractures have the same symptoms and signs. Recent Findings: A reliable imaging modality is always required, as appropriate management and early diagnosis of malignant vertebral fractures are important. Summary: The aim of this study was to highlight the importance of recognizing, through diagnostic techniques, the presence of vertebral fractures, since the majority of vertebral fractures do not come to clinical attention and to determine the parameters that are useful for the early differentiation between benign versus malignant vertebral fractures.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2023

A deep learning upscaling framework: Reactive transport and mineral precipitation in fracture-matrix systems

Pore-scale modeling has limited applicability at large scales due to its high computational cost. One common approach to upscale pore-scale models is the use of effective medium theories, which homogenize small-scale features in a porous structure and characterize the medium by macroscale properties (e.g., permeability) and equations (e.g., Darcy’s law). However, there are classes of physical processes for which effective medium approximations may become inaccurate, e.g., mineral precipitation and clogging during reactive transport. We have developed a deep learning upscaling framework, in which pore-scale modeling is directly employed in macroscale systems, without relying on effective medium approximations. The upscaling framework is first developed for general multiscale systems and then applied to modeling reactive transport with mineral precipitation in the altered layer in fracture-matrix structures. Solute transport from the fractures to the matrix is modeled as a wall boundary condition for the fractures, which, in turn, is predicted by recurrent neural networks using the concentration histories at the fracture-matrix boundary. Specifically, we consider a meter-scale fracture network embedded in sandstones, where the smallest feature is at the micron scale. Here the proposed framework allows us to span five orders of magnitude in length scales by capturing mineral precipitation in the altered layer of the rock matrix at the pore scale across the entire meter-scale fracture network.

42 ENGINEERING↗