Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase field fracture”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Connecting In Situ Stress and Wellbore Deviation to Near-Well Fracture Complexity Using Phase-Field Simulations

The interactions among in situ stress, rock fabric, wellbore geometry, natural fractures, and other natural or man-made defects create highly complex fracture trajectories in the near-wellbore region, far more intricate than those in the far-field. These near-wellbore complexities are critical for the Utah FORGE project and Enhanced Geothermal Systems (EGS) in general. Frictional pressure loss in the near-wellbore region during stimulation can significantly influence the growth of far-field fractures, while pressure losses during circulation serve as a major source of energy dissipation. Near-wellbore fracture complexities are often observable through image logs, offering valuable insights into in situ stress characteristics. However, leveraging this information requires a high-fidelity model capable of capturing the interplay among the diverse factors influencing fracture behavior.

58 GEOSCIENCES↗
DOE OSTI2026

Predicting fracture behavior in single crystal nickel using a coupled crystal plasticity phase field damage approach

Understanding the fracture behavior of single crystal metals is critical for predicting material performance under mechanical loading. Here, in this study, we investigate the fracture characteristics of single crystal nickel tensile bars using a crystal plasticity coupled phase field damage (CP-PFD) model. Experimental tensile tests were conducted on 15 specimens spanning five crystallographic orientations and three thickness variants per orientation. The results revealed two distinct fracture modes: brittle fractures with 45-degree angled surfaces and ductile fractures characterized by significant necking. The CP-PFD model successfully replicated these fracture behaviors, demonstrating strong agreement with experimental observations. The model effectively predicted the strain at which necking and fracture occurred, as well as the orientation-dependent fracture mechanisms. By comparing experimental and simulated fracture surfaces, we establish the CP-PFD model as a robust tool for predicting single crystal behavior and damage evolution. This work provides insight into the microstructural dependence of fracture behavior and establishes a predictive framework for modeling orientation-dependent damage evolution in single-crystal nickel.

Crystal plasticity↗
DOE OSTI2025

Phase-field modeling and experiments of dynamic fracture in single crystal quartz

Predicting the onset and characteristics of brittle fracture is important for a wide range of engineering and geological material applications. In this paper, we study important aspects of brittle fracture in α-quartz by phase-field modeling and experiments using a top-down approach. In the modeling framework, the work term in the Griffith energy balance is replaced with internal energy contributions that represent surface energy, thermal energy, and elastic strain energy stored in defects. This allows parametrization of individual energy contributions in terms of internal state variables and keeps track of energy partitioning after the onset of fracture. The path and history dependence of fracture is included in evolution laws for internal state variables, e.g., entropy evolution, while the energy remains a true potential. In the experimental part, dynamic compression experiments coupled with X-ray phase contrast imaging are performed on cube-like samples with a hole. In the top-down analysis, dynamic compression and three point bending experiments from the literature are simulated with the developed phase-field damage model. In conclusion, the fitted model highlights the strain rate, size, and stress state dependence of damage nucleation and evolution in single crystal α-quartz.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Coupled phase field damage and crystal plasticity analysis of intragranular fracture: The role of crystallographic orientation and voids

Damage evolution in engineering metal alloys at the grain scale exhibits significant microstructural heterogeneity and anisotropy. These heterogeneities create local hotspots for stress and strain localization, leading to void nucleation. Crystal orientation influences the active slip systems around voids, affecting lattice rotation and potentially forming discontinuities. At low triaxiality, voids may change shape due to lower stress, rotation, elongation, and coalescence. At high triaxiality, the correlation between crystal orientation and void growth rate becomes stronger, resembling the behavior observed in isolated single crystals. Therefore, understanding the effects of crystal orientation, heterogeneous strain, and defect evolution is crucial for single crystal fracture characterization. Here, in this work, a coupled phase-field damage (PFD) and crystal plasticity (CP) model is implemented within a finite element framework to analyze crystal deformation and failure. The CP method employs a dislocation density-based constitutive model, while intragranular failure is modeled using an anisotropic PFD method. The PFD model considers both the stored energy due to elastic stretching and the energy release due to defect formation and crack formation. A single crystal Al2219 with an intracrystalline spherical void is chosen to analyze fracture. The study finds that fracture propagation is strongly correlated with crystal orientations. This coupled CP-PFD model provides accurate failure prediction in crystalline materials by incorporating the effects of crystal orientations and existing voids. This study demonstrates how the local microstructure and defects influence plastic deformation and failure mechanisms in metal alloys.

Aluminum alloy↗
DOE OSTI2025

A chain stretch-based gradient-enhanced model for damage and fracture in elastomers

Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.

Elastomers↗
DOE OSTI2024

Sprain energy consequences for damage localization and fracture mechanics

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

Science & Technology - Other Topics↗
DOE OSTI2024

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

GEOS-DEV/FORGE

This dataset is a repository that provides the input deck for numerical models associated with two Utah FORGE research projects: "Closing the Loop Between In situ Stress Complexity and EGS Fracture Compexity" (Project Number: 2-2446) and "Coupled Investigation of Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior" (Project Number: 5-2428). Three kinds of simulation decks are included: (i) hydraulic fracturing simulation (HydroFrac folder), (ii) phase-field simulation (PhaseField folder), and (iii) thermo-hydro-mechanical simulation (ThermoHydroMech folder).

Cusini, Matteo [Lawrence Livermore National Labora↗
DOE OSTI2025

Material Fracturing and Failure Simulation Datasets

Fracturing is a fundamental physics phenomena with broad relevance across multiple domains, ranging from infrastructure integrity, aerospace durability, reservoir production, and seismic events. We present a diverse dataset of simulated fracture evolution and material failure generated from two numerical solvers: the phase-field method and the combined finite-discrete element method (FDEM). These solvers differ in formulation, physical fidelity, and computational efficiency. The dataset includes five materials: PBX, anisotropic shale, tungsten, aluminum, and steel. For each, phase-field simulations span 400,000 cases: 200,000 under uniaxial tension and 200,000 under biaxial tension. The computationally expensive FDEM simulations include 90,000 split evenly among PBX, shale, and tungsten under uniaxial loading. All simulations begin with randomized initial fracture patterns. Each entry includes temporal data capturing fracture propagation dynamics. This comprehensive dataset is designed to support the development of foundational or surrogate machine learning approaches for predicting material failure. While no such models are introduced here, the dataset lays a robust foundation for advancing future research and innovation in these areas.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phase-field modeling of thermally-grown oxide and damage evolution in environmental barrier coatings

Silicon carbide-based ceramic matrix composites protected by environmental barrier coatings (EBCs) present a promising materials solution for next-generation gas turbines. Developing more robust and efficient EBCs is therefore of significant technological importance. During the service in high-temperature oxidative environments, there is a thermally grown oxide (TGO) layer, spontaneously formed in the EBC system. TGO is recognized as a critical factor for the degradation and failure of EBCs, yet the detailed mechanisms of TGO growth and its effect on EBC failure remain unclear. In this study we develop a comprehensive chemo-mechano-phase-field model to simulate growth of the TGO in EBCs, factoring in creep and deformation, and especially the cracking behaviors. The volume expansion due to TGO growth and the resulting large inelastic deformation are addressed by using our recently developed, so-called incremental realization of inelastic deformation (IRID) algorithm, in combination with an adapted Hu-Chen spectral solver for elasticity. Simulations of TGO growth are performed considering different growth modes of TGOs determined mainly by the ratio of oxidant permeability in the topcoat to that in the TGO itself. Large-scale three-dimensional (3D) simulations are performed to model the formation of interconnecting vertical/channel cracks (often called ‘mud cracks’). The simulated crack morphology are in excellent agreement with the experimental observations from the literature. The simulations also provide insights into the cracking of EBCs and its dependence on the structure and constituent properties of the coating system. Furthermore, these results demonstrate the developed damage model can be a useful tool for design of more durable EBCs.

Cracking↗
DOE OSTI2026

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

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

Materials science↗
DOE OSTI2026

The Effects of Realistic Impact-processed Asteroid Topology on Atmospheric Airbursts

When small asteroids (∼10–200 m) enter Earth’s atmosphere they can break up before hitting the ground, releasing a considerable amount of energy in what is known as an atmospheric airburst. To date, numerous studies have examined how an asteroid’s strength, size, shape, and entry characteristics affect the breakup process; but few have assessed the effects of heterogenous internal topologies. We present a two-step high-fidelity simulation approach, which we use to model the atmospheric breakup of asteroids with realistic impact-processed internal topologies and their associated ground effects. In the first step, we model the asteroid as it is impacted to generate an “impact-processed” topology. We then map the resulting fractures (i.e., damage field) from that simulation to the initial state of an atmospheric breakup simulation. For both phases, we use the smoothed particle hydrodynamics code Spheral++. As a demonstrative example, we apply our new simulation pipeline to a Tunguska-scale impactor (50 m diameter, 20 km s −1 entry velocity, and 45° angle) in a case study comparing impact-processed and fully intact asteroids in both planar 2D and full 3D. For sufficiently strong asteroids, we find that impact processing generally causes the asteroid to break up at higher altitudes. We also observe that impact-processed asteroids break up into smaller fragments than objects with equivalent properties that have not been impact-processed.

Astronomy and AstroPhysics↗
DOE OSTI2024

Development of Ti 3 SiC 2 MAX phase tubular structures for solar receiver applications

Solar receiver tubes are key components of concentrating solar-thermal power (CSP) systems that harvest solar energy. For better efficiency, the Gen3 CSP receivers, which collect heat into a heat transfer fluid, require a temperature exceeding 700 °C during operation and need to perform under extreme conditions of high temperature and high thermal stress. Operators are seeking CSP designs using new high-temperature structural materials with high thermal conductivity and high creep resistance to achieve a design life of 30 years and thus help recover the plant capital cost sooner. MAX phase materials, which consist of an early transition metal element, an A-group element, and carbon or nitrogen, are expected to exhibit high creep resistance as well as high fracture toughness. Here, in this paper, we describe fabricating both (1) dense Ti 3 SiC 2 MAX phase disks and (2) short-length tubes using field-assisted sintering technology (FAST). First, the disk samples that we fabricated are fully dense and contain ≈90 % Ti 3 SiC 2 MAX phase materials and ≈10 % TiC phase materials. We determined a flexure strength of 519 ± 32 MPa by conducting a four-point bending test at room temperature with rectangular bar samples of ≈100 % density. The thermal conductivity of the Ti 3 SiC 2 MAX phase samples, measured by light flashing analysis, decreases linearly from a value of 41 W . m -1 . K -1 at room temperature to a value of 36 W . m -1 . K -1 at 650 °C. A solar reflectance measurement of the Ti 3 SiC 2 MAX phase revealed that, temperature increases from 400 to 1400 °C, thermal emittance increases from 0.39 to 0.49, while selectivity decreases from 1.8 to 1.4, respectively. Whereas the surface oxidized MAX phase samples after 100 h exposure to air at 1000 °C exhibit that of SiC. Next, we discuss fabrication of the crack-free Ti 3 SiC 2 MAX phase tubular structures accomplished by using FAST processing in graphite bedding. A Ti 3 SiC 2 MAX phase content of > 95 % with traceable ≈3% remaining TiC phase and ≈15 % porosity were demonstrated after high-temperature annealing. An average fracture strength of ≈250 MPa was determined with Ti 3 SiC 2 MAX phase tubes of ≈85 % density by diametral compression testing at room temperature. Our work demonstrated that using FAST processing to produce Ti 3 SiC 2 MAX phase tubular structures for CSP receiver applications is a viable approach.

14 SOLAR ENERGY↗
DOE OSTI2024

Fracture Dynamics in Silicon Anode Solid-State Batteries

Solid-state batteries (SSBs) with silicon anodes could enable improved safety and energy density compared to lithium-ion batteries. However, degradation arising from the massive volumetric changes of silicon anodes during cycling is not well understood in solid-state systems. Here, we use operando X-ray computed microtomography to reveal micro- to macro-scale chemo-mechanical degradation processes of silicon anodes in SSBs. Mud-type channel cracks driven by biaxial tensile stress form across the electrode during delithiation. We also find detrimental cracks at the silicon/ solid electrolyte interface that form due to local reaction competition between neighboring domains of different sizes. Continuum phase-field damage modeling quantifies stress-driven channel cracking and shows that the lithiated silicon stress state is critical for determining the extent of interfacial fracture. This work reveals the mechanisms that govern SSBs compared to conventional lithium-ion batteries and provides guidelines for engineering chemomechanically resilient electrodes for high-energy batteries.

25 ENERGY STORAGE↗
NASA NTRS1983

Analogy between fluid cavitation and fracture mechanics

When the stresses imposed on a fluid are sufficiently large, rupture or cavitation can occur. Such conditions can exist in many two-phase flow applications, such as the choked flows, which can occur in seals and bearings. Nonspherical bubbles with large aspect ratios have been observed in fluids under rapid acceleration and high shear fields. These bubbles are geometrically similar to fracture surface patterns (Griffith crack model) existing in solids. Analogies between crack growth in solid and fluid cavitation are proposed and supported by analysis and observation (photographs). Healing phenomena (void condensation), well accepted in fluid mechanics, have been observed in some polymers and hypothesized in solid mechanics. By drawing on the strengths of the theories of solid mechanics and cavitation, a more complete unified theory can be developed.

Hendricks, R. C.↗
DOE OSTI2025

The Alamo multiphysics solver for phase field simulations with strong-form mechanics and block structured adaptive mesh refinement

Alamo is a high-performance scientific code that uses block-structured adaptive mesh refinement to solve such problems as: the ignition and burn of solid rocket propellant, plasticity, damage and fracture in materials undergoing loading, and the interaction of compressible flow with eroding solid materials. Alamo is powered by AMReX, and provides a set of unique methods, models, and algorithms that enable it to solve solid-mechanics problems (coupled to other physical behavior such as fluid flow or thermal diffusion) using the power of block-structured adaptive mesh refinement.

36 MATERIALS SCIENCE↗
NASA NTRS2009

Recent Developments in Ultra High Temperature Ceramics at NASA Ames

NASA Ames is pursuing a variety of approaches to modify and control the microstructure of UHTCs with the goal of improving fracture toughness, oxidation resistance and controlling thermal conductivity. The overall goal is to produce materials that can perform reliably as sharp leading edges or nose tips in hypersonic reentry vehicles. Processing approaches include the use of preceramic polymers as the SiC source (as opposed to powder techniques), the addition of third phases to control grain growth and oxidation, and the use of processing techniques to produce high purity materials. Both hot pressing and field assisted sintering have been used to make UHTCs. Characterization of the mechanical and thermal properties of these materials is ongoing, as is arcjet testing to evaluate performance under simulated reentry conditions. The preceramic polymer approach has generated a microstructure in which elongated SiC grains grow in the form of an in-situ composite. This microstructure has the advantage of improving fracture toughness while potentially improving oxidation resistance by reducing the amount and interconnectivity of SiC in the material. Addition of third phases, such as Ir, results in a very fine-grained microstructure, even in hot-pressed samples. The results of processing and compositional changes on microstructure and properties are reported, along with selected arcjet results.

Johnson, Sylvia M.↗