Shape optimization of open-boundary axisymmetric induction devices