A parallel solver for large scale DFN flow simulations

Flows in fractured media have been modeled with many different approaches in order to get reliable and efficient simulations for many critical applications. The common issues to be tackled are the wide range of scales involved in the phenomenon, the complexity of the domain and the huge computational cost. In the present paper we propose a parallel implementation of the PDE-constrained optimization method presented in previous works of the authors for dealing with arbitrary Discrete Fracture Networks on non conforming grids. We show the scalability performances and the efficiency of the parallel algorithm and we also test the robustness of the method on complex and strongly connected DFN configurations, that would result very diffcult to mesh with conventional approaches relying on some kind of mesh conformity at the interfaces.