This class represents transient incompressible flow problem. More...

#include <supg.h>

Inheritance diagram for oofem::SUPG:

Collaboration diagram for oofem::SUPG:

Public Member Functions
	SUPG (int i, EngngModel *_master=NULL)

virtual	~SUPG ()

virtual void	solveYourselfAt (TimeStep *tStep)
	Solves problem for given time step. More...

virtual void	updateYourself (TimeStep *tStep)
	Updates internal state after finishing time step. More...

virtual double	giveUnknownComponent (ValueModeType mode, TimeStep tStep, Domain d, Dof *dof)
	Returns requested unknown. More...

virtual void	updateComponent (TimeStep tStep, NumericalCmpn cmpn, Domain d)
	Updates components mapped to numerical method if necessary during solution process. More...

virtual double	giveReynoldsNumber ()

virtual contextIOResultType	saveContext (DataStream &stream, ContextMode mode)
	Stores the state of model to output stream. More...

virtual contextIOResultType	restoreContext (DataStream &stream, ContextMode mode)
	Restores the state of model from output stream. More...

virtual void	updateDomainLinks ()
	Updates domain links after the domains of receiver have changed. More...

virtual TimeStep *	giveNextStep ()
	Returns next time step (next to current step) of receiver. More...

virtual TimeStep *	giveSolutionStepWhenIcApply (bool force=false)
	Returns the solution step when Initial Conditions (IC) apply. More...

virtual NumericalMethod *	giveNumericalMethod (MetaStep *mStep)
	Returns reference to receiver's numerical method. More...

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

virtual int	checkConsistency ()
	Allows programmer to test some receiver's internal data, before computation begins. More...

virtual const char *	giveClassName () const
	Returns class name of the receiver. More...

virtual const char *	giveInputRecordName () const

virtual fMode	giveFormulation ()
	Indicates type of non linear computation (total or updated formulation). More...

virtual void	printDofOutputAt (FILE stream, Dof iDof, TimeStep *tStep)
	DOF printing routine. More...

virtual int	requiresUnknownsDictionaryUpdate ()
	Indicates if EngngModel requires Dofs dictionaries to be updated. More...

virtual bool	giveEquationScalingFlag ()
	Returns the Equation scaling flag, which is used to indicate that governing equation(s) are scaled, or non-dimensionalized. More...

virtual double	giveVariableScale (VarScaleType varId)
	Returns the scale factor for given variable type. More...

virtual void	updateDofUnknownsDictionary (DofManager dman, TimeStep tStep)
	Updates necessary values in Dofs unknown dictionaries. More...

virtual int	giveUnknownDictHashIndx (ValueModeType mode, TimeStep *tStep)
	This method is responsible for computing unique dictionary id (ie hash value) from given valueModeType and time step. More...

virtual MaterialInterface *	giveMaterialInterface (int n)
	Returns material interface representation for given domain. More...

Public Member Functions inherited from oofem::FluidModel
	FluidModel (int i, EngngModel *master)

virtual int	forceEquationNumbering (int id)
	Forces equation renumbering on given domain. More...

virtual int	forceEquationNumbering ()
	Forces equation renumbering on all domains associated to engng model. More...

Public Member Functions inherited from oofem::EngngModel
	EngngModel (int i, EngngModel *_master=NULL)
	Constructor. More...

virtual	~EngngModel ()
	Destructor. More...

	EngngModel (const EngngModel &)=delete

EngngModel &	operator= (const EngngModel &)=delete

Domain *	giveDomain (int n)
	Service for accessing particular problem domain. More...

void	setDomain (int i, Domain *ptr, bool iDeallocateOld=true)
	Sets i-th domain of receiver. More...

int	giveNumberOfDomains ()
	Returns number of domains in problem. More...

const std::string &	giveDescription () const

const time_t &	giveStartTime ()

bool	giveSuppressOutput () const

virtual ErrorEstimator *	giveDomainErrorEstimator (int n)
	Service for accessing ErrorEstimator corresponding to particular domain. More...

void	setNumberOfEquations (int id, int neq)

FILE *	giveOutputStream ()
	Returns file descriptor of output file. More...

std::string	giveOutputBaseFileName ()
	Returns base output file name to which extensions, like .out .vtu .osf should be added. More...

std::string	giveReferenceFileName ()
	Returns reference file name. More...

void	letOutputBaseFileNameBe (const std::string &src)
	Sets the base output file name. More...

ContextOutputMode	giveContextOutputMode ()
	Returns domain context output mode. More...

int	giveContextOutputStep ()
	Returns domain context output step. More...

void	setContextOutputMode (ContextOutputMode contextMode)
	Sets context output mode of receiver. More...

void	setUDContextOutputMode (int cStep)
	Sets user defined context output mode (it sets contextOutputMode to contextOutputMode), setting contextOutputStep to given value. More...

void	setProblemMode (problemMode pmode)
	Sets domain mode to given mode. More...

void	setParallelMode (bool newParallelFlag)
	Sets the problem to run in parallel (or not). More...

problemMode	giveProblemMode ()
	Returns domain mode. More...

void	setProblemScale (problemScale pscale)
	Sets scale in multiscale simulation. More...

problemScale	giveProblemScale ()
	Returns scale in multiscale simulation. More...

virtual void	setRenumberFlag ()
	Sets the renumber flag to true. More...

virtual void	resetRenumberFlag ()
	Sets the renumber flag to false. More...

double	giveSolutionStepTime ()
	Returns the user time of the current simulation step in seconds. More...

void	giveAnalysisTime (int &rhrs, int &rmin, int &rsec, int &uhrs, int &umin, int &usec)
	Returns the real and user time for the analysis. More...

void	terminateAnalysis ()
	Performs analysis termination after finishing analysis. More...

virtual void	solveYourself ()
	Starts solution process. More...

virtual void	terminate (TimeStep *tStep)
	Terminates the solution of time step. More...

virtual void	doStepOutput (TimeStep *tStep)
	Prints the ouput of the solution step (using virtual this->printOutputAtservice) to the stream detemined using this->giveOutputStream() method and calls exportModuleManager to do output. More...

void	saveStepContext (TimeStep *tStep, ContextMode mode)
	Saves context of given solution step, if required (determined using this->giveContextOutputMode() method). More...

virtual void	initializeYourself (TimeStep *tStep)
	Provides the opportunity to initialize state variables stored in element integration points according to initial conditions using function initializeYourself() on element level. More...

virtual int	initializeAdaptive (int tStepNumber)
	Initializes the newly generated discretization state according to previous solution. More...

virtual int	giveNumberOfDomainEquations (int di, const UnknownNumberingScheme &num)
	Returns number of equations for given domain in active (current time step) time step. More...

virtual FieldPtr	giveField (FieldType key, TimeStep *)
	Returns the smart pointer to requested field, Null otherwise. More...

EngngModel *	giveMasterEngngModel ()
	Returns the master engnmodel. More...

virtual double	giveLoadLevel ()
	Returns the current load level. More...

virtual double	giveEigenValue (int eigNum)
	Only relevant for eigen value analysis. Otherwise returns zero. More...

virtual void	setActiveVector (int i)
	Only relevant for eigen value analysis. Otherwise does noting. More...

int	updateSharedDofManagers (FloatArray &answer, const UnknownNumberingScheme &s, int ExchangeTag)
	Exchanges necessary remote DofManagers data. More...

int	exchangeRemoteElementData (int ExchangeTag)
	Exchanges necessary remote element data with remote partitions. More...

virtual int	giveCurrentNumberOfIterations ()
	Returns number of iterations that was required to reach equilibrium - used for adaptive step length in staggered problem. More...

MPI_Comm	giveParallelComm ()
	Returns the communication object of reciever. More...

int	packRemoteElementData (ProcessCommunicator &processComm)
	Packs data of local element to be received by their remote counterpart on remote partitions. More...

int	unpackRemoteElementData (ProcessCommunicator &processComm)
	Unpacks data for remote elements (which are mirrors of remote partition's local elements). More...

int	packDofManagers (ArrayWithNumbering *src, ProcessCommunicator &processComm)
	Packing function for vector values of DofManagers. More...

int	unpackDofManagers (ArrayWithNumbering *dest, ProcessCommunicator &processComm)
	Unpacking function for vector values of DofManagers . More...

ProblemCommunicator *	giveProblemCommunicator (EngngModelCommType t)

void	initializeCommMaps (bool forceInit=false)

virtual int	instanciateYourself (DataReader &dr, InputRecord ir, const char outFileName, const char *desc)
	Initializes whole problem according to its description stored in inputStream. More...

void	Instanciate_init ()
	Initialization of the receiver state (opening the default output stream, empty domain creation, initialization of parallel context, etc) before Initialization form DataReader. More...

int	instanciateDomains (DataReader &dr)
	Instanciate problem domains by calling their instanciateYourself() service. More...

int	instanciateMetaSteps (DataReader &dr)
	Instanciate problem meta steps by calling their instanciateYourself() service. More...

virtual int	instanciateDefaultMetaStep (InputRecord *ir)
	Instanciate default metastep, if nmsteps is zero. More...

virtual void	updateAttributes (MetaStep *mStep)
	Update receiver attributes according to step metaStep attributes. More...

void	initMetaStepAttributes (MetaStep *mStep)
	Update e-model attributes attributes according to step metaStep attributes. More...

MetaStep *	giveCurrentMetaStep ()
	Returns current meta step. More...

virtual TimeStep *	giveCurrentStep (bool force=false)
	Returns current time step. More...

virtual TimeStep *	givePreviousStep (bool force=false)
	Returns previous time step. More...

TimeStep *	generateNextStep ()
	Generate new time step (and associate metastep). More...

virtual void	preInitializeNextStep ()
	Does a pre-initialization of the next time step (implement if necessarry) More...

virtual int	giveNumberOfFirstStep (bool force=false)
	Returns number of first time step used by receiver. More...

int	giveNumberOfMetaSteps ()
	Return number of meta steps. More...

MetaStep *	giveMetaStep (int i)
	Returns the i-th meta step. More...

int	giveNumberOfSteps (bool force=false)
	Returns total number of steps. More...

virtual double	giveEndOfTimeOfInterest ()
	Returns end of time interest (time corresponding to end of time integration). More...

int	giveNumberOfTimeStepWhenIcApply ()
	Returns the time step number, when initial conditions should apply. More...

ExportModuleManager *	giveExportModuleManager ()
	Returns receiver's export module manager. More...

EngngModelTimer *	giveTimer ()
	Returns reference to receiver timer (EngngModelTimer). More...

virtual int	giveNewEquationNumber (int domain, DofIDItem)
	Increases number of equations of receiver's domain and returns newly created equation number. More...

virtual int	giveNewPrescribedEquationNumber (int domain, DofIDItem)
	Increases number of prescribed equations of receiver's domain and returns newly created equation number. More...

std::string	giveContextFileName (int tStepNumber, int stepVersion) const
	Returns the filename for the context file for the given step and version. More...

std::string	giveDomainFileName (int domainNum, int domainSerNum) const
	Returns the filename for the given domain (used by adaptivity and restore) More...

virtual void	initStepIncrements ()
	Initializes solution of new time step. More...

virtual bool	requiresEquationRenumbering (TimeStep *tStep)
	Returns true if equation renumbering is required for given solution step. More...

virtual ParallelContext *	giveParallelContext (int n)
	Returns the parallel context corresponding to given domain (n) and unknown type Default implementation returns i-th context from parallelContextList. More...

virtual void	initParallelContexts ()
	Creates parallel contexts. More...

virtual void	assemble (SparseMtrx &answer, TimeStep tStep, const MatrixAssembler &ma, const UnknownNumberingScheme &s, Domain domain)
	Assembles characteristic matrix of required type into given sparse matrix. More...

virtual void	assemble (SparseMtrx &answer, TimeStep tStep, const MatrixAssembler &ma, const UnknownNumberingScheme &r_s, const UnknownNumberingScheme &c_s, Domain domain)
	Assembles characteristic matrix of required type into given sparse matrix. More...

void	assembleVector (FloatArray &answer, TimeStep tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain domain, FloatArray *eNorms=NULL)
	Assembles characteristic vector of required type from dofManagers, element, and active boundary conditions, into given vector. More...

void	assembleVectorFromDofManagers (FloatArray &answer, TimeStep tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain domain, FloatArray *eNorms=NULL)
	Assembles characteristic vector of required type from dofManagers into given vector. More...

void	assembleVectorFromElements (FloatArray &answer, TimeStep tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain domain, FloatArray *eNorms=NULL)
	Assembles characteristic vector of required type from elements into given vector. More...

void	assembleVectorFromBC (FloatArray &answer, TimeStep tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain domain, FloatArray *eNorms=NULL)
	Assembles characteristic vector of required type from boundary conditions. More...

void	assembleExtrapolatedForces (FloatArray &answer, TimeStep tStep, CharType type, Domain domain)
	Assembles the extrapolated internal forces vector, useful for obtaining a good initial guess in nonlinear analysis with Dirichlet boundary conditions. More...

void	assemblePrescribedExtrapolatedForces (FloatArray &answer, TimeStep tStep, CharType type, Domain domain)

void	assembleVectorFromContacts (FloatArray &answer, TimeStep tStep, CharType type, ValueModeType mode, const UnknownNumberingScheme &s, Domain domain, FloatArray *eNorms=NULL)

virtual int	checkProblemConsistency ()
	Allows programmer to test problem its internal data, before computation begins. More...

virtual void	init ()
	Initializes the receiver state. More...

virtual void	postInitialize ()
	Performs post-initialization for all the problem contents (which is called after initializeFrom). More...

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

virtual void	printOutputAt (FILE file, TimeStep tStep, const IntArray &nodeSets, const IntArray &elementSets)

void	outputNodes (FILE file, Domain &domain, TimeStep tStep, int setNum)
	Outputs all nodes in the given set. More...

void	outputElements (FILE file, Domain &domain, TimeStep tStep, int setNum)
	Outputs all elements in the given set. More...

void	printYourself ()
	Prints state of receiver. Useful for debugging. More...

virtual int	useNonlocalStiffnessOption ()
	Returns nonzero if nonlocal stiffness option activated. More...

bool	isParallel () const
	Returns true if receiver in parallel mode. More...

int	giveRank () const
	Returns domain rank in a group of collaborating processes (0..groupSize-1) More...

int	giveNumberOfProcesses () const
	Returns the number of collaborating processes. More...

EngngModelContext *	giveContext ()
	Context requesting service. More...

virtual int	giveNumberOfSlaveProblems ()
	Returns number of slave problems. More...

virtual EngngModel *	giveSlaveProblem (int i)
	Returns i-th slave problem. More...

virtual int	estimateMaxPackSize (IntArray &commMap, DataStream &buff, int packUnpackType)
	Determines the space necessary for send/receive buffer. More...

virtual void	balanceLoad (TimeStep *tStep)
	Recovers the load balance between processors, if needed. More...

virtual LoadBalancer *	giveLoadBalancer ()
	Returns reference to receiver's load balancer. More...

virtual LoadBalancerMonitor *	giveLoadBalancerMonitor ()
	Returns reference to receiver's load balancer monitor. More...

void	initParallel ()
	Request domain rank and problem size. More...

EngngModel *	giveEngngModel ()
	Returns reference to itself -> required by communicator.h. More...

virtual bool	isElementActivated (int elemNum)

virtual bool	isElementActivated (Element *e)

virtual void	drawYourself (oofegGraphicContext &gc)

virtual void	drawElements (oofegGraphicContext &gc)

virtual void	drawNodes (oofegGraphicContext &gc)

virtual void	showSparseMtrxStructure (int type, oofegGraphicContext &gc, TimeStep *tStep)
	Shows the sparse structure of required matrix, type == 1 stiffness. More...

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

Protected Member Functions
void	updateInternalState (TimeStep *tStep)

void	applyIC (TimeStep *tStep)

void	evaluateElementStabilizationCoeffs (TimeStep *tStep)

void	updateElementsForNewInterfacePosition (TimeStep *tStep)

void	updateDofUnknownsDictionary_predictor (TimeStep *tStep)

void	updateDofUnknownsDictionary_corrector (TimeStep *tStep)

void	updateSolutionVectors (FloatArray &solutionVector, FloatArray &accelerationVector, FloatArray &incrementalSolutionVector, TimeStep *tStep)

void	updateSolutionVectors_predictor (FloatArray &solutionVector, FloatArray &accelerationVector, TimeStep *tStep)

void	updateDofManVals (TimeStep *tStep)

Protected Member Functions inherited from oofem::EngngModel
virtual void	packMigratingData (TimeStep *tStep)
	Packs receiver data when rebalancing load. More...

virtual void	unpackMigratingData (TimeStep *tStep)
	Unpacks receiver data when rebalancing load. More...

Protected Attributes
std::unique_ptr< SparseLinearSystemNM >	nMethod
	Numerical method used to solve the problem. More...

LinSystSolverType	solverType

SparseMtrxType	sparseMtrxType

std::unique_ptr< SparseMtrx >	lhs

std::unique_ptr< PrimaryField >	VelocityPressureField

FloatArray	accelerationVector

FloatArray	incrementalSolutionVector

FloatArray	internalForces

FloatArray	eNorm

double	deltaT

int	deltaTF

double	atolv
	Convergence tolerance. More...

double	rtolv

int	maxiter
	Max number of iterations. More...

bool	stopmaxiter
	Flag if set to true (default), then when max number of iteration reached, computation stops otherwise computation continues with next step. More...

double	alpha
	Integration constant. More...

int	initFlag

int	consistentMassFlag

bool	equationScalingFlag

double	lscale

double	uscale

double	dscale

double	Re
	Reynolds number. More...

std::unique_ptr< MaterialInterface >	materialInterface

Protected Attributes inherited from oofem::EngngModel
int	ndomains
	Number of receiver domains. More...

std::vector< std::unique_ptr< Domain > >	domainList
	List of problem domains. More...

int	numberOfSteps
	Total number of time steps. More...

int	numberOfEquations
	Total number of equation in current time step. More...

int	numberOfPrescribedEquations
	Total number or prescribed equations in current time step. More...

IntArray	domainNeqs
	Number of equations per domain. More...

IntArray	domainPrescribedNeqs
	Number of prescribed equations per domain. More...

bool	renumberFlag
	Renumbering flag (renumbers equations after each step, necessary if Dirichlet BCs change). More...

bool	profileOpt
	Profile optimized numbering flag (using Sloan's algorithm). More...

int	equationNumberingCompleted
	Equation numbering completed flag. More...

int	nMetaSteps
	Number of meta steps. More...

std::vector< MetaStep >	metaStepList
	List of problem metasteps. More...

std::unique_ptr< TimeStep >	stepWhenIcApply
	Solution step when IC (initial conditions) apply. More...

std::unique_ptr< TimeStep >	currentStep
	Current time step. More...

std::unique_ptr< TimeStep >	previousStep
	Previous time step. More...

int	number
	Receivers id. More...

std::string	dataOutputFileName
	Path to output stream. More...

std::string	coreOutputFileName
	String with core output file name. More...

FILE *	outputStream
	Output stream. More...

std::string	referenceFileName
	String with reference file name. More...

ContextOutputMode	contextOutputMode
	Domain context output mode. More...

int	contextOutputStep

ExportModuleManager *	exportModuleManager
	Export module manager. More...

InitModuleManager *	initModuleManager
	Initialization module manager. More...

problemMode	pMode
	Domain mode. More...

problemScale	pScale
	Multiscale mode. More...

time_t	startTime
	Solution start time. More...

EngngModel *	master
	Master e-model; if defined receiver is in maintained (slave) mode. More...

EngngModelContext *	context
	Context. More...

EngngModelTimer	timer
	E-model timer. More...

int	parallelFlag
	Flag indicating that the receiver runs in parallel. More...

enum fMode	nonLinFormulation
	Type of non linear formulation (total or updated formulation). More...

ErrorEstimator *	defaultErrEstimator
	Error estimator. Useful for adaptivity, or simply printing errors output. More...

int	rank
	Domain rank in a group of collaborating processes (0..groupSize-1). More...

int	numProcs
	Total number of collaborating processes. More...

int	nonlocalExt
	Flag indicating if nonlocal extension active, which will cause data to be sent between shared elements before computing the internal forces. More...

char	processor_name [PROCESSOR_NAME_LENGTH]
	Processor name. More...

MPI_Comm	comm
	Communication object for this engineering model. More...

CommunicatorBuff *	commBuff
	Common Communicator buffer. More...

ProblemCommunicator *	communicator
	Communicator. More...

ProblemCommunicator *	nonlocCommunicator
	NonLocal Communicator. Necessary when nonlocal constitutive models are used. More...

std::vector< ParallelContext >	parallelContextList
	List where parallel contexts are stored. More...

bool	suppressOutput
	Flag for suppressing output to file. More...

std::string	simulationDescription

LoadBalancer *	lb
	Load Balancer. More...

LoadBalancerMonitor *	lbm

bool	loadBalancingFlag
	If set to true, load balancing is active. More...

bool	force_load_rebalance_in_first_step
	Debug flag forcing load balancing after first step. More...

Additional Inherited Members
Public Types inherited from oofem::EngngModel
enum	EngngModel_UpdateMode { EngngModel_Unknown_Mode, EngngModel_SUMM_Mode, EngngModel_SET_Mode }

enum	EngngModelCommType { PC_default, PC_nonlocal }

enum	InitialGuess { IG_None = 0, IG_Tangent = 1 }
	Means to choose methods for finding a good initial guess. More...

Protected Types inherited from oofem::EngngModel
enum	{ InternalForcesExchangeTag, MassExchangeTag, LoadExchangeTag, ReactionExchangeTag, RemoteElementExchangeTag }
	Message tags. More...

Detailed Description

This class represents transient incompressible flow problem.

Solution is based on algorithm with SUPG/PSPG stabilization.

Definition at line 104 of file supg.h.

Constructor & Destructor Documentation

oofem::SUPG::SUPG	(	int	i,
		EngngModel *	_master = `NULL`
	)

Definition at line 210 of file supg.C.

References consistentMassFlag, dscale, equationScalingFlag, initFlag, lscale, oofem::EngngModel::ndomains, and uscale.

oofem::SUPG::~SUPG ( )

virtual

Definition at line 220 of file supg.C.

Member Function Documentation

void oofem::SUPG::applyIC ( TimeStep * tStep )

protected

Definition at line 935 of file supg.C.

References accelerationVector, oofem::FloatArray::at(), oofem::Domain::giveDofManager(), oofem::EngngModel::giveDomain(), oofem::Domain::giveElements(), oofem::Domain::giveNumberOfDofManagers(), oofem::EngngModel::giveNumberOfDomainEquations(), materialInterface, OOFEM_LOG_INFO, requiresUnknownsDictionaryUpdate(), oofem::FloatArray::resize(), updateElementsForNewInterfacePosition(), oofem::SUPGElement::updateInternalState(), oofem::Element::updateYourself(), VelocityPressureField, and oofem::FloatArray::zero().

Referenced by solveYourselfAt().

int oofem::SUPG::checkConsistency ( )

virtual

Allows programmer to test some receiver's internal data, before computation begins.

Returns: Nonzero if receiver check is o.k.

Reimplemented from oofem::EngngModel.

Definition at line 862 of file supg.C.

References oofem::EngngModel::checkConsistency(), equationScalingFlag, oofem::ForceLoadBVT, oofem::Domain::giveBcs(), oofem::EngngModel::giveDomain(), oofem::Domain::giveElements(), oofem::Domain::giveIcs(), giveVariableScale(), OOFEM_ERROR, OOFEM_WARNING, oofem::PressureBVT, uscale, oofem::VelocityBVT, oofem::VST_Force, and oofem::VST_Pressure.

void oofem::SUPG::evaluateElementStabilizationCoeffs ( TimeStep * tStep )

protected

Definition at line 1020 of file supg.C.

References oofem::EngngModel::giveDomain(), oofem::Domain::giveElements(), and oofem::FMElement::updateStabilizationCoeffs().

Referenced by solveYourselfAt().

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

inlinevirtual

Returns class name of the receiver.

Implements oofem::EngngModel.

Definition at line 177 of file supg.h.

virtual bool oofem::SUPG::giveEquationScalingFlag ( )

inlinevirtual

Returns the Equation scaling flag, which is used to indicate that governing equation(s) are scaled, or non-dimensionalized.

Reimplemented from oofem::EngngModel.

Definition at line 186 of file supg.h.

virtual fMode oofem::SUPG::giveFormulation ( )

inlinevirtual

Indicates type of non linear computation (total or updated formulation).

This is used for example on Nodal level to update coordinates if updated formulation is done, or on element level, when non linear contributions are computed.

Reimplemented from oofem::EngngModel.

Definition at line 180 of file supg.h.

References oofem::TL.

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

inlinevirtual

Definition at line 178 of file supg.h.

References _IFT_SUPG_Name.

virtual MaterialInterface* oofem::SUPG::giveMaterialInterface ( int n )

inlinevirtual

Returns material interface representation for given domain.

Reimplemented from oofem::EngngModel.

Definition at line 192 of file supg.h.

TimeStep * oofem::SUPG::giveNextStep ( )

virtual

Returns next time step (next to current step) of receiver.

Reimplemented from oofem::EngngModel.

Definition at line 415 of file supg.C.

References oofem::EngngModel::currentStep, deltaT, deltaTF, oofem::Function::evaluateAtTime(), oofem::EngngModel::giveDomain(), oofem::Domain::giveElements(), oofem::Domain::giveFunction(), giveSolutionStepWhenIcApply(), giveVariableScale(), materialInterface, oofem::min(), OOFEM_LOG_INFO, oofem::EngngModel::previousStep, and oofem::VST_Time.

NumericalMethod * oofem::SUPG::giveNumericalMethod ( MetaStep * mStep )

virtual

Returns reference to receiver's numerical method.

Reimplemented from oofem::EngngModel.

Definition at line 225 of file supg.C.

References oofem::classFactory, oofem::ClassFactory::createSparseLinSolver(), nMethod, OOFEM_ERROR, and solverType.

Referenced by solveYourselfAt(), and updateDomainLinks().

double oofem::SUPG::giveReynoldsNumber ( )

virtual

Implements oofem::FluidModel.

Definition at line 387 of file supg.C.

References equationScalingFlag, and Re.

TimeStep * oofem::SUPG::giveSolutionStepWhenIcApply ( bool force = false )

virtual

Returns the solution step when Initial Conditions (IC) apply.

Parameters

force when set to true then receiver reply is returned instead of master (default)

Reimplemented from oofem::EngngModel.

Definition at line 398 of file supg.C.

References deltaT, oofem::EngngModel::giveNumberOfTimeStepWhenIcApply(), oofem::EngngModel::giveSolutionStepWhenIcApply(), giveVariableScale(), oofem::EngngModel::master, oofem::EngngModel::stepWhenIcApply, and oofem::VST_Time.

Referenced by giveNextStep().

double oofem::SUPG::giveUnknownComponent	(	ValueModeType	,
		TimeStep *	,
		Domain *	,
		Dof *
	)

virtual

Returns requested unknown.

Unknown at give time step is characterized by its type and mode and by its equation number. This function is used by Dofs, when they are requested for their associated unknowns.

See also: Dof::giveUnknown

Reimplemented from oofem::EngngModel.

Definition at line 319 of file supg.C.

References oofem::Dof::__giveEquationNumber(), oofem::__ValueModeTypeToString(), accelerationVector, oofem::FloatArray::at(), giveUnknownDictHashIndx(), oofem::Dof::giveUnknowns(), OOFEM_ERROR, requiresUnknownsDictionaryUpdate(), and VelocityPressureField.

int oofem::SUPG::giveUnknownDictHashIndx	(	ValueModeType	mode,
		TimeStep *	tStep
	)

virtual

This method is responsible for computing unique dictionary id (ie hash value) from given valueModeType and time step.

This function is used by particular dofs to access unknown identified by given parameters from its dictionary using computed index. Usually the hash algorithm should produce index that depend on time step relatively to actual one to avoid storage of complete history.

Reimplemented from oofem::EngngModel.

Definition at line 1128 of file supg.C.

References oofem::EngngModel::giveCurrentStep(), oofem::TimeStep::giveNumber(), oofem::EngngModel::givePreviousStep(), and OOFEM_ERROR.

Referenced by giveUnknownComponent().

double oofem::SUPG::giveVariableScale ( VarScaleType varId )

virtual

Returns the scale factor for given variable type.

Reimplemented from oofem::EngngModel.

Definition at line 996 of file supg.C.

References dscale, lscale, OOFEM_ERROR, uscale, oofem::VST_Density, oofem::VST_Force, oofem::VST_Length, oofem::VST_Pressure, oofem::VST_Time, oofem::VST_Velocity, and oofem::VST_Viscosity.

Referenced by checkConsistency(), giveNextStep(), and giveSolutionStepWhenIcApply().

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

virtual

Initializes receiver according to object description in input reader.

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

Reimplemented from oofem::EngngModel.

Definition at line 237 of file supg.C.

References _IFT_EngngModel_lstype, _IFT_EngngModel_smtype, _IFT_SUPG_alpha, _IFT_SUPG_atolv, _IFT_SUPG_cmflag, _IFT_SUPG_deltat, _IFT_SUPG_deltatFunction, _IFT_SUPG_dscale, _IFT_SUPG_lscale, _IFT_SUPG_maxiter, _IFT_SUPG_miflag, _IFT_SUPG_rtolv, _IFT_SUPG_scaleflag, _IFT_SUPG_stopmaxiter, _IFT_SUPG_uscale, alpha, atolv, consistentMassFlag, deltaT, deltaTF, dscale, equationScalingFlag, oofem::EngngModel::giveContext(), oofem::EngngModel::giveDomain(), oofem::EngngModelContext::giveFieldManager(), oofem::InputRecord::hasField(), oofem::EngngModel::initializeFrom(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, oofem::IRRT_OK, lscale, materialInterface, maxiter, Re, oofem::FieldManager::registerField(), requiresUnknownsDictionaryUpdate(), rtolv, solverType, sparseMtrxType, stopmaxiter, uscale, and VelocityPressureField.

void oofem::SUPG::printDofOutputAt	(	FILE *	stream,
		Dof *	iDof,
		TimeStep *	tStep
	)

virtual

DOF printing routine.

Called by DofManagers to print Dof specific part. Dof class provides component printing routines, but emodel is responsible for what will be printed at DOF level.

Parameters

stream	output stream
iDof	dof to be processed
tStep	solution step

Reimplemented from oofem::EngngModel.

Definition at line 920 of file supg.C.

References dscale, oofem::Dof::giveDofID(), OOFEM_ERROR, oofem::Dof::printSingleOutputAt(), and uscale.

virtual int oofem::SUPG::requiresUnknownsDictionaryUpdate ( )

inlinevirtual

Indicates if EngngModel requires Dofs dictionaries to be updated.

If EngngModel does not support changes of static system, the dof forwards the requests for its unknowns to EngngModel, where unknowns are naturally kept. This is possible, because dof equation number is same during whole solution. But when changes of static system are allowed, several problem arise. For example by solving simple incremental static with allowed static changes, the incremental displacement vector of structure can not be added to total displacement vector of structure, because equation numbers may have changed, and one can not simply add these vector to obtain new total displacement vector, because incompatible displacement will be added. To solve this problem, unknown dictionary at dof level has been assumed. Dof then keeps its unknowns in its own private dictionary. After computing increment of solution, engngModel updates for each dof its unknowns in its dictionary (using updateUnknownsDictionary function). For aforementioned example engngModel updates incremental values but also total value by asking dof for previous total value (dof will use its dictionary, does not asks back EngngModel) adds corresponding increment and updates total value in dictionary.

Reimplemented from oofem::EngngModel.

Definition at line 184 of file supg.h.

Referenced by applyIC(), giveUnknownComponent(), initializeFrom(), solveYourselfAt(), updateDofUnknownsDictionary_corrector(), updateDofUnknownsDictionary_predictor(), and updateInternalState().

contextIOResultType oofem::SUPG::restoreContext	(	DataStream &	stream,
		ContextMode	mode
	)

virtual

Restores the state of model from output stream.

Restores not only the receiver state, but also same function is invoked for all DofManagers and Elements in associated domain. Note that by restoring element context also contexts of all associated integration points (and material statuses) are restored. Each context is associated with unique time step. Only one context per time step is allowed. Restore context function will restore such context, which is related (through its step number) to time step number and version given in obj parameter. Restoring context will change current time step in order to correspond to newly restored context.

Parameters

stream	Context file.
mode	Determines amount of info in stream.

Returns: contextIOResultType.

Exceptions

ContextIOERR exception if error encountered.

Reimplemented from oofem::EngngModel.

Definition at line 835 of file supg.C.

References accelerationVector, oofem::CIO_OK, materialInterface, oofem::EngngModel::restoreContext(), oofem::FloatArray::restoreYourself(), THROW_CIOERR, and VelocityPressureField.

contextIOResultType oofem::SUPG::saveContext	(	DataStream &	stream,
		ContextMode	mode
	)

virtual

Stores the state of model to output stream.

Stores not only the receiver state, but also same function is invoked for all DofManagers and Elements in associated domain. Note that by storing element context also contexts of all associated integration points (and material statuses) are stored.

Parameters

stream	Context stream.
mode	Determines amount of info in stream.

Returns: contextIOResultType.

Exceptions

ContextIOERR If error encountered.

Reimplemented from oofem::EngngModel.

Definition at line 808 of file supg.C.

References accelerationVector, oofem::CIO_OK, materialInterface, oofem::EngngModel::saveContext(), oofem::FloatArray::storeYourself(), THROW_CIOERR, and VelocityPressureField.

void oofem::SUPG::solveYourselfAt ( TimeStep * tStep )

virtual

Solves problem for given time step.

Should assemble characteristic matrices and vectors if necessary and solve problem using appropriate numerical method. After finishing solution, this->updateYourself function for updating solution state and then this->terminate function (for updating nodal and element values) should be called.

Reimplemented from oofem::EngngModel.

Definition at line 453 of file supg.C.

void oofem::SUPG::updateComponent	(	TimeStep *	tStep,
		NumericalCmpn	cmpn,
		Domain *	d
	)

virtual

Updates components mapped to numerical method if necessary during solution process.

Some numerical methods may require updating mapped components during solution process (e.g., updating of tangent stiffness when using updated Newton-Raphson method).

Parameters

tStep	Time when component is updated.
cmpn	Numerical component to update.
d	Domain.

Reimplemented from oofem::EngngModel.

Definition at line 357 of file supg.C.

References alpha, oofem::EngngModel::assemble(), oofem::EngngModel::assembleVector(), dscale, eNorm, oofem::Domain::giveElements(), internalForces, oofem::EngngModel::InternalForcesExchangeTag, oofem::InternalRhs, lhs, lscale, oofem::NonLinearLhs, OOFEM_ERROR, oofem::EngngModel::updateSharedDofManagers(), uscale, and oofem::FloatArray::zero().

Referenced by solveYourselfAt().

void oofem::SUPG::updateDofManVals ( TimeStep * tStep )

protected

void oofem::SUPG::updateDofUnknownsDictionary	(	DofManager *	,
		TimeStep *
	)

virtual

Updates necessary values in Dofs unknown dictionaries.

See also: EngngModel::requiresUnknownsDictionaryUpdate; Dof::updateUnknownsDictionary

Reimplemented from oofem::EngngModel.

Definition at line 1046 of file supg.C.

Referenced by updateInternalState().

void oofem::SUPG::updateDofUnknownsDictionary_corrector ( TimeStep * tStep )

protected

Definition at line 1096 of file supg.C.

References alpha, oofem::FloatArray::at(), deltaT, oofem::Domain::giveDofManager(), oofem::EngngModel::giveDomain(), oofem::Domain::giveNumberOfDofManagers(), oofem::TimeStep::giveTimeIncrement(), incrementalSolutionVector, and requiresUnknownsDictionaryUpdate().

Referenced by solveYourselfAt().

void oofem::SUPG::updateDofUnknownsDictionary_predictor ( TimeStep * tStep )

protected

Definition at line 1060 of file supg.C.

References deltaT, oofem::Domain::giveDofManager(), oofem::EngngModel::giveDomain(), oofem::Domain::giveNumberOfDofManagers(), oofem::TimeStep::givePreviousStep(), oofem::TimeStep::giveTimeIncrement(), and requiresUnknownsDictionaryUpdate().

Referenced by solveYourselfAt().

void oofem::SUPG::updateDomainLinks ( )

virtual

Updates domain links after the domains of receiver have changed.

Used mainly after restoring context - the domains may change and this service is then used to update domain variables in all components belonging to receiver like error estimators, solvers, etc, having domains as attributes.

Reimplemented from oofem::EngngModel.

Definition at line 912 of file supg.C.

References oofem::EngngModel::giveCurrentMetaStep(), giveNumericalMethod(), oofem::NumericalMethod::reinitialize(), and oofem::EngngModel::updateDomainLinks().

void oofem::SUPG::updateElementsForNewInterfacePosition ( TimeStep * tStep )

protected

Definition at line 1031 of file supg.C.

References oofem::EngngModel::giveDomain(), oofem::Domain::giveElements(), OOFEM_LOG_DEBUG, and oofem::SUPGElement::updateElementForNewInterfacePosition().

Referenced by applyIC(), and solveYourselfAt().

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

protected

Definition at line 791 of file supg.C.

References oofem::EngngModel::domainList, requiresUnknownsDictionaryUpdate(), and updateDofUnknownsDictionary().

Referenced by updateYourself().

void oofem::SUPG::updateSolutionVectors	(	FloatArray &	solutionVector,
		FloatArray &	accelerationVector,
		FloatArray &	incrementalSolutionVector,
		TimeStep *	tStep
	)

protected

Definition at line 1186 of file supg.C.

References oofem::FloatArray::add(), alpha, oofem::FloatArray::at(), deltaT, oofem::Domain::giveDofManagers(), oofem::EngngModel::giveDomain(), oofem::Domain::giveElements(), and oofem::TimeStep::giveTimeIncrement().

Referenced by solveYourselfAt().

void oofem::SUPG::updateSolutionVectors_predictor	(	FloatArray &	solutionVector,
		FloatArray &	accelerationVector,
		TimeStep *	tStep
	)

protected

Definition at line 1140 of file supg.C.

References oofem::FloatArray::at(), deltaT, oofem::Domain::giveDofManagers(), oofem::EngngModel::giveDomain(), oofem::Domain::giveElements(), and oofem::TimeStep::giveTimeIncrement().

Referenced by solveYourselfAt().

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

virtual

Updates internal state after finishing time step.

(for example total values may be updated according to previously solved increments). Then element values are also updated (together with related integration points and material statuses).

Reimplemented from oofem::EngngModel.

Definition at line 778 of file supg.C.

References materialInterface, updateInternalState(), and oofem::EngngModel::updateYourself().