Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture behavior”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Transverse Rupture Strength of Uranium Dioxide

Uranium dioxide (UO 2 ) fuel is used as fuel in light water reactors (LWRs). While the fuel pellet is technically the first engineering barrier for radionuclide release, pellet fracturing at intermediate- to high-burnup values releases fission gases into the fuel rod plenum [1, 2]. Therefore, the true engineering barrier is the fuel cladding, which performs very well in LWR environments [3]. The extreme temperature gradients generated by fission energy and the low thermal conductivity of UO 2 quickly induce radial cracking in UO 2 during operation [4]. Cracks in the fuel provide opportunities for fuel relocation, increased fission gas release, and pellet-cladding mechanical interaction (PCMI) [5]. The ability to predict and engineer the fracture of UO 2 fuel pellets using modern computational tools is therefore a key engineering goal that has been the focus of ongoing experimental and computational efforts [6, 7]. Accurate predictions of fuel pellet cracking during operation requires knowledge of more complex phenomena, but improved understanding of the fundamental fracture behavior of unirradiated UO 2 is first necessary.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Connecting In Situ Stress and Wellbore Deviation to Near-Well Fracture Complexity Using Phase-Field Simulations

The interactions among in situ stress, rock fabric, wellbore geometry, natural fractures, and other natural or man-made defects create highly complex fracture trajectories in the near-wellbore region, far more intricate than those in the far-field. These near-wellbore complexities are critical for the Utah FORGE project and Enhanced Geothermal Systems (EGS) in general. Frictional pressure loss in the near-wellbore region during stimulation can significantly influence the growth of far-field fractures, while pressure losses during circulation serve as a major source of energy dissipation. Near-wellbore fracture complexities are often observable through image logs, offering valuable insights into in situ stress characteristics. However, leveraging this information requires a high-fidelity model capable of capturing the interplay among the diverse factors influencing fracture behavior.

58 GEOSCIENCES↗
DOE OSTI2025

Microstructure & Mechanical Properties of Scaled Thermally Aged LPBF 316H SS Builds

As part of the US Department of Energy’s Advanced Materials and Manufacturing Technologies program’s mission to accelerate qualification of advanced manufacturing pathways for nuclear applications, laser powder bed fusion (LPBF) 316H stainless steel (SS) has been selected as a model system to develop a rapid code case framework. This effort directly addresses the grand challenges of (1)expanding the limited portfolio of materials currently codified for elevated-temperature nuclear structural service under Section III, Division 5of the American Society of Mechanical Engineers’ Boiler and Pressure Vessel Code; and (2) significantly reducing qualification timelines that traditionally exceed a decade. The strategic importance of LPBF 316H lies in its immediate industrial relevance, existing data foundation from wrought 316H, and alignment with ongoing code case development for LPBF 316L.Prior work revealed accelerated precipitation of deleterious secondary phases and reduced creep ductility in as-printed LPBF 316H. Building on that prior research, FY2025activities focused on establishing an understanding of the key failure mechanisms of crept 316H specimens to aid in code case development and on evaluating stress relief (SR) parameters on the high-temperature performance and thermal aging-induced degradation of tensile and fracture behavior in LPBF316H.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microstructure and Mechanical Properties of Scaled Thermally Aged LPBF 316H SS Builds

As part of the US Department of Energy’s Advanced Materials and Manufacturing Technologies program’s mission to accelerate qualification of advanced manufacturing pathways for nuclear applications, laser powder bed fusion (LPBF) 316H stainless steel (SS) has been selected as a model system to develop a rapid code case framework. This effort directly addresses the grand challenges of (1) expanding the limited portfolio of materials currently codified for elevated-temperature nuclear structural service under Section III, Division 5 of the American Society of Mechanical Engineers’ Boiler and Pressure Vessel Code; and (2) significantly reducing qualification timelines that traditionally exceed a decade. The strategic importance of LPBF 316H lies in its immediate industrial relevance, existing data foundation from wrought 316H, and alignment with ongoing code case development for LPBF 316L. Prior work revealed accelerated precipitation of deleterious secondary phases and reduced creep ductility in as-printed LPBF 316H. Building on that prior research, FY 2025 activities focused on establishing an understanding of the key failure mechanisms of crept 316H specimens to aid in code case development and on evaluating stress relief (SR) parameters on the high-temperature performance and thermal aging induced degradation of tensile and fracture behavior in LPBF 316H.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Creating Hydride Rim Structures in Zircaloy 4 Cladding Tubes to Mimic High Burnup

During their long service in commercial light water reactors, zirconium alloy cladding tubes are subject to oxidation and subsequent hydrogen pickup which can degrade the materials mechanical properties [1]. These degradations can lead to premature cladding ruptures during reactor power transients [2] as well as negatively impact the claddings post quench ductility following high temperature steam oxidation in thermal hydraulic transients [3]. Hydrogen has low solid solubility in zirconium and usually precipitates as brittle zirconium hydride platelets in the cladding metal matrix. Due to temperature driven diffusion (Soret effect) and the grain texture of stress relieved zirconium metal tubes, these hydrides precipitate with a circumferential orientation near the outer surface of the cladding tubes forming a hydride rim which can serve as a site for crack initiation during circumferential ‘hoop’ loading of the cladding [4rim formation][5]. In order to better understand the effect of hydride rims on the fracture behavior of stress relieved zirconium alloy cladding tubes a procedure to produce these rim-like structures is being developed at Idaho National Laboratory (INL).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Quantitative Assessment of Parent Well Effect on Hydraulic Fracture Propagation at HFTS2: Insights from Cross-Well Strain Measurements and Microseismic Data

Understanding fracture propagation behavior is essential for optimizing hydraulic fracturing in unconventional reservoirs. This study demonstrates the value of integrating Low-Frequency Distributed Acoustic Sensing (LF-DAS) and microseismic data, which together provide a more complete picture of fracture growth. Using data from Hydraulic Fracturing Test Site 2 (HFTS2), we identify stress changes in depletion zones induced by parent wells as a key factor influencing fracture propagation. This result is shown by new measurements of in-situ fracture propagation velocity and fracture-hit volume (fluid volume at fracture hit?) from LF-DAS and event density from microseismic. These findings highlight the importance of considering parent well effects, well spacing, and stimulation sequencing in completion design to improve reservoir development and production efficiency.

depletion zones↗
DOE OSTI2021

Developing Upscaling Approach for Swarming Hydraulic Fractures Observed at Hydraulic Fracturing Test Site through Multiscale Simulations

We show that this work aims to address a challenge posed by recent observations of tightly spaced hydraulic fractures in core samples from the hydraulic fracturing test site (HFTS) in the Middle Wolfcamp Formation. Many fractures in retrieved cores have subfoot spacing, which is at odds with conventional models in which usually one hydraulic fracture is initiated per cluster. Models assuming a single fracture at each cluster, although a common practice, often predict excessive fracture propagation that is inconsistent with microseismic observation. Here, we aim to develop a numerical approach to effectively account for densely spaced hydraulic fractures in field-scale simulations. Because it is impractical to explicitly model all aforementioned fractures, we develop a new upscaling law that enables existing simulation tools to predict reservoir response to fracture swarms. The upscaling law is derived based on an energy equivalence argument and validated through multiscale simulations using a high-fidelity code, GEOS. The swarming fractures are first modeled with a spacing that is much smaller than the cluster spacing; these fractures are then approximated by an upscaled, single fracture based on the proposed upscaling law. The upscaled fracture is shown to successfully match the energy input rate and produce the total fracture aperture and average propagation length of the explicitly simulated swarm. Afterward, the upscaling approach is further implemented in 3D field-scale simulations and validated against the HFTS microseismic data of a horizontal well. Our results show that hydraulic fracture swarming can significantly affect fracture propagation behaviors compared with the propagation of single fractures as assumed by conventional modeling approaches. Under the considered situations, the conventional treatment yields fast propagation speed that far exceeds that indicated by the microseismic data. We also illustrate that this discrepancy can be reduced readily through the implementation of the upscaling law. Our results demonstrate the importance of accounting for the fracture swarming effect in field-scale simulations and the efficacy of this approach to enable realistic predictions of reservoir responses to fracture swarms, without the need to model tightly spaced fractures individually.

02 PETROLEUM↗
DOE OSTI2020

Three-Dimensional Evaluation of Sand Particle Fracture Using Discrete-Element Method and Synchrotron Microcomputed Tomography Images

Recent research showed that fracture of sand particles plays a significant role in determining the plastic bulk volumetric changes of granular materials under different loading conditions. One of the major tools used to better understand the influence of particle fracture on the behavior of granular materials is discrete-element modeling (DEM). This paper employed the bonded block model (BBM) to simulate the fracture behavior of sand. Each sand particle is modeled as an agglomerate of rigid blocks bonded at their contacts using the linear-parallel contact model, which can transmit both moment and force. DEM simulated particles closely matched the actual three-dimensional (3D) shape of sand particles acquired using high-resolution 3D synchrotron microcomputed tomography (SMT). Results from unconfined one-dimensional (1D) compression of a single synthetic silica cube were used to calibrate the model parameters. Particle fracture was investigated for specimens composed of three sand particles that were loaded under confined 1D compression. Breakage energy measured from DEM models matched well with that measured experimentally. The paper studied the effects of contact loading condition and particle interaction on the fracture mode of particles using BBM that can closely capture the 3D shape of real sand particles.

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 OSTI2022

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗
DOE OSTI2024

Quartz Dissolution Effects on Flow Channelization and Transport Behavior in Three‐Dimensional Fracture Networks

We perform a set of reactive transport simulations in three-dimensional fracture networks to characterize the impact of geochemical reactions on flow channelization. Flow channelization, a frequently observed phenomenon in porous and fractured subsurface rock formations, results from the spatially variable hydraulic resistance offered by a geological structure. In addition to geo-structural features such as network connectivity, geometry, and hydraulic resistance, geochemical reactions, for example, dissolution and precipitation, can dynamically inhibit or enhance flow channelization. These geochemical processes can change the fracture permeability leading to increased flow channelization, which are localized connected regions of high volumetric flow rates that are seemingly ubiquitous in the subsurface. In our simulations, fractures partially filled with quartz are gradually dissolved until quasi-steady state conditions are obtained. We compare the flow field's initial unreacted and final dissolved states in terms of flow and transport observations. We observe that the dissolved fracture networks provide less resistance to flow and exhibit increased flow channelization when compared to their unreacted counterparts. However, there is substantial variability in the magnitude of these changes which implies that the channelization strongly depends on the network structure. In turn, we identify the interplay between the particular network structure and the impact of geochemical dissolution on flow channelization. The presented results indicate that geological systems that have been weathering or reactive for longer times in older landscapes are likely to have increased flow channelization compared to their equivalent but younger counterparts, which implies a time dependence on flow channelization in fractured media.

Hyman, Jeffrey D.↗
DOE OSTI2022

Experimentally validated multiphysics modeling of fracture induced by thermal shocks in sintered UO 2 pellets

Uranium Dioxide (UO 2 ) fuel powers almost all commercial Nuclear Power Plants (NPPs) worldwide, generating carbon-free energy and contributing to the fight against climate change. UO 2 fuel incurs damage and fractures due to large thermal gradients that develop across the fuel pellet during normal and transient operating conditions. A comprehensive understanding of the underlying mechanisms by which these processes take place is still lacking. A combined experimental and computational approach is utilized here to quantify the behavior of UO 2 fuel fracture induced by thermal shock. Here, this work introduces both (1) an experimental study to understand the fuel fracturing behavior of sintered UO 2 pellets when exposed to thermal shock, and (2) a Multiphysics phase-field fracture model capable of simulating this process. Parametric studies were conducted to evaluate the effects of uncertainties in fracture properties on the fracture behavior of UO 2 due to thermal shocking. A set of energy release rate (or equivalently fracture toughness) and contract area (the part of the fuel pellet in direct contact with the cold bath) were able to capture the overall fracture trends of the corresponding experimental data. Our combined approach presents a new method for accounting for the effects of microstructure and sample size on the energy release rate/fracture toughness. The experimental data were collected from multiple experiments that exposed UO 2 pellets to high-temperature conditions (589–676 °C) followed by a quench in sub-zero water. This work demonstrates that joint experimental and computational efforts are able to advance the understanding of thermal fracture in the primary fuel source for existing and future NPPs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Controlling Sustainability of Hydraulic Fracture Permeability in Ductile Shales

Hydraulic fracturing technology, along with horizontal drilling and associated completion technologies, dramatically increased US oil and gas production, by increasing the permeability of tight shales and allowing for cost-effective extraction of trapped hydrocarbons. For efficient and economical production, the targeted shale has to be fractured easily (i.e., good ‘frackability’), but also the permeability of induced fractures have to be sustained during production (good ‘sustainability’). For these reasons, currently, ductile shale with high clay content (>~40%) is difficult to exploit as a hydrocarbon resource, although hydrocarbons can still be found in it. Good frackability, high-TOC, low-clay-content resource shales—the ‘low-hanging fruit’—will be depleted, and we need to develop tools that allow tapping into currently underutilized, clay-rich ductile shales. This research project aims (1) to understand the behavior of fractures in clay-rich, ductile (and sometimes swelling) shales and (2) to begin to develop technologies for efficient and economical production from such shales. With these objectives, we examine time-dependent, coupled mechanical-hydrological behavior of open and proppant-filled fractures within different types of shales. In the preceding project, the focus of the research was to understand the behavior of the fractures and to predict their permeability changes due to fracture closure from shale matrix deformation and proppant embedment. In the current project, we investigate the possibility of chemically manipulating the sustainability of hydraulic fractures in ductile shales—specifically by altering the proppant-embedding behavior—via chemical means, at different stages of oil and gas production.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2020

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗
DOE OSTI2020

Deterministic Discrete Fracture Network (DFN) Model for the EGS Collab Project on the 4850 Level of the Sanford Underground Research Facility (SURF)

The EGS Collab is conducting hydraulic fracture stimulation and fluid circulation experiments in the Sanford Underground Research Facility (SURF) located in Lead, South Dakota. A total of eight ~60m-long subhorizontal boreholes were drilled from the 4850 Level (~1.5 km below the ground surface) into the crystalline rock of this former mine. Six of these holes are used for geophysical monitoring, one is used for hydraulic fracture stimulation, and the remaining hole was designed as a production borehole that receives water from the injection well via the induced and natural fracture system. The primary goal of creating the discrete fracture network model is to show that these modeling methods are critical for the development of enhanced geothermal systems (EGS). This includes the prediction of rock behavior during fracturing and during an extended period of water flow between the parallel injection and production boreholes. Understanding the results from the induced fracturing and flow is complicated by the presence of significant natural fractures that interact with the stimulation and/or flow pathways. The delineation and characterization of natural fractures is thus an important part of the project, and therefore a model of the Discrete Fracture Network (DFN) was developed on a deterministic basis. The DFN was populated using observations and interpretations integrated from drift (horizontal passageways that allow access in the underground) fracture mapping, analysis of core recovered from the eight boreholes, borehole televiewer logs and videos, and observations of flow between and within boreholes and in the drift. The natural fracture system is dominated by a pervasive northwest-trending, steeply dipping shear system that is identifiable in the drifts and the core. Hydraulic fracture stimulation, flow/tracer circulation tests, and geophysical monitoring revealed that the behavior of the injected water, and perhaps the growth of induced fractures, has been significantly influenced by the existing fractures identified in the DFN.

40 EE - Geothermal Technologies Office (EE-4G)↗
DOE OSTI2022

Investigation of fracture in porous materials: a phase-field fracture study informed by ReaxFF

Microscopic features (e.g, pore shapes, sizes, and distribution) in porous material substantially affect the overall mechanical properties such as stiffness and strength. In turn, these material properties determine the macroscopic behaviors of fracture in the porous material. In certain cases, macroscopic properties can be derived from the porous skeleton and void ratio (i.e., porosity), but in many other cases, such derivation is a challenging task. This paper presents a numerical investigation of microporosity and micropore shapes effect on the macrofracture behavior in porous amorphous silica. For this study, we extend the recently-proposed combined molecular dynamic (MD) and phase-field (PF) fracture modeling approach by including different pore shapes in the atomistic domain. In the MD simulations, we adopt ReaxFF to evaluate the material properties, where four different micropore cases are considered. Based on the material properties derived from MD simulations, the macrofracture propagation of porous media is studied using hybrid PF simulation. In the characterization of the pore structure, the concept of pore ligament is proposed to relate the pore shape and the critical energy release rate. Two classical fracture problems were used to evaluate the effect of pore shape on the macrofracture behavior. The results of the case study show that although the micropore shapes change the macrofracture behaviors, these effects vary with the geometry and loading conditions of macroscopic boundary value problems. Furthermore, the case study also shows that the influence of micropore structure can be captured at the macroscopic level through the material properties derived from the MD simulations.

42 ENGINEERING↗