Dumux::MimeticTwoPLocalStiffness< TypeTag > Class Template Reference

compute local stiffness matrix for conforming finite elements for the full 2-phase pressure equation More...

#include <dumux/porousmediumflow/2p/sequential/diffusion/mimetic/mimetic.hh>

Inheritance diagram for Dumux::MimeticTwoPLocalStiffness< TypeTag >:

Description

template<class TypeTag>
class Dumux::MimeticTwoPLocalStiffness< TypeTag >

compute local stiffness matrix for conforming finite elements for the full 2-phase pressure equation

Public Types
enum	{ m = 1 }
enum	{ size = 2 * dim }
using	VBlockType = Dune::FieldVector< Scalar, m >
using	MBlockType = Dune::FieldMatrix< Scalar, m, m >
using	BCBlockType = std::array< BoundaryConditions::Flags, m >

Public Member Functions
	MimeticTwoPLocalStiffness (Problem &problem, bool levelBoundaryAsDirichlet, const GridView &gridView, bool procBoundaryAsDirichlet=true)
	Constructor. More...
void	initialize ()
void	reset ()
void	setErrorInfo (Scalar maxErr, Scalar dt)
void	assemble (const Element &element, int k=1)
	Assembles local stiffness matrix for given element and order. More...
void	assemble (const Element &cell, const Dune::BlockVector< VBlockType > &localSolution, int orderOfShapeFns=1)
void	assembleBoundaryCondition (const Element &element, int k=1)
	Assembles only boundary conditions for given element. More...
template<class Vector >
void	completeRHS (const Element &element, Dune::FieldVector< int, 2 *dim > &local2Global, Vector &f)
Scalar	constructPressure (const Element &element, Dune::FieldVector< Scalar, 2 *dim > &pressTrace)
void	constructVelocity (const Element &element, Dune::FieldVector< Scalar, 2 *dim > &vel, Dune::FieldVector< Scalar, 2 *dim > &pressTrace, Scalar press)
void	constructVelocity (const Element &element, int fIdx, Scalar &vel, Dune::FieldVector< Scalar, 2 *dim > &pressTrace, Scalar press)
void	computeReconstructionMatrices (const Element &element, const Dune::FieldMatrix< Scalar, 2 *dim, 2 *dim > &W, Dune::FieldVector< Scalar, 2 *dim > &F, Scalar &dInv)
void	assembleElementMatrices (const Element &element, Dune::FieldVector< Scalar, 2 *dim > &faceVol, Dune::FieldVector< Scalar, 2 *dim > &F, Dune::FieldMatrix< Scalar, 2 *dim, 2 *dim > &Pi, Scalar &dInv, Scalar &qmean)
const Problem &	problem () const
virtual void	assemble (const Entity &e, int k=1)=0
	Assembles local stiffness matrix including boundary conditions for given element and order. More...
virtual void	assemble (const Entity &e, const Dune::BlockVector< VBlockType > &localSolution, int k=1)=0
	assemble local stiffness matrix including boundary conditions for given element and order More...
virtual void	assembleBoundaryCondition (const Entity &e, int k=1)=0
	assemble only boundary conditions for given element and order More...
void	print (std::ostream &s, int width, int precision)
	Prints contents of local stiffness matrix. More...
const MBlockType &	mat (int i, int j) const
	Accesses the local stiffness matrix. More...
const VBlockType &	rhs (int i) const
	Accesses right hand side. More...
const BoundaryTypes &	bc (int i) const
	Accesses boundary condition for each dof. More...
void	setcurrentsize (int s)
	Sets the current size of the local stiffness matrix. More...
int	currentsize ()
	Gets the current size of the local stiffness matrix. More...

Protected Attributes
Dune::Matrix< MBlockType >	A
std::vector< VBlockType >	b
std::vector< BoundaryTypes >	bctype

Member Typedef Documentation

BCBlockType

using Dumux::LocalStiffness< TypeTag, m >::BCBlockType = std::array<BoundaryConditions::Flags, m>

inherited

MBlockType

using Dumux::LocalStiffness< TypeTag, m >::MBlockType = Dune::FieldMatrix<Scalar, m, m>

inherited

VBlockType

template<class TypeTag >

using Dumux::MimeticTwoPLocalStiffness< TypeTag >::VBlockType = Dune::FieldVector<Scalar, m>

Member Enumeration Documentation

anonymous enum

template<class TypeTag >

anonymous enum

Enumerator
m

anonymous enum

template<class TypeTag >

anonymous enum

Enumerator
size

Constructor & Destructor Documentation

MimeticTwoPLocalStiffness()

template<class TypeTag >

Dumux::MimeticTwoPLocalStiffness< TypeTag >::MimeticTwoPLocalStiffness ( Problem &  problem, bool  levelBoundaryAsDirichlet, const GridView &  gridView, bool  procBoundaryAsDirichlet = true  )

inline

Constructor.

Member Function Documentation

assemble() [1/4]

template<class TypeTag >

void Dumux::MimeticTwoPLocalStiffness< TypeTag >::assemble ( const Element &  cell, const Dune::BlockVector< VBlockType > &  localSolution, int  orderOfShapeFns = 1  )

inline

assemble() [2/4]

template<class TypeTag >

void Dumux::MimeticTwoPLocalStiffness< TypeTag >::assemble ( const Element &  element, int  k = 1  )

inline

Assembles local stiffness matrix for given element and order.

On exit the following things have been done:

• The stiffness matrix for the given entity and polynomial degree has been assembled and is accessible with the mat() method.
• The boundary conditions have been evaluated and are accessible with the bc() method
• The right hand side has been assembled. It contains either the value of the essential boundary condition or the assembled source term and Neumann boundary condition. It is accessible via the rhs() method.

Parameters

element	a codim 0 entity reference
k	order of CR basis (only k = 1 is implemented)

assemble() [3/4]

virtual void Dumux::LocalStiffness< TypeTag, m >::assemble ( const Entity &  e, const Dune::BlockVector< VBlockType > &  localSolution, int  k = 1  )

pure virtualinherited

assemble local stiffness matrix including boundary conditions for given element and order

Unlike the method with only two arguments, this one additionally takes the local solution in order to allow assembly of nonlinear operators.

On exit the following things have been done:

• The stiffness matrix for the given entity and polynomial degree has been assembled and is accessible with the mat() method.
• The boundary conditions have been evaluated and are accessible with the bc() method. The boundary conditions are either Neumann, process or Dirichlet. Neumann indicates that the corresponding node (assuming a nodal basis) is at the Neumann boundary, process indicates that the node is at a process boundary (arising from the parallel decomposition of the mesh). Process boundaries are treated as homogeneous Dirichlet conditions, i.e. the corresponding value in the right hand side is set to 0. Finally, Dirichlet indicates that the node is at the Dirichlet boundary.
• The right hand side has been assembled. It contains either the value of the essential boundary condition or the assembled source term and Neumann boundary condition. It is accessible via the rhs() method.

Parameters

e	a codim 0 entity reference
localSolution	The current solution on the entity, which is needed by nonlinear assemblers
k	order of Lagrange basis (default is 1)

assemble() [4/4]

virtual void Dumux::LocalStiffness< TypeTag, m >::assemble ( const Entity &  e, int  k = 1  )

pure virtualinherited

Assembles local stiffness matrix including boundary conditions for given element and order.

On exit the following things have been done:

• The stiffness matrix for the given entity and polynomial degree has been assembled and is accessible with the mat() method.
• The boundary conditions have been evaluated and are accessible with the bc() method. The boundary conditions are either neumann, process or dirichlet. Neumann indicates that the corresponding node (assuming a nodal basis) is at the Neumann boundary, process indicates that the node is at a process boundary (arising from the parallel decomposition of the mesh). Process boundaries are treated as homogeneous Dirichlet conditions, i.e. the corresponding value in the right hand side is set to 0. Finally, Dirichlet indicates that the node is at the Dirichlet boundary.
• The right hand side has been assembled. It contains either the value of the essential boundary condition or the assembled source term and neumann boundary condition. It is accessible via the rhs() method.

Parameters

e	a codim 0 entity reference
k	order of Lagrange basis (default is 1)

assembleBoundaryCondition() [1/2]

template<class TypeTag >

void Dumux::MimeticTwoPLocalStiffness< TypeTag >::assembleBoundaryCondition ( const Element &  element, int  k = 1  )

inline

Assembles only boundary conditions for given element.

On exit the following things have been done:

• The boundary conditions have been evaluated and are accessible with the bc() method
• The right hand side contains either the value of the essential boundary condition or the assembled neumann boundary condition. It is accessible via the rhs() method.

Parameters

element	a codim 0 entity reference
k	order of CR basis

assembleBoundaryCondition() [2/2]

virtual void Dumux::LocalStiffness< TypeTag, m >::assembleBoundaryCondition ( const Entity &  e, int  k = 1  )

pure virtualinherited

assemble only boundary conditions for given element and order

On exit the following things have been done:

• The boundary conditions have been evaluated and are accessible with the bc() method. The boundary conditions are either Neumann, process or Dirichlet. Neumann indicates that the corresponding node (assuming a nodal basis) is at the Neumann boundary, process indicates that the node is at a process boundary (arising from the parallel decomposition of the mesh). Process boundaries are treated as homogeneous Dirichlet conditions, i.e. the corresponding value in the right hand side is set to 0. Finally, Dirichlet indicates that the node is at the Dirichlet boundary.
• The right hand side has been assembled as far as boundary conditions are concerned. It contains either the value of the essential boundary condition or the assembled Neumann boundary condition. It is accessible via the rhs() method.

Parameters

e	a codim 0 entity reference
k	order of Lagrange basis (default is 1)

assembleElementMatrices()

template<class TypeTag >

void Dumux::MimeticTwoPLocalStiffness< TypeTag >::assembleElementMatrices	(	const Element &	element,
		Dune::FieldVector< Scalar, 2 *dim > &	faceVol,
		Dune::FieldVector< Scalar, 2 *dim > &	F,
		Dune::FieldMatrix< Scalar, 2 *dim, 2 *dim > &	Pi,
		Scalar &	dInv,
		Scalar &	qmean
	)

bc()

const BoundaryTypes & Dumux::LocalStiffness< TypeTag, m >::bc ( int  i ) const

inlineinherited

Accesses boundary condition for each dof.

Access boundary condition type for each degree of freedom. Elements are undefined without prior call to the assemble method.

Parameters

i

index i

completeRHS()

template<class TypeTag >

template<class Vector >

void Dumux::MimeticTwoPLocalStiffness< TypeTag >::completeRHS ( const Element &  element, Dune::FieldVector< int, 2 *dim > &  local2Global, Vector &  f  )

inline

computeReconstructionMatrices()

template<class TypeTag >

void Dumux::MimeticTwoPLocalStiffness< TypeTag >::computeReconstructionMatrices ( const Element &  element, const Dune::FieldMatrix< Scalar, 2 *dim, 2 *dim > &  W, Dune::FieldVector< Scalar, 2 *dim > &  F, Scalar &  dInv  )

inline

constructPressure()

template<class TypeTag >

Scalar Dumux::MimeticTwoPLocalStiffness< TypeTag >::constructPressure ( const Element &  element, Dune::FieldVector< Scalar, 2 *dim > &  pressTrace  )

inline

constructVelocity() [1/2]

template<class TypeTag >

void Dumux::MimeticTwoPLocalStiffness< TypeTag >::constructVelocity ( const Element &  element, Dune::FieldVector< Scalar, 2 *dim > &  vel, Dune::FieldVector< Scalar, 2 *dim > &  pressTrace, Scalar  press  )

inline

constructVelocity() [2/2]

template<class TypeTag >

void Dumux::MimeticTwoPLocalStiffness< TypeTag >::constructVelocity ( const Element &  element, int  fIdx, Scalar &  vel, Dune::FieldVector< Scalar, 2 *dim > &  pressTrace, Scalar  press  )

inline

currentsize()

int Dumux::LocalStiffness< TypeTag, m >::currentsize ( )

inlineinherited

Gets the current size of the local stiffness matrix.

initialize()

template<class TypeTag >

void Dumux::MimeticTwoPLocalStiffness< TypeTag >::initialize ( )

inline

mat()

const MBlockType & Dumux::LocalStiffness< TypeTag, m >::mat ( int  i, int  j  ) const

inlineinherited

Accesses the local stiffness matrix.

Access elements of the local stiffness matrix. Elements are undefined without prior call to the assemble method.

Parameters

i	index i
j	index j

print()

void Dumux::LocalStiffness< TypeTag, m >::print ( std::ostream &  s, int  width, int  precision  )

inlineinherited

Prints contents of local stiffness matrix.

Parameters

s	output stream
width	the width
precision	the precision

problem()

template<class TypeTag >

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

inline

reset()

template<class TypeTag >

void Dumux::MimeticTwoPLocalStiffness< TypeTag >::reset ( )

inline

rhs()

const VBlockType & Dumux::LocalStiffness< TypeTag, m >::rhs ( int  i ) const

inlineinherited

Accesses right hand side.

Access elements of the right hand side vector. Elements are undefined without prior call to the assemble method.

Parameters

i

index i

setcurrentsize()

void Dumux::LocalStiffness< TypeTag, m >::setcurrentsize ( int  s )

inlineinherited

Sets the current size of the local stiffness matrix.

Parameters

s	size which is to be set

setErrorInfo()

template<class TypeTag >

void Dumux::MimeticTwoPLocalStiffness< TypeTag >::setErrorInfo ( Scalar  maxErr, Scalar  dt  )

inline

Member Data Documentation

A

Dune::Matrix<MBlockType> Dumux::LocalStiffness< TypeTag, m >::A

protectedinherited

b

std::vector<VBlockType> Dumux::LocalStiffness< TypeTag, m >::b

protectedinherited

bctype

std::vector<BoundaryTypes> Dumux::LocalStiffness< TypeTag, m >::bctype

protectedinherited

---

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

• mimetic.hh