In this paper a Lagrangian finite element approach for the simulation of incompressible flows is presented, based on the so-called Particle Finite Element Method (PFEM). The spatial discretization and the definition of the boundary terms are discussed in detail with a specific focus on free-surface flows. Additionally, some problems that can arise from the use of such a method are pointed out. Some numerical examples are given and discussed in the last section of the paper