Discrete particle method