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)
 
virtual ~SUPGElement ()
 
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 giveCharacteristicMatrix (FloatMatrix &answer, CharType type, TimeStep *tStep)
 Computes characteristic matrix of receiver of requested type in given time step. More...
 
virtual void giveCharacteristicVector (FloatArray &answer, CharType type, ValueModeType mode, TimeStep *tStep)
 Computes characteristic vector of receiver of requested type in given time step. More...
 
virtual void updateStabilizationCoeffs (TimeStep *tStep)
 Updates the stabilization coefficients used for CBS and SUPG algorithms. More...
 
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). More...
 
virtual void computeAdvectionTerm_MB (FloatArray &answer, TimeStep *tStep)=0
 Computes nonlinear advection terms for momentum balance equations(s). More...
 
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. More...
 
virtual void computeDiffusionTerm_MB (FloatArray &answer, TimeStep *tStep)=0
 Computes diffusion terms for momentum balance equations(s). More...
 
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. More...
 
virtual void computePressureTerm_MB (FloatMatrix &answer, TimeStep *tStep)=0
 Computes pressure terms for momentum balance equations(s). More...
 
virtual void computeBCLhsPressureTerm_MC (FloatMatrix &answer, TimeStep *tStep)
 Computes Lhs terms due to boundary conditions - pressure. More...
 
virtual void computeLSICStabilizationTerm_MB (FloatMatrix &answer, TimeStep *tStep)=0
 Computes SLIC stabilization term for momentum balance equation(s). More...
 
virtual void computeLinearAdvectionTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0
 Computes the linear advection term for mass conservation equation. More...
 
virtual void computeAdvectionTerm_MC (FloatArray &answer, TimeStep *tStep)=0
 Computes advection terms for mass conservation equation. More...
 
virtual void computeAdvectionDerivativeTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0
 Computes the derivative of advection terms for mass conservation equation with respect to nodal velocities. More...
 
virtual void computeDiffusionDerivativeTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0
 Computes diffusion derivative terms for mass conservation equation. More...
 
virtual void computeDiffusionTerm_MC (FloatArray &answer, TimeStep *tStep)=0
 Computes diffusion terms for mass conservation equation. More...
 
virtual void computeAccelerationTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0
 Computes acceleration terms for mass conservation equation. More...
 
virtual void computePressureTerm_MC (FloatMatrix &answer, TimeStep *tStep)=0
 Computes pressure terms for mass conservation equation. 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 computeBCRhsTerm_MB (FloatArray &answer, TimeStep *tStep)=0
 Computes Rhs terms due to boundary conditions. More...
 
virtual void computeBCRhsTerm_MC (FloatArray &answer, TimeStep *tStep)=0
 Computes Rhs terms due to boundary conditions. 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)
 
virtual double computeCriticalTimeStep (TimeStep *tStep)=0
 Computes the critical time increment. More...
 
virtual void updateInternalState (TimeStep *tStep)
 Updates element state after equilibrium in time step has been reached. More...
 
virtual int checkConsistency ()
 Performs consistency check. 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 giveLocalVelocityDofMap (IntArray &map)
 
virtual void giveLocalPressureDofMap (IntArray &map)
 
- 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
 
Elementoperator= (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 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. 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 IntArraygivePartitionList () 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...
 
IntArraygiveBodyLoadArray ()
 Returns array containing load numbers of loads acting on element. More...
 
IntArraygiveBoundaryLoadArray ()
 Returns array containing load numbers of boundary loads acting on element. 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 void printOutputAt (FILE *file, TimeStep *tStep)
 Prints output of receiver to stream, for given time step. More...
 
virtual const char * giveClassName () const
 
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 int giveNumberOfInternalDofManagers () const
 
virtual DofManagergiveInternalDofManager (int i) const
 Returns i-th internal element dof manager of the receiver. More...
 
virtual double giveCharacteristicValue (CharType type, TimeStep *tStep)
 Computes characteristic value of receiver of requested type in given time step. 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 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 IntArraygiveBodyLoadList () const
 Returns receiver list of bodyloads. More...
 
const IntArraygiveBoundaryLoadList () 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 computeNumberOfDofs ()
 Computes or simply returns total number of element's local DOFs. 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 giveDofManDofIDMask (int inode, IntArray &answer) const
 Returns dofmanager dof mask for node. More...
 
virtual void giveInternalDofManDofIDMask (int inode, IntArray &answer) const
 Returns internal dofmanager dof mask for node. 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 computeVolumeAround (GaussPoint *gp)
 Returns volume related to given integration point. 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 IntegrationRulegiveBoundaryEdgeIntegrationRule (int order, int boundary)
 Returns boundary edge integration rule. More...
 
virtual IntegrationRulegiveBoundarySurfaceIntegrationRule (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 IntArraygiveDofManArray () const
 
void addDofManager (DofManager *dMan)
 
DofManagergiveDofManager (int i) const
 
NodegiveNode (int i) const
 Returns reference to the i-th node of element. More...
 
virtual ElementSidegiveSide (int i) const
 Returns reference to the i-th side of element. More...
 
virtual FEInterpolationgiveInterpolation () const
 
virtual FEInterpolationgiveInterpolation (DofIDItem id) const
 Returns the interpolation for the specific dof id. More...
 
virtual MaterialgiveMaterial ()
 
int giveMaterialNumber () const
 
CrossSectiongiveCrossSection ()
 
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 MaterialMode giveMaterialMode ()
 Returns material mode for receiver integration points. 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 updateYourself (TimeStep *tStep)
 Updates element state after equilibrium in time step has been reached. 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 IntegrationRulegiveDefaultIntegrationRulePtr ()
 Access method for default integration rule. More...
 
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. More...
 
virtual int giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep)
 Returns the integration point corresponding value in full form. 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...
 
virtual const char * giveInputRecordName () const =0
 
DomaingiveDomain () 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...
 
virtual InterfacegiveInterface (InterfaceType t)
 Interface requesting service. More...
 
std::string errorInfo (const char *func) const
 Returns string for prepending output (used by error reporting macros). More...
 

Protected Member Functions

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

Protected Attributes

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...
 
Domaindomain
 Link to domain object, useful for communicating with other FEM components. More...
 

Detailed Description

General stabilized SUPG/PSPG element for CFD analysis.

Definition at line 58 of file supgelement.h.

Constructor & Destructor Documentation

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

Definition at line 55 of file supgelement.C.

oofem::SUPGElement::~SUPGElement ( )
virtual

Definition at line 60 of file supgelement.C.

Member Function Documentation

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

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

Definition at line 339 of file supgelement.C.

Referenced by oofem::TR21_2D_SUPG::checkConsistency(), oofem::Quad10_2D_SUPG::checkConsistency(), oofem::TR1_2D_SUPG::checkConsistency(), and computeHomogenizedReinforceTerm_MC().

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::TR1_2D_SUPG2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2_AXI, oofem::SUPGElement2, and oofem::TR1_2D_SUPG_AXI.

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

virtual void oofem::SUPGElement::computeAccelerationTerm_MC ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual
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::TR1_2D_SUPG2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2_AXI, oofem::SUPGElement2, and oofem::TR1_2D_SUPG_AXI.

Referenced by oofem::SUPGTangentAssembler::matrixFromElement(), and updateElementForNewInterfacePosition().

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::TR1_2D_SUPG2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2_AXI, oofem::SUPGElement2, and oofem::TR1_2D_SUPG_AXI.

Referenced by oofem::SUPGTangentAssembler::matrixFromElement(), and updateElementForNewInterfacePosition().

virtual void oofem::SUPGElement::computeAdvectionTerm_MB ( FloatArray answer,
TimeStep tStep 
)
pure virtual
virtual void oofem::SUPGElement::computeAdvectionTerm_MC ( FloatArray answer,
TimeStep tStep 
)
pure virtual
virtual void oofem::SUPGElement::computeBCRhsTerm_MB ( FloatArray answer,
TimeStep tStep 
)
pure virtual
virtual void oofem::SUPGElement::computeBCRhsTerm_MC ( FloatArray answer,
TimeStep tStep 
)
pure virtual
virtual double oofem::SUPGElement::computeCriticalTimeStep ( TimeStep tStep)
pure virtual
virtual void oofem::SUPGElement::computeDeviatoricStrain ( FloatArray answer,
GaussPoint gp,
TimeStep tStep 
)
protectedpure virtual
void oofem::SUPGElement::computeDeviatoricStress ( FloatArray answer,
GaussPoint gp,
TimeStep tStep 
)
protectedvirtual
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::TR1_2D_SUPG2, oofem::TR1_2D_SUPG, oofem::TR1_2D_SUPG2_AXI, oofem::SUPGElement2, and oofem::TR1_2D_SUPG_AXI.

Referenced by giveCharacteristicMatrix(), oofem::SUPGTangentAssembler::matrixFromElement(), and updateElementForNewInterfacePosition().

virtual void oofem::SUPGElement::computeDiffusionDerivativeTerm_MC ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual
virtual void oofem::SUPGElement::computeDiffusionTerm_MB ( FloatArray answer,
TimeStep tStep 
)
pure virtual
virtual void oofem::SUPGElement::computeDiffusionTerm_MC ( FloatArray answer,
TimeStep tStep 
)
pure virtual
virtual void oofem::SUPGElement::computeHomogenizedReinforceTerm_MB ( FloatMatrix answer,
Load load,
TimeStep tStep 
)
inlinevirtual

Reimplemented in oofem::TR1_2D_SUPG.

Definition at line 186 of file supgelement.h.

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

Referenced by computeBCLhsTerm_MB().

virtual void oofem::SUPGElement::computeHomogenizedReinforceTerm_MC ( FloatMatrix answer,
Load load,
TimeStep tStep 
)
inlinevirtual
virtual void oofem::SUPGElement::computeLinearAdvectionTerm_MC ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual
virtual void oofem::SUPGElement::computeLSICStabilizationTerm_MB ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual
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 180 of file supgelement.h.

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

Referenced by computeBCLhsPressureTerm_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 173 of file supgelement.h.

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

Referenced by computeBCLhsTerm_MB().

virtual void oofem::SUPGElement::computePressureTerm_MB ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual
virtual void oofem::SUPGElement::computePressureTerm_MC ( FloatMatrix answer,
TimeStep tStep 
)
pure virtual
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 166 of file supgelement.h.

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

Referenced by computeBCLhsTerm_MB().

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

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 91 of file supgelement.C.

References oofem::FloatMatrix::assemble(), computeBCLhsPressureTerm_MB(), computeBCLhsPressureTerm_MC(), computeBCLhsTerm_MB(), computeDiffusionDerivativeTerm_MB(), computeDiffusionDerivativeTerm_MC(), computeLinearAdvectionTerm_MC(), oofem::Element::computeNumberOfDofs(), computePressureTerm_MB(), computePressureTerm_MC(), oofem::FEMComponent::giveDomain(), oofem::Domain::giveEngngModel(), giveLocalPressureDofMap(), giveLocalVelocityDofMap(), oofem::EngngModel::giveVariableScale(), oofem::FloatMatrix::isNotEmpty(), OOFEM_ERROR, oofem::FloatMatrix::resize(), oofem::FloatMatrix::times(), oofem::VST_Density, oofem::VST_Velocity, and oofem::FloatMatrix::zero().

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

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 151 of file supgelement.C.

References oofem::FloatArray::assemble(), oofem::FloatArray::beProductOf(), oofem::FloatArray::clear(), computeAdvectionTerm_MC(), computeBCRhsTerm_MB(), computeBCRhsTerm_MC(), computeDiffusionTerm_MB(), computeDiffusionTerm_MC(), computeLinearAdvectionTerm_MC(), oofem::Element::computeNumberOfDofs(), computePressureTerm_MB(), computePressureTerm_MC(), oofem::FMElement::computeVectorOfPressures(), oofem::FMElement::computeVectorOfVelocities(), giveLocalPressureDofMap(), giveLocalVelocityDofMap(), OOFEM_ERROR, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

int oofem::SUPGElement::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
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::TR1_2D_SUPG, oofem::TR1_2D_SUPG2, oofem::TR21_2D_SUPG, oofem::TR1_2D_SUPG2_AXI, and oofem::SUPGElement2.

Definition at line 366 of file supgelement.C.

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

Referenced by computeHomogenizedReinforceTerm_MC(), oofem::TR1_2D_SUPG2_AXI::giveInternalStateAtNode(), oofem::TR21_2D_SUPG::giveInternalStateAtNode(), oofem::Quad10_2D_SUPG::giveInternalStateAtNode(), oofem::TR1_2D_SUPG2::giveInternalStateAtNode(), and oofem::TR1_2D_SUPG::giveInternalStateAtNode().

IRResultType oofem::SUPGElement::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
irInput record to initialize from.
Returns
IRResultType

Reimplemented from oofem::Element.

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

Definition at line 65 of file supgelement.C.

References _IFT_SUPGElement_bcodes, _IFT_SUPGElement_bsides, boundaryCodes, boundarySides, oofem::Element::initializeFrom(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, and oofem::IntArray::isEmpty().

Referenced by oofem::SUPGElement2::initializeFrom(), oofem::TR1_2D_SUPG2_AXI::initializeFrom(), oofem::TR1_2D_SUPG2::initializeFrom(), and oofem::TR1_2D_SUPG::initializeFrom().

void oofem::SUPGElement::updateInternalState ( 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
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 351 of file supgelement.C.

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

Referenced by oofem::SUPG::applyIC(), and computeHomogenizedReinforceTerm_MC().

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

Updates the stabilization coefficients used for CBS and SUPG algorithms.

Parameters
tStepActive time step.

Reimplemented from oofem::FMElement.

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

Definition at line 81 of file supgelement.h.

Member Data Documentation

IntArray oofem::SUPGElement::boundaryCodes
protected
IntArray oofem::SUPGElement::boundarySides
protected
double oofem::SUPGElement::t_pspg
protected

Definition at line 70 of file supgelement.h.

Referenced by oofem::TR1_2D_SUPG_AXI::computeAccelerationTerm_MC(), oofem::SUPGElement2::computeAccelerationTerm_MC(), oofem::TR1_2D_SUPG2_AXI::computeAccelerationTerm_MC(), oofem::TR1_2D_SUPG::computeAccelerationTerm_MC(), oofem::TR1_2D_SUPG_AXI::computeAdvectionDerivativeTerm_MC(), oofem::SUPGElement2::computeAdvectionDerivativeTerm_MC(), oofem::TR1_2D_SUPG2_AXI::computeAdvectionDerivativeTerm_MC(), oofem::TR1_2D_SUPG::computeAdvectionDerivativeTerm_MC(), oofem::TR1_2D_SUPG2::computeAdvectionDerivativeTerm_MC(), oofem::TR1_2D_SUPG_AXI::computeAdvectionTerm_MC(), oofem::SUPGElement2::computeAdvectionTerm_MC(), oofem::TR1_2D_SUPG2_AXI::computeAdvectionTerm_MC(), oofem::TR1_2D_SUPG::computeAdvectionTerm_MC(), oofem::TR1_2D_SUPG2::computeAdvectionTerm_MC(), oofem::TR1_2D_SUPG_AXI::computeBCRhsTerm_MC(), oofem::SUPGElement2::computeBCRhsTerm_MC(), oofem::TR1_2D_SUPG2_AXI::computeBCRhsTerm_MC(), oofem::TR1_2D_SUPG::computeBCRhsTerm_MC(), oofem::TR1_2D_SUPG_AXI::computeDiffusionDerivativeTerm_MC(), oofem::SUPGElement2::computeDiffusionDerivativeTerm_MC(), oofem::TR1_2D_SUPG2_AXI::computeDiffusionDerivativeTerm_MC(), oofem::TR1_2D_SUPG_AXI::computeDiffusionTerm_MC(), oofem::TR1_2D_SUPG2_AXI::computeDiffusionTerm_MC(), oofem::TR1_2D_SUPG::computeHomogenizedReinforceTerm_MC(), oofem::TR1_2D_SUPG_AXI::computeLoadVector(), oofem::SUPGElement2::computeLoadVector(), oofem::TR1_2D_SUPG::computeLoadVector(), oofem::TR1_2D_SUPG_AXI::computePressureTerm_MC(), oofem::SUPGElement2::computePressureTerm_MC(), oofem::TR1_2D_SUPG2_AXI::computePressureTerm_MC(), oofem::TR1_2D_SUPG::computePressureTerm_MC(), oofem::TR1_2D_SUPG2::computePressureTerm_MC(), oofem::Tet1_3D_SUPG::updateStabilizationCoeffs(), oofem::TR1_2D_SUPG_AXI::updateStabilizationCoeffs(), oofem::TR1_2D_SUPG2_AXI::updateStabilizationCoeffs(), oofem::TR1_2D_SUPG2::updateStabilizationCoeffs(), oofem::TR1_2D_SUPG::updateStabilizationCoeffs(), oofem::TR21_2D_SUPG::updateStabilizationCoeffs(), and oofem::Quad10_2D_SUPG::updateStabilizationCoeffs().

double oofem::SUPGElement::t_supg
protected

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

Definition at line 70 of file supgelement.h.

Referenced by oofem::TR1_2D_SUPG_AXI::computeAccelerationTerm_MB(), oofem::SUPGElement2::computeAccelerationTerm_MB(), oofem::TR1_2D_SUPG2_AXI::computeAccelerationTerm_MB(), oofem::TR1_2D_SUPG::computeAccelerationTerm_MB(), oofem::TR1_2D_SUPG2::computeAccelerationTerm_MB(), oofem::TR1_2D_SUPG_AXI::computeAdvectionDerivativeTerm_MB(), oofem::SUPGElement2::computeAdvectionDerivativeTerm_MB(), oofem::TR1_2D_SUPG2_AXI::computeAdvectionDerivativeTerm_MB(), oofem::TR1_2D_SUPG::computeAdvectionDerivativeTerm_MB(), oofem::TR1_2D_SUPG2::computeAdvectionDerivativeTerm_MB(), oofem::TR1_2D_SUPG_AXI::computeAdvectionTerm_MB(), oofem::SUPGElement2::computeAdvectionTerm_MB(), oofem::TR1_2D_SUPG2_AXI::computeAdvectionTerm_MB(), oofem::TR1_2D_SUPG::computeAdvectionTerm_MB(), oofem::TR1_2D_SUPG2::computeAdvectionTerm_MB(), oofem::TR1_2D_SUPG_AXI::computeBCRhsTerm_MB(), oofem::SUPGElement2::computeBCRhsTerm_MB(), oofem::TR1_2D_SUPG2_AXI::computeBCRhsTerm_MB(), oofem::TR1_2D_SUPG::computeBCRhsTerm_MB(), oofem::TR1_2D_SUPG2::computeBCRhsTerm_MB(), oofem::TR1_2D_SUPG_AXI::computeDiffusionDerivativeTerm_MB(), oofem::SUPGElement2::computeDiffusionDerivativeTerm_MB(), oofem::TR1_2D_SUPG2_AXI::computeDiffusionDerivativeTerm_MB(), oofem::TR1_2D_SUPG_AXI::computeDiffusionTerm_MB(), oofem::SUPGElement2::computeDiffusionTerm_MB(), oofem::TR1_2D_SUPG2_AXI::computeDiffusionTerm_MB(), oofem::TR1_2D_SUPG::computeHomogenizedReinforceTerm_MB(), oofem::TR1_2D_SUPG_AXI::computeLoadVector(), oofem::SUPGElement2::computeLoadVector(), oofem::TR1_2D_SUPG::computeLoadVector(), oofem::TR1_2D_SUPG_AXI::computePressureTerm_MB(), oofem::SUPGElement2::computePressureTerm_MB(), oofem::TR1_2D_SUPG2_AXI::computePressureTerm_MB(), oofem::TR1_2D_SUPG::computePressureTerm_MB(), oofem::TR1_2D_SUPG2::computePressureTerm_MB(), oofem::Tet1_3D_SUPG::updateStabilizationCoeffs(), oofem::TR1_2D_SUPG_AXI::updateStabilizationCoeffs(), oofem::TR1_2D_SUPG2_AXI::updateStabilizationCoeffs(), oofem::TR1_2D_SUPG2::updateStabilizationCoeffs(), oofem::TR1_2D_SUPG::updateStabilizationCoeffs(), oofem::TR21_2D_SUPG::updateStabilizationCoeffs(), and oofem::Quad10_2D_SUPG::updateStabilizationCoeffs().


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 Jan 2 2018 20:07:42 for OOFEM by doxygen 1.8.11 written by Dimitri van Heesch, © 1997-2011