A new numerical scheme is proposed for flow computation in complex discrete fracture networks. The method is based on a three-field domain decomposition framework, in which independent variables are introduced at the interfaces generated in the process of decoupling the original problem on the whole network into a set of fracture-local problems. A PDE-constrained formulation is then used to enforce compatibility conditions at the interfaces. The combination of the three-field domain decomposition and of the optimization based coupling strategy results in a novel method which can handle non-conforming meshes, independently built on each geometrical object of the computational domain, and ensures local mass conservation property at fracture intersections, which is of paramount importance for hydro-geological applications. An iterative solver is devised for the method, suitable for parallel implementation on parallel computing architectures.
A three-field based optimization formulation for flow simulations in networks of fractures on non-conforming meshes / Berrone, S; Grappein, D; Pieraccini, S; Scialò, S. - In: SIAM JOURNAL ON SCIENTIFIC COMPUTING. - ISSN 1064-8275. - (In corso di stampa).
|Titolo:||A three-field based optimization formulation for flow simulations in networks of fractures on non-conforming meshes|
|Data di pubblicazione:||Being printed|
|Appare nelle tipologie:||1.1 Articolo in rivista|
File in questo prodotto: