Shape optimization for Navier-Stokes variational-hemivariational inequality with slip boundary condition