Numerical method for multi-phase inverse Stefan design problems