Hybrid numerical scheme for nonlinear two-dimensional phase-change problems with the irregular geometry