projection method used on incompressible Navier-Stokes (with boundary force)