We present a numerical framework for the simulation of three-dimensional multi-fluid flows based on a finite element/level-set approach. The method allows a full Eulerian ” tracking ” of the interfaces between the fluids, and the properties of the interfaces can be directly taken into account as surface forces. The resolution of the fluid equations and the advection of the interface can be easily decoupled, which enables the use of efficient solving strategies. We also present a 3D benchmark of the rise of a drop in a viscous fluid. We use two different setups and compare our results to previous results obtained with other approaches to validate our method.