The optimization of water distribution networks (WDN) has evolved, requiring approaches that seek to reduce capital costs and maximize the reliability of the system simultaneously. Hence, several evolutionary algorithms, such as the non-dominated sorting-based multi-objective evolutionary algorithm (NSGA-II), have been widely used despite the high computational costs required to achieve an acceptable solution. Alternatively, energy-based methods have been used to reach near-optimal solutions with reduced computational requirements. This paper presents a method to combine the domain knowledge given by energy-based methods with an evolutionary algorithm, in a way that improves the convergence rate and reduces the overall computational requirements to find near-optimal Pareto fronts (PFs). This method is divided into three steps: parameters calibration, preprocessing of the optimal power use surface (OPUS) results, and periodic feedback using OPUS in NSGA II. The method was tested in four benchmark networks with different characteristics, seeking to minimize the costs of the WDN and maximizing its reliability. Then the results were compared with a generic implementation of NSGA-II, and the performance and quality of the solutions were evaluated using two metrics: hypervolume (HV) and modified inverted generational distance (IGD+). The results showed that the feedback procedure increases the efficiency of the algorithm, particularly the first time the algorithm is retrofitted.