Add Continuum Geometry class for direct numerical simulations

For carrying out direct numerical simulations, it would be handy if we had a DNS geometry tailored for that purpose. It would simply be a bunch of cubic models (volume, area, etc.).

[jgostick]

I like this idea, a geometry class called Continuum which does all the stuff needed to create a ‘pore-less’ region. This could also be used to do porous continuum using different pore-scale models to get effective diffusivity of the region.

[ma-sadeghi]

import openpnm as op
import openpnm.models as mods
import openpnm.models.geometry as gmods
from openpnm.geometry import GenericGeometry


class Continuum(GenericGeometry):
    r"""
    Continuum subclass of GenericGeometry.

    This subclass is meant to be used for regions that are treated as
    continuum rather than an actual pore network.

    Parameters
    ----------
    network : GenericNetwork
        The network with which this Geometry should be associated
    pores : array_like
        The pores in the domain where this Geometry applies
    throats : array_like
        The throats in the domain where this Geometry applies
    name : str
        The name of the object, which is also used as the label where this
        geometry is defined.

    """

    def __init__(self, **kwargs):
        super().__init__(**kwargs)

        spacing = op.topotools.get_spacing(self.network).mean()

        self.add_model(propname='pore.diameter',
                       model=mods.misc.constant,
                       value=spacing)

        self.add_model(propname='pore.volume',
                       model=mods.geometry.pore_volume.cube,
                       pore_diameter='pore.diameter')

        self.add_model(propname='throat.diameter',
                       model=mods.misc.from_neighbor_pores,
                       prop='pore.diameter', mode='mean')

        self.add_model(propname='throat.length',
                       model=mods.misc.constant,
                       value=1e-50)

        self.add_model(propname='throat.cross_sectional_area',
                       model=mods.geometry.throat_cross_sectional_area.cuboid,
                       throat_diameter='throat.diameter')

        self.add_model(propname='throat.volume',
                       model=mods.geometry.throat_volume.cuboid,
                       throat_diameter='throat.diameter',
                       throat_length='throat.length')

        self.add_model(propname='throat.diffusive_size_factors',
                       model=gmods.diffusive_size_factors.cubes_and_cuboids,
                       pore_diameter="pore.diameter",
                       throat_diameter="throat.diameter")

        self.add_model(propname='throat.hydraulic_size_factors',
                       model=gmods.hydraulic_size_factors.cubes_and_cuboids,
                       pore_diameter="pore.diameter",
                       throat_diameter="throat.diameter")