“…Subsequent efforts have been directed to extend this type of model with regards to non-matched densities [57], divergence-free mixture velocities [25], thermodynamic consistency [1,40], and flows with more than two fluids [9,15,26,27,48,49,50]. Regarding the inclusion of surfactants we refer to [30,55,56,74,76,81], all of which are restricted to two fluids. Our phase field model builds up on [33], where surfactants in two-phase flow are studied within a free energy framework.…”