Fracture in mode I using a conserved phase-field model.

We present a continuum phase-field model of crack propagation. It includes a phase-field that is proportional to the mass density and a displacement field that is governed by linear elastic theory. Generic macroscopic crack growth laws emerge naturally from this model. In contrast to classical continuum fracture mechanics simulations, our model avoids numerical front tracking. The added phase-field smooths the sharp interface, enabling us to use equations of motion for the material (grounded in basic physical principles) rather than for the interface (which often are deduced from complicated theories or empirical observations). The interface dynamics thus emerges naturally. In this paper, we look at stationary solutions of the model, mode I fracture, and also discuss numerical issues. We find that the Griffith's threshold underestimates the critical value at which our system fractures due to long wavelength modes excited by the fracture process.