The Cubic Interpolated Pseudo-Particle Navier Stokes equation was applied to investigate the two-dimensional laminar square lid driven cavity flow of water at Reynolds number 1000. CIP-NSE scheme was used to solve hyperbolic term of the vorticity transport equation. In the CIP-NSE, the gradient and the value of the vorticity at the nodes is determined and the stream function is then determined using the vorticity equation. The Runge-Kutta method has been used to calculate the velocity and position of the particle with the effects of drag force and gravitational forces. The hard sphere model has been applied to show the collisions effect on particles in the lid-driven cavity. It is discovered that the numerical simulation of CIP-NSE provided a very good agreement with the established benchmark results by previous researchers.