In this paper we derive rates of convergence for regularizations of the multidimensional two-phase Stefan problem and use the regularized problems to define backward-difference in time and C° piecewise-linear in space Galerkin approximations. We find an L2 rate of convergence of order \^ in the e-regularization and an L rate of convergence of order (h2/e + Ai/ \6F) in the Galerkin estimates whi...