A stabilized finite element method for shape optimization in low Reynolds number flows