A numerical procedure for the calculation of solidification and melting phase change using PISO algorithm is presented. In case of phase change problem, the coupling between velocity/pressure/temperature and liquid fraction is important. The converged temperature and liquid fraction solution which satisfies the energy balance is acquired by applying enthalpy method into inner iteration in matrix solver. And a modified PISO algorithm version is introduced to properly solve the coupling between velocity/pressure/temperature and liquid fraction. A comparison of the proposed procedure with a standard iterative method shows improvement both in terms of computing speed and robustness.