Definition of the spatial parameters for the heterogeneous problem which uses the non-isothermal or isothermal CO2 fully implicit model. More...

#include <test/porousmediumflow/co2/implicit/spatialparams.hh>

Inheritance diagram for Dumux::HeterogeneousSpatialParams< GridGeometry, Scalar >:

Description

template<class GridGeometry, class Scalar>
class Dumux::HeterogeneousSpatialParams< GridGeometry, Scalar >

Definition of the spatial parameters for the heterogeneous problem which uses the non-isothermal or isothermal CO2 fully implicit model.

Public Types
using	MaterialLaw = EffToAbsLaw< EffectiveLaw >

using	MaterialLawParams = typename MaterialLaw::Params

using	PermeabilityType = Scalar

Public Member Functions
	HeterogeneousSpatialParams (std::shared_ptr< const GridGeometry > gridGeometry, std::shared_ptr< const GridData< Grid > > gridData)

void	getParamsFromGrid ()
	Reads layer information from the grid. More...

template<class ElementSolution >
PermeabilityType	permeability (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const
	Function for defining the (intrinsic) permeability \([m^2]\). More...

PermeabilityType	permeability (std::size_t eIdx) const
	Function for defining the (intrinsic) permeability \([m^2]\). More...

template<class SolidSystem , class ElementSolution >
Scalar	inertVolumeFraction (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol, int compIdx) const
	Returns the volume fraction of the inert component with index compIdx \([-]\). More...

Scalar	inertVolumeFraction (std::size_t eIdx) const
	Returns the porosity \([-]\). More...

const MaterialLawParams &	materialLawParamsAtPos (const GlobalPosition &globalPos) const
	Function for defining the parameters needed by constitutive relationships (kr-sw, pc-sw, etc.). More...

template<class FluidSystem >
int	wettingPhaseAtPos (const GlobalPosition &globalPos) const
	Function for defining which phase is to be considered as the wetting phase. More...

decltype(auto)	materialLawParams (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const
	Function for defining the parameters needed by constitutive relationships (kr-sw, pc-sw, etc.). More...

int	wettingPhase (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const
	Function for defining which phase is to be considered as the wetting phase. More...

int	wettingPhaseAtPos (const GlobalPosition &globalPos) const
	Function for defining which phase is to be considered as the wetting phase. More...

const GlobalPosition &	gravity (const GlobalPosition &pos) const
	Returns the acceleration due to gravity \(\mathrm{[m/s^2]}\). More...

Scalar	harmonicMean (const Scalar T1, const Scalar T2, const GlobalPosition &normal) const
	Harmonic average of a discontinuous scalar field at discontinuity interface (for compatibility reasons with the function below) More...

DimWorldMatrix	harmonicMean (const DimWorldMatrix &T1, const DimWorldMatrix &T2, const GlobalPosition &normal) const
	Harmonic average of a discontinuous tensorial field at discontinuity interface. More...

decltype(auto)	permeability (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const
	Function for defining the (intrinsic) permeability \([m^2]\). More...

Scalar	porosity (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const
	Function for defining the porosity. That is possibly solution dependent. More...

Scalar	inertVolumeFraction (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol, int compIdx) const
	Function for defining the solid volume fraction. That is possibly solution dependent. More...

Scalar	beaversJosephCoeffAtPos (const GlobalPosition &globalPos) const
	Function for defining the Beavers-Joseph coefficient for multidomain problems \(\mathrm{[-]}\). More...

Scalar	forchCoeff (const SubControlVolumeFace &scvf) const
	Apply the Forchheimer coefficient for inertial forces calculation. More...

const GridGeometry &	fvGridGeometry () const
	The finite volume grid geometry. More...

const GridGeometry &	gridGeometry () const
	The finite volume grid geometry. More...

Static Public Member Functions
static constexpr bool	evaluatePermeabilityAtScvfIP ()
	If the permeability should be evaluated directly at the scvf integration point (for convergence tests with analytical and continuous perm functions) or is evaluated at the scvs (for permeability fields with discontinuities) -> default. More...

Protected Member Functions
HeterogeneousSpatialParams< GridGeometry, Scalar > &	asImp_ ()

const HeterogeneousSpatialParams< GridGeometry, Scalar > &	asImp_ () const

Member Typedef Documentation

◆ MaterialLaw

template<class GridGeometry , class Scalar >

using Dumux::HeterogeneousSpatialParams< GridGeometry, Scalar >::MaterialLaw = EffToAbsLaw<EffectiveLaw>

◆ MaterialLawParams

template<class GridGeometry , class Scalar >

using Dumux::HeterogeneousSpatialParams< GridGeometry, Scalar >::MaterialLawParams = typename MaterialLaw::Params

◆ PermeabilityType

template<class GridGeometry , class Scalar >

using Dumux::HeterogeneousSpatialParams< GridGeometry, Scalar >::PermeabilityType = Scalar

Constructor & Destructor Documentation

◆ HeterogeneousSpatialParams()

template<class GridGeometry , class Scalar >

Dumux::HeterogeneousSpatialParams< GridGeometry, Scalar >::HeterogeneousSpatialParams	(	std::shared_ptr< const GridGeometry >	gridGeometry,
		std::shared_ptr< const GridData< Grid > >	gridData
	)

inline

Member Function Documentation

◆ asImp_() [1/2]

HeterogeneousSpatialParams< GridGeometry, Scalar > & Dumux::FVSpatialParamsOneP< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::asImp_ ( )

inlineprotectedinherited

◆ asImp_() [2/2]

const HeterogeneousSpatialParams< GridGeometry, Scalar > & Dumux::FVSpatialParamsOneP< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::asImp_ ( ) const

inlineprotectedinherited

◆ beaversJosephCoeffAtPos()

Scalar Dumux::FVSpatialParamsOneP< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::beaversJosephCoeffAtPos ( const GlobalPosition & globalPos ) const

inlineinherited

Function for defining the Beavers-Joseph coefficient for multidomain problems \(\mathrm{[-]}\).

Returns: Beavers-Joseph coefficient \(\mathrm{[-]}\)

Parameters

globalPos The global position

◆ evaluatePermeabilityAtScvfIP()

static constexpr bool Dumux::FVSpatialParamsOneP< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::evaluatePermeabilityAtScvfIP ( )

inlinestaticconstexprinherited

If the permeability should be evaluated directly at the scvf integration point (for convergence tests with analytical and continuous perm functions) or is evaluated at the scvs (for permeability fields with discontinuities) -> default.

◆ forchCoeff()

Scalar Dumux::FVSpatialParamsOneP< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::forchCoeff ( const SubControlVolumeFace & scvf ) const

inlineinherited

Apply the Forchheimer coefficient for inertial forces calculation.

Parameters

scvf	The sub-control volume face where the intrinsic velocity ought to be calculated.

◆ fvGridGeometry()

const GridGeometry & Dumux::FVSpatialParamsOneP< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::fvGridGeometry ( ) const

inlineinherited

The finite volume grid geometry.

◆ getParamsFromGrid()

template<class GridGeometry , class Scalar >

void Dumux::HeterogeneousSpatialParams< GridGeometry, Scalar >::getParamsFromGrid ( )

inline

Reads layer information from the grid.

◆ gravity()

const GlobalPosition & Dumux::FVSpatialParamsOneP< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::gravity ( const GlobalPosition & pos ) const

inlineinherited

Returns the acceleration due to gravity \(\mathrm{[m/s^2]}\).

The default behaviour is a constant gravity vector; if the Problem.EnableGravity parameter is true, \(\boldsymbol{g} = ( 0,\dots,\ -9.81)^T \), else \(\boldsymbol{g} = ( 0,\dots, 0)^T \).

Parameters

pos	the spatial position at which to evaulate the gravity vector

◆ gridGeometry()

const GridGeometry & Dumux::FVSpatialParamsOneP< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::gridGeometry ( ) const

inlineinherited

The finite volume grid geometry.

◆ harmonicMean() [1/2]

DimWorldMatrix Dumux::FVSpatialParamsOneP< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::harmonicMean	(	const DimWorldMatrix &	T1,
		const DimWorldMatrix &	T2,
		const GlobalPosition &	normal
	)		const

inlineinherited

Harmonic average of a discontinuous tensorial field at discontinuity interface.

Note: We do a harmonic average of the part normal to the interface (alpha*I) and an arithmetic average of the tangential part (T - alpha*I).

Returns: the averaged tensor

Parameters

T1	first tensor
T2	second tensor
normal	The unit normal vector of the interface

◆ harmonicMean() [2/2]

Scalar Dumux::FVSpatialParamsOneP< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::harmonicMean	(	const Scalar	T1,
		const Scalar	T2,
		const GlobalPosition &	normal
	)		const

inlineinherited

Harmonic average of a discontinuous scalar field at discontinuity interface (for compatibility reasons with the function below)

Returns: the averaged scalar

Parameters

T1	first scalar parameter
T2	second scalar parameter
normal	The unit normal vector of the interface

◆ inertVolumeFraction() [1/3]

Scalar Dumux::FVSpatialParamsOneP< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::inertVolumeFraction	(	const Element &	element,
		const SubControlVolume &	scv,
		const ElementSolution &	elemSol,
		int	compIdx
	)		const

inlineinherited

Function for defining the solid volume fraction. That is possibly solution dependent.

Parameters

element	The current element
scv	The sub-control volume inside the element.
elemSol	The solution at the dofs connected to the element.
compIdx	The solid component index

Returns: the volume fraction of the inert solid component with index compIdx

Note: this overload is enable if there is only one inert solid component and the user didn't choose to implement a inertVolumeFractionAtPos overload. It then forwards to the simpler porosity interface. With more than one solid components or active solid components (i.e. dissolution) please overload the more general inertVolumeFraction/inertVolumeFractionAtPos interface.

◆ inertVolumeFraction() [2/3]

template<class GridGeometry , class Scalar >

template<class SolidSystem , class ElementSolution >

Scalar Dumux::HeterogeneousSpatialParams< GridGeometry, Scalar >::inertVolumeFraction	(	const Element &	element,
		const SubControlVolume &	scv,
		const ElementSolution &	elemSol,
		int	compIdx
	)		const

inline

Returns the volume fraction of the inert component with index compIdx \([-]\).

Parameters

element	The current element
scv	The sub-control volume inside the element.
elemSol	The element solution
compIdx	The solid component index

Returns: The solid volume fraction

◆ inertVolumeFraction() [3/3]

template<class GridGeometry , class Scalar >

Scalar Dumux::HeterogeneousSpatialParams< GridGeometry, Scalar >::inertVolumeFraction ( std::size_t eIdx ) const

inline

Returns the porosity \([-]\).

Parameters

eIdx	The element index

◆ materialLawParams()

decltype(auto) Dumux::FVSpatialParams< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::materialLawParams	(	const Element &	element,
		const SubControlVolume &	scv,
		const ElementSolution &	elemSol
	)		const

inlineinherited

Function for defining the parameters needed by constitutive relationships (kr-sw, pc-sw, etc.).

Parameters

element	The current element
scv	The sub-control volume inside the element.
elemSol	The solution at the dofs connected to the element.

Returns: the material parameters object

◆ materialLawParamsAtPos()

template<class GridGeometry , class Scalar >

const MaterialLawParams & Dumux::HeterogeneousSpatialParams< GridGeometry, Scalar >::materialLawParamsAtPos ( const GlobalPosition & globalPos ) const

inline

Function for defining the parameters needed by constitutive relationships (kr-sw, pc-sw, etc.).

Parameters

globalPos The position of the center of the element

Returns: The material parameters object

◆ permeability() [1/3]

decltype(auto) Dumux::FVSpatialParamsOneP< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::permeability	(	const Element &	element,
		const SubControlVolume &	scv,
		const ElementSolution &	elemSol
	)		const

inlineinherited

Function for defining the (intrinsic) permeability \([m^2]\).

Note: It is possibly solution dependent.

Parameters

element	The current element
scv	The sub-control volume inside the element.
elemSol	The solution at the dofs connected to the element.

Returns: permeability

◆ permeability() [2/3]

template<class GridGeometry , class Scalar >

template<class ElementSolution >

PermeabilityType Dumux::HeterogeneousSpatialParams< GridGeometry, Scalar >::permeability	(	const Element &	element,
		const SubControlVolume &	scv,
		const ElementSolution &	elemSol
	)		const

inline

Function for defining the (intrinsic) permeability \([m^2]\).

Note: It is possibly solution dependent.

Parameters

element	The current element
scv	The sub-control volume inside the element.
elemSol	The solution at the dofs connected to the element.

Returns: The instrinsic permeability

◆ permeability() [3/3]

template<class GridGeometry , class Scalar >

PermeabilityType Dumux::HeterogeneousSpatialParams< GridGeometry, Scalar >::permeability ( std::size_t eIdx ) const

inline

Function for defining the (intrinsic) permeability \([m^2]\).

Note: It is possibly solution dependent.

Parameters

eIdx	The element index

◆ porosity()

Scalar Dumux::FVSpatialParamsOneP< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::porosity	(	const Element &	element,
		const SubControlVolume &	scv,
		const ElementSolution &	elemSol
	)		const

inlineinherited

Function for defining the porosity. That is possibly solution dependent.

Note: this can only be used for solids with one inert component (see inertVolumeFraction for the more general interface)

Parameters

element	The current element
scv	The sub-control volume inside the element.
elemSol	The solution at the dofs connected to the element.

Returns: the porosity

◆ wettingPhase()

int Dumux::FVSpatialParams< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::wettingPhase	(	const Element &	element,
		const SubControlVolume &	scv,
		const ElementSolution &	elemSol
	)		const

inlineinherited

Function for defining which phase is to be considered as the wetting phase.

Parameters

element	The current element
scv	The sub-control volume inside the element.
elemSol	The solution at the dofs connected to the element.

Returns: the wetting phase index

◆ wettingPhaseAtPos() [1/2]

int Dumux::FVSpatialParams< GridGeometry, Scalar, HeterogeneousSpatialParams< GridGeometry, Scalar > >::wettingPhaseAtPos ( const GlobalPosition & globalPos ) const

inlineinherited

Function for defining which phase is to be considered as the wetting phase.

Returns: the wetting phase index

Parameters

globalPos The global position

◆ wettingPhaseAtPos() [2/2]

template<class GridGeometry , class Scalar >

template<class FluidSystem >

int Dumux::HeterogeneousSpatialParams< GridGeometry, Scalar >::wettingPhaseAtPos ( const GlobalPosition & globalPos ) const

inline

Function for defining which phase is to be considered as the wetting phase.

Parameters

globalPos The position of the center of the element

Returns: The wetting phase index

The documentation for this class was generated from the following file:

porousmediumflow/co2/implicit/spatialparams.hh

Description

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Member Typedef Documentation

◆ MaterialLaw

◆ MaterialLawParams

◆ PermeabilityType

Constructor & Destructor Documentation

◆ HeterogeneousSpatialParams()

Member Function Documentation

◆ asImp_() [1/2]

◆ asImp_() [2/2]

◆ beaversJosephCoeffAtPos()

◆ evaluatePermeabilityAtScvfIP()

◆ forchCoeff()

◆ fvGridGeometry()

◆ getParamsFromGrid()

◆ gravity()

◆ gridGeometry()

◆ harmonicMean() [1/2]

◆ harmonicMean() [2/2]

◆ inertVolumeFraction() [1/3]

◆ inertVolumeFraction() [2/3]

◆ inertVolumeFraction() [3/3]

◆ materialLawParams()

◆ materialLawParamsAtPos()

◆ permeability() [1/3]

◆ permeability() [2/3]

◆ permeability() [3/3]

◆ porosity()

◆ wettingPhase()

◆ wettingPhaseAtPos() [1/2]

◆ wettingPhaseAtPos() [2/2]