Calculates the element-wise residual for cell-centered discretization schemes in models where coupling occurs across the element facets. We only overwrite the function for the computation of a flux across a single sub-control volume face, as we need to additionally check if a boundary face couples to a facet element. More...

#include <dumux/multidomain/facet/cellcentered/localresidual.hh>

Inheritance diagram for Dumux::CCFacetCouplingLocalResidual< TypeTag >:

Description

template<class TypeTag>
class Dumux::CCFacetCouplingLocalResidual< TypeTag >

Calculates the element-wise residual for cell-centered discretization schemes in models where coupling occurs across the element facets. We only overwrite the function for the computation of a flux across a single sub-control volume face, as we need to additionally check if a boundary face couples to a facet element.

Public Types
using	ElementResidualVector = typename ParentType::ElementResidualVector
	export the type used for element residuals More...

Public Member Functions
template<class Problem >
NumEqVector	evalFlux (const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFluxVariablesCache &elemFluxVarsCache, const SubControlVolumeFace &scvf) const

void	evalFlux (ElementResidualVector &residual, const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementBoundaryTypes &elemBcTypes, const ElementFluxVariablesCache &elemFluxVarsCache, const SubControlVolumeFace &scvf) const
	evaluate the flux residual for a sub control volume face More...

NumEqVector	evalFlux (const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFluxVariablesCache &elemFluxVarsCache, const SubControlVolumeFace &scvf) const
	evaluate the flux residual for a sub control volume face More...

Discretization specific interface
Note The following method are the discretization specific wrapper methods
void	evalFlux (ElementResidualVector &residual, const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementBoundaryTypes &elemBcTypes, const ElementFluxVariablesCache &elemFluxVarsCache, const SubControlVolumeFace &scvf) const
	Compute the source local residual, i.e. the deviation of the source term from zero and add to the residual. More...

NumEqVector	evalFlux (const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFluxVariablesCache &elemFluxVarsCache, const SubControlVolumeFace &scvf) const
	Compute the flux local residual, i.e. the deviation of the flux term from zero. More...

void	evalStorage (ElementResidualVector &residual, const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &prevElemVolVars, const ElementVolumeVariables &curElemVolVars, const SubControlVolume &scv) const
	Compute the storage local residual, i.e. the deviation of the storage term from zero for instationary problems. More...

void	evalSource (ElementResidualVector &residual, const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &curElemVolVars, const SubControlVolume &scv) const
	Compute the source local residual, i.e. the deviation of the source term from zero. More...

User interface
Note The following methods are usually expensive to evaluate They are useful for outputting / postprocessing residual information.
ElementResidualVector	evalStorage (const Problem &problem, const Element &element, const GridGeometry &gridGeometry, const GridVariables &gridVariables, const SolutionVector &sol) const
	Compute the storage term for the current solution. More...

Main interface
Note Methods used by the assembler to compute derivatives and residual
ElementResidualVector	evalStorage (const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &prevElemVolVars, const ElementVolumeVariables &curElemVolVars) const
	Compute the storage local residual, i.e. the deviation of the storage term from zero for instationary problems. More...

ElementResidualVector	evalFluxAndSource (const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFluxVariablesCache &elemFluxVarsCache, const ElementBoundaryTypes &bcTypes) const

Model specific interface
Note The following method are the model specific implementations of the local residual
NumEqVector	computeStorage (const Problem &problem, const SubControlVolume &scv, const VolumeVariables &volVars) const
	Calculate the source term of the equation. More...

NumEqVector	computeSource (const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const SubControlVolume &scv) const
	Calculate the source term of the equation. More...

NumEqVector	computeFlux (const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const SubControlVolumeFace &scvf, const ElementFluxVariablesCache &elemFluxVarsCache) const
	Calculate the source term of the equation. More...

Interfaces for analytic Jacobian computation
template<class PartialDerivativeMatrix >
void	addStorageDerivatives (PartialDerivativeMatrix &partialDerivatives, const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const VolumeVariables &curVolVars, const SubControlVolume &scv) const
	Compute the derivative of the storage residual. More...

template<class PartialDerivativeMatrix >
void	addSourceDerivatives (PartialDerivativeMatrix &partialDerivatives, const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const VolumeVariables &curVolVars, const SubControlVolume &scv) const
	Compute the derivative of the source residual. More...

template<class PartialDerivativeMatrices , class T = TypeTag>
std::enable_if_t< GetPropType< T, Properties::GridGeometry >::discMethod !=DiscretizationMethods::box, void >	addFluxDerivatives (PartialDerivativeMatrices &derivativeMatrices, const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &curElemVolVars, const ElementFluxVariablesCache &elemFluxVarsCache, const SubControlVolumeFace &scvf) const
	Compute the derivative of the flux residual. More...

template<class JacobianMatrix , class T = TypeTag>
std::enable_if_t< GetPropType< T, Properties::GridGeometry >::discMethod==DiscretizationMethods::box, void >	addFluxDerivatives (JacobianMatrix &A, const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &curElemVolVars, const ElementFluxVariablesCache &elemFluxVarsCache, const SubControlVolumeFace &scvf) const
	Compute the derivative of the flux residual for the box method. More...

template<class PartialDerivativeMatrices >
void	addCCDirichletFluxDerivatives (PartialDerivativeMatrices &derivativeMatrices, const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &curElemVolVars, const ElementFluxVariablesCache &elemFluxVarsCache, const SubControlVolumeFace &scvf) const
	Compute the derivative of the Dirichlet flux residual for cell-centered schemes. More...

template<class PartialDerivativeMatrices >
void	addRobinFluxDerivatives (PartialDerivativeMatrices &derivativeMatrices, const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &curElemVolVars, const ElementFluxVariablesCache &elemFluxVarsCache, const SubControlVolumeFace &scvf) const
	Compute the derivative of Robin type boundary conditions ("solution dependent Neumann") More...

Interfaces accessed by local residual implementations
const Problem &	problem () const
	the problem More...

const TimeLoop &	timeLoop () const

bool	isStationary () const
	returns true if the residual is stationary More...

Implementation &	asImp ()

const Implementation &	asImp () const

Member Typedef Documentation

◆ ElementResidualVector

template<class TypeTag >

using Dumux::CCFacetCouplingLocalResidual< TypeTag >::ElementResidualVector = typename ParentType::ElementResidualVector

export the type used for element residuals

Member Function Documentation

◆ addCCDirichletFluxDerivatives()

template<class TypeTag >

template<class PartialDerivativeMatrices >

void Dumux::FVLocalResidual< TypeTag >::addCCDirichletFluxDerivatives	(	PartialDerivativeMatrices &	derivativeMatrices,
		const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	curElemVolVars,
		const ElementFluxVariablesCache &	elemFluxVarsCache,
		const SubControlVolumeFace &	scvf
	)		const

inlineinherited

Compute the derivative of the Dirichlet flux residual for cell-centered schemes.

◆ addFluxDerivatives() [1/2]

template<class TypeTag >

template<class JacobianMatrix , class T = TypeTag>

std::enable_if_t< GetPropType< T, Properties::GridGeometry >::discMethod==DiscretizationMethods::box, void > Dumux::FVLocalResidual< TypeTag >::addFluxDerivatives	(	JacobianMatrix &	A,
		const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	curElemVolVars,
		const ElementFluxVariablesCache &	elemFluxVarsCache,
		const SubControlVolumeFace &	scvf
	)		const

inlineinherited

Compute the derivative of the flux residual for the box method.

◆ addFluxDerivatives() [2/2]

template<class TypeTag >

template<class PartialDerivativeMatrices , class T = TypeTag>

std::enable_if_t< GetPropType< T, Properties::GridGeometry >::discMethod !=DiscretizationMethods::box, void > Dumux::FVLocalResidual< TypeTag >::addFluxDerivatives	(	PartialDerivativeMatrices &	derivativeMatrices,
		const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	curElemVolVars,
		const ElementFluxVariablesCache &	elemFluxVarsCache,
		const SubControlVolumeFace &	scvf
	)		const

inlineinherited

Compute the derivative of the flux residual.

◆ addRobinFluxDerivatives()

template<class TypeTag >

template<class PartialDerivativeMatrices >

void Dumux::FVLocalResidual< TypeTag >::addRobinFluxDerivatives	(	PartialDerivativeMatrices &	derivativeMatrices,
		const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	curElemVolVars,
		const ElementFluxVariablesCache &	elemFluxVarsCache,
		const SubControlVolumeFace &	scvf
	)		const

inlineinherited

Compute the derivative of Robin type boundary conditions ("solution dependent Neumann")

◆ addSourceDerivatives()

template<class TypeTag >

template<class PartialDerivativeMatrix >

void Dumux::FVLocalResidual< TypeTag >::addSourceDerivatives	(	PartialDerivativeMatrix &	partialDerivatives,
		const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const VolumeVariables &	curVolVars,
		const SubControlVolume &	scv
	)		const

inlineinherited

Compute the derivative of the source residual.

◆ addStorageDerivatives()

template<class TypeTag >

template<class PartialDerivativeMatrix >

void Dumux::FVLocalResidual< TypeTag >::addStorageDerivatives	(	PartialDerivativeMatrix &	partialDerivatives,
		const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const VolumeVariables &	curVolVars,
		const SubControlVolume &	scv
	)		const

inlineinherited

Compute the derivative of the storage residual.

◆ asImp() [1/2]

template<class TypeTag >

Implementation & Dumux::FVLocalResidual< TypeTag >::asImp ( )

inlineprotectedinherited

◆ asImp() [2/2]

template<class TypeTag >

const Implementation & Dumux::FVLocalResidual< TypeTag >::asImp ( ) const

inlineprotectedinherited

◆ computeFlux()

template<class TypeTag >

NumEqVector Dumux::FVLocalResidual< TypeTag >::computeFlux	(	const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const SubControlVolumeFace &	scvf,
		const ElementFluxVariablesCache &	elemFluxVarsCache
	)		const

inlineinherited

Calculate the source term of the equation.

Parameters

problem	The problem to solve
element	The DUNE Codim<0> entity for which the residual ought to be calculated
fvGeometry	The finite-volume geometry of the element
elemVolVars	The volume variables associated with the element stencil
scvf	The sub-control volume over which we integrate the flux
elemFluxVarsCache	the flux variable caches for the element's flux stencils

Note: has to be implemented by the model specific residual class

◆ computeSource()

template<class TypeTag >

NumEqVector Dumux::FVLocalResidual< TypeTag >::computeSource	(	const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const SubControlVolume &	scv
	)		const

inlineinherited

Calculate the source term of the equation.

Parameters

problem	The problem to solve
element	The DUNE Codim<0> entity for which the residual ought to be calculated
fvGeometry	The finite-volume geometry of the element
elemVolVars	The volume variables associated with the element stencil
scv	The sub-control volume over which we integrate the source term

Note: This is the default implementation for all models as sources are computed in the user interface of the problem

◆ computeStorage()

template<class TypeTag >

NumEqVector Dumux::FVLocalResidual< TypeTag >::computeStorage	(	const Problem &	problem,
		const SubControlVolume &	scv,
		const VolumeVariables &	volVars
	)		const

inlineinherited

Calculate the source term of the equation.

Parameters

problem	The problem to solve
scv	The sub-control volume over which we integrate the storage term
volVars	The volume variables associated with the scv

Note: has to be implemented by the model specific residual class

◆ evalFlux() [1/5]

template<class TypeTag >

NumEqVector Dumux::FVLocalResidual< TypeTag >::evalFlux	(	const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const ElementFluxVariablesCache &	elemFluxVarsCache,
		const SubControlVolumeFace &	scvf
	)		const

inlineinherited

Compute the flux local residual, i.e. the deviation of the flux term from zero.

Parameters

problem	The problem to solve
element	The DUNE Codim<0> entity for which the residual ought to be calculated
fvGeometry	The finite-volume geometry of the element
elemVolVars	The volume averaged variables for all sub-control volumes of the element at the current time level
elemFluxVarsCache	The flux variable caches for the element stencil
scvf	The sub control volume face the flux term is integrated over

◆ evalFlux() [2/5]

template<class TypeTag >

NumEqVector Dumux::FVLocalResidual< TypeTag >::evalFlux	(	const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const ElementFluxVariablesCache &	elemFluxVarsCache,
		const SubControlVolumeFace &	scvf
	)		const

inline

evaluate the flux residual for a sub control volume face

◆ evalFlux() [3/5]

template<class TypeTag >

template<class Problem >

NumEqVector Dumux::CCFacetCouplingLocalResidual< TypeTag >::evalFlux	(	const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const ElementFluxVariablesCache &	elemFluxVarsCache,
		const SubControlVolumeFace &	scvf
	)		const

inline

◆ evalFlux() [4/5]

template<class TypeTag >

void Dumux::FVLocalResidual< TypeTag >::evalFlux	(	ElementResidualVector &	residual,
		const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const ElementBoundaryTypes &	elemBcTypes,
		const ElementFluxVariablesCache &	elemFluxVarsCache,
		const SubControlVolumeFace &	scvf
	)		const

inlineinherited

Compute the source local residual, i.e. the deviation of the source term from zero and add to the residual.

Note: This is implemented for different discretization schemes

Parameters

residual	The residual vector to fill
problem	The problem to solve
element	The DUNE Codim<0> entity for which the residual ought to be calculated
fvGeometry	The finite-volume geometry of the element
elemVolVars	The volume averaged variables for all sub-control volumes of the element at the current time level
elemBcTypes	the boundary types for the boundary entities of an elements
elemFluxVarsCache	The flux variable caches for the element stencil
scvf	The sub control volume face the flux term is integrated over

◆ evalFlux() [5/5]

template<class TypeTag >

void Dumux::FVLocalResidual< TypeTag >::evalFlux	(	ElementResidualVector &	residual,
		const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const ElementBoundaryTypes &	elemBcTypes,
		const ElementFluxVariablesCache &	elemFluxVarsCache,
		const SubControlVolumeFace &	scvf
	)		const

inline

evaluate the flux residual for a sub control volume face

◆ evalFluxAndSource()

template<class TypeTag >

ElementResidualVector Dumux::FVLocalResidual< TypeTag >::evalFluxAndSource	(	const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const ElementFluxVariablesCache &	elemFluxVarsCache,
		const ElementBoundaryTypes &	bcTypes
	)		const

inlineinherited

◆ evalSource()

template<class TypeTag >

void Dumux::FVLocalResidual< TypeTag >::evalSource	(	ElementResidualVector &	residual,
		const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	curElemVolVars,
		const SubControlVolume &	scv
	)		const

inlineinherited

Compute the source local residual, i.e. the deviation of the source term from zero.

Parameters

residual	The residual vector to fill
problem	The problem to solve
element	The DUNE Codim<0> entity for which the residual ought to be calculated
fvGeometry	The finite-volume geometry of the element
curElemVolVars	The volume averaged variables for all sub-control volumes of the element at the current time level
scv	The sub control volume the source term is integrated over

Compute source with the model specific storage residual

subtract source from local rate (sign convention in user interface)

◆ evalStorage() [1/3]

template<class TypeTag >

ElementResidualVector Dumux::FVLocalResidual< TypeTag >::evalStorage	(	const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	prevElemVolVars,
		const ElementVolumeVariables &	curElemVolVars
	)		const

inlineinherited

Compute the storage local residual, i.e. the deviation of the storage term from zero for instationary problems.

Parameters

element	The DUNE Codim<0> entity for which the residual ought to be calculated
fvGeometry	The finite-volume geometry of the element
prevElemVolVars	The volume averaged variables for all sub-control volumes of the element at the previous time level
curElemVolVars	The volume averaged variables for all sub-control volumes of the element at the current time level

◆ evalStorage() [2/3]

template<class TypeTag >

ElementResidualVector Dumux::FVLocalResidual< TypeTag >::evalStorage	(	const Problem &	problem,
		const Element &	element,
		const GridGeometry &	gridGeometry,
		const GridVariables &	gridVariables,
		const SolutionVector &	sol
	)		const

inlineinherited

Compute the storage term for the current solution.

This can be used to figure out how much of each conservation quantity is inside the element.

Parameters

problem	The problem to solve
element	The DUNE Codim<0> entity for which the storage term ought to be calculated
gridGeometry	The finite-volume grid geometry
gridVariables	The grid variables (volume and flux variables)
sol	The solution vector

◆ evalStorage() [3/3]

template<class TypeTag >

void Dumux::FVLocalResidual< TypeTag >::evalStorage	(	ElementResidualVector &	residual,
		const Problem &	problem,
		const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	prevElemVolVars,
		const ElementVolumeVariables &	curElemVolVars,
		const SubControlVolume &	scv
	)		const

inlineinherited

Compute the storage local residual, i.e. the deviation of the storage term from zero for instationary problems.

Parameters

residual	The residual vector to fill
problem	The problem to solve
element	The DUNE Codim<0> entity for which the residual ought to be calculated
fvGeometry	The finite-volume geometry of the element
prevElemVolVars	The volume averaged variables for all sub-control volumes of the element at the previous time level
curElemVolVars	The volume averaged variables for all sub-control volumes of the element at the current time level
scv	The sub control volume the storage term is integrated over

Compute storage with the model specific storage residual

◆ isStationary()

template<class TypeTag >

bool Dumux::FVLocalResidual< TypeTag >::isStationary ( ) const

inlineinherited

returns true if the residual is stationary

◆ problem()

template<class TypeTag >

const Problem & Dumux::FVLocalResidual< TypeTag >::problem ( ) const

inlineinherited

the problem

◆ timeLoop()

template<class TypeTag >

const TimeLoop & Dumux::FVLocalResidual< TypeTag >::timeLoop ( ) const

inlineinherited

the timeloop for instationary problems calling this for stationary leads to undefined behaviour

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

multidomain/facet/cellcentered/localresidual.hh

Description

Public Types

Public Member Functions

Interfaces accessed by local residual implementations

Member Typedef Documentation

◆ ElementResidualVector

Member Function Documentation

◆ addCCDirichletFluxDerivatives()

◆ addFluxDerivatives() [1/2]

◆ addFluxDerivatives() [2/2]

◆ addRobinFluxDerivatives()

◆ addSourceDerivatives()

◆ addStorageDerivatives()

◆ asImp() [1/2]

◆ asImp() [2/2]

◆ computeFlux()

◆ computeSource()

◆ computeStorage()

◆ evalFlux() [1/5]

◆ evalFlux() [2/5]

◆ evalFlux() [3/5]

◆ evalFlux() [4/5]

◆ evalFlux() [5/5]

◆ evalFluxAndSource()

◆ evalSource()

◆ evalStorage() [1/3]

◆ evalStorage() [2/3]

◆ evalStorage() [3/3]

◆ isStationary()

◆ problem()

◆ timeLoop()