Calculus with combinatorial differential forms for fluid flow analysis in porous and fractured media

The fabric of porous and fractured media contains solid regions (grains) and voids. The space conducting fluids is a system of connected voids with variable geometries. Relative to the grain sizes, the voids can be voluminous with three comparable large extensions, narrow expansive with two comparable large extensions and one smaller extension, and thin long with one comparable large extension and two smaller extensions. The widely used representation of void spaces by systems of spheres connected by cylinders (pore network models) is an acceptable approximation for some special cases, but not for most porous and fractured media. We propose a flexible method for modelling such media by mapping their measured fabric's characteristics - void and grain volume distributions and shapes - onto polyhedral tessellations of space. The map assigns voluminous voids and grains to polyhedrons (3D), narrow expansive voids to some polyhedral faces (2D), and thin long voids to some polyhedral edges (1D), as dictated by experimental data. The analysis of transport through such discrete structures with components of different dimensions is performed by a novel mathematical method, which uses combinatorial differential forms to represent physical properties and their fluxes, as well as structure-preserving operators on such forms to formulate the conservation laws exactly and directly in matrix form, ready for computation. The method allows for individual material properties, such as conductivity, to be assigned to voids of all dimensions, so that the three types of voids are suitably represented. Publicly available XCT images of four different rocks are used to test the method.