In the present work a message passing interface (MPI) parallel implementation of an optimization-based approach for the simulation of underground flows in large discrete fracture networks is proposed. The software is capable of parallel execution of meshing, discretization, resolution, and postprocessing of the solution. We describe how optimal scalability performances are achieved combining high efficiency in computations with an optimized use of MPI communication protocols. Also, a novel graph-topology for communications, called the multi-Master approach, is tested, allowing for high scalability performances on massive fracture networks. Strong scalability and weak scalability simulations on random networks counting order of 105 fractures are reported. [ABSTRACT FROM AUTHOR]