Development of a PISO-SPH method for computing incompressible flows