This paper investigates control and design-for-control strategies to improve the resilience of sectorized water distribution networks (WDN), while minimizing pressure induced pipe stress and leakage. Both evolutionary algorithms (EA) and gradient-based mathematical optimization approaches are investigated for the solution of the resulting large-scale non-linear (NLP) and bi-objective mixed-integer non-linear programs (BOMINLP). While EAs have been successfully applied to solve discrete network design problems for large-scale WDNs, gradient-based mathematical optimization methods are more computationally efficient when dealing with large search spaces associated with continuous variables in optimal network control problems. Considering the advantages of each method, we propose a sequential hybrid method for the optimal design-for-control of large-scale WDNs, where a refinement stage relying on gradient-based mathematical optimization is used to solve continuous optimal control problems corresponding to design solutions returned by an initial EA search. The proposed method is applied to compute the Pareto front of a bi-objective design-for-control problem for the operational network BWPnet, where we consider reopening closed connections between isolated supply areas. The results show that the considered design-for-control strategy increases the resilience of BWPnet while minimizing pressure induced leakage. Moreover, the refinement stage of the proposed hybrid method efficiently improves the coarse approximation computed by the initial EA search, returning a continuous and even Pareto front approximation.