Parallel solution, adaptivity, computational convergence, and open-source code of 2d and 3d pressurized phase-field fracture problems
For researchers in computational fracture mechanics, this provides an efficient open-source solver with proven convergence and scalability.
The paper presents a scalable parallel solver for phase-field fracture problems, achieving optimal linear solver scaling and convergence of functionals like crack opening displacement. Adaptive mesh refinement recovers optimal convergence rates for non-smooth solutions.
We present a scalable, parallel implementation of a solver for the solution of a phase-field model for quasi-static brittle fracture. The code is available as open source. Numerical solutions in 2d and 3d with adaptive mesh refinement show optimal scaling of the linear solver based on algebraic multigrid, and convergence of the phase-field model towards exact values of functionals of interests such as the crack opening displacement or the total crack volume. In contrast to uniform refinement, adaptive mesh refinement allows us to recover optimal convergence rates for the non-smooth solutions encountered in typical test problems. We also present numerical studies of the influence of the finite domain size on functional evaluations used to approximate the infinite domain.