The phase field method uses a unified approach to model nucleation and propagation of cracks. Here, cracks are represented by means of a scalar field variable, which interpolates smoothly between the values 0 (blue) and 1 (red) indicating fractured and undamaged material, respectively (see Fig. 1).