Adaptive regularization, linearization, and discretization and a posteriori error control for the two-phase Stefan problem

We consider in this paper the time-dependent two-phase Stefan problem and derive a posteriori error estimates and adaptive strategies for its conforming spatial and backward Euler temporal discretizations. Regularization of the enthalpy-temperature function and iterative linearization of the arising systems of nonlinear algebraic equations are considered. Our estimators yield a guaranteed and fully computable upper bound on the dual norm of the residual, as well as on the L2(L2) error of the temperature and the L2(H1) error of the enthalpy. Moreover, they allow to distinguish the space, time, regularization, and linearization error components. An adaptive algorithm is proposed, which ensures computational savings through the online choice of a sufficient regularization parameter, a stopping criterion for the linearization iterations, local space mesh refinement, time step adjustment, and equilibration of the spatial and temporal errors. We also prove the efficiency of our estimate. Our analysis is quite general and is not focused on a specific choice of the space discretization and of the linearization. As an example, we apply it to the vertex-centered finite volume (finite element with mass lumping and quadrature) and Newton methods. Numerical results illustrate the effiectiveness of our estimates and the performance of the adaptive algorithm.