OOFEM 2.4
OOFEM.org - Object Oriented Finite Element Solver
oofem::SUPGElement Class Referenceabstract

General stabilized SUPG/PSPG element for CFD analysis. More...

#include <supgelement.h>

+ Inheritance diagram for oofem::SUPGElement:
+ Collaboration diagram for oofem::SUPGElement:

Public Member Functions

 SUPGElement (int n, Domain *aDomain)
 
void initializeFrom (InputRecord &ir, int priority) override
 
void giveInputRecord (DynamicInputRecord &input) override
 Setups the input record string of receiver.
 
void giveCharacteristicMatrix (FloatMatrix &answer, CharType type, TimeStep *tStep) override
 Computes characteristic matrix of receiver of requested type in given time step.
 
void giveCharacteristicVector (FloatArray &answer, CharType type, ValueModeType mode, TimeStep *tStep) override
 Computes characteristic vector of receiver of requested type in given time step.
 
void updateStabilizationCoeffs (TimeStep *tStep) override
 Updates the stabilization coefficients used for CBS and SUPG algorithms.
 
virtual void updateElementForNewInterfacePosition (TimeStep *tStep)
 
virtual void computeAccelerationTerm_MB (FloatMatrix &answer, TimeStep *tStep)=0
 Computes acceleration terms (generalized mass matrix with stabilization terms) for momentum balance equations(s).
 
virtual void computeAdvectionTerm_MB (FloatArray &answer, TimeStep *tStep)=0
 Computes nonlinear advection terms for momentum balance equations(s).
 
virtual void computeAdvectionDerivativeTerm_MB (FloatMatrix &answer, TimeStep *tStep)=0
 Computes the derivative of advection terms for momentum balance equations(s) with respect to nodal velocities.
 
virtual void computeDiffusionTerm_MB (FloatArray &answer, TimeStep *tStep)=0
 Computes diffusion terms for momentum balance equations(s).
 
virtual void computeDiffusionDerivativeTerm_MB (FloatMatrix &answer, MatResponseMode mode, TimeStep *tStep)=0
 Computes the derivative of diffusion terms for momentum balance equations(s) with respect to nodal velocities.
 
virtual void computePressureTerm_MB (FloatMatrix &answer, TimeStep *tStep)=0
 Computes pressure terms for momentum balance equations(s).
 
virtual void computeBCLhsPressureTerm_MC (FloatMatrix &answer, TimeStep *tStep)
 Computes Lhs terms due to boundary conditions - pressure.
 
virtual void computeLSICStabilizationTerm_MB (FloatMatrix &answer, TimeStep *tStep)=0
 Computes SLIC stabilization term for momentum balance equation(s).
 
virtual void computeLinearAdvectionTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0
 Computes the linear advection term for mass conservation equation.
 
virtual void computeAdvectionTerm_MC (FloatArray &answer, TimeStep *tStep)=0
 Computes advection terms for mass conservation equation.
 
virtual void computeAdvectionDerivativeTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0
 Computes the derivative of advection terms for mass conservation equation with respect to nodal velocities.
 
virtual void computeDiffusionDerivativeTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0
 Computes diffusion derivative terms for mass conservation equation.
 
virtual void computeDiffusionTerm_MC (FloatArray &answer, TimeStep *tStep)=0
 Computes diffusion terms for mass conservation equation.
 
virtual void computeAccelerationTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0
 Computes acceleration terms for mass conservation equation.
 
virtual void computePressureTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0
 Computes pressure terms for mass conservation equation.
 
virtual void computeBCLhsTerm_MB (FloatMatrix &answer, TimeStep *tStep)
 Computes Lhs terms due to boundary conditions - velocity.
 
virtual void computeBCLhsPressureTerm_MB (FloatMatrix &answer, TimeStep *tStep)
 Computes Lhs terms due to boundary conditions - pressure.
 
virtual void computeBCRhsTerm_MB (FloatArray &answer, TimeStep *tStep)=0
 Computes Rhs terms due to boundary conditions.
 
virtual void computeBCRhsTerm_MC (FloatArray &answer, TimeStep *tStep)=0
 Computes Rhs terms due to boundary conditions.
 
virtual void computeSlipWithFrictionBCTerm_MB (FloatMatrix &answer, Load *load, int side, TimeStep *tStep)
 Computes Lhs term due to applied slip with friction bc.
 
virtual void computePenetrationWithResistanceBCTerm_MB (FloatMatrix &answer, Load *load, int side, TimeStep *tStep)
 Computes Lhs contribution due to applied Penetration bc.
 
virtual void computeOutFlowBCTerm_MB (FloatMatrix &answer, int side, TimeStep *tStep)
 Computes Lhs contribution due to outflow BC.
 
virtual void computeHomogenizedReinforceTerm_MB (FloatMatrix &answer, Load *load, TimeStep *tStep)
 
virtual void computeHomogenizedReinforceTerm_MC (FloatMatrix &answer, Load *load, TimeStep *tStep)
 
virtual double computeCriticalTimeStep (TimeStep *tStep)=0
 Computes the critical time increment.
 
void updateInternalState (TimeStep *tStep) override
 Updates element state after equilibrium in time step has been reached.
 
int checkConsistency () override
 Performs consistency check.
 
int giveInternalStateAtNode (FloatArray &answer, InternalStateType type, InternalStateMode mode, int node, TimeStep *tStep) override
 Returns internal state variable (like stress,strain) at node of element in Reduced form, the way how is obtained is dependent on InternalValueType.
 
virtual void giveLocalVelocityDofMap (IntArray &map)
 
virtual void giveLocalPressureDofMap (IntArray &map)
 
- 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)
 Constructor.
 
 Element (const Element &src)=delete
 
Elementoperator= (const Element &src)=delete
 
virtual ~Element ()
 Virtual destructor.
 
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.
 
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.
 
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 DofManagergiveInternalDofManager (int i) const
 Returns i-th internal element dof manager of the receiver.
 
virtual void setInternalDofManager (int num, std::unique_ptr< DofManager > dm)
 
virtual double giveCharacteristicValue (CharType type, TimeStep *tStep)
 Computes characteristic value of receiver of requested type in given time step.
 
virtual void computeLoadVector (FloatArray &answer, BodyLoad *load, CharType type, ValueModeType mode, TimeStep *tStep)
 Computes the contribution of the given body load (volumetric).
 
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.
 
virtual void computeTangentFromSurfaceLoad (FloatMatrix &answer, BoundaryLoad *load, int boundary, MatResponseMode rmode, TimeStep *tStep)
 Computes the tangent contribution of the given load at the given boundary.
 
virtual void computeTangentFromEdgeLoad (FloatMatrix &answer, BoundaryLoad *load, int boundary, MatResponseMode rmode, TimeStep *tStep)
 Computes the tangent contribution of the given load at the given boundary.
 
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.
 
const IntArraygiveBodyLoadList () const
 Returns receiver list of bodyloads.
 
const IntArraygiveBoundaryLoadList () const
 Returns receiver list of boundary loads.
 
void computeVectorOf (ValueModeType u, TimeStep *tStep, FloatArray &answer)
 Returns local vector of unknowns.
 
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.
 
void computeVectorOf (PrimaryField &field, const IntArray &dofIDMask, ValueModeType u, TimeStep *tStep, FloatArray &answer, bool padding=false)
 Returns local vector of unknowns.
 
void computeVectorOfPrescribed (ValueModeType u, TimeStep *tStep, FloatArray &answer)
 Returns local vector of prescribed unknowns.
 
void computeVectorOfPrescribed (const IntArray &dofIDMask, ValueModeType type, TimeStep *tStep, FloatArray &answer)
 Returns local vector of prescribed unknowns.
 
virtual int computeNumberOfDofs ()
 Computes or simply returns total number of element's local DOFs.
 
virtual int computeNumberOfGlobalDofs ()
 Computes the total number of element's global dofs.
 
int computeNumberOfPrimaryMasterDofs ()
 Computes the total number of element's primary master DOFs.
 
virtual bool computeGtoLRotationMatrix (FloatMatrix &answer)
 Returns transformation matrix from global c.s.
 
virtual bool giveRotationMatrix (FloatMatrix &answer)
 Transformation matrices updates rotation matrix between element-local and primary DOFs, taking into account nodal c.s.
 
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.
 
virtual void giveDofManDofIDMask (int inode, IntArray &answer) const
 Returns dofmanager dof mask for node.
 
virtual void giveInternalDofManDofIDMask (int inode, IntArray &answer) const
 Returns internal dofmanager dof mask for node.
 
virtual void giveElementDofIDMask (IntArray &answer) const
 Returns element dof mask for node.
 
virtual void computeField (ValueModeType mode, TimeStep *tStep, const FloatArray &lcoords, FloatArray &answer)
 Computes the unknown vector interpolated at the specified local coordinates.
 
virtual double computeVolumeAround (GaussPoint *gp)
 Returns volume related to given integration point.
 
virtual double computeVolumeAreaOrLength ()
 Computes the volume, area or length of the element depending on its spatial dimension.
 
double computeMeanSize ()
 Computes the size of the element defined as its length.
 
virtual double computeVolume ()
 Computes the volume.
 
virtual double computeArea ()
 Computes the area (zero for all but 2d geometries).
 
virtual double computeLength ()
 Computes the length (zero for all but 1D geometries)
 
virtual IntArray giveBoundaryEdgeNodes (int boundary, bool includeHierarchical=false) const
 Returns list of receiver boundary nodes for given edge.
 
virtual IntArray giveBoundarySurfaceNodes (int boundary, bool includeHierarchical=false) const
 Returns list of receiver boundary nodes for given surface.
 
virtual IntArray giveBoundaryNodes (int boundary) const
 Returns list of receiver boundary nodes for given boundary entity.
 
virtual std::unique_ptr< IntegrationRulegiveBoundaryEdgeIntegrationRule (int order, int boundary)
 Returns boundary edge integration rule.
 
virtual std::unique_ptr< IntegrationRulegiveBoundarySurfaceIntegrationRule (int order, int boundary)
 Returns boundary surface integration rule.
 
int giveDofManagerNumber (int i) const
 Translates local to global indices for dof managers.
 
const IntArraygiveDofManArray () const
 
void addDofManager (DofManager *dMan)
 
DofManagergiveDofManager (int i) const
 
NodegiveNode (int i) const
 Returns reference to the i-th node of element.
 
virtual ElementSidegiveSide (int i) const
 Returns reference to the i-th side of element.
 
virtual FEInterpolationgiveInterpolation () const
 
virtual FEInterpolationgiveInterpolation (DofIDItem id) const
 Returns the interpolation for the specific dof id.
 
virtual const FEInterpolationgetGeometryInterpolation () const
 Returns the geometry interpolation for the receiver.
 
virtual MaterialgiveMaterial ()
 
int giveMaterialNumber () const
 
CrossSectiongiveCrossSection ()
 
int getActivityTimeFunctionNumber ()
 Gets the activity time function number.
 
void setActivityTimeFunctionNumber (int funcIndx)
 Sets the activity time function of receiver.
 
void setMaterial (int matIndx)
 Sets the material of receiver.
 
virtual void setCrossSection (int csIndx)
 Sets the cross section model of receiver.
 
virtual int giveNumberOfDofManagers () const
 
void setNumberOfDofManagers (int i)
 Sets number of element dof managers.
 
virtual int giveNumberOfNodes () const
 Returns number of nodes of receiver.
 
void setDofManagers (const IntArray &dmans)
 Sets receiver dofManagers.
 
void setDofManager (int id, int dm)
 Sets receiver specific dofManager.
 
void setBodyLoads (const IntArray &bodyLoads)
 Sets receiver bodyLoadArray.
 
void setIntegrationRules (std ::vector< std ::unique_ptr< IntegrationRule > > irlist)
 Sets integration rules.
 
virtual integrationDomain giveIntegrationDomain () const
 Returns integration domain for receiver, used to initialize integration point over receiver volume.
 
virtual MaterialMode giveMaterialMode ()
 Returns material mode for receiver integration points.
 
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.
 
int giveRegionNumber ()
 
virtual void updateYourself (TimeStep *tStep)
 Updates element state after equilibrium in time step has been reached.
 
virtual void initializeYourself (TimeStep *timeStepWhenICApply)
 Initialization according to state given by initial conditions.
 
virtual bool isActivated (TimeStep *tStep)
 
virtual bool isCast (TimeStep *tStep)
 
virtual void initForNewStep ()
 Initializes receivers state to new time step.
 
virtual Element_Geometry_Type giveGeometryType () const =0
 Returns the element geometry type.
 
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 ()
 Returns the element spatial dimension (1, 2, or 3).
 
virtual int giveNumberOfBoundarySides ()
 Returns number of boundaries (entities of element_dimension-1: points, edges, surfaces)
 
virtual int giveNumberOfEdges () const
 
virtual int giveNumberOfSurfaces () const
 
virtual int giveDefaultIntegrationRule () const
 Returns id of default integration rule.
 
virtual IntegrationRulegiveDefaultIntegrationRulePtr ()
 Access method for default integration rule.
 
int giveNumberOfIntegrationRules ()
 
virtual IntegrationRulegiveIntegrationRule (int i)
 
std::vector< std ::unique_ptr< IntegrationRule > > & giveIntegrationRulesArray ()
 
virtual int testElementExtension (ElementExtension ext)
 Tests if the element implements required extension.
 
virtual int giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep)
 Returns the integration point corresponding value in full form.
 
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.
 
virtual double giveCharacteristicLength (const FloatArray &normalToCrackPlane)
 Returns the size of element in the given direction, in some cases adjusted (e.g.
 
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.
 
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.
 
virtual double giveCharacteristicSize (GaussPoint *gp, FloatArray &normalToCrackPlane, ElementCharSizeMethod method)
 Returns characteristic element size for a given integration point and given direction.
 
virtual double giveParentElSize () const
 Returns the size (length, area or volume depending on element type) of the parent element.
 
virtual void updateBeforeNonlocalAverage (TimeStep *tStep)
 Updates internal element state (in all integration points of receiver) before nonlocal averaging takes place.
 
virtual int computeGlobalCoordinates (FloatArray &answer, const FloatArray &lcoords)
 Computes the global coordinates from given element's local coordinates.
 
virtual bool computeLocalCoordinates (FloatArray &answer, const FloatArray &gcoords)
 Computes the element local coordinates from given global coordinates.
 
virtual int giveLocalCoordinateSystem (FloatMatrix &answer)
 Returns local coordinate system of receiver Required by material models with ortho- and anisotrophy.
 
virtual void giveLocalCoordinateSystemVector (InternalStateType isttype, FloatArray &answer)
 Returns unit vector of local coordinate system of receiver If local system is equal to global one, returns unit vectors.
 
virtual void computeMidPlaneNormal (FloatArray &answer, const GaussPoint *gp)
 Computes mid-plane normal of receiver at integration point.
 
virtual int adaptiveMap (Domain *oldd, TimeStep *tStep)
 Initializes the internal state variables stored in all IPs according to state in given domain.
 
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.
 
virtual int adaptiveUpdate (TimeStep *tStep)
 Updates the internal state variables stored in all IPs according to already mapped state.
 
virtual int adaptiveFinish (TimeStep *tStep)
 Finishes the mapping for given time step.
 
void updateLocalNumbering (EntityRenumberingFunctor &f) override
 Local renumbering support.
 
template<class T >
void ipEvaluator (T *src, void(T ::*f)(GaussPoint *gp))
 Integration point evaluator, loops over receiver IP's and calls given function (passed as f parameter) on them. The IP is parameter to function f.
 
template<class T , class S >
void ipEvaluator (T *src, void(T ::*f)(GaussPoint *, S &), S &_val)
 Integration point evaluator, loops over receiver IP's and calls given function (passed as f parameter) on them. The IP is parameter to function f as well as additional array.
 
virtual void drawYourself (oofegGraphicContext &gc, TimeStep *tStep)
 
virtual void drawAnnotation (oofegGraphicContext &gc, TimeStep *tStep)
 
virtual void drawRawGeometry (oofegGraphicContext &gc, TimeStep *tStep)
 
virtual void drawDeformedGeometry (oofegGraphicContext &gc, TimeStep *tStep, UnknownType)
 
virtual void drawScalar (oofegGraphicContext &gc, TimeStep *tStep)
 
virtual void drawSpecial (oofegGraphicContext &gc, TimeStep *tStep)
 
virtual void giveLocalIntVarMaxMin (oofegGraphicContext &gc, TimeStep *tStep, double &emin, double &emax)
 
virtual int giveInternalStateAtSide (FloatArray &answer, InternalStateType type, InternalStateMode mode, int side, TimeStep *tStep)
 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.
 
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)
 Sets receiver globally unique number.
 
elementParallelMode giveParallelMode () const
 Return elementParallelMode of receiver.
 
void setParallelMode (elementParallelMode _mode)
 Sets parallel mode of element.
 
virtual elementParallelMode giveKnotSpanParallelMode (int) const
 Returns the parallel mode for particular knot span of the receiver.
 
int packUnknowns (DataStream &buff, TimeStep *tStep)
 Pack all necessary data of element (according to its parallel_mode) integration points into given communication buffer.
 
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.
 
int estimatePackSize (DataStream &buff)
 Estimates the necessary pack size to hold all packed data of receiver.
 
const IntArraygivePartitionList () const
 Returns partition list of receiver.
 
void setPartitionList (IntArray &pl)
 Sets partition list of receiver.
 
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.
 
virtual double giveRelativeSelfComputationalCost ()
 Returns the weight representing relative computational cost of receiver The reference element is triangular plane stress element.
 
virtual double predictRelativeRedistributionCost ()
 Returns the relative redistribution cost of the receiver.
 
IntArraygiveBodyLoadArray ()
 Returns array containing load numbers of loads acting on element.
 
IntArraygiveBoundaryLoadArray ()
 Returns array containing load numbers of boundary loads acting on element.
 
void initializeFinish () override
 Finishes the initialization.
 
void postInitialize () override
 Performs post initialization steps.
 
void saveContext (DataStream &stream, ContextMode mode) override
 Stores receiver state to output stream.
 
void restoreContext (DataStream &stream, ContextMode mode) override
 Restores the receiver state previously written in stream.
 
const char * giveClassName () const override
 
void printOutputAt (FILE *file, TimeStep *tStep) override
 Prints output of receiver to stream, for given time step.
 
virtual const IntArray giveLocation ()
 
virtual void recalculateCoordinates (int nodeNumber, FloatArray &coords)
 
void setSharedEdgeID (int iedge, int globalID)
 Set global edge ID corresponding to given local edge id.
 
void setSharedSurfaceID (int isurf, int globalID)
 
const IntArraygiveSharedEdgeIDs () const
 
const IntArraygiveSharedSurfaceIDs () const
 
- Public Member Functions inherited from oofem::FEMComponent
 FEMComponent (int n, Domain *d)
 Regular constructor, creates component with given number and belonging to given domain.
 
virtual ~FEMComponent ()=default
 Virtual destructor.
 
virtual const char * giveInputRecordName () const =0
 
DomaingiveDomain () const
 
virtual void setDomain (Domain *d)
 Sets associated Domain.
 
int giveNumber () const
 
void setNumber (int num)
 Sets number of receiver.
 
virtual void initializeFrom (InputRecord &ir)
 Initializes receiver according to object description stored in input record.
 
virtual void printYourself ()
 Prints receiver state on stdout. Useful for debugging.
 
virtual InterfacegiveInterface (InterfaceType t)
 Interface requesting service.
 
std::string errorInfo (const char *func) const
 Returns string for prepending output (used by error reporting macros).
 

Protected Member Functions

virtual void computeDeviatoricStrain (FloatArray &answer, GaussPoint *gp, TimeStep *tStep)=0
 
virtual void computeDeviatoricStress (FloatArray &answer, const FloatArray &eps, GaussPoint *gp, TimeStep *tStep)=0
 
virtual void computeTangent (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep)=0
 
- Protected Member Functions inherited from oofem::Element
virtual void computeGaussPoints ()
 Initializes the array of integration rules member variable.
 

Protected Attributes

IntArray boundarySides
 Array of boundary sides.
 
IntArray boundaryCodes
 Boundary sides codes.
 
double t_supg = 0.
 Stabilization coefficients, updated for each solution step in updateStabilizationCoeffs()
 
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
 Array containing indexes of loads (body loads and boundary loads are kept separately), that apply on receiver.
 
IntArray boundaryLoadArray
 
std::vector< std ::unique_ptr< IntegrationRule > > integrationRulesArray
 List of integration rules of receiver (each integration rule contains associated integration points also).
 
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
 In parallel mode, globalNumber contains globally unique DoFManager number.
 
int numberOfGaussPoints
 Number of integration points as specified by nip.
 
elementParallelMode parallel_mode
 Determines the parallel mode of the element.
 
IntArray partitions
 List of partition sharing the shared element or remote partition containing remote element counterpart.
 
IntArray globalEdgeIDs
 List of global element edge IDs.
 
IntArray globalSurfaceIDs
 List of global element surface IDs.
 
- Protected Attributes inherited from oofem::FEMComponent
int number
 Component number.
 
Domaindomain
 Link to domain object, useful for communicating with other FEM components.
 

Static Protected Attributes

static ParamKey IPK_SUPGElement_bsides
 
static ParamKey IPK_SUPGElement_bcodes
 

Additional Inherited Members

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

Detailed Description

General stabilized SUPG/PSPG element for CFD analysis.

Definition at line 53 of file supgelement.h.

Constructor & Destructor Documentation

◆ SUPGElement()

oofem::SUPGElement::SUPGElement ( int  n,
Domain aDomain 
)

Definition at line 61 of file supgelement.C.

Member Function Documentation

◆ checkConsistency()

int oofem::SUPGElement::checkConsistency ( void  )
overridevirtual

Performs consistency check.

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

Returns
Zero value if check fail, otherwise nonzero.

Reimplemented from oofem::Element.

Reimplemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, and oofem::TR21_2D_SUPG.

Definition at line 327 of file supgelement.C.

◆ computeAccelerationTerm_MB()

virtual void oofem::SUPGElement::computeAccelerationTerm_MB ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual

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

Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.

Referenced by oofem::SUPGTangentAssembler::matrixFromElement(), and oofem::SUPGInternalForceAssembler::vectorFromElement().

◆ computeAccelerationTerm_MC()

virtual void oofem::SUPGElement::computeAccelerationTerm_MC ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual

◆ computeAdvectionDerivativeTerm_MB()

virtual void oofem::SUPGElement::computeAdvectionDerivativeTerm_MB ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual

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

Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.

Referenced by oofem::SUPGTangentAssembler::matrixFromElement().

◆ computeAdvectionDerivativeTerm_MC()

virtual void oofem::SUPGElement::computeAdvectionDerivativeTerm_MC ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual

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

Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.

Referenced by oofem::SUPGTangentAssembler::matrixFromElement().

◆ computeAdvectionTerm_MB()

virtual void oofem::SUPGElement::computeAdvectionTerm_MB ( FloatArray answer,
TimeStep tStep 
)
pure virtual

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

Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.

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

◆ computeAdvectionTerm_MC()

virtual void oofem::SUPGElement::computeAdvectionTerm_MC ( FloatArray answer,
TimeStep tStep 
)
pure virtual

◆ computeBCLhsPressureTerm_MB()

void oofem::SUPGElement::computeBCLhsPressureTerm_MB ( FloatMatrix answer,
TimeStep tStep 
)
virtual

◆ computeBCLhsPressureTerm_MC()

void oofem::SUPGElement::computeBCLhsPressureTerm_MC ( FloatMatrix answer,
TimeStep tStep 
)
virtual

◆ computeBCLhsTerm_MB()

◆ computeBCRhsTerm_MB()

virtual void oofem::SUPGElement::computeBCRhsTerm_MB ( FloatArray answer,
TimeStep tStep 
)
pure virtual

Computes Rhs terms due to boundary conditions.

Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.

◆ computeBCRhsTerm_MC()

virtual void oofem::SUPGElement::computeBCRhsTerm_MC ( FloatArray answer,
TimeStep tStep 
)
pure virtual

Computes Rhs terms due to boundary conditions.

Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.

◆ computeCriticalTimeStep()

virtual double oofem::SUPGElement::computeCriticalTimeStep ( TimeStep tStep)
pure virtual

◆ computeDeviatoricStrain()

virtual void oofem::SUPGElement::computeDeviatoricStrain ( FloatArray answer,
GaussPoint gp,
TimeStep tStep 
)
protectedpure virtual

◆ computeDeviatoricStress()

virtual void oofem::SUPGElement::computeDeviatoricStress ( FloatArray answer,
const FloatArray eps,
GaussPoint gp,
TimeStep tStep 
)
protectedpure virtual

◆ computeDiffusionDerivativeTerm_MB()

virtual void oofem::SUPGElement::computeDiffusionDerivativeTerm_MB ( FloatMatrix answer,
MatResponseMode  mode,
TimeStep tStep 
)
pure virtual

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

Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.

Referenced by oofem::SUPGTangentAssembler::matrixFromElement().

◆ computeDiffusionDerivativeTerm_MC()

virtual void oofem::SUPGElement::computeDiffusionDerivativeTerm_MC ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual

Computes diffusion derivative terms for mass conservation equation.

Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.

Referenced by oofem::SUPGTangentAssembler::matrixFromElement().

◆ computeDiffusionTerm_MB()

virtual void oofem::SUPGElement::computeDiffusionTerm_MB ( FloatArray answer,
TimeStep tStep 
)
pure virtual

Computes diffusion terms for momentum balance equations(s).

Implemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR1_2D_SUPG_AXI.

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

◆ computeDiffusionTerm_MC()

virtual void oofem::SUPGElement::computeDiffusionTerm_MC ( FloatArray answer,
TimeStep tStep 
)
pure virtual

◆ computeHomogenizedReinforceTerm_MB()

virtual void oofem::SUPGElement::computeHomogenizedReinforceTerm_MB ( FloatMatrix answer,
Load load,
TimeStep tStep 
)
inlinevirtual

Reimplemented in oofem::TR1_2D_SUPG.

Definition at line 183 of file supgelement.h.

References oofem::FloatMatrix::clear(), and OOFEM_WARNING.

◆ computeHomogenizedReinforceTerm_MC()

virtual void oofem::SUPGElement::computeHomogenizedReinforceTerm_MC ( FloatMatrix answer,
Load load,
TimeStep tStep 
)
inlinevirtual

Reimplemented in oofem::TR1_2D_SUPG.

Definition at line 187 of file supgelement.h.

References oofem::FloatMatrix::clear(), and OOFEM_WARNING.

◆ computeLinearAdvectionTerm_MC()

virtual void oofem::SUPGElement::computeLinearAdvectionTerm_MC ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual

◆ computeLSICStabilizationTerm_MB()

virtual void oofem::SUPGElement::computeLSICStabilizationTerm_MB ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual

◆ computeOutFlowBCTerm_MB()

virtual void oofem::SUPGElement::computeOutFlowBCTerm_MB ( FloatMatrix answer,
int  side,
TimeStep tStep 
)
inlinevirtual

Computes Lhs contribution due to outflow BC.

Reimplemented in oofem::TR1_2D_SUPG, and oofem::TR1_2D_SUPG_AXI.

Definition at line 177 of file supgelement.h.

References oofem::FloatMatrix::clear(), and OOFEM_WARNING.

◆ computePenetrationWithResistanceBCTerm_MB()

virtual void oofem::SUPGElement::computePenetrationWithResistanceBCTerm_MB ( FloatMatrix answer,
Load load,
int  side,
TimeStep tStep 
)
inlinevirtual

Computes Lhs contribution due to applied Penetration bc.

Reimplemented in oofem::TR1_2D_SUPG, and oofem::TR1_2D_SUPG_AXI.

Definition at line 170 of file supgelement.h.

References oofem::FloatMatrix::clear(), and OOFEM_WARNING.

◆ computePressureTerm_MB()

virtual void oofem::SUPGElement::computePressureTerm_MB ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual

◆ computePressureTerm_MC()

virtual void oofem::SUPGElement::computePressureTerm_MC ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual

◆ computeSlipWithFrictionBCTerm_MB()

virtual void oofem::SUPGElement::computeSlipWithFrictionBCTerm_MB ( FloatMatrix answer,
Load load,
int  side,
TimeStep tStep 
)
inlinevirtual

Computes Lhs term due to applied slip with friction bc.

Reimplemented in oofem::TR1_2D_SUPG, and oofem::TR1_2D_SUPG_AXI.

Definition at line 163 of file supgelement.h.

References oofem::FloatMatrix::clear(), and OOFEM_WARNING.

◆ computeTangent()

virtual void oofem::SUPGElement::computeTangent ( FloatMatrix answer,
MatResponseMode  mode,
GaussPoint gp,
TimeStep tStep 
)
protectedpure virtual

◆ giveCharacteristicMatrix()

void oofem::SUPGElement::giveCharacteristicMatrix ( FloatMatrix answer,
CharType  type,
TimeStep tStep 
)
overridevirtual

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

Parameters
answerRequested characteristic matrix (stiffness, tangent, ...). If element has no capability to compute requested type of characteristic matrix error function is invoked.
typeId of characteristic component requested.
tStepTime step when answer is computed.

Reimplemented from oofem::Element.

Reimplemented in oofem::SUPGElement2.

Definition at line 95 of file supgelement.C.

References oofem::FloatMatrix::assemble(), oofem::FloatMatrix::isNotEmpty(), OOFEM_ERROR, oofem::FloatMatrix::resize(), oofem::FloatMatrix::times(), oofem::VST_Density, oofem::VST_Velocity, and oofem::FloatMatrix::zero().

◆ giveCharacteristicVector()

void oofem::SUPGElement::giveCharacteristicVector ( FloatArray answer,
CharType  type,
ValueModeType  mode,
TimeStep tStep 
)
overridevirtual

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

If element has no capability to compute requested type of characteristic vector error function is invoked.

Parameters
answerRequested characteristic vector.
typeId of characteristic component requested.
modeDetermines mode of answer.
tStepTime step when answer is computed.

Reimplemented from oofem::Element.

Reimplemented in oofem::SUPGElement2.

Definition at line 155 of file supgelement.C.

References oofem::FloatArray::assemble(), oofem::FloatArray::beProductOf(), oofem::FloatArray::clear(), OOFEM_ERROR, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

◆ giveInputRecord()

void oofem::SUPGElement::giveInputRecord ( DynamicInputRecord input)
overridevirtual

Setups the input record string of receiver.

Parameters
inputDynamic input record to be filled by receiver.

Reimplemented from oofem::Element.

Reimplemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, and oofem::TR1_2D_SUPG2_AXI.

Definition at line 84 of file supgelement.C.

References oofem::DynamicInputRecord::setField().

◆ giveInternalStateAtNode()

int oofem::SUPGElement::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
answerContains result, zero sized if not supported.
typeDetermines the internal variable requested (physical meaning).
modeDetermines the mode of variable (recovered, local, ...).
nodeNode number, for which variable is required.
tStepTime step.
Returns
Nonzero if o.k, zero otherwise.

Reimplemented from oofem::Element.

Reimplemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, and oofem::TR21_2D_SUPG.

Definition at line 355 of file supgelement.C.

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

◆ giveLocalPressureDofMap()

virtual void oofem::SUPGElement::giveLocalPressureDofMap ( IntArray map)
inlinevirtual

◆ giveLocalVelocityDofMap()

virtual void oofem::SUPGElement::giveLocalVelocityDofMap ( IntArray map)
inlinevirtual

◆ initializeFrom()

void oofem::SUPGElement::initializeFrom ( InputRecord ir,
int  priority 
)
overridevirtual

Reimplemented from oofem::Element.

Reimplemented in oofem::SUPGElement2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, and oofem::TR1_2D_SUPG2_AXI.

Definition at line 67 of file supgelement.C.

References OOFEM_ERROR, and PM_UPDATE_PARAMETER_AND_REPORT.

◆ updateElementForNewInterfacePosition()

virtual void oofem::SUPGElement::updateElementForNewInterfacePosition ( TimeStep tStep)
inlinevirtual

◆ updateInternalState()

void oofem::SUPGElement::updateInternalState ( TimeStep tStep)
overridevirtual

Updates element state after equilibrium in time step has been reached.

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

Parameters
tStepTime step for newly reached state.
See also
Material::updateYourself
IntegrationRule::updateYourself
GaussPoint::updateYourself
Element::updateInternalState

Reimplemented from oofem::Element.

Reimplemented in oofem::SUPGElement2.

Definition at line 339 of file supgelement.C.

Referenced by oofem::SUPG::applyIC().

◆ updateStabilizationCoeffs()

void oofem::SUPGElement::updateStabilizationCoeffs ( TimeStep tStep)
inlineoverridevirtual

Updates the stabilization coefficients used for CBS and SUPG algorithms.

Parameters
tStepActive time step.

Reimplemented from oofem::FMElement.

Reimplemented in oofem::Tet1_3D_SUPG, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR1_2D_SUPG2_AXI, oofem::TR1_2D_SUPG_AXI, and oofem::TR21_2D_SUPG.

Definition at line 78 of file supgelement.h.

Member Data Documentation

◆ boundaryCodes

IntArray oofem::SUPGElement::boundaryCodes
protected

Boundary sides codes.

Definition at line 59 of file supgelement.h.

◆ boundarySides

IntArray oofem::SUPGElement::boundarySides
protected

Array of boundary sides.

Definition at line 57 of file supgelement.h.

◆ IPK_SUPGElement_bcodes

ParamKey oofem::SUPGElement::IPK_SUPGElement_bcodes
staticprotected

Definition at line 68 of file supgelement.h.

◆ IPK_SUPGElement_bsides

ParamKey oofem::SUPGElement::IPK_SUPGElement_bsides
staticprotected

Definition at line 67 of file supgelement.h.

◆ t_lsic

double oofem::SUPGElement::t_lsic = 0.
protected

Definition at line 65 of file supgelement.h.

◆ t_pspg

double oofem::SUPGElement::t_pspg = 0.
protected

Definition at line 65 of file supgelement.h.

◆ t_supg

double oofem::SUPGElement::t_supg = 0.
protected

Stabilization coefficients, updated for each solution step in updateStabilizationCoeffs()

Definition at line 65 of file supgelement.h.


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

This page is part of the OOFEM documentation. Copyright (c) 2011 Borek Patzak
Project e-mail: info@oofem.org
Generated at Tue Dec 23 2025 23:54:19 for OOFEM by doxygen 1.9.8 written by Dimitri van Heesch, © 1997-2011