Performances of the Mixed Virtual Element Method on Complex Grids for Underground Flow

The numerical simulation of physical processes in the underground frequently entails challenges related to the geometry and/or data. The former are mainly due to the shape of sedimentary layers and the presence of fractures and faults, while the latter are connected to the properties of the rock matrix which might vary abruptly in space. The development of approximation schemes has recently focused on the overcoming of such difficulties with the objective of obtaining numerical schemes with good approximation properties. In this work we carry out a numerical study on the performance of the Mixed Virtual Element Method (MVEM) for the solution of a single-phase flow model in fractured porous media. This method is able to handle grid cells of polytopal type and treat hybrid dimensional problems. It has been proven to be robust with respect to the variation of the permeability field and of the shape of the elements. Our numerical experiments focus on two test cases that cover several of the aforementioned critical aspects.