The design of pumpjet propulsors (PJP) is addressed through a simulation-based design optimization approach built on a parametric description of the main geometrical characteristics of the system, a RANSE solver with mixing plane interface capabilities and a genetic algorithm. Both Rotor/Stator (PJP-R/S) and Stator/Rotor (PJP-S/R) configurations are considered and optimal designs from a multi-objective optimization process aimed at maximizing the propulsive efficiency at the lowest possible cavitation inception index are compared to a reference ducted propeller with decelerating nozzle (Gaggero et al., 2012), which served as baseline during the design. Detached Eddy Simulations (DES) were finally carried out to highlight, in addition to the performance improvements provided by the PJPs, also the influence of the rotor/stator/nozzle interactions on the vortical structures shed by the propulsor.