The article deals with the shape optimization of a supersonic ejector for propulsion of an experimental supersonic wind tunnel. This ejector contains several primary nozzles arranged around the mixing chamber wall. CFD software Fluent was used to compute the flow in the ejector. A dynamic mesh method was applied to find an optimal shape of the three-dimensional geometry. During the work it was found out that the previously developed optimization method for subsonic ejectors must be modified. The improved method is more stable and the solution requires fewer optimization steps. The shapes of the mixing chamber, the diffuser, inlet parts and the optimal declination of the primary nozzles are obtained as the optimization results.