A finite element formulation for solving multidimensional phase-change problems is presented. The formulation considers the temperature as the unique state variable, it is conservative in the weak form sense and it preserves the moving interface condition. In this work, a consistent Jacobian matrix that ensures numerical convergence and stability is derived. Also, a comparative analysis with other different phase-change finite element techniques is performed. Finally, two numerical examples are analyzed in order to show the performance of the proposed methodology.