Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture propagation”

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 163 records · Page 9

Finite domain solution of a hydraulic fracture in a permeable rock

In this work, we present a domain-based algorithm to simulate the propagation of a plane-strain hydraulic fracture in a zero-toughness permeable elastic medium. The algorithm utilizes a domain-based method to solve the elasticity equation and integrates a multi-scale tip asymptote, which is particular to hydraulic fractures, into this framework. This integration is key to accurately model the energy dissipation and the fluid leak-off in the fracture tip region. The algorithm combines a 2D finite volume method (FVM) for solving the elasticity equation with a 1D FVM for solving the nonlinear lubrication equation. Incorporating the far-field asymptotics and using a moving-mesh scheme reduces the computational burden while improving the accuracy of the scheme. The paper concludes with an analysis of the numerical results. This study demonstrates the potential of this domain-based approach for modeling hydraulic fractures in poroelastic media.

Domain-based method↗

Hydro-chemo-mechanical phase field formulation for corrosion induced cracking in reinforced concrete

Steel corrosion induced cracking in reinforced concrete structures is often caused by volumetric straining at the steel-concrete interface. This volumetric strain is mainly controlled by the rate of generation of corrosion product, which depends on the diffusion of Cl{sup −} and dissolved O{sub 2} as well as the chemical reaction at the steel-concrete interface. While many numerical modeling techniques have been proposed to tackle this complex phenomenon, most have been limited to using staggered solution schemes, where the diffusion, chemical reaction, mechanical strain, and crack fields are solved for separately. Unfortunately, such approaches introduce numerical errors that can lead to unrealistic predictions. Therefore, a coupled solution, is needed to remove the formulation induced errors and understand this complex phenomenon. To provide a unifying solution framework for corrosion induced cracking in reinforced concrete, a phase field formulation is presented and implemented with FEM algorithms. The governing equations are derived from the system's free energy and Fick's law of diffusion. The proposed formulation was first validated in 1D by comparing the modeling results with the semi-analytical solution approximated via Laurent series to avoid singularities in the complex domain. 2D plane strain solutions to the embedded rebar undergoing corrosion induced cracking with pre-cracks are also presented and compared with the approximated solutions from linear elastic fracture mechanics theory (LEFM). The modeling of crack propagation enabled by the proposed formulation is also presented along with parametric studies to reveal the roles of the fracture toughness and permeability of the concrete on the corrosion induced cracking. A case study was also conducted on H-pile steel with concrete jackets to demonstrate the feasibility of modeling corrosion induced fracture in reinforced concrete structures with complex steel-concrete interfaces.

36 MATERIALS SCIENCE↗

Investigating Grain Structure and Microcracking in SiCf-SiCm Composites Using 4D-STEM

Silicon carbide (SiC) fiber-reinforced SiC matrix composites (SiCf-SiCm) are promising materials for accident-tolerant fuel (ATF) cladding in light water reactors. This paper, using four-dimensional scanning transmission electron microscopy (4D-STEM), studied the microstructure of SiCf-SiCm at nanoscale and provided grain statistics at the CVI/fiber region. Here, the results reveal that CVI region mainly contain irregular sub-micron grains (115 nm to 1800 nm) with a preferred orientation, while the fiber region has nano equiaxed grains (21 nm to 78 nm) with random orientations. Centered dark field imaging over the fiber regions verified the grain size measurement by 4D-STEM. Through 4D-STEM virtual dark field imaging and orientation mapping, intragranular propagation was identified as the fracture mechanism for a microcrack observed within the CVI region. The amorphous pyrolytic carbon layer was also shown to effectively arrest the microcrack.

4D-STEM↗

Dynamic fracture processes in hydrogen embrittled iron

Mechanistic understanding of hydrogen embrittlement is critical to the transport and storage of hydrogen as a sustainable fuel. Experimental observation of dynamic fracture processes in hydrogen environments has been challenging. Here, in-situ transmission X-ray microscopy is used to image polycrystalline iron foils that are concurrently strained and electrochemically charged with hydrogen. Using this technique, voids and microcracks down to ∼30 nm in size are identified at or near the crack tip and tracked as strain is increased. We observe both hydrogen-induced transgranular failure and hydrogen-induced intergranular failure with distinct fracture processes. We find that transgranular fracture occurs through the formation of sharp branched structures at the main crack tip that are reminiscent of slip band formation. This is followed by void growth and coalescence, and interaction with secondary cracks that lead to propagation of the main crack. Intergranular fracture is found to occur through void coalescence along grain boundaries, along with significant secondary cracking near the crack tip. The crack growth rate is greater for void-mediated intergranular failure than for slip-mediated transgranular failure. Additionally, the removal of hydrogen results in recovered ductility for a crack initially propagated under hydrogen charging.

36 MATERIALS SCIENCE↗

Impacts of Mode Mixity on Controlled Spalling of (100)-Oriented Germanium

Controlled spalling is a technology to prepare single-crystal thin films of semiconductors by fracture with a subsurface crack propagating nearly parallel to the substrate surface. Practical applications require uniform thickness and a smooth surface across the whole film. Both wafer-scale and patterned-stressor-defined small-area spalling of germanium substrates are conducted experimentally and numerically. River line features are observed on spalled surfaces close to lateral edges of the spall, regardless of the spall direction and the size of the spalled area. Three-dimensional finite element method modeling shows the river lines are caused by mixed mode I?+?III loading near the lateral edges of spall and predicts a spall depth variation near the lateral edges of spall due to mixed mode I?+?II loading. The absolute range of river lines increases with lateral size of spall, while the relative range of river lines decreases, consistent with variations in mode mixity.

36 MATERIALS SCIENCE↗

Understanding geopolymer binder-aggregate interfacial characteristics at molecular level

The interfacial characteristics of geopolymer binder to aggregate composites are poorly understood, especially at molecular level. Herein, molecular models are developed to study, for the first time, the geopolymer-aggregate interface. Chemically, various forms of interfacial bonding are characterized, including Al-O-Si bonding through condensation reactions, NaO and H-bonding. An atomic-level interfacial transition zone (ITZ) is identified, attributed to the concentration of –OH groups. Increasing the Si/Al ratio of geopolymer is found to decrease the ITZ density, but have limited effect on the ITZ width. A heterogeneous diffusion characteristic occurs in geopolymer, due to the weak interfacial interaction. Mechanically, lowering the Si/Al ratio promotes the interfacial strength due to the stronger interfacial interaction and higher cross-linking degree in geopolymer. Under loading the interfacial fracture undergoes three stages: crack propagation, chain bridging (including aluminosilicate and ionic bridging) and breakage. The above atomic-level findings may facilitate a better design of geopolymer concrete in engineering.

36 MATERIALS SCIENCE↗

Residual fracture energy of natural and recycled aggregate concrete after exposure to high temperatures

Abstract Fracture energy ( G F ) represents the average total energy consumed in the whole process of crack propagation. Using one parameter of G F , it is not sufficient to describe the fracture behavior during crack stable and unstable propagation periods. In this study, the G F , stable fracture energy ( G FS ) and unstable fracture energy ( G FU ) of natural and recycled aggregate concrete (RAC) were investigated by using three‐point bending tests. Totally, 128 cubic specimens (100 mm × 100 mm × 100 mm) and 64 notched beams with size of 510 mm × 100 mm × 100 mm were manufactured to test the compressive strength ( f cu ), splitting tensile strength ( f t ), load–deflection (P‐δ), and load‐cracking mouth opening displacement (P‐CMOD) curves, considering different replacement levels (0%, 30%, 70%, and 100%) of recycled coarse aggregate (RA) and high temperatures exposure (20, 100, 200, and 300°C). The effects of RA substitutions and temperatures on basic mechanical properties, G F , G FS , and G FU were analyzed. The results demonstrated that the f cu , f t , and elastic modulus ( E ) decreased with the incremental temperature. The G FS of concrete with different RA substitutions first decreased and then increased, and the minimum value appeared at 200°C. The G F and G FU of concrete with different RA substitutions increased with the elevated temperature, and the increase rate decreased with RA substitutions. When the temperature is less than 100°C, the G F and G FU of concrete remained mostly unchanged with the increase of RA substitutions, while the temperature reached 200 and 300°C, the G F and G FU of concrete increased and decreased with the increase of RA substitutions, respectively.

Chen, Wei↗

Mechanism of Microseismic Generation During Hydraulic Fracturing - With Evidence From HFTS 2 Observations

The objective of this study is to understand how microseismic events are generated during hydraulic fracturing, as well as the role of geomechanical conditions (i.e., stress and mechanical stratigraphy) in this process. In the industry, microseismic event clouds have been generally used as an “outer-boundary” of the “stimulated reservoir volume” (SRV). However, by comparing with other surveillance data (Low frequency Distributed Acoustic Sensing, or LF-DAS strain) in the Hydraulic Fracturing Test Site (HFTS) 2 experiment, we show that this assumption is fundamentally flawed. The HFTS 2 data has three unique observations that have not been commonly observed in other datasets: 1. Due to influence of offset pad depletion, microseismic data shows that hydraulic fractures from the child well can propagate over 3000 feet into the depleted low stress zone. 2. By comparing microseismic and horizontal fiber LF-DAS strain data, we observe that microseismic event cloud does not necessarily reflect the created hydraulic fracture volume. Particularly, the extent of microseismic event clouds near heel stages are much shorter than what is shown with LF-DAS strain data. 3. Microseismic event magnitudes are larger in the depleted regions. Through geomechanical analysis, we demonstrate that the “bedding-plane-slip” model is likely the mechanism for microseismic generation during hydraulic fracturing. Furthermore, this model successfully explains the above field observations from HFTS 2 experiment. We also provide a quantitative relationship connecting the microseismic event magnitude with fracture width increment and layer mechanical property contrast.

58 GEOSCIENCES↗

Interpretation of HFTS 2 Microseismic Data Using Bedding-Plane-Slip Mechanism

Summary The objective of this study is to understand how microseismic events are generated during hydraulic fracturing, as well as the role of geomechanical conditions (i.e., stress and mechanical stratigraphy) in this process. In the industry, microseismic event clouds have been generally used as an “outer boundary” of the “stimulated reservoir volume” (SRV). However, by comparing with other surveillance data (low-frequency distributed acoustic sensing, or LF-DAS, strain) in the Hydraulic Fracturing Test Site (HFTS) 2 experiment, we show that this assumption is fundamentally flawed. The HFTS 2 data have three unique observations that have not been commonly observed in other datasets: (1) Due to the influence of offset pad depletion, microseismic data show that hydraulic fractures from the child well can propagate over 3,000 ft into the depleted low-stress zone. (2) By comparing microseismic and horizontal fiber LF-DAS strain data, we observe that the microseismic event cloud does not necessarily reflect the created hydraulic fracture volume. Particularly, the extent of microseismic event clouds near heel stages is much shorter than what is shown with LF-DAS strain data. (3) Microseismic event magnitudes are larger in the depleted regions. Through geomechanical analysis, we demonstrate that the “bedding-plane-slip” model is likely the mechanism for microseismic generation during hydraulic fracturing. This model successfully explains the above field observations from the HFTS 2 experiment. We also provide a quantitative relationship connecting the microseismic event magnitude with fracture width increment and layer mechanical property contrast.

Engineering↗

Anisotropic thermal-conductivity degradation in the phase-field method accounting for crack directionality

Dynamic loads applied to metals may lead to brittle or ductile fracture depending on the imposed strain rates, material properties, and specimen geometry. With an increase in the applied velocity, brittle to ductile failure mode transition may also be observed. Furthermore, at high strain rates, ductile fracture may be preceded by shear bands, which are narrow bands of intense plastic deformation, typically accompanied by a significant rise in temperature. Reliable models are needed to predict the response of metals subject to dynamic loads. In particular, capturing the interplay between heat conduction and crack propagation using the phase-field fracture method is still an open research field. To accurately capture the heat transfer physics across crack surfaces, damage models degrading thermal-conductivity are necessary. While isotropic thermal-conductivity degradation models were proposed, they may lead to errors as they are indifferent to crack directionality. In our current work, a new anisotropic approach is proposed in which thermal-conductivity, which depends on the phase-field gradient, is degraded solely across the crack. Here, it is shown that this approach improves the near-field approximation of temperature and heat flux compared with isotropic degradation, when taking the discontinuous crack solutions as reference. Additionally, the anisotropic degradation approach is implemented in a unified dynamic fracture model and a couple of examples demonstrate the viability of this approach.

42 ENGINEERING↗

Impact of position and density of nanoscale voids on fracture initiation in iron from phase field fracture simulation

Understanding the impact of these bubbles on crack propagation, like that of helium bubble-induced cracking in irradiated materials is incredibly complex. A useful first study towards understanding bubble effects on fracture is to examine how voids impact fracture first. In this work, we used phase-field fracture simulations to examine the influence of voids and their distribution on Mode I fracture in Fe. Assuming brittle fracture, two simulation configurations were considered: (1) nanoscale systems with one or two voids, and (2) nanoscale systems with an experimentally relevant distribution of voids, with up to 20% void area. Results from simulations with one and two voids showed that voids within 10 nm of a crack tip reduce the stress required for crack growth, with the magnitude of reduction depending on void-to-crack orientation. Comparisons with linear elastic fracture mechanics and evaluation of one versus two void systems revealed deviations from linear superposition, implying complex interactions between void and crack tip stress fields. In multi-void simulations, as void sizes increase, the nearest void to the crack tip exerts a greater influence on fracture stress than the overall porosity. Furthermore this study provides valuable insights into the relationship between void size and concentration, and the stress necessary for crack growth, marking a step forward towards understanding He bubble-induced fracture in ferrous materials.

36 MATERIALS SCIENCE↗

Utah FORGE: Direct Shear Test Data for Investigating Seismic Precursors to Shear Failure of Fractures

This dataset includes results of direct shear tests to investigate the mechanical and geophysical response of dry and saturated fractures in Indiana limestone and Sierra White granite. Direct shear tests were performed on tensile-induced fractures in Indiana limestone and Sierra White granite in a custom water-pressurized chamber. The provided Excel files include the representative seismic wave signals and the normalized wave amplitudes of ultrasonic wave transducers. A link to the published journal article presenting the data and describing the experiment in detail is provided as well.

15 GEOTHERMAL ENERGY↗

Large scale polymer toughening of two-dimensional materials revealed by in situ TEM fracture tests and multiscale simulations

Two-dimensional (2D) materials offer significant potential for applications in energy-harvesting devices, batteries, sensors, and transistors. However, their intrinsic brittleness makes them prone to mechanical failure, limiting their practical use. In this work, we perform in situ transmission electron microscopy (TEM) fracture tests on monolayer MoSe2 and uncover an extrinsic toughening effect induced by an ultrathin adsorbed polystyrene adlayer. This adlayer substantially enhances the fracture resistance of the 2D flakes. Through a combination of molecular dynamics simulations and finite element analysis, we elucidate the molecular mechanism behind this toughening effect. Further, it arises from the active crack-bridging behavior of entangled polymer chains and the formation of a fracture process zone that stabilizes crack propagation and increases the energy required for crack extension. The proposed toughening mechanism offers a pathway to improving the mechanical reliability of 2D material-based devices by mitigating the risk of sudden failure.

2D materials↗

Phase-field modeling of crack growth under coupled creep-fatigue

Crack growth under coupled creep-fatigue strongly influences service life of metallic components operating at high temperatures. However, there lack numerical models that enable direct simulations of crack growth and plasticity development under coupled creep-fatigue. Here, in this study, a unified phase-field model is developed for simulating crack growth under fatigue, creep, and coupled creep-fatigue. This model is able to reproduce the Paris laws for both brittle and ductile materials in pure fatigue mode, which shows that the degradation rate due to cyclic stress in ductile materials is generally larger than that in brittle materials. Fatigue simulations with different yield strengths demonstrate that plasticity facilitates crack growth in typical ductile materials. High-throughput simulations are performed under coupled creep-fatigue conditions, with varying hold-stress levels and hold-time. By analyzing crack growth rate under different hold-time, a creep-fatigue interaction term is obtained. Simulation results demonstrate that the creep-fatigue interaction is caused by smaller stress gradient and enhanced degradation rate ahead of the crack tips under coupled creep-fatigue, compared to the situation of pure fatigue. This work reveals the origin of different Paris-law exponents in brittle and ductile materials and demonstrates how the interplay between creep and fatigue affects crack growth.

36 MATERIALS SCIENCE↗

Radiation damage effects in beryllium for next generation neutrino beam targetry (Final Technical Report)

Current and future high-power accelerators put severe requirements on materials used for target and beam windows and target facilities have been recognized as a critical challenge in development of future particle accelerators. In accelerators, window and target materials are exposed to extreme conditions, which include bombardment with very high energy protons (1- 100 GeV) and thermomechanical shock waves. Radiation can cause direct damage in the material, and it leads to production of transmutation products (especially helium), both phenomena having a potential adverse effect on the stability and durability of the target/window material. At high enough temperatures, He can aggregate to form gas bubbles, which in turn cause significant dimensional changes (swelling), enable easy crack propagation, and eventually cause failure by fracture. On the other hand, if the temperature is too low, radiation damage accumulates in the form of internal defects (e.g., dislocations), leading to hardening and a decreased ductility of the material. In this project, we will focus on beryllium since it is considered to be one of the candidate materials for beam windows and targets in the next-generation proton accelerators, e.g., the Long Baseline Neutrino Facility (LBNF). Radiation effects in Be have been studied in the context of nuclear fusion reactor applications. However, key differences exist between reactor and accelerator conditions, including neutron vs. proton irradiation, continuous vs. pulsed beam flux, much higher energies of bombarding particles in accelerators, and higher operating temperatures for typical reactors. For example, the impact of beam pulsing on the radiation damage and the He bubble kinetics is largely unknown. While results obtained on Be from fusion research might not be directly transferrable to understanding target materials, there is an opportunity to bring state-of-the-art tools from materials research in nuclear reactors to aid design of target and beam window materials in high-power accelerators. To this end, the overarching goal of this project are to develop an experimentally-validated computational framework capable of predicting radiation damage evolution in beryllium relevant to beam window and target conditions, focusing on He bubble formation and growth as a function of irradiation temperature. Our model will be based on the cluster dynamics formalism, where size distribution of defects and He bubbles is simulated as a function of time, temperature, and radiation dose. Parameters for the model will be taken from published experiments and from high-fidelity atomistic simulations proposed in this project. In addition, we will carry out a series of targeted ex-situ and in-situ dual-beam experiments using low-energy protons to provide critical data for validation of the model on the effects of radiation on He clustering, He bubble distribution, and dislocation loop density/size in proton irradiated Be.

36 MATERIALS SCIENCE↗

On the Mechanical Behavior of LP-DED C103 Thin-Wall Structures

Laser Powder Directed Energy Deposition (LP-DED) can produce thin-wall features on the order of 1 mm. These features are essential for large structures operating in extreme environments such as regeneratively cooled nozzles and heat exchangers, which often make use of refractory metals. In this work, the mechanical behavior of LP-DED C103 was investigated via quasi-static tensile testing and low cycle fatigue (LCF) testing. The effects of vacuum stress relief (SR) and hot isostatic pressing (HIP) heat treatments were investigated for specimens in the vertical and horizontal build orientations during tensile testing. The AB and SR properties were lower than literature values for wrought and laser powder bed fusion (L-PBF) bulk components but higher than electron beam powder bed fusion (EB-PBF). The application of a HIP cycle improved strength by 7% and ductility by 27% past the initial as-built condition. Fracture images reveal that interlayer stress concentration sites are responsible for fracture in specimens in the vertical orientation. Meanwhile, fracture in the horizontal specimens mainly propagates at a slanted angle typical of plane stress conditions. The LCF results show cycles to failure ranging from 100 cycles to 8000 cycles for max strain levels of 2% and 0.5%, respectively. Fractography on the fatigue specimens reveals an increasing propagation zone as max strain levels are increased. The impact of these findings and future work are discussed in detail.

36 MATERIALS SCIENCE↗

Investigation of acoustic waves under subsurface conditions to improve the predictions of rock mechanical properties and natural fracture characteristics

Mechanical properties and natural fracture characteristics are critical to investigate for subsurface engineering applications, including carbon storage, well drilling, and stimulation, as they govern rock stability, fluid flow, and mechanical behavior under stress. This dissertation integrates experimental and machine learning approaches to enhance the prediction and understanding of these properties by analyzing acoustic wave behavior under varied subsurface conditions. First, the influence of temperature, pore pressure, and supercritical CO2 (scCO2) saturation on poroelastic properties is examined using Gray Berea sandstone samples. The results show that temperature and pore pressure significantly affect the bulk modulus and Biot’s coefficient, while scCO2 saturation impacts rock compressibility, informing strategies for effective geological carbon storage. The study extends this understanding by experimentally evaluating the impact of reservoir depletion on the dynamic mechanical properties of the emerging Caney shale in South Oklahoma with the employment of unsupervised machine learning to predict static mechanical properties across the Caney shale. Integrating petrophysical data and chemostratigraphy, the workflow—featuring K-means clustering, principal component analysis (PCA), and inverse distance weighting (IDW)—improves stratigraphic characterization and the estimation of static-to-dynamic modulus ratios, which is vital for optimizing drilling and stimulation strategies. Finally, the work explores how natural fracture characteristics in shale influence acoustic waveforms and shear wave splitting (SWS) analysis. Experimental data on fractured samples under different stress and temperature conditions, combined with machine learning models such as K-nearest neighbors (KNN) and extreme gradient boosting (XGBoost), reveal key fracture properties impacting SWS and wave propagation. Together, these studies provide a comprehensive framework for linking acoustic wave behavior with rock properties, advancing the methods for monitoring and predicting geomechanical changes. The insights offered valuable implications for safer, more efficient CO2 injection, hydrocarbon extraction, and subsurface management.

Elkholy, Sherif↗

Finite deformation implementation of a mixed-mode single-integral type cohesive zone with reorienting surfaces of separation

To model material ductile failure and crack propagation, cohesive zone elements can be embedded along potential fracture paths in a finite element simulation. When damage criteria are met, elements in the mesh decohere, simulating the formation and propagation of a crack. In this paper, we present a novel computational algorithm based on finite deformation theory, essential to modeling crack initiation and growth in solids undergoing large deformations. This new algorithm was formulated within a Lagrangian frame of reference to extend previous cohesive zone algorithms to include modeling crack growth in finite deformation contexts. The local coordinate system, necessary for defining an embedded cohesive zone, is constructed based upon the current configuration and is updated within the nonlinear iteration process, thereby resulting in the convergence of the solution for a growing crack in a large deformation quasi-static setting. The model’s accuracy was demonstrated by comparing finite element model simulation results with the analytic case of a constant surface separation, as shown in the verification examples. The power and efficacy of the algorithm to capture large deformations during crack growth were then demonstrated with a double cantilever beam example case. It indicates that the model can be applied to a variety of physical circumstances for predicting crack initiation and growth with delamination and fracture.

42 ENGINEERING↗