Computation of sharp phase boundaries by spreading: the planar and spherically symmetric cases