Simulation of fluid flows of multi-materials is an intriguing topic in Computational Mechanics. Capturing the physics of the interface between different materials poses a challenge because of the discontinuities that may occur on the interface. Several methods have been proposed in the literature to deal with this issue. In this paper, a technique based on Nitsche's Method has been employed on a fixed mesh combined with the PFEM-2 strategy for the solution of Navier-Stokes Equations on multi-fluid flows. The novelty of this technique is its capability of capturing the strong and weak discontinuities, and its compatibility for the application of various types of boundary conditions on the interface.