Search NASA⌕ Search

SEARCH · Search NASA

Results for “FRACTURE MECHANICS”

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

A coupled thermo-hydro-mechanical model for simulating leakoff-dominated hydraulic fracturing with application to geologic carbon storage

A potential risk of injecting CO2 into storage reservoirs with marginal permeability (≲ 10 mD (1 mD = 10 –15 m 2 )) is that commercial injection rates could induce fracturing of the reservoir and/or the caprock. Such fracturing is essentially fluid-driven fracturing in the leakoff-dominated regime. Recent studies suggested that fracturing, if contained within the lower portion of the caprock complex, could substantially improve the injectivity without compromising the overall seal integrity. Modeling this phenomenon entails complex coupled interactions among the fluids, the fracture, the reservoir, and the caprock. Here, we develop a simple method to capture all these interplays in high fidelity by sequentially coupling a hydraulic fracturing module with a coupled thermal-hydrological-mechanical (THM) model for nonisothermal multiphase flow. The model was made numerically tractable by taking advantage of self-stabilizing features of leakoff-dominated fracturing. The model is validated against the PKN solution in the leakoff-dominated regime. Moreover, we employ the model to study thermo-poromechanical responses of a fluid-driven fracture in a field-scale carbon storage reservoir that is loosely based on the In Salah project's Krechba reservoir. The model reveals complex yet intriguing behaviors of the reservoir-caprock-fluid system with fracturing induced by cold CO 2 injection. We also study the effects of the in situ stress contrast between the reservoir and caprock and thermal contraction on the vertical containment of the fracture. The proposed model proves effective in simulating practical problems on length and time scales relevant to geological carbon storage.

58 GEOSCIENCES↗
DOE OSTI2021

Hydro-mechanical measurements and x-ray imaging of sheared crystalline rock fractures from EGS Collab experiments 1 & 2

Measuring hydro-mechanical properties of natural fractures is a prerequisite for optimizing hydraulic stimulation design and well placement. We completed experiments to characterize shear on natural fractures in schist, amphibolite, and rhyolite specimens drilled from EGS Collab Project’s field sites at the Sanford Underground Research Facility (SURF) in South Dakota. A triaxial direct shear method and coupled x-ray imaging were used to perform hydroshearing and mechanical shearing at the site’s in-situ stress conditions. This produced simultaneous measurements of fracture and matrix strength, permeability, stress-dependent aperture, dilation, and friction strength. Our results identified that only a subset of the natural fractures was weak enough for hydroshearing. Generally, hydroshearing increases fracture permeability by a factor of 10 or more and the enhancement is retainable over time. However, the shear slip does not always result in permeability enhancement. High content of phyllosilicates was found to associate with exceptionally weak fractures that also exhibited poor or even negative enhancement after stimulation. Combining our measurements with site data, we can predict that most observable fractures at the two EGS Collab sites do not meet the criteria for hydroshearing before tensile opening. In some cases, the visible fractures are low permeability and as strong as the adjacent rock. To induce hydroshearing before tensile opening, injection must target known weak and favorably oriented fractures with confirmed pre-existing permeability.

58 GEOSCIENCES↗
DOE OSTI2023

The effects of carbonaceous inclusions and their distributions on dynamic failure processes in boron carbide ceramics

Commercially available boron carbide ceramics typically have heterogeneous microstructures that contain distributions of processing-induced inclusions. The inclusions that are rich in carbon (i.e., carbonaceous) govern the underlying mechanisms of brittle fracture through wing crack formation, and thus dictate the mechanical response of the ceramic. In this study, we investigate the dynamic failure of five boron carbide ceramic materials with different inclusion populations. All of the materials were prepared by hot-pressing; four of these boron carbides contained different sizes and concentrations of carbonaceous inclusions, while one contained no carbonaceous inclusions. The heterogeneity distributions were characterized in some detail for statistical analysis using scanning electron microscopy and quantitative image analysis. A modified compression Kolsky bar setup with in situ ultra-high-speed microscopic imaging (10 million frames per second) was then used to study the influence of the inclusion distributions on the dynamic failure processes in these materials, at nominal high strain rates of 10 2 -10 3 s -1 . The in situ ultra-high-speed microscopy highlighted the link between micro and macroscale failure processes and demonstrated that the carbonaceous inclusions are indeed the preferential sites for nucleation of wing cracks, as previously hypothesized based on post-mortem observations. The relative orientation of an inclusion with respect to the compression axis was shown to affect the likelihood that it would participate in crack nucleation. In conclusion, all of the ceramics were also found to have orientation-dependent peak compressive stress, regardless of the presence of carbonaceous inclusions, suggesting that grain orientation distributions are also important.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Review—The Lithiation/Delithiation Behavior of Si-Based Electrodes: A Connection between Electrochemistry and Mechanics

Silicon is a promising alternative anode material to graphite because of its high gravimetric and volumetric energy densities. However, severe capacity fading is observed in Si electrodes, and it is a result of mechanical changes of Si, such as volume changes, stress or fracture. Furthermore, these mechanical behaviors are strongly coupled with the electrochemistry of the Li-Si alloying reaction in Si-based electrodes, including both thermodynamics and kinetics. Therefore, the electrochemical properties of Si-based electrodes are strongly dependent on the control of the mechanics of Si during lithiation/delithiation. As such, it is very important to understand the correlation between electrochemistry and mechanics. Here, we review lithiation/delithiation behaviors of various types of Si-based electrodes, applying a fundamental understanding of electrochemistry and mechanics and the correlation between them.

25 ENERGY STORAGE↗
DOE OSTI2024

Impacts of irregularly-distributed acidified brine flow on geo-chemo-mechanical alteration in an artificial shale fracture under differential stress

The efficacy of geological carbon sequestration is reliant on the integrity of the caprock and its resistance to physical and chemical alteration. Caprocks with high abundance of reactive carbonates like calcite are susceptible to acid-promoted dissolution and can result in structural weakening. This work investigates the effect of acidified brine flow through an artificially fractured, high-carbonate (30 % by XRD) shale under differential compressive stress. Cylindrical samples were cut in half vertically and milled to create an artificial fracture with interlocking asperities and open channels. Samples were sheared with a single applied stress in a custom flow cell housed within an industrial CT scanner. Further, either acidic (pH 4) or reservoir-simulated (pH 9.5) brine was flowed through the artificial fracture for 7–8 days under reservoir pressure and room temperature. Model simulations indicate flow mainly occurred in open channels, with limited flow between overlapping asperities. Analysis of fracture surfaces by optical and scanning electron microscopy show increased surface alteration and roughness after exposure to pH 4 versus pH 9.5 brine indicating mineral dissolution/loss, and this effect is greater in areas that receive the highest brine flows. Similarly, surface analysis by scratch testing shows fracture toughness decreases more after exposure to acidic versus reservoir-simulated brine, with the greatest alteration in areas of highest acidic brine flows. Despite weakening, no shear slip occurred. Overall, the results indicate that acidified brine can result in significant physical and geomechanical alteration of irregular fracture surfaces in shale caprock, with greatest effects in preferential flow regions.

58 GEOSCIENCES↗
DOE OSTI2026

In situ quantification of fracture slip induced by hydraulic injections in a deep borehole: A comparison of two different borehole techniques

In situ measurements of fracture deformation during fluid injection are rare, yet essential for understanding the mechanical response of fractured rock. In this study, we evaluate the reliability of two methods by comparing their slip vector estimates: high-resolution borehole acoustic televiewer images captured before and after injection tests, and displacement data from a three-component borehole deformation probe recorded during the injections. Acoustic televiewer images capture only final in-plane displacement, whereas three-component borehole deformation measurements provide full 3D, transient fracture movement. Four injection tests in a fractured granitic rock mass along an inclined borehole at the Bedretto Underground Laboratory (Switzerland), beneath more than 1100 m of overburden, were analyzed. The two methods yielded consistent kinematics and comparable slip magnitudes, typically in the range of 0.2 – 0.6 mm. Angular differences between estimated slip directions ranged from 0.2° to 28°. For two test intervals with pre-existing fractures optimally oriented for slip, the measured slip vectors agreed with forward-modeled slip vectors based on the regional in situ stress field, supporting the Wallace–Bott hypothesis. Measurable shear displacement was identified in two test intervals containing fractures oriented nearly perpendicular to the least principal stress, including a newly formed hydraulic fracture. This demonstrates that even fractures unfavorably oriented for slip can undergo shear due to local stress perturbations around an inclined borehole.

Acoustic televiewer↗
DOE OSTI2020

Two-Phase Fluid Flow Properties of Rough Fractures With Heterogeneous Wettability: Analysis With Lattice Boltzmann Simulations

Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.

58 GEOSCIENCES↗
DOE OSTI2025

Thermo-Mechanical Phase-Field Modeling of Fracture in High-Burnup UO2 Fuels Under Transient Conditions

This study presents a novel multiphysics phase-field fracture model to analyze high-burnup uranium dioxide (UO2) fuel behavior under transient reactor conditions. Fracture is treated as a stochastic phase transition, which inherently accounts for the random microstructural effects that lead to variations in the value of fracture strength. Moreover, the model takes into consideration the effects of temperature and burnup on thermal conductivity. Therefore, the model is able to predict crack initiation, propagation, and complex morphologies in response to thermal gradients and stress distributions. Several simulations were conducted to investigate the effects of operational and transient conditions on fracture behavior and the resulting cracking patterns. High-burnup fuels exhibit reduced thermal conductivity, elevating temperature gradients and resulting in extensive radial and circumferential cracks. Transient heating rates and temperatures significantly affect fracture patterns, with higher heating rates generating steeper gradients and more irregular crack trajectories. This approach provides critical insights into fuel integrity during accident scenarios and supports the safety evaluation of extended burnup limits.

Chemistry↗
DOE OSTI2022

Dynamic aspects of graphene deformation and fracture from approximate density functional theory

Graphene is one of the most intriguing two-dimensional carbon materials. Its mechanical strength and failure are key concerns for materials engineering and applications. Despite the success of fracture mechanics, the mechanism of how pristine materials fail remains an elusive problem. While many theoretical studies based on molecular dynamics using empirical forcefields have tried to address this question, atomic-scale mechanics are not clearly understood. Especially, a widely employed bond-breaking approach based on the critical bond length has not been rigorously tested. Here, utilizing molecular dynamics simulations with density functional based tight binding, we investigate how the failure of the pristine material initiates. The Wiberg bond order (W BO ) to estimate the change of chemical bonds shows a transition from sp 2 (W BO ~ 1.33) to sp 3 (WBO < 1.0) during the deformation. However, it reveals that a single threshold value for either the WBO or bond length is insufficient to decide failure of pristine material without free edges or defects. Instead, collective behaviors of the local atomic group govern the fracture initiation of pristine graphene. Overall, our study provides dynamic mechanical responses based on quantum mechanics, which have not been captured by empirical forcefields, opening opportunities to design properties by precisely coupling the mechanics and quantum mechanics.

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

Tailoring the Toughening Effects in Two-Dimensional Nanomaterial-Reinforced Ceramic Matrix Composites

Abstract Ceramic matrix composites (CMCs) reinforced by two-dimensional (2D) nanomaterials have shown extraordinary load-carrying capacities, even in the harsh environments required by emerging applications. Their exceptional mechanical performance, especially fracture toughness, primarily arises from their heterogeneous microstructures. The deliberate dispersion of 2D reinforcements enables toughening mechanisms that are extrinsic to the matrix and thus endows the composites with substantial resistance to catastrophic failure. However, the incomplete understanding of the fracture behavior of such nanocomposites, especially the complex energy dissipation process of the matrix/reinforcement interface, limits the development of stronger and tougher CMCs. To overcome these limitations, we investigate crack deflection and energy dissipation in nanocomposites using an extended cohesive shear-lag model. This new model accounts for interfacial debonding and friction, which critically control the toughening of nanocomposites. Our analysis provides mechanistic insights for optimizing the toughening effects of CMCs.

Mechanics↗
DOE OSTI2021

Wide-blocky veins explained by dependency of crystal growth rate on fracture surface type: Insights from phase-field modeling

Vein microstructures contain a wealth of information on coupled chemical and mechanical processes of fracturing, fluid transport, and crystal growth. Numerical simulations have been used for exploring the factors controlling the development of vein microstructures; however, they have not been quantitatively validated against natural veins. Here we combined phase-field modeling with microtextural analysis of previously unexplained wide-blocky calcite veins in natural limestone and of the fresh fracture surface in this limestone. Results show that the wide-blocky vein textures can only be reproduced if ~10%–20% of crystals grow faster than the rest. This fraction corresponds to the amount of transgranularly broken grains that were observed on the experimental fracture surfaces, which are dominantly intergranular. We hypothesize that transgranular fractures allow faster growth of vein minerals due to the lack of clay coatings and other nucleation discontinuities that are common along intergranular cracks. Our simulation results show remarkable similarity to the natural veins and reproduce the nonlinear relationship between vein crystal width and vein aperture. This allows accurate simulations of crystal growth processes and related permeability evolution in fractured rocks.

Geology↗
DOE OSTI2020

Fluid-driven fractures in granular media: Insights from numerical investigations

We investigate the mechanisms of opening-mode fracture initiation in granular media. The study is based on a simulation of grain-scale fluid-grain interactions through a coupled numerical approach in which the discrete element method is used to solve for the mechanics of a solid granular medium, and computational fluid dynamics is used to model fluid flow and drag forces. In this work, we present benchmark problems with analytical solutions and validate this numerical model against experiments on a viscous-drag-driven cavity in the literature. Additional simulation results show fracture initiation mechanisms in a random granular packing subjected to constant boundary stresses and to fluid injection with a localized source. The dimensionless variable $F_s/F_{sk}$ (ratio of seepage force $F_s$ and skeletal force $F_{sk}$) incorporates the impacts of physical properties and injection parameters including fluid viscosity, injection velocity, grain size, and effective stresses, and it has been used as a criterion separating regimes of fluid invasion and drag-driven fracture opening. Our simulation results show that $F_s/F_{sk}$ in combination with $τ_1$ (ratio of diffusion time from hydromechanical coupling and injection time) serves as a prediction of fracture opening within granular packing. We suggest a simple criterion ($F_s/F_{sk}$ > 1 or $τ_1$ > 0.17) that is valid for various types of granular media and injection conditions to determine if fracture opening will occur. Among other applications, this study is useful to predict the initiation and propagation of fractures in natural sediments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Multilayer Architecture Barriers for High Temperature Nuclear Applications

Innovative thin film multilayer metal/ceramic barrier coatings have been developed by Argonne National Laboratory to protect fuel cladding against different phenomena that can be encountered under the harsh operating conditions of advanced nuclear reactors. One phenomenon is Fuel-Cladding Chemical Interaction (FCCI), which can represent a challenge to the performance of steel-clad metallic fuels in sodium-cooled fast reactors (SFRs), especially when operating to high burnups at elevated cladding temperatures, and long fuel residence time. Another phenomenon is the corrosion of zirconium-based cladding under prototypic boiling water reactor (BWR) conditions and associated cladding during off normal conditions, including loss of coolant accidents (LOCA). This work provides details of the development and testing of a thin film (few microns) multilayer metal/ceramic (CrY/Y 2 O 3 ) coatings on steel cladding and its performance against FCCI. The coating architecture evolved through systematic mechanical testing, radiation tolerance assessments, and high-temperature diffusion studies from a preliminary multilayer concept to an optimized high-performance design. The optimized structure leverages stoichiometric Y 2 O 3 for its exceptional radiation stability, chemical inertness, and ultra-low fission product diffusivity (e.g., Ce diffusivity ~10 -20 m 2 /s in irradiated Y 2 O 3 at 650 °C). This is paired with ductile chromium-rich interlayers (Cr 95 Y 5 ), which significantly enhance fracture toughness and mechanical resilience. Design refinements, including increased ceramic layer thickness and optimized Cr–Y interlayer composition, eliminated fracture propagation and preserved coating integrity under extreme mechanical stress (>10× typical reactor conditions) and prolonged radiation exposure (~300 dpa). Those conditions are well beyond expected conditions during normal operations in SFR, and cover some of the transient conditions as well. In addition, the coating deposition was demonstrated to be applied uniformly to the inner surface of steel tubes representing prototypic SFR cladding. Another phenomenon under consideration in this work is the corrosion of zirconium-based cladding under prototypic boiling water reactor (BWR) conditions and associated off normal conditions including loss of coolant accident (LOCA). Another multilayer thin film coating is considered for mitigation of this phenomenon, where a coating of yttrium aluminum and yttrium aluminum ceramic/metal thin film coating. Preliminary results on the coating performance during high temperature thermal cycling (up to 1000 C) have shown promising performance of the coating. On-going optimization of the coating will continue in FY26.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Peridynamics and surrogate modeling of pressure-driven well stimulation

In this work we use the peridynamics theory of solid mechanics to simulate fracture in an annular rock domain subject to an in-situ stress and create surrogate models that predict the area of the resulting cracks. Peridynamics is a non-local formulation of continuum mechanics that naturally accommodates material discontinuities. Furthermore, unlike other fracture modeling techniques there is no need to provide information about the crack path. We utilize the peridynamics code Peridigm and take a two-stage approach to fracture modeling. First an implicit solve is performed to compute the in-situ stress state. We then execute an explicit solve where a pressure loading designed to emulate fluid-driven hydraulic fracture is applied at the borehole and transmitted to the pre-stressed rock. We present results from polynomial and single and multi-level Gaussian process surrogate models constructed from a sampling study of the peridynamics model. The surrogates predict crack area given a measure of the in-situ stress anisotropy and rise time and amplitude of the pressure loading. These surrogates take a minuscule fraction of peridynamics model's running time to evaluate and are a step towards enabling advanced optimization and uncertainty quantification workflows that require many model evaluations.

42 ENGINEERING↗
DOE OSTI2025

Influence of Initial fabric and water-wetting on particle fracture and force chain evolution in natural silica sand

Particle fracture has a significant influence on the engineering behavior of granular materials. However, the combined influence of the initial fabric and wetting on particle fracture and force chain evolution in granular materials, particularly in sands, remains insufficiently explored in the current literature. Here, in this paper, one-dimensional (1D) confined compression experiments were conducted on dry and wet specimens, and particle-level fracture was captured using three-dimensional (3D) in-situ Synchrotron Micro-Computed Tomography (SMT). A quantitative assessment of particle fracture revealed that wet specimens, which have disturbed fabric, exhibited a higher fracture percentage in comparison to dry specimens. Complementary 3D finite element (FE) simulations were performed using dry and wet (only solid fabric was considered and the influence of pore water is not considered) sand assemblies to assess interparticle contact forces and particle-scale stress distributions within the specimens. Force chain analysis was subsequently conducted, encompassing quantification of interparticle forces, characterization of force network structures, and monitoring the dynamic evolution of force chains under different strain levels. The results show that the specimens with disturbed fabric led to a more dynamic and less persistent force network, more fabric instability, and thus more reorganization of force chain structures. In addition, the frequent rearrangement of the force network in the presence of water (with reduced inter-particle contact friction) likely exacerbates localized stress concentrations, promoting failure in previously unengaged or weakly connected particles. The results reported in this paper offer a new insight into how the initial fabric and wetting cause different fracture behavior. The findings can also pave the way for more in-depth future investigations into the mechanics governing particle fracture in granular assemblies.

Finite element analysis↗
DOE OSTI2022

Temperature-Dependent Fracture Resistance of Silicon Nanopillars during Electrochemical Lithiation

During the lithation of silicon anodes, the solid-state diffusion of lithium into Li x Si follows the Arrhenius law, the resulting morphology and fracture behavior are determined by the silicon anode operation temperature. Here, we reveal the temperature dependence of the lithiation mechanics of crystalline silicon nanopillars (SiNPs) via microscopic observations of the anisotropic growth and fracture behavior. In this work, we fabricated 1D SiNP structures with various orientations ($\langle100\rangle, \langle110\rangle$, and $\langle111\rangle$) as working electrodes and operated them at temperatures ranging from –20 to 40 °C. The lithiation of crystalline silicon at low temperatures exhibited preferential volume expansion along $\langle110\rangle$ and decreased fracture resistance. Furthermore, low temperatures caused the catastrophic fracture of amorphous silicon after the second lithiation. Our findings demonstrate the importance of silicon anode temperature control to prevent mechanical fracture during the cycle of lithium-ion batteries in harsh environments (e.g., electric vehicles in winter).

25 ENERGY STORAGE↗
DOE OSTI2025

Ductility mechanisms in complex concentrated refractory alloys from atomistic fracture simulations

The striking variation in damage tolerance among refractory complex concentrated alloys is examined through the analysis of atomistic fracture simulations, contrasting behavior in elemental Nb with that in brittle NbMoTaW and ductile Nb45Ta25Ti15Hf15. We employ machine-learning interatomic potentials (MLIPs), including a new MLIP developed for NbTaTiHf, in atomistic simulations of crack tip extension mechanisms based on analyses of atomistic fracture resistance curves. While the initial behavior of sharp cracks shows good correspondence with the Rice theory, fracture resistance curves reveal marked changes in fracture modes for the complex alloys as crack extension proceeds. In NbMoTaW, compositional complexity appears to promote dislocation nucleation relative to pure Nb, despite theoretical predictions that the alloy should be relatively more brittle. In Nb45Ta25Ti15Hf15, alloying alters the fracture mode compared to elemental Nb, promoting crack tip blunting and enhancing resistance to crack propagation.

Computational methods↗
DOE OSTI2020

Mechanical testing of PH13-8Mo H950 Steel for Xue-Wierzbicki Fracture Criterion Determination at 20° C and -40° C

Mechanical testing was conducted to collect the data needed to build a Xue-Wierzbicki (XW) fracture model for PH13-8 Mo H950 stainless steel (PH 13-8 SS). This model is intended for use in structural analysis of this material between room temperature and -40° C. Tests were performed on four different specimen geometries such that a range of stress states were characterized at room temperature and -40° C. Tensile tests on R5 tensile specimens were also performed to assess material anisotropy. Fracture toughness test were also conducted. The fracture toughness of this material at -40° C was 68% of the room-temperature value. Material strength generally increased with decreasing temperature while the opposite trend was observed for ductility. These trends were most pronounced for specimens with the largest stress triaxialities.

36 MATERIALS SCIENCE↗