Esta tesis trata sobre la implementación de una nueva herramienta para resolver problemas fluidodinámicos tratando, en particular, los problemas de interfaces fluido-fluido y sólido-fluido. Esta herramienta está basada en el "Semi-Lagrangian Particle Finite Element Method" (SL-PFEM) para resolver las ecuaciones de Navier-Stokes. De manera similar al "Particle Finite Element Method" (PFEM) y su segunda versión (PFEM-2), partículas Lagrangianas son usadas para resolver el problema de convección mientras que una malla fija se usa para resolver la parte elíptica de las ecuaciones de Navier-Stokes. Uno de los objetivos es que la herramienta debe ser capaz de resolver los problemas de manera eficiente. Para ello, primero se desarrolla el esquema numérico para que sea un esquema de segundo orden tanto en espacio y tiempo y, luego, se optimiza la implementación para obtener un buen rendimiento en entornos de supercomputación.Para resolver problemas de interfaces como los de fluido-fluido, solido-fluido o la combinación de ambas, un método de enriquecimiento es aplicado al esquema numérico. Esto, junto con una variación del método "Level-Set" aplicado al SL-PFEM permite capturar y calcular de forma precisa estos tipos de problemas. Además, esta implementación lleva a la capacidad de resolver problemas con más de un cuerpo que puedan interactuar con varios fluidos.