Modeling Fatigue Crack Growth Using Symbolic Regression
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Brittle fracture is difficult to study in situ due to the speed of a growing crack and the often-catastrophic nature of failure in brittle materials. As a result, the influence of microstructural considerations, such as orientation, grain boundary locations, and strain field, on the crack path remains poorly understood. Presented in this study is a method addressing this knowledge gap, which utilizes the double-cleavage drilled compression geometry to achieve quasi-stable fracture in aluminum oxynitride (AlON). Synchrotron X-ray micro-computed tomography is used to characterize the crack shape and length, while high-energy diffraction microscopy provides information on the strains, orientations, and shapes of grains in the microstructure surrounding the crack tip. During testing, the crack grew in discrete and irregular jumps while the fracture toughness falls within reported ranges. The crack in AlON is found to have no greater tendency to crack intergranularly as compared to transgranularly, and grains which are cracked transgranularly do not display a trend in orientation or stress when compared to those around which the crack followed a grain boundary. The high resolution of the crack path and microstructural data provides a path forward for modeling and understanding 3D brittle fracture.
Solidification cracking remains a persistent barrier to laser powder bed fusion (LPBF) processing of solid-solution-strengthened (SSS) Ni-based superalloys and is commonly attributed to carbide formation, liquation-type failure, or elemental segregation. In this work, the cracking behavior of three SSS Ni-based superalloys—Inconel 625 (625), Inconel 617 (617), and Haynes 230 (230)—is examined under comparable LPBF conditions to evaluate the origins of their cracking susceptibility. Carbides were present in both 625 and 230, yet cracking behavior differed significantly. MC-type carbides in 625 remain discrete and preserve liquid connectivity, whereas Cr-rich M23C6 carbides in 230 form at cellular triple junctions that bottleneck the interdendritic liquid network, influencing cracking through liquid connectivity rather than as an independent cause. Alloy 617 exhibited anomalous segregation without resolvable secondary phases yet cracked mildly, indicating that extensive carbide formation is not a prerequisite for cracking. To rationalize these observations, a transport-limited framework is proposed in which cracking occurs when terminal-liquid feeding cannot accommodate solidification-induced strain; among the mechanisms evaluated, this framework is most consistent with the physical constraints imposed by LPBF solidification. Transport-based crack-susceptibility indexes incorporating permeability, viscosity, and solidification kinetics distinguish the alloy hierarchy; permeability ratios Krat increased from approximately 7 (625) to 11.5 (617) and 14.5 (230), with corresponding increases in the μ-weighted mushy zone risk index (CSIMZRM). These results are consistent with transport-limited liquid feeding as a rate-limiting contributor to solidification cracking in LPBF-processed SSS Ni-based superalloys, providing a physically grounded basis for alloy and process design.
Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.
Surface crack detection and dimensional measurement at active mining sites present significant safety and operational challenges. Manual inspection methods are labor-intensive, spatially incomplete, and expose personnel to hazardous environments, while existing automated approaches have been developed primarily for concrete civil infrastructure and have not been validated on the complex, variable surfaces characteristic of mining environments. This dissertation presents an automated pipeline that integrates deep learning semantic segmentation with Structure-from-Motion photogrammetry to detect surface cracks and measure their aperture, length, and vertical displacement from standard RGB imagery acquired during routine Uncrewed Aerial Vehicle (UAV) survey operations, without requiring additional sensor hardware or manual measurement. The pipeline combines a U-Net architecture with an EfficientNet-B0 encoder, pretrained on the SDNET2018 concrete crack dataset and fine-tuned on a mining-specific dataset spanning laboratory concrete specimens, coal refuse impoundment embankments, and post-blast limestone quarry benches. Photogrammetric reconstruction is performed using COLMAP Structure-from-Motion and Multi-View Stereo, with crack segmentation masks projected into the reconstructed point cloud to enable three-dimensional vertical displacement measurement through local plane fitting and bimodal surface detection. The pipeline was validated across 36 controlled laboratory specimens at three imaging distances and four vertical displacement levels, achieving aperture measurement RMSE of 0.047 cm and R² of 0.954, and vertical displacement RMSE of 0.140 cm and R² of 0.966, against independent caliper measurements. Field application at a coal refuse impoundment in southwestern Pennsylvania detected 71 crack components across the embankment crest, with a dominant longitudinal crack exhibiting aperture values reaching 28 cm and a 95th percentile vertical displacement of 35.53 cm, consistent in magnitude and spatial distribution with simultaneously acquired LiDAR-derived estimates. Application across four post-blast limestone quarry bench datasets in California successfully characterized blast-induced fracture networks at ground sampling distances ranging from 0.59 to 1.23 cm/pixel, with detected crack geometries physically consistent with observable surface conditions at each site. The results demonstrate that deep learning-based crack detection and photogrammetric measurement can be integrated into routine UAV inspection workflows at mining sites, providing repeatable, scalable, and quantitative crack characterization across surface types, crack scales, and displacement magnitudes not previously addressed in the literature. The pipeline requires no dedicated surveying equipment beyond the UAV platforms already deployed at mine sites for survey and monitoring purposes, supporting practical adoption within existing operational workflows.
Nuclear power plant components operate in environments that promote multiple degradation mecha- nisms, several of which involve crack initiation and growth. An ongoing effort in the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is developing a general capability within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework for simulating three-dimensional crack growth under a range of driving conditions, including fatigue, stress corrosion cracking (SCC), brittle fracture, and stress-relaxation cracking. This report demonstrates an end-to-end workflow that uses this capability to model weld-residual-stress-driven SCC in the J-groove weld of a pressurized-water reactor control rod drive mechanism penetration in the vessel head. The workflow consists of a thermomechanical welding simulation with temperature-dependent plasticity, followed by cooldown to ambient conditions, and a restart of the simulation using the MOOSE extended finite element method (XFEM) module to propagate a three-dimensional crack through the residual stress field. New welding capabilities were developed to properly initialize newly activated elements in the weld region, and robustness improvements were made to the mesh-based algorithm for defining cutting planes in the 3D XFEM algorithm, allowing it to handle complex crack fronts and stress fields. Together these advances allowed the simulated SCC crack to grow from an initial elliptical flaw in the weld, across the weld, through the tube wall, and almost to the triple point (where the weld, tube, and reactor pressure vessel head intersect) over roughly 36 years of simulated service. These results demonstrate a workflow that can be extended to fully three-dimensional welding simulations and more complex crack interaction problems.
Abstract Single‐crystalline layered oxide (LiNi 1‐ x ‐ y Mn x Co y O 2 ) cathodes have been found to exhibit exceptional electrochemical properties when coupled with various inorganic solid electrolytes (ISEs) in all‐solid‐state batteries (ASSBs). Their advantages stem from the robust morphological integrity with the absence of grain boundaries and the high electrochemical oxidative stability. Here, ASSBs featuring single‐crystalline LiNiO 2 (LNO) with the highest Ni content are reported, offering a high theoretical specific capacity of 275 mAh g ‐1 alongside a high average discharge voltage (3.7 V vs Li + /Li). Through a careful investigation, it is demonstrated that micron‐sized single‐crystalline LNO (µSC‐LNO) composite cathodes with a halide ISE exhibit a high initial discharge capacity of 205 mAh g ‐1 with an outstanding cycle performance over 200 cycles in room‐temperature ASSBs. The significance of engineering parameters is emphasized, such as particle size and specific density, in promoting a homogeneous and fast Li + transport within the composite cathodes. Furthermore, the formation of undesirable interphase between the halide ISE in the cathode and sulfide ISE separator is elucidated, which may be a critical factor impeding long‐term cyclability of ASSBs. This work provides insights into the design of composite cathodes for high‐energy‐density ASSBs.
Ceramic materials are known for their high hardness and strength but are limited by low toughness and sudden failure. Inspired by the microstructure of dental enamel, which features undulating rods that promote crack deflection and energy absorption, architected specimens were designed and manufactured. A bespoke ball-on-ring (BoR) testing apparatus measured the static biaxial rupture strength of these bioinspired materials. Additionally, impact testing with spherical steel projectiles assessed their behavior under dynamic conditions. Monolithic alumina specimens were compared to architected alumina reinforced with yttria-stabilized zirconia (YSZ) rods produced via direct ink write 3D printing. BoR tests revealed that monolithic alumina had a Weibull modulus of 10.53 and a characteristic strength of 372.3 MPa, while the composite specimens showed a Weibull modulus of 5.46 and a characteristic strength of 213.3 MPa. Although the composite failed at lower stresses and projectile velocities, it exhibited notable crack deflection and fracture resistance, with cracks being effectively interrupted by the rod inclusions. The composite specimens also resulted in fewer and larger fragments upon impact. Finite element simulations confirmed the effectiveness of the rod structures in enhancing fracture resistance.
Fatigue analysis tools can vary across industries. For example, automotive engineers often use the Wöhler (S-N) approach to design for safe-life, while aerospace engineers prioritize damage tolerance and inspection intervals, relying instead on crack growth models such as Paris’ law. Although both approaches may deal with the control of cracks in similar materials, their analysis tools and material characterizations are fundamentally distinct. This divide mirrors the classic split between stress-based strength analysis and linear elastic fracture mechanics. However, modern nonlinear models that incorporate material softening, such as cohesive laws, blur this boundary and capture fracture behaviors across scales. This presentation describes the CF23 fatigue model, which uses cohesive softening to link S-N crack initiation with crack propagation rates. CF23 spans the full fatigue spectrum, from initial propagation transients to steady-state growth and threshold conditions, offering a unified framework that bridges Wöhler and Paris-based methodologies. Example applications include fatigue crack propagation transients in adhesive interfaces and skin/stiffener separation.
A fracture study of coupon-scale composite cylinders with embedded defects was conducted with an objective to assess and validate a modeling approach using two available cohesive material models. The study included experimental and simulation evaluations of initiation of crack growth and progression. Interrupted thermal experiments used acoustic emissions monitoring to identify the onset of crack progression during each cooling interval and ultrasonic scanning provided images of defect growth. Verification, validation, and uncertainty quantification (VVUQ) processes were performed in the assessment of the simulation predicted temperature at which crack propagation begins (quantity of interest). The Sobol sensitivity analysis identified the hoop direction elastic modulus in the carbon fiber reinforced polymer (CFRP) plies as the most influential parameter for simulations using both cohesive models, accounting for at least 70% of the variation in the temperature at crack propagation. The UQ temperature range for the Tvergaard-Hutchinson model was higher (more conservative) than the experimental acoustic measurement indicators of crack progression, while the temperature range for the Thouless-Parmigiani model enveloped the experimental data points for the primary defect size of 0.75 x 1 in. The simulations could not capture the stable crack growth indicated in the experiments. This is likely due to the models’ inability to represent anisotropic fracture toughness attributed to the structure of the orthotropic fiber weave in a woven composite laminate.
The present report follows on the cohesive fatigue damage model methodology proposed in NASA-TP-2018-219838. In that report, an empirical function describing the incremental damage due to cyclic loading was used to calculate fatigue damage within a cohesive formulation. The form of the function was developed such that, when integrated at a constant stress amplitude from no damage to failure, it produces a life versus load response that is consistent with an S-N diagram. Therefore, the parameters of the damage model could be obtained by fitting the model predictions to an S-N diagram. The finite element analyses performed demonstrate that the cohesive fatigue accumulation function provides a link between the S-N diagram that describes crack initiation, and the Paris law that characterizes the rate of crack propagation. However, when the model was proposed, it was not known whether the form of the damage accumulation function associated with a desired S-N diagram is unique and, if not, if the link between S-N and the Paris law is unique and independent of the fatigue function selected. In the effort described herein, several alternative forms of the damage function that reproduce the desired features of S-N diagrams were found and evaluated. The effects of each of these functions on the predicted parameters of the Paris law and the propagation threshold are discussed. The results indicate that the predicted exponent m of the Paris law is indeed independent of the damage accumulation function. However, different functions predict different values for the pre-factor C of the Paris law. Therefore, the proper damage accumulation function must be selected by comparison with experiments. One of the new damage accumulation functions proposed herein was found to be particularly useful for analysis because of the ease with which the model parameters can be determined with a minimal amount of experimental information. The effectiveness of the proposed methodology and damage function was demonstrated by conducting analyses of a double cantilever beam test, a mixed-mode bending test, and a three-point bending test of a skin/doubler specimen. The results indicate that the same set of model parameters can provide accurate predictions of the rate of fatigue crack propagation for a variety of material interfaces, mode mixities, load levels, and stress ratios.
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.
In this study, we investigated the effects of combined intense particle and heat flux exposure on advanced tungsten plasma-facing materials within the DIII-D fusion facility. Our test matrix included two types of dispersoid-strengthened tungsten (containing either 100 nm diameter TiO 2 or Ni particles), along with high-purity polycrystalline tungsten as a reference. This experiment relied on a sample geometry angled at 15° relative to the divertor surface, thereby allowing the surfaces to intercept steady-state perpendicular heat fluxes (q ⟂ ) ranging from 10.1 to 19.6 MW/m 2 . During each shot, the samples were exposed to 42 Hz edge-localized modes (ELMs), allowing us to test the material response to transient heating. We correlated the exposure conditions with extensive post-test surface composition analysis and microscopy to determine how the plasma modified each surface. The angled specimens closest to the strike point received the highest combined heat and particle flux and melted midway through the experiment. EBSD analysis revealed they were completely recrystallized throughout, with an average grain size >100 µm. On the other hand, the specimens that received a lower steady state heat flux survived with more superficial surface damage. Whereas the high-purity polycrystalline tungsten exhibited a higher surface roughness, the dispersoid-strengthened material exhibited more extensive shallow inter-granular cracking. In addition, the surface was depleted of dispersoids following plasma exposure, possibly because of evaporation and/or sputtering. The results described here provide insights into the performance of these materials in a fusion environment which can guide further optimization for use in long-pulse devices.
After the rollout of Space Shuttle Discovery in April 2005 in preparation for return-to-flight, there was a failure of the Orbiter (OV-103) helium signature leak test in the gaseous hydrogen (GH2) system. Leakage was attributed to the Flow Control Valve (FCV) in Main Engine 3. The FCV determined to be the source of the leak for OV-103 is designated as LV-58. The nitrile/Buna N rubber O-ring seal was removed from LV-58, and failure analysis indicated radial cracks providing leak paths in one quadrant. Cracks were eventually found in 6 of 9 FCV O-rings among the three Shuttle Orbiters, though none were as severe as those for LV-58, OV-103. Testing by EM10 at MSFC on all 9 FCV O- rings included: laser dimensional, Shore A hardness and properties from a dynamic mechanical analyzer (DMA) and an Instron tensile machine. The following test data was obtained on the cracked quadrant of the LV-58, OV-103 O-ring: (1) the estimated compression set was only 9.5%, compared to none for the rest of the O-ring; (2) Shore A hardness for the O.D. was higher by almost 4 durometer points than for the rest of the O-ring; and (3) DMA data showed that the storage/elastic modulus E was almost 25% lower than for the rest of the O-ring. Of the 8 FCV O-rings tested on an Instron, 4 yielded tensile strengths that were below the MIL spec requirement of 1350 psi-a likely influence of rubber cracking. Comparisons were made between values of modulus determined by DNA (elastic) and Instron (Young s). Each nitrile/Buna N O-ring used in the FCV conforms to the MIL-P-25732C specification. A number of such O-rings taken from shelf storage at MSFC and Kennedy Space Center (KSC) were used to generate a reference curve of DMA glass transition temperature (Tg) vs. shelf storage time ranging from 8 to 26 years. A similar reference curve of TGA onset temperature (of rubber weight loss) vs. shelf storage time was also generated. The DMA and TGA data for the used FCV O-rings were compared to the reference curves. Correlations were also made between the DMA modulus (at 22 C) and Shore A hardness for all 9 of the FCV O-rings used among the three Shuttle Orbiters. The radial cracking in the FCV O-rings was determined to be due to ozone attack, as nitrile/Buna N rubber is susceptible to such attack. Nitrile/Buna N material under MIL-P25732C should be used in a hydraulic fluid environment to help protect it from cracking. However, the FCV O-rings were used in an air only environment. The FCV design has as much as a 9-mil gap that allows the O.D. of the O-ring to be directly exposed to ozone, pressurized air and some elevated temperatures, accelerating the weathering process that leads to O-ring cracking. Space Shuttle flights will likely not continue past 2010. Therefore, Shuttle management decided to continue using the nitrile/Buna N material for the FCVs, but have each O-ring replaced after 3 years to minimize any chances for crack initiation.
Conventional solid electrolyte interphases (SEIs) strongly adhere to micro-silicon (µ-Si) and crack under volume changes, causing poor cycling performance. Nano-silicon improves cycling performance but remains costly with limited calendar life. Here potentiostatic ageing tests demonstrate that both calendar and cycle ageing are governed by SEI cracking and dissolution with different relative contributions. When the system is not dominated by SEI dissolution, the relative calendar life of Si anodes could correlates positively with their cycle life. LiF-rich SEI that enables long cycle life in µ-Si is therefore expected to enhance calendar life as well. Using this framework, we screened electrolytes, SEIs and electrodes and validated them with full-cell storage. LiF-rich SEI minimizes cracking and dissolution, enabling μ-Si to achieve excellent calendar life, whereas nano-silicon suffers from SEI dissolution and needs reduced electrolyte–electrode contact for better calendar life. This work clarifies calendar-ageing behaviour and accelerates electrolytes and SEI development for long-life Si anodes.
Ni-based structural alloys in molten salt environments often experience simultaneous mechanical loading and corrosive attack, yet the mechanisms governing stress-corrosion interactions remain unclear. Prior studies largely emphasize tensile stress, while the role of compressive stress has received limited attention. Here, reactive molecular dynamics simulations are used to investigate the coupled effects of applied strain and corrosion in Ni 0.75 Cr 0.25 exposed to molten FLiNaK at 800 °C. A Σ5(210) grain boundary model is subjected to tensile (+4%) to compressive (−4%) uniaxial strains, and corrosion behavior is evaluated through fluorine adsorption, charge redistribution, and grain boundary evolution. Tensile strain accelerates intergranular corrosion susceptibility by reducing local atomic packing through elastic dilation and increasing excess free volume at the grain boundary, which enhances atomic mobility and salt infiltration. In contrast, compressive strain can suppress corrosion by promoting the formation of a ridge-like surface layer along the grain boundary, limiting salt access to the underlying alloy. These results provide atomistic insight into how stress states influence grain boundary corrosion in molten salts.
Following the successful Artemis I skip reentry on December 11, 2022, unexpected char liberation from the Orion heat- shield prompted formation of an Anomaly Response Team to determine root cause and establish flight rationale and corrective actions for subsequent missions. Through a comprehensive investigation including detailed hardware analysis, modeling & simulation, and ground testing, the team determined that the low permeability Avcoat experienced extreme internal gas pressure buildup during the skip entry that could not adequately outgas, leading to crack formation and char liberation. Key discoveries included significant performance differences between permeable and impermeable heatshield regions and successful replication of the anomaly through ground testing. Based on extensive ground testing and analysis, Artemis II flew a modified trajectory without the skip entry to minimize crack-inducing conditions, and successfully splashed down on April 10, 2026 with significantly reduced char loss. For Artemis III and beyond, a corrective action was implemented to produce a more permeable version of Avcoat that will substantially reduce internal pressure accumulation during entry.