#include <transportelement.h>

Inheritance diagram for oofem::TransportElement:

Collaboration diagram for oofem::TransportElement:

Public Types
enum	ElementMode { HeatTransferEM , HeatMass1TransferEM , Mass1TransferEM }

Public Member Functions
	TransportElement (int n, Domain *d, ElementMode em=HeatTransferEM)
void	giveCharacteristicMatrix (FloatMatrix &answer, CharType type, TimeStep *tStep) override
void	giveCharacteristicVector (FloatArray &answer, CharType type, ValueModeType mode, TimeStep *tStep) override
void	giveDofManDofIDMask (int inode, IntArray &answer) const override
virtual void	computeInternalForcesVector (FloatArray &answer, TimeStep *tStep)
virtual void	computeExternalForcesVector (FloatArray &answer, TimeStep *tStep, ValueModeType mode)
virtual void	computeInertiaForcesVector (FloatArray &answer, TimeStep *tStep)
virtual void	computeLumpedCapacityVector (FloatArray &answer, TimeStep *tStep)
void	computeLoadVector (FloatArray &answer, BodyLoad load, CharType type, ValueModeType mode, TimeStep tStep) override
void	computeBoundarySurfaceLoadVector (FloatArray &answer, BoundaryLoad load, int boundary, CharType type, ValueModeType mode, TimeStep tStep, bool global=true) override
void	computeTangentFromSurfaceLoad (FloatMatrix &answer, BoundaryLoad load, int boundary, MatResponseMode rmode, TimeStep tStep) override
void	computeTangentFromEdgeLoad (FloatMatrix &answer, BoundaryLoad load, int boundary, MatResponseMode rmode, TimeStep tStep) override
void	computeBoundaryEdgeLoadVector (FloatArray &answer, BoundaryLoad load, int edge, CharType type, ValueModeType mode, TimeStep tStep, bool global=true) override
void	computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer) override
TransportCrossSection *	giveTransportCrossSection ()
Material *	giveMaterial () override
void	initializeFrom (InputRecord &ir, int priority) override
void	postInitialize () override
	Performs post initialization steps.
const char *	giveClassName () const override
virtual double	giveThicknessAt (const FloatArray &gcoords)
virtual void	computeCapacityMatrix (FloatMatrix &answer, TimeStep *tStep)
virtual void	computeConductivityMatrix (FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
virtual void	computeBCVectorAt (FloatArray &answer, TimeStep *tStep, ValueModeType mode)
virtual void	computeBCMtrxAt (FloatMatrix &answer, TimeStep *tStep, ValueModeType mode)
virtual void	computeInternalSourceRhsVectorAt (FloatArray &answer, TimeStep *tStep, ValueModeType mode)
virtual void	computeIntSourceLHSMatrix (FloatMatrix &answer, TimeStep *tStep)
virtual void	computeIntSourceLHSSubMatrix (FloatMatrix &answer, MatResponseMode rmode, int iri, TimeStep *tStep)
virtual void	computeFlow (FloatArray &answer, GaussPoint gp, TimeStep tStep)
virtual double	computeVof (TimeStep *tStep)
void	updateInternalState (TimeStep *tStep) override
int	checkConsistency () override
virtual int	EIPrimaryFieldI_evaluateFieldVectorAt (FloatArray &answer, PrimaryField &pf, const FloatArray &coords, IntArray &dofId, ValueModeType mode, TimeStep *tStep) override
int	giveInternalStateAtNode (FloatArray &answer, InternalStateType type, InternalStateMode mode, int node, TimeStep *tStep) override
virtual void	computeConstitutiveMatrixAt (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep)
Public Member Functions inherited from oofem::Element
	Element (int n, Domain *aDomain)
	Element (const Element &src)=delete
Element &	operator= (const Element &src)=delete
virtual	~Element ()
	Virtual destructor.
void	giveLocationArray (IntArray &locationArray, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) const
void	giveLocationArray (IntArray &locationArray, const IntArray &dofIDMask, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) const
virtual void	giveBoundaryLocationArray (IntArray &locationArray, const IntArray &bNodes, const UnknownNumberingScheme &s, IntArray *dofIds=NULL)
virtual void	giveBoundaryLocationArray (IntArray &locationArray, const IntArray &bNodes, const IntArray &dofIDMask, const UnknownNumberingScheme &s, IntArray *dofIds=NULL)
virtual int	giveNumberOfDofs ()
virtual int	giveNumberOfInternalDofManagers () const
virtual DofManager *	giveInternalDofManager (int i) const
virtual void	setInternalDofManager (int num, std::unique_ptr< DofManager > dm)
virtual double	giveCharacteristicValue (CharType type, TimeStep *tStep)
const IntArray &	giveBodyLoadList () const
const IntArray &	giveBoundaryLoadList () const
void	computeVectorOf (ValueModeType u, TimeStep *tStep, FloatArray &answer)
void	computeVectorOf (const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false)
void	computeBoundaryVectorOf (const IntArray &bNodes, const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false)
void	computeVectorOf (PrimaryField &field, const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false)
void	computeVectorOfPrescribed (ValueModeType u, TimeStep *tStep, FloatArray &answer)
void	computeVectorOfPrescribed (const IntArray &dofIDMask, ValueModeType type, TimeStep *tStep, FloatArray &answer)
virtual int	computeNumberOfDofs ()
virtual int	computeNumberOfGlobalDofs ()
int	computeNumberOfPrimaryMasterDofs ()
virtual bool	computeGtoLRotationMatrix (FloatMatrix &answer)
virtual bool	giveRotationMatrix (FloatMatrix &answer)
virtual bool	computeDofTransformationMatrix (FloatMatrix &answer, const IntArray &nodes, bool includeInternal)
virtual void	giveInternalDofManDofIDMask (int inode, IntArray &answer) const
virtual void	giveElementDofIDMask (IntArray &answer) const
virtual double	computeVolumeAround (GaussPoint *gp)
virtual double	computeVolumeAreaOrLength ()
	Computes the volume, area or length of the element depending on its spatial dimension.
double	computeMeanSize ()
virtual double	computeVolume ()
virtual double	computeArea ()
virtual double	computeLength ()
virtual IntArray	giveBoundaryEdgeNodes (int boundary, bool includeHierarchical=false) const
virtual IntArray	giveBoundarySurfaceNodes (int boundary, bool includeHierarchical=false) const
virtual IntArray	giveBoundaryNodes (int boundary) const
virtual std::unique_ptr< IntegrationRule >	giveBoundaryEdgeIntegrationRule (int order, int boundary)
virtual std::unique_ptr< IntegrationRule >	giveBoundarySurfaceIntegrationRule (int order, int boundary)
int	giveDofManagerNumber (int i) const
const IntArray &	giveDofManArray () const
void	addDofManager (DofManager *dMan)
DofManager *	giveDofManager (int i) const
Node *	giveNode (int i) const
virtual ElementSide *	giveSide (int i) const
virtual FEInterpolation *	giveInterpolation () const
virtual FEInterpolation *	giveInterpolation (DofIDItem id) const
virtual const FEInterpolation *	getGeometryInterpolation () const
int	giveMaterialNumber () const
CrossSection *	giveCrossSection ()
int	getActivityTimeFunctionNumber ()
void	setActivityTimeFunctionNumber (int funcIndx)
void	setMaterial (int matIndx)
virtual void	setCrossSection (int csIndx)
virtual int	giveNumberOfDofManagers () const
void	setNumberOfDofManagers (int i)
	Sets number of element dof managers.
virtual int	giveNumberOfNodes () const
void	setDofManagers (const IntArray &dmans)
void	setDofManager (int id, int dm)
void	setBodyLoads (const IntArray &bodyLoads)
void	setIntegrationRules (std ::vector< std ::unique_ptr< IntegrationRule > > irlist)
virtual integrationDomain	giveIntegrationDomain () const
virtual MaterialMode	giveMaterialMode ()
virtual int	giveIntegrationRuleLocalCodeNumbers (IntArray &answer, IntegrationRule &ie)
int	giveRegionNumber ()
virtual void	updateYourself (TimeStep *tStep)
virtual void	initializeYourself (TimeStep *timeStepWhenICApply)
virtual bool	isActivated (TimeStep *tStep)
virtual bool	isCast (TimeStep *tStep)
virtual void	initForNewStep ()
virtual Element_Geometry_Type	giveGeometryType () const =0
virtual Element_Geometry_Type	giveEdgeGeometryType (int id) const
	Returns the receiver edge geometry type.
virtual Element_Geometry_Type	giveSurfaceGeometryType (int id) const
	Returns the receiver surface geometry type.
virtual int	giveSpatialDimension ()
virtual int	giveNumberOfBoundarySides ()
	Returns number of boundaries (entities of element_dimension-1: points, edges, surfaces).
virtual int	giveNumberOfEdges () const
virtual int	giveNumberOfSurfaces () const
virtual int	giveDefaultIntegrationRule () const
virtual IntegrationRule *	giveDefaultIntegrationRulePtr ()
int	giveNumberOfIntegrationRules ()
virtual IntegrationRule *	giveIntegrationRule (int i)
std::vector< std ::unique_ptr< IntegrationRule > > &	giveIntegrationRulesArray ()
virtual int	testElementExtension (ElementExtension ext)
virtual int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep)
int	giveGlobalIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep)
virtual double	giveLengthInDir (const FloatArray &normalToCrackPlane)
virtual double	giveCharacteristicLength (const FloatArray &normalToCrackPlane)
double	giveCharacteristicLengthForPlaneElements (const FloatArray &normalToCrackPlane)
double	giveCharacteristicLengthForAxisymmElements (const FloatArray &normalToCrackPlane)
virtual double	giveCharacteristicSize (GaussPoint *gp, FloatArray &normalToCrackPlane, ElementCharSizeMethod method)
virtual double	giveParentElSize () const
virtual void	updateBeforeNonlocalAverage (TimeStep *tStep)
virtual int	computeGlobalCoordinates (FloatArray &answer, const FloatArray &lcoords)
virtual bool	computeLocalCoordinates (FloatArray &answer, const FloatArray &gcoords)
virtual int	giveLocalCoordinateSystem (FloatMatrix &answer)
virtual void	giveLocalCoordinateSystemVector (InternalStateType isttype, FloatArray &answer)
virtual void	computeMidPlaneNormal (FloatArray &answer, const GaussPoint *gp)
virtual int	adaptiveMap (Domain oldd, TimeStep tStep)
virtual int	mapStateVariables (Domain &iOldDom, const TimeStep &iTStep)
virtual int	adaptiveUpdate (TimeStep *tStep)
virtual int	adaptiveFinish (TimeStep *tStep)
void	updateLocalNumbering (EntityRenumberingFunctor &f) override
template<class T>
void	ipEvaluator (T src, void(T ::f)(GaussPoint *gp))
	Integration point evaluator, loops over receiver IP's and calls given function (passed as f parameter) on them. The IP is parameter to function f.
template<class T, class S>
void	ipEvaluator (T src, void(T ::f)(GaussPoint *, S &), S &_val)
	Integration point evaluator, loops over receiver IP's and calls given function (passed as f parameter) on them. The IP is parameter to function f as well as additional array.
virtual void	drawYourself (oofegGraphicContext &gc, TimeStep *tStep)
virtual void	drawAnnotation (oofegGraphicContext &gc, TimeStep *tStep)
virtual void	drawRawGeometry (oofegGraphicContext &gc, TimeStep *tStep)
virtual void	drawDeformedGeometry (oofegGraphicContext &gc, TimeStep *tStep, UnknownType)
virtual void	drawScalar (oofegGraphicContext &gc, TimeStep *tStep)
virtual void	drawSpecial (oofegGraphicContext &gc, TimeStep *tStep)
virtual void	giveLocalIntVarMaxMin (oofegGraphicContext &gc, TimeStep *tStep, double &emin, double &emax)
virtual int	giveInternalStateAtSide (FloatArray &answer, InternalStateType type, InternalStateMode mode, int side, TimeStep *tStep)
virtual void	showSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep)
	Shows sparse structure.
virtual void	showExtendedSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep)
	Shows extended sparse structure (for example, due to nonlocal interactions for tangent stiffness).
int	giveLabel () const
int	giveGlobalNumber () const
void	setGlobalNumber (int num)
elementParallelMode	giveParallelMode () const
void	setParallelMode (elementParallelMode _mode)
	Sets parallel mode of element.
virtual elementParallelMode	giveKnotSpanParallelMode (int) const
int	packUnknowns (DataStream &buff, TimeStep *tStep)
int	unpackAndUpdateUnknowns (DataStream &buff, TimeStep *tStep)
int	estimatePackSize (DataStream &buff)
const IntArray *	givePartitionList () const
void	setPartitionList (IntArray &pl)
virtual double	predictRelativeComputationalCost ()
virtual double	giveRelativeSelfComputationalCost ()
virtual double	predictRelativeRedistributionCost ()
IntArray *	giveBodyLoadArray ()
	Returns array containing load numbers of loads acting on element.
IntArray *	giveBoundaryLoadArray ()
	Returns array containing load numbers of boundary loads acting on element.
void	initializeFinish () override
void	giveInputRecord (DynamicInputRecord &input) override
void	saveContext (DataStream &stream, ContextMode mode) override
void	restoreContext (DataStream &stream, ContextMode mode) override
void	printOutputAt (FILE file, TimeStep tStep) override
virtual const IntArray	giveLocation ()
virtual void	recalculateCoordinates (int nodeNumber, FloatArray &coords)
void	setSharedEdgeID (int iedge, int globalID)
void	setSharedSurfaceID (int isurf, int globalID)
const IntArray *	giveSharedEdgeIDs () const
const IntArray *	giveSharedSurfaceIDs () const
Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
virtual	~FEMComponent ()=default
	Virtual destructor.
virtual const char *	giveInputRecordName () const =0
Domain *	giveDomain () const
virtual void	setDomain (Domain *d)
int	giveNumber () const
void	setNumber (int num)
virtual void	initializeFrom (InputRecord &ir)
virtual void	printYourself ()
	Prints receiver state on stdout. Useful for debugging.
virtual Interface *	giveInterface (InterfaceType t)
std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros).
Public Member Functions inherited from oofem::EIPrimaryFieldInterface
	EIPrimaryFieldInterface ()
Public Member Functions inherited from oofem::Interface
	Interface ()
	Constructor.
virtual	~Interface ()

Protected Member Functions
virtual void	computeCapacitySubMatrix (FloatMatrix &answer, MatResponseMode rmode, int iri, TimeStep *tStep)
virtual void	computeConductivitySubMatrix (FloatMatrix &answer, int iri, MatResponseMode rmode, TimeStep *tStep)
void	computeBCSubVectorAt (FloatArray &answer, TimeStep *tStep, ValueModeType mode, int indx)
void	computeBCSubMtrxAt (FloatMatrix &answer, TimeStep *tStep, ValueModeType mode, int indx)
void	computeBodyBCSubVectorAt (FloatArray &answer, Load load, TimeStep tStep, ValueModeType mode, int indx)
void	computeEdgeBCSubVectorAt (FloatArray &answer, Load load, int iEdge, TimeStep tStep, ValueModeType mode, int indx)
void	computeSurfaceBCSubVectorAt (FloatArray &answer, Load load, int iSurf, TimeStep tStep, ValueModeType mode, int indx)
virtual void	computeNAt (FloatArray &answer, const FloatArray &lcoord)
virtual void	computeNmatrixAt (FloatMatrix &answer, const FloatArray &lcoords)
virtual void	computeBmatrixAt (FloatMatrix &answer, const FloatArray &lcoords)
virtual void	computeGradientMatrixAt (FloatMatrix &answer, const FloatArray &lcoords)
virtual void	computeInternalSourceRhsSubVectorAt (FloatArray &answer, TimeStep *, ValueModeType mode, int indx)
virtual void	computeEgdeNAt (FloatArray &answer, int iEdge, const FloatArray &lcoord)
virtual void	giveEdgeDofMapping (IntArray &mask, int iEdge)
virtual double	computeEdgeVolumeAround (GaussPoint *gp, int iEdge)
virtual void	computeSurfaceNAt (FloatArray &answer, int iSurf, const FloatArray &lcoord)
virtual double	computeSurfaceVolumeAround (GaussPoint *gp, int iSurf)
virtual void	giveSurfaceDofMapping (IntArray &mask, int iSurf)
virtual int	giveApproxOrder (int unknownIndx)
void	assembleLocalContribution (FloatMatrix &answer, FloatMatrix &src, int ndofs, int rdof, int cdof)
void	assembleLocalContribution (FloatArray &answer, FloatArray &src, int ndofs, int rdof)
virtual void	computeGaussPoints ()

Protected Attributes
ElementMode	emode
int	vofFunction = 0
	Fuction determining the relative volume of reference material in element.
Protected Attributes inherited from oofem::Element
int	numberOfDofMans
	Number of dofmanagers.
IntArray	dofManArray
	Array containing dofmanager numbers.
int	material
	Number of associated material.
int	crossSection
	Number of associated cross section.
IntArray	bodyLoadArray
IntArray	boundaryLoadArray
std::vector< std ::unique_ptr< IntegrationRule > >	integrationRulesArray
FloatMatrix	elemLocalCS
	Transformation material matrix, used in orthotropic and anisotropic materials, global->local transformation.
int	activityTimeFunction
	Element activity time function. If defined, nonzero value indicates active receiver, zero value inactive element.
int	globalNumber
int	numberOfGaussPoints
elementParallelMode	parallel_mode
	Determines the parallel mode of the element.
IntArray	partitions
IntArray	globalEdgeIDs
IntArray	globalSurfaceIDs
Protected Attributes inherited from oofem::FEMComponent
int	number
	Component number.
Domain *	domain
	Link to domain object, useful for communicating with other FEM components.

Static Protected Attributes
static const double	stefanBoltzmann = 5.67e-8
	Stefan–Boltzmann constant W/m2/K4.
static ParamKey	IPK_TransportElement_vof_function

Additional Inherited Members
Static Public Attributes inherited from oofem::Element
static ParamKey	IPK_Element_mat
static ParamKey	IPK_Element_crosssect
static ParamKey	IPK_Element_nodes
static ParamKey	IPK_Element_bodyload
static ParamKey	IPK_Element_boundaryload
static ParamKey	IPK_Element_lcs
static ParamKey	IPK_Element_partitions
static ParamKey	IPK_Element_remote
static ParamKey	IPK_Element_activityTimeFunction
static ParamKey	IPK_Element_nip

Detailed Description

This abstract class represent a general base element class for transport problems. In actual implementation, the same approximation order of all unknowns is assumed, but this can be easily implemented.

Definition at line 54 of file transportelement.h.

Member Enumeration Documentation

◆ ElementMode

enum oofem::TransportElement::ElementMode

Enumerator
HeatTransferEM
HeatMass1TransferEM
Mass1TransferEM

Definition at line 57 of file transportelement.h.

Constructor & Destructor Documentation

◆ TransportElement()

oofem::TransportElement::TransportElement	(	int	n,
		Domain *	d,
		ElementMode	em = HeatTransferEM )

Definition at line 75 of file transportelement.C.

References oofem::Element::Element(), emode, and vofFunction.

Referenced by oofem::Brick1_ht::Brick1_ht(), oofem::LatticeTransportElement::LatticeTransportElement(), oofem::Line1_ht::Line1_ht(), oofem::QBrick1_ht::QBrick1_ht(), oofem::QQuad1_ht::QQuad1_ht(), oofem::Quad1_ht::Quad1_ht(), oofem::QWedge_ht::QWedge_ht(), oofem::Tetrah1_ht::Tetrah1_ht(), oofem::Tr1_ht::Tr1_ht(), oofem::Tr1Darcy::Tr1Darcy(), and oofem::Wedge_ht::Wedge_ht().

Member Function Documentation

◆ assembleLocalContribution() [1/2]

void oofem::TransportElement::assembleLocalContribution	(	FloatArray &	answer,
		FloatArray &	src,
		int	ndofs,
		int	rdof )

protected

Assembles the given source vector of size (ndofs) into target answer. The coefficients in src vector are assumed to represent local sub-vector for specific DOF. The answer vector is element global one, and it typically contains the values for different DOFs. In this routine is assumed, that the number of dofs per node is the same. The DOFs in answer are ordered in a same way as the DOFS at the element level are (all Dofs corresponding to given node subsequently, followed by next node DOFS).

Parameters

answer	Receiver vector.
src	Source vector containing local subvector corresponding to single DOF.
ndofs	Number of DOFs per node (assumed same for all nodes).
rdof	Rows of src are localized into rows of answer corresponding to rdof-th dof.

Definition at line 1253 of file transportelement.C.

References oofem::FloatArray::at(), and oofem::Element::giveNumberOfDofManagers().

◆ assembleLocalContribution() [2/2]

void oofem::TransportElement::assembleLocalContribution	(	FloatMatrix &	answer,
		FloatMatrix &	src,
		int	ndofs,
		int	rdof,
		int	cdof )

protected

Assembles the given source matrix of size (ndofs, ndofs) into target matrix answer. The coefficients in src matrix are assumed to represent local sub-matrix for specific DOF. The answer matrix is element global one, and it typically contains the values for different DOFs. In this routine is assumed, that the number of dofs per node is the same. The DOFs in answer are ordered in a same way as the DOFS at the element level are (all Dofs corresponding to given node subsequently, followed by next node DOFS).

Parameters

answer	Receiver matrix.
src	Source matrix containing local sub matrix corresponding to single DOF.
ndofs	Number of DOFs per node (assumed same for all nodes).
rdof	Rows of src are localized into rows of answer corresponding to rdof-th dof.
cdof	Columns of src are localized into columns of answer corresponding to cdof-th dof.

Definition at line 1237 of file transportelement.C.

References oofem::FloatMatrix::at(), and oofem::Element::giveNumberOfDofManagers().

Referenced by computeBCMtrxAt(), computeBCVectorAt(), computeCapacityMatrix(), computeConductivityMatrix(), computeInternalSourceRhsVectorAt(), and computeIntSourceLHSMatrix().

◆ checkConsistency()

int oofem::TransportElement::checkConsistency ( void )

overridevirtual

Performs consistency check. This method is called at startup for all elements in particular domain. This method is intended to check data compatibility. Particular element types should test if compatible material and crossSection both with required capabilities are specified. Derived classes should provide their own analysis specific tests. Some printed input if incompatibility is found should be provided (error or warning member functions). Method can be also used to initialize some variables, since this is invoked after all domain components are instanciated.

Returns: Zero value if check fail, otherwise nonzero.

Reimplemented from oofem::Element.

Definition at line 167 of file transportelement.C.

References oofem::Element::giveCrossSection(), giveMaterial(), oofem::Element::material, and OOFEM_WARNING.

◆ computeBCMtrxAt()

void oofem::TransportElement::computeBCMtrxAt	(	FloatMatrix &	answer,
		TimeStep *	tStep,
		ValueModeType	mode )

virtual

Computes the LHS contribution to balance equation(s) due to boundary conditions

Definition at line 894 of file transportelement.C.

References assembleLocalContribution(), computeBCSubMtrxAt(), oofem::Element::computeNumberOfDofs(), emode, HeatMass1TransferEM, HeatTransferEM, oofem::FloatMatrix::isNotEmpty(), Mass1TransferEM, OOFEM_ERROR, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by giveCharacteristicMatrix().

◆ computeBCSubMtrxAt()

void oofem::TransportElement::computeBCSubMtrxAt	(	FloatMatrix &	answer,
		TimeStep *	tStep,
		ValueModeType	mode,
		int	indx )

protected

Computes the part of LHS due to applied BCs. the result should be localized.

Parameters

answer	Contains the result.
tStep	Solution step.
mode	Mode of value (incremental, total, ...).
indx	Unknown index.

Todo: Deal with coupled fields

Todo: Deal with coupled fields

Definition at line 1125 of file transportelement.C.

References oofem::FloatMatrix::assemble(), oofem::Element::boundaryLoadArray, oofem::FloatMatrix::clear(), oofem::Element::computeBoundaryVectorOf(), computeEdgeVolumeAround(), computeEgdeNAt(), computeSurfaceNAt(), computeSurfaceVolumeAround(), oofem::ConvectionBC, oofem::FEMComponent::domain, oofem::FloatArray::dotProduct(), oofem::EdgeLoadBGT, emode, giveApproxOrder(), oofem::GeneralBoundaryCondition::giveBCGeoType(), oofem::Element::giveBoundaryLoadArray(), oofem::Element::giveBoundarySurfaceIntegrationRule(), giveEdgeDofMapping(), oofem::Element::giveElementDofIDMask(), oofem::Element::giveNumberOfDofManagers(), oofem::BoundaryLoad::giveProperty(), giveSurfaceDofMapping(), oofem::GeneralBoundaryCondition::giveType(), HeatTransferEM, oofem::ScalarFunction::isDefined(), oofem::Load::isDofExcluded(), oofem::GeneralBoundaryCondition::isImposed(), Mass1TransferEM, OOFEM_ERROR, oofem::FloatMatrix::plusDyadSymmUpper(), oofem::BoundaryLoad::propertyMultExpr, oofem::FloatMatrix::resize(), oofem::GaussIntegrationRule::SetUpPointsOnLine(), oofem::SurfaceLoadBGT, oofem::FloatMatrix::symmetrized(), and oofem::FloatMatrix::zero().

Referenced by computeBCMtrxAt().

◆ computeBCSubVectorAt()

void oofem::TransportElement::computeBCSubVectorAt	(	FloatArray &	answer,
		TimeStep *	tStep,
		ValueModeType	mode,
		int	indx )

protected

Computes the part of RHS due to applied BCs. The part corresponds to unknown identified by indx param. The result should be localized. The part corresponding to Dirichlet BC is not taken into account.

Parameters

answer	Contains the result.
tStep	Solution step.
mode	Mode of value (incremental, total, ...).
indx	Unknown index.

Definition at line 918 of file transportelement.C.

References oofem::FloatArray::add(), oofem::Element::bodyLoadArray, oofem::Element::boundaryLoadArray, computeBodyBCSubVectorAt(), computeEdgeBCSubVectorAt(), computeSurfaceBCSubVectorAt(), oofem::FEMComponent::domain, oofem::EdgeLoadBGT, oofem::GeneralBoundaryCondition::giveBCGeoType(), oofem::Element::giveBoundaryLoadArray(), oofem::Element::giveNumberOfDofManagers(), oofem::Load::isDofExcluded(), oofem::GeneralBoundaryCondition::isImposed(), OOFEM_ERROR, oofem::FloatArray::resize(), oofem::SurfaceLoadBGT, and oofem::FloatArray::zero().

Referenced by computeBCVectorAt().

◆ computeBCVectorAt()

void oofem::TransportElement::computeBCVectorAt	(	FloatArray &	answer,
		TimeStep *	tStep,
		ValueModeType	mode )

virtual

Computes the RHS contribution to balance equation(s) due to boundary conditions

Definition at line 873 of file transportelement.C.

References assembleLocalContribution(), computeBCSubVectorAt(), oofem::Element::computeNumberOfDofs(), emode, HeatMass1TransferEM, HeatTransferEM, Mass1TransferEM, OOFEM_ERROR, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by computeExternalForcesVector(), and oofem::TransportExternalForceAssembler::vectorFromElement().

◆ computeBmatrixAt()

void oofem::TransportElement::computeBmatrixAt	(	FloatMatrix &	answer,
		const FloatArray &	lcoords )

protectedvirtual

Todo: We should change the transposition in evaldNdx;

Reimplemented in oofem::Lattice2d_mt, and oofem::Lattice3d_mt.

Definition at line 258 of file transportelement.C.

References oofem::FloatMatrix::beTranspositionOf(), emode, oofem::Element::giveInterpolation(), oofem::FloatMatrix::giveNumberOfColumns(), oofem::FloatMatrix::giveNumberOfRows(), HeatMass1TransferEM, HeatTransferEM, Mass1TransferEM, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by computeInternalForcesVector().

◆ computeBodyBCSubVectorAt()

void oofem::TransportElement::computeBodyBCSubVectorAt	(	FloatArray &	answer,
		Load *	load,
		TimeStep *	tStep,
		ValueModeType	mode,
		int	indx )

protected

Definition at line 958 of file transportelement.C.

References oofem::FloatArray::add(), oofem::FloatArray::at(), oofem::Element::computeGlobalCoordinates(), computeNAt(), oofem::Load::computeValueAt(), oofem::Element::computeVolumeAround(), oofem::Load::giveApproxOrder(), oofem::Element::giveIntegrationRule(), oofem::Element::giveInterpolation(), oofem::Element::giveNumberOfDofManagers(), and oofem::FloatArray::resize().

Referenced by computeBCSubVectorAt().

◆ computeBoundaryEdgeLoadVector()

void oofem::TransportElement::computeBoundaryEdgeLoadVector	(	FloatArray &	answer,
		BoundaryLoad *	load,
		int	edge,
		CharType	type,
		ValueModeType	mode,
		TimeStep *	tStep,
		bool	global = true )

overridevirtual

Computes the contribution of the given load at the given boundary edge. In general, the answer should include only relevant DOFs at the edge. The related is giveBoundaryLocationArray method, which should return corresponding code numbers..

Note: Elements which do not have an contribution should resize the vector to be empty.

Parameters

answer	Requested contribution of load (in Global c.s.).
load	Load to compute contribution from.
edge	Edge number.
type	Type of the contribution.
mode	Determines mode of answer.
tStep	Time step when answer is computed.
global	if true (default) then contribution is in global c.s., when false then contribution is in element local c.s.

Todo: Very similar to function computeBoundaryLoadVector and computeEdgeBCSubVectorAt. They should be all merged together.

Todo: Deal with coupled fields (I think they should be another class of problems completely).

Todo: Change the load type to "BoundaryEdgeLoad" maybe?

Reimplemented from oofem::Element.

Definition at line 787 of file transportelement.C.

◆ computeBoundarySurfaceLoadVector()

void oofem::TransportElement::computeBoundarySurfaceLoadVector	(	FloatArray &	answer,
		BoundaryLoad *	load,
		int	boundary,
		CharType	type,
		ValueModeType	mode,
		TimeStep *	tStep,
		bool	global = true )

overridevirtual

Computes the contribution of the given load at the given boundary surface in global coordinate system. In general, the answer should include only relevant DOFs at the edge. The related is giveBoundaryLocationArray method, which should return corresponding code numbers.

Note: Elements which do not have an contribution should resize the vector to be empty.

Parameters

answer	Requested contribution of load.
load	Load to compute contribution from.
boundary	Boundary number.
type	Type of the contribution.
mode	Determines mode of answer.
tStep	Time step when answer is computed.
global	if true (default) then contribution is in global c.s., when false then contribution is in element local c.s.

Todo: Deal with coupled fields (I think they should be another class of problems completely).

Reimplemented from oofem::Element.

Definition at line 607 of file transportelement.C.

◆ computeCapacityMatrix()

void oofem::TransportElement::computeCapacityMatrix	(	FloatMatrix &	answer,
		TimeStep *	tStep )

virtual

Computes the capacity matrix of the receiver

Reimplemented in oofem::Lattice2d_mt, and oofem::Lattice3d_mt.

Definition at line 191 of file transportelement.C.

References assembleLocalContribution(), computeCapacitySubMatrix(), oofem::Element::computeNumberOfDofs(), emode, HeatMass1TransferEM, HeatTransferEM, Mass1TransferEM, OOFEM_ERROR, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by computeInertiaForcesVector(), computeLumpedCapacityVector(), and giveCharacteristicMatrix().

◆ computeCapacitySubMatrix()

void oofem::TransportElement::computeCapacitySubMatrix	(	FloatMatrix &	answer,
		MatResponseMode	rmode,
		int	iri,
		TimeStep *	tStep )

protectedvirtual

Computes the matrix \( P=\int_\Omega N^{\mathrm{T}} N\;\mathrm{d}\Omega \). The result should be multiplied by corresponding coefficient and localized.

Parameters

answer	Contains the result.
rmode	Determines the capacity coefficient to be used.
iri	Index of integration rule to use.
tStep	Time step.

Definition at line 338 of file transportelement.C.

References oofem::FloatMatrix::clear(), computeNAt(), oofem::Element::computeVolumeAround(), oofem::TransportMaterial::giveCharacteristicValue(), giveMaterial(), oofem::Element::integrationRulesArray, oofem::FloatMatrix::plusDyadSymmUpper(), and oofem::FloatMatrix::symmetrized().

Referenced by computeCapacityMatrix().

◆ computeConductivityMatrix()

void oofem::TransportElement::computeConductivityMatrix	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		TimeStep *	tStep )

virtual

Computes the conductivity matrix of the receiver

Reimplemented in oofem::Lattice2d_mt, and oofem::Lattice3d_mt.

Definition at line 215 of file transportelement.C.

References assembleLocalContribution(), computeConductivitySubMatrix(), oofem::Element::computeNumberOfDofs(), emode, HeatMass1TransferEM, HeatTransferEM, Mass1TransferEM, OOFEM_ERROR, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by giveCharacteristicMatrix().

◆ computeConductivitySubMatrix()

void oofem::TransportElement::computeConductivitySubMatrix	(	FloatMatrix &	answer,
		int	iri,
		MatResponseMode	rmode,
		TimeStep *	tStep )

protectedvirtual

Computes the matrix \( C=\int_\Omega \sum_i \left( \frac{\partial N}{\partial x_i}\right)^{\mathrm{T}} \frac{\partial N}{\partial x_i} \;\mathrm{d}\Omega \). The result should be multiplied by corresponding coefficient and localized.

Todo: Is this expression really correct?

Parameters

answer	Contains the result.
iri	Index of integration rule to use.
rmode	Determines the constitutive sub matrix to be used.
tStep	Solution step.

Definition at line 356 of file transportelement.C.

References oofem::FloatMatrix::beProductOf(), computeConstitutiveMatrixAt(), computeGradientMatrixAt(), oofem::Element::computeVolumeAround(), oofem::Element::giveNumberOfDofManagers(), oofem::Element::integrationRulesArray, oofem::FloatMatrix::plusProductSymmUpper(), oofem::FloatMatrix::resize(), oofem::FloatMatrix::symmetrized(), and oofem::FloatMatrix::zero().

Referenced by computeConductivityMatrix().

◆ computeConstitutiveMatrixAt()

void oofem::TransportElement::computeConstitutiveMatrixAt	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		GaussPoint *	gp,
		TimeStep *	tStep )

virtual

Computes constitutive matrix of receiver.

Parameters

answer	Computed answer.
rMode	Material response mode of answer.
gp	Integration point for which constitutive matrix is computed.
tStep	Time step.

Definition at line 480 of file transportelement.C.

References giveCharacteristicMatrix(), and giveMaterial().

Referenced by computeConductivitySubMatrix(), computeFlow(), and oofem::TransportGradientNeumann::computeTangent().

◆ computeEdgeBCSubVectorAt()

void oofem::TransportElement::computeEdgeBCSubVectorAt	(	FloatArray &	answer,
		Load *	load,
		int	iEdge,
		TimeStep *	tStep,
		ValueModeType	mode,
		int	indx )

protected

Computes the part of RHS due to applied BCs on particular edge. The part corresponds to unknown identified by indx param. The result should be localized. The part corresponding to Dirichlet BC is not taken into account.

Parameters

answer	Contains the result.
load	Edge load object.
iEdge	Element edge number subjected to bc.
tStep	Solution step.
mode	Mode of value (incremental, total, ...).
indx	Unknown index.

Todo: Deal with coupled fields

Definition at line 976 of file transportelement.C.

References oofem::FloatArray::add(), oofem::FloatArray::assemble(), oofem::FloatArray::at(), oofem::Element::computeBoundaryVectorOf(), computeEdgeVolumeAround(), computeEgdeNAt(), oofem::BoundaryLoad::computeValueAt(), oofem::ConvectionBC, oofem::FEMComponent::domain, oofem::FloatArray::dotProduct(), emode, oofem::BoundaryLoad::giveApproxOrder(), giveApproxOrder(), giveEdgeDofMapping(), oofem::Element::giveElementDofIDMask(), oofem::Load::giveFormulationType(), oofem::Element::giveInterpolation(), oofem::Element::giveNumberOfDofManagers(), oofem::BoundaryLoad::giveProperty(), oofem::GeneralBoundaryCondition::giveType(), HeatTransferEM, oofem::ScalarFunction::isDefined(), Mass1TransferEM, OOFEM_ERROR, oofem::BoundaryLoad::propertyMultExpr, oofem::RadiationBC, oofem::FloatArray::resize(), oofem::GaussIntegrationRule::SetUpPointsOnLine(), oofem::TransmissionBC, and oofem::FloatArray::zero().

Referenced by computeBCSubVectorAt().

◆ computeEdgeVolumeAround()

virtual double oofem::TransportElement::computeEdgeVolumeAround	(	GaussPoint *	gp,
		int	iEdge )

inlineprotectedvirtual

Computes the length around a integration point on a edge.

Reimplemented in oofem::Brick1_ht, oofem::Lattice2d_mt, oofem::Lattice3d_mt, oofem::QBrick1_ht, oofem::QQuad1_ht, oofem::Quad1_ht, oofem::QuadAxisym1_ht, oofem::QWedge_ht, oofem::Tetrah1_ht, oofem::Tr1_ht, oofem::Tr1Darcy, oofem::TrAxisym1_ht, and oofem::Wedge_ht.

Definition at line 262 of file transportelement.h.

Referenced by computeBCSubMtrxAt(), computeBoundaryEdgeLoadVector(), and computeEdgeBCSubVectorAt().

◆ computeEgdeNAt()

void oofem::TransportElement::computeEgdeNAt	(	FloatArray &	answer,
		int	iEdge,
		const FloatArray &	lcoord )

protectedvirtual

Computes the basis functions at the edge for one unknown.

Definition at line 291 of file transportelement.C.

References oofem::Element::giveInterpolation().

Referenced by computeBCSubMtrxAt(), and computeEdgeBCSubVectorAt().

◆ computeExternalForcesVector()

void oofem::TransportElement::computeExternalForcesVector	(	FloatArray &	answer,
		TimeStep *	tStep,
		ValueModeType	mode )

virtual

Reimplemented in oofem::Tr1Darcy.

Definition at line 867 of file transportelement.C.

References computeBCVectorAt().

◆ computeField()

void oofem::TransportElement::computeField	(	ValueModeType	mode,
		TimeStep *	tStep,
		const FloatArray &	lcoords,
		FloatArray &	answer )

overridevirtual

Computes the unknown vector interpolated at the specified local coordinates. Used for exporting data and mapping fields.

See also: giveElementDofIDMask The unknown vector should match the element field as specified by the element dof IDs.

Parameters

mode	Mode (total, increment, etc) of the output
tStep	Time step to evaluate at
lcoords	Local coordinates to evaluate at
answer	Results

Reimplemented from oofem::Element.

Definition at line 1265 of file transportelement.C.

References oofem::FloatArray::beProductOf(), computeNmatrixAt(), and oofem::Element::computeVectorOf().

◆ computeFlow()

void oofem::TransportElement::computeFlow	(	FloatArray &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep )

virtual

Computes a flow vector in an integration point.

Parameters

answer	Flow vector.
gp	Integration point.
tStep	Time step.

Reimplemented in oofem::Lattice3d_mt.

Definition at line 1275 of file transportelement.C.

References oofem::FloatArray::beProductOf(), computeConstitutiveMatrixAt(), computeGradientMatrixAt(), oofem::Element::computeVectorOf(), emode, oofem::Element::giveElementDofIDMask(), oofem::GaussPoint::giveNaturalCoordinates(), oofem::FloatMatrix::giveNumberOfColumns(), oofem::FloatMatrix::giveNumberOfRows(), oofem::FloatArray::giveSize(), HeatMass1TransferEM, HeatTransferEM, oofem::Element::isActivated(), Mass1TransferEM, oofem::FloatArray::negated(), OOFEM_ERROR, oofem::FloatArray::resize(), oofem::FloatMatrix::resize(), oofem::FloatMatrix::setSubMatrix(), and oofem::FloatArray::zero().

Referenced by oofem::TransportMaterial::giveIPValue().

◆ computeGradientMatrixAt()

void oofem::TransportElement::computeGradientMatrixAt	(	FloatMatrix &	answer,
		const FloatArray &	lcoords )

protectedvirtual

Computes the gradient matrix corresponding to one unknown.

Todo: We should change the transposition in evaldNdx;

Reimplemented in oofem::Lattice2d_mt, and oofem::Lattice3d_mt.

Definition at line 281 of file transportelement.C.

References oofem::FloatMatrix::beTranspositionOf(), and oofem::Element::giveInterpolation().

Referenced by computeConductivitySubMatrix(), computeFlow(), and updateInternalState().

◆ computeInertiaForcesVector()

void oofem::TransportElement::computeInertiaForcesVector	(	FloatArray &	answer,
		TimeStep *	tStep )

virtual

Definition at line 533 of file transportelement.C.

References oofem::FloatArray::beProductOf(), computeCapacityMatrix(), and oofem::Element::computeVectorOf().

Referenced by giveCharacteristicVector().

◆ computeInternalForcesVector()

void oofem::TransportElement::computeInternalForcesVector	(	FloatArray &	answer,
		TimeStep *	tStep )

virtual

Todo: Can/should this part not be built into the flux itself? Probably not? / Mikael

Reimplemented in oofem::Tr1Darcy.

Definition at line 488 of file transportelement.C.

References oofem::FloatArray::beProductOf(), oofem::FloatArray::clear(), computeBmatrixAt(), oofem::TransportMaterial::computeInternalSourceVector(), computeNmatrixAt(), oofem::Element::computeVectorOf(), oofem::Element::computeVolumeAround(), oofem::TransportMaterial::giveFluxVector(), giveMaterial(), oofem::TransportMaterial::hasInternalSource(), oofem::Element::integrationRulesArray, N, and oofem::FloatArray::plusProduct().

Referenced by giveCharacteristicVector().

◆ computeInternalSourceRhsSubVectorAt()

void oofem::TransportElement::computeInternalSourceRhsSubVectorAt	(	FloatArray &	answer,
		TimeStep *	tStep,
		ValueModeType	mode,
		int	indx )

protectedvirtual

Computes the contribution to balance equation(s) due to internal sources

Definition at line 375 of file transportelement.C.

References oofem::FloatArray::add(), oofem::FloatArray::at(), oofem::FloatArray::clear(), oofem::TransportMaterial::computeInternalSourceVector(), computeNAt(), oofem::Element::computeVolumeAround(), oofem::Element::giveDefaultIntegrationRulePtr(), giveMaterial(), and oofem::TransportMaterial::hasInternalSource().

Referenced by computeInternalSourceRhsVectorAt().

◆ computeInternalSourceRhsVectorAt()

void oofem::TransportElement::computeInternalSourceRhsVectorAt	(	FloatArray &	answer,
		TimeStep *	tStep,
		ValueModeType	mode )

virtual

Computes the contribution to balance equation(s) due to internal sources

Reimplemented in oofem::Lattice2d_mt, and oofem::Lattice3d_mt.

Definition at line 402 of file transportelement.C.

References assembleLocalContribution(), computeInternalSourceRhsSubVectorAt(), emode, oofem::FloatArray::giveSize(), HeatMass1TransferEM, HeatTransferEM, oofem::FloatArray::isEmpty(), oofem::FloatArray::isNotEmpty(), Mass1TransferEM, OOFEM_ERROR, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by oofem::TransportExternalForceAssembler::vectorFromElement().

◆ computeIntSourceLHSMatrix()

void oofem::TransportElement::computeIntSourceLHSMatrix	(	FloatMatrix &	answer,
		TimeStep *	tStep )

virtual

Computes the LHS contribution to balance equation(s) due to material internal source

Definition at line 427 of file transportelement.C.

References assembleLocalContribution(), oofem::FloatMatrix::clear(), computeIntSourceLHSSubMatrix(), oofem::Element::computeNumberOfDofs(), emode, giveMaterial(), oofem::TransportMaterial::hasInternalSource(), HeatMass1TransferEM, HeatTransferEM, Mass1TransferEM, OOFEM_ERROR, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

◆ computeIntSourceLHSSubMatrix()

void oofem::TransportElement::computeIntSourceLHSSubMatrix	(	FloatMatrix &	answer,
		MatResponseMode	rmode,
		int	iri,
		TimeStep *	tStep )

virtual

Computes the part of internal source LHS contribution corresponding to unknown identified by rmode parameter

Definition at line 456 of file transportelement.C.

References oofem::FloatMatrix::clear(), computeNAt(), oofem::Element::computeVolumeAround(), oofem::TransportMaterial::giveCharacteristicValue(), giveMaterial(), oofem::Element::integrationRulesArray, oofem::FloatMatrix::plusDyadSymmUpper(), and oofem::FloatMatrix::symmetrized().

Referenced by computeIntSourceLHSMatrix().

◆ computeLoadVector()

void oofem::TransportElement::computeLoadVector	(	FloatArray &	answer,
		BodyLoad *	load,
		CharType	type,
		ValueModeType	mode,
		TimeStep *	tStep )

overridevirtual

Computes the contribution of the given body load (volumetric).

Parameters

answer	Requested contribution of load.
load	Load to compute contribution from.
type	Type of the contribution.
mode	Determines mode of answer.
tStep	Time step when answer is computed.

Todo: We could technically support other types of volume loads. Even a bulk load is a type of transmission b.c.

Todo: Deal with coupled fields (I think they should be another class of problems completely).

Reimplemented from oofem::Element.

Definition at line 561 of file transportelement.C.

References oofem::FloatArray::clear(), oofem::Load::computeValueAt(), oofem::Element::computeVectorOf(), oofem::Element::computeVolumeAround(), oofem::ConvectionBC, emode, oofem::FEInterpolation::evalN(), oofem::Load::giveApproxOrder(), oofem::Element::giveElementDofIDMask(), oofem::FEInterpolation::giveIntegrationRule(), oofem::Element::giveInterpolation(), oofem::FEInterpolation::giveInterpolationOrder(), oofem::GeneralBoundaryCondition::giveType(), HeatMass1TransferEM, N, oofem::FloatArray::plusProduct(), oofem::RadiationBC, and oofem::TransmissionBC.

◆ computeLumpedCapacityVector()

void oofem::TransportElement::computeLumpedCapacityVector	(	FloatArray &	answer,
		TimeStep *	tStep )

virtual

Definition at line 544 of file transportelement.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), computeCapacityMatrix(), oofem::FloatMatrix::giveNumberOfRows(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by giveCharacteristicVector().

◆ computeNAt()

void oofem::TransportElement::computeNAt	(	FloatArray &	answer,
		const FloatArray &	lcoord )

protectedvirtual

Computes the basis functions.

Parameters

answer	The basis functions evaluated at lcoord.
lcoord	The local coordinate.

Definition at line 238 of file transportelement.C.

References oofem::Element::giveInterpolation().

Referenced by computeBodyBCSubVectorAt(), computeCapacitySubMatrix(), computeInternalSourceRhsSubVectorAt(), computeIntSourceLHSSubMatrix(), and computeNmatrixAt().

◆ computeNmatrixAt()

void oofem::TransportElement::computeNmatrixAt	(	FloatMatrix &	answer,
		const FloatArray &	lcoords )

protectedvirtual

Computes the interpolation matrix corresponding to all unknowns. In the default implementation the same approximation order is assumed, but it can be extended.

Reimplemented in oofem::Lattice2d_mt, and oofem::Lattice3d_mt.

Definition at line 245 of file transportelement.C.

References oofem::FloatMatrix::beNMatrixOf(), computeNAt(), emode, HeatTransferEM, and Mass1TransferEM.

Referenced by computeField(), computeInternalForcesVector(), EIPrimaryFieldI_evaluateFieldVectorAt(), and updateInternalState().

◆ computeSurfaceBCSubVectorAt()

void oofem::TransportElement::computeSurfaceBCSubVectorAt	(	FloatArray &	answer,
		Load *	load,
		int	iSurf,
		TimeStep *	tStep,
		ValueModeType	mode,
		int	indx )

protected

Computes the part of RHS due to applied BCs on particular surface. The part corresponds to unknown identified by indx param. The result should be localized. The part corresponding to Dirichlet BC is not taken into account.

Parameters

answer	Contains the result.
load	Surface load object.
iSurf	Element surface number subjected to bc.
tStep	Solution step.
mode	Mode of value (incremental, total, ...).
indx	Unknown index.

Todo: Deal with coupled fields

Definition at line 1049 of file transportelement.C.

References oofem::FloatArray::add(), oofem::FloatArray::assemble(), oofem::FloatArray::at(), oofem::Element::computeBoundaryVectorOf(), computeSurfaceNAt(), computeSurfaceVolumeAround(), oofem::BoundaryLoad::computeValueAt(), oofem::ConvectionBC, oofem::FEMComponent::domain, oofem::FloatArray::dotProduct(), oofem::Element_SurfaceLoadSupport, emode, oofem::BoundaryLoad::giveApproxOrder(), oofem::Element::giveBoundarySurfaceIntegrationRule(), oofem::Element::giveElementDofIDMask(), oofem::Load::giveFormulationType(), oofem::Element::giveInterpolation(), oofem::Element::giveNumberOfDofManagers(), oofem::BoundaryLoad::giveProperty(), giveSurfaceDofMapping(), oofem::GeneralBoundaryCondition::giveType(), HeatTransferEM, oofem::ScalarFunction::isDefined(), Mass1TransferEM, OOFEM_ERROR, oofem::BoundaryLoad::propertyMultExpr, oofem::RadiationBC, oofem::FloatArray::resize(), oofem::Element::testElementExtension(), oofem::TransmissionBC, and oofem::FloatArray::zero().

Referenced by computeBCSubVectorAt().

◆ computeSurfaceNAt()

void oofem::TransportElement::computeSurfaceNAt	(	FloatArray &	answer,
		int	iSurf,
		const FloatArray &	lcoord )

protectedvirtual

Definition at line 314 of file transportelement.C.

References oofem::Element::giveInterpolation().

Referenced by computeBCSubMtrxAt(), and computeSurfaceBCSubVectorAt().

◆ computeSurfaceVolumeAround()

virtual double oofem::TransportElement::computeSurfaceVolumeAround	(	GaussPoint *	gp,
		int	iSurf )

inlineprotectedvirtual

Reimplemented in oofem::Brick1_ht, oofem::QBrick1_ht, and oofem::Tetrah1_ht.

Definition at line 265 of file transportelement.h.

Referenced by computeBCSubMtrxAt(), and computeSurfaceBCSubVectorAt().

◆ computeTangentFromEdgeLoad()

void oofem::TransportElement::computeTangentFromEdgeLoad	(	FloatMatrix &	answer,
		BoundaryLoad *	load,
		int	boundary,
		MatResponseMode	rmode,
		TimeStep *	tStep )

overridevirtual

Computes the tangent contribution of the given load at the given boundary.

Note: Elements which do not have an contribution should resize the vector to be empty.

Parameters

answer	Requested contribution of load.
load	Load to compute contribution from.
boundary	Surface number.
rmode	Mode of the contribution.
tStep	Time step when answer is computed.

Todo: Deal with coupled fields (I think they should be another class of problems completely).

Reimplemented from oofem::Element.

Definition at line 733 of file transportelement.C.

References oofem::FloatArray::at(), oofem::FloatArray::beProductOf(), oofem::FEInterpolation::boundaryEdgeGiveNodes(), oofem::FEInterpolation::boundaryEvalN(), oofem::FEInterpolation::boundaryGiveTransformationJacobian(), oofem::FEInterpolation::boundaryLocal2Global(), oofem::FloatMatrix::clear(), oofem::Element::computeBoundaryVectorOf(), oofem::ConvectionBC, emode, oofem::BoundaryLoad::giveApproxOrder(), oofem::Element::giveBoundaryEdgeIntegrationRule(), oofem::Element::giveElementDofIDMask(), oofem::Element::giveGeometryType(), oofem::Element::giveInterpolation(), oofem::FEInterpolation::giveInterpolationOrder(), oofem::BoundaryLoad::giveProperty(), giveThicknessAt(), oofem::BoundaryLoad::giveType(), HeatMass1TransferEM, oofem::ScalarFunction::isDefined(), N, OOFEM_ERROR, oofem::FloatMatrix::plusProductSymmUpper(), oofem::BoundaryLoad::propertyMultExpr, oofem::FloatArray::subtract(), and oofem::FloatMatrix::symmetrized().

◆ computeTangentFromSurfaceLoad()

void oofem::TransportElement::computeTangentFromSurfaceLoad	(	FloatMatrix &	answer,
		BoundaryLoad *	load,
		int	boundary,
		MatResponseMode	rmode,
		TimeStep *	tStep )

overridevirtual

Computes the tangent contribution of the given load at the given boundary.

Note: Elements which do not have an contribution should resize the vector to be empty.

Parameters

answer	Requested contribution of load.
load	Load to compute contribution from.
boundary	Surface number.
rmode	Mode of the contribution.
tStep	Time step when answer is computed.

Todo: Deal with coupled fields (I think they should be another class of problems completely).

Reimplemented from oofem::Element.

Definition at line 682 of file transportelement.C.

References oofem::FEInterpolation::boundaryEvalN(), oofem::FEInterpolation::boundaryGiveNodes(), oofem::FEInterpolation::boundaryGiveTransformationJacobian(), oofem::FEInterpolation::boundaryLocal2Global(), oofem::FloatMatrix::clear(), oofem::Element::computeBoundaryVectorOf(), oofem::ConvectionBC, oofem::FloatArray::dotProduct(), emode, oofem::BoundaryLoad::giveApproxOrder(), oofem::Element::giveBoundarySurfaceIntegrationRule(), oofem::Element::giveElementDofIDMask(), oofem::Element::giveGeometryType(), oofem::Element::giveInterpolation(), oofem::FEInterpolation::giveInterpolationOrder(), oofem::BoundaryLoad::giveProperty(), giveThicknessAt(), oofem::BoundaryLoad::giveType(), HeatMass1TransferEM, oofem::ScalarFunction::isDefined(), N, OOFEM_ERROR, oofem::FloatMatrix::plusProductSymmUpper(), oofem::BoundaryLoad::propertyMultExpr, and oofem::FloatMatrix::symmetrized().

◆ computeVof()

double oofem::TransportElement::computeVof ( TimeStep * tStep )

virtual

Computes the relative volume of reference material in element. Default implementation uses user provided function.

Parameters

tStep Time step

Definition at line 1431 of file transportelement.C.

References oofem::Function::evaluateAtTime(), oofem::FEMComponent::giveDomain(), oofem::Domain::giveFunction(), oofem::TimeStep::giveTargetTime(), and vofFunction.

Referenced by oofem::TwoPhaseMaterial::giveVof().

◆ EIPrimaryFieldI_evaluateFieldVectorAt()

int oofem::TransportElement::EIPrimaryFieldI_evaluateFieldVectorAt	(	FloatArray &	answer,
		PrimaryField &	pf,
		const FloatArray &	coords,
		IntArray &	dofId,
		ValueModeType	mode,
		TimeStep *	tStep )

overridevirtual

Evaluates the value of field at given point of interest (should be located inside receiver's volume) using element interpolation.

Todo

This should use local coordinates instead of having all elements search for it manually.

Shouldn't this just be part of Element? It's very much the core of functionality for elements.

Parameters

answer	Field evaluated at coordinate.
pf	Field to use for evaluation.
coords	Coordinate.
dofId	IDs of DOFs to evaluate.
mode	Mode of field.
tStep	Time step to evaluate at.

Returns: Zero if ok, nonzero when error encountered.

Implements oofem::EIPrimaryFieldInterface.

Definition at line 1372 of file transportelement.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::IntArray::at(), oofem::Element::computeLocalCoordinates(), computeNmatrixAt(), oofem::Element::computeVectorOf(), oofem::IntArray::findFirstIndexOf(), oofem::Element::giveElementDofIDMask(), oofem::FloatArray::giveSize(), oofem::IntArray::giveSize(), OOFEM_ERROR, oofem::FloatArray::resize(), oofem::sum(), and oofem::FloatArray::zero().

◆ giveApproxOrder()

int oofem::TransportElement::giveApproxOrder ( int unknownIndx )

protectedvirtual

Reimplemented in oofem::Lattice2d_mt, oofem::Lattice3d_mt, oofem::QuadAxisym1_ht, and oofem::TrAxisym1_ht.

Definition at line 332 of file transportelement.C.

References oofem::Element::giveInterpolation().

Referenced by computeBCSubMtrxAt(), and computeEdgeBCSubVectorAt().

◆ giveCharacteristicMatrix()

void oofem::TransportElement::giveCharacteristicMatrix	(	FloatMatrix &	answer,
		CharType	type,
		TimeStep *	tStep )

overridevirtual

Computes characteristic matrix of receiver of requested type in given time step.

Parameters

answer	Requested characteristic matrix (stiffness, tangent, ...). If element has no capability to compute requested type of characteristic matrix error function is invoked.
type	Id of characteristic component requested.
tStep	Time step when answer is computed.

Reimplemented from oofem::Element.

Definition at line 112 of file transportelement.C.

References oofem::__CharTypeToString(), oofem::FloatMatrix::add(), oofem::FloatMatrix::at(), computeBCMtrxAt(), computeCapacityMatrix(), computeConductivityMatrix(), oofem::FloatMatrix::giveNumberOfColumns(), oofem::FloatMatrix::giveNumberOfRows(), and OOFEM_ERROR.

Referenced by computeConstitutiveMatrixAt().

◆ giveCharacteristicVector()

void oofem::TransportElement::giveCharacteristicVector	(	FloatArray &	answer,
		CharType	type,
		ValueModeType	mode,
		TimeStep *	tStep )

overridevirtual

Computes characteristic vector of receiver of requested type in given time step. If element has no capability to compute requested type of characteristic vector error function is invoked.

Parameters

answer	Requested characteristic vector.
type	Id of characteristic component requested.
mode	Determines mode of answer.
tStep	Time step when answer is computed.

Reimplemented from oofem::Element.

Definition at line 144 of file transportelement.C.

References oofem::__CharTypeToString(), computeInertiaForcesVector(), computeInternalForcesVector(), computeLumpedCapacityVector(), OOFEM_ERROR, and oofem::FloatArray::resize().

◆ giveClassName()

const char * oofem::TransportElement::giveClassName ( ) const

inlineoverridevirtual

Returns: Class name of the receiver.

Reimplemented from oofem::Element.

Reimplemented in oofem::TrAxisym1_hmt, oofem::TrAxisym1_ht, oofem::TrAxisym1_mt, oofem::Wedge_hmt, oofem::Wedge_ht, and oofem::Wedge_mt.

Definition at line 97 of file transportelement.h.

◆ giveDofManDofIDMask()

void oofem::TransportElement::giveDofManDofIDMask	(	int	inode,
		IntArray &	answer ) const

overridevirtual

Returns dofmanager dof mask for node. This mask defines the dofs which are used by element in node. Mask influences the code number ordering for particular node. Code numbers are ordered according to node order and dofs belonging to particular node are ordered according to this mask. If element requests dofs using node mask which are not in node then error is generated. This masking allows node to be shared by different elements with different dofs in same node. Elements local code numbers are extracted from node using this mask. Must be defined by particular element.

Parameters

inode	Mask is computed for local dofmanager with inode number.
answer	Mask for node.

Reimplemented from oofem::Element.

Definition at line 97 of file transportelement.C.

References emode, HeatMass1TransferEM, HeatTransferEM, Mass1TransferEM, and OOFEM_ERROR.

◆ giveEdgeDofMapping()

void oofem::TransportElement::giveEdgeDofMapping	(	IntArray &	mask,
		int	iEdge )

protectedvirtual

Gives the node indexes for given edge.

Definition at line 303 of file transportelement.C.

References oofem::Element::giveInterpolation().

Referenced by computeBCSubMtrxAt(), and computeEdgeBCSubVectorAt().

◆ giveInternalStateAtNode()

int oofem::TransportElement::giveInternalStateAtNode	(	FloatArray &	answer,
		InternalStateType	type,
		InternalStateMode	mode,
		int	node,
		TimeStep *	tStep )

overridevirtual

Returns internal state variable (like stress,strain) at node of element in Reduced form, the way how is obtained is dependent on InternalValueType. The value may be local, or smoothed using some recovery technique. Returns zero if element is unable to respond to request.

Parameters

answer	Contains result, zero sized if not supported.
type	Determines the internal variable requested (physical meaning).
mode	Determines the mode of variable (recovered, local, ...).
node	Node number, for which variable is required.
tStep	Time step.

Returns: Nonzero if o.k, zero otherwise.

Reimplemented from oofem::Element.

Definition at line 1442 of file transportelement.C.

References oofem::FloatArray::at(), oofem::DofManager::end(), oofem::DofManager::findDofWithDofId(), oofem::Element::giveNode(), and oofem::FloatArray::resize().

Referenced by oofem::Brick1_ht::drawScalar(), oofem::Quad1_ht::drawScalar(), and oofem::Tetrah1_ht::drawScalar().

◆ giveMaterial()

Material * oofem::TransportElement::giveMaterial ( )

overridevirtual

Returns: Reference to the associated material of element.

Reimplemented from oofem::Element.

Definition at line 1421 of file transportelement.C.

References oofem::FEMComponent::domain, oofem::Element::giveCrossSection(), giveMaterial(), and oofem::Element::material.

Referenced by oofem::NonStationaryTransportProblem::applyIC(), checkConsistency(), oofem::Lattice2d_mt::computeCapacityMatrix(), oofem::Lattice3d_mt::computeCapacityMatrix(), computeCapacitySubMatrix(), oofem::Lattice2d_mt::computeConductivityMatrix(), oofem::Lattice3d_mt::computeConductivityMatrix(), computeConstitutiveMatrixAt(), oofem::Lattice3d_mt::computeFlow(), oofem::Tr1Darcy::computeInternalForcesVector(), computeInternalForcesVector(), computeInternalSourceRhsSubVectorAt(), oofem::Lattice2d_mt::computeInternalSourceRhsVectorAt(), oofem::Lattice3d_mt::computeInternalSourceRhsVectorAt(), computeIntSourceLHSMatrix(), computeIntSourceLHSSubMatrix(), oofem::Tr1Darcy::computeStiffnessMatrix(), giveMaterial(), oofem::Lattice2d_mt::updateInternalState(), oofem::Lattice3d_mt::updateInternalState(), updateInternalState(), and oofem::NonStationaryTransportProblem::updateYourself().

◆ giveSurfaceDofMapping()

void oofem::TransportElement::giveSurfaceDofMapping	(	IntArray &	mask,
		int	iSurf )

protectedvirtual

Definition at line 323 of file transportelement.C.

References oofem::Element::giveInterpolation().

Referenced by computeBCSubMtrxAt(), and computeSurfaceBCSubVectorAt().

◆ giveThicknessAt()

virtual double oofem::TransportElement::giveThicknessAt ( const FloatArray & gcoords )

inlinevirtual

Gives the thickness at some global coordinate. For solid elements, the value returned is 1.0 (which is the default implementation).

Parameters

gcoords Global coordinates.

Returns: Thickness of element at given coordinate.

Todo: Move this into the base element?

Reimplemented in oofem::QQuad1_ht, oofem::Quad1_ht, oofem::QuadAxisym1_ht, oofem::Tr1_ht, oofem::Tr1Darcy, and oofem::TrAxisym1_ht.

Definition at line 106 of file transportelement.h.

Referenced by computeBoundarySurfaceLoadVector(), computeTangentFromEdgeLoad(), and computeTangentFromSurfaceLoad().

◆ giveTransportCrossSection()

TransportCrossSection * oofem::TransportElement::giveTransportCrossSection ( )

Definition at line 1414 of file transportelement.C.

References oofem::Element::giveCrossSection().

◆ initializeFrom()

void oofem::TransportElement::initializeFrom	(	InputRecord &	ir,
		int	priority )

overridevirtual

Reimplemented from oofem::Element.

Reimplemented in oofem::Wedge_ht.

Definition at line 83 of file transportelement.C.

References oofem::FEMComponent::giveDomain(), oofem::Element::initializeFrom(), IPK_TransportElement_vof_function, oofem::FEMComponent::number, PM_UPDATE_PARAMETER, and vofFunction.

◆ postInitialize()

void oofem::TransportElement::postInitialize ( )

overridevirtual

Performs post initialization steps.

Reimplemented from oofem::Element.

Definition at line 91 of file transportelement.C.

References oofem::Element::postInitialize().

◆ updateInternalState()

void oofem::TransportElement::updateInternalState ( TimeStep * tStep )

overridevirtual

Updates element state after equilibrium in time step has been reached. Default implementation updates all integration rules defined by integrationRulesArray member variable. Doing this, all integration points and their material statuses are updated also. All temporary history variables, which now describe equilibrium state are copied into equilibrium ones. The existing internal state is used for update.

Parameters

tStep Time step for newly reached state.

See also: Material::updateYourself; IntegrationRule::updateYourself; GaussPoint::updateYourself; Element::updateInternalState

Todo: This method should only used when applying IC (but nltransienttransport and nonstationarytransport uses it wrong, and they should be removed)

Todo: Why is the state vector the unknown solution at the gauss point? / Mikael

Todo: We need to sort out multiple materials for coupled (heat+mass) problems

Reimplemented from oofem::Element.

Definition at line 1327 of file transportelement.C.

References oofem::FloatArray::beProductOf(), computeGradientMatrixAt(), computeNmatrixAt(), oofem::Element::computeVectorOf(), oofem::Element::giveElementDofIDMask(), oofem::TransportMaterial::giveFluxVector(), giveMaterial(), oofem::CemhydMat::initMaterial(), oofem::Element::integrationRulesArray, oofem::TimeStep::isIcApply(), and oofem::TransportMaterial::updateInternalState().

Referenced by oofem::NLTransientTransportProblem::applyIC(), oofem::NonStationaryTransportProblem::applyIC(), oofem::TransientTransportProblem::applyIC(), and oofem::StationaryTransportProblem::solveYourselfAt().

Member Data Documentation

◆ emode

ElementMode oofem::TransportElement::emode

protected

Definition at line 60 of file transportelement.h.

Referenced by oofem::Brick1_hmt::Brick1_hmt(), oofem::Brick1_mt::Brick1_mt(), computeBCMtrxAt(), computeBCSubMtrxAt(), computeBCVectorAt(), computeBmatrixAt(), computeBoundaryEdgeLoadVector(), computeBoundarySurfaceLoadVector(), computeCapacityMatrix(), computeConductivityMatrix(), computeEdgeBCSubVectorAt(), computeFlow(), computeInternalSourceRhsVectorAt(), computeIntSourceLHSMatrix(), computeLoadVector(), computeNmatrixAt(), oofem::Line1_ht::computeNumberOfDofs(), oofem::QBrick1_ht::computeNumberOfDofs(), oofem::Tetrah1_ht::computeNumberOfDofs(), oofem::Tr1_ht::computeNumberOfDofs(), computeSurfaceBCSubVectorAt(), computeTangentFromEdgeLoad(), computeTangentFromSurfaceLoad(), oofem::Quad1_ht::drawScalar(), giveDofManDofIDMask(), oofem::Line1_hmt::Line1_hmt(), oofem::Line1_mt::Line1_mt(), oofem::QBrick1_hmt::QBrick1_hmt(), oofem::QBrick1_mt::QBrick1_mt(), oofem::QQuad1_hmt::QQuad1_hmt(), oofem::QQuad1_mt::QQuad1_mt(), oofem::Quad1_hmt::Quad1_hmt(), oofem::Quad1_mt::Quad1_mt(), oofem::QuadAxisym1_hmt::QuadAxisym1_hmt(), oofem::QuadAxisym1_mt::QuadAxisym1_mt(), oofem::QWedge_hmt::QWedge_hmt(), oofem::QWedge_mt::QWedge_mt(), oofem::Tetrah1_hmt::Tetrah1_hmt(), oofem::Tr1_hmt::Tr1_hmt(), oofem::Tr1_mt::Tr1_mt(), TransportElement(), oofem::TrAxisym1_hmt::TrAxisym1_hmt(), oofem::TrAxisym1_mt::TrAxisym1_mt(), oofem::Wedge_hmt::Wedge_hmt(), and oofem::Wedge_mt::Wedge_mt().

◆ IPK_TransportElement_vof_function

ParamKey oofem::TransportElement::IPK_TransportElement_vof_function

staticprotected

Definition at line 66 of file transportelement.h.

Referenced by initializeFrom().

◆ stefanBoltzmann

const double oofem::TransportElement::stefanBoltzmann = 5.67e-8

staticprotected

Stefan–Boltzmann constant W/m2/K4.

Definition at line 62 of file transportelement.h.

Referenced by computeBoundaryEdgeLoadVector(), and computeBoundarySurfaceLoadVector().

◆ vofFunction

int oofem::TransportElement::vofFunction = 0

protected

Fuction determining the relative volume of reference material in element.

Definition at line 64 of file transportelement.h.

Referenced by computeVof(), initializeFrom(), and TransportElement().

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

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Static Protected Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ ElementMode

Constructor & Destructor Documentation

◆ TransportElement()

Member Function Documentation

◆ assembleLocalContribution() [1/2]

◆ assembleLocalContribution() [2/2]

◆ checkConsistency()

◆ computeBCMtrxAt()

◆ computeBCSubMtrxAt()

◆ computeBCSubVectorAt()

◆ computeBCVectorAt()

◆ computeBmatrixAt()

◆ computeBodyBCSubVectorAt()

◆ computeBoundaryEdgeLoadVector()

◆ computeBoundarySurfaceLoadVector()

◆ computeCapacityMatrix()

◆ computeCapacitySubMatrix()

◆ computeConductivityMatrix()

◆ computeConductivitySubMatrix()

◆ computeConstitutiveMatrixAt()

◆ computeEdgeBCSubVectorAt()

◆ computeEdgeVolumeAround()

◆ computeEgdeNAt()

◆ computeExternalForcesVector()

◆ computeField()

◆ computeFlow()

◆ computeGradientMatrixAt()

◆ computeInertiaForcesVector()

◆ computeInternalForcesVector()

◆ computeInternalSourceRhsSubVectorAt()

◆ computeInternalSourceRhsVectorAt()

◆ computeIntSourceLHSMatrix()

◆ computeIntSourceLHSSubMatrix()

◆ computeLoadVector()

◆ computeLumpedCapacityVector()

◆ computeNAt()

◆ computeNmatrixAt()

◆ computeSurfaceBCSubVectorAt()

◆ computeSurfaceNAt()

◆ computeSurfaceVolumeAround()

◆ computeTangentFromEdgeLoad()

◆ computeTangentFromSurfaceLoad()

◆ computeVof()

◆ EIPrimaryFieldI_evaluateFieldVectorAt()

◆ giveApproxOrder()

◆ giveCharacteristicMatrix()

◆ giveCharacteristicVector()

◆ giveClassName()

◆ giveDofManDofIDMask()

◆ giveEdgeDofMapping()

◆ giveInternalStateAtNode()

◆ giveMaterial()

◆ giveSurfaceDofMapping()

◆ giveThicknessAt()

◆ giveTransportCrossSection()

◆ initializeFrom()

◆ postInitialize()

◆ updateInternalState()

Member Data Documentation

◆ emode

◆ IPK_TransportElement_vof_function

◆ stefanBoltzmann

◆ vofFunction