We derive a posteriori estimates for a discretization in space of the standard
Cahn–Hilliard equation with a double obstacle free energy.
The derived estimates are robust and efficient, and in practice are combined
with a heuristic time step adaptation.
We present numerical experiments in two and three space dimensions and compare
our method with an existing heuristic spatial mesh adaptation algorithm.