How to introduce an initial crack in phase field simulations to accurately predict the linear elastic fracture propagation threshold?

Variational phase field fracture models are now widely used to simulate crack propagation in structures. A critical aspect of these simulations is the correct determination of the propagation threshold of pre-existing cracks, as it highly relies on how the initial cracks are implemented. While prior studies briefly discuss initial crack implementation techniques, we present here a systematic investigation. Various techniques to introduce initial cracks in phase field fracture simulations are tested, from the crack explicit meshing to the replacement by a fully damaged phase field, including different variants for the boundary conditions. Our focus here is on phase field models aiming to approximate, in the $\Gamma$-convergence limit, Griffith quasi-static propagation in the framework of Linear Elastic Fracture Mechanics. Therefore, a sharp crack model from classic linear elastic fracture mechanics based on Griffith criterion is the reference in this work. To assess the different techniques to introduce initial cracks, we rely on path-following methods to compute the sharp crack and the phase field smeared crack solutions. The underlying idea is that path-following ensures staying at equilibrium at each instant so that any difference between phase field and sharp crack models can be attributed to numerical artifacts. Thus, by comparing the results from both models, we can provide practical recommendations for reliably incorporating initial cracks in phase field fracture simulations. The comparison shows that an improper initial crack implementation often requires the smeared crack to transition to a one-element-wide phase band to adequately represent a displacement jump along a crack. This transition increases the energy required to propagate the crack, leading to a significant overshoot in the force-displacement response. The take-home message is that to predict the propagation threshold accurately and avoid artificial toughening; the crack must be initialized either setting the phase field to its damage state over a one-element-wide band or meshing the crack explicitly as a one-element-wide slit and imposing the fully cracked state on the crack surface.