#include <tr1_2d_supg2.h>

Inheritance diagram for oofem::TR1_2D_SUPG2:

Collaboration diagram for oofem::TR1_2D_SUPG2:

Public Member Functions
	TR1_2D_SUPG2 (int n, Domain *d)
void	computeAccelerationTerm_MB (FloatMatrix &answer, TimeStep *tStep) override
void	computeAdvectionTerm_MB (FloatArray &answer, TimeStep *tStep) override
void	computeAdvectionDerivativeTerm_MB (FloatMatrix &answer, TimeStep *tStep) override
void	computeDiffusionTerm_MB (FloatArray &answer, TimeStep *tStep) override
void	computeDiffusionDerivativeTerm_MB (FloatMatrix &answer, MatResponseMode mode, TimeStep *tStep) override
void	computePressureTerm_MB (FloatMatrix &answer, TimeStep *tStep) override
void	computeLSICStabilizationTerm_MB (FloatMatrix &answer, TimeStep *tStep) override
void	computeLinearAdvectionTerm_MC (FloatMatrix &answer, TimeStep *tStep) override
void	computeAdvectionTerm_MC (FloatArray &answer, TimeStep *tStep) override
void	computeAdvectionDerivativeTerm_MC (FloatMatrix &answer, TimeStep *tStep) override
void	computeDiffusionDerivativeTerm_MC (FloatMatrix &answer, TimeStep *tStep) override
void	computeDiffusionTerm_MC (FloatArray &answer, TimeStep *tStep) override
void	computeAccelerationTerm_MC (FloatMatrix &answer, TimeStep *tStep) override
void	computePressureTerm_MC (FloatMatrix &answer, TimeStep *tStep) override
void	computeBCRhsTerm_MB (FloatArray &answer, TimeStep *tStep) override
void	computeBCRhsTerm_MC (FloatArray &answer, TimeStep *tStep) override
void	updateStabilizationCoeffs (TimeStep *tStep) override
void	updateElementForNewInterfacePosition (TimeStep *tStep) override
double	computeCriticalTimeStep (TimeStep *tStep) override
	Computes the critical time increment.
const char *	giveClassName () const override
const char *	giveInputRecordName () const override
void	giveDofManDofIDMask (int inode, IntArray &answer) const override
int	computeNumberOfDofs () override
void	initializeFrom (InputRecord &ir, int priority) override
void	giveInputRecord (DynamicInputRecord &input) override
void	updateYourself (TimeStep *tStep) override
void	saveContext (DataStream &stream, ContextMode mode) override
void	restoreContext (DataStream &stream, ContextMode mode) override
Interface *	giveInterface (InterfaceType) override
int	EIPrimaryFieldI_evaluateFieldVectorAt (FloatArray &answer, PrimaryField &pf, const FloatArray &coords, IntArray &dofId, ValueModeType mode, TimeStep *tStep) override
double	computeLEPLICVolumeFraction (const FloatArray &n, const double p, LEPlic *matInterface, bool updFlag) override
	Computes corresponding volume fraction to given interface position.
void	formMaterialVolumePoly (Polygon &matvolpoly, LEPlic *matInterface, const FloatArray &normal, const double p, bool updFlag) override
	Assembles the true element material polygon (takes receiver vof into accout).
void	formVolumeInterfacePoly (Polygon &matvolpoly, LEPlic *matInterface, const FloatArray &normal, const double p, bool updFlag) override
	Assembles receiver material polygon based solely on given interface line.
double	truncateMatVolume (const Polygon &matvolpoly, double &volume) override
	Truncates given material polygon to receiver.
void	giveElementCenter (LEPlic *mat_interface, FloatArray &center, bool updFlag) override
	Computes the receiver center (in updated Lagrangian configuration).
void	formMyVolumePoly (Polygon &myPoly, LEPlic *mat_interface, bool updFlag) override
	Assembles receiver volume.
Element *	giveElement () override
	Return number of receiver's element.
double	computeMyVolume (LEPlic *matInterface, bool updFlag) override
	Computes the volume of receiver.
void	NodalAveragingRecoveryMI_computeNodalValue (FloatArray &answer, int node, InternalStateType type, TimeStep *tStep) override
void	SPRNodalRecoveryMI_giveSPRAssemblyPoints (IntArray &pap) override
void	SPRNodalRecoveryMI_giveDofMansDeterminedByPatch (IntArray &answer, int pap) override
int	SPRNodalRecoveryMI_giveNumberOfIP () override
SPRPatchType	SPRNodalRecoveryMI_givePatchType () override
int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep) override
int	giveDefaultIntegrationRule () const override
IntegrationRule *	giveDefaultIntegrationRulePtr () override
int	giveInternalStateAtNode (FloatArray &answer, InternalStateType type, InternalStateMode mode, int node, TimeStep *tStep) override
void	drawRawGeometry (oofegGraphicContext &gc, TimeStep *tStep) override
void	drawScalar (oofegGraphicContext &gc, TimeStep *tStep) override
void	printOutputAt (FILE file, TimeStep tStep) override
Public Member Functions inherited from oofem::TR1_2D_SUPG
	TR1_2D_SUPG (int n, Domain *d)
FEInterpolation *	giveInterpolation () const override
void	computeLoadVector (FloatArray &answer, BodyLoad load, CharType type, ValueModeType mode, TimeStep tStep) override
void	computeSlipWithFrictionBCTerm_MB (FloatMatrix &answer, Load load, int side, TimeStep tStep) override
void	computePenetrationWithResistanceBCTerm_MB (FloatMatrix &answer, Load load, int side, TimeStep tStep) override
void	computeOutFlowBCTerm_MB (FloatMatrix &answer, int side, TimeStep *tStep) override
void	computeHomogenizedReinforceTerm_MB (FloatMatrix &answer, Load load, TimeStep tStep) override
void	computeHomogenizedReinforceTerm_MC (FloatMatrix &answer, Load load, TimeStep tStep) override
MaterialMode	giveMaterialMode () override
Element_Geometry_Type	giveGeometryType () const override
int	checkConsistency () override
	Used to check consistency and initialize some element geometry data (area,b,c).
double	computeVolumeAround (GaussPoint *gp) override
double	computeCriticalLEPlicTimeStep (TimeStep *tStep) override
	Computes critical time step.
double	LS_PCS_computeF (LevelSetPCS ls, TimeStep tStep) override
void	LS_PCS_computedN (FloatMatrix &answer) override
double	LS_PCS_computeVolume () override
	Returns receiver's volume.
double	LS_PCS_computeS (LevelSetPCS ls, TimeStep tStep) override
void	LS_PCS_computeVOFFractions (FloatArray &answer, FloatArray &fi) override
Public Member Functions inherited from oofem::SUPGElement
	SUPGElement (int n, Domain *aDomain)
void	giveCharacteristicMatrix (FloatMatrix &answer, CharType type, TimeStep *tStep) override
void	giveCharacteristicVector (FloatArray &answer, CharType type, ValueModeType mode, TimeStep *tStep) override
virtual void	computeBCLhsPressureTerm_MC (FloatMatrix &answer, TimeStep *tStep)
virtual void	computeBCLhsTerm_MB (FloatMatrix &answer, TimeStep *tStep)
virtual void	computeBCLhsPressureTerm_MB (FloatMatrix &answer, TimeStep *tStep)
void	updateInternalState (TimeStep *tStep) override
Public Member Functions inherited from oofem::FMElement
	FMElement (int n, Domain *aDomain)
void	computeVectorOfVelocities (ValueModeType mode, TimeStep *tStep, FloatArray &velocities)
void	computeVectorOfPressures (ValueModeType mode, TimeStep *tStep, FloatArray &pressures)
FloatArray	computeVectorOfVelocities (ValueModeType mode, TimeStep *tStep)
FloatArray	computeVectorOfPressures (ValueModeType mode, 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)
virtual void	computeBoundarySurfaceLoadVector (FloatArray &answer, BoundaryLoad load, int boundary, CharType type, ValueModeType mode, TimeStep tStep, bool global=true)
virtual void	computeTangentFromSurfaceLoad (FloatMatrix &answer, BoundaryLoad load, int boundary, MatResponseMode rmode, TimeStep tStep)
virtual void	computeTangentFromEdgeLoad (FloatMatrix &answer, BoundaryLoad load, int boundary, MatResponseMode rmode, TimeStep tStep)
virtual void	computeBoundaryEdgeLoadVector (FloatArray &answer, BoundaryLoad load, int edge, CharType type, ValueModeType mode, TimeStep tStep, bool global=true)
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	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 void	computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer)
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 (DofIDItem id) const
virtual const FEInterpolation *	getGeometryInterpolation () const
virtual Material *	giveMaterial ()
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 int	giveIntegrationRuleLocalCodeNumbers (IntArray &answer, IntegrationRule &ie)
int	giveRegionNumber ()
virtual void	initializeYourself (TimeStep *timeStepWhenICApply)
virtual bool	isActivated (TimeStep *tStep)
virtual bool	isCast (TimeStep *tStep)
virtual void	initForNewStep ()
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
int	giveNumberOfIntegrationRules ()
virtual IntegrationRule *	giveIntegrationRule (int i)
std::vector< std ::unique_ptr< IntegrationRule > > &	giveIntegrationRulesArray ()
virtual int	testElementExtension (ElementExtension ext)
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	drawDeformedGeometry (oofegGraphicContext &gc, TimeStep *tStep, UnknownType)
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
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.
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.
std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros).
Public Member Functions inherited from oofem::SpatialLocalizerInterface
	SpatialLocalizerInterface (Element *element)
virtual int	SpatialLocalizerI_containsPoint (const FloatArray &coords)
int	SpatialLocalizerI_BBoxContainsPoint (const FloatArray &coords)
virtual void	SpatialLocalizerI_giveBBox (FloatArray &bb0, FloatArray &bb1)
virtual double	SpatialLocalizerI_giveClosestPoint (FloatArray &lcoords, FloatArray &closest, const FloatArray &gcoords)
Public Member Functions inherited from oofem::Interface
	Interface ()
	Constructor.
virtual	~Interface ()
Public Member Functions inherited from oofem::EIPrimaryFieldInterface
	EIPrimaryFieldInterface ()
Public Member Functions inherited from oofem::ZZNodalRecoveryModelInterface
	ZZNodalRecoveryModelInterface (Element *element)
	Constructor.
virtual bool	ZZNodalRecoveryMI_computeNValProduct (FloatMatrix &answer, InternalStateType type, TimeStep *tStep)
virtual void	ZZNodalRecoveryMI_computeNNMatrix (FloatArray &answer, InternalStateType type)
Public Member Functions inherited from oofem::NodalAveragingRecoveryModelInterface
	NodalAveragingRecoveryModelInterface ()
	Constructor.
Public Member Functions inherited from oofem::SPRNodalRecoveryModelInterface
	SPRNodalRecoveryModelInterface ()
	Constructor.
Public Member Functions inherited from oofem::LEPlicElementInterface
	LEPlicElementInterface ()
bool	isBoundary ()
	Returns true if cell is boundary.
void	setTempLineConstant (double tp)
void	setTempInterfaceNormal (FloatArray tg)
void	setTempVolumeFraction (double v)
void	setPermanentVolumeFraction (double v)
void	addTempVolumeFraction (double v)
double	giveVolumeFraction ()
double	giveTempVolumeFraction ()
void	giveTempInterfaceNormal (FloatArray &n)
double	giveTempLineConstant ()
void	updateYourself (TimeStep *tStep)
void	saveContext (DataStream &stream, ContextMode mode)
void	restoreContext (DataStream &stream, ContextMode mode)
Public Member Functions inherited from oofem::LevelSetPCSElementInterface
	LevelSetPCSElementInterface ()

Protected Member Functions
void	computeGaussPoints () override
void	postInitialize () override
	Performs post initialization steps.
void	computeDeviatoricStress (FloatArray &answer, const FloatArray &eps, GaussPoint gp, TimeStep tStep) override
void	computeTangent (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep) override
void	computeNVector (FloatArray &answer, GaussPoint *gp)
void	updateVolumePolygons (Polygon &referenceFluidPoly, Polygon &secondFluidPoly, int &rfPoints, int &sfPoints, const FloatArray &normal, const double p, bool updFlag)
double	computeVolumeAroundID (GaussPoint *gp, integrationDomain id, const std::vector< FloatArray > &idpoly)
void	updateIntegrationRules ()
Material *	_giveMaterial (int indx)
Protected Member Functions inherited from oofem::TR1_2D_SUPG
void	giveLocalVelocityDofMap (IntArray &map) override
void	giveLocalPressureDofMap (IntArray &map) override
void	computeNMtrx (FloatArray &answer, GaussPoint *gp)
void	computeDeviatoricStrain (FloatArray &answer, GaussPoint gp, TimeStep tStep) override
virtual void	initGeometry ()

Protected Attributes
Polygon	myPoly [2]
std::vector< FloatArray >	vcoords [2]
integrationDomain	id [2]
int	mat [2]
std ::unique_ptr< IntegrationRule >	defaultIRule
Protected Attributes inherited from oofem::TR1_2D_SUPG
double	b [3]
double	c [3]
double	area = 0.
Protected Attributes inherited from oofem::SUPGElement
IntArray	boundarySides
	Array of boundary sides.
IntArray	boundaryCodes
	Boundary sides codes.
double	t_supg = 0.
double	t_pspg = 0.
double	t_lsic = 0.
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.
Protected Attributes inherited from oofem::LEPlicElementInterface
bool	permanentVofFlag
double	vof
	Volume fraction of reference fluid in element.
double	temp_vof
double	p
	Line constant of line segment representing interface.
double	temp_p
FloatArray	normal
	Interface segment normal.
FloatArray	temp_normal

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
Static Protected Attributes inherited from oofem::TR1_2D_SUPG
static FEI2dTrLin	interp
static ParamKey	IPK_TR1_2D_SUPG_vof
static ParamKey	IPK_TR1_2D_SUPG_pvof
static ParamKey	IPK_TR1_2D_SUPG_mat0
static ParamKey	IPK_TR1_2D_SUPG_mat1
Static Protected Attributes inherited from oofem::SUPGElement
static ParamKey	IPK_SUPGElement_bsides
static ParamKey	IPK_SUPGElement_bcodes

Detailed Description

Class representing 2d linear triangular element for solving incompressible fluid with SUPG solver

This class is similar to TR1_2D_SUPG, but the difference is in handling multiple fluids. This class uses the interface position within an element to perform an integration for each fluid separately when evaluating contributing terms. It does not rely on rule of mixture which interpolates the properties using VOF value, but uses separate integration on each fluid volume.

Todo: Save & restore context will not work (some way how to save/restore dynamic integration rules have to be found: element has to restore these rules based on restored data, and then rules will be restored, or integration rule can store more info about gauss points if desired (when dynamic).

Todo: Integrate sub_IPRule as ordinary element integration rule.

Definition at line 65 of file tr1_2d_supg2.h.

Constructor & Destructor Documentation

◆ TR1_2D_SUPG2()

oofem::TR1_2D_SUPG2::TR1_2D_SUPG2	(	int	n,
		Domain *	d )

Definition at line 74 of file tr1_2d_supg2.C.

References oofem::Element::numberOfDofMans, and oofem::TR1_2D_SUPG::TR1_2D_SUPG().

Member Function Documentation

◆ _giveMaterial()

Material * oofem::TR1_2D_SUPG2::_giveMaterial ( int indx )

inlineprotected

Definition at line 182 of file tr1_2d_supg2.h.

References oofem::FEMComponent::domain, and mat.

Referenced by computeAccelerationTerm_MB(), computeAdvectionDerivativeTerm_MB(), computeAdvectionTerm_MB(), computeBCRhsTerm_MB(), computeDeviatoricStress(), computeDiffusionDerivativeTerm_MB(), computeDiffusionTerm_MB(), computeLSICStabilizationTerm_MB(), computePressureTerm_MC(), computeTangent(), printOutputAt(), and updateStabilizationCoeffs().

◆ computeAccelerationTerm_MB()

void oofem::TR1_2D_SUPG2::computeAccelerationTerm_MB	(	FloatMatrix &	answer,
		TimeStep *	tStep )

overridevirtual

Computes acceleration terms (generalized mass matrix with stabilization terms) for momentum balance equations(s).

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 186 of file tr1_2d_supg2.C.

References _giveMaterial(), oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::TR1_2D_SUPG::b, oofem::TR1_2D_SUPG::c, computeNVector(), oofem::FMElement::computeVectorOfVelocities(), computeVolumeAroundID(), oofem::TimeStep::givePreviousStep(), oofem::Element::integrationRulesArray, OOFEM_ERROR, oofem::FloatMatrix::resize(), oofem::SUPGElement::t_supg, vcoords, and oofem::FloatMatrix::zero().

◆ computeAccelerationTerm_MC()

void oofem::TR1_2D_SUPG2::computeAccelerationTerm_MC	(	FloatMatrix &	answer,
		TimeStep *	tStep )

overridevirtual

Computes acceleration terms for mass conservation equation.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 756 of file tr1_2d_supg2.C.

◆ computeAdvectionDerivativeTerm_MB()

void oofem::TR1_2D_SUPG2::computeAdvectionDerivativeTerm_MB	(	FloatMatrix &	answer,
		TimeStep *	tStep )

overridevirtual

Computes the derivative of advection terms for momentum balance equations(s) with respect to nodal velocities.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 363 of file tr1_2d_supg2.C.

References _giveMaterial(), oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::TR1_2D_SUPG::b, oofem::TR1_2D_SUPG::c, computeNVector(), oofem::FMElement::computeVectorOfVelocities(), computeVolumeAroundID(), oofem::TimeStep::givePreviousStep(), oofem::Element::integrationRulesArray, OOFEM_ERROR, oofem::FloatMatrix::resize(), oofem::SUPGElement::t_supg, vcoords, and oofem::FloatMatrix::zero().

Referenced by computeAdvectionTerm_MB().

◆ computeAdvectionDerivativeTerm_MC()

void oofem::TR1_2D_SUPG2::computeAdvectionDerivativeTerm_MC	(	FloatMatrix &	answer,
		TimeStep *	tStep )

overridevirtual

Computes the derivative of advection terms for mass conservation equation with respect to nodal velocities.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 715 of file tr1_2d_supg2.C.

References oofem::TR1_2D_SUPG::area, oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::TR1_2D_SUPG::b, oofem::TR1_2D_SUPG::c, oofem::FMElement::computeVectorOfVelocities(), oofem::TimeStep::givePreviousStep(), oofem::FloatMatrix::resize(), oofem::SUPGElement::t_pspg, and oofem::FloatMatrix::zero().

◆ computeAdvectionTerm_MB()

void oofem::TR1_2D_SUPG2::computeAdvectionTerm_MB	(	FloatArray &	answer,
		TimeStep *	tStep )

overridevirtual

Computes nonlinear advection terms for momentum balance equations(s).

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 282 of file tr1_2d_supg2.C.

References _giveMaterial(), oofem::FloatArray::at(), oofem::TR1_2D_SUPG::b, oofem::FloatArray::beProductOf(), oofem::TR1_2D_SUPG::c, computeAdvectionDerivativeTerm_MB(), computeNVector(), oofem::FMElement::computeVectorOfVelocities(), computeVolumeAroundID(), oofem::TimeStep::givePreviousStep(), oofem::Element::integrationRulesArray, oofem::FEMComponent::number, OOFEM_ERROR, oofem::FloatArray::resize(), oofem::SUPGElement::t_supg, vcoords, and oofem::FloatArray::zero().

◆ computeAdvectionTerm_MC()

void oofem::TR1_2D_SUPG2::computeAdvectionTerm_MC	(	FloatArray &	answer,
		TimeStep *	tStep )

overridevirtual

Computes advection terms for mass conservation equation.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 688 of file tr1_2d_supg2.C.

References oofem::TR1_2D_SUPG::area, oofem::FloatArray::at(), oofem::TR1_2D_SUPG::b, oofem::TR1_2D_SUPG::c, oofem::FMElement::computeVectorOfVelocities(), oofem::TimeStep::givePreviousStep(), oofem::FloatArray::resize(), and oofem::SUPGElement::t_pspg.

◆ computeBCRhsTerm_MB()

void oofem::TR1_2D_SUPG2::computeBCRhsTerm_MB	(	FloatArray &	answer,
		TimeStep *	tStep )

overridevirtual

Computes Rhs terms due to boundary conditions.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 802 of file tr1_2d_supg2.C.

References _giveMaterial(), oofem::FloatArray::at(), oofem::TR1_2D_SUPG::b, oofem::Element::bodyLoadArray, oofem::BodyLoadBGT, oofem::Element::boundaryLoadArray, oofem::TR1_2D_SUPG::c, oofem::Load::computeComponentArrayAt(), computeNVector(), oofem::BoundaryLoad::computeValueAt(), oofem::FMElement::computeVectorOfVelocities(), computeVolumeAroundID(), oofem::FEMComponent::domain, oofem::ForceLoadBVT, oofem::GeneralBoundaryCondition::giveBCGeoType(), oofem::GeneralBoundaryCondition::giveBCValType(), oofem::Element::giveBodyLoadArray(), oofem::Element::giveBoundaryLoadArray(), oofem::Element::giveNode(), oofem::TimeStep::givePreviousStep(), oofem::FloatArray::giveSize(), id, oofem::Element::integrationRulesArray, length(), OOFEM_ERROR, oofem::FloatArray::resize(), oofem::SUPGElement::t_supg, vcoords, and oofem::FloatArray::zero().

◆ computeBCRhsTerm_MC()

void oofem::TR1_2D_SUPG2::computeBCRhsTerm_MC	(	FloatArray &	answer,
		TimeStep *	tStep )

overridevirtual

Computes Rhs terms due to boundary conditions.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 901 of file tr1_2d_supg2.C.

◆ computeCriticalTimeStep()

double oofem::TR1_2D_SUPG2::computeCriticalTimeStep ( TimeStep * tStep )

overridevirtual

Computes the critical time increment.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1168 of file tr1_2d_supg2.C.

◆ computeDeviatoricStress()

void oofem::TR1_2D_SUPG2::computeDeviatoricStress	(	FloatArray &	answer,
		const FloatArray &	eps,
		GaussPoint *	gp,
		TimeStep *	tStep )

overrideprotectedvirtual

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1196 of file tr1_2d_supg2.C.

References _giveMaterial(), and oofem::LEPlicElementInterface::temp_vof.

◆ computeDiffusionDerivativeTerm_MB()

void oofem::TR1_2D_SUPG2::computeDiffusionDerivativeTerm_MB	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		TimeStep *	tStep )

overridevirtual

Computes the derivative of diffusion terms for momentum balance equations(s) with respect to nodal velocities.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 506 of file tr1_2d_supg2.C.

References _giveMaterial(), oofem::FloatMatrix::at(), oofem::TR1_2D_SUPG::b, oofem::FloatMatrix::beProductOf(), oofem::TR1_2D_SUPG::c, computeVolumeAroundID(), oofem::FEMComponent::domain, oofem::FluidModel::giveReynoldsNumber(), oofem::Element::integrationRulesArray, mat, OOFEM_ERROR, oofem::FloatMatrix::plusProductSymmUpper(), oofem::FloatMatrix::resize(), oofem::FloatMatrix::symmetrized(), oofem::FloatMatrix::times(), vcoords, and oofem::FloatMatrix::zero().

◆ computeDiffusionDerivativeTerm_MC()

void oofem::TR1_2D_SUPG2::computeDiffusionDerivativeTerm_MC	(	FloatMatrix &	answer,
		TimeStep *	tStep )

inlineoverridevirtual

Computes diffusion derivative terms for mass conservation equation.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 100 of file tr1_2d_supg2.h.

References oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

◆ computeDiffusionTerm_MB()

void oofem::TR1_2D_SUPG2::computeDiffusionTerm_MB	(	FloatArray &	answer,
		TimeStep *	tStep )

overridevirtual

Computes diffusion terms for momentum balance equations(s).

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 444 of file tr1_2d_supg2.C.

References _giveMaterial(), oofem::FloatArray::at(), oofem::TR1_2D_SUPG::b, oofem::TR1_2D_SUPG::c, oofem::FMElement::computeVectorOfVelocities(), computeVolumeAroundID(), oofem::FEMComponent::domain, oofem::FluidModel::giveReynoldsNumber(), oofem::Element::integrationRulesArray, mat, OOFEM_ERROR, oofem::FloatArray::resize(), vcoords, and oofem::FloatArray::zero().

◆ computeDiffusionTerm_MC()

void oofem::TR1_2D_SUPG2::computeDiffusionTerm_MC	(	FloatArray &	answer,
		TimeStep *	tStep )

inlineoverridevirtual

Computes diffusion terms for mass conservation equation.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 104 of file tr1_2d_supg2.h.

References oofem::FloatArray::resize(), and oofem::FloatArray::zero().

◆ computeGaussPoints()

void oofem::TR1_2D_SUPG2::computeGaussPoints ( )

overrideprotectedvirtual

Initializes the array of integration rules member variable. Element can have multiple integration rules for different tasks. For example structural element family class uses this feature to implement transparent support for reduced and selective integration of some strain components. Must be defined by terminator classes.

See also: IntegrationRule

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 154 of file tr1_2d_supg2.C.

References defaultIRule, oofem::Element::giveCrossSection(), and oofem::Element::integrationRulesArray.

Referenced by postInitialize().

◆ computeLEPLICVolumeFraction()

double oofem::TR1_2D_SUPG2::computeLEPLICVolumeFraction	(	const FloatArray &	n,
		const double	p,
		LEPlic *	matInterface,
		bool	updFlag )

overridevirtual

Computes corresponding volume fraction to given interface position.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1250 of file tr1_2d_supg2.C.

References computeMyVolume(), oofem::Polygon::computeVolume(), formVolumeInterfacePoly(), OOFEM_WARNING, and oofem::LEPlicElementInterface::p.

◆ computeLinearAdvectionTerm_MC()

void oofem::TR1_2D_SUPG2::computeLinearAdvectionTerm_MC	(	FloatMatrix &	answer,
		TimeStep *	tStep )

overridevirtual

Computes the linear advection term for mass conservation equation.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 682 of file tr1_2d_supg2.C.

◆ computeLSICStabilizationTerm_MB()

void oofem::TR1_2D_SUPG2::computeLSICStabilizationTerm_MB	(	FloatMatrix &	answer,
		TimeStep *	tStep )

overridevirtual

Computes SLIC stabilization term for momentum balance equation(s).

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 639 of file tr1_2d_supg2.C.

References _giveMaterial(), oofem::FloatMatrix::at(), oofem::TR1_2D_SUPG::b, oofem::TR1_2D_SUPG::c, computeVolumeAroundID(), oofem::Element::integrationRulesArray, OOFEM_ERROR, oofem::FloatMatrix::resize(), oofem::SUPGElement::t_lsic, vcoords, and oofem::FloatMatrix::zero().

◆ computeMyVolume()

double oofem::TR1_2D_SUPG2::computeMyVolume	(	LEPlic *	matInterface,
		bool	updFlag )

overridevirtual

Computes the volume of receiver.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1555 of file tr1_2d_supg2.C.

References oofem::TR1_2D_SUPG::area, oofem::Element::giveNode(), oofem::FEMComponent::giveNumber(), oofem::LEPlic::giveUpdatedXCoordinate(), and oofem::LEPlic::giveUpdatedYCoordinate().

Referenced by computeLEPLICVolumeFraction().

◆ computeNumberOfDofs()

int oofem::TR1_2D_SUPG2::computeNumberOfDofs ( )

overridevirtual

Computes or simply returns total number of element's local DOFs. Must be defined by particular element.

Returns: Number of local DOFs of element.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 89 of file tr1_2d_supg2.C.

◆ computeNVector()

void oofem::TR1_2D_SUPG2::computeNVector	(	FloatArray &	answer,
		GaussPoint *	gp )

protected

Definition at line 82 of file tr1_2d_supg2.C.

References oofem::GaussPoint::giveNaturalCoordinates(), and oofem::TR1_2D_SUPG::interp.

Referenced by computeAccelerationTerm_MB(), computeAdvectionDerivativeTerm_MB(), computeAdvectionTerm_MB(), computeBCRhsTerm_MB(), and updateStabilizationCoeffs().

◆ computePressureTerm_MB()

void oofem::TR1_2D_SUPG2::computePressureTerm_MB	(	FloatMatrix &	answer,
		TimeStep *	tStep )

overridevirtual

Computes pressure terms for momentum balance equations(s).

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 584 of file tr1_2d_supg2.C.

References oofem::TR1_2D_SUPG::area, oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::TR1_2D_SUPG::b, oofem::TR1_2D_SUPG::c, oofem::FMElement::computeVectorOfPressures(), oofem::FMElement::computeVectorOfVelocities(), oofem::TimeStep::givePreviousStep(), oofem::LEPlicElementInterface::p, oofem::FloatMatrix::resize(), oofem::SUPGElement::t_supg, and oofem::FloatMatrix::zero().

◆ computePressureTerm_MC()

void oofem::TR1_2D_SUPG2::computePressureTerm_MC	(	FloatMatrix &	answer,
		TimeStep *	tStep )

overridevirtual

Computes pressure terms for mass conservation equation.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 762 of file tr1_2d_supg2.C.

References _giveMaterial(), oofem::FloatMatrix::at(), oofem::TR1_2D_SUPG::b, oofem::TR1_2D_SUPG::c, computeVolumeAroundID(), oofem::Element::integrationRulesArray, OOFEM_ERROR, oofem::FloatMatrix::resize(), oofem::SUPGElement::t_pspg, vcoords, and oofem::FloatMatrix::zero().

◆ computeTangent()

void oofem::TR1_2D_SUPG2::computeTangent	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep )

overrideprotectedvirtual

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1206 of file tr1_2d_supg2.C.

References _giveMaterial(), and oofem::LEPlicElementInterface::temp_vof.

◆ computeVolumeAroundID()

double oofem::TR1_2D_SUPG2::computeVolumeAroundID	(	GaussPoint *	gp,
		integrationDomain	id,
		const std::vector< FloatArray > &	idpoly )

protected

◆ drawRawGeometry()

void oofem::TR1_2D_SUPG2::drawRawGeometry	(	oofegGraphicContext &	gc,
		TimeStep *	tStep )

overridevirtual

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1897 of file tr1_2d_supg2.C.

References gc, oofem::Element::giveNode(), OOFEG_RAW_GEOMETRY_LAYER, OOFEG_RAW_GEOMETRY_WIDTH, and oofem::LEPlicElementInterface::p.

◆ drawScalar()

void oofem::TR1_2D_SUPG2::drawScalar	(	oofegGraphicContext &	gc,
		TimeStep *	tStep )

overridevirtual

◆ EIPrimaryFieldI_evaluateFieldVectorAt()

int oofem::TR1_2D_SUPG2::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.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1595 of file tr1_2d_supg2.C.

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

◆ formMaterialVolumePoly()

void oofem::TR1_2D_SUPG2::formMaterialVolumePoly	(	Polygon &	matvolpoly,
		LEPlic *	matInterface,
		const FloatArray &	normal,
		const double	p,
		bool	updFlag )

overridevirtual

Assembles the true element material polygon (takes receiver vof into accout).

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1265 of file tr1_2d_supg2.C.

References oofem::Polygon::addVertex(), oofem::Polygon::clear(), formVolumeInterfacePoly(), oofem::Element::giveNode(), oofem::FEMComponent::giveNumber(), oofem::LEPlic::giveUpdatedXCoordinate(), oofem::LEPlic::giveUpdatedYCoordinate(), oofem::LEPlicElementInterface::normal, oofem::LEPlicElementInterface::p, oofem::Vertex::setCoords(), TRSUPG_ZERO_VOF, and oofem::LEPlicElementInterface::vof.

Referenced by drawScalar().

◆ formMyVolumePoly()

void oofem::TR1_2D_SUPG2::formMyVolumePoly	(	Polygon &	myPoly,
		LEPlic *	mat_interface,
		bool	updFlag )

overridevirtual

Assembles receiver volume.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1532 of file tr1_2d_supg2.C.

References oofem::Polygon::addVertex(), oofem::Polygon::clear(), oofem::Element::giveNode(), oofem::FEMComponent::giveNumber(), oofem::LEPlic::giveUpdatedXCoordinate(), oofem::LEPlic::giveUpdatedYCoordinate(), and oofem::Vertex::setCoords().

Referenced by truncateMatVolume().

◆ formVolumeInterfacePoly()

void oofem::TR1_2D_SUPG2::formVolumeInterfacePoly	(	Polygon &	matvolpoly,
		LEPlic *	matInterface,
		const FloatArray &	normal,
		const double	p,
		bool	updFlag )

overridevirtual

Assembles receiver material polygon based solely on given interface line.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1297 of file tr1_2d_supg2.C.

References oofem::Polygon::addVertex(), oofem::Polygon::clear(), oofem::Element::giveNode(), oofem::FEMComponent::giveNumber(), oofem::LEPlic::giveUpdatedXCoordinate(), oofem::LEPlic::giveUpdatedYCoordinate(), oofem::LEPlicElementInterface::normal, oofem::LEPlicElementInterface::p, and oofem::Vertex::setCoords().

Referenced by computeLEPLICVolumeFraction(), and formMaterialVolumePoly().

◆ giveClassName()

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

inlineoverridevirtual

Returns: Class name of the receiver.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 118 of file tr1_2d_supg2.h.

◆ giveDefaultIntegrationRule()

int oofem::TR1_2D_SUPG2::giveDefaultIntegrationRule ( ) const

inlineoverridevirtual

Returns id of default integration rule. Various element types can use different integration rules for implementation of selective or reduced integration of selected components. One particular integration rule from defined integration rules is default. There may be some operations (defined by parent analysis type class) which use default integration rule.

Returns: Id of default integration rule. (index into integrationRulesArray).

Reimplemented from oofem::Element.

Definition at line 156 of file tr1_2d_supg2.h.

◆ giveDefaultIntegrationRulePtr()

IntegrationRule * oofem::TR1_2D_SUPG2::giveDefaultIntegrationRulePtr ( )

inlineoverridevirtual

Access method for default integration rule.

Returns: Pointer to default integration rule.

See also: giveDefaultIntegrationRule

Reimplemented from oofem::Element.

Definition at line 157 of file tr1_2d_supg2.h.

References defaultIRule.

◆ giveDofManDofIDMask()

void oofem::TR1_2D_SUPG2::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::TR1_2D_SUPG.

Definition at line 95 of file tr1_2d_supg2.C.

◆ giveElement()

Element * oofem::TR1_2D_SUPG2::giveElement ( )

inlineoverridevirtual

Return number of receiver's element.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 144 of file tr1_2d_supg2.h.

◆ giveElementCenter()

void oofem::TR1_2D_SUPG2::giveElementCenter	(	LEPlic *	mat_interface,
		FloatArray &	center,
		bool	updFlag )

overridevirtual

Computes the receiver center (in updated Lagrangian configuration).

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1573 of file tr1_2d_supg2.C.

References oofem::FloatArray::add(), oofem::FloatArray::at(), oofem::Element::giveNode(), oofem::FEMComponent::giveNumber(), oofem::LEPlic::giveUpdatedCoordinate(), oofem::FloatArray::resize(), oofem::FloatArray::times(), and oofem::FloatArray::zero().

◆ giveInputRecord()

void oofem::TR1_2D_SUPG2::giveInputRecord ( DynamicInputRecord & input )

overridevirtual

Setups the input record string of receiver.

Parameters

input Dynamic input record to be filled by receiver.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 130 of file tr1_2d_supg2.C.

References oofem::TR1_2D_SUPG::IPK_TR1_2D_SUPG_mat0, oofem::TR1_2D_SUPG::IPK_TR1_2D_SUPG_mat1, oofem::TR1_2D_SUPG::IPK_TR1_2D_SUPG_pvof, oofem::TR1_2D_SUPG::IPK_TR1_2D_SUPG_vof, mat, oofem::LEPlicElementInterface::permanentVofFlag, oofem::DynamicInputRecord::setField(), and oofem::LEPlicElementInterface::vof.

◆ giveInputRecordName()

const char * oofem::TR1_2D_SUPG2::giveInputRecordName ( ) const

inlineoverridevirtual

Returns: Input record name of the receiver.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 119 of file tr1_2d_supg2.h.

References _IFT_TR1_2D_SUPG2_Name.

◆ giveInterface()

Interface * oofem::TR1_2D_SUPG2::giveInterface ( InterfaceType t )

overridevirtual

Interface requesting service.

See also: InterfaceType

Returns: Requested interface if implemented, otherwise NULL.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1175 of file tr1_2d_supg2.C.

References oofem::EIPrimaryFieldInterfaceType, oofem::LEPlicElementInterfaceType, oofem::NodalAveragingRecoveryModelInterfaceType, oofem::SpatialLocalizerInterfaceType, oofem::SPRNodalRecoveryModelInterfaceType, and oofem::ZZNodalRecoveryModelInterfaceType.

◆ giveInternalStateAtNode()

int oofem::TR1_2D_SUPG2::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::TR1_2D_SUPG.

Definition at line 1879 of file tr1_2d_supg2.C.

Referenced by drawScalar().

◆ giveIPValue()

int oofem::TR1_2D_SUPG2::giveIPValue	(	FloatArray &	answer,
		GaussPoint *	gp,
		InternalStateType	type,
		TimeStep *	tStep )

overridevirtual

Returns the integration point corresponding value in full form.

Parameters

answer	Contain corresponding integration point value, zero sized if not available.
gp	Integration point to check.
type	Determines the type of internal variable.
tStep	Time step.

Returns: Nonzero if o.k, zero otherwise.

Todo: Which "error type" should be used? Why are there several? I don't see the point of this enum when there could be different function calls just as well (and different IST values)

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1781 of file tr1_2d_supg2.C.

References oofem::FloatArray::at(), oofem::LEPlicElementInterface::giveTempVolumeFraction(), and oofem::FloatArray::resize().

Referenced by drawScalar(), and NodalAveragingRecoveryMI_computeNodalValue().

◆ initializeFrom()

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

overridevirtual

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 101 of file tr1_2d_supg2.C.

References oofem::FEMComponent::giveDomain(), oofem::TR1_2D_SUPG::IPK_TR1_2D_SUPG_mat0, oofem::TR1_2D_SUPG::IPK_TR1_2D_SUPG_mat1, oofem::TR1_2D_SUPG::IPK_TR1_2D_SUPG_pvof, oofem::TR1_2D_SUPG::IPK_TR1_2D_SUPG_vof, mat, oofem::Element::material, oofem::FEMComponent::number, PM_UPDATE_PARAMETER, PM_UPDATE_PARAMETER_AND_REPORT, oofem::LEPlicElementInterface::setPermanentVolumeFraction(), oofem::LEPlicElementInterface::temp_vof, and oofem::LEPlicElementInterface::vof.

◆ NodalAveragingRecoveryMI_computeNodalValue()

void oofem::TR1_2D_SUPG2::NodalAveragingRecoveryMI_computeNodalValue	(	FloatArray &	answer,
		int	node,
		InternalStateType	type,
		TimeStep *	tStep )

overridevirtual

Computes the element value in given node.

Parameters

answer	Contains the result.
node	Element node number.
type	Determines the type of internal variable to be recovered.
tStep	Time step.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1793 of file tr1_2d_supg2.C.

References giveIPValue(), and oofem::Element::integrationRulesArray.

◆ postInitialize()

void oofem::TR1_2D_SUPG2::postInitialize ( )

overrideprotectedvirtual

Performs post initialization steps.

Reimplemented from oofem::Element.

Definition at line 145 of file tr1_2d_supg2.C.

References computeGaussPoints(), oofem::TR1_2D_SUPG::initGeometry(), and updateIntegrationRules().

◆ printOutputAt()

void oofem::TR1_2D_SUPG2::printOutputAt	(	FILE *	file,
		TimeStep *	tStep )

overridevirtual

Prints output of receiver to stream, for given time step. This is used for output into the standard output file.

Parameters

file	File pointer to print to.
tStep	Time step to write for.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1843 of file tr1_2d_supg2.C.

References _giveMaterial(), oofem::LEPlicElementInterface::giveVolumeFraction(), oofem::Element::integrationRulesArray, and oofem::LEPlicElementInterface::vof.

◆ restoreContext()

void oofem::TR1_2D_SUPG2::restoreContext	(	DataStream &	stream,
		ContextMode	mode )

overridevirtual

Restores the receiver state previously written in stream.

See also: saveContext

Parameters

stream	Input stream.
mode	Determines amount of info available in stream (state, definition, ...).

Exceptions

throws an ContextIOERR exception if error encountered.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1870 of file tr1_2d_supg2.C.

◆ saveContext()

void oofem::TR1_2D_SUPG2::saveContext	(	DataStream &	stream,
		ContextMode	mode )

overridevirtual

Stores receiver state to output stream.

Parameters

stream	Output stream.
mode	Determines amount of info required in stream (state, definition, ...).

Exceptions

throws an ContextIOERR exception if error encountered.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1863 of file tr1_2d_supg2.C.

◆ SPRNodalRecoveryMI_giveDofMansDeterminedByPatch()

void oofem::TR1_2D_SUPG2::SPRNodalRecoveryMI_giveDofMansDeterminedByPatch	(	IntArray &	answer,
		int	pap )

overridevirtual

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1816 of file tr1_2d_supg2.C.

References oofem::IntArray::at(), oofem::Element::giveNode(), oofem::FEMComponent::giveNumber(), OOFEM_ERROR, and oofem::IntArray::resize().

◆ SPRNodalRecoveryMI_giveNumberOfIP()

int oofem::TR1_2D_SUPG2::SPRNodalRecoveryMI_giveNumberOfIP ( )

overridevirtual

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1829 of file tr1_2d_supg2.C.

◆ SPRNodalRecoveryMI_givePatchType()

SPRPatchType oofem::TR1_2D_SUPG2::SPRNodalRecoveryMI_givePatchType ( )

overridevirtual

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1836 of file tr1_2d_supg2.C.

References oofem::SPRPatchType_2dxy.

◆ SPRNodalRecoveryMI_giveSPRAssemblyPoints()

void oofem::TR1_2D_SUPG2::SPRNodalRecoveryMI_giveSPRAssemblyPoints ( IntArray & pap )

overridevirtual

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1807 of file tr1_2d_supg2.C.

References oofem::IntArray::at(), oofem::Element::giveNode(), and oofem::IntArray::resize().

◆ truncateMatVolume()

double oofem::TR1_2D_SUPG2::truncateMatVolume	(	const Polygon &	matvolpoly,
		double &	volume )

overridevirtual

Truncates given material polygon to receiver.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1514 of file tr1_2d_supg2.C.

References oofem::TR1_2D_SUPG::area, oofem::Graph::clip(), oofem::Polygon::computeVolume(), oofem::Polygon::draw(), formMyVolumePoly(), gc, and OOFEG_DEBUG_LAYER.

◆ updateElementForNewInterfacePosition()

void oofem::TR1_2D_SUPG2::updateElementForNewInterfacePosition ( TimeStep * tStep )

inlineoverridevirtual

Reimplemented from oofem::SUPGElement.

Definition at line 114 of file tr1_2d_supg2.h.

References updateIntegrationRules().

◆ updateIntegrationRules()

void oofem::TR1_2D_SUPG2::updateIntegrationRules ( )

protected

◆ updateStabilizationCoeffs()

void oofem::TR1_2D_SUPG2::updateStabilizationCoeffs ( TimeStep * tStep )

overridevirtual

Updates the stabilization coefficients used for CBS and SUPG algorithms.

Parameters

tStep Active time step.

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 908 of file tr1_2d_supg2.C.

References _giveMaterial(), oofem::TR1_2D_SUPG::area, oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::TR1_2D_SUPG::b, oofem::FloatArray::beProductOf(), oofem::TR1_2D_SUPG::c, oofem::FloatMatrix::computeFrobeniusNorm(), oofem::FloatArray::computeNorm(), computeNVector(), oofem::FMElement::computeVectorOfVelocities(), computeVolumeAroundID(), oofem::FEMComponent::domain, oofem::TimeStep::givePreviousStep(), oofem::TimeStep::giveTimeIncrement(), oofem::Element::integrationRulesArray, oofem::max(), oofem::FloatMatrix::resize(), oofem::sum(), oofem::SUPGElement::t_lsic, oofem::SUPGElement::t_pspg, oofem::SUPGElement::t_supg, oofem::FloatArray::times(), vcoords, oofem::LEPlicElementInterface::vof, oofem::VST_Density, oofem::VST_Length, oofem::VST_Time, oofem::VST_Velocity, oofem::VST_Viscosity, and oofem::FloatMatrix::zero().

◆ updateVolumePolygons()

void oofem::TR1_2D_SUPG2::updateVolumePolygons	(	Polygon &	referenceFluidPoly,
		Polygon &	secondFluidPoly,
		int &	rfPoints,
		int &	sfPoints,
		const FloatArray &	normal,
		const double	p,
		bool	updFlag )

protected

Definition at line 1406 of file tr1_2d_supg2.C.

References oofem::Polygon::addVertex(), oofem::Polygon::clear(), oofem::Element::giveNode(), oofem::LEPlicElementInterface::normal, oofem::LEPlicElementInterface::p, and oofem::Vertex::setCoords().

Referenced by updateIntegrationRules().

◆ updateYourself()

void oofem::TR1_2D_SUPG2::updateYourself ( 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

Reimplemented from oofem::TR1_2D_SUPG.

Definition at line 1638 of file tr1_2d_supg2.C.

Member Data Documentation

◆ defaultIRule

std :: unique_ptr< IntegrationRule > oofem::TR1_2D_SUPG2::defaultIRule

protected

default integration rule over element volume. the standard integrationRulesArray contains two rules on element subvolumes.

Definition at line 85 of file tr1_2d_supg2.h.

Referenced by computeGaussPoints(), and giveDefaultIntegrationRulePtr().

◆ id

integrationDomain oofem::TR1_2D_SUPG2::id[2]

protected

Definition at line 75 of file tr1_2d_supg2.h.

Referenced by computeBCRhsTerm_MB().

◆ mat

int oofem::TR1_2D_SUPG2::mat[2]

protected

mat[0] reference fluid. mat[1] second fluid.

Definition at line 80 of file tr1_2d_supg2.h.

Referenced by _giveMaterial(), computeDiffusionDerivativeTerm_MB(), computeDiffusionTerm_MB(), giveInputRecord(), and initializeFrom().

◆ myPoly

Polygon oofem::TR1_2D_SUPG2::myPoly[2]

protected

myPoly[0] occupied by reference fluid. myPoly[1] occupied by second fluid (air).

Definition at line 72 of file tr1_2d_supg2.h.

Referenced by updateIntegrationRules().

◆ vcoords

std::vector< FloatArray > oofem::TR1_2D_SUPG2::vcoords[2]

protected

Definition at line 73 of file tr1_2d_supg2.h.

Referenced by computeAccelerationTerm_MB(), computeAdvectionDerivativeTerm_MB(), computeAdvectionTerm_MB(), computeBCRhsTerm_MB(), computeDiffusionDerivativeTerm_MB(), computeDiffusionTerm_MB(), computeLSICStabilizationTerm_MB(), computePressureTerm_MC(), updateIntegrationRules(), and updateStabilizationCoeffs().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TR1_2D_SUPG2()

Member Function Documentation

◆ _giveMaterial()

◆ computeAccelerationTerm_MB()

◆ computeAccelerationTerm_MC()

◆ computeAdvectionDerivativeTerm_MB()

◆ computeAdvectionDerivativeTerm_MC()

◆ computeAdvectionTerm_MB()

◆ computeAdvectionTerm_MC()

◆ computeBCRhsTerm_MB()

◆ computeBCRhsTerm_MC()

◆ computeCriticalTimeStep()

◆ computeDeviatoricStress()

◆ computeDiffusionDerivativeTerm_MB()

◆ computeDiffusionDerivativeTerm_MC()

◆ computeDiffusionTerm_MB()

◆ computeDiffusionTerm_MC()

◆ computeGaussPoints()

◆ computeLEPLICVolumeFraction()

◆ computeLinearAdvectionTerm_MC()

◆ computeLSICStabilizationTerm_MB()

◆ computeMyVolume()

◆ computeNumberOfDofs()

◆ computeNVector()

◆ computePressureTerm_MB()

◆ computePressureTerm_MC()

◆ computeTangent()

◆ computeVolumeAroundID()

◆ drawRawGeometry()

◆ drawScalar()

◆ EIPrimaryFieldI_evaluateFieldVectorAt()

◆ formMaterialVolumePoly()

◆ formMyVolumePoly()

◆ formVolumeInterfacePoly()

◆ giveClassName()

◆ giveDefaultIntegrationRule()

◆ giveDefaultIntegrationRulePtr()

◆ giveDofManDofIDMask()

◆ giveElement()

◆ giveElementCenter()

◆ giveInputRecord()

◆ giveInputRecordName()

◆ giveInterface()

◆ giveInternalStateAtNode()

◆ giveIPValue()

◆ initializeFrom()

◆ NodalAveragingRecoveryMI_computeNodalValue()

◆ postInitialize()

◆ printOutputAt()

◆ restoreContext()

◆ saveContext()

◆ SPRNodalRecoveryMI_giveDofMansDeterminedByPatch()

◆ SPRNodalRecoveryMI_giveNumberOfIP()

◆ SPRNodalRecoveryMI_givePatchType()

◆ SPRNodalRecoveryMI_giveSPRAssemblyPoints()

◆ truncateMatVolume()

◆ updateElementForNewInterfacePosition()

◆ updateIntegrationRules()

◆ updateStabilizationCoeffs()

◆ updateVolumePolygons()

◆ updateYourself()

Member Data Documentation

◆ defaultIRule

◆ id

◆ mat

◆ myPoly

◆ vcoords