CFD based shape optimization of axisymmetric cavitators in supercavitating flows