English

Parallel matrix-free higher-order finite element solvers for phase-field fracture problems

Numerical Analysis 2020-05-04 v1 Numerical Analysis

Abstract

Phase-field fracture models lead to variational problems that can be written as a coupled variational equality and inequality system. Numerically, such problems can be treated with Galerkin finite elements and primal-dual active set methods. Specifically, low-order and high-order finite elements may be employed, where, for the latter, only few studies exist to date. The most time-consuming part in the discrete version of the primal-dual active set (semi-smooth Newton) algorithm consists in the solutions of changing linear systems arising at each semi-smooth Newton step. We propose a new parallel matrix-free monolithic multigrid preconditioner for these systems. We provide two numerical tests, and discuss the performance of the parallel solver proposed in the paper. Furthermore, we compare our new preconditioner with a block-AMG preconditioner available in the literature.

Keywords

Cite

@article{arxiv.2005.00331,
  title  = {Parallel matrix-free higher-order finite element solvers for phase-field fracture problems},
  author = {Daniel Jodlbauer and Ulrich Langer and Thomas Wick},
  journal= {arXiv preprint arXiv:2005.00331},
  year   = {2020}
}

Comments

20 pages

R2 v1 2026-06-23T15:14:19.014Z