Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracturing materials”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Artificial neural network potentials for mechanics and fracture dynamics of two-dimensional crystals **

Understanding the mechanics and failure of materials at the nanoscale is critical for their engineering and applications. The accurate atomistic modeling of brittle failure with crack propagation in covalent crystals requires a quantum mechanics-based description of individual bond-breaking events. Artificial neural network potentials (NNPs) have emerged to overcome the traditional, physics-based modeling tradeoff between accuracy and accessible time and length scales. Previous studies have shown successful applications of NNPs for describing the structure and dynamics of molecular systems and amorphous or liquid phases of materials. However, their application to deformation and failure processes in materials is still uncommon. In this study, we discuss the apparent limitations of NNPs for the description of deformation and fracture under loadings and propose a way to generate and select training data for their employment in simulations of deformation and fracture simulations of crystals. We applied the proposed approach to 2D crystalline graphene, utilizing the density-functional tight-binding method for more efficient and extensive data generation in place of density functional theory. Then, we explored how the data selection affects the accuracy of the developed artificial NNPs. It revealed that NNP’s reliability should not only be measured based on the total energy and atomic force comparisons for reference structures but also utilize comparisons for physical properties, e.g. stress–strain curves and geometric deformation. In sharp contrast to popular reactive bond order potentials, our optimized NNP predicts straight crack propagation in graphene along both armchair and zigzag (ZZ) lattice directions, as well as higher fracture toughness of ZZ edge direction. Our study provides significant insight into crack propagation mechanisms on atomic scales and highlights strategies for NNP developments of broader materials.

2D materials↗
DOE OSTI2021

Connectivity-informed drainage network generation using deep convolution generative adversarial networks

Abstract Stochastic network modeling is often limited by high computational costs to generate a large number of networks enough for meaningful statistical evaluation. In this study, Deep Convolutional Generative Adversarial Networks (DCGANs) were applied to quickly reproduce drainage networks from the already generated network samples without repetitive long modeling of the stochastic network model, Gibb’s model. In particular, we developed a novel connectivity-informed method that converts the drainage network images to the directional information of flow on each node of the drainage network, and then transforms it into multiple binary layers where the connectivity constraints between nodes in the drainage network are stored. DCGANs trained with three different types of training samples were compared; (1) original drainage network images, (2) their corresponding directional information only, and (3) the connectivity-informed directional information. A comparison of generated images demonstrated that the novel connectivity-informed method outperformed the other two methods by training DCGANs more effectively and better reproducing accurate drainage networks due to its compact representation of the network complexity and connectivity. This work highlights that DCGANs can be applicable for high contrast images common in earth and material sciences where the network, fractures, and other high contrast features are important.

42 ENGINEERING↗
DOE OSTI2025

FY25 Status Update on A709 Code Case Testing at ANL

This report provides an update on the status of the creep rupture testing on the precipitation-treated (PT) Alloy 709 samples fabricated from the first, the second and the third commercial heats, in support of the first, 100,000-hour Alloy 709 Code Case development under the American Society of Mechanical Engineers (ASME) Section III, Division 5 for high-temperature reactor construction. In Fiscal Year (FY) 25, no new tests were initiated, while 5 previously loaded tests were completed. This report presents an overview of the creep-rupture data gathered at Argonne National Laboratory (ANL) since FY21 that totaled 52 completed tests. Metallographic images showing crack morphologies in some representative samples ruptured in FY25 are also included. Key conclusions include: 1) the creep behavior of the three heats of Alloy 709 can be well described by the Larson-Miller relationship; 2) higher temperatures and larger stresses lead to shorter creep life and higher minimum creep rate; 3) Alloy 709 is overall very ductile under creep; 4) the crack morphologies near fracture surfaces are dependent on material’s grain size and test conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Fatigue and Fracture Behavior of Additively Manufactured Austenitic Stainless Steel

Additive manufacturing (AM) includes a diverse suite of innovative manufacturing processes for producing near-net shape components, typically from powder or wire feedstock. Reported mechanical properties of AM materials vary significantly depending on the details of the manufacturing process and the characteristics of the processing defects (namely, lack of fusion defects). However, an excellent combination of strength, ductility and fracture resistance can be achieved in AM type 304L and 316L austenitic stainless steels by minimizing processing defects. It is also important to recognize that localized solidification processing during AIVI produces microstructures more analogous to weld microstructures than wrought microstructures. Consequently, the mechanical behavior of AM austenitic stainless steels in harsh environments can diverge from the performance of wrought materials. This report gives an overview of the fracture and fatigue response of type 304L materials from both directed energy deposition (DED) and powder bed fusion (PBF) techniques. In particular, the mechanical performance of these materials is considered for high-pressure hydrogen applications by evaluating fatigue and fracture resistance after thermally precharging of test specimens in high-pressure gaseous hydrogen. The mechanical behaviors are considered with respect to previous reports on hydrogen-assisted fracture of austenitic stainless steel welds and the unique characteristics of the AM microstructures. Fatigue crack growth can be relatively insensitive to processing defects, displaying similar behavior as wrought materials. Fracture resistance of dense AM austenitic stainless steel, on the other hand, is more consistent with weld metal than with compositionally-similar wrought materials. Hydrogen effects in the AM materials are generally more severe than in wrought materials, but comparable to measurements on welded austenitic stainless steels in hydrogen environments. While hydrogenassisted fracture manifests differently in welded and AM austenitic stainless steel, the fracture process appears to have a common origin in the compositional microsegregation intrinsic to solidification processes.

additive manufacturing, austenitic stainless steel↗
DOE OSTI2024

Full-field quantitative visualization of shock-driven pore collapse and failure modes in PMMA

The dynamic collapse of pores under shock loading is thought to be directly related to hot spot generation and material failure, which is critical to the performance of porous energetic and structural materials. However, the shock compression response of porous materials at the local, individual pore scale is not well understood. This study examines, quantitatively, the collapse phenomenon of a single spherical void in PMMA at shock stresses ranging from 0.4 to 1.0 GPa. Using a newly developed internal digital image correlation technique in conjunction with plate impact experiments, full-field quantitative deformation measurements are conducted in the material surrounding the collapsing pore for the first time. The experimental results reveal two failure mode transitions as shock stress is increased: (i) the first in situ evidence of shear localization via adiabatic shear banding and (ii) dynamic fracture initiation at the pore surface. Numerical simulations using thermo-viscoplastic dynamic finite element analysis provide insights into the formation of adiabatic shear bands (ASBs) and stresses at which failure mode transitions occur. Further numerical and theoretical modeling indicates the dynamic fracture to occur along the weakened material inside an adiabatic shear band. Finally, analysis of the evolution of pore asymmetry and models for ASB spacing elucidate the mechanisms for the shear band initiation sites, and elastostatic theory explains the experimentally observed ASB and fracture paths based on the directions of maximum shear.

42 ENGINEERING↗
DOE OSTI2025

Computer Vision Pipeline for Image Analysis for Freeze‐Fracture Electron Microscopy: Rosette Cellulose Synthase Complexes Case

In materials science, plant biology, agriculture, and environmental research, the automated analysis of high-magnification, complex microscopy images, such as those generated by freeze-fracture electron microscopy (FF-TEM), remains a critical challenge that limits the scalability of data interpretation. We present a deep learning computer vision pipeline for high-throughput detection and morphological characterization analysis of cellulose synthase complexes (CSCs, or rosettes) in FF-TEM images. The pipeline integrates preprocessing, detection, human-in-the-loop verification, and semantic segmentation to quantify features such as rosette diameter and inter-lobe spacing. The approach was trained and tested on a curated dataset of high-resolution FF-TEM micrographs of Physcomitrium patens, expanded via strategic tiling and augmentation to over 650 images. We compare YOLOv8 and YOLOv9 architectures and demonstrate that YOLOv9 achieves superior performance in both localization accuracy (mAP50-95 = 0.854) and inference speed. The resulting distributions revealed biological variability consistent with prior manual studies, validating the approach for high-throughput applications. Our results show that the pipeline achieves human-expert level accuracy while dramatically reducing analysis time, enabling scalable, reproducible structural characterization of intramembrane protein complexes. The pipeline is broadly applicable to other domains requiring precise interpretation of complex microscopy data and establishes a foundation for future artificial intelligence (AI)-assisted workflows in biological imaging.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Meshfree Methods

Meshfree methods have undergone substantial development and have received much attention in the last two decades. This new family of numerical methods is designed to inherit the main advantages of the finite element method such as compact supports of shape functions and good approximation properties while, at the same time, overcome the main disadvantages of the finite element method caused by the mesh dependence. The meshfree methods share a common feature that no mesh is needed and shape functions are constructed from sets of points, thus eliminating the need for time consuming mesh generation. The most significant advantage of meshfree methods is the flexibility in customizing approximation functions for desired regularity and for capturing essential physics and features of the particular problems of interest. Adaptivity formulation and multiple-scale solution strategies also can be implemented with relative ease. It has become clear that the meshfree methods provide considerable advantages over the conventional finite element methods in solving problems involving moving discontinuities, evolving material interfaces, multiple-scale phenomena, large material distortion and structural deformation, and fracture and damage processes. This Chapter gives an overview of many classes of meshfree methods, with more detailed discussions on Smoothed Particle Hydrodynamics (SPH), the Reproducing Kernel Particle Method (RKPM), Peridynamics (PD), the Material Point Method (MPM), as well as their applications in various challenging engineering problems.2

Chen, Jiun-Shyan↗
DOE OSTI2021

In Situ TEM Tensile Testing of Bicrystals with Tailored Misorientation Angles

Grain boundaries have complex structural features that influence strength, ductility and fracture in metals and alloys. Grain boundary misorientation angle has been identified as a key parameter that controls their mechanical behavior, but the effect of misorientation angle has been challenging to isolate in polycrystalline materials. Here, we describe the use of bicrystal Au thin films made using a rapid melt growth process to study deformation at a single grain boundary. Tensile testing is performed on bicrystals with different misorientation angles using in situ TEM, as well as on a single crystalline sample. Plastic deformation is initiated through dislocation nucleation from free surfaces. Grain boundary sliding is not observed, and failure occurs away from the grain boundary through plastic collapse in all cases. The failure behavior in these nanoscale bicrystals does not appear to depend on the misorientation angle or grain boundary energy but instead has a more complex dependence on sample surface structure and dislocation activity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Composition-dependent fracture energy in metallic glasses

The interplay between metallic glasses (MGs) mechanical properties, fracture energy (G), and glass-forming ability (GFA) and their dependence on alloy composition remain poorly understood. Here, we perform molecular dynamics simulations to investigate the intrinsic composition-dependence of G in Cu x Zr 100-x MGs (x = 20, 30, 40, 50, 64). The results indicate that the value of G increases with Cu content. In addition, it is revealed that MGs with higher G values display higher Poisson’s ratio (ν) and GFA, suggesting a close correlation between fracture toughness, mechanical properties, and GFA. This correlation between G, ν and GFA can be understood based on the fragility (m) of supercooled liquids, which is directly related to the structural heterogeneity in MGs. Larger m values are related to dynamic slowdown and supercooled liquid stabilization, which enhance GFA and the formation of pronounced structural heterogeneity, comprised of loosely packed regions that favor β-relaxation and the activation of shear transformation zones. Those concurrently promote the expansion of the plastic zone at the crack tip, enhancing the observed value of G. Furthermore, these simulation results shed light on the intrinsic relationship between fracture toughness, mechanical properties, and alloy composition in MGs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

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↗
DOE OSTI2020

Scale bridging damage model for quasi-brittle metals informed with crack evolution statistics

Computationally efficient methods for bridging length scales, from highly resolved micro/meso-scale models that can explicitly model crack growth, to macro-scale continuum models that are more suitable for modeling large parts, have been of interest to researchers for decades. In this work, an improved brittle damage model is presented for the simulation of dynamic fracture in continuum scale quasi-brittle metal components. Crack evolution statistics, including the number, length, and orientation of individual cracks, are extracted from high-fidelity, finite discrete element method (FDEM) simulations and used to generate effective material moduli that reflect the material’s damaged state over time. This strategy allows for the retention of small-scale physical behaviors such as crack growth and coalescence in continuum scale hydrodynamic simulations. However, the high-fidelity simulations required to generate the crack statistics are computationally expensive. Thus, steps were taken to produce a flexible constitutive model to reduce the number of costly high-fidelity simulations needed to produce accurate results. A new stress based degradation criterion is introduced for the degradation of individual material zones. This allows for the development of a heterogeneous damage distribution within the bulk material. Then a flow stress model is added to the hydrodynamic simulation to account for plasticity in quasi-brittle materials. As a result, the effective moduli model can be applied to a larger range of materials. The effective moduli constitutive model is used to simulate beryllium flyer plate experiments. The results from the continuum scale simulations using statistics from a single high-fidelity simulation are found to be in excellent agreement with numerical and experimental velocity interferometer data. The same set of crack statistics are used to extrapolate the results of a higher rate flyer plate case using the effective moduli model. In conclusion, the extension of this model to higher rate cases shows promise for further reducing the number of costly high-fidelity simulations needed to generate crack statistics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Constitutive modeling of δ-phase zircaloy hydride based on strain rate dependent nanoindentation and nano-scale impact dataset

In this work, a constitutive model for zircaloy-4 embedded with δ-phase hydride rim/blister structures (δ-phase zircaloy hydride) is presented. Nanoindentation and nano-scale impact experiments were performed on the δ-hydride blisters as a function of temperature in order to obtain strain rate dependent location specific constitutive response. A power law viscoplastic constitutive model was chosen to represent the strain rate dependent behavior. The strength coefficient and the strain hardening exponent in the power law relation were found to decrease with increase in temperature indicating softening behavior at higher temperatures. The strain rate exponent remained constant with increase in temperature indicating that the rate dependent ductility does not change significantly at higher temperatures. A finite element method based model was created to analyze the failure behavior and the ductile to brittle transition in zircaloy hydride. The model predicted that fracture strength of δ-phase zircaloy hydride decreases linearly with increase in temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Grain-boundary fracture mechanisms in Li 7 La 3 Zr 2 O 12 (LLZO) solid electrolytes: When phase transformation acts as a temperature-dependent toughening mechanism

Garnet-type, solid electrolytes, such as Li 7 La 3 Zr 2 O 12 (LLZO), are a promising alternative to liquid electrolytes for lithium-metal batteries. However, such solid-electrolyte materials frequently exhibit undesirable lithium (Li) metal plating and fracture along grain boundaries. In this work, we employ atomistic simulations to investigate the mechanisms and key fracture properties associated with intergranular fracture along one such boundary. Our results show that, in the case of a Σ5 (310) grain boundary, this boundary exhibits brittle fracture behavior, i.e. the absence of dislocation activity ahead of the propagating crack tip, accompanied with a decrease in work of separation, peak stress, and maximum stress intensity factor as the temperature increases from 300 K to 1500 K. As the crack propagates, we predict two temperature-dependent Li clustering regimes. For temperatures at or below 900 K, Li tends to cluster in the bulk region away from the crack plane driven by a void-coalescence mechanism concomitant a simultaneous cubic-to-tetragonal phase transition. The tetragonalization of LLZO in this temperature regime acts as an emerging toughening mechanism. At higher temperatures, this phase transition mechanism is suppressed leading to a more uniform distribution of Li throughout the grain-boundary system and lower fracture properties as compared to lower temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental and theoretical study of the restrained shrinkage cracking of early age well cement

The shrinkage cracking occurring in oil & gas well cement material is a major factor that reduces the cement sheath integrity and impairs the durability of an underground well. Shrinkage mainly takes place in the early hydration period of the cement, and the resulting restrained shrinkage cracking can form a leaking network for oil & gas or carbon dioxide stored in an underground reservoir. It is thus important to understand the early age behavior of well cement and the leaking pathway formation in the well. This study describes the application of the digital image correlation (DIC) method as a tool to measure the strain development and crack distributions of well cement in the lab during early phase changes. A theoretical model is presented to calculate the stress distribution built up in the cement sheath during early ages of well cement hydration, which helps to interpret the experimental observations and determine the major factors influencing the cracking sequences. This study presents the whole hydration – drying shrinkage – and cracking mechanism of early age well cement under restrained ring conditions experimentally and theoretically. Finally, the results offer useful insights on the key material and geometric parameters which control the fracture of a well cementing structure.

36 MATERIALS SCIENCE↗
DOE OSTI2026

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↗
DOE OSTI2023

A phase-field model for hydraulic fracture nucleation and propagation in porous media

Many geo-engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase-field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase-field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy-based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase-field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase-field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. In conclusion, the comparison of numerical results for two-dimensional test cases with existing analytical solutions demonstrates that the proposed phase-field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three-dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Use of Mini-CT Specimens for Fracture Toughness Characterization of Irradiated Highly Embrittled Weld

In the 1990’s, the Heavy Section Steel Irradiation Program at the ORNL performed investigation of the shape of the fracture toughness master curve for reactor pressure vessel steel highly embrittled because of irradiation exposure. A radiation-sensitive reactor pressure vessel (RPV) weld with intentionally enhanced copper content, designated KS-01, has been characterized in terms of static initiation (KJc) and Charpy impact toughness in the unirradiated and irradiated conditions. The objective of this original project was to investigate the ability of highly embrittled material to maintain the shape of the unirradiated transition fracture toughness curve, as well as to examine the ability of the Charpy 41-J shift to predict the fracture toughness shift at such high level of embrittlement. Irradiation of this weld was performed at the University of Michigan Ford Reactor in the new HSSI irradiation-anneal-reirradiation (IAR) facility. Broken specimens from that project have been saved in ORNL storage. To verify applicability of Mini-CT specimens for fracture toughness characterization of RPV materials as part of the DOE Light Water Reactor Sustainability Program, Mini-CT specimens were machined from broken Charpy specimens and tested according to ASTM E1921 standard. As result of this study, the fracture toughness of this weld derived by testing Mini-CT specimens in the unirradiated and irradiated conditions is compared to previously reported fracture toughness derived by large number of conventional specimens.

Sokolov, Mikhail A.↗
DOE OSTI2021

Crack tip microplasticity mediated by microstructure gradients

Abstract Traditional fracture theories infer damage at cracks (local field) by surveying loading conditions away from cracks (far field). This approach has been successful in predicting ductile fracture, but it normally assumes isotropic and homogeneous materials. However, myriads of manufacturing procedures induce heterogeneous microstructural gradients that can affect the accuracy of traditional fracture models. This work presents a microstructure‐sensitive finite element approach to explore the shielding effects of grain size and crystallographic orientation gradients on crack tip microplasticity and blunting. A dislocation density‐based crystal plasticity model conveys texture evolution, grain size effects, and directional hardening by computing the constraint from dislocation structures. The results demonstrate that the microstructure can act as a buffer between the local and far fields that affects the crack tip microplasticity variability. For nominal opening loading, grain size and texture affect the local ductility and induce a non‐negligible multiaxial plastic deformation. Furthermore, driving forces based on measuring displacements away from the crack tip are less affected by the microstructure, which suggests that traditional experimental methods smear out important crack tip variability.

Castelluccio, Gustavo M.↗