Use of algebraic dual spaces in domain decomposition methods for Darcy flow in 3D domains

In this work we use algebraic dual spaces with a domain decomposition method to solve the Darcy equations. We define the broken Sobolev spaces and their finite dimensional counterparts. A global trace space is defined that connects the solution between the broken spaces. Use of algebraic dual spaces results in a sparse, metric-free representation of the incompressibility constraint, the pressure gradient term, and on the continuity constraint between the sub domains. To demonstrate this, we solve two test cases: (i) a manufactured solution case, and (ii) an industrial benchmark reservoir modelling problem SPE10. The results demonstrate that the dual spaces can be used for domain decomposition formulation, and despite having more unknowns, requires less simulation time compared to the continuous Galerkin formulation, without compromising on the accuracy of the solution.