Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fractures”

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 217 records · Page 12

Geologic stress modulates fluid mixing at fracture intersections

Fracture intersections are critical links that enable flow and transport in subsurface fracture networks, and their behavior strongly influences fluid mixing in a network. Although all subsurface fractures are subjected to geological stress, we lack a fundamental understanding of how fracture intersection geometry evolves under stress and how these changes influence fluid mixing. Here, we combine 3D printing, 3D X-ray tomographic imaging, and 3D pore-scale numerical simulations to reveal stress-induced changes in intersection geometry and their impact on mixing. Mixing is found to be strongly affected by partial closure of an intersection under stress. As an intersection closes, the void area for fluid flow and diffusion decreases leading to substantial deviations between conventional mixing models and full pore-scale modeling. To address this, we propose a modified mixing model that accounts for intersection deformation, which is essential for accurate modeling of solute transport and mixing through fracture networks.

15 GEOTHERMAL ENERGY↗

Dynamic Binary Complexes (DBC) as Super-Adjustable Viscosity Modifiers for Hydraulic Fracturing Fluids

In the preceding project year two, we refined three DBC formulations from a selection of over 50 different chemistries. The optimization study primarily encompassed testing for (i) reversibility extent, (ii) performance in the presence of chemical additives, (iii) adhesion and friction behavior during displacement in wellbores and pipelines, (iv) corrosion protection performance, and (v) injection performance with model fracture systems at the laboratory scale. Highly promising results obtained from all these tests signify the significant potential of DBCs in enhancing hydrocarbon recovery from unconventional reservoirs. The primary activities in the third project year included publishing experimental findings across multiple articles and conducting outreach initiatives. Throughout the year, we undertook tasks such as replicating experimental results, further optimizing various formulations and their associated experimental sets, and conducting additional tests to address missing components based on reviewer feedback and suggestions. We also explored the surfactant and friction-reduction aspects of selected formulations through drag reduction tests. In addition, we constructed an improved fracturing performance setup and performed flow injection tests. The specific DBC formulations focused on during this project period were A8/B1, A12/B5, and A10/B12. We also obtained results for additional DBC formulations and a select few commercial fracturing fluids for the purpose of comparison. Within the project's scope, we aim to enhance the experimental findings with the development of various models. The first two years focused on two key aspects: (i) the creation of a high-fidelity hydraulic fracturing model for non-Newtonian fluids to gain insights into the implementation of DBC fluids in fracking environments, and (ii) the development of a multiphase flow simulator for estimating total production, fluid saturation in the reservoir, and the creation of a fracture propagation model and kinetic Monte Carlo (kMC) models for diverse applications. In the third year, we delved into the fundamental nanostructural properties of DBCs, exploring aspects such as material chemistry, pH tunability, and control of DBC formation and stability. Subsequently, in the extension year, we conducted a systematic investigation of various building blocks containing primary, secondary, and tertiary amine functional groups to understand their impact on rheological and viscoelastic properties. Furthermore, we explored a Dissipative Particle Dynamics (DPD) model to simulate self-assembly processes with precision, creating a high-fidelity representation of relevant nanostructures. The tasks performed this year with the significant results obtained have been discussed in Section 2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

FracML: A Machine Learning Based Tool to Quantify Reservoir Scale Fracture Network for CO2 Storage

Poster on “FRACML: A Machine Learning Based Tool to Quantify Reservoir Scale Fracture Network for CO2 Storage” for the CCUS 2025 conference held in Houston, Texas March 3-5, 2025. The accurate characterization of subsurface fracture networks is essential for the secure operation of carbon capture, utilization, and storage (CCUS) projects. A thorough understanding of the spatial distribution of subsurface faults and fractures is crucial for predicting CO2 plume evolution and minimizing risks such as potential leakage into overlying formations or induced seismicity. In this context, robust fracture network quantification plays a pivotal role in reservoir management, providing the data necessary to fine-tune operational parameters, and ensure the environmental and economic viability of CCUS projects. As part of the U.S. Department of Energy’s SMART (Science-informed Machine Learning for Accelerating Real-time Decisions in Subsurface Applications) initiative, we focused on the development and application of a machine learning-based tool (FRACML) designed to quantify and map fracture networks using real-world (non-synthetic) data from an active CO2 injection site. Our objective is to demonstrate the utility of this tool in improving operational efficiency and safety across CCUS sites.

artifical intelligence / machine learning (AI/ML)↗

Fracture resistance of a TiB2 particle/SiC matrix composite at elevated temperature

The fracture resistance of a commercial TiB 2 particle/SiC matrix composite was evaluated at temperatures ranging from 20 to 1400 °C. A laser interferometric strain gauge (LISG) was used to continuously monitor the crack mouth opening displacement (CMOD) of the chevron-notched and straight-notched, three-point bend specimens used. Crack growth resistance curves (R-curves) were determined from the load versus displacement curves and displacement calibrations. Fracture toughness, work-of-fracture, and R-curve levels were found to decrease with increasing temperature. Microstructure, fracture surface, and oxidation coat were examined to explain the fracture behavior.

Fracture↗

A review of thermo-hydro-mechanical modeling of coupled processes in fractured rock: From continuum to discontinuum perspective

Coupled thermo-hydro-mechanical (THM) processes in fractured rock are playing a crucial role in geoscience and geoengineering applications. Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling. This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions. By comparing various models, one can better assess the uncertainties in predictions, particularly those related to the conceptual models. The review explores how these methodologies have significantly enhanced the fundamental understanding of how fractures respond to fluid injection and production, and improved predictive capabilities pertaining to coupled processes within fractured systems. It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience. The selection and calibration of model parameters should be based on typical ranges and applied to the specific conditions of applications. The challenges arising from inherent heterogeneity and uncertainties, nonlinear THM coupled processes, scale dependence, and computational limitations in representing field scale fractures are discussed. Realizing potential advances on computational capacity calls for methodical conceptualization, mathematical modeling, selection of numerical solution strategies, implementation, and calibration to foster simulation outcomes that intricately reflect the nuanced complexities of geological phenomena. Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical field of study.

Coupling scheme↗

3D probabilistic fracture mechanics / computational fluid dynamics simulation of a reactor pressure vessel under transient conditions

Reactor pressure vessels (RPVs) are safety-critical light-water-reactor components that, under irradiation, experience long-term material degradation in the form of embrittlement. This can increase their susceptibility to fracture under thermal-shock conditions, which could occur during off-normal transients such as loss-of-coolant accidents (LOCAs). During a LOCA, the most severe conditions for the RPV occur when emergency core cooling water is injected through the cold legs into the water-and-steam-filled RPV. The rapid cooling of the downcomer and internal RPV surface causes decreased temperature and elevated thermally driven tensile stresses in the RPV wall. This, combined with long-term material embrittlement, may cause fracture initiation at pre-existing flaws, challenging the integrity of the RPV. Assessing RPV integrity during transients with a large spatial variation in the coolant temperature requires a modeling approach that considers the effects of spatially varying coolant temperature on the fracture probability of a population of flaws distributed throughout the RPV, accounting for spatially varying embrittlement. Here, the present study addresses this need by demonstrating first-of-its-kind coupling of a high-fidelity 3D computational fluid dynamics code with 3D probabilistic fracture mechanics This was accomplished using representative models of a pressurized-water reactor subjected to small- and medium-break LOCA conditions, both of which can result in large spatial temperature variations. While the observed impact of accounting for 3D effects was minimal under the small-break LOCA this study indicates a significant increase in the probability of fracture initiation under the medium-break LOCA when 3D effects are considered, relative to a spatially uniform cooling scenario.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

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↗

Phase instability-coupled fracture behavior in garnet LLZO solid electrolytes: a machine learning-enabled atomistic study

Fracture in the garnet-type solid electrolyte Li 7 La 3 Zr 2 O 12 (LLZO) poses a critical threat to both the performance and safety of solid-state batteries. To unravel the coupled chemomechanical processes that govern fracture in LLZO under external loading, we carry out large-scale molecular-dynamics simulations with a validated machine-learning force field. Our results demonstrate that triaxial stresses at crack flanks trigger a localized cubic-to-tetragonal phase transformation, which is accompanied by Li-ion rearrangement. The emergent tetragonal domains feature lattice contraction normal to the fracture plane, imposing coherent misfit strains that provide an additional driving force for further crack propagation. Crucially, introducing Li deficiencies stabilizes the cubic phase, postponing the phase transition and thereby delaying fracture initiation. These findings highlight the role of intrinsic phase instability in dictating LLZO's fracture resistance and its critical connection to local Li concentration. This chemomechanical coupling points toward targeted strategies to enhance the mechanical robustness of garnet electrolytes, including tuning Li content, ensuring dopant homogeneity, and refining processing protocols.

Materials science↗

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗

Demonstration of thermal fracture limits of gas-cooled Nd:APG-1 and Tm:YLF slabs

High-energy short-pulse laser systems demonstrated to date have been limited to low repetition rates due to significant thermal management challenges associated with scaling these systems to high average power, including the need to understand the thermal fracture limits of laser gain media. In this paper, we demonstrate the thermal fracture limits of two high-average-power gain media in a realistic diode-pumped amplifier geometry: a 5 mm thick Nd:APG-1 slab and a 2 mm thick Tm:YLF slab, both cooled by room-temperature helium gas. Finite element analysis revealed maximum tensile stress at fracture of 25 ± 3 MPa for Nd:APG-1 and 30 ± 3 MPa for Tm:YLF. These stress levels correspond to average volumetric heat loads of 13 ± 1 W/cm 3 in Nd:APG-1 and 247 ± 12 W/cm 3 in Tm:YLF. Additional simulations show geometric scaling: thinner slabs tolerate higher heat loads before fracture, while larger pump spots reach the fracture stress at lower heat loads.

Lasers↗

Engineered Water for Improved Oil Recovery from Fractured Reservoirs (Final Scientific/Technical Report)

The goal of this proposal is to field test the novel technology of engineered water (modifying ionic composition and/or adding surfactants or nanoparticles) injection to improve oil recovery from the carbonate matrix in fractured reservoirs. Carbonate reservoirs tend to be oil-wet/mixed-wet due to positively charged carbonate mineral surfaces and negatively charged acidic/asphaltic components in the oil. Many carbonate reservoirs are also highly fractured. Water injection does not work efficiently in fractured oil-wet reservoirs because water does not imbibe into the matrix. The water composition can be engineered with ionic modification, surfactants and nanoparticles which can change the wettability of carbonate rocks towards more water-wet conditions which can imbibe water into the matrix and push the oil into the fractures. The engineered water can also lower interfacial tension and minimize capillary forces so that gravity can push the oil up and out of the matrix. Wettability alteration also increases water relative permeability which can also improve oil recovery in non-fractured carbonate reservoirs. This technology, if successful, can be a game changer for West Texas carbonates.

02 PETROLEUM↗

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↗

Final Report for Closing the Loop Between In Situ Stress Complexity and EGS Fracture Complexity

Characterizing in situ stress is essential in Enhanced Geothermal System (EGS) for both risk assessment and operation design. Standard stress measurement techniques such as diagnostic fracture injection tests (DFITs) require creating hydraulic fractures and apply analytical methods such as the G -function analysis to estimate the stress. These classical approaches rely on several assumptions, including but not limited to: (i) the induced fracture is planar and persistent, and (ii) the fracture plane is normal to the minimum horizontal stress. However, these assumptions may not hold in EGS environments like Utah FORGE, where thermo-hydro-mechanical (THM) coupling, wellbore deviation, and rock heterogeneity can alter the local stress condition and consequently lead to fracture complexities in both near-well and far-field regions.

15 GEOTHERMAL ENERGY↗

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↗

Investigation of flaw geometry and loading effects on plane strain fracture in metallic structures

The effects on fracture and flaw growth of weld-induced residual stresses, combined bending and tension stresses, and stress fields adjacent to circular holes in 2219-T87 aluminum and 5AI-2.5Sn(ELI) titanium alloys were evaluated. Static fracture tests were conducted in liquid nitrogen; fatigue tests were performed in room air, liquid nitrogen, and liquid hydrogen. Evaluation of results was based on linear elastic fracture mechanics concepts and was directed to improving existing methods of estimating minimum fracture strength and fatigue lives for pressurized structure in spacecraft and booster systems. Effects of specimen design in plane-strain fracture toughness testing were investigated. Four different specimen types were tested in room air, liquid nitrogen and liquid hydrogen environments using the aluminum and titanium alloys. Interferometry and holograph were used to measure crack-opening displacements in surface-flawed plexiglass test specimens. Comparisons were made between stress intensities calculated using displacement measurements, and approximate analytical solutions.

Hall, L. R.↗

Fundamental Mechanisms of Tensile Fracture in Aluminum Sheet Unidirectionally Reinforced with Boron Filament

Results are presented from an experimental research effort to gain a more complete understanding of the physics of tensile fracture in unidirectionally reinforced B-Al composite sheet. By varying the degree of filament degradation resulting from fabrication, composite specimens were produced which failed in tension by the cumulative mode, the noncumulative mode, or by any desired combination of the two modes. Radiographic and acoustic emission techniques were combined to identify and physically describe a previously unrecognized fundamental fracture mechanism which was responsible for the noncumulative mode. The tensile strength of the composite was found to be severely limited by the noncumulative mechanism which involved the initiation and sustenance of a chain reaction of filament fractures at a relatively low stress level followed by ductile fracture of the matrix. The minimum average filament stress required for initiation of the fracture mechanism was shown to be approximately 170 ksi, and appeared to be independent of filament diameter, number of filament layers, and the identity of the matrix alloy.

Herring, H. W.↗

The fracture energy and some mechanical properties of a polyurethane elastomer.

The energy required to form a unit of new surface in the fracture of a polyurethane elastomer is determined. The rate sensitivity of the material has been reduced by swelling it in toluene. This paper primarily describes the experimental work of measuring the lower limit of the fracture energy. With this value and the creep compliance as a basis, the rate dependence of fracture energy for the unswollen material has been determined. It is thus shown that the dependence of the fracture energy on the rate of crack propagation can be explained by energy dissipation around the tip of the crack. Good agreement between the theoretically and experimentally determined relationships for the rate-sensitive fracture energy is demonstrated.

Mueller, H. K.↗

Effect of specimen thickness of fatigue-crack-growth behavior and fracture toughness of 7075-T6 and 7178-T6 aluminum alloys

A study was made to determine the effects of specimen thickness on fatigue crack growth and fracture behavior of 7075-T6 and 7178-T6 aluminum alloy sheet and plate. Specimen thicknesses ranged from 5.1 to 12.7 mm (0.20 to 0.50 in.) for 7075-T6 and from 1.3 to 6.4 mm (0.05 to 0.25 in.) for 7178-T6. The stress ratios R used in the crack growth experiments were 0.02 and 0.50. For 7075-T6, specimen thickness had relatively little effect on fatigue-crack growth. However, the fracture toughness of the thickness of the thickest gage of 7075-T6 was about two-thirds of the fracture toughness of the thinner gages of 7075-T6. For 7178-T6, fatigue cracks generally grew somewhat faster in the thicker gages than in the thinnest gage. The fracture toughness of the thickest gage of 7178-T6 was about two-thirds of the fracture toughness of the thinner gages of 7178-T6. Stress intensity methods were used to analyze the experimental results. For a given thickness and value of R, the rate of fatigue crack growth was essentially a single-valued function of the stress intensity range for 7075-T6 and 7178-T6. An empirical equation developed by Forman, Kearney, and Engle fit the 7075-T6 and 7178-T6 crack growth data reasonably well.

Hudson, C. M.↗