A mesh-less method for solving incompressible, multi-phase flow problems has been developed and is discussed along with the presentation of benchmark results showing good agreement with theoretical and experimental results. Results of a systematic, parametric study of the single phase flow around a 2D circular cylinder at Reynolds numbers up to 1000 are presented and discussed. Simulation results show good agreement with experimental results. Extension of the method to deal with multiphase flow including liquid-to-vapor phase transition along with applications to cavitating flow are discussed. Insight gleaned from numerical experiments of the cavity closure problem are discussed along with recommendations for additional research. Several conclusions regarding the use of the method are made.