Parallel meshing, discretization and computation of flow in massive Discrete Fracure Networks

In the present work an 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 post-processing of the solution. We describe how optimal scalability performances are achieved combining high efficiency in computations to an optimized use of MPI communication protocols. Also, a novel graph-topology for communications, called 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 10^5 fractures are reported.