I guess you can’t avoid the diffusion - this stuff exactly forms your interface. That’s why it’s called diffusion interface. Moreover, to specify diffusion to be equal zero you need to put omega=1/2 which is explicitly unstable.
I now solved the issue by simply not using two fluids: I only need advection, so just saving and importing the velocity field proved to be easier…
I now let particles move step by step on the velocity field. It’s not perfect, but it allows me to directly use the velocity field more easily, which opens up some new options for me.
What I do to get a 1024x1024x1024 velocity field without filling up my disk (or having to code compression ):