Discretization schemes

The discretization schemes available in DuMu^x More...

Description

The discretization schemes available in DuMu^x

Modules
	Box FV scheme
	The box method is a collocated finite volume scheme with control volumes centered at grid nodes.
	Cell-centered FV scheme
	Finite volume schemes with degrees of freedom located at grid cell centers.
	Staggered FV scheme
	A staggered finite volume scheme with degrees of freedom at cell-centers and facets.
	Finite element method
	The finite element method.

Files
file	basegridgeometry.hh
	Base class for grid geometries.
file	box.hh
	Defines a type tag and some properties for models using the box scheme.
file	ccmpfa.hh
	Properties for all models using cell-centered finite volume scheme with mpfa.
file	cctpfa.hh
	Properties for all models using cell-centered finite volume scheme with TPFA.
file	checkoverlapsize.hh
	Check the overlap size for different discretization methods.
file	elementsolution.hh
	Element solution classes and factory functions.
file	evalgradients.hh
	free functions for the evaluation of primary variable gradients inside elements.
file	evalsolution.hh
	free functions for the evaluation of primary variables inside elements.
file	extrusion.hh
	Helper classes to compute the integration elements.
file	fluxstencil.hh
	The flux stencil specialized for different discretization schemes.
file	functionspacebasis.hh
	Provides helper aliases and functionality to obtain the types and instances of Dune::Functions function space bases that underlie different discretization schemes.
file	discretization/fvgridvariables.hh
	The grid variable class for finite volume schemes storing variables on scv and scvf (volume and flux variables)
file	fvproperties.hh
	Declares properties required for finite-volume models models.
file	localview.hh
	Free function to get the local view of a grid cache object.
file	method.hh
	The available discretization methods in Dumux.
file	projector.hh
	Contains functionality for L2-projections from one function space into another, which can live both on the same or different grids of potentially different dimensionality.
file	rotationpolicy.hh
	Rotation policy for defining rotational symmetric grid geometries.
file	rotationsymmetricgridgeometrytraits.hh
	Rotation policy for defining rotational symmetric grid geometries.
file	rotationsymmetricscv.hh
	Wrapper to make a sub control volume rotation symmetric.
file	rotationsymmetricscvf.hh
	Wrapper to make a sub control volume face rotation symmetric.
file	scvandscvfiterators.hh
	Class providing iterators over sub control volumes and sub control volume faces of an element.
file	staggered.hh
	Defines a type tag and some properties for models using the staggered scheme. This scheme features degrees of freedom at the elements' centers and intersections (faces). TODO: detailed documentation and figures.
file	subcontrolvolumebase.hh
	Base class for a sub control volume.
file	subcontrolvolumefacebase.hh
	Base class for a sub control volume face.

Classes
class	Dumux::BaseGridGeometry< GV, Traits >
	Base class for all finite volume grid geometries. More...
struct	Dumux::CheckOverlapSize< discMethod >
	Check if the overlap size is valid for a given discretization method. More...
struct	Dumux::NoExtrusion
	Default implementation that performs no extrusion (extrusion with identity) More...
struct	Dumux::RotationalExtrusion< radAx >
	Rotation symmetric extrusion policy for rotating about an external axis. More...
struct	Dumux::SphericalExtrusion
	Rotation symmetric extrusion policy for spherical rotation. More...
class	Dumux::FluxStencil< FVElementGeometry, discMethod >
	The flux stencil specialized for different discretization schemes. More...
class	Dumux::FVGridVariables< GG, GVV, GFVC >
	The grid variable class for finite volume schemes storing variables on scv and scvf (volume and flux variables) More...
class	Dumux::Projector< ScalarType >
	Does an L2-projection from one discrete function space into another. The convenience functions makeProjectorPair or makeProjector can be used to create such a projection. More...
struct	Dumux::RotationSymmetricGridGeometryTraits< BaseTraits, rotPolicy >
	Traits for rotation symmetric grid geometries. More...
class	Dumux::RotationSymmetricSubControlVolume< SubControlVolume, rotationPolicy >
	Wrapper to make a sub control volume rotation symmetric. More...
class	Dumux::RotationSymmetricSubControlVolume< SubControlVolume, RotationPolicy::disc >
	Wrapper to make a sub control volume rotation symmetric. More...
class	Dumux::RotationSymmetricSubControlVolume< SubControlVolume, RotationPolicy::ball >
	Wrapper to make a sub control volume rotation symmetric. More...
class	Dumux::RotationSymmetricSubControlVolume< SubControlVolume, RotationPolicy::toroid >
	Wrapper to make a sub control volume rotation symmetric. More...
class	Dumux::RotationSymmetricSubControlVolumeFace< SubControlVolumeFace, rotationPolicy >
	Wrapper to make a sub control volume face rotation symmetric. More...
class	Dumux::RotationSymmetricSubControlVolumeFace< SubControlVolumeFace, RotationPolicy::disc >
	Wrapper to make a sub control volume face rotation symmetric. More...
class	Dumux::RotationSymmetricSubControlVolumeFace< SubControlVolumeFace, RotationPolicy::ball >
	Wrapper to make a sub control volume face rotation symmetric. More...
class	Dumux::RotationSymmetricSubControlVolumeFace< SubControlVolumeFace, RotationPolicy::toroid >
	Wrapper to make a sub control volume face rotation symmetric. More...
class	Dumux::ScvIterator< SubControlVolume, Vector, FVElementGeometry >
	Iterators over sub control volumes. More...
class	Dumux::ScvfIterator< SubControlVolumeFace, Vector, FVElementGeometry >
	Iterators over sub control volume faces of an fv geometry. More...
class	Dumux::SubControlVolumeBase< Imp, ScvGeometryTraits >
	Base class for a sub control volume, i.e a part of the control volume we are making the balance for. Defines the general interface. More...
class	Dumux::SubControlVolumeFaceBase< Imp, ScvfGeometryTraits >
	Base class for a sub control volume face, i.e a part of the boundary of a sub control volume we computing a flux on. More...

Enumerations
enum class	Dumux::DiscretizationMethod {   Dumux::DiscretizationMethod::none , Dumux::DiscretizationMethod::box , Dumux::DiscretizationMethod::cctpfa , Dumux::DiscretizationMethod::ccmpfa ,   Dumux::DiscretizationMethod::staggered , Dumux::DiscretizationMethod::fem }
	The available discretization methods in Dumux. More...
enum class	Dumux::RotationPolicy { Dumux::RotationPolicy::disc , Dumux::RotationPolicy::ball , Dumux::RotationPolicy::toroid }
	Rotation policies. More...

Functions
	Dumux::BaseGridGeometry< GV, Traits >::BaseGridGeometry (const GridView &gridView)
	Constructor computes the bouding box of the entire domain, for e.g. setting boundary conditions. More...
template<class Element , class FVElementGeometry , class PrimaryVariables >
auto	Dumux::evalGradients (const Element &element, const typename Element::Geometry &geometry, const typename FVElementGeometry::GridGeometry &gridGeometry, const BoxElementSolution< FVElementGeometry, PrimaryVariables > &elemSol, const typename Element::Geometry::GlobalCoordinate &globalPos, bool ignoreState=false)
	Evaluates the gradient of a given box element solution to a given global position. More...
template<class Element , class FVElementGeometry , class PrimaryVariables >
Dune::FieldVector< typename Element::Geometry::GlobalCoordinate, PrimaryVariables::dimension >	Dumux::evalGradients (const Element &element, const typename Element::Geometry &geometry, const typename FVElementGeometry::GridGeometry &gridGeometry, const CCElementSolution< FVElementGeometry, PrimaryVariables > &elemSol, const typename Element::Geometry::GlobalCoordinate &globalPos)
	Evaluates the gradient of a given CCElementSolution to a given global position. This function is only here for (compilation) compatibility reasons with the box scheme. The solution within the control volumes is constant and thus gradients are zero. One can compute gradients towards the sub-control volume faces after reconstructing the solution on the faces. More...
template<class Element , class FVElementGeometry , class PrimaryVariables >
PrimaryVariables	Dumux::evalSolution (const Element &element, const typename Element::Geometry &geometry, const typename FVElementGeometry::GridGeometry &gridGeometry, const BoxElementSolution< FVElementGeometry, PrimaryVariables > &elemSol, const typename Element::Geometry::GlobalCoordinate &globalPos, bool ignoreState=false)
	Interpolates a given box element solution at a given global position. Uses the finite element cache of the grid geometry. More...
template<class Element , class FVElementGeometry , class PrimaryVariables >
PrimaryVariables	Dumux::evalSolution (const Element &element, const typename Element::Geometry &geometry, const BoxElementSolution< FVElementGeometry, PrimaryVariables > &elemSol, const typename Element::Geometry::GlobalCoordinate &globalPos, bool ignoreState=false)
	Interpolates a given box element solution at a given global position. More...
template<class Element , class FVElementGeometry , class PrimaryVariables >
PrimaryVariables	Dumux::evalSolution (const Element &element, const typename Element::Geometry &geometry, const typename FVElementGeometry::GridGeometry &gridGeometry, const CCElementSolution< FVElementGeometry, PrimaryVariables > &elemSol, const typename Element::Geometry::GlobalCoordinate &globalPos, bool ignoreState=false)
	Interpolates a given cell-centered element solution at a given global position. More...
template<class Element , class FVElementGeometry , class PrimaryVariables >
PrimaryVariables	Dumux::evalSolution (const Element &element, const typename Element::Geometry &geometry, const CCElementSolution< FVElementGeometry, PrimaryVariables > &elemSol, const typename Element::Geometry::GlobalCoordinate &globalPos, bool ignoreState=false)
	Interpolates a given cell-centered element solution at a given global position. Overload of the above evalSolution() function without a given gridGeometry. For compatibility reasons with the box scheme. More...
template<class GridCache >
GridCache::LocalView	Dumux::localView (const GridCache &gridCache)
	Free function to get the local view of a grid cache object. More...
template<class FEBasisDomain , class FEBasisTarget , class GlueType >
auto	Dumux::makeProjectorPair (const FEBasisDomain &feBasisDomain, const FEBasisTarget &feBasisTarget, GlueType glue)
	Creates a pair of projectors between the space with basis feBasisDomain to the space with basis feBasisTarget. More...
template<class FEBasisDomain , class FEBasisTarget , class GlueType >
auto	Dumux::makeProjector (const FEBasisDomain &feBasisDomain, const FEBasisTarget &feBasisTarget, GlueType glue)
	Creates a forward projector from the space feBasisDomain to the space with basis feBasisTarget. More...

Enumeration Type Documentation

DiscretizationMethod

enum class Dumux::DiscretizationMethod

strong

The available discretization methods in Dumux.

Note

Use none if specifying a discretization method is required but the class in question is not specific to a a discretization method or the classification is non-applicable

Enumerator
none
box
cctpfa
ccmpfa
staggered
fem

RotationPolicy

enum class Dumux::RotationPolicy

strong

Rotation policies.

• disc (or annulus): rotate a segment around an axis perpendicular to the line through the segment (we choose )
• ball (or shell): rotate a segment around a point on the line through the segment
• toroid: rotate a polygon around one of its axis (we rotate about the second axis)

Enumerator
disc
ball
toroid

Function Documentation

BaseGridGeometry()

template<class GV , class Traits >

Dumux::BaseGridGeometry< GV, Traits >::BaseGridGeometry ( const GridView &  gridView )

inline

Constructor computes the bouding box of the entire domain, for e.g. setting boundary conditions.

Parameters

gridView

the grid view on which to construct the grid geometry

evalGradients() [1/2]

template<class Element , class FVElementGeometry , class PrimaryVariables >

auto Dumux::evalGradients	(	const Element &	element,
		const typename Element::Geometry &	geometry,
		const typename FVElementGeometry::GridGeometry &	gridGeometry,
		const BoxElementSolution< FVElementGeometry, PrimaryVariables > &	elemSol,
		const typename Element::Geometry::GlobalCoordinate &	globalPos,
		bool	ignoreState = false
	)

Evaluates the gradient of a given box element solution to a given global position.

Parameters

element	The element
geometry	The element geometry
gridGeometry	The finite volume grid geometry
elemSol	The primary variables at the dofs of the element
globalPos	The global position
ignoreState	If true, the state of primary variables is ignored

Returns

Dune::FieldVector with as many entries as dimension of the PrimaryVariables object (i.e. numEq). Each entry is a GlobalCoordinate object holding the priVar gradient.

evalGradients() [2/2]

template<class Element , class FVElementGeometry , class PrimaryVariables >

Dune::FieldVector< typename Element::Geometry::GlobalCoordinate, PrimaryVariables::dimension > Dumux::evalGradients	(	const Element &	element,
		const typename Element::Geometry &	geometry,
		const typename FVElementGeometry::GridGeometry &	gridGeometry,
		const CCElementSolution< FVElementGeometry, PrimaryVariables > &	elemSol,
		const typename Element::Geometry::GlobalCoordinate &	globalPos
	)

Evaluates the gradient of a given CCElementSolution to a given global position. This function is only here for (compilation) compatibility reasons with the box scheme. The solution within the control volumes is constant and thus gradients are zero. One can compute gradients towards the sub-control volume faces after reconstructing the solution on the faces.

Parameters

element	The element
geometry	The element geometry
gridGeometry	The finite volume grid geometry
elemSol	The primary variables at the dofs of the element
globalPos	The global position

Exceptions

Dune::NotImplemented

Returns

Dune::FieldVector with as many entries as dimension of the PrimaryVariables object (i.e. numEq). Each entry is a GlobalCoordinate object holding the priVar gradient.

evalSolution() [1/4]

template<class Element , class FVElementGeometry , class PrimaryVariables >

PrimaryVariables Dumux::evalSolution	(	const Element &	element,
		const typename Element::Geometry &	geometry,
		const BoxElementSolution< FVElementGeometry, PrimaryVariables > &	elemSol,
		const typename Element::Geometry::GlobalCoordinate &	globalPos,
		bool	ignoreState = false
	)

Interpolates a given box element solution at a given global position.

Overload of the above evalSolution() function without a given gridGeometry. The local basis is computed on the fly.

Returns

the interpolated primary variables

Parameters

element	The element
geometry	The element geometry
elemSol	The primary variables at the dofs of the element
globalPos	The global position
ignoreState	If true, the state of primary variables is ignored

The box scheme always uses linear Ansatz functions

evalSolution() [2/4]

template<class Element , class FVElementGeometry , class PrimaryVariables >

PrimaryVariables Dumux::evalSolution	(	const Element &	element,
		const typename Element::Geometry &	geometry,
		const CCElementSolution< FVElementGeometry, PrimaryVariables > &	elemSol,
		const typename Element::Geometry::GlobalCoordinate &	globalPos,
		bool	ignoreState = false
	)

Interpolates a given cell-centered element solution at a given global position. Overload of the above evalSolution() function without a given gridGeometry. For compatibility reasons with the box scheme.

Returns

the primary variables (constant over the element)

Parameters

element	The element
geometry	The element geometry
elemSol	The primary variables at the dofs of the element
globalPos	The global position
ignoreState	If true, the state of primary variables is ignored

evalSolution() [3/4]

template<class Element , class FVElementGeometry , class PrimaryVariables >

PrimaryVariables Dumux::evalSolution	(	const Element &	element,
		const typename Element::Geometry &	geometry,
		const typename FVElementGeometry::GridGeometry &	gridGeometry,
		const BoxElementSolution< FVElementGeometry, PrimaryVariables > &	elemSol,
		const typename Element::Geometry::GlobalCoordinate &	globalPos,
		bool	ignoreState = false
	)

Interpolates a given box element solution at a given global position. Uses the finite element cache of the grid geometry.

Returns

the interpolated primary variables

Parameters

element	The element
geometry	The element geometry
gridGeometry	The finite volume grid geometry
elemSol	The primary variables at the dofs of the element
globalPos	The global position
ignoreState	If true, the state of primary variables is ignored

evalSolution() [4/4]

template<class Element , class FVElementGeometry , class PrimaryVariables >

PrimaryVariables Dumux::evalSolution	(	const Element &	element,
		const typename Element::Geometry &	geometry,
		const typename FVElementGeometry::GridGeometry &	gridGeometry,
		const CCElementSolution< FVElementGeometry, PrimaryVariables > &	elemSol,
		const typename Element::Geometry::GlobalCoordinate &	globalPos,
		bool	ignoreState = false
	)

Interpolates a given cell-centered element solution at a given global position.

Returns

the primary variables (constant over the element)

Parameters

element	The element
geometry	The element geometry
gridGeometry	The finite volume grid geometry
elemSol	The primary variables at the dofs of the element
globalPos	The global position
ignoreState	If true, the state of primary variables is ignored

localView()

template<class GridCache >

GridCache::LocalView Dumux::localView ( const GridCache &  gridCache )

inline

Free function to get the local view of a grid cache object.

Note

A local object is only functional after calling its bind/bindElement method.

Template Parameters

GridCache

the grid caching type (such as GridGeometry)

Parameters

gridCache

the grid caching object we want to localView from

makeProjector()

template<class FEBasisDomain , class FEBasisTarget , class GlueType >

auto Dumux::makeProjector	(	const FEBasisDomain &	feBasisDomain,
		const FEBasisTarget &	feBasisTarget,
		GlueType	glue
	)

Creates a forward projector from the space feBasisDomain to the space with basis feBasisTarget.

Parameters

feBasisDomain	The domain finite element space basis
feBasisTarget	The target finite element space basis
glue	The glue object containing the intersections between the grids.

Returns

The forward projector from the space with basis feBasisDomain to the space with basis feBasisTarget.

makeProjectorPair()

template<class FEBasisDomain , class FEBasisTarget , class GlueType >

auto Dumux::makeProjectorPair	(	const FEBasisDomain &	feBasisDomain,
		const FEBasisTarget &	feBasisTarget,
		GlueType	glue
	)

Creates a pair of projectors between the space with basis feBasisDomain to the space with basis feBasisTarget.

Parameters

feBasisDomain	The domain finite element space basis
feBasisTarget	The target finite element space basis
glue	The glue object containing the intersections between the grids.

Returns

An std::pair of projectors where the first is the forward projector from the space with basis feBasisDomain to the space with basis feBasisTarget and the second does the backward projection.