Test problem for the 1pnc (Navier-) Stokes problem. More...

#include <test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/problem_stokes.hh>

Inheritance diagram for Dumux::StokesSubProblem< TypeTag >:

Description

template<class TypeTag>
class Dumux::StokesSubProblem< TypeTag >

Test problem for the 1pnc (Navier-) Stokes problem.

Test problem for the one-phase (Navier-) Stokes problem.

Horizontal flow from left to right with a parabolic velocity profile.

Vertical flow from top to bottom with a parabolic velocity profile.

Public Member Functions
	StokesSubProblem (std::shared_ptr< const GridGeometry > gridGeometry, std::shared_ptr< CouplingManager > couplingManager)

	StokesSubProblem (std::shared_ptr< const GridGeometry > gridGeometry, std::shared_ptr< CouplingManager > couplingManager)

const std::string &	name () const
	The problem name. More...

	StokesSubProblem (std::shared_ptr< const GridGeometry > gridGeometry, std::shared_ptr< CouplingManager > couplingManager)

const std::string &	name () const
	The problem name. More...

	StokesSubProblem (std::shared_ptr< const GridGeometry > gridGeometry, std::shared_ptr< CouplingManager > couplingManager)

	StokesSubProblem (std::shared_ptr< const GridGeometry > gridGeometry, std::shared_ptr< CouplingManager > couplingManager)

const std::string &	name () const
	The problem name. More...

	StokesSubProblem (std::shared_ptr< const GridGeometry > gridGeometry, std::shared_ptr< CouplingManager > couplingManager)

const std::string &	name () const
	The problem name. More...

Scalar	temperatureAtPos (const GlobalPosition &globalPos) const
	Returns the temperature \(\mathrm{[K]}\) at a given global position. More...

const GravityVector &	gravity () const
	Returns the acceleration due to gravity. More...

bool	enableInertiaTerms () const
	Returns whether interia terms should be considered. More...

template<class SolutionVector , class G = GridGeometry>
std::enable_if< G::discMethod==DiscretizationMethod::staggered, void >::type	applyInitialFaceSolution (SolutionVector &sol, const SubControlVolumeFace &scvf, const PrimaryVariables &initSol) const
	Applys the initial face solution (velocities on the faces). Specialization for staggered grid discretization. More...

Scalar	pseudo3DWallFriction (const Scalar velocity, const Scalar viscosity, const Scalar height, const Scalar factor=8.0) const
	An additional drag term can be included as source term for the momentum balance to mimic 3D flow behavior in 2D: More...

template<class ElementVolumeVariables , class ElementFaceVariables , class G = GridGeometry>
std::enable_if< G::discMethod==DiscretizationMethod::staggered, Scalar >::type	pseudo3DWallFriction (const SubControlVolumeFace &scvf, const ElementVolumeVariables &elemVolVars, const ElementFaceVariables &elemFaceVars, const Scalar height, const Scalar factor=8.0) const
	Convenience function for staggered grid implementation. More...

Scalar	betaBJ (const Element &element, const SubControlVolumeFace &scvf) const
	Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition. More...

Scalar	velocityPorousMedium (const Element &element, const SubControlVolumeFace &scvf) const
	Returns the velocity in the porous medium (which is 0 by default according to Saffmann). More...

const Scalar	bjsVelocity (const Element &element, const SubControlVolume &scv, const SubControlVolumeFace &faceOnPorousBoundary, const Scalar velocitySelf) const
	helper function to evaluate the slip velocity on the boundary when the Beavers-Joseph-Saffman condition is used More...

const Scalar	beaversJosephVelocity (const Element &element, const SubControlVolume &scv, const SubControlVolumeFace &faceOnPorousBoundary, const Scalar velocitySelf, const Scalar tangentialVelocityGradient) const
	helper function to evaluate the slip velocity on the boundary when the Beavers-Joseph condition is used More...

Problem parameters
bool	shouldWriteRestartFile () const

Scalar	temperature () const
	Returns the temperature within the domain in [K]. More...

NumEqVector	sourceAtPos (const GlobalPosition &globalPos) const
	Returns the sources within the domain. More...

Scalar	temperature () const
	Returns the temperature within the domain in [K]. More...

NumEqVector	sourceAtPos (const GlobalPosition &globalPos) const
	Returns the sources within the domain. More...

Scalar	temperature () const
	Returns the temperature within the domain in [K]. More...

NumEqVector	sourceAtPos (const GlobalPosition &globalPos) const
	Returns the sources within the domain. More...

bool	shouldWriteRestartFile () const

Scalar	temperature () const
	Returns the temperature within the domain in [K]. More...

NumEqVector	sourceAtPos (const GlobalPosition &globalPos) const
	Returns the sources within the domain. More...

Scalar	temperature () const
	Returns the temperature within the domain in [K]. More...

NumEqVector	sourceAtPos (const GlobalPosition &globalPos) const
	Returns the sources within the domain. More...

bool	shouldWriteRestartFile () const

Scalar	temperature () const
	Returns the temperature within the domain in [K]. More...

NumEqVector	sourceAtPos (const GlobalPosition &globalPos) const
	Returns the sources within the domain. More...

Boundary conditions
BoundaryTypes	boundaryTypes (const Element &element, const SubControlVolumeFace &scvf) const
	Specifies which kind of boundary condition should be used for which equation on a given boundary segment. More...

PrimaryVariables	dirichletAtPos (const GlobalPosition &globalPos) const
	Evaluates the boundary conditions for a Dirichlet control volume. More...

template<class ElementVolumeVariables , class ElementFaceVariables >
NumEqVector	neumann (const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFaceVariables &elemFaceVars, const SubControlVolumeFace &scvf) const
	Evaluates the boundary conditions for a Neumann control volume. More...

const CouplingManager &	couplingManager () const
	Get the coupling manager. More...

BoundaryTypes	boundaryTypes (const Element &element, const SubControlVolumeFace &scvf) const
	Specifies which kind of boundary condition should be used for which equation on a given boundary segment. More...

PrimaryVariables	dirichletAtPos (const GlobalPosition &globalPos) const
	Evaluates the boundary conditions for a Dirichlet control volume. More...

template<class ElementVolumeVariables , class ElementFaceVariables >
NumEqVector	neumann (const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFaceVariables &elemFaceVars, const SubControlVolumeFace &scvf) const
	Evaluates the boundary conditions for a Neumann control volume. More...

const CouplingManager &	couplingManager () const
	Get the coupling manager. More...

BoundaryTypes	boundaryTypes (const Element &element, const SubControlVolumeFace &scvf) const
	Specifies which kind of boundary condition should be used for which equation on a given boundary segment. More...

PrimaryVariables	dirichletAtPos (const GlobalPosition &pos) const
	Evaluates the boundary conditions for a Dirichlet control volume. More...

NumEqVector	neumann (const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFaceVariables &elemFaceVars, const SubControlVolumeFace &scvf) const
	Evaluates the boundary conditions for a Neumann control volume. More...

const CouplingManager &	couplingManager () const
	Get the coupling manager. More...

BoundaryTypes	boundaryTypes (const Element &element, const SubControlVolumeFace &scvf) const
	Specifies which kind of boundary condition should be used for which equation on a given boundary segment. More...

PrimaryVariables	dirichletAtPos (const GlobalPosition &globalPos) const
	Evaluates the boundary conditions for a Dirichlet control volume. More...

template<class ElementVolumeVariables , class ElementFaceVariables >
NumEqVector	neumann (const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFaceVariables &elemFaceVars, const SubControlVolumeFace &scvf) const
	Evaluates the boundary conditions for a Neumann control volume. More...

const CouplingManager &	couplingManager () const
	Get the coupling manager. More...

BoundaryTypes	boundaryTypes (const Element &element, const SubControlVolumeFace &scvf) const
	Specifies which kind of boundary condition should be used for which equation on a given boundary segment. More...

PrimaryVariables	dirichletAtPos (const GlobalPosition &globalPos) const
	Evaluates the boundary conditions for a Dirichlet control volume. More...

template<class ElementVolumeVariables , class ElementFaceVariables >
NumEqVector	neumann (const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFaceVariables &elemFaceVars, const SubControlVolumeFace &scvf) const
	Evaluates the boundary conditions for a Neumann control volume. More...

const CouplingManager &	couplingManager () const
	Get the coupling manager. More...

BoundaryTypes	boundaryTypes (const Element &element, const SubControlVolumeFace &scvf) const
	Specifies which kind of boundary condition should be used for which equation on a given boundary segment. More...

PrimaryVariables	dirichletAtPos (const GlobalPosition &pos) const
	Evaluates the boundary conditions for a Dirichlet control volume. More...

NumEqVector	neumann (const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFaceVariables &elemFaceVars, const SubControlVolumeFace &scvf) const
	Evaluates the boundary conditions for a Neumann control volume. More...

const CouplingManager &	couplingManager () const
	Get the coupling manager. More...

Volume terms
PrimaryVariables	initialAtPos (const GlobalPosition &globalPos) const
	Evaluates the initial value for a control volume. More...

Scalar	permeability (const Element &element, const SubControlVolumeFace &scvf) const
	Returns the intrinsic permeability of required as input parameter for the Beavers-Joseph-Saffman boundary condition. More...

Scalar	alphaBJ (const SubControlVolumeFace &scvf) const
	Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition. More...

void	setInjectionState (const bool yesNo)

bool	isInjectionPeriod () const

PrimaryVariables	initialAtPos (const GlobalPosition &globalPos) const
	Evaluates the initial value for a control volume. More...

Scalar	permeability (const Element &element, const SubControlVolumeFace &scvf) const
	Returns the intrinsic permeability of required as input parameter for the Beavers-Joseph-Saffman boundary condition. More...

Scalar	alphaBJ (const SubControlVolumeFace &scvf) const
	Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition. More...

void	setTimeLoop (TimeLoopPtr timeLoop)
	Sets the time loop pointer. More...

Scalar	time () const
	Returns the time. More...

PrimaryVariables	initialAtPos (const GlobalPosition &globalPos) const
	Evaluates the initial value for a control volume. More...

const Scalar	refVelocity () const
	Returns the reference velocity. More...

const Scalar	refPressure () const
	Returns the reference pressure. More...

const Scalar	refMoleFrac () const
	Returns the reference mass fraction. More...

const Scalar	refTemperature () const
	Returns the reference temperature. More...

void	setTimeLoop (TimeLoopPtr timeLoop)

Scalar	time () const

Scalar	permeability (const Element &element, const SubControlVolumeFace &scvf) const
	Returns the intrinsic permeability of required as input parameter for the Beavers-Joseph-Saffman boundary condition. More...

Scalar	alphaBJ (const SubControlVolumeFace &scvf) const
	Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition. More...

PrimaryVariables	initialAtPos (const GlobalPosition &globalPos) const
	Evaluates the initial value for a control volume. More...

Scalar	permeability (const Element &element, const SubControlVolumeFace &scvf) const
	Returns the intrinsic permeability of required as input parameter for the Beavers-Joseph-Saffman boundary condition. More...

Scalar	alphaBJ (const SubControlVolumeFace &scvf) const
	Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition. More...

PrimaryVariables	initialAtPos (const GlobalPosition &globalPos) const
	Evaluates the initial value for a control volume. More...

Scalar	permeability (const Element &element, const SubControlVolumeFace &scvf) const
	Returns the intrinsic permeability of required as input parameter for the Beavers-Joseph-Saffman boundary condition. More...

Scalar	alphaBJ (const SubControlVolumeFace &scvf) const
	Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition. More...

PrimaryVariables	initialAtPos (const GlobalPosition &globalPos) const
	Evaluates the initial value for a control volume. More...

Scalar	permeability (const Element &element, const SubControlVolumeFace &scvf) const
	Returns the intrinsic permeability of required as input parameter for the Beavers-Joseph-Saffman boundary condition. More...

Scalar	alphaBJ (const SubControlVolumeFace &scvf) const
	Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition. More...

Constructor & Destructor Documentation

◆ StokesSubProblem() [1/6]

template<class TypeTag >

Dumux::StokesSubProblem< TypeTag >::StokesSubProblem	(	std::shared_ptr< const GridGeometry >	gridGeometry,
		std::shared_ptr< CouplingManager >	couplingManager
	)

inline

◆ StokesSubProblem() [2/6]

template<class TypeTag >

Dumux::StokesSubProblem< TypeTag >::StokesSubProblem	(	std::shared_ptr< const GridGeometry >	gridGeometry,
		std::shared_ptr< CouplingManager >	couplingManager
	)

inline

◆ StokesSubProblem() [3/6]

template<class TypeTag >

Dumux::StokesSubProblem< TypeTag >::StokesSubProblem	(	std::shared_ptr< const GridGeometry >	gridGeometry,
		std::shared_ptr< CouplingManager >	couplingManager
	)

inline

◆ StokesSubProblem() [4/6]

template<class TypeTag >

Dumux::StokesSubProblem< TypeTag >::StokesSubProblem	(	std::shared_ptr< const GridGeometry >	gridGeometry,
		std::shared_ptr< CouplingManager >	couplingManager
	)

inline

◆ StokesSubProblem() [5/6]

template<class TypeTag >

Dumux::StokesSubProblem< TypeTag >::StokesSubProblem	(	std::shared_ptr< const GridGeometry >	gridGeometry,
		std::shared_ptr< CouplingManager >	couplingManager
	)

inline

◆ StokesSubProblem() [6/6]

template<class TypeTag >

Dumux::StokesSubProblem< TypeTag >::StokesSubProblem	(	std::shared_ptr< const GridGeometry >	gridGeometry,
		std::shared_ptr< CouplingManager >	couplingManager
	)

inline

Member Function Documentation

◆ alphaBJ() [1/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::alphaBJ ( const SubControlVolumeFace & scvf ) const

inline

Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition.

◆ alphaBJ() [2/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::alphaBJ ( const SubControlVolumeFace & scvf ) const

inline

Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition.

◆ alphaBJ() [3/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::alphaBJ ( const SubControlVolumeFace & scvf ) const

inline

Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition.

◆ alphaBJ() [4/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::alphaBJ ( const SubControlVolumeFace & scvf ) const

inline

Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition.

◆ alphaBJ() [5/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::alphaBJ ( const SubControlVolumeFace & scvf ) const

inline

Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition.

◆ alphaBJ() [6/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::alphaBJ ( const SubControlVolumeFace & scvf ) const

inline

Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition.

◆ applyInitialFaceSolution()

template<class TypeTag >

template<class SolutionVector , class G = GridGeometry>

std::enable_if< G::discMethod==DiscretizationMethod::staggered, void >::type Dumux::NavierStokesProblem< TypeTag >::applyInitialFaceSolution	(	SolutionVector &	sol,
		const SubControlVolumeFace &	scvf,
		const PrimaryVariables &	initSol
	)		const

inlineinherited

Applys the initial face solution (velocities on the faces). Specialization for staggered grid discretization.

◆ beaversJosephVelocity()

template<class TypeTag >

const Scalar Dumux::NavierStokesProblem< TypeTag >::beaversJosephVelocity	(	const Element &	element,
		const SubControlVolume &	scv,
		const SubControlVolumeFace &	faceOnPorousBoundary,
		const Scalar	velocitySelf,
		const Scalar	tangentialVelocityGradient
	)		const

inlineinherited

helper function to evaluate the slip velocity on the boundary when the Beavers-Joseph condition is used

◆ betaBJ()

template<class TypeTag >

Scalar Dumux::NavierStokesProblem< TypeTag >::betaBJ	(	const Element &	element,
		const SubControlVolumeFace &	scvf
	)		const

inlineinherited

Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition.

This member function must be overloaded in the problem implementation, if the BJS boundary condition is used.

◆ bjsVelocity()

template<class TypeTag >

const Scalar Dumux::NavierStokesProblem< TypeTag >::bjsVelocity	(	const Element &	element,
		const SubControlVolume &	scv,
		const SubControlVolumeFace &	faceOnPorousBoundary,
		const Scalar	velocitySelf
	)		const

inlineinherited

helper function to evaluate the slip velocity on the boundary when the Beavers-Joseph-Saffman condition is used

Deprecated:: "Use beaversJosephVelocity(element, scv, faceOnPorousBoundary, velocitySelf, tangentialVelocityGradient) instead"

◆ boundaryTypes() [1/6]

template<class TypeTag >

BoundaryTypes Dumux::StokesSubProblem< TypeTag >::boundaryTypes	(	const Element &	element,
		const SubControlVolumeFace &	scvf
	)		const

inline

Specifies which kind of boundary condition should be used for which equation on a given boundary segment.

Parameters

element	The finite element
scvf	The sub control volume face

◆ boundaryTypes() [2/6]

template<class TypeTag >

BoundaryTypes Dumux::StokesSubProblem< TypeTag >::boundaryTypes	(	const Element &	element,
		const SubControlVolumeFace &	scvf
	)		const

inline

Specifies which kind of boundary condition should be used for which equation on a given boundary segment.

Parameters

element	The finite element
scvf	The sub control volume face

◆ boundaryTypes() [3/6]

template<class TypeTag >

BoundaryTypes Dumux::StokesSubProblem< TypeTag >::boundaryTypes	(	const Element &	element,
		const SubControlVolumeFace &	scvf
	)		const

inline

Specifies which kind of boundary condition should be used for which equation on a given boundary segment.

Parameters

element	The finite element
scvf	The sub control volume face

◆ boundaryTypes() [4/6]

template<class TypeTag >

BoundaryTypes Dumux::StokesSubProblem< TypeTag >::boundaryTypes	(	const Element &	element,
		const SubControlVolumeFace &	scvf
	)		const

inline

Specifies which kind of boundary condition should be used for which equation on a given boundary segment.

Parameters

element	The finite element
scvf	The sub control volume face

◆ boundaryTypes() [5/6]

template<class TypeTag >

BoundaryTypes Dumux::StokesSubProblem< TypeTag >::boundaryTypes	(	const Element &	element,
		const SubControlVolumeFace &	scvf
	)		const

inline

Specifies which kind of boundary condition should be used for which equation on a given boundary segment.

Parameters

element	The finite element
scvf	The sub control volume face

◆ boundaryTypes() [6/6]

template<class TypeTag >

BoundaryTypes Dumux::StokesSubProblem< TypeTag >::boundaryTypes	(	const Element &	element,
		const SubControlVolumeFace &	scvf
	)		const

inline

Specifies which kind of boundary condition should be used for which equation on a given boundary segment.

Parameters

element	The finite element
scvf	The sub control volume face

◆ couplingManager() [1/6]

template<class TypeTag >

const CouplingManager & Dumux::StokesSubProblem< TypeTag >::couplingManager ( ) const

inline

Get the coupling manager.

◆ couplingManager() [2/6]

template<class TypeTag >

const CouplingManager & Dumux::StokesSubProblem< TypeTag >::couplingManager ( ) const

inline

Get the coupling manager.

◆ couplingManager() [3/6]

template<class TypeTag >

const CouplingManager & Dumux::StokesSubProblem< TypeTag >::couplingManager ( ) const

inline

Get the coupling manager.

◆ couplingManager() [4/6]

template<class TypeTag >

const CouplingManager & Dumux::StokesSubProblem< TypeTag >::couplingManager ( ) const

inline

Get the coupling manager.

◆ couplingManager() [5/6]

template<class TypeTag >

const CouplingManager & Dumux::StokesSubProblem< TypeTag >::couplingManager ( ) const

inline

Get the coupling manager.

◆ couplingManager() [6/6]

template<class TypeTag >

const CouplingManager & Dumux::StokesSubProblem< TypeTag >::couplingManager ( ) const

inline

Get the coupling manager.

◆ dirichletAtPos() [1/6]

template<class TypeTag >

PrimaryVariables Dumux::StokesSubProblem< TypeTag >::dirichletAtPos ( const GlobalPosition & globalPos ) const

inline

Evaluates the boundary conditions for a Dirichlet control volume.

◆ dirichletAtPos() [2/6]

template<class TypeTag >

PrimaryVariables Dumux::StokesSubProblem< TypeTag >::dirichletAtPos ( const GlobalPosition & globalPos ) const

inline

Evaluates the boundary conditions for a Dirichlet control volume.

◆ dirichletAtPos() [3/6]

template<class TypeTag >

PrimaryVariables Dumux::StokesSubProblem< TypeTag >::dirichletAtPos ( const GlobalPosition & globalPos ) const

inline

Evaluates the boundary conditions for a Dirichlet control volume.

◆ dirichletAtPos() [4/6]

template<class TypeTag >

PrimaryVariables Dumux::StokesSubProblem< TypeTag >::dirichletAtPos ( const GlobalPosition & globalPos ) const

inline

Evaluates the boundary conditions for a Dirichlet control volume.

◆ dirichletAtPos() [5/6]

template<class TypeTag >

PrimaryVariables Dumux::StokesSubProblem< TypeTag >::dirichletAtPos ( const GlobalPosition & pos ) const

inline

Evaluates the boundary conditions for a Dirichlet control volume.

◆ dirichletAtPos() [6/6]

template<class TypeTag >

PrimaryVariables Dumux::StokesSubProblem< TypeTag >::dirichletAtPos ( const GlobalPosition & pos ) const

inline

Evaluates the boundary conditions for a Dirichlet control volume.

◆ enableInertiaTerms()

template<class TypeTag >

bool Dumux::NavierStokesProblem< TypeTag >::enableInertiaTerms ( ) const

inlineinherited

Returns whether interia terms should be considered.

◆ gravity()

template<class TypeTag >

const GravityVector & Dumux::NavierStokesProblem< TypeTag >::gravity ( ) const

inlineinherited

Returns the acceleration due to gravity.

If the Problem.EnableGravity parameter is true, this means \(\boldsymbol{g} = ( 0,\dots,\ -9.81)^T \), else \(\boldsymbol{g} = ( 0,\dots, 0)^T \)

◆ initialAtPos() [1/6]

template<class TypeTag >

PrimaryVariables Dumux::StokesSubProblem< TypeTag >::initialAtPos ( const GlobalPosition & globalPos ) const

inline

Evaluates the initial value for a control volume.

Parameters

globalPos The global position

◆ initialAtPos() [2/6]

template<class TypeTag >

PrimaryVariables Dumux::StokesSubProblem< TypeTag >::initialAtPos ( const GlobalPosition & globalPos ) const

inline

Evaluates the initial value for a control volume.

Parameters

globalPos The global position

◆ initialAtPos() [3/6]

template<class TypeTag >

PrimaryVariables Dumux::StokesSubProblem< TypeTag >::initialAtPos ( const GlobalPosition & globalPos ) const

inline

Evaluates the initial value for a control volume.

Parameters

globalPos The global position

◆ initialAtPos() [4/6]

template<class TypeTag >

PrimaryVariables Dumux::StokesSubProblem< TypeTag >::initialAtPos ( const GlobalPosition & globalPos ) const

inline

Evaluates the initial value for a control volume.

Parameters

globalPos The global position

◆ initialAtPos() [5/6]

template<class TypeTag >

PrimaryVariables Dumux::StokesSubProblem< TypeTag >::initialAtPos ( const GlobalPosition & globalPos ) const

inline

Evaluates the initial value for a control volume.

Parameters

globalPos The global position

◆ initialAtPos() [6/6]

template<class TypeTag >

PrimaryVariables Dumux::StokesSubProblem< TypeTag >::initialAtPos ( const GlobalPosition & globalPos ) const

inline

Evaluates the initial value for a control volume.

Parameters

globalPos The global position

◆ isInjectionPeriod()

template<class TypeTag >

bool Dumux::StokesSubProblem< TypeTag >::isInjectionPeriod ( ) const

inline

◆ name() [1/4]

template<class TypeTag >

const std::string & Dumux::StokesSubProblem< TypeTag >::name ( ) const

inline

The problem name.

◆ name() [2/4]

template<class TypeTag >

const std::string & Dumux::StokesSubProblem< TypeTag >::name ( ) const

inline

The problem name.

◆ name() [3/4]

template<class TypeTag >

const std::string & Dumux::StokesSubProblem< TypeTag >::name ( ) const

inline

The problem name.

◆ name() [4/4]

template<class TypeTag >

const std::string & Dumux::StokesSubProblem< TypeTag >::name ( ) const

inline

The problem name.

◆ neumann() [1/6]

template<class TypeTag >

template<class ElementVolumeVariables , class ElementFaceVariables >

NumEqVector Dumux::StokesSubProblem< TypeTag >::neumann	(	const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const ElementFaceVariables &	elemFaceVars,
		const SubControlVolumeFace &	scvf
	)		const

inline

Evaluates the boundary conditions for a Neumann control volume.

Parameters

element	The element for which the Neumann boundary condition is set
fvGeometry	The fvGeometry
elemVolVars	The element volume variables
elemFaceVars	The element face variables
scvf	The boundary sub control volume face

◆ neumann() [2/6]

template<class TypeTag >

template<class ElementVolumeVariables , class ElementFaceVariables >

NumEqVector Dumux::StokesSubProblem< TypeTag >::neumann	(	const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const ElementFaceVariables &	elemFaceVars,
		const SubControlVolumeFace &	scvf
	)		const

inline

Evaluates the boundary conditions for a Neumann control volume.

Parameters

element	The element for which the Neumann boundary condition is set
fvGeometry	The fvGeometry
elemVolVars	The element volume variables
elemFaceVars	The element face variables
scvf	The boundary sub control volume face

◆ neumann() [3/6]

template<class TypeTag >

NumEqVector Dumux::StokesSubProblem< TypeTag >::neumann	(	const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const ElementFaceVariables &	elemFaceVars,
		const SubControlVolumeFace &	scvf
	)		const

inline

Evaluates the boundary conditions for a Neumann control volume.

Parameters

element	The element for which the Neumann boundary condition is set
fvGeometry	The fvGeometry
elemVolVars	The element volume variables
elemFaceVars	The element face variables
scvf	The boundary sub control volume face

◆ neumann() [4/6]

template<class TypeTag >

template<class ElementVolumeVariables , class ElementFaceVariables >

NumEqVector Dumux::StokesSubProblem< TypeTag >::neumann	(	const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const ElementFaceVariables &	elemFaceVars,
		const SubControlVolumeFace &	scvf
	)		const

inline

Evaluates the boundary conditions for a Neumann control volume.

Parameters

element	The element for which the Neumann boundary condition is set
fvGeometry	The fvGeometry
elemVolVars	The element volume variables
elemFaceVars	The element face variables
scvf	The boundary sub control volume face

◆ neumann() [5/6]

template<class TypeTag >

template<class ElementVolumeVariables , class ElementFaceVariables >

NumEqVector Dumux::StokesSubProblem< TypeTag >::neumann	(	const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const ElementFaceVariables &	elemFaceVars,
		const SubControlVolumeFace &	scvf
	)		const

inline

Evaluates the boundary conditions for a Neumann control volume.

Parameters

element	The element for which the Neumann boundary condition is set
fvGeometry	The fvGeometry
elemVolVars	The element volume variables
elemFaceVars	The element face variables
scvf	The boundary sub control volume face

◆ neumann() [6/6]

template<class TypeTag >

NumEqVector Dumux::StokesSubProblem< TypeTag >::neumann	(	const Element &	element,
		const FVElementGeometry &	fvGeometry,
		const ElementVolumeVariables &	elemVolVars,
		const ElementFaceVariables &	elemFaceVars,
		const SubControlVolumeFace &	scvf
	)		const

inline

Evaluates the boundary conditions for a Neumann control volume.

Parameters

element	The element for which the Neumann boundary condition is set
fvGeometry	The fvGeometry
elemVolVars	The element volume variables
elemFaceVars	The element face variables
scvf	The boundary sub control volume face

For this method, the values variable stores primary variables.

◆ permeability() [1/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::permeability	(	const Element &	element,
		const SubControlVolumeFace &	scvf
	)		const

inline

Returns the intrinsic permeability of required as input parameter for the Beavers-Joseph-Saffman boundary condition.

◆ permeability() [2/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::permeability	(	const Element &	element,
		const SubControlVolumeFace &	scvf
	)		const

inline

Returns the intrinsic permeability of required as input parameter for the Beavers-Joseph-Saffman boundary condition.

◆ permeability() [3/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::permeability	(	const Element &	element,
		const SubControlVolumeFace &	scvf
	)		const

inline

Returns the intrinsic permeability of required as input parameter for the Beavers-Joseph-Saffman boundary condition.

◆ permeability() [4/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::permeability	(	const Element &	element,
		const SubControlVolumeFace &	scvf
	)		const

inline

Returns the intrinsic permeability of required as input parameter for the Beavers-Joseph-Saffman boundary condition.

◆ permeability() [5/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::permeability	(	const Element &	element,
		const SubControlVolumeFace &	scvf
	)		const

inline

Returns the intrinsic permeability of required as input parameter for the Beavers-Joseph-Saffman boundary condition.

◆ permeability() [6/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::permeability	(	const Element &	element,
		const SubControlVolumeFace &	scvf
	)		const

inline

Returns the intrinsic permeability of required as input parameter for the Beavers-Joseph-Saffman boundary condition.

◆ pseudo3DWallFriction() [1/2]

template<class TypeTag >

Scalar Dumux::NavierStokesProblem< TypeTag >::pseudo3DWallFriction	(	const Scalar	velocity,
		const Scalar	viscosity,
		const Scalar	height,
		const Scalar	factor = `8.0`
	)		const

inlineinherited

An additional drag term can be included as source term for the momentum balance to mimic 3D flow behavior in 2D:

\[ f_{drag} = -(8 \mu / h^2)v \]

Here, \(h\) corresponds to the extruded height that is bounded by the imaginary walls. See Flekkoy et al. (1995) [25]
A value of 8.0 is used as a default factor, corresponding to the velocity profile at the center plane of the virtual height (maximum velocity). Setting this value to 12.0 corresponds to an depth-averaged velocity (Venturoli and Boek, 2006) [74].

◆ pseudo3DWallFriction() [2/2]

template<class TypeTag >

template<class ElementVolumeVariables , class ElementFaceVariables , class G = GridGeometry>

std::enable_if< G::discMethod==DiscretizationMethod::staggered, Scalar >::type Dumux::NavierStokesProblem< TypeTag >::pseudo3DWallFriction	(	const SubControlVolumeFace &	scvf,
		const ElementVolumeVariables &	elemVolVars,
		const ElementFaceVariables &	elemFaceVars,
		const Scalar	height,
		const Scalar	factor = `8.0`
	)		const

inlineinherited

Convenience function for staggered grid implementation.

◆ refMoleFrac()

template<class TypeTag >

const Scalar Dumux::StokesSubProblem< TypeTag >::refMoleFrac ( ) const

inline

Returns the reference mass fraction.

◆ refPressure()

template<class TypeTag >

const Scalar Dumux::StokesSubProblem< TypeTag >::refPressure ( ) const

inline

Returns the reference pressure.

◆ refTemperature()

template<class TypeTag >

const Scalar Dumux::StokesSubProblem< TypeTag >::refTemperature ( ) const

inline

Returns the reference temperature.

◆ refVelocity()

template<class TypeTag >

const Scalar Dumux::StokesSubProblem< TypeTag >::refVelocity ( ) const

inline

Returns the reference velocity.

◆ setInjectionState()

template<class TypeTag >

void Dumux::StokesSubProblem< TypeTag >::setInjectionState ( const bool yesNo )

inline

◆ setTimeLoop() [1/2]

template<class TypeTag >

void Dumux::StokesSubProblem< TypeTag >::setTimeLoop ( TimeLoopPtr timeLoop )

inline

Sets the time loop pointer.

◆ setTimeLoop() [2/2]

template<class TypeTag >

void Dumux::StokesSubProblem< TypeTag >::setTimeLoop ( TimeLoopPtr timeLoop )

inline

◆ shouldWriteRestartFile() [1/3]

template<class TypeTag >

bool Dumux::StokesSubProblem< TypeTag >::shouldWriteRestartFile ( ) const

inline

◆ shouldWriteRestartFile() [2/3]

template<class TypeTag >

bool Dumux::StokesSubProblem< TypeTag >::shouldWriteRestartFile ( ) const

inline

◆ shouldWriteRestartFile() [3/3]

template<class TypeTag >

bool Dumux::StokesSubProblem< TypeTag >::shouldWriteRestartFile ( ) const

inline

◆ sourceAtPos() [1/6]

template<class TypeTag >

NumEqVector Dumux::StokesSubProblem< TypeTag >::sourceAtPos ( const GlobalPosition & globalPos ) const

inline

Returns the sources within the domain.

Parameters

globalPos The global position

◆ sourceAtPos() [2/6]

template<class TypeTag >

NumEqVector Dumux::StokesSubProblem< TypeTag >::sourceAtPos ( const GlobalPosition & globalPos ) const

inline

Returns the sources within the domain.

Parameters

globalPos The global position

◆ sourceAtPos() [3/6]

template<class TypeTag >

NumEqVector Dumux::StokesSubProblem< TypeTag >::sourceAtPos ( const GlobalPosition & globalPos ) const

inline

Returns the sources within the domain.

Parameters

globalPos The global position

◆ sourceAtPos() [4/6]

template<class TypeTag >

NumEqVector Dumux::StokesSubProblem< TypeTag >::sourceAtPos ( const GlobalPosition & globalPos ) const

inline

Returns the sources within the domain.

Parameters

globalPos The global position

◆ sourceAtPos() [5/6]

template<class TypeTag >

NumEqVector Dumux::StokesSubProblem< TypeTag >::sourceAtPos ( const GlobalPosition & globalPos ) const

inline

Returns the sources within the domain.

Parameters

globalPos The global position

◆ sourceAtPos() [6/6]

template<class TypeTag >

NumEqVector Dumux::StokesSubProblem< TypeTag >::sourceAtPos ( const GlobalPosition & globalPos ) const

inline

Returns the sources within the domain.

Parameters

globalPos The global position

◆ temperature() [1/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::temperature ( ) const

inline

Returns the temperature within the domain in [K].

This problem assumes a temperature of 10 degrees Celsius.

◆ temperature() [2/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::temperature ( ) const

inline

Returns the temperature within the domain in [K].

This problem assumes a temperature of 10 degrees Celsius.

◆ temperature() [3/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::temperature ( ) const

inline

Returns the temperature within the domain in [K].

◆ temperature() [4/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::temperature ( ) const

inline

Returns the temperature within the domain in [K].

This problem assumes a temperature of 10 degrees Celsius.

◆ temperature() [5/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::temperature ( ) const

inline

Returns the temperature within the domain in [K].

This problem assumes a temperature of 10 degrees Celsius.

◆ temperature() [6/6]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::temperature ( ) const

inline

Returns the temperature within the domain in [K].

◆ temperatureAtPos()

template<class TypeTag >

Scalar Dumux::NavierStokesProblem< TypeTag >::temperatureAtPos ( const GlobalPosition & globalPos ) const

inlineinherited

Returns the temperature \(\mathrm{[K]}\) at a given global position.

This is not specific to the discretization. By default it just calls temperature().

Parameters

globalPos The position in global coordinates where the temperature should be specified.

◆ time() [1/2]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::time ( ) const

inline

Returns the time.

◆ time() [2/2]

template<class TypeTag >

Scalar Dumux::StokesSubProblem< TypeTag >::time ( ) const

inline

◆ velocityPorousMedium()

template<class TypeTag >

Scalar Dumux::NavierStokesProblem< TypeTag >::velocityPorousMedium	(	const Element &	element,
		const SubControlVolumeFace &	scvf
	)		const

inlineinherited

Returns the velocity in the porous medium (which is 0 by default according to Saffmann).

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

Description

Public Member Functions

Volume terms

Constructor & Destructor Documentation

◆ StokesSubProblem() [1/6]

◆ StokesSubProblem() [2/6]

◆ StokesSubProblem() [3/6]

◆ StokesSubProblem() [4/6]

◆ StokesSubProblem() [5/6]

◆ StokesSubProblem() [6/6]

Member Function Documentation

◆ alphaBJ() [1/6]

◆ alphaBJ() [2/6]

◆ alphaBJ() [3/6]

◆ alphaBJ() [4/6]

◆ alphaBJ() [5/6]

◆ alphaBJ() [6/6]

◆ applyInitialFaceSolution()

◆ beaversJosephVelocity()

◆ betaBJ()

◆ bjsVelocity()

◆ boundaryTypes() [1/6]

◆ boundaryTypes() [2/6]

◆ boundaryTypes() [3/6]

◆ boundaryTypes() [4/6]

◆ boundaryTypes() [5/6]

◆ boundaryTypes() [6/6]

◆ couplingManager() [1/6]

◆ couplingManager() [2/6]

◆ couplingManager() [3/6]

◆ couplingManager() [4/6]

◆ couplingManager() [5/6]

◆ couplingManager() [6/6]

◆ dirichletAtPos() [1/6]

◆ dirichletAtPos() [2/6]

◆ dirichletAtPos() [3/6]

◆ dirichletAtPos() [4/6]

◆ dirichletAtPos() [5/6]

◆ dirichletAtPos() [6/6]

◆ enableInertiaTerms()

◆ gravity()

◆ initialAtPos() [1/6]

◆ initialAtPos() [2/6]

◆ initialAtPos() [3/6]

◆ initialAtPos() [4/6]

◆ initialAtPos() [5/6]

◆ initialAtPos() [6/6]

◆ isInjectionPeriod()

◆ name() [1/4]

◆ name() [2/4]

◆ name() [3/4]

◆ name() [4/4]

◆ neumann() [1/6]

◆ neumann() [2/6]

◆ neumann() [3/6]

◆ neumann() [4/6]

◆ neumann() [5/6]

◆ neumann() [6/6]

◆ permeability() [1/6]

◆ permeability() [2/6]

◆ permeability() [3/6]

◆ permeability() [4/6]

◆ permeability() [5/6]

◆ permeability() [6/6]

◆ pseudo3DWallFriction() [1/2]

◆ pseudo3DWallFriction() [2/2]

◆ refMoleFrac()

◆ refPressure()

◆ refTemperature()

◆ refVelocity()

◆ setInjectionState()

◆ setTimeLoop() [1/2]

◆ setTimeLoop() [2/2]

◆ shouldWriteRestartFile() [1/3]

◆ shouldWriteRestartFile() [2/3]

◆ shouldWriteRestartFile() [3/3]

◆ sourceAtPos() [1/6]

◆ sourceAtPos() [2/6]

◆ sourceAtPos() [3/6]

◆ sourceAtPos() [4/6]