- Check the relation between $\int_{\Omega} \alpha(x) \mathrm{d}x$ and the crack length (the relation most likely depends on $\ell$). - Create an end checker that stop the simulation when the crack length reaches certain value.