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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Stress effects on flow and transport in three-dimensional fracture networks

We investigate the effects of various external stress regimes on fracture apertures, fluid flow, and solute transport in three-dimensional fracture networks. We use well-established geomechanics equations coupled with discrete fracture network modeling to characterize changes in primary flow paths within a complex network as a function of stress magnitude and orientation. These changes manifest in the alterations of the fluid flow field and are measured in terms of Eulerian and Lagrangian flow observables including solute transport, which is a key problem in many hydrologic applications. Changes in primary flow paths affect the solute transport in the network by promoting anomalously early arrival or long tailing behavior. However, early time arrival is not ubiquitous in anisotropically stressed networks and in most cases there is a delayed arrival of solute, which is attributed to (i) the presence of low-velocity zones normal to the flow direction, (ii) the angle between flow direction and major compressive principal stress directions, and (iii) changes in the primary flow paths (i.e., increases in tortuosity and active network structure). Overall, flows become more channelized in anisotropically stressed fracture networks than in unstressed and isotropically stressed networks.

58 GEOSCIENCES↗
DOE OSTI2021

Thermal Experiments for Fractured Rock Characterization: Theoretical Analysis and Inverse Modeling

Abstract Field‐scale properties of fractured rocks play a crucial role in many subsurface applications, yet methodologies for identification of the statistical parameters of a discrete fracture network (DFN) are scarce. We present an inversion technique to infer two such parameters, fracture density and fractal dimension, from cross‐borehole thermal experiments data. It is based on a particle‐based heat‐transfer model, whose evaluation is accelerated with a deep neural network (DNN) surrogate that is integrated into a grid search. The DNN is trained on a small number of the heat‐transfer model runs and predicts the cumulative density function of the thermal field. The latter is used to compute fine posterior distributions of the (to be estimated) parameters. Our synthetic experiments reveal that fracture density is well constrained by data, while fractal dimension is harder to determine. Adding nonuniform prior information related to the DFN connectivity improves the inference of this parameter.

Zhou, Zitong↗
DOE OSTI2025

Co‐Evolution of Specific Stiffness and Permeability of Rock Fractures Offset in Shear

Abstract Fractures and faults represent planes of weakness and compliance in rock masses that serve as focal points for both microearthquakes and fluid transport, with seismicity and permeability evolution closely linked. Contact stiffness is highly stress‐sensitive and directly influences permeability. We explore the co‐evolution of specific stiffness and permeability of rough fractures under normal stress and shear offset using numerical simulations. Individual rough fractures are represented by variable amplitude (Root mean square) and wavelength ( λ ) using a granular mechanics model. Contacting rough surfaces are mated, offset in shear, and then compacted in displacement mode. The compacting fractures generate stress‐dependent changes in contact porosity, which govern both permeability and stiffness evolution. We establish a universal dimensionless relationship linking specific stiffness and permeability that inherently incorporates the effects of surface roughness, shear offset, and microcracking. The observed cracking effect—where local stress redistribution and pressure‐driven microcrack propagation dynamically alter the aperture field—introduces a nonlinear permeability response at high stress. Increased roughness amplitude and larger shear offsets reduce stiffness while dampening permeability sensitivity to stress, demonstrating a strong interplay between surface texture and hydro‐mechanical behavior. While the model captures this behavior effectively, deviations emerge at very low porosities due to extreme aperture sensitivity in this limit.

He, Xinxin [Energy and Mineral Engineering Pennsyl↗
DOE OSTI2025

Predicting Flow in Fracture Networks With Quantum Algorithms

Uncertainty quantification plays a crucial role in the modeling of subsurface flow. For instance, uncertainties in the properties of geologic fracture networks significantly impact flow, requiring numerous simulations to accurately estimate quantities of interest. However, each simulation is computationally expensive because it requires solving a large linear system to capture features that involve both small and large fractures. An example is in percolation, where the interaction of many small fractures (which cumulatively can have a large surface area) with the rock matrix must be modeled precisely. Quantum computing is an emerging tool with the potential to address this issue. Quantum algorithms offer a significant speedup in solving linear systems, achieving efficiencies that are challenging to match with classical approaches. These classical approaches include direct solvers, such as LU decomposition, and iterative methods, notably preconditioned conjugate gradient, commonly used in subsurface modeling to solve large sparse systems. However, applying quantum algorithms to geologic fracture flow requires careful attention to algorithmic and problem-specific constraints to fully realize this quantum advantage. In this work we describe a quantum algorithm for generalized Monte Carlo applications with a quadratic speedup over the classical approaches which can be combined with the quantum speedup, currently under investigation, for solving quantum linear systems for subsurface flow. We show that for quantum algorithms the computational cost of estimating a quantity of interest for a statistical ensemble of networks is roughly the same as that of a single realization, essentially implying that one can get uncertainty quantification for free.

58 GEOSCIENCES↗
DOE OSTI2023

River thorium concentrations can record bedrock fracture processes including some triggered by distant seismic events

Fractures are integral to the hydrology and geochemistry of watersheds, but our understanding of fracture dynamics is very limited because of the challenge of monitoring the subsurface. Here we provide evidence that long-term, high-frequency measurements of the river concentration of the ultra-trace element thorium (Th) can provide a signature of bedrock fracture processes spanning neighboring watersheds in Colorado. River Th concentrations show abrupt (subdaily) excursions and biexponential decay with approximately 1-day and 1-week time constants, concentration patterns that are distinct from all other solutes except beryllium and arsenic. The patterns are uncorrelated with daily precipitation records or seasonal trends in atmospheric deposition. Groundwater Th analyses are consistent with bedrock release and dilution upon mixing with river water. Most Th excursions have no seismic signatures that are detectable 50 km from the site, suggesting the Th concentrations can reveal aseismic fracture or fault events. We find, however, a weak statistical correlation between Th and seismic motion caused by distant earthquakes, possibly the first chemical signature of dynamic earthquake triggering, a phenomenon previously identified only through geophysical methods.

58 GEOSCIENCES↗
DOE OSTI2021

StressNet - Deep learning to predict stress with fracture propagation in brittle materials

Abstract Catastrophic failure in brittle materials is often due to the rapid growth and coalescence of cracks aided by high internal stresses. Hence, accurate prediction of maximum internal stress is critical to predicting time to failure and improving the fracture resistance and reliability of materials. Existing high-fidelity methods, such as the Finite-Discrete Element Model (FDEM), are limited by their high computational cost. Therefore, to reduce computational cost while preserving accuracy, a deep learning model, StressNet, is proposed to predict the entire sequence of maximum internal stress based on fracture propagation and the initial stress data. More specifically, the Temporal Independent Convolutional Neural Network (TI-CNN) is designed to capture the spatial features of fractures like fracture path and spall regions, and the Bidirectional Long Short-term Memory (Bi-LSTM) Network is adapted to capture the temporal features. By fusing these features, the evolution in time of the maximum internal stress can be accurately predicted. Moreover, an adaptive loss function is designed by dynamically integrating the Mean Squared Error (MSE) and the Mean Absolute Percentage Error (MAPE), to reflect the fluctuations in maximum internal stress. After training, the proposed model is able to compute accurate multi-step predictions of maximum internal stress in approximately 20 seconds, as compared to the FDEM run time of 4 h, with an average MAPE of 2% relative to test data.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Material Fracturing and Failure Simulation Datasets

Fracturing is a fundamental physics phenomena with broad relevance across multiple domains, ranging from infrastructure integrity, aerospace durability, reservoir production, and seismic events. We present a diverse dataset of simulated fracture evolution and material failure generated from two numerical solvers: the phase-field method and the combined finite-discrete element method (FDEM). These solvers differ in formulation, physical fidelity, and computational efficiency. The dataset includes five materials: PBX, anisotropic shale, tungsten, aluminum, and steel. For each, phase-field simulations span 400,000 cases: 200,000 under uniaxial tension and 200,000 under biaxial tension. The computationally expensive FDEM simulations include 90,000 split evenly among PBX, shale, and tungsten under uniaxial loading. All simulations begin with randomized initial fracture patterns. Each entry includes temporal data capturing fracture propagation dynamics. This comprehensive dataset is designed to support the development of foundational or surrogate machine learning approaches for predicting material failure. While no such models are introduced here, the dataset lays a robust foundation for advancing future research and innovation in these areas.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mineral Fabric as a Hidden Variable in Fracture Formation in Layered Media

Two longstanding goals in subsurface science are to induce fractures with a desired geometry and to adaptively control the interstitial geometry of existing fractures in response to changing subsurface conditions. Here, we demonstrate that microscopic mineral fabric and structure interact with macroscopic strain fields to generate emergent meso-scale geometries of induced fractures. These geometries define preferential directions of flow. Using additively manufactured rock, we demonstrate that highly conductive flow paths can be formed in tensile fractures by creating corrugated surfaces. Generation, suppression and enhancement of corrugations depend on the relative orientation between mineral fabric and layering. These insights into the role of micro-scale structure on macro-scale flow provide a new method for designing subsurface strategies to maximize potential production or to inhibit flow.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Incorporating a naphthalene diimide polymer into a fullerene electron-transport layer to improve the fracture energy of perovskite solar cells

By blending a naphthalene diimide polymer into C 60 , we made a solution-processed electron-transport layer (ETL) for perovskite solar cells with fracture energies of 1.25 J m −2 , over 3× higher than that of thermally evaporated C 60 . Fracture energies were measured in a double cantilever beam configuration, and fracture surface images showed a fracture location near the ETL/perovskite interface, indicating a toughening of the interface between the ETL and Ag. We show that this modification to the ETL has no adverse effect on solar cell performance, and highlight the additional benefit of reduced parasitic absorption; a finding relevant for tandem solar cells.

electron transport layer↗
DOE OSTI2024

An integrated experimental–modeling approach to identify key processes for carbon mineralization in fractured mafic and ultramafic rocks

Abstract Controlling atmospheric warming requires immediate reduction of carbon dioxide (CO2) emissions, as well as the active removal and sequestration of CO2 from current point sources. One promising proposed strategy to reduce atmospheric CO2 levels is geologic carbon sequestration (GCS), where CO2 is injected into the subsurface and reacts with the formation to precipitate carbonate minerals. Rapid mineralization has recently been reported for field tests in mafic and ultramafic rocks. However, unlike saline aquifers and depleted oil and gas reservoirs historically considered for GCS, these formations can have extremely low porosities and permeabilities, limiting storage volumes and reactive mineral surfaces to the preexisting fracture network. As a result, coupling between geochemical interactions and the fracture network evolution is a critical component of long-term, sustainable carbon storage. In this paper, we summarize recent advances in integrating experimental and modeling approaches to determine the first-order processes for carbon mineralization in a fractured mafic/ultramafic rock system. We observe the critical role of fracture aperture, flow, and surface characteristics in controlling the quantity, identity, and morphology of secondary precipitates and present where the influence of these factors can be reflected in newly developed thermo-hydro-mechanical–chemical models. Our findings provide a roadmap for future work on carbon mineralization, as we present the most important system components and key challenges that we are overcoming to enable GCS in mafic and ultramafic rocks.

58 GEOSCIENCES↗
DOE OSTI2022

Fracturing in Wet Granular Media Illuminated by Photoporomechanics

Here, we study fluid-induced deformation and fracture of granular media and apply photoporomechanics to uncover the underpinning grain-scale mechanics. We fabricate spherical photoelastic particles of 2-mm diameter to form a monolayer granular pack in a circular Hele-Shaw cell that is initially filled with a viscous fluid. The key distinct feature of our system is that, with spherical particles, the granular pack has a connected pore space, thus allowing for pore-pressure diffusion and the study of effective stress in coupled poromechanical processes. We inject air into the fluid-filled photoelastic granular pack, varying the initial packing density and confining weight. With our recently developed experimental technique, photoporomechanics, we find two different modes of fluid invasion: fracturing in fluid-filled elastic media (with strong photoelastic response) and viscous fingering in frictional fluids (with weak or negligible photoelastic response). We directly visualize the evolving effective stress field and discover an effective stress shadow behind the propagating fracture tips, where the granular pack exhibits undrained behavior. We conceptualize the behavior of the system by means of a mechanistic model for a wedge of the granular pack bounded by two growing fractures. The model captures the pore-pressure build-up inside the stress shadow region and the grain compaction in the annular region outside. Our model reveals that a jamming transition determines the distinct rheological behavior of the wet granular pack, from a friction-dominated to an elasticity-dominated response.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Fracture Failure Analysis for MQXFA Magnet Aluminum Shells

The High-Luminosity LHC Accelerator Upgrade Project (HL-LHC AUP) is approaching the production phase of the US-contributed Q1 and Q3 Interaction Region Quadrupoles (MQXFA). The structures for the MQXFA prototypes were design and inspected by the US-LARP (LHC Accelerator Research Program), AUP developed criteria, which will be used for the pre-series structures. As the first two full-length prototypes with 4.2 m magnetic length, MQXFAP1 and MQXFAP2, were designed and assembled at Lawrence Berkeley National Laboratory (LBNL), and tested at Brookhaven National Laboratory (BNL). The end aluminum short shell of MQXFAP2 was fractured along the shell length during the test, and tests were stopped. Analytical and Finite Element analysis were performed in light of the graded procedure defined in the Structure Design Criteria to investigate the fracture failure for MQXFAP2. In this paper, we report the fracture analysis of the current shell design, including the elasto-plastic simulations with submodel technique, and calculations with Linear Elastic Fracture Mechanics (LEFM). Test material properties are also presented. The results of this analysis explain why the end shell of MQXFAP2 failed, and suggest fillets on the end shell notches to meet the margin specified in the Structural Design Criteria.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

“Figure‐of‐Eight” Suture‐Button Technique for Fixation of Displaced Anterior Cruciate Ligament Avulsion Fracture

Objective To assess the clinical results of the “figure‐of‐eight” suture‐button technique in the arthroscopic treatment of anterior cruciate ligament (ACL) tibial avulsion fractures. Methods This was a retrospective study reviewing data from September 2013 to June 2019. A total of 27 patients (13 males and 14 females) who underwent arthroscopic “figure‐of‐eight” suture‐button fixation for displaced ACL avulsion fractures were analyzed. The mean age of the patients in the sample was 15.8 years (10–29 years), with a mean follow‐up of 24 months (6–48 months). According to Meyers–McKeever classification, 11 patients were classified as type III and 16 as type IV. All patients were evaluated following the guidelines of the radiological union, the Lysholm knee scoring scale, and the International Knee Documentation Committee (IKDC). Results Fractures were united within 3 months after surgery in all 27 cases. During the last follow‐up, all the anterior drawer and Lachman tests were negative, except in 1 patient whose anterior drawer test was 1° positive. The range of motion was improved from 72.22° ± 27.92° before surgery to 137.78° ± 7.38° at the last follow‐up ( P < 0.05); the Lysholm score was improved from 45.81 ± 10.94 before surgery to 93.04 ± 5.66 at the last follow‐up ( P < 0.05); and the IKDC score was increased from 43.89 ± 11.16 before surgery to 90.26 ± 5.86 at the last follow‐up ( P < 0.05). In 1 patient, an inflammatory reaction was observed at the medial incision of the tibial tubercle; the symptoms disappeared with administration of antibiotics for 1 week. All patients returned to their preinjury physical activities at the last follow‐up. Conclusion The “figure‐of‐eight” suture‐button technique achieves a satisfactory clinical outcome and provides an effective method for the treatment of displaced ACL avulsion fractures.

Kuang, Shi‐da↗
DOE OSTI2022

Exceptional fracture toughness of CrCoNi-based medium- and high-entropy alloys at 20 kelvin

Medium- and high-entropy alloys based on the CrCoNi-system have been shown to display outstanding strength, tensile ductility and fracture toughness (damage-tolerance properties), especially at cryogenic temperatures. Here we examine the J Ic and (back-calculated) K JIc fracture toughness values of the face-centered cubic, equiatomic CrCoNi and CrMnFeCoNi alloys at 20 K. At flow stress values of ~1.5 GPa, crack-initiation K JIc toughnesses were found to be exceptionally high, respectively 235 and 415 MPa√m for CrMnFeCoNi and CrCoNi, with the latter displaying a crack-growth toughness K ss exceeding 540 MPa√m after 2.25 mm of stable cracking, which to our knowledge is the highest such value ever reported. Characterization of the crack-tip regions in CrCoNi by scanning electron and transmission electron microscopy reveal deformation structures at 20 K that are quite distinct from those at higher temperatures and involve heterogeneous nucleation, but restricted growth, of stacking faults and fine nanotwins, together with transformation to the hexagonal closed-packed phase. The coherent interfaces of these features can promote both the arrest and transmission of dislocations to generate respectively strength and ductility which strongly contributes to sustained strain hardening. Indeed, we believe that these nominally single-phase, concentrated solid-solution alloys develop their fracture resistance through a progressive synergy of deformation mechanisms, including dislocation glide, stacking-fault formation, nano-twinning and eventually in situ phase transformation, all of which serve to extend continuous strain hardening which simultaneously elevates strength and ductility (by delaying plastic instability), leading to truly exceptional resistance to fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessment of four strain energy decomposition methods for phase field fracture models using quasi-static and dynamic benchmark cases

Abstract Strain energy decomposition methods in phase field fracture models separate strain energy that contributes to fracture from that which does not. However, various decomposition methods have been proposed in the literature, and it can be difficult to determine an appropriate method for a given problem. The goal of this work is to facilitate the choice of strain decomposition method by assessing the performance of three existing methods (spectral decomposition of the stress or the strain and deviatoric decomposition of the strain) and one new method (deviatoric decomposition of the stress) with several benchmark problems. In each benchmark problem, we compare the performance of the four methods using both qualitative and quantitative metrics. In the first benchmark, we compare the predicted mechanical behavior of cracked material. We then use four quasi-static benchmark cases: a single edge notched tension test, a single edge notched shear test, a three-point bending test, and a L-shaped panel test. Finally, we use two dynamic benchmark cases: a dynamic tensile fracture test and a dynamic shear fracture test. All four methods perform well in tension, the two spectral methods perform better in compression and with mixed mode (though the stress spectral method performs the best), and all the methods show minor issues in at least one of the shear cases. In general, whether the strain or the stress is decomposed does not have a significant impact on the predicted behavior.

Zhang, Shuaifang↗
DOE OSTI2021

Understanding Discrete Fracture Networks Through Spectral Graph Theory

Discrete Fracture Network models (DFNs) are used to simulate fluid flow and particle transport through fracture networks in low permeability rock. Understanding these processes are essential in many subsurface applications, such as environmental restoration of contaminated fractured media, CO 2 sequestration, detection of low-level nuclear tests, and hydrocarbon extraction. Compared with other models, DFNs allow for incorporation of a wider range of network characteristics but have substantially greater computation cost. These networks can be represented with graphs, allowing the use of graph theory tools to study the networks. I used Python to simulate flow and transport on a range of DFNs and analyzed these networks using methods from network analysis and spectral graph theory. My purpose was to find ways to gain insight about flow and transport on DFNs using these graph representations, bypassing the computationally intensive meshing typically required. My work is still in progress, but I have discovered several interesting trends and patterns that I believe could be useful towards my goal. If I am able to bring these results to fruition, they will aid subsurface geologists in extracting flow and transport information about fracture networks more efficiently.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Fracture Toughness Characterization of Generation II FeCrAl Alloys after ~18 dpa Irradiation

FeCrAl alloys are promising candidate materials for the accident tolerant fuel (ATF) cladding applications due to their excellent corrosion resistance to the elevated temperature steam environment. Currently, the handbook on FeCrAl material properties contains only limited data regarding the fracture toughness properties of any FeCrAl alloy. This includes alloys currently under investigation within the Advanced Fuels Campaign (AFC) at Oak Ridge National Laboratory (ORNL). In this project, a series of irradiation capsules have been irradiated in the High Flux Isotope Reactor (HFIR) at ORNL with two Generation II FeCrAl candidate alloys, i.e., C06M and C36M, to assess the fracture response of these alloys after neutron irradiation. These alloys represent the “book-end” compositions for C26M, the alloy currently being developed as the leading candidate for LWR cladding. A total of six irradiation capsules were irradiated in HFIR at target temperatures of 200°C, 330°C, and 500°C up to target damage doses of 8 displacements per atom (dpa) and 16 dpa. These damage doses represent the expected middle and end of life damage levels for typical LWR cladding while the irradiation temperature regimes will provide insight into the role of varying microstructural features on the fracture toughness properties of neutron irradiated FeCrAl alloys. To date, irradiation of all capsules has been completed in HFIR. This report summarizes the latest results of microhardness and fracture toughness PIE for the 16 dpa capsules (FCAB2, FCAB4, and FCAB6), for which the measured irradiation conditions were: 204°C/17.6dpa, 343°C/18.3dpa, and 507°C/18.6dpa. The main conclusions of this study can be summarized as follows: 1) After the 204°C/17.6dpa irradiation, both C06M and C36M exhibited significant irradiation hardening and embrittlement 2) After the 343°C/18.3dpa irradiation, both C06M and C36M exhibited small irradiation hardening without irradiation embrittlement 3) After the 507°C/18.6dpa irradiation, both C06M and C36M exhibited irradiation softening without irradiation embrittlement 4) Comparing the microhardness and Master Curve reference temperature T 0q before and after neutron irradiation, we did not observe a linear correlation between the two parameters for both C06M and C36M. This should be mainly due to a flat response of the Master Curve reference temperature T 0q to the irradiations at 166-204°C and 315-343°C ranges 5) C06M showed a lower T 0q , meaning better toughness, than C36M at the unirradiated condition and such trend was kept even after neutron irradiation except for the 166-204°C irradiation where both materials had similar T 0q . 6) In terms of hardening and embrittlement, the irradiation effect on both C06M and C36M appeared to saturate after an irradiation dose of 7 dpa.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Integrated Multiphysics Modeling of Environmentally Assisted Brittle Fracture

Brittle materials, such as cement, compose major portions of built infrastructure and are vulnerable to degradation and fracture from chemo-mechanical effects. Currently, methods of modeling infrastructure do not account for the presence of a reactive environment, such as water, on the acceleration of failure. Here, we have developed methodologies and models of concrete and cement fracture that account for varying material properties, such as strength, shrinkage, and fracture toughness due to degradation or hydration. The models have been incorporated into peridynamics, non-local continuum mechanics methodology, that can model intersecting and branching brittle fracture that occurs in multicomponent brittle materials, such as concrete. Through development of new peridynamic capabilities, decalcification of cement and differential shrinkage in clay-cement composites have been evaluated, along with exemplar problems in nuclear waste cannisters and wellbores. We have developed methods to simulate multiphase phenomena in cement and cement-composite materials for energy and infrastructure applications.

36 MATERIALS SCIENCE↗