A double-layer reduced model for fault flow on slipping domains with hybrid finite volume scheme

In this work we are interested in dealing with single-phase flows in fractured porous media for underground processes. We focus our attention on domains where the presence of faults, with thickness several orders of magnitude smaller than other characteristic sizes, can allow one part of the domain to slide past to the other. We propose a mathematical scheme where a reduced model for the fault flows is employed yielding a problem of co-dimension one. The hybrid finite volume method is used to obtain the discretized problem, which employs two different meshes for each fault, one associated with the porous-medium domain on each side of the fault. These two meshes can move with the corresponding domain, resulting in non-matching grids between the two parts of the fault. In an earlier paper a mathematical scheme was proposed where the numerical discretization considers the hybrid finite volume method. In this paper we focus on the well-posedness of the continuous problem, the convergence of the discretized problem, and with several numerical tests we support the theoretical findings.