Just fall back from simulating water surface to correctly implementing the particle levelset (PLS) method itself. The test is the classic Zalesak's disk problem which is the rigid body rotation in a constant vorticity velocity field.
The PLS seems working correctly on this simple problem.
The PLS seems working correctly on this simple problem.