This formulation requires well defined boundary conditions; and solvers to those systems are, so far, computationally expensive. The existing SPH-algorithm solves the dynamic in a Projection-Correction Scheme extrapolate the ˆ uk+1, the velocity-projection, without pressure term and later correct it via solving the p as an explicit function of density. Moreover, since field variables is approxi...