A gas-liquid two-phase model for the simulation of a power-law fluid mold filling process with the consideration of phase change is proposed, in which the governing equations for the melt and air in the cavity, including the mass conservation, momentum conservation and energy conservation equations, are unified into one system of equation. A revised Enthalpy method, which can be used for both the melt and air in the mold cavity, is proposed to describe the phase change during the mold filling. Finite volume method on non-staggered grid is used to solve the system. The level set method is used to capture the interface evolution during the mold filling process. The interface evolution and the distributions of physical quantities such as velocity, pressure and temperature and so on are given. The “frozen skin” layers under different temperature and velocities conditions are discussed in detail. Numerical results show that increasing the temperatures of the melt and cavity is a better way to get rid of the “frozen skin” layer than increasing the injection velocity.