For debugging purposes, I created several modes to visualize the simulation.

normal

simulation cell types: fluid (dark blue), surface (light blue), and empty(white)

Fluid velocities: magnitude and direction of velocities denoted by red vectors

Marker particles: negative particles(magenta), positive particles(yellow)

Level Set Visualization: Value of sin(phi) shows that the level set is correctly a signed distact function