A general formula for the second-order projection method combined with the level set method is developed to simulate unsteady, incompressible multifluid flow with phase change. A subcell conception is introduced in a modified mass transfer model to accurately calculate the mass transfer across the interface. Teh third-order essentially nonoscillatory(ENO) scheme an..