Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “crack nucleation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

An Experimental Study of Incremental Surface Loading of an Elastic Plate: Application to Volcano Tectonics

Models of surface fractures due to volcanic loading an elastic plate are commonly used to constrain thickness of planetary lithospheres, but discrepancies exist in predictions of the style of initial failure and in the nature of subsequent fracture evolution. In this study, we perform an experiment to determine the mode of initial failure due to the incremental addition of a conical load to the surface of an elastic plate and compare the location of initial failure with that predicted by elastic theory. In all experiments, the mode of initial failure was tension cracking at the surface of the plate, with cracks oriented circumferential to the load. The cracks nucleated at a distance from load center that corresponds the maximum radial stress predicted by analytical solutions, so a tensile failure criterion is appropriate for predictions of initial failure. With continued loading of the plate, migration of tensional cracks was observed. In the same azimuthal direction as the initial crack, subsequent cracks formed at a smaller radial distance than the initial crack. When forming in a different azimuthal direction, the subsequent cracks formed at a distance greater than the radial distance of the initial crack. The observed fracture pattern may explain the distribution of extensional structures in annular bands around many large scale, circular volcanic features.

Williams, K. K.↗
DOE OSTI2023

On the convexity of phase-field fracture formulations: Analytical study and comparison of various degradation functions

Efficient and accurate fracture modeling is of great importance in applications where catastrophic outcomes under extreme scenarios are possible. The phase-field (PF) approach to fracture received significant attention over the past decade, due to its capability to capture complicated fracture patterns (e.g., crack merging and branching). Specifically, crack initiation and propagation are modeled via minimization of the total energy functional, which is regularized with the aid of a phase field. Despite the promising results and modeling capabilities of the PF method in many applications, the solution of fracture problems remains computationally challenging mainly due to the non-convexity of the total energy functional with respect to the combined unknown (phase field and displacement) fields. Understanding the effects of their coupling on convexity is crucial in order to address frequently encountered hurdles in fracture modeling (e.g., inefficient solvers and non-physical crack nucleation). In this paper, we develop convexity criteria for a wide class of PF fracture formulations. For this class of formulations, the second variation of the total energy functional is expressed in terms of Hessian matrices (evaluated at individual material points). Depending on the choice of geometric crack functions and degradation functions, we classify the formulations into three categories and analytically study each one separately. To study the sign of the second variation, we derive inequalities which are satisfied at material points when the Hessian matrix is locally positive semi-definite. These inequalities provide objective criteria for comparing degradation functions. Finally, the applicability of the proposed convexity criteria is demonstrated in the context of a one-dimensional problem, solved using a conventional monolithic solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Resolving electrochemically triggered topological defect dynamics and structural degradation in layered oxides

Understanding topological defects-controlled structural degradation of layered oxides—a key cathode material for high-performance lithium-ion batteries—plays a critical role in developing next-generation cathode materials. Here, by constructing a nanobattery in an electron microscope enabling atomic-scale monitoring of electrochemcial reactions, we captured the electrochemically driven atomistic dynamics and evolution of dislocations—a most important topological defect in material. We deciphered how dislocations nucleate, move, and annihilate within layered cathodes at the atomic scale. Specifically, we found two types of dislocation configurations, i.e., single dislocations and dislocation dipoles. Both pure dislocation glide/climb and mixed motions were captured, and the dislocation glide and climb velocities were first experimentally measured. Moreover, dislocation activity-mediated structural degradation such as crack nucleation, phase transformation, and lattice reorientation was unraveled. Our work provides deep insights into the atomistic dynamics of electrochemically driven dislocation activities in layered oxides.

Wang, Chunyang (ORCID:0000000184613952)↗
NASA NTRS1982

Cyclic-inelastic deformation and fatigue resistance of notched-thin aluminum plates

The results of experiments designed to develop data to assess the accuracy and utility of the critical location concept in applications of fatigue-crack nucleation analysis at notch roots are presented and discussed. Fully reversed and nonzero mean-stress data are presented over a range of lives which encompass both elastic and inelastic deformations for thin-notched specimens and smooth specimens made of 2024 T351 aluminum-alloy sheet. Notch-root strains were measured via an extensometer, whereas the formation of small cracks was detected via an eddy-current transducer. Data reported indicated the validity of the assumption that smooth and notched specimens form cracks at the same cycle number when identical deformation histories are imposed at their respective critical locations. They also serve to demonstrate the accuracy and utility of the critical-location approach in analysis to predict the formation of small cracks at notches in coupons and components.

Leis, B. N.↗
DOE OSTI2020

Role of Material Condition in Precursor Corrosion Damage and Stress Corrosion Crack Initiation Behavior of Alloy 600 in PWR Primary Water

Stress corrosion crack (SCC) initiation testing has been performed on a 15% cold-worked UNS N06600 (Alloy 600) heat in the mill-annealed (MA), solution annealed (SA), and thermally treated (TT) conditions to assess the role of grain boundary (GB) carbides on stress-assisted intergranular attack (IGA) and short crack nucleation and growth. The SCC initiation tests were conducted at constant load equivalent to materials yield stress in 360oC simulated PWR primary water with dissolved hydrogen at the Ni/NiO stability line. Results revealed the highest SCC initiation susceptibility in the Alloy 600 MA material followed by the TT and SA materials, suggesting the GB carbide distribution did not have a controlling effect on SCC initiation resistance. Quantitative assessments were conducted on IGA and short cracks to help understand this phenomenon and the role of GB carbides on precursor damage development that lead to differences in macroscopic SCC initiation behavior are discussed.

Zhai, Ziqing↗
DOE OSTI2025

Towards robust surrogate models: Benchmarking machine learning approaches to expediting phase field simulations of brittle fracture

Data-driven approaches have the potential to make modeling complex, nonlinear physical phenomena significantly more computationally tractable. For example, computational modeling of fracture is a core challenge where machine learning techniques have the potential to provide a much needed speedup that would enable progress in areas such as multi-scale modeling and uncertainty quantification. Currently, phase field modeling (PFM) of fracture is one such approach that offers a convenient variational formulation to model crack nucleation, branching and propagation. To date, machine learning techniques have shown promise in approximating PFM simulations. While standard fracture benchmarks represent realistic scenarios frequently observed in practice, they typically do not provide sufficiently challenging tests for data-driven methods. Here, to address this gap, we introduce a challenging dataset based on PFM simulations designed to benchmark and advance ML methods for fracture modeling. This dataset includes three energy decomposition methods, two boundary conditions, and 1000 random initial crack configurations for a total of 6000 simulations. Each sample contains 100 time steps capturing the temporal evolution of the crack field. Alongside this dataset, we also implement and evaluate Physics Informed Neural Networks (PINN), Fourier Neural Operators (FNO), and UNet models as baselines, and explore the impact of ensembling strategies on prediction accuracy. With this combination of our dataset and baseline models drawn from the literature we aim to provide a standardized and challenging benchmark for evaluating machine learning approaches to solid mechanics. Our results highlight both the promise and limitations of popular current models, and demonstrate the utility of this dataset as a testbed for advancing machine learning in fracture mechanics research.

Benchmark dataset↗
DOE OSTI2021

Thermomechanical Microstructural Predictions of Fracture Nucleation of Zircaloy-4 Alloys With δ and ε Hydride Distributions

Abstract A crystalline dislocation-density formulation that was incorporated with a nonlinear finite-element (FE) method was utilized to understand and to predict the thermomechanical behavior of an hexagonal closest packed (h.c.p.) zircaloy system with hydrides with either face-centered cubic (f.c.c.) or body-centered cubic (b.c.c.) hydrides. This formulation was then used with a recently developed fracture methodology that is adapted for finite inelastic strains and multiphase crystalline systems to understand how different microstructurally based fracture modes nucleate and propagate. The interrelated microstructural characteristics of the different crystalline hydride and matrix phases with the necessary orientation relationships (ORs) have been represented, such that a detailed physical understanding of fracture nucleation and propagation can be predicted for the simultaneous thermomechanical failure modes of hydride populations and the matrix. The effects of volume fraction, morphology, crystalline structure, and orientation and distribution of the hydrides on simultaneous and multiple fracture modes were investigated for radial, circumferential, and mixed distributions. Another key aspect was accounting for temperatures changes due to the effects of thermal conduction and dissipated plastic work and their collective effects on fracture. For hydrided aggregates subjected to high temperatures, thermal softening resulted in higher ductility due to increased dislocation-density activity, which led to higher shear strain accumulation and inhibited crack nucleation and growth. The predictions provide validated insights into why circumferential hydrides are more fracture-resistant than radial hydrides for different volume fractions and thermomechanical loading conditions.

Engineering↗
NASA NTRS2005

Prediction of Fretting Crack Location and Orientation in a Single Crystal Nickel Alloy

Fretting is a structural damage mechanism arising between two nominally clamped surfaces subjected to an oscillatory loading. A critical location for fretting induced damage has been identified at the blade/disk and blade/damper interfaces of gas turbine engine turbomachinery and space propulsion components. The high- temperature, high-frequency loading environment seen by these components lead to severe stress gradients at the edge-of-contact that could potentially foster crack growth leading to component failure. These contact stresses drive crack nucleation in fretting and are very sensitive to the geometry of the contacting bodies, the contact loads, materials, temperature, and contact surface tribology (friction). Recently, a high-frequency, high-temperature load frame has been designed for experimentally investigating fretting damage of single crystal nickel materials employed in aircraft and spacecraft turbomachinery. A modeling method for characterizing the fretting stresses of the spherical fretting contact stress behavior in this experiment is developed and described. The calculated fretting stresses for a series of experiments are then correlated to the observed fretting damage. Results show that knowledge of the normal stresses and resolved shear stresses on each crystal plane can aid in predicting crack locations and orientations.

Matlik, J. F.↗
DOE OSTI2024

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Characterization of damage processes in Montney siltstone under triaxial compression using acoustic emission and diagnostic imaging

SUMMARY Crack nucleation and rock failure processes in a fine-grained siltstone (Montney Formation) under triaxial compression are investigated using combined diagnostic techniques, including ultrasonic-wave measurement, acoustic-emission (AE) monitoring, computed tomography (CT) scanning, and thin-section imaging. The sample displays a weak-to-moderate inherent seismic anisotropy and noticeable stress-induced anisotropy prior to failure. No AE event was detected until the applied axial stress reached 95 per cent of the peak value. The signal-to-noise ratio is relatively low, however, and detectable AE events are more diffuse than those observed in highly brittle rocks. The AE locations correlate with a shear fracture zone imaged by CT scanning. AE moment-tensor analysis reveals that events with larger relative magnitudes are characterized by high volumetric (tensile or compressive) components, and the initiation of the failure zone is dominated by combined shear-tensile (ST) failure. Stress inversion of the AE events with high tensile components is in good agreement with the known applied stress. Microscopic imaging of thin sections from the failed sample shows that the failure zone is an en echelon structure consisting of a major fracture with branching micro and minor cracks. This failure mechanism is consistent with a ST source mechanism and is interpreted to be associated with the fine granular structure and mineral composition of Montney siltstone.

Jia, Suzie Q. (ORCID:0000000326496190)↗
DOE OSTI2023

Variations in Creep Performance of P92 Steel Within Its Composition Range: Influence of N Solubility

Underspecified composition ranges often lead to alloys with unpredictable mechanical performance. To better understand the changes in microstructure and mechanical performance associated with variations of key elements, three versions of P92 are formulated within, or close to, the specified allowable for N, B, and C ranges. Chromium and Si are also varied to influence N solubility. Different service conditions (i.e., temperature and stress) are explored. It is observed that >80% decrease in creep life occurs at 625 °C and 155 MPa for the highest B and N containing alloy. Multiscale characterization reveals key changes due to the trace element variation. The high B and N containing alloy forms deleterious BN precipitates with morphology that promotes crack nucleation and damage accumulation, but this alloy additionally forms higher fractions of beneficial MX precipitates. The alloy with the lowest B and N concentrations but greater C content shows the best creep performance—a consequence of the refined M 23 C 6 carbide precipitate population and the absence of large-scale inclusions or BN precipitates. Calculations of creep activation energy reveal that the high B and N containing alloy is more prone to damage accumulation which causes an early onset of accelerating creep and greater minimum creep rate.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A variational phase-field model For ductile fracture with coalescence dissipation

A novel phase-field model for ductile fracture is presented. The model is developed within a consistent variational framework in the context of finite-deformation kinematics. A novel coalescence dissipation introduces a new coupling mechanism between plasticity and fracture by degrading the fracture toughness as the equivalent plastic strain increases. The proposed model is compared with a recent alternative where plasticity and fracture are strongly coupled. Several representative numerical examples motivate specific modeling choices. In particular, a linear crack geometric function provides an “unperturbed” ductile response prior to crack initiation, and Lorentz-type degradation functions ensure that the critical fracture strength remains independent of the phase-field regularization length. In addition, the response of the model is demonstrated to converge with a vanishing phase-field regularization length. Additionally, the model is then applied to calibrate and simulate a three-point bending experiment of an aluminum alloy specimen with a complex geometry. The effect of the proposed coalescence dissipation coupling on simulations of the experiment is first investigated in a two-dimensional plane strain setting. The calibrated model is then applied to a three-dimensional calculation, where the calculated load-deflection curves and the crack trajectory show excellent agreement with experimental observations. Finally, the model is applied to simulate crack nucleation and growth in a specimen from a recent Sandia Fracture Challenge.

42 ENGINEERING↗
DOE OSTI2021

Computational study of the effect of core–skin structure on the mechanical properties of carbon nanofibers

The effect of the core–skin structure on the mechanical properties of carbon nanofibers is investigated in large-scale molecular dynamics simulations of tensile deformation of carbon nanofibers with the core–skin and homogeneous structures. Contrary to an established notion of the deleterious effect of the skin layer on the strength of carbon fibers, the presence of a high-quality skin layer is found to increase both the Young’s modulus and tensile strength of the nanofiber. A detailed analysis of the fracture process indicates that the nanofiber strengthening is related to the ability of skin layer to suppress crack nucleation at the core–skin interface. Furthermore, the computational predictions suggest that the design of new approaches to carbon fiber manufacturing and processing leading to the generation of a high-quality skin layer while avoiding the introduction of structural defects at the core–skin interface may yield a significant enhancement of the mechanical properties of carbon fibers.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Additively manufactured novel Al-Cu-Sc-Zr alloy: Microstructure and mechanical properties

An in-depth understanding of microstructure and resultant properties is paramount in the design of a novel alloy system, especially for additive manufacturing (AM). The present investigation aims to characterize a prototypical AM Al alloy with great potential for structural applications. An Al-1.5Cu-0.8Sc-0.4Zr alloy designed using integrated computational material engineering was printed using the laser powder bed fusion AM process. This novel alloy shows promising combination of strength and ductility in as-built and peak-aged conditions. This improvement in the tensile properties is attributed to the presence of both coherent L1 2 Al 3 Sc/Al 3 (Sc,Zr) precipitates and Cu-rich regions. The microstructures were studied via extensive microscopy at different length scales using X-ray microscopy, scanning electron microscopy, and transmission electron microscopy. Fractography revealed that the columnar grain boundaries in as-built condition allow easy slip transfer as compared to the equiaxed grains, with the apex of the melt pool acting as the crack nucleation site. Furthermore, the peak aged condition resulted in improved strength while marginally sacrificing ductility due to precipitates decorating dislocations, grain boundaries and melt pool boundaries thus acting as obstacles to slip transfer.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Insights into the failure mode of the Chelyabinsk meteor from high-fidelity simulation

In 2013, an asteroid with a diameter of approximate 20 m broke up in the atmosphere above the Russian city of Chelyabinsk producing what is known as an airburst event. The event was recorded on video from multiple angles and a 500-kg fragment was recovered from Lake Chebarkul shortly thereafter. These aspects, coupled with the asteroid’s large size, distinguish the Chelyabinsk event as unique. In this study, we use the smoothed particle hydrodynamics method to examine the breakup mode of a Chelyabinsk-sized monolithic asteroid and probe whether a monolithic structure is consistent with observations from the Chelyabinsk event. We model the object as a homogeneous sphere with properties consistent with Chelyabinsk meteorites. We use planar-2D simulations to test the sensitivity of results to variations in the material model and to qualitatively inspect the cause of crack nucleation and fragmentation. We then compare results from a 3D simulation to observations from the event. Our 3D simulation predicts a burst height within 3 km of the observed peak luminosity suggesting the responsible asteroid may have been monolithic.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

In situ X-ray imaging of pore formation mechanisms and dynamics in laser powder-blown directed energy deposition additive manufacturing

Directed energy deposition (DED) additive manufacturing (AM) is receiving growing attention in many applications, such as repair, remanufacturing, and fabrication of functionally graded structures. However, the laser-matter interactions and melt pool dynamics in laser DED with powder flow are still unclear, particularly in how pores form and flow inside the melt pool during the process. Understanding the porosity formation mechanisms is critical in the qualification, certification, and overall properties of a DED AM part. Porosity is a common phenomenon and can significantly hinder the quality of DED fabricated parts, as the pores can act as sites of crack nucleation and propagation. Here in this paper, we reveal four types of pore formation mechanisms through in-situ and operando high-speed high-resolution X-ray imaging in the DED AM process. Our results confirm that porosity within the feedstock powder induces pores in the process. We also observed pore formation mechanisms unique to the laser-based, powder-blown DED process as a result of powder delivery, keyhole dynamics, melt pool dynamics, and shield gas. High-speed X-ray images provide direct evidence for pore formation mechanisms and show that the pores related to the interaction between the delivered powder and melt pool are the largest in size in laser-based powder-blown DED AM. These results will guide porosity mitigation, elimination, and control in DED AM.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On formulations for modeling pressurized cracks within phase-field methods for fracture

Over the past few decades, the phase-field method for fracture has seen widespread appeal due to the many benefits associated with its ability to regularize a sharp crack geometry. Along the way, several different models for including the effects of pressure loads on the crack faces have been developed. Further, this work investigates the performance of these models and compares them to a relatively new formulation for incorporating crack-face pressure loads. It is shown how the new formulation can be obtained either by modifying the trial space in the traditional variational principle or by postulating a new functional that is dependent on the rates of the primary variables. The key differences between the new formulation and existing models for pressurized cracks in a phase-field setting are highlighted. Model-based simulations developed with discretized versions of the new formulation and existing models are then used to illustrate the advantages and differences. In order to analyze the results, a domain form of the J-integral is developed for diffuse cracks subjected to pressure loads. Results are presented for a one-dimensional cohesive crack, steady crack growth, and crack nucleation from a pressurized enclosure.

42 ENGINEERING↗
DOE OSTI2022

Transition from sub-Rayleigh anticrack to supershear crack propagation in snow avalanches

Snow slab avalanches, characterized by a distinct, broad fracture line, are released following anticrack propagation in highly porous weak snow layers buried below cohesive slabs. The anticrack mechanism is driven by the volumetric collapse of the weak layer, which leads to the closure of crack faces and to the onset of frictional contact. Here, on the basis of snow fracture experiments, full-scale avalanche measurements and numerical simulations, we report the existence of a transition from sub-Rayleigh anticrack to supershear crack propagation. This transition follows the Burridge–Andrews mechanism, in which a supershear daughter crack nucleates ahead of the main fracture front and eventually propagates faster than the shear wave speed. Furthermore, we show that the supershear propagation regime can exist even if the shear-to-normal stress ratio is lower than the static friction coefficient as a result of the loss of frictional resistance during collapse. This finding shows that snow slab avalanches have fundamental similarities with strike-slip earthquakes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗