Class representing 2d quadrilateral element with linear velocity and constant pressure approximations for solving incompressible fluid problems with SUPG solver. More...

#include <quad10_2d_supg.h>

Inheritance diagram for oofem::Quad10_2D_SUPG:

Collaboration diagram for oofem::Quad10_2D_SUPG:

Public Member Functions
	Quad10_2D_SUPG (int n, Domain *d)

virtual	~Quad10_2D_SUPG ()

virtual FEInterpolation *	giveInterpolation () const

virtual FEInterpolation *	giveInterpolation (DofIDItem id) const
	Returns the interpolation for the specific dof id. More...

virtual const char *	giveClassName () const

virtual const char *	giveInputRecordName () const

virtual MaterialMode	giveMaterialMode ()
	Returns material mode for receiver integration points. More...

virtual void	giveInternalDofManDofIDMask (int i, IntArray &answer) const
	Returns internal dofmanager dof mask for node. More...

virtual void	giveDofManDofIDMask (int inode, IntArray &answer) const
	Returns dofmanager dof mask for node. More...

virtual int	computeNumberOfDofs ()
	Computes or simply returns total number of element's local DOFs. More...

virtual IRResultType	initializeFrom (InputRecord *ir)
	Initializes receiver according to object description stored in input record. More...

virtual void	giveInputRecord (DynamicInputRecord &input)
	Setups the input record string of receiver. More...

virtual void	updateYourself (TimeStep *tStep)
	Updates element state after equilibrium in time step has been reached. More...

virtual int	checkConsistency ()
	Performs consistency check. More...

virtual contextIOResultType	saveContext (DataStream &stream, ContextMode mode, void *obj=NULL)
	Stores receiver state to output stream. More...

virtual contextIOResultType	restoreContext (DataStream &stream, ContextMode mode, void *obj=NULL)
	Restores the receiver state previously written in stream. More...

virtual double	LS_PCS_computeF (LevelSetPCS ls, TimeStep tStep)
	Evaluates F in level set equation of the form $\phi_t + F(\nabla\phi, x) \|\nabla\phi\| = 0$ where for interface position driven by flow with speed u: $F = u\cdot \frac{\nabla\phi}{\|\nabla\phi\|}$ . More...

virtual void	LS_PCS_computedN (FloatMatrix &answer)
	Returns gradient of shape functions. More...

virtual double	LS_PCS_computeVolume ()
	Returns receiver's volume. More...

virtual double	LS_PCS_computeS (LevelSetPCS ls, TimeStep tStep)
	Evaluates S in level set equation of the form $\phi_t = S(\phi) (1-\|\nabla\phi\|) = 0$ where $S=\frac{\phi}{\sqrt{\phi^2+\epsilon^2}}$ . More...

virtual void	LS_PCS_computeVOFFractions (FloatArray &answer, FloatArray &fi)
	Returns VOF fractions for each material on element according to nodal values of level set function (passed as parameter) More...

virtual Interface *	giveInterface (InterfaceType t)
	Interface requesting service. More...

virtual void	NodalAveragingRecoveryMI_computeNodalValue (FloatArray &answer, int node, InternalStateType type, TimeStep *tStep)
	Computes the element value in given node. More...

virtual int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep)
	Returns the integration point corresponding value in full form. More...

int	giveInternalStateAtNode (FloatArray &answer, InternalStateType type, InternalStateMode mode, int node, TimeStep *tStep)
	Returns internal state variable (like stress,strain) at node of element in Reduced form, the way how is obtained is dependent on InternalValueType. More...

virtual void	drawRawGeometry (oofegGraphicContext &gc, TimeStep *tStep)

virtual void	drawScalar (oofegGraphicContext &gc, TimeStep *tStep)

virtual double	computeCriticalTimeStep (TimeStep *tStep)
	Computes the critical time increment. More...

virtual void	computeAdvectionTerm (FloatMatrix &answer, TimeStep *tStep)

virtual void	computeAdvectionDeltaTerm (FloatMatrix &answer, TimeStep *tStep)

virtual void	computeMassDeltaTerm (FloatMatrix &answer, TimeStep *tStep)

virtual void	computeLSICTerm (FloatMatrix &answer, TimeStep *tStep)

virtual void	computeAdvectionEpsilonTerm (FloatMatrix &answer, TimeStep *tStep)

virtual void	computeMassEpsilonTerm (FloatMatrix &answer, TimeStep *tStep)

virtual int	giveNumberOfInternalDofManagers () const

virtual DofManager *	giveInternalDofManager (int i) const
	Returns i-th internal element dof manager of the receiver. More...

Helping functions for computing VOFFractions.
void	computeIntersection (int iedge, FloatArray &intcoords, FloatArray &fi)

void	computeMiddlePointOnParabolicArc (FloatArray &answer, int iedge, FloatArray borderpoints)

void	computeCenterOf (FloatArray &C, FloatArray c, int dim)

void	computeQuadraticRoots (FloatArray Coeff, double &r1, double &r2)

void	computeCoordsOfEdge (FloatArray &answer, int iedge)

void	computeQuadraticFunct (FloatArray &answer, int iedge)

void	computeQuadraticFunct (FloatArray &answer, FloatArray line)

Public Member Functions inherited from oofem::SUPGElement2
	SUPGElement2 (int n, Domain *aDomain)

virtual	~SUPGElement2 ()

virtual void	giveCharacteristicMatrix (FloatMatrix &answer, CharType, TimeStep *tStep)
	Computes characteristic matrix of receiver of requested type in given time step. More...

virtual void	giveCharacteristicVector (FloatArray &answer, CharType, ValueModeType, TimeStep *tStep)
	Computes characteristic vector of receiver of requested type in given time step. More...

virtual void	updateElementForNewInterfacePosition (TimeStep *tStep)

virtual void	computeAccelerationTerm_MB (FloatMatrix &answer, TimeStep *tStep)
	Computes acceleration terms (generalized mass matrix with stabilization terms) for momentum balance equations(s). More...

virtual void	computeAdvectionTerm_MB (FloatArray &answer, TimeStep *tStep)
	Computes nonlinear advection terms for momentum balance equations(s). More...

virtual void	computeAdvectionDerivativeTerm_MB (FloatMatrix &answer, TimeStep *tStep)
	Computes the derivative of advection terms for momentum balance equations(s) with respect to nodal velocities. More...

virtual void	computeDiffusionTerm_MB (FloatArray &answer, TimeStep *tStep)
	Computes diffusion terms for momentum balance equations(s). More...

virtual void	computeDiffusionDerivativeTerm_MB (FloatMatrix &answer, MatResponseMode mode, TimeStep *tStep)
	Computes the derivative of diffusion terms for momentum balance equations(s) with respect to nodal velocities. More...

virtual void	computePressureTerm_MB (FloatMatrix &answer, TimeStep *tStep)
	Computes pressure terms for momentum balance equations(s). More...

virtual void	computeLSICStabilizationTerm_MB (FloatMatrix &answer, TimeStep *tStep)
	Computes SLIC stabilization term for momentum balance equation(s). More...

virtual void	computeLinearAdvectionTerm_MC (FloatMatrix &answer, TimeStep *tStep)
	Computes the linear advection term for mass conservation equation. More...

virtual void	computeAdvectionTerm_MC (FloatArray &answer, TimeStep *tStep)
	Computes advection terms for mass conservation equation. More...

virtual void	computeAdvectionDerivativeTerm_MC (FloatMatrix &answer, TimeStep *tStep)
	Computes the derivative of advection terms for mass conservation equation with respect to nodal velocities. More...

virtual void	computeDiffusionDerivativeTerm_MC (FloatMatrix &answer, TimeStep *tStep)
	Computes diffusion derivative terms for mass conservation equation. More...

virtual void	computeDiffusionTerm_MC (FloatArray &answer, TimeStep *tStep)
	Computes diffusion terms for mass conservation equation. More...

virtual void	computeAccelerationTerm_MC (FloatMatrix &answer, TimeStep *tStep)
	Computes acceleration terms for mass conservation equation. More...

virtual void	computePressureTerm_MC (FloatMatrix &answer, TimeStep *tStep)
	Computes pressure terms for mass conservation equation. More...

virtual void	computeBCRhsTerm_MB (FloatArray &answer, TimeStep *tStep)
	Computes Rhs terms due to boundary conditions. More...

virtual void	computeBCRhsTerm_MC (FloatArray &answer, TimeStep *tStep)
	Computes Rhs terms due to boundary conditions. More...

virtual void	computeLoadVector (FloatArray &answer, BodyLoad load, CharType type, ValueModeType mode, TimeStep tStep)
	Computes the contribution of the given body load (volumetric). More...

virtual void	updateInternalState (TimeStep *tStep)
	Updates element state after equilibrium in time step has been reached. More...

int	giveInternalStateAtNode (FloatArray &answer, InternalStateType type, InternalStateMode mode, int node, TimeStep *tStep)
	Returns internal state variable (like stress,strain) at node of element in Reduced form, the way how is obtained is dependent on InternalValueType. More...

Public Member Functions inherited from oofem::SUPGElement
	SUPGElement (int n, Domain *aDomain)

virtual	~SUPGElement ()

virtual void	computeBCLhsPressureTerm_MC (FloatMatrix &answer, TimeStep *tStep)
	Computes Lhs terms due to boundary conditions - pressure. More...

virtual void	computeBCLhsTerm_MB (FloatMatrix &answer, TimeStep *tStep)
	Computes Lhs terms due to boundary conditions - velocity. More...

virtual void	computeBCLhsPressureTerm_MB (FloatMatrix &answer, TimeStep *tStep)
	Computes Lhs terms due to boundary conditions - pressure. More...

virtual void	computeSlipWithFrictionBCTerm_MB (FloatMatrix &answer, Load load, int side, TimeStep tStep)
	Computes Lhs term due to applied slip with friction bc. More...

virtual void	computePenetrationWithResistanceBCTerm_MB (FloatMatrix &answer, Load load, int side, TimeStep tStep)
	Computes Lhs contribution due to applied Penetration bc. More...

virtual void	computeOutFlowBCTerm_MB (FloatMatrix &answer, int side, TimeStep *tStep)
	Computes Lhs contribution due to outflow BC. More...

virtual void	computeHomogenizedReinforceTerm_MB (FloatMatrix &answer, Load load, TimeStep tStep)

virtual void	computeHomogenizedReinforceTerm_MC (FloatMatrix &answer, Load load, TimeStep tStep)

Public Member Functions inherited from oofem::FMElement
	FMElement (int n, Domain *aDomain)

virtual	~FMElement ()

virtual void	computeVectorOfVelocities (ValueModeType mode, TimeStep *tStep, FloatArray &velocities)

virtual void	computeVectorOfPressures (ValueModeType mode, TimeStep *tStep, FloatArray &pressures)

Public Member Functions inherited from oofem::Element
	Element (int n, Domain *aDomain)
	Constructor. More...

	Element (const Element &src)=delete

Element &	operator= (const Element &src)=delete

virtual	~Element ()
	Virtual destructor. More...

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)
	Returns internal state variable (like stress,strain) at side of element in Reduced form If side is possessing DOFs, otherwise recover techniques will not work due to absence of side-shape functions. More...

virtual void	showSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep)
	Shows sparse structure. More...

virtual void	showExtendedSparseMtrxStructure (CharType mtrx, oofegGraphicContext &gc, TimeStep *tStep)
	Shows extended sparse structure (for example, due to nonlocal interactions for tangent stiffness) More...

int	giveLabel () const

int	giveGlobalNumber () const

void	setGlobalNumber (int num)
	Sets receiver globally unique number. More...

elementParallelMode	giveParallelMode () const
	Return elementParallelMode of receiver. More...

void	setParallelMode (elementParallelMode _mode)
	Sets parallel mode of element. More...

virtual elementParallelMode	giveKnotSpanParallelMode (int) const
	Returns the parallel mode for particular knot span of the receiver. More...

int	packUnknowns (DataStream &buff, TimeStep *tStep)
	Pack all necessary data of element (according to its parallel_mode) integration points into given communication buffer. More...

int	unpackAndUpdateUnknowns (DataStream &buff, TimeStep *tStep)
	Unpack and updates all necessary data of element (according to its parallel_mode) integration points into given communication buffer. More...

int	estimatePackSize (DataStream &buff)
	Estimates the necessary pack size to hold all packed data of receiver. More...

const IntArray *	givePartitionList () const
	Returns partition list of receiver. More...

void	setPartitionList (IntArray &pl)
	Sets partition list of receiver. More...

virtual double	predictRelativeComputationalCost ()
	Returns the weight representing relative computational cost of receiver The reference element is triangular plane stress element with linear approximation, single integration point and linear isotropic material. More...

virtual double	giveRelativeSelfComputationalCost ()
	Returns the weight representing relative computational cost of receiver The reference element is triangular plane stress element. More...

virtual double	predictRelativeRedistributionCost ()
	Returns the relative redistribution cost of the receiver. More...

IntArray *	giveBodyLoadArray ()
	Returns array containing load numbers of loads acting on element. More...

IntArray *	giveBoundaryLoadArray ()
	Returns array containing load numbers of boundary loads acting on element. More...

virtual void	printOutputAt (FILE file, TimeStep tStep)
	Prints output of receiver to stream, for given time step. More...

void	giveLocationArray (IntArray &locationArray, const UnknownNumberingScheme &s, IntArray *dofIds=NULL) const
	Returns the location array (array of code numbers) of receiver for given numbering scheme. More...

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)
	Returns the location array for the boundary of the element. More...

virtual void	giveBoundaryLocationArray (IntArray &locationArray, const IntArray &bNodes, const IntArray &dofIDMask, const UnknownNumberingScheme &s, IntArray *dofIds=NULL)

virtual int	giveNumberOfDofs ()

virtual double	giveCharacteristicValue (CharType type, TimeStep *tStep)
	Computes characteristic value of receiver of requested type in given time step. More...

virtual void	computeBoundarySurfaceLoadVector (FloatArray &answer, BoundaryLoad load, int boundary, CharType type, ValueModeType mode, TimeStep tStep, bool global=true)
	Computes the contribution of the given load at the given boundary surface in global coordinate system. More...

virtual void	computeTangentFromSurfaceLoad (FloatMatrix &answer, SurfaceLoad load, int boundary, MatResponseMode rmode, TimeStep tStep)
	Computes the tangent contribution of the given load at the given boundary. More...

virtual void	computeTangentFromEdgeLoad (FloatMatrix &answer, EdgeLoad load, int boundary, MatResponseMode rmode, TimeStep tStep)
	Computes the tangent contribution of the given load at the given boundary. More...

virtual void	computeBoundaryEdgeLoadVector (FloatArray &answer, BoundaryLoad load, int edge, CharType type, ValueModeType mode, TimeStep tStep, bool global=true)
	Computes the contribution of the given load at the given boundary edge. More...

const IntArray &	giveBodyLoadList () const
	Returns receiver list of bodyloads. More...

const IntArray &	giveBoundaryLoadList () const
	Returns receiver list of boundary loads. More...

void	computeVectorOf (ValueModeType u, TimeStep *tStep, FloatArray &answer)
	Returns local vector of unknowns. More...

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)
	Boundary version of computeVectorOf. More...

void	computeVectorOf (PrimaryField &field, const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false)
	Returns local vector of unknowns. More...

void	computeVectorOfPrescribed (ValueModeType u, TimeStep *tStep, FloatArray &answer)
	Returns local vector of prescribed unknowns. More...

void	computeVectorOfPrescribed (const IntArray &dofIDMask, ValueModeType type, TimeStep *tStep, FloatArray &answer)
	Returns local vector of prescribed unknowns. More...

virtual int	computeNumberOfGlobalDofs ()
	Computes the total number of element's global dofs. More...

int	computeNumberOfPrimaryMasterDofs ()
	Computes the total number of element's primary master DOFs. More...

virtual bool	computeGtoLRotationMatrix (FloatMatrix &answer)
	Returns transformation matrix from global c.s. More...

virtual bool	giveRotationMatrix (FloatMatrix &answer)
	Transformation matrices updates rotation matrix between element-local and primary DOFs, taking into account nodal c.s. More...

virtual bool	computeDofTransformationMatrix (FloatMatrix &answer, const IntArray &nodes, bool includeInternal)
	Returns transformation matrix for DOFs from global coordinate system to local coordinate system in nodes. More...

virtual void	giveElementDofIDMask (IntArray &answer) const
	Returns element dof mask for node. More...

virtual void	computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer)
	Computes the unknown vector interpolated at the specified local coordinates. More...

virtual double	computeVolumeAreaOrLength ()
	Computes the volume, area or length of the element depending on its spatial dimension. More...

double	computeMeanSize ()
	Computes the size of the element defined as its length. More...

virtual double	computeVolume ()
	Computes the volume. More...

virtual double	computeArea ()
	Computes the area (zero for all but 2d geometries). More...

virtual double	computeLength ()
	Computes the length (zero for all but 1D geometries) More...

virtual void	giveBoundaryEdgeNodes (IntArray &bNodes, int boundary)
	Returns list of receiver boundary nodes for given edge. More...

virtual void	giveBoundarySurfaceNodes (IntArray &bNodes, int boundary)
	Returns list of receiver boundary nodes for given surface. More...

virtual IntegrationRule *	giveBoundaryEdgeIntegrationRule (int order, int boundary)
	Returns boundary edge integration rule. More...

virtual IntegrationRule *	giveBoundarySurfaceIntegrationRule (int order, int boundary)
	Returns boundary surface integration rule. More...

int	giveDofManagerNumber (int i) const
	Translates local to global indices for dof managers. More...

const IntArray &	giveDofManArray () const

void	addDofManager (DofManager *dMan)

DofManager *	giveDofManager (int i) const

Node *	giveNode (int i) const
	Returns reference to the i-th node of element. More...

virtual ElementSide *	giveSide (int i) const
	Returns reference to the i-th side of element. More...

virtual Material *	giveMaterial ()

int	giveMaterialNumber () const

CrossSection *	giveCrossSection ()

void	setMaterial (int matIndx)
	Sets the material of receiver. More...

virtual void	setCrossSection (int csIndx)
	Sets the cross section model of receiver. More...

virtual int	giveNumberOfDofManagers () const

virtual int	giveNumberOfNodes () const
	Returns number of nodes of receiver. More...

void	setDofManagers (const IntArray &dmans)
	Sets receiver dofManagers. More...

void	setBodyLoads (const IntArray &bodyLoads)
	Sets receiver bodyLoadArray. More...

void	setIntegrationRules (std::vector< std::unique_ptr< IntegrationRule > > irlist)
	Sets integration rules. More...

virtual integrationDomain	giveIntegrationDomain () const
	Returns integration domain for receiver, used to initialize integration point over receiver volume. More...

virtual int	giveIntegrationRuleLocalCodeNumbers (IntArray &answer, IntegrationRule &ie)
	Assembles the code numbers of given integration element (sub-patch) This is done by obtaining list of nonzero shape functions and by collecting the code numbers of nodes corresponding to these shape functions. More...

int	giveRegionNumber ()

virtual void	postInitialize ()
	Performs post initialization steps. More...

virtual void	initializeYourself (TimeStep *timeStepWhenICApply)
	Initialization according to state given by initial conditions. More...

virtual bool	isActivated (TimeStep *tStep)

virtual bool	isCast (TimeStep *tStep)

virtual void	initForNewStep ()
	Initializes receivers state to new time step. More...

virtual Element_Geometry_Type	giveGeometryType () const
	Returns the element geometry type. More...

virtual int	giveSpatialDimension ()
	Returns the element spatial dimension (1, 2, or 3). More...

virtual int	giveNumberOfBoundarySides ()

virtual int	giveDefaultIntegrationRule () const
	Returns id of default integration rule. More...

virtual IntegrationRule *	giveDefaultIntegrationRulePtr ()
	Access method for default integration rule. More...

int	giveNumberOfIntegrationRules ()

virtual IntegrationRule *	giveIntegrationRule (int i)

std::vector< std::unique_ptr< IntegrationRule > > &	giveIntegrationRulesArray ()

virtual int	testElementExtension (ElementExtension ext)
	Tests if the element implements required extension. More...

int	giveGlobalIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep)

virtual double	giveLengthInDir (const FloatArray &normalToCrackPlane)
	Default implementation returns length of element projection into specified direction. More...

virtual double	giveCharacteristicLength (const FloatArray &normalToCrackPlane)
	Returns the size of element in the given direction, in some cases adjusted (e.g. More...

double	giveCharacteristicLengthForPlaneElements (const FloatArray &normalToCrackPlane)
	Returns the size of element in the given direction if the direction is in the XY plane, otherwise gives the mean size defined as the square root of the element area. More...

double	giveCharacteristicLengthForAxisymmElements (const FloatArray &normalToCrackPlane)
	Returns the size of an axisymmetric element in the given direction if the direction is in the XY plane, otherwise gives the mean distance vrom the symmetry axis multiplied by pi. More...

virtual double	giveCharacteristicSize (GaussPoint *gp, FloatArray &normalToCrackPlane, ElementCharSizeMethod method)
	Returns characteristic element size for a given integration point and given direction. More...

virtual double	giveParentElSize () const
	Returns the size (length, area or volume depending on element type) of the parent element. More...

virtual void	updateBeforeNonlocalAverage (TimeStep *tStep)
	Updates internal element state (in all integration points of receiver) before nonlocal averaging takes place. More...

virtual int	computeGlobalCoordinates (FloatArray &answer, const FloatArray &lcoords)
	Computes the global coordinates from given element's local coordinates. More...

virtual bool	computeLocalCoordinates (FloatArray &answer, const FloatArray &gcoords)
	Computes the element local coordinates from given global coordinates. More...

virtual int	giveLocalCoordinateSystem (FloatMatrix &answer)
	Returns local coordinate system of receiver Required by material models with ortho- and anisotrophy. More...

virtual void	computeMidPlaneNormal (FloatArray &answer, const GaussPoint *gp)
	Computes mid-plane normal of receiver at integration point. More...

virtual int	adaptiveMap (Domain oldd, TimeStep tStep)
	Initializes the internal state variables stored in all IPs according to state in given domain. More...

virtual int	mapStateVariables (Domain &iOldDom, const TimeStep &iTStep)
	Maps the internal state variables stored in all IPs from the old domain to the new domain. More...

virtual int	adaptiveUpdate (TimeStep *tStep)
	Updates the internal state variables stored in all IPs according to already mapped state. More...

virtual int	adaptiveFinish (TimeStep *tStep)
	Finishes the mapping for given time step. More...

virtual void	updateLocalNumbering (EntityRenumberingFunctor &f)
	Local renumbering support. More...

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. More...

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. More...

Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
	Regular constructor, creates component with given number and belonging to given domain. More...

virtual	~FEMComponent ()
	Virtual destructor. More...

Domain *	giveDomain () const

virtual void	setDomain (Domain *d)
	Sets associated Domain. More...

int	giveNumber () const

void	setNumber (int num)
	Sets number of receiver. More...

virtual void	printYourself ()
	Prints receiver state on stdout. Useful for debugging. More...

std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros). More...

Public Member Functions inherited from oofem::LevelSetPCSElementInterface
	LevelSetPCSElementInterface ()

Public Member Functions inherited from oofem::Interface
	Interface ()
	Constructor. More...

virtual	~Interface ()

Public Member Functions inherited from oofem::ZZNodalRecoveryModelInterface
	ZZNodalRecoveryModelInterface (Element *element)
	Constructor. More...

virtual bool	ZZNodalRecoveryMI_computeNValProduct (FloatMatrix &answer, InternalStateType type, TimeStep *tStep)
	Computes the element contribution to $\int_\Omega N^{\mathrm{T}}\alpha\;\mathrm{d}\Omega$ , where $\alpha$ is quantity to be recovered (for example stress or strain vector). More...

virtual void	ZZNodalRecoveryMI_computeNNMatrix (FloatArray &answer, InternalStateType type)
	Computes the element contribution to $\int_\Omega N^{\mathrm{T}} \cdot N\;\mathrm{d}\Omega$ term. More...

Public Member Functions inherited from oofem::NodalAveragingRecoveryModelInterface
	NodalAveragingRecoveryModelInterface ()
	Constructor. More...

Protected Member Functions
virtual void	giveLocalVelocityDofMap (IntArray &map)

virtual void	giveLocalPressureDofMap (IntArray &map)

virtual void	computeGaussPoints ()
	Initializes the array of integration rules member variable. More...

virtual void	computeNuMatrix (FloatMatrix &answer, GaussPoint *gp)

virtual void	computeUDotGradUMatrix (FloatMatrix &answer, GaussPoint gp, TimeStep tStep)

virtual void	computeBMatrix (FloatMatrix &anwer, GaussPoint *gp)

virtual void	computeDivUMatrix (FloatMatrix &answer, GaussPoint *gp)

virtual void	computeNpMatrix (FloatMatrix &answer, GaussPoint *gp)

virtual void	computeGradPMatrix (FloatMatrix &answer, GaussPoint *gp)

virtual void	computeDivTauMatrix (FloatMatrix &answer, GaussPoint gp, TimeStep tStep)

virtual void	computeGradUMatrix (FloatMatrix &answer, GaussPoint gp, TimeStep tStep)

virtual int	giveNumberOfSpatialDimensions ()

virtual double	computeVolumeAround (GaussPoint *gp)
	Returns volume related to given integration point. More...

virtual void	updateStabilizationCoeffs (TimeStep *tStep)
	Updates the stabilization coefficients used for CBS and SUPG algorithms. More...

Protected Member Functions inherited from oofem::SUPGElement2
virtual void	computeDeviatoricStrain (FloatArray &answer, GaussPoint gp, TimeStep tStep)

virtual void	computeDeviatoricStress (FloatArray &answer, GaussPoint gp, TimeStep tStep)

virtual void	computeEdgeLoadVector_MB (FloatArray &answer, Load load, int id, TimeStep tStep)

virtual void	computeSurfaceLoadVector_MB (FloatArray &answer, Load load, int id, TimeStep tStep)

virtual void	computeEdgeLoadVector_MC (FloatArray &answer, Load load, int id, TimeStep tStep)

virtual void	computeSurfaceLoadVector_MC (FloatArray &answer, Load load, int id, TimeStep tStep)

Protected Attributes
ElementDofManager	pressureNode

Protected Attributes inherited from oofem::SUPGElement
IntArray	boundarySides
	Array of boundary sides. More...

IntArray	boundaryCodes
	Boundary sides codes. More...

double	t_supg
	Stabilization coefficients, updated for each solution step in updateStabilizationCoeffs() More...

double	t_pspg

double	t_lsic

Protected Attributes inherited from oofem::Element
int	numberOfDofMans
	Number of dofmanagers. More...

IntArray	dofManArray
	Array containing dofmanager numbers. More...

int	material
	Number of associated material. More...

int	crossSection
	Number of associated cross section. More...

IntArray	bodyLoadArray
	Array containing indexes of loads (body loads and boundary loads are kept separately), that apply on receiver. More...

IntArray	boundaryLoadArray

std::vector< std::unique_ptr< IntegrationRule > >	integrationRulesArray
	List of integration rules of receiver (each integration rule contains associated integration points also). More...

FloatMatrix	elemLocalCS
	Transformation material matrix, used in orthotropic and anisotropic materials, global->local transformation. More...

int	activityTimeFunction
	Element activity time function. If defined, nonzero value indicates active receiver, zero value inactive element. More...

int	globalNumber
	In parallel mode, globalNumber contains globally unique DoFManager number. More...

int	numberOfGaussPoints
	Number of integration points as specified by nip. More...

elementParallelMode	parallel_mode
	Determines the parallel mode of the element. More...

IntArray	partitions
	List of partition sharing the shared element or remote partition containing remote element counterpart. More...

Protected Attributes inherited from oofem::FEMComponent
int	number
	Component number. More...

Domain *	domain
	Link to domain object, useful for communicating with other FEM components. More...

Static Protected Attributes
static FEI2dQuadLin	velocityInterpolation

static FEI2dQuadConst	pressureInterpolation

Detailed Description

Class representing 2d quadrilateral element with linear velocity and constant pressure approximations for solving incompressible fluid problems with SUPG solver.

Definition at line 57 of file quad10_2d_supg.h.

Constructor & Destructor Documentation

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

Definition at line 67 of file quad10_2d_supg.C.

References oofem::Element::numberOfDofMans.

oofem::Quad10_2D_SUPG::~Quad10_2D_SUPG ( )

virtual

Definition at line 73 of file quad10_2d_supg.C.

Member Function Documentation

int oofem::Quad10_2D_SUPG::checkConsistency ( )

virtual

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::SUPGElement2.

Definition at line 529 of file quad10_2d_supg.C.

References oofem::SUPGElement::checkConsistency().

Referenced by giveMaterialMode().

void oofem::Quad10_2D_SUPG::computeAdvectionDeltaTerm	(	FloatMatrix &	answer,
		TimeStep *	tStep
	)

virtual

Definition at line 377 of file quad10_2d_supg.C.

References oofem::FloatMatrix::clear(), computeNuMatrix(), computeUDotGradUMatrix(), computeVolumeAround(), oofem::Element::giveCrossSection(), oofem::Element::integrationRulesArray, and oofem::FloatMatrix::plusProductUnsym().

Referenced by LS_PCS_computeVolume(), and updateStabilizationCoeffs().

void oofem::Quad10_2D_SUPG::computeAdvectionEpsilonTerm	(	FloatMatrix &	answer,
		TimeStep *	tStep
	)

virtual

Definition at line 433 of file quad10_2d_supg.C.

References oofem::FloatMatrix::clear(), computeGradPMatrix(), computeUDotGradUMatrix(), computeVolumeAround(), oofem::Element::integrationRulesArray, and oofem::FloatMatrix::plusProductUnsym().

Referenced by LS_PCS_computeVolume(), and updateStabilizationCoeffs().

void oofem::Quad10_2D_SUPG::computeAdvectionTerm	(	FloatMatrix &	answer,
		TimeStep *	tStep
	)

virtual

Definition at line 359 of file quad10_2d_supg.C.

References oofem::FloatMatrix::clear(), computeNuMatrix(), computeUDotGradUMatrix(), computeVolumeAround(), oofem::Element::giveCrossSection(), oofem::Element::integrationRulesArray, and oofem::FloatMatrix::plusProductUnsym().

Referenced by LS_PCS_computeVolume(), and updateStabilizationCoeffs().

void oofem::Quad10_2D_SUPG::computeBMatrix	(	FloatMatrix &	anwer,
		GaussPoint *	gp
	)

protectedvirtual

Implements oofem::SUPGElement2.

Definition at line 187 of file quad10_2d_supg.C.

References oofem::FloatMatrix::at(), oofem::FEI2dQuadLin::evaldNdx(), oofem::GaussPoint::giveNaturalCoordinates(), oofem::FloatMatrix::resize(), velocityInterpolation, and oofem::FloatMatrix::zero().

Referenced by giveNumberOfInternalDofManagers().

void oofem::Quad10_2D_SUPG::computeCenterOf	(	FloatArray &	C,
		FloatArray	c,
		int	dim
	)

Referenced by LS_PCS_computeVolume().

void oofem::Quad10_2D_SUPG::computeCoordsOfEdge	(	FloatArray &	answer,
		int	iedge
	)

Referenced by LS_PCS_computeVolume().

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

virtual

Computes the critical time increment.

Implements oofem::SUPGElement2.

Definition at line 513 of file quad10_2d_supg.C.

Referenced by LS_PCS_computeVolume().

void oofem::Quad10_2D_SUPG::computeDivTauMatrix	(	FloatMatrix &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

protectedvirtual

Implements oofem::SUPGElement2.

Definition at line 269 of file quad10_2d_supg.C.

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

Referenced by giveNumberOfInternalDofManagers().

void oofem::Quad10_2D_SUPG::computeDivUMatrix	(	FloatMatrix &	answer,
		GaussPoint *	gp
	)

protectedvirtual

Implements oofem::SUPGElement2.

Definition at line 204 of file quad10_2d_supg.C.

References oofem::FloatMatrix::at(), oofem::FEI2dQuadLin::evaldNdx(), oofem::GaussPoint::giveNaturalCoordinates(), oofem::FloatMatrix::resize(), velocityInterpolation, and oofem::FloatMatrix::zero().

Referenced by computeLSICTerm(), and giveNumberOfInternalDofManagers().

void oofem::Quad10_2D_SUPG::computeGaussPoints ( )

protectedvirtual

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::Element.

Definition at line 138 of file quad10_2d_supg.C.

References oofem::Element::giveCrossSection(), oofem::Element::integrationRulesArray, and oofem::CrossSection::setupIntegrationPoints().

Referenced by giveNumberOfInternalDofManagers().

void oofem::Quad10_2D_SUPG::computeGradPMatrix	(	FloatMatrix &	answer,
		GaussPoint *	gp
	)

protectedvirtual

Implements oofem::SUPGElement2.

Definition at line 259 of file quad10_2d_supg.C.

References oofem::FloatMatrix::beTranspositionOf(), oofem::FEI2dQuadConst::evaldNdx(), oofem::GaussPoint::giveNaturalCoordinates(), and pressureInterpolation.

Referenced by computeAdvectionEpsilonTerm(), computeMassEpsilonTerm(), and giveNumberOfInternalDofManagers().

void oofem::Quad10_2D_SUPG::computeGradUMatrix	(	FloatMatrix &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

protectedvirtual

Implements oofem::SUPGElement2.

Definition at line 232 of file quad10_2d_supg.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FMElement::computeVectorOfVelocities(), oofem::FloatArray::dotProduct(), oofem::FEI2dQuadLin::evaldNdx(), oofem::GaussPoint::giveNaturalCoordinates(), oofem::FloatMatrix::resize(), velocityInterpolation, and oofem::FloatMatrix::zero().

Referenced by giveNumberOfInternalDofManagers().

void oofem::Quad10_2D_SUPG::computeIntersection	(	int	iedge,
		FloatArray &	intcoords,
		FloatArray &	fi
	)

Referenced by LS_PCS_computeVolume().

void oofem::Quad10_2D_SUPG::computeLSICTerm	(	FloatMatrix &	answer,
		TimeStep *	tStep
	)

virtual

Definition at line 414 of file quad10_2d_supg.C.

References oofem::FloatMatrix::clear(), computeDivUMatrix(), computeVolumeAround(), oofem::Element::giveCrossSection(), oofem::Element::integrationRulesArray, oofem::FloatMatrix::plusProductSymmUpper(), and oofem::FloatMatrix::symmetrized().

Referenced by LS_PCS_computeVolume(), and updateStabilizationCoeffs().

void oofem::Quad10_2D_SUPG::computeMassDeltaTerm	(	FloatMatrix &	answer,
		TimeStep *	tStep
	)

virtual

Definition at line 396 of file quad10_2d_supg.C.

References oofem::FloatMatrix::clear(), computeNuMatrix(), computeUDotGradUMatrix(), computeVolumeAround(), oofem::Element::giveCrossSection(), oofem::Element::integrationRulesArray, and oofem::FloatMatrix::plusProductUnsym().

Referenced by LS_PCS_computeVolume(), and updateStabilizationCoeffs().

void oofem::Quad10_2D_SUPG::computeMassEpsilonTerm	(	FloatMatrix &	answer,
		TimeStep *	tStep
	)

virtual

Definition at line 450 of file quad10_2d_supg.C.

References oofem::FloatMatrix::clear(), computeGradPMatrix(), computeNuMatrix(), computeVolumeAround(), oofem::Element::integrationRulesArray, and oofem::FloatMatrix::plusProductUnsym().

Referenced by LS_PCS_computeVolume(), and updateStabilizationCoeffs().

void oofem::Quad10_2D_SUPG::computeMiddlePointOnParabolicArc	(	FloatArray &	answer,
		int	iedge,
		FloatArray	borderpoints
	)

Referenced by LS_PCS_computeVolume().

void oofem::Quad10_2D_SUPG::computeNpMatrix	(	FloatMatrix &	answer,
		GaussPoint *	gp
	)

protectedvirtual

Implements oofem::SUPGElement2.

Definition at line 219 of file quad10_2d_supg.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FEI2dQuadConst::evalN(), oofem::GaussPoint::giveNaturalCoordinates(), pressureInterpolation, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by giveNumberOfInternalDofManagers().

void oofem::Quad10_2D_SUPG::computeNuMatrix	(	FloatMatrix &	answer,
		GaussPoint *	gp
	)

protectedvirtual

Implements oofem::SUPGElement2.

Definition at line 159 of file quad10_2d_supg.C.

References oofem::FloatMatrix::beNMatrixOf(), oofem::FEI2dQuadLin::evalN(), oofem::GaussPoint::giveNaturalCoordinates(), and velocityInterpolation.

Referenced by computeAdvectionDeltaTerm(), computeAdvectionTerm(), computeMassDeltaTerm(), computeMassEpsilonTerm(), computeUDotGradUMatrix(), and giveNumberOfInternalDofManagers().

int oofem::Quad10_2D_SUPG::computeNumberOfDofs ( )

virtual

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::Element.

Definition at line 100 of file quad10_2d_supg.C.

Referenced by giveMaterialMode().

void oofem::Quad10_2D_SUPG::computeQuadraticFunct	(	FloatArray &	answer,
		int	iedge
	)

Referenced by LS_PCS_computeVolume().

void oofem::Quad10_2D_SUPG::computeQuadraticFunct	(	FloatArray &	answer,
		FloatArray	line
	)

void oofem::Quad10_2D_SUPG::computeQuadraticRoots	(	FloatArray	Coeff,
		double &	r1,
		double &	r2
	)

Referenced by LS_PCS_computeVolume().

void oofem::Quad10_2D_SUPG::computeUDotGradUMatrix	(	FloatMatrix &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

protectedvirtual

Implements oofem::SUPGElement2.

Definition at line 168 of file quad10_2d_supg.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatArray::beProductOf(), computeNuMatrix(), oofem::FMElement::computeVectorOfVelocities(), oofem::FEI2dQuadLin::evaldNdx(), oofem::GaussPoint::giveNaturalCoordinates(), oofem::FloatMatrix::resize(), velocityInterpolation, and oofem::FloatMatrix::zero().

Referenced by computeAdvectionDeltaTerm(), computeAdvectionEpsilonTerm(), computeAdvectionTerm(), computeMassDeltaTerm(), and giveNumberOfInternalDofManagers().

double oofem::Quad10_2D_SUPG::computeVolumeAround ( GaussPoint * gp )

protectedvirtual

Returns volume related to given integration point.

Used typically in subroutines, that perform integration over element volume. Should be implemented by particular elements.

See also: GaussPoint

Parameters

gp	Integration point for which volume is computed.

Returns: Volume for integration point.

Reimplemented from oofem::Element.

Definition at line 596 of file quad10_2d_supg.C.

References oofem::FEMComponent::domain, oofem::GaussPoint::giveNaturalCoordinates(), oofem::FEInterpolation::giveTransformationJacobian(), oofem::GaussPoint::giveWeight(), and velocityInterpolation.

Referenced by computeAdvectionDeltaTerm(), computeAdvectionEpsilonTerm(), computeAdvectionTerm(), computeLSICTerm(), computeMassDeltaTerm(), computeMassEpsilonTerm(), and giveNumberOfInternalDofManagers().

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

virtual

Reimplemented from oofem::Element.

Definition at line 664 of file quad10_2d_supg.C.

References oofem::oofegGraphicContext::getElementColor(), oofem::oofegGraphicContext::getElementEdgeColor(), oofem::Node::giveCoordinate(), oofem::Element::giveNode(), OOFEG_RAW_GEOMETRY_LAYER, OOFEG_RAW_GEOMETRY_WIDTH, and oofem::oofegGraphicContext::testElementGraphicActivity().

Referenced by LS_PCS_computeVolume().

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

virtual

Reimplemented from oofem::Element.

Definition at line 695 of file quad10_2d_supg.C.

Referenced by LS_PCS_computeVolume().

virtual const char* oofem::Quad10_2D_SUPG::giveClassName ( ) const

inlinevirtual

Implements oofem::Interface.

Definition at line 72 of file quad10_2d_supg.h.

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

virtual

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 107 of file quad10_2d_supg.C.

Referenced by giveMaterialMode().

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

virtual

Setups the input record string of receiver.

Parameters

input Dynamic input record to be filled by receiver.

Reimplemented from oofem::SUPGElement2.

Definition at line 130 of file quad10_2d_supg.C.

References oofem::SUPGElement2::giveInputRecord(), oofem::DofManager::giveInputRecord(), and pressureNode.

Referenced by giveMaterialMode().

virtual const char* oofem::Quad10_2D_SUPG::giveInputRecordName ( ) const

inlinevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 73 of file quad10_2d_supg.h.

References _IFT_Quad10_2D_SUPG_Name.

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

virtual

Interface requesting service.

See also: InterfaceType

Returns: Requested interface if implemented, otherwise NULL.

Reimplemented from oofem::FEMComponent.

Definition at line 627 of file quad10_2d_supg.C.

References oofem::LevelSetPCSElementInterfaceType, oofem::NodalAveragingRecoveryModelInterfaceType, and oofem::ZZNodalRecoveryModelInterfaceType.

Referenced by LS_PCS_computeVolume().

DofManager * oofem::Quad10_2D_SUPG::giveInternalDofManager ( int i ) const

virtual

Returns i-th internal element dof manager of the receiver.

Parameters

i	Internal number of DOF.

Returns: DOF number i.

Reimplemented from oofem::Element.

Definition at line 93 of file quad10_2d_supg.C.

References pressureNode.

Referenced by giveNumberOfInternalDofManagers().

void oofem::Quad10_2D_SUPG::giveInternalDofManDofIDMask	(	int	inode,
		IntArray &	answer
	)		const

virtual

Returns internal 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 114 of file quad10_2d_supg.C.

Referenced by giveMaterialMode().

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

virtual

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 657 of file quad10_2d_supg.C.

References oofem::SUPGElement::giveInternalStateAtNode().

Referenced by drawScalar(), and LS_PCS_computeVolume().

FEInterpolation * oofem::Quad10_2D_SUPG::giveInterpolation ( ) const

virtual

Returns: Interpolation of the element geometry, or NULL if none exist.

Reimplemented from oofem::Element.

Definition at line 77 of file quad10_2d_supg.C.

References velocityInterpolation.

FEInterpolation * oofem::Quad10_2D_SUPG::giveInterpolation ( DofIDItem id ) const

virtual

Returns the interpolation for the specific dof id.

Special elements which uses a mixed interpolation should reimplement this method.

Parameters

id	ID of the dof for the for the requested interpolation.

Returns: Appropriate interpolation, or NULL if none exists.

Reimplemented from oofem::Element.

Definition at line 83 of file quad10_2d_supg.C.

References pressureInterpolation, and velocityInterpolation.

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

virtual

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::Element.

Definition at line 542 of file quad10_2d_supg.C.

References oofem::FloatArray::at(), oofem::FEMComponent::domain, oofem::GaussPoint::giveElement(), oofem::MaterialInterface::giveElementMaterialMixture(), oofem::Domain::giveEngngModel(), oofem::Element::giveIPValue(), oofem::EngngModel::giveMaterialInterface(), oofem::FEMComponent::giveNumber(), and oofem::Domain::giveNumber().

Referenced by drawScalar(), LS_PCS_computeVolume(), and NodalAveragingRecoveryMI_computeNodalValue().

void oofem::Quad10_2D_SUPG::giveLocalPressureDofMap ( IntArray & map )

protectedvirtual

Reimplemented from oofem::SUPGElement.

Definition at line 649 of file quad10_2d_supg.C.

Referenced by giveNumberOfInternalDofManagers().

void oofem::Quad10_2D_SUPG::giveLocalVelocityDofMap ( IntArray & map )

protectedvirtual

Reimplemented from oofem::SUPGElement.

Definition at line 642 of file quad10_2d_supg.C.

References oofem::IntArray::enumerate().

Referenced by giveNumberOfInternalDofManagers().

virtual MaterialMode oofem::Quad10_2D_SUPG::giveMaterialMode ( )

inlinevirtual

Returns material mode for receiver integration points.

Should be specialized.

Returns: Material mode of element.

Reimplemented from oofem::Element.

Definition at line 74 of file quad10_2d_supg.h.

References checkConsistency(), computeNumberOfDofs(), giveDofManDofIDMask(), giveInputRecord(), giveInternalDofManDofIDMask(), initializeFrom(), LS_PCS_computedN(), LS_PCS_computeF(), restoreContext(), saveContext(), and updateYourself().

virtual int oofem::Quad10_2D_SUPG::giveNumberOfInternalDofManagers ( ) const

inlinevirtual

Returns: Number of internal DOF managers of element.

Reimplemented from oofem::Element.

Definition at line 137 of file quad10_2d_supg.h.

References computeBMatrix(), computeDivTauMatrix(), computeDivUMatrix(), computeGaussPoints(), computeGradPMatrix(), computeGradUMatrix(), computeNpMatrix(), computeNuMatrix(), computeUDotGradUMatrix(), computeVolumeAround(), giveInternalDofManager(), giveLocalPressureDofMap(), giveLocalVelocityDofMap(), giveNumberOfSpatialDimensions(), and updateStabilizationCoeffs().

int oofem::Quad10_2D_SUPG::giveNumberOfSpatialDimensions ( )

protectedvirtual

Implements oofem::SUPGElement2.

Definition at line 467 of file quad10_2d_supg.C.

Referenced by giveNumberOfInternalDofManagers().

IRResultType oofem::Quad10_2D_SUPG::initializeFrom ( InputRecord * ir )

virtual

Initializes receiver according to object description stored in input record.

This function is called immediately after creating object using constructor. Input record can be imagined as data record in component database belonging to receiver. Receiver may use value-name extracting functions to extract particular field from record.

See also: IR_GIVE_FIELD; IR_GIVE_OPTIONAL_FIELD

Parameters

ir	Input record to initialize from.

Returns: IRResultType

Reimplemented from oofem::SUPGElement2.

Definition at line 121 of file quad10_2d_supg.C.

References oofem::SUPGElement2::initializeFrom(), oofem::ElementDofManager::initializeFrom(), and pressureNode.

Referenced by giveMaterialMode().

void oofem::Quad10_2D_SUPG::LS_PCS_computedN ( FloatMatrix & answer )

virtual

Returns gradient of shape functions.

Implements oofem::LevelSetPCSElementInterface.

Definition at line 503 of file quad10_2d_supg.C.

Referenced by giveMaterialMode().

double oofem::Quad10_2D_SUPG::LS_PCS_computeF	(	LevelSetPCS *	ls,
		TimeStep *	tStep
	)

virtual

Evaluates F in level set equation of the form

$\phi_t + F(\nabla\phi, x) |\nabla\phi| = 0$

where for interface position driven by flow with speed u:

$F = u\cdot \frac{\nabla\phi}{|\nabla\phi|}$

.

Implements oofem::LevelSetPCSElementInterface.

Definition at line 474 of file quad10_2d_supg.C.

References oofem::IntArray::at(), oofem::FloatArray::at(), oofem::FMElement::computeVectorOfVelocities(), oofem::Element::dofManArray, oofem::LevelSetPCS::giveLevelSetDofManValue(), and oofem::norm().

Referenced by giveMaterialMode().

double oofem::Quad10_2D_SUPG::LS_PCS_computeS	(	LevelSetPCS *	ls,
		TimeStep *	tStep
	)

virtual

Evaluates S in level set equation of the form

$\phi_t = S(\phi) (1-|\nabla\phi|) = 0$

where

$S=\frac{\phi}{\sqrt{\phi^2+\epsilon^2}}$

.

Implements oofem::LevelSetPCSElementInterface.

Definition at line 496 of file quad10_2d_supg.C.

Referenced by LS_PCS_computeVolume().

void oofem::Quad10_2D_SUPG::LS_PCS_computeVOFFractions	(	FloatArray &	answer,
		FloatArray &	fi
	)

virtual

Returns VOF fractions for each material on element according to nodal values of level set function (passed as parameter)

Implements oofem::LevelSetPCSElementInterface.

Definition at line 508 of file quad10_2d_supg.C.

Referenced by LS_PCS_computeVolume().

virtual double oofem::Quad10_2D_SUPG::LS_PCS_computeVolume ( )

inlinevirtual

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

virtual

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.

Implements oofem::NodalAveragingRecoveryModelInterface.

Definition at line 520 of file quad10_2d_supg.C.

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

Referenced by LS_PCS_computeVolume().

contextIOResultType oofem::Quad10_2D_SUPG::restoreContext	(	DataStream &	stream,
		ContextMode	mode,
		void *	obj = `NULL`
	)

virtual

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, ...).
obj	Special parameter for sending extra information.

Returns: contextIOResultType.

Exceptions

throws an ContextIOERR exception if error encountered.

Reimplemented from oofem::Element.

Definition at line 579 of file quad10_2d_supg.C.

References oofem::CIO_OK, oofem::Element::restoreContext(), and THROW_CIOERR.

Referenced by giveMaterialMode().

contextIOResultType oofem::Quad10_2D_SUPG::saveContext	(	DataStream &	stream,
		ContextMode	mode,
		void *	obj = `NULL`
	)

virtual

Stores receiver state to output stream.

Parameters

stream	Output stream.
mode	Determines amount of info required in stream (state, definition, ...).
obj	Special parameter, used only to send particular integration point to material class version of this method.

Returns: contextIOResultType.

Exceptions

throws an ContextIOERR exception if error encountered.

Reimplemented from oofem::Element.

Definition at line 562 of file quad10_2d_supg.C.

References oofem::CIO_OK, oofem::Element::saveContext(), and THROW_CIOERR.

Referenced by giveMaterialMode().

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

protectedvirtual

Updates the stabilization coefficients used for CBS and SUPG algorithms.

Parameters

tStep Active time step.

Reimplemented from oofem::SUPGElement.

Definition at line 277 of file quad10_2d_supg.C.

References computeAdvectionDeltaTerm(), computeAdvectionEpsilonTerm(), computeAdvectionTerm(), oofem::FloatMatrix::computeFrobeniusNorm(), oofem::SUPGElement2::computeLinearAdvectionTerm_MC(), computeLSICTerm(), computeMassDeltaTerm(), computeMassEpsilonTerm(), oofem::FloatArray::computeNorm(), oofem::FMElement::computeVectorOfVelocities(), oofem::Element::giveCrossSection(), oofem::TimeStep::giveTimeIncrement(), oofem::Element::integrationRulesArray, N, oofem::SUPGElement::t_lsic, oofem::SUPGElement::t_pspg, and oofem::SUPGElement::t_supg.

Referenced by giveNumberOfInternalDofManagers().

void oofem::Quad10_2D_SUPG::updateYourself ( TimeStep * tStep )

virtual

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.