Shape optimization for delay of laminar-turbulent transition