This study proposed a kind of optimization design for a reversible axial-flow pump based on an ordinary one-way pump. Three-dimensional (3D) Reynolds-averaged Navier−Stokes (RANS) equations was used to predict the pump performance, and the optimized design was validated by an external characteristic test. Six main geometry parameters of an impeller and diffuser based on an orthogonal experiment were set as design variables. The efficiency and head under forward and reverse design conditions were set as the optimization objective. Based on 120 groups of sample designs obtained from Latin hypercube sampling (LHS), a two-layer artificial neural network (ANN) was used to build a non-linear function with high accuracy between the design variables and optimization objective. The optimized design was obtained from 300 groups of Pareto-optimal solutions using the non-dominated based genetic algorithm (NSGA) for multiobjective optimization. After optimization, there was a slight decrease in the forward pump efficiency and head. The reverse pump efficiency and head on the other hand was largely improved and the high efficiency range was also widened.