In this research, the propulsion of the proposed jellyfish-inspired mantle undulated propulsion robot (MUPRo) is optimized. To reliably predict the hydrodynamic forces acting on the robot, the proposed non-intrusive reduced-order model (NIROM) based on proper orthogonal decomposition (POD) additionally considers the POD basis that has an important contribution to the features on the specified boundary. The proposed model establishes a mapping between the parameter-driven motion of the mantle and the evolution of the fluid characteristics around the swimmer. Moreover, to predict new cases where the input needs to be updated, the input of the proposed model is taken from the kinematics of the robot rather than extracted from full-order high-fidelity models. In this paper, it takes approximately 950s to perform a simulation using the full-order high-fidelity model. However, the computational cost for one prediction with the proposed POD-NIROM is around 0.54s, of which about 0.2s is contributed by preprocessing. Compared with the NIROM based on the classic POD method, the proposed POD-NIROM can effectively update the input and reasonably predict the characteristics on the boundary. The analysis of the hydrodynamic performance of the MUPRo pinpoints that, under a certain period and a certain undulation amplitude, the hydrodynamic force generated by swinging-like mantle motion (k<0.5) is greater, outperforming Aequorea victoria in startup acceleration. It is demonstrated that considering a certain power loss and a certain tail beat amplitude, the wave-like mantle motion (k>0.5) can produce greater propulsion, which means higher propulsion efficiency.