Numerical Methods for Phase Change