On the numerical approach to a two-phase Stefan problem with nonlinear flux