Finite element approximation of the enthalpy formulation for Stefan problems on evolving surfaces

We propose, and analyse, a spatially discrete evolving surface finite element method for the approximation of the enthalpy formulation of the two-phase Stefan problem posed on an evolving surface. Our approach does not rely on mass-lumping and discrete maximum principles. We prove this numerical method is numerically stable, and prove $\mathcal{O}(\sqrt{h})$ error bounds for the temperature in the $L^2_{L^2}$ norm under minimal regularity assumptions by introducing a new projection-type operator. We complement our analysis with discussion on the implementation of this numerical method, where we propose a novel implementation that avoids errors due to numerical quadrature and which has not previously been considered in the literature even in the stationary, flat setting. We also include numerical experiments and experimental order of convergence demonstrations.