Adaptive Mesh Coarsening for Efficient Phase-Field Fracture Simulations
Abstract
Phase-field models provide a versatile framework for simulating fracture nucleation and propagation without explicit crack tracking. Their principal computational challenge is the need to resolve a small regularization length with a sufficiently fine finite element mesh. This requirement can render uniform-mesh simulations prohibitively expensive, particularly for three-dimensional problems, problems involving distributed crack nucleation, and soft nearly incompressible materials. Adaptive mesh refinement offers a natural means of reducing this cost, but existing approaches often rely on heuristic refinement indicators, are primarily designed for problems containing pre-existing cracks, and retain increasingly large refined regions as cracks grow in size. This work presents an adaptive mesh coarsening framework along with adaptive mesh refinement for phase-field fracture. The coarsening strategy replaces the fractured region in the crack wake with a coarse crack band that preserves the essential mechanical behavior of a crack, making it necessary to have a refined mesh only in a small region near the tip of the growing crack. The method introduces a physics-based refinement indicator derived from the violation of the material strength surface, which is a necessary condition for fracture evolution. The indicator therefore robustly identifies regions where crack nucleation or propagation is imminent and can be applied across arbitrary materials, geometries, and loading conditions. The framework is implemented in parallel within FEniCSx; the supporting finite element codes are made available. Its generality and substantial computational benefit are demonstrated through benchmark problems involving crack propagation and nucleation under quasi-static and dynamic loading, for soft and hard materials, and for thermomechanical fracture.
Related
- Matrix-free phase-field modeling of fracture in micromechanical testing simulations of inelastic materials
- Numerical benchmark to assess sensitivity to spatial discretization in quasi-static fracture simulations : Application to phase-field models
- Adaptive Crack-Tip Localization of Phase-Field Fracture Computations
- Simulations of Particle-Laden Flows with Large Dispersed-Phase Size Disparities Using Scalable Parallel Adaptive Methods
- Efficient quantum phase estimation with adaptive entanglement-assisted Hadamard test
BibTeX
@misc{dahal2026adaptive,
title = {{Adaptive Mesh Coarsening for Efficient Phase-Field Fracture Simulations}},
author = {Aarosh Dahal and Abhinav Gupta and Aditya Kumar},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.21201},
eprint = {2609.21201},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.21201},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - UNPB
TI - Adaptive Mesh Coarsening for Efficient Phase-Field Fracture Simulations
AU - Aarosh Dahal
AU - Abhinav Gupta
AU - Aditya Kumar
PY - 2026
DA - 2026/09/18
PB - arXiv
DO - 10.48550/arxiv.2609.21201
UR - https://arxiv.org/abs/2609.21201
AB - Phase-field models provide a versatile framework for simulating fracture nucleation and propagation without explicit crack tracking. Their principal computational challenge is the need to resolve a small regularization length with a sufficiently fine finite element mesh. This requirement can render uniform-mesh simulations prohibitively expensive, particularly for three-dimensional problems, problems involving distributed crack nucleation, and soft nearly incompressible materials. Adaptive mesh refinement offers a natural means of reducing this cost, but existing approaches often rely on heuristic refinement indicators, are primarily designed for problems containing pre-existing cracks, and retain increasingly large refined regions as cracks grow in size. This work presents an adaptive mesh coarsening framework along with adaptive mesh refinement for phase-field fracture. The coarsening strategy replaces the fractured region in the crack wake with a coarse crack band that preserves the essential mechanical behavior of a crack, making it necessary to have a refined mesh only in a small region near the tip of the growing crack. The method introduces a physics-based refinement indicator derived from the violation of the material strength surface, which is a necessary condition for fracture evolution. The indicator therefore robustly identifies regions where crack nucleation or propagation is imminent and can be applied across arbitrary materials, geometries, and loading conditions. The framework is implemented in parallel within FEniCSx; the supporting finite element codes are made available. Its generality and substantial computational benefit are demonstrated through benchmark problems involving crack propagation and nucleation under quasi-static and dynamic loading, for soft and hard materials, and for thermomechanical fracture.
ER -