A multigrid method for steady incompressible Navier-Stokes equations based on partial flux splitting