A fully eulerian method for shape optimization, with application to navier-stokes flows